estructura gravimétrica y magnética de la corteza del suroeste ...

Nieves Sánchez Jiménez
del mapa. En la esquina inferior izquierda se alcanzan los máximos valores de anomalía de
Bouguer, 88 mGal. En la parte suroriental del mapa se observa más claramente el mínimo
que ya se apreciaba en el Mapa de Anomalías presentado en el capítulo 4 (fig. 4.9). En la
figura 6.3 se aprecia más claramente la magnitud de este mínimo, alcanzando los -122 mGal.
Se puede observar también cómo la isolínea de valor 0 mGal, con una dirección NE-SO en la
parte meridional del mapa, gira en las cercanías de Córdoba y toma la orientación de las
estructuras variscas en dirección NO-SE al adentrarse en el Macizo Ibérico. Esta isolínea
sigue, al menos en la parte sur, los afloramientos paleozoicos del Macizo Ibérico siendo
ligeramente oblicua al contacto de éstos con los sedimentos de la cuenca del Guadalquivir.
Este mapa presenta la dificultad, frente al análisis espectral, de no ser cuadrado, lo
que implica que un mismo número de armónicos va a representar longitudes de onda
diferentes en las direcciones x e y. El programa OASIS utilizado para el cálculo del espectro
de potencia radial del mapa, permite expandir el mapa en los bordes de forma que la malla
sea periódica en ellos, es decir, de forma que un punto en el borde izquierdo de la malla
encaje con el borde derecho y que un punto en la parte superior de la malla encaje con un
punto situado en la parte inferior, tanto en valor como en pendiente. Para ello, añade nodos
en la malla hasta hacerla cuadrada, o rectangular, y rellena estos nodos con valores
interpolados utilizando datos de zonas válidas de la malla. Así, consideramos todas las
longitudes de onda de las anomalías presentes en el mapa contenidas totalmente en él.
Para comprobar la influencia de la forma del mapa en la estimación de profundidades
a partir de la gráfica del espectro radial, se ha realizado el análisis espectral en estos dos
casos: expandiendo la malla hasta conseguir un mapa cuadrado y expandiendo la malla un 10
% en todas las direcciones manteniendo un mapa rectangular.
ln espectro de potencia radial
110
6
4
2
0
-2
-4
-6
A. Expansión a una
superfice cuadrada
Regional = f (0-0,016)
Prof.= 34,53 ± 2,21 km
Residual = f (0,016-0,045)
Prof.= 11,59 ± 0,50 km
Ruido blanco
0 0.02 0.04 0.06 0.08 0.1
frecuencia (km -1 )
ln espectro de potencia radial
6
4
2
0
-2
-4
-6
B. Expansión a una
superficie rectangular
Regional = f (0-0,016)
Prof.= 33,17 ± 1,62 km
Residual = f (0,016-0,045)
Prof.= 12,66 ± 0,36 km
Ruido blanco
0 0.02 0.04 0.06 0.08 0.1
frecuencia (km -1 )
Figura 6. 4. Representación del espectro de potencia radial para el mapa gravimétrico con una expansión a una
superficie cuadrada (450 x 450 km) (A) y a una superficie rectangular (450 x 240 km) (B).