Descarga la aplicación para disfrutar aún más
Esta es una vista previa del archivo. Inicie sesión para ver el archivo original
Av. Hidalgo 935, Colonia Centro, C.P. 44100, Guadalajara, Jalisco, México bibliotecadigital@redudg.udg.mx - Tel. 31 34 22 77 ext. 11959 UNIVERSIDAD DE GUADALAJARA COORDINACIÓN GENERAL ACADÉMICA Coordinación de Bibliotecas Biblioteca Digital La presente tesis es publicada a texto completo en virtud de que el autor ha dado su autorización por escrito para la incorporación del documento a la Biblioteca Digital y al Repositorio Institucional de la Universidad de Guadalajara, esto sin sufrir menoscabo sobre sus derechos como autor de la obra y los usos que posteriormente quiera darle a la misma. UNIVERSIDAD DE GUADALAJARA CENTRO UNIVERSITARIO DE CIENCIAS EXACTAS E INGENIERIAS “ANÁLISIS EN LA FRONTERA ENTRE LA SIERRA MADRE OCCIDENTAL Y EL EJE NEOVOLCÁNICO (W DE MÉXICO) MEDIANTE INTERPRETACIÓN 2D DE ANOMALÍAS AEROMAGNÉTICAS.” T E S I S Q U E P A R A O B T E N E R E L T I T U L O D E L I C E N C I A D O E N F Í S I C A P R E S E N T A R O D R I G O N A V A R R O L E Ó N Director: Dr. Miguel Angel Alatorre Zamora INVIERNO 2019 UNIVERSIDAD DE GUADALAJARA CENTRO UNIVERSITARIO DE CIENCIAS EXACTAS E INGENIERIAS “ANÁLISIS EN LA FRONTERA ENTRE LA SIERRA MADRE OCCIDENTAL Y EL EJE NEOVOLCÁNICO (W DE MÉXICO) MEDIANTE INTERPRETACIÓN 2D DE ANOMALÍAS AEROMAGNÉTICAS.” T E S I S Q U E P A R A O B T E N E R E L T I T U L O D E L I C E N C I A D O E N F Í S I C A P R E S E N T A R O D R I G O N A V A R R O L E Ó N Director: Dr. Miguel Angel Alatorre Zamora INVIERNO 2019 1 Comité de Titulación Licenciatura en Física Centro Universitario de Ciencias Exactas e Ingenierías Universidad de Guadalajara Estimados Miembros del Comité Por medio de la presente hago de su conocimiento que, después de revisar la tesis profesional del pasante de la Licenciatura en Física C. RODRIGO NAVARRO LEÓN que lleva por título “ANALISIS EN LA FRONTERA ENTRE LA SIERRA MADRE OCCIDENTAL Y EL EJE NEOVOLCÁNICO (W DE MÉXICO) MEDIANTE INTERPRETACIÓN 2D DE ANOMALÍAS AEROMAGNÉTICAS”, considero que la misma cumple con los objetivos presentados en el Art. 13 del reglamento de titulación, así como las formas exigidas por tan respetable Comité, por lo que otorgo mi autorización para su impresión y defensa. A T E N T A M E N T E DR. MIGUEL ANGEL ALATORRE ZAMORA DIRECTOR DE TRABAJO DE TITULACIÓN AGRADECIMIENTOS Al Dr. Miguel Angel Alatorre por todo el apoyo y disponibilidad que me brindo durante los últimos años, por su gran esfuerzo y dedicación hacia todos sus estudiantes. Gracias por ser un ejemplo a seguir. A mis padres por toda la motivación y apoyo, por enseñarme desde pequeño la importancia del conocimiento, pero sobre todo gracias por el gran sacrificio que han hecho para ayudarme a llegar aquí. A los sinodales por el tiempo que dedicaron a revisar este trabajo, así como por todas sus sugerencias, las cuales ayudan al enriquecimiento de mi tesis. A la Universidad de Guadalajara por todo lo que eh aprendido a lo largo de mi formación Universitaria. INDICE RESUMEN ........................................................................................................................... 1 1.- INTRODUCCIÓN ............................................................................................................ 2 1.1 OBJETIVOS ................................................................................................................ 3 2. - ZONA DE ESTUDIO ...................................................................................................... 4 2.1 GEOLOGÍA DEL ÁREA DE ESTUDIO ........................................................................ 5 2.2 PROVINCIAS FISIOGRÁFICAS .................................................................................. 8 2.3 ESTRUCTURAS TECTÓNICAS ASOCIADAS A LA ZONA DE ESTUDIO ................ 9 2.3.1 SIERRA MADRE OCCIDENTAL ......................................................................... 10 2.3.2 EJE NEOVOLCANICO ........................................................................................ 11 2.3.3 BLOQUE JALISCO.............................................................................................. 12 2.3.4 RIFT TEPIC-ZACOALCO .................................................................................... 13 3.- ANTECEDENTES ......................................................................................................... 14 3.1 FUNDAMENTO FÍSICO ............................................................................................ 14 3.2 EL CAMPO GEOMAGNÉTICO ................................................................................. 15 3.3 MAGNETOMETRIA ................................................................................................... 18 3.4 ANOMALIAS MAGNETICAS ..................................................................................... 20 3.5 FILTRADO DE ANOMALIAS MAGNETICA ............................................................... 21 3.5.1 CONTINUACIÓN ASCENDENTE ....................................................................... 21 3.5.2 CONTINUACIÓN DESCENDENTE ..................................................................... 22 3.5.3 REDUCCIÓN AL POLO ...................................................................................... 23 3.6 SUCEPTIBILIDAD MAGNETICA ............................................................................... 23 4.- DESARROLLO DEL TEMA .......................................................................................... 25 4.1 METODOLOGÍA ........................................................................................................ 25 4.2 DATOS ...................................................................................................................... 25 4.2.1 INFORMACIÓN DE LA CARTA MAGNÉTICA DE CAMPO TOTAL .................... 25 4.2.2 INFORMACIÓN DE LA ELEVACIÓN .................................................................. 26 4.3 PERFILES ................................................................................................................. 26 4.5 INTERPRETACIÓN DIRECTA .................................................................................. 28 5.- RESULTADOS Y DISCUSIÓN ..................................................................................... 30 5.1 CORRELACION VISUAL ENTRE ANOMALÍAS MAGNÉTICAS Y TOPOGRAFÍA ... 30 5.2 CONTINUACIONES ASCENDENTES DE PERFILES AEROMAGNÉTICOS ........... 39 5.3 PERFILES GEOLOGICOS 2D OBTENIDOS POR INTERPRETACIÓN DIRECTA .. 44 5.2.1 MODELO DEL PERFIL P-01 ............................................................................... 46 5.2.2 MODELO DEL PERFIL P-02 ............................................................................... 50 5.2.3 MODELO DEL PERFIL P-04 ............................................................................... 53 5.2.4 MODELO DEL PERFIL P-06 ............................................................................... 56 5.2.5 MODELO DEL PERFIL P-08 ............................................................................... 59 6.- CONCLUSIONES Y PERSPECTIVA ............................................................................ 63 BIBLIOGRAFÍA ................................................................................................................. 66 INDICE DE FIGURAS Figura 1: Ubicación de la Zona de estudio. ......................................................................... 4 Figura 2: Mapa divisional donde se indica el área de estudio (Obtenida del Servicio Geológico Mexicano [5]). ...................................................................................................... 5 Figura 3: Carta Geológica-Minera Tepic F13-8 (Obtenida del Servicio Geológico Mexicano [4]), escala 1:250,000. El trazo de líneas corresponde al sistema coordenado UTM. La simbología sobre la geología es detallada en la Figura 4. ................................................... 6 Figura 4: Explicación de la Carta Geológica-Minera Tepic F13-8 (Obtenida del Servicio Geológico Mexicano [4]) donde se muestra la Simbología y la Columna Geológica de las rocas presentes en la zona de estudio. ................................................................................ 7 Figura 5: Principales provincias fisiográficas de México que experimentaron cambios significativos en su relieve o configuración durante el Cenozoico, el recuadro rojo representa la zona de estudio (Modificado de [7]). ............................................................... 9 Figura 6: Configuración tectónica del centro-Oeste de México. Se observan las zonas de los Rift Tepic-Zacoalco, Chapala y Colima. El recuadro rojo representa la zona de estudio (Modificado de [8]). ............................................................................................................. 10 Figura 7: Campo geomagnético. En primera aproximación este puede ser considerado como el producido por un dipolo magnético. 1: Plano de la órbita de la Tierra, 2: Eje de rotación, 3: Eje magnético. (Obtenida de [23]). .................................................................. 16 Figura 8: Componente del campo geomagnético. (Obtenida de [23]). .............................. 18 Figura 9: Principio de la magnetometría (Tomado de [28]). .............................................. 19 Figura 10: Mapa de Campo Magnético Total (Modificado del Servicio Geológico Mexicano [5]) donde se muestra la zona de interés; así como los perfiles seleccionados. ................ 27 Figura 11: Valores de Intensidad de Campo Magnético, en nano Teslas, de la Figura 10 (Modificado del Servicio Geológico Mexicano [5]). ............................................................. 27 Figura 12: Mapa Geológico-Minero (modificado del Servicio Geológico Mexicano [4]) donde se muestra la zona de interés; así como los perfiles seleccionados y la geología presente. ............................................................................................................................ 28 Figura 13: Sección P-01; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 31 Figura 14: Sección P-02; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 32 Figura 15: Sección P-03; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 33 Figura 16: Sección P-04; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 34 Figura 17: Sección P-05; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 35 Figura 18: Sección P-06; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 36 Figura 19: Sección P-07; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 37 Figura 20: Sección P-08; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. ....... 38 Figura 21: Continuación Ascendente para el Perfil P-01; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 40 Figura 22: Continuación Ascendente para el Perfil P-02; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 40 Figura 23: Continuación Ascendente para el Perfil P-03; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 41 Figura 24: Continuación Ascendente para el Perfil P-04; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 41 Figura 25: Continuación Ascendente para el Perfil P-05; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 42 Figura 26: Continuación Ascendente para el Perfil P-06; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 42 Figura 27: Continuación Ascendente para el Perfil P-07; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 43 Figura 28: Continuación Ascendente para el Perfil P-08; cada 500 m en el intervalo de 200 a 3700 m. ........................................................................................................................... 43 Figura 29: P-01; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. .............................................................................................................. 46 Figura 30: P-01’; Sección meridional del perfil P-01 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ............... 47 Figura 31: P-01”; Sección septentrional del perfil P-01 (distancia de 40 - 81.6 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ........................................................................................................................................... 49 Figura 32: P-02; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. .............................................................................................................. 50 Figura 33: P-02’; Sección meridional del perfil P-02 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ............... 51 Figura 34: P-02”; Sección septentrional del perfil P-02 (distancia de 40 - 80.6 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ........................................................................................................................................... 52 Figura 35: P-04; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. .............................................................................................................. 53 Figura 36: P-04’; Sección meridional del perfil P-04 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ............... 54 Figura 37: P-04”; Sección septentrional del perfil P-04 (distancia de 40 - 81.6 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ........................................................................................................................................... 55 Figura 38: P-06; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. .............................................................................................................. 56 Figura 39: P-06’; Sección meridional del perfil P-06 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ............... 57 Figura 40: P-06”; Sección septentrional del perfil P-06 (distancia de 40 - 80.2 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ........................................................................................................................................... 58 Figura 41: P-08; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. .............................................................................................................. 60 Figura 42: P-08’; Sección meridional del perfil P-08 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ............... 60 Figura 43: P-08”; Sección septentrional del perfil P-08 (distancia de 40 - 80.6 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. ........................................................................................................................................... 61 INDICE DE TABLAS Tabla 1: Simbología y valores de susceptibilidad magnética de cada material presente la zona de estudio. ................................................................................................................. 45 1 RESUMEN Palabras clave: Eje Neovolcánico, Sierra Madre, Bloque Jalisco, magnetometría, anomalías magnéticas, modelo magnético. Este trabajo presenta un análisis de secciones en 2D de anomalías aeromagnéticas obtenidas directamente de la Carta Magnética de Campo Total Tepic F13-8. El área de estudio abarca la estructura conocida como Rift Tepic-Zacoalco, que corre desde el extremo oeste de Guadalajara, hasta el Pacifico, en la región de Nayarit. El estudio pretende comprender y caracterizar las estructuras subsuperficiales a lo largo de la unión entre tres grandes rasgos fisiográficos de la República Mexicana, esto es, La Sierra Madre Occidental, el Bloque Jalisco y en medio de ellos el Eje Neovolcánico. Una mejor comprensión de los rasgos tectónicos y estructurales de esta región puede ayudar a posteriores estudios sísmicos y geológicos. Para el estudio se realizaron ocho perfiles paralelos con dirección casi Norte - Sur en la zona de estudio, con una longitud aproximada de 80 kilómetros cada uno. La obtención de los datos a interpretar se ha hecho mediante el trazado de los perfiles directamente sobre la Carta Magnética mencionada anteriormente, cortando de forma normal las isolíneas de anomalías. Se desarrollara un algoritmo para prolongar la altitud de las anomalías (que son tomadas a 300 metros sobre el terreno), con la intención de filtrar anomalías con un alto número de onda. Finalmente, los perfiles prolongados se interpretaran de forma directa con un software comercial, que aplica un algoritmo clásico sumamente vigente. El análisis por el método geofísico aplicado en este trabajo, es decir, el aeromagnetismo, demostró ser útil para lograr los objetivos propuestos, con la investigación se obtiene una interpretación geológica de la zona, mediante los modelos de los diferentes perfiles y sentando así las bases para futuras investigaciones. 2 1.- INTRODUCCIÓN “La Geofísica es la ciencia cuyo objeto es el estudio científico de la Tierra y su aplicación a la búsqueda de recursos naturales, reducción de efectos causados por desastres naturales y la preservación del medio ambiente”, (fuente, Ministerio de Fomento del gobierno Español). [1] Los métodos Geofísicos son en general la forma más rápida, económica y de menor impacto ambiental, que debe emplearse para obtener información detallada del subsuelo; profundidad, extensión, naturaleza, estructura de las formaciones, etc. La magnetometría es una disciplina geofísica muy antigua y sus inicios se remontan al año 1600 con la publicación del libro "De Magnete" por el inglés William Gilbert; dicho libro es una compilación de los conocimientos ya existentes sobre magnetismo, además en el mismo se estableció el uso del concepto campo geomagnético general. Pero no fue hasta la invención del primer magnetómetro por Carl Friedrich Gauss en 1833 que fue posible medir la intensidad magnética absoluta. Los primeros estudios de magnetometría se realizan con exploración aérea (aeromagnetometría), posteriormente en áreas específicas, se realiza la exploración magnética terrestre, enfocada a la localización puntual de yacimientos minerales ferríferos ocultos o cualquier yacimiento metálico, que contenga minerales magnéticos (como magnetita, pirrotita, etc.). [2] La magnetometría es un método geofísico relativamente simple; el campo magnético de la Tierra es perturbado por afloramientos rocosos o yacimientos minerales que contienen magnetita, lo que se refleja en anomalías magnéticas susceptibles de ser medidas en la superficie terrestre con un magnetómetro. La prospección magnética es una técnica basada en la medida y estudio de las variaciones del campo magnético terrestre, obteniéndose medidas del valor total del campo magnético o bien, opcionalmente, del gradiente de dicho campo magnético. Estas variaciones son debidas a la presencia de cuerpos susceptibles de ser magnetizados y 3 que, por tal motivo, contribuyen a modificar el campo magnético terrestre en su entorno. [3] Las anomalías magnéticas son cualquier desviación de la intensidad o la inclinación magnética del valor normal. Estas pueden deberse a la actividad solar, la presencia de objetos magnéticos o debido a estructuras geológicas formadas por rocas con minerales ferromagnéticos, siendo este último la causa más común de la existencia de variaciones espaciales del campo magnético. 1.1 OBJETIVOS Establecer una correlación entre los datos topográficos y las anomalías aeromagnéticas en la zona de estudio, realizando una interpretación cualitativa de los resultados. Modelar por medio del programa comercial GM-SYS varios perfiles con una dirección casi Norte – Sur. Generar modelos geológicos 2-D y realizar una interpretación cualitativa de los mismos. Utilizar los modelos geológicos propuestos para probar la existencia de cuerpos magnéticos bajo la superficie, lo cual es de vital importancia para posteriores estudios sísmicos y geológicos de la zona. 4 2. - ZONA DE ESTUDIO La zona de estudio está comprendida por la unión entre tres grandes rasgos fisiográficos del occidente de la República Mexicana, esto es, La Sierra Madre Occidental, el Bloque Jalisco y en medio de estos dos el Eje Neovolcánico. La Figura 1 señala la ubicación de la zona de estudio en un mapa satelital, la cual incluye la zona central del estado de Nayarit y la zona limítrofe occidental del estado de Jalisco. Además de ser un área de confluencia de tres grandes rasgos fisiográficos del occidente de México, el área es también tectónica y estructuralmente importante por incluir al Rift Tepic-Zacoalco, que se describirá en líneas posteriores. Figura 1: Ubicación de la Zona de estudio. El área de estudio se localiza en la zona central del estado de Nayarit, comprendiendo los paralelos 21o y 22o de latitud Norte y 106o y 104o de longitud Oeste respecto al meridiano de Greenwich (Figura 2). 5 Figura 2: Mapa divisional donde se indica el área de estudio (Obtenida del Servicio Geológico Mexicano [5]). 2.1 GEOLOGÍA DEL ÁREA DE ESTUDIO Debido a que este trabajo busca caracterizar la zona mencionada, es decir, generar información geológica de gran relevancia, por ejemplo, sobre las estructuras y los cuerpos geológicos que se encuentran en la sub-superficie de la zona, para ello es necesario conocer la geología presente en la superficie, razón por la cual el estudio se apoyará en la información que proporciona una Carta Geológico-Minera de la zona. La Figura 3, corresponde a la Carta Geológico-Minera Tepic F13-8 [4], que incluye al área de estudio, y muestra los diferentes tipos de rocas emplazadas en la región. 6 Figura 3: Carta Geológica-Minera Tepic F13-8 (Obtenida del Servicio Geológico Mexicano [4]), escala 1:250,000. El trazo de líneas corresponde al sistema coordenado UTM. La simbología sobre la geología es detallada en la Figura 4. La explicación de la Figura 3; donde se observa la Geología de la zona de interés. Se describe a detalle en la Figura 4; mostrando Simbología por edades y Columna Geológica de las rocas presentes en la zona de estudio. Para más información sobre los Símbolos mineros, Depósitos minerales, Símbolos topográficos, así como un breve resumen, se recomienda ver la Carta Geológica-Minera [4]. 7 Figura 4: Explicación de la Carta Geológica-Minera Tepic F13-8 (Obtenida del Servicio Geológico Mexicano [4]) donde se muestra la Simbología y la Columna Geológica de las rocas presentes en la zona de estudio. 8 Fisiográficamente se ubica dentro de la Provincia de la Sierra Madre Occidental; dentro de la Provincia Llanura Costera del Pacifico; la porción centro-sur se encuentra dentro de la Provincia del Eje Neovolcánico y el extremo sur de la carta dentro de la Provincia Sierra Madre del Sur. [6] Las rocas más antiguas son de carácter volcanoclástico y consisten de intercalación de lavas y tobas andesiticas, areniscas y lutitas aparentemente coronadas por un miembro calcáreo. La secuencia descrita puede correlacionarse con rocas similares al sur de la carta, cuya edad se estima del Jurásico superior-Cretásico inferior. Al sur de Compostela, aflora la secuencia calcárea. Intrusionando a la Secuencia Volcánico sedimentaria y cubierta discordantemente por las unidades volcánicas terciarias, aflora una roca que varía de granito a granodiorita. Por posición estratigráfica y correlación se considera de edad Cretácico superior. Este intrusivo aflora en la región sur de Tepic. En la zona del graben de Tepic-Zacoalco se tienen derrames de basaltos e ignimbritas de composición Riolítica datadas del Mioceno superior-Plioceno inferior. [4] 2.2 PROVINCIAS FISIOGRÁFICAS Las principales provincias fisiográficas conformadas por procesos geológicos ocurridos en la era Cenozoica se pueden observar en la Figura 5. Las provincias son: Sierra de La Giganta, Golfo de California, Llanura Costera de Sonora y Sinaloa, Sierra Madre Occidental, Sierras y Cuencas de Chihuahua, Sierras y Cuencas de Coahuila, Mesa Central, Sierra Madre Oriental, Faja Volcánica Mexicana o Eje Neovolcánico, Llanura Costera del Golfo, Sierra Madre del Sur, Sierra Juárez, Sierra de Chiapas y Península de Yucatán. Entre las más extensas se encuentran las provincias Golfo de California, Sierra Madre Occidental, Sierra Madre Oriental, Sierras y Cuencas de Coahuila y Chihuahua, Mesa Central, Faja Volcánica Transmexicana, Llanura Costera del Golfo de México y Península de Yucatán, Sierra de La Giganta, Sierra de Juárez y Sierra Madre de Chiapas. Los procesos que definieron la evolución y la configuración de estas provincias varían desde la apertura incipiente de cuencas oceánicas (Golfo y Península de California), hasta 9 apilamientos volcánicos en zonas de magmatismo de arco (Faja Volcánica Transmexicana). Otros procesos incluyen la deformación por compresión, el acortamiento horizontal y el fallamiento lateral en la Sierra Madre de Chiapas y la emersión de planicies que previamente pertenecían a plataformas marinas como las franjas litorales de la llanura costera del Golfo de México y la plataforma de Yucatán. [7] Figura 5: Principales provincias fisiográficas de México que experimentaron cambios significativos en su relieve o configuración durante el Cenozoico, el recuadro rojo representa la zona de estudio (Modificado de [7]). 2.3 ESTRUCTURAS TECTÓNICAS ASOCIADAS A LA ZONA DE ESTUDIO La zona de estudio se caracteriza por la interacción del Rift Tepic-Zacoalco, la Sierra Madre Occidental, el Bloque Jalisco y el Eje Neovolcánico, lo cual se puede observar en la Figura 6 la cual muestra la configuración tectónica del centro-Oeste de México y los rifts que delimitan el Bloque Jalisco. Dado que cada una de estas estructuras tectónicas presenta características diversas, a continuación se abordan cada una de ellas. 10 Figura 6: Configuración tectónica del centro-Oeste de México. Se observan las zonas de los Rift Tepic-Zacoalco, Chapala y Colima. El recuadro rojo representa la zona de estudio (Modificado de [8]). 2.3.1 SIERRA MADRE OCCIDENTAL La Sierra Madre Occidental junto al Eje Neovolcánico forma parte de las Principales provincias fisiográficas de México (Figura 5). La Sierra Madre Occidental es el sistema montañoso más espacioso del territorio nacional, abarcando todo el oeste mexicano y el extremo suroccidental de los Estados Unidos, siendo una continuación de las Montañas Rocallosas en Canadá y los Estados Unidos. Cubre una extensión de 289.000 km2, lo que representa la sexta parte del territorio de México. La Sierra Madre Occidental se extiende en dirección noroeste a Sureste casi en forma paralela a las costas del océano Pacífico, inicia en el límite internacional con el estado de Arizona, E.U.A., y termina aproximadamente en el río Santiago, a la altura del 11 estado de Nayarit, en donde se conecta con la Sierra Volcánica Transversal o Eje Neovolcánico. [9] Entre las características Fisiográficas de la Sierra Madre Occidental, está la separación del golfo de California por una amplia llanura costera que se ensancha hacia el desierto de Sonora. Sus montañas presentan una altitud media de 2,440 m y una máxima que alcanza los 3,500 m, constituyen el borde occidental de la árida altiplanicie Mexicana. Originalmente era una gran meseta, pero millones de años de erosión crearon un paisaje con picos, mesetas, grandes cañones y barrancas. Su altura media es de 2 250 m. por su continuidad, sirve de barrera occidental a la Mesa del Centro. Siendo una región relativamente alta. Compuesta en su mayor parte por una compleja secuencia de rocas ígneas intrusivas del cenozoico cubiertas en algunos lugares por delgados depósitos de suelos residuales. [9, 10] 2.3.2 EJE NEOVOLCANICO De manera similar a la Sierra Madre Occidental el Eje Neovolcánico forma parte de las principales provincias fisiográficas de México (Figura 5). También se le conoce como Sierra Volcánica Transversal. Eje Volcánico Transversal, Cordillera Neovolcánica ó Faja Volcánica Transmexicana, además suele dividirse en tres sectores: occidental, central y oriental. El Eje Neovolcánico junto con la Sierra Madre del Sur es una de las provincias con mayor variación de relieve y de tipos de rocas. Se extiende desde el Océano Pacifico hasta el Golfo de México, constituyendo una ancha faja de 130 km. Inicia en la Costa Occidental en la desembocadura del río Grande Santiago a la Bahía de Banderas, continua hacia el sureste hasta encontrar el volcán de Colima para después continuar aproximadamente sobre el paralelo 19º N, hasta llegar al pico de Orizaba, alcanzando 880 km. de longitud. Esta estructura determina el límite físico entre el Norte del continente y Centroamérica. [9] 12 Es un área de montañas altas compuesta por una secuencia compleja de rocas volcánicas del cenozoico, con fosas profundas rellenadas con material lacustre del Terciario y Cuaternario. Siendo está la cordillera más alta del país, con algunas cimas coronadas de nieve permanentemente. [9,10] El Eje Neovolcánico se desarrolla sobre el margen sudoccidental de la Placa de Norteamérica como resultado de la subducción de las Placas Rivera y Cocos a lo largo de la trinchera de Acapulco. [11] Los datos gravimétricos y modelos propuestos en [12] sugieren que el espesor cortical de la placa superior es máximo en la parte oriente del Eje Neovolcánico alcanzando los 50 km, mientras que hacia el occidente el espesor disminuye hasta los 28 km. 2.3.3 BLOQUE JALISCO El Bloque Jalisco (Figura 6) es la estructura bajo la que subduce la placa Rivera, está delimitado por la Trinchera Mesoamericana hacia occidente, la zona del Rift Colima en la parte meridional y la zona del Rift Tepic-Zacoalco en la parte septentrional oriental [13], de manera que representa lo que se conoce como un bloque tectónico, o micro placa. Los Rifts presentan orientaciones que van desde NW-SE a NE-SW y son asociados con rocas volcánicas recientes del Eje Neovolcánico, que son delimitadas por lineamientos regionales con orientaciones WNW (Figura 6) [14]. El Bloque Jalisco se mueve de manera independiente con respecto a las placas Rivera y Norteamericana, a través del Rift Tepic-Zacoalco y el Rift Colima y a lo largo de una zona de subducción en su límite costero con la placa Rivera. La evolución del bloque Jalisco parece estar relacionado con la geometría y dinámica de la placa Rivera. [15] Los Rift Tepic-Zacoalco y de Colima se unen con el Rift de Chapala, en el límite NE del Bloque Jalisco, formando una unión tipo Rift-Rift-Rift, lo que es esencialmente un punto triple continental, aproximadamente a 50 km al suroeste de la ciudad de Guadalajara, 13 formado por la unión del Bloque Jalisco, el Bloque de Michoacán y la placa Norteamericana. [13, 15, 16] 2.3.4 RIFT TEPIC-ZACOALCO El Rift Tepic-Zacoalco (Figura 6) tiene una extensión aproximada de 250 km y aproximadamente 50 km de ancho. Es una región amplia, topográficamente baja que consta de varios grabenes y semigrabenes desarrollados en distintos episodios a partir del Mioceno Tardio [14,18]. Dichos grabenes y semigrabenes están delimitados por fallas escalonadas NNW que se extienden desde el punto triple hasta la costa Oeste de México. [17] En el área de los tres Rifts son dominantes las andesitas, rocas calcoalcalinas y rocas relacionadas, típicas de los arcos volcánicos de subducción. [13] La parte norte del Rift Tepic-Zacoalco parece inactiva, mientras que la parte sur muestra evidencia geológica e indicadores geomórficos de fallas recientes [18]. Precisamente es en la parte norte de este Rift donde se ha realizado este estudio. 14 3.- ANTECEDENTES 3.1 FUNDAMENTO FÍSICO Como sabemos de diversas fuentes ([19], [20]) que la Fuerza eléctrica entre dos cargas puntuales q1 y q2 separadas por una distancia r se puede calcular gracias a la Ley de Coulomb: 𝐹 = 𝑞1𝑞2 4𝜋𝜀𝑟2 donde 𝜀 es la permitividad. De manera similar se puede calcular la fuerza magnética (Fm) entre dos polos magnéticos de carga m1 y m2 también dada por la ley de Coulomb: 𝐹𝑚 = 𝑚1𝑚2 4𝜋𝜇𝑟2 donde r es la distancia que separa m1 y m2, mientras que µ es la permeabilidad magnética. La fuerza es atractiva en polos de signo opuesto y repulsiva entre polos del mismo signo. La magnetización (M) se define como la influencia que ejerce un campo magnético sobre un material. En un campo magnético como el de la Tierra la magnetización se define por: 𝑀 = 𝑘𝐻 donde H es la Intensidad del Campo Magnético, mientras que k es la susceptibilidad magnética, la cual se explica más adelante. La magnetización ocurre al colocar un cuerpo en un campo magnético y este adquiere una magnetización inducida o polarización, y como se puede observar en la ecuación anterior, su valor dependerá tanto del campo como de la susceptibilidad del cuerpo. La Intensidad de Campo Magnético es una cantidad de campo fundamental, descrita como: 𝐻 = 𝐵 𝜇0 − 𝑀 donde B es la densidad de flujo magnético, 𝜇𝑜 es la permeabilidad magnética en el vacío. De la ecuación anterior es posible llegar a una relación más sencilla, para ello se requiere despejar B: 𝐵 = (𝐻 + 𝑀)𝜇0 15 Sustituir la Magnetización: 𝐵 = (𝐻 + 𝑘𝐻)𝜇0 Despejar la Intensidad del Campo Magnético: 𝐵 = (1 + 𝑘)𝜇0𝐻 Donde la permeabilidad magnética relativa se puede definir como 𝜇𝑅 = 1 + 𝑘 = 𝜇 𝜇𝑜 . Entonces: 𝐵 = 𝜇𝑅 𝜇0 𝐻 Dado que la permeabilidad magnética 𝜇 es el producto de la permeabilidad magnética en el vacío y la permeabilidad magnética relativa, finalmente obtenemos: 𝐵 = 𝜇 𝐻 De esta ecuación podemos observar como los campos magnéticos se relacionan directamente; por ejemplo, en el aire y agua la permeabilidad magnética es prácticamente 1, entonces es posible considerar B = H. 3.2 EL CAMPO GEOMAGNÉTICO El campo geomagnético lambien llamado campo magnético terrestre, es el campo magnético que se extiende desde el núcleo interno de la Tierra hasta el límite en el que se encuentra con el viento solar, debido a que la Tierra se comporta como un enorme imán. La Tierra posee un poderoso campo magnético, como si el planeta fuera una enorme imán en forma de barra, cuyo polo sur se considera cerca del polo norte geográfico y polo norte cerca del polo sur geográfico. Aunque los polos magnéticos terrestres reciben el nombre de polo norte magnético y polo sur magnético, su magnetismo real es el opuesto al que indican sus nombres. Precisamente por este comportamiento de imán, es que presenta líneas de campo o líneas flujo que salen del polo sur magnético y entran por el polo norte magnético. Las posiciones de los polos magnéticos no son constantes y muestran notables cambios de año en año. El campo magnético de la Tierra tiene tendencia a trasladarse hacia el oeste a razón de 19 a 24 km por año. [21] 16 Como ya se mencionó en primera aproximación, el campo de la Tierra puede ser considerado como el producido por un imán de barra localizado cerca del centro de la Tierra, y cuyo momento magnético es de 8.0×1022 Am2. El campo en la superficie tiene una magnitud que va desde unos 30000nT cerca del ecuador hasta unos 60000nT cerca de los polos. [22] El eje a lo largo del imán se conoce como eje magnético, que extendido hasta la superficie de la Tierra emerge en los puntos llamados polos magnéticos. El eje magnético está inclinado aproximadamente 11.5o con respecto al eje de rotación de la Tierra (que a su vez forma un ángulo de 23.5o con la normal al plano de la órbita de la Tierra con respecto al Sol), de manera que el polo norte magnético se encuentra ubicado sobre la Tierra a 79.5o de latitud norte y 71.6o de longitud oeste, mientras que el polo sur magnético está a 79.5o de latitud sur y 108.4o de longitud este. En la Figura 7 se muestra la forma de este campo magnético. El eje magnético pasa por un punto que se encuentra directamente por debajo del Océano Pacífico medio, a unos 340km del centro de la Tierra. [23] Figura 7: Campo geomagnético. En primera aproximación este puede ser considerado como el producido por un dipolo magnético. 1: Plano de la órbita de la Tierra, 2: Eje de rotación, 3: Eje magnético. (Obtenida de [23]). 17 Actualmente se sabe que aproximadamente el 95% del campo magnético de la Tierra es producido en su interior. De forma que, dentro de este existe un mecanismo que es capaz de generar y mantener el campo de un simple dipolo. La explicación está en las propiedades y movimientos del núcleo de la Tierra, un cuerpo esférico de material metálico, probablemente de una composición níquel-hierro, cuya parte externa se encuentra en estado líquido. El diámetro del núcleo es de alrededor de 6920 km, poco mayor que el radio terrestre. El núcleo se encuentra rodeado por el manto, consistente de roca semisólida y de un espesor de 2860 km. El magnetismo interno de la Tierra se explica por la llamada teoría del dinamo, según la cual el núcleo externo (de consistencia líquida) está dando vueltas lentamente, con respecto al manto semisólido, generando de esta manera corrientes eléctricas que rodean al núcleo. Estas corrientes generan a su vez un campo magnético, parte del cual escapa a la superficie de la Tierra, dándonos el campo magnético que observamos, y otra parte interacciona con el núcleo líquido que está en movimiento, sosteniendo de esta manera la acción del dínamo. [23, 24] La intensidad del campo geomagnético se representa de dos maneras [23, 25]: 1. Componentes XYZ: formado por tres componentes ortogonales entre sí, siendo X la componente que apunta hacia el Norte, mientras que Y es la componente que apunta al Este y la componente Z es la que apunta a el centro de la Tierra. 2. Componentes HDI: H es la magnitud de la componente horizontal que es considerada positiva cualquiera sea su dirección, mientras que D representa la declinación magnética, es decir, el ángulo desde X hacia H, finalmente I es el ángulo formado entre la intensidad total y la componente horizontal. En la Figura 8 es posible apreciar de manera gráfica la relación existente entre las componentes antes mencionas. Por simple geómetra es posible calcular las siguientes relaciones entre variables: X =H Cos D Y =H Sen D F2 = X2 +Y2 +Z2 Tan I = Z / H 18 Figura 8: Componente del campo geomagnético. (Obtenida de [23]). El campo geomagnético sufre variaciones en tiempo y espacio debido a distintas causas por elementos internos y externos. Una prueba visible del campo magnético son las auroras, pues la interacción del campo magnético de la tierra con las partículas del viento solar crea las condiciones de generación de los fenómenos conocidos como auroras cerca de los polos. 3.3 MAGNETOMETRIA El método magnético es el método geofísico de prospección más antiguo y sus aplicaciones principales son variadas, pero entre las principales se encuentran: Exploración Geológica. Exploración Minera. Medir las propiedades de las rocas. Estudios de contaminación de suelos. Búsqueda y localización de objetos metálicos sepultados. Detección de minas, granadas, etc. Arqueología. Búsqueda de agua subterránea. 19 Método de campos potenciales económico. Interpretar el tipo de fallas que afectan el basamento magnético. Las técnicas magnéticas miden el campo magnético remanente asociado con cambios en el campo magnético de la Tierra asociada a la estructura geológica. [26] El estudio de los métodos magnéticos aplicados a problemas geológicos ha ido avanzando en paralelo con el desarrollo de los magnetómetros. En la industria petrolera se ha utilizado como método de reconocimiento para estudios regionales desde principios del siglo pasado y con el desarrollo del magnetómetro aéreo en 1941 se facilitaron las mediciones magnéticas a lugares de difícil acceso por tierra; además, se empezó a cubrir mayor superficie en menor tiempo, lo que representa una gran ventaja de esta técnica exploratoria. [26, 27] El principio de la magnetometría se describe en la Figura 9. Un cuerpo mineralizado genera anomalías magnéticas al cambiar los valores que se esperan para el Campo magnético de la tierra, es decir, al existir un cuerpo mineralizado y colocar un magnetómetro sobre el mediríamos el campo magnético de la tierra y también del cuerpo mineralizado. De no existir dicho cuerpo solo mediríamos el campo magnético de la tierra. Figura 9: Principio de la magnetometría (Tomado de [28]). 20 En la prospección petrolífera el método magnético entrega informaciones acerca de la profundidad de las rocas pertenecientes al basamento. A partir de estos conocimientos se puede localizar y definir la extensión de las cuencas sedimentarias ubicadas encima del basamento, que posiblemente contienen reservas de hidrocarburos. [28] El método se emplea en la delineación de depósitos magnéticos intrasedimentarios como rocas subvolcánicas e intrusiones emplazadas a poca profundidad, que cortan la secuencia sedimentaria normal. Como las rocas sedimentarias generalmente ejercen un efecto magnético despreciable en comparación con el efecto magnético generado por las rocas ígneas. La mayoría de las variaciones de la intensidad magnética medidas en la superficie terrestre resulta de cambios litológicos o topográficos asociados con rocas ígneas o con rocas del basamento. El desarrollo reciente de magnetómetros de alta precisión posibilita ahora la definición de pequeñas repuestas magnéticas de alta frecuencia y la detección de variaciones muy pequeñas de la intensidad magnética, que podrían ser relacionadas con variaciones diminutas en el carácter magnético de rocas sedimentarias someras. [28] Las exploraciones mineras también aplican el método magnético en la búsqueda directa de minerales magnéticos y minerales no magnéticos, pues estos ejercen un efecto magnético que es medible en la superficie terrestre. También gracias a la magnetometría es posible localizar zonas de fallas, de cizallamiento y de fracturas, las cuales es común que contengan una gran cantidad de minerales. El conocimiento de los sistemas de fallas y fracturas es de vital importancia para la explotación de materiales subterráneos. Por ejemplo Héctor López-Loera [26] utiliza este método para la prospección de agua subterránea. 3.4 ANOMALIAS MAGNETICAS Como se ha mencionado previamente, en algunas zonas de la corteza terrestre se producen grandes desviaciones del valor teórico del campo magnético principal. Estas 21 desviaciones se conocen en general como anomalías magnéticas. En buena parte son debidas a la magnetización remanente de rocas volcánicas. [29] Las anomalías que abarcan enormes áreas se llaman regionales, mientras que las que ocupan áreas de decenas de metros a unos pocos kilómetros cuadrados se llaman anomalías locales. De manera que grandes volúmenes de materiales magnéticos cambiarán la intensidad del campo terrestre. Los Teslas son muy grandes para expresar las anomalías causadas por la geología, así que comúnmente se utiliza el nano tesla (nT = 10-9T), el sistema CGS utiliza el gamma (𝛾) pero este es equivalente al nT. 3.5 FILTRADO DE ANOMALIAS MAGNETICA Las anomalías magnéticas contienen la suma de todos los efectos magnéticos de distintas fuentes, como cuerpos superficiales, profundos, intermedios, ruido, entre otros. Razón por la cual es de gran utilidad realizar algún proceso de filtrado a los datos, eliminando así el ruido presente en los datos magnéticos, separando los efectos magnéticos profundos y superficiales, enfatizando las anomalías de interés y facilitando la interpretación de los datos que se desean. Existen diferentes procesos de filtrado para distintos usos, a continuación se explican brevemente algunos de ellos: 3.5.1 CONTINUACIÓN ASCENDENTE La continuación ascendente transforma el campo potencial medido en una superficie en el campo que se mediría en otra superficie hacia arriba de la fuente. Esta transformación atenúa las anomalías con respecto a la longitud de onda; y cuanto más corta es la longitud de onda, mayor es la atenuación. 22 En este sentido, el proceso de continuación ascendente degrada los datos medidos, y podríamos preguntarnos por qué ese proceso tendría alguna aplicación, existen dos ejemplos útiles. Primero, a veces es necesario comparar o fusionar levantamientos aéreos medidos a altitudes diferentes, y la continuación ascendente proporciona una manera de transformar los levantamientos individuales en una superficie consistente. En segundo lugar, la continuación ascendente tiende a acentuar las anomalías causadas por fuentes profundas a expensas de las anomalías causadas por fuentes poco profundas. Un levantamiento magnético sobre terreno volcánico joven, por ejemplo, puede estar dominado por anomalías de longitud de onda corta debido a rocas volcánicas cercanas a la superficie; La continuación hacia arriba se puede utilizar para atenuar las anomalías de fuentes poco profundas con el fin de enfatizar fuentes más profundas. [27] Para una Continuación Ascendente se puede considerar la siguiente ecuación: 𝐹 ↑= ⅇ−𝛥𝑧|𝑘| donde 𝐹 ↑ es la exponencial requerida para realizar una Continuación Ascendente, ∆z es el cambio que se desea hacer en la elevación, es decir, el valor donde se desea evaluar los datos, mientras que k es el número de onda. El proceso de continuación hacia arriba atenúa todos los números de onda excepto |𝑘| = 0, y el grado de atenuación aumenta con el aumento de z. Para lograr una continuación ascendente es necesario: aplicar la transformada de Fourier a los datos magnéticos medidos previamente, después multiplicar la trasformada por el término exponencial 𝐹 ↑, y finalmente aplicar la transformada de Fourier inversa a este producto. De esta manera las longitudes de onda más cortas de las anomalías originales son esencialmente eliminadas, mientras que las anomalías fundamentales se mantienen en una forma más suave. 3.5.2 CONTINUACIÓN DESCENDENTE La Continuación Descendente es muy similar a la anterior, pero ésta se realiza en regiones más cercanas a la fuente, bajando el punto de observación. Cabe mencionar que en este caso se realzan las anomalías más superficiales. Este proceso se utiliza 23 comúnmente para el estudio de anomalías específicas, o también como herramienta para unificar levantamientos magnéticos realizados a distintas altitudes. [27] Para una Continuación Descendente se puede emplear la siguiente ecuación: 𝐹 ↓= ⅇ𝛥𝑧|𝑘| donde 𝐹 ↓ es la exponencial requerida para realizar la Continuación Descendente, ∆𝑧 es el cambio que se desea hacer en la elevación, mientras que k es el número de onda. Realizar una Continuación Descendente es muy similar a una Continuación Ascendente, primeramente se aplica la transformada de Fourier a los datos magnéticos, después se multiplica la trasformada por el término exponencial 𝐹 ↓, y finalmente se aplica la transformada de Fourier inversa a este producto. 3.5.3 REDUCCIÓN AL POLO Las anomalías de gravedad positiva tienden a ubicarse sobre concentraciones de masa de forma vertical, Lo que no necesariamente ocurre con las anomalías magnéticas, pues la magnetización y el campo ambiental no están dirigidos verticalmente. Una distribución uniforme del material magnético puede desplazar la anomalía lateralmente, distorsionar su forma e incluso cambiar su signo si el campo magnético inductor no es vertical. [27] La Reducción al Polo es un procedimiento que recalcula los datos de Intensidad de Campo Magnético Total a un conjunto de datos que tiene una inclinación de 90°, es decir, transforma las anomalías medidas con cualquier dirección a la que sería medida bajo un campo vertical. De esta manera las anomalías asimétricas debidas a la inclinación se eliminan y se reubica las anomalías sobre los cuerpos causantes, de modo que los cuerpos simétricos producirán anomalías simétricas. 3.6 SUCEPTIBILIDAD MAGNETICA Por definición, el valor de la susceptibilidad magnética proporciona una medida del grado en que una sustancia puede ser magnetizada. Esta ocurre cuando un cuerpo está situado 24 en un campo magnético adquiriendo una magnetización, es decir, nos permite describir la forma y grado de la magnetización de cualquier sustancia al ser colocada en un campo magnético. [31, 32] La susceptibilidad magnética es una propiedad inherente de la materia y se define dada la siguiente relación: kH = M Donde k es la susceptibilidad magnética, H es el campo y M la magnetización. Despejando k podemos observar que la susceptibilidad es adimensional. Los materiales pueden clasificarse magnéticamente, según el valor de su susceptibilidad magnética, en tres grandes grupos: Diamagnéticos (k<0), Paramagnéticos (k>0) y Ferromagnéticos (k>>0). Ferromagnéticos: se refiere a materiales fuertemente atraídos por un imán, como una pieza de hierro. En particular, minerales ferromagnéticos conocidos incluyen magnetita, maghemita, pirrotita y pentlandita. [33] Paramagnéticos: se refiere a materiales menos magnéticos que ferromagnéticos, y separables en separadores magnéticos con corrientes de campo de hasta 1.70 amperios. [33] Diamagnéticos: se refiere a materiales que no son magnéticos a 1.70 amperios, y serían repelidos por electroimanes con corrientes más fuertes. [33] La determinación de k es un procedimiento experimental que requiere cuidado, debido a la magnitud de los valores de k para la mayoría de las sustancias. El valor de k para sustancias paramagnéticas y diamagnéticas es |k| << 1, típicamente k ≈10−6 (SI), para la mayoría de estas sustancias. Del mismo modo, es importante enfatizar que a pesar de que k es una variable adimensional, su valor depende del sistema de unidades. [32] 25 4.- DESARROLLO DEL TEMA 4.1 METODOLOGÍA Con la finalidad de entender las estructuras subsuperficiales pertenecientes a los sistemas geológicos ya mencionados y cubrir los objetivos planteados, se modelarán de forma directa ocho perfiles o secciones que atraviesan al campo de anomalías aeromagnéticas en sentido casi norte - sur. Para cada perfil se obtienen manualmente los valores de anomalías, que no son equidistantes. Estos valores se han interpolado con un algoritmo de splines cúbicos implementado en Matlab, para obtener datos equidistantes y modelar en una etapa futura estructuras más regionales. Esto último requerirá elevaciones también equiespaciadas. Anomalías y elevaciones se muestran en figuras de apartados posteriores. 4.2 DATOS 4.2.1 INFORMACIÓN DE LA CARTA MAGNÉTICA DE CAMPO TOTAL Para la realización de este estudio, la información Magnética se obtuvo de la Carta Magnética de Campo Total (Tepic F13-8, escala 1:250,000 [5]) con base cartográfica INEGI, realizada por el Servicio Geológico Mexicano (SGM) con datos obtenidos en una fecha de vuelo promedio del 27 de Marzo de 1984, con una altura de vuelo sobre el nivel de terreno de 300m y distancia entre líneas de vuelo de 1,000 m. Los parámetros del campo Geomagnético al centro de la carta son: Altura: 1.8 Km Intensidad: 43,594 nT Inclinación: 48o 15’ Declinación: 08o 56’ Las características del equipo son: Avión: ISLANDER BN2-B27 Magnetómetro: GEOMETRICS G-803 Sensibilidad: 0.25 nT 26 Sistema de adquisición de datos: GEOMETRICS G-704 Cámara: FOTOGRÁFICA 35mm, AUTOMAX G-2 Magnetómetro estación base: GEOMETRICS G-826A Sensibilidad: 1 nT Radar Altímetro: SPERRY Sistema de navegación: VISUAL 4.2.2 INFORMACIÓN DE LA ELEVACIÓN En este estudio se obtuvieron los datos de elevación colocando los perfiles deseados en el programa de Google llamado Google Earth Pro, guardándolos de tal manera que puedan ser analizados utilizando el SOFTWARE de Garmin (GPS) llamado MapSource, que está disponible para descargar de forma gratuita. El MapSource tiene una gran precisión y es capaz de proporcionar datos de elevación y posición de rutas o puntos específicos. Por estas razones se utilizó para este trabajo. 4.3 PERFILES Los perfiles que se interpretan en este trabajo son ocho, y se han elegido siguiendo una convención común: deben ser normales a las isolíneas de rasgos anómalos de tendencia bidimensional. Cada perfil tiene una distancia aproximada de ochenta kilómetros de largo, además estos son paralelos unos de los otros, presentando una dirección SSO-NNE y siendo nombrados P-01, P-02, P-03, P-04, P-05, P-06, P-07 y P-08 de oeste a este, como se puede apreciar en la Figuras 10 y 12. En la Figura 10 se pueden observar los distintos perfiles seleccionados, esbozados sobre la Carta mencionada. Los perfiles intersecan las principales anomalías magnéticas (debido a su gran intensidad o abundancia en la zona). 27 Figura 10: Mapa de Campo Magnético Total (Modificado del Servicio Geológico Mexicano [5]) donde se muestra la zona de interés; así como los perfiles seleccionados. Figura 11: Valores de Intensidad de Campo Magnético, en nano Teslas, de la Figura 10 (Modificado del Servicio Geológico Mexicano [5]). Para interpretar de manera adecuada la Figura 10, se requiere el uso de la escala de valores de Intensidad de Campo magnético localizada en la Figura 11, la cual se encuentra en nano Teslas. De manera similar en la Figura 12 se pueden observar los ocho perfiles seleccionados sobrepuestos en esta ocasión sobre la Carta Geológica-Minera ya mencionada. La explicación de la Figura 12 se halla en la Figura 3, donde se muestra la Simbología y una Columna Geológica de las rocas de la zona. 28 Figura 12: Mapa Geológico-Minero (modificado del Servicio Geológico Mexicano [4]) donde se muestra la zona de interés; así como los perfiles seleccionados y la geología presente. 4.5 INTERPRETACIÓN DIRECTA Para preparar los perfiles que se habrán de interpretar, se interpolaron los valores magnéticos y de elevación para una misma posición. La interpolación se realizó mediante un programa desarrollado en Matlab, empleándose la función SPLINE. Esta función interpola aplicando el método de spline cúbico (Véase [34]). Dicho método ajusta un polinomio de tercer grado en cada intervalo de datos, brindando los valores intermedios entre los puntos ó valores conocidos. La finalidad de obtener los valores previos es para correlacionar las anomalías observadas con un modelo geológico estructural propuesto. Esto se realiza a partir de la interpretación de las señales observadas en la Figura 10 y del conocimiento geológico (Figura 12); primero se propone un modelo inicial, que después será modificado mediante prueba y error; con la finalidad de reducir los errores entre la anomalía observada y la generada. Sin embargo, para una señal existe más de una solución, así que este problema de ambigüedad solo se podrá reducir mediante el conocimiento geológico de la zona. 29 Este proceso de modelación se realizó utilizando el paquete comercial GM-SYS (1999) de NGA (Northwest Geophysical Associates, Inc.), a partir de las secciones magnéticas. Dicho programa cuenta con una interfaz gráfica con diversas opciones que permiten el cálculo para la respuesta magnética y gravimétrica, manipulación de un modelo geométrico complejo, así como el error entre la señal calculada y la observada. El programa emplea el algoritmo bidimensional de Talwani (1965) [35] y los métodos descritos por Talwani y Heirtzler (1964) [36] para calcular los resultados de los modelos gravimétricos y magnéticos utilizando los algoritmos desarrollados en [37]. 30 5.- RESULTADOS Y DISCUSIÓN 5.1 CORRELACION VISUAL ENTRE ANOMALÍAS MAGNÉTICAS Y TOPOGRAFÍA Los más altos contrastes magnéticos ocurren en la interface entre la superficie terrestre y la atmósfera ([38]). La susceptibilidad magnética del aire es de 0.0 SI, mientras que las rocas de la corteza, donde se colectan los datos magnéticos, muestran susceptibilidades en el rango de 0.0 a 0.2 SI ([39]). Variaciones en la topografía superficial producen anomalías magnéticas que se relacionan a los contrastes magnéticos entre la superficie del terreno y el aire, y no a los cuerpos magnéticos sepultados que son los objetivos de los levantamientos magnéticos ([38]). En esta sección se muestran correlaciones visuales entre la topografía y las anomalías aeromagnéticas en los perfiles investigados, buscando precisamente los rasgos anómalos que pudiesen ser debidos a la topografía irregular de la zona de estudio. Las interpolaciones de los perfiles de anomalías magnéticas y elevaciones obtenidas mediante el uso de la función spline cúbico se muestran en las Figuras 13, 14, 15, 16, 17, 18, 19 y 20. Las figuras mencionadas se han elaborado con el fin de buscar, mediante correlación visual, un efecto que podría tener la topografía local en las anomalías aeromagnéticas. Una interpretación cualitativa de los perfiles magnéticos en función de los rasgos predominantes y la topografía es elaborada en esta sección. En los primeros 9 km el perfil P-01 (Figura 13) presenta una topografía compleja que varía desde 48 metros sobre el nivel del mar (m.s.m.) hasta los 509 m.s.m. Después se tiene una topografía más suave hasta el kilómetro 36 variando de 9 m.s.m. a 121 m.s.m. y posteriormente se incrementa nuevamente pero de manera más abrupta hasta alcanzar los 795 m.s.m. Esta topografía se incluye en el modelado. La anomalía magnética (Figura 13) tiene una amplitud absoluta de aproximadamente 1000 nT. Varios máximos y mínimos caracterizan al perfil. Destaca un máximo de aproximadamente 500 nT hacia los 50 kms del perfil. En este Perfil es posible observar que existe una ligera correlación directa entre el valor magnético y la elevación en varios puntos a lo largo del perfil. Esto se aprecia en los 31 primeros 9 Km. En el tramo entre los 10 y los 35 kilómetros se puede observar una ligera correlación positiva (figura 13). En este tramo las anomalías magnéticas pueden deberse a una profundización del basamento intrusivo, perteneciente al Bloque Jalisco. La correlación es más o menos clara también en el tramo entre 50 y 60 kilómetros, entre el máximo de casi 500 nT ya mencionado, con una elevación que llega casi a los 800 metros sobre el nivel del mar. La correlación en estos tramos del perfil puede deberse a la presencia de grandes y espesos cuerpos basálticos. Figura 13: Sección P-01; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. Los primeros 11 kilómetros del perfil P-02 (Figura 14) presentan una topografía compleja formada por diferentes picos alcanzando una elevación de 600.95 m.s.m. en el kilómetro 5.2, posteriormente disminuye de manera abrupta a una topografía más suave hasta el kilómetro 41 donde la topografía aumenta progresivamente debido a una montaña cuyo punto más alto se encuentra en el kilómetro 60 y alcanza los 660.94 m.s.m. 32 En este perfil la correlación no es muy clara. Una posible correlación positiva parece darse en los primeros 7 kilómetros. En el perfil magnético resalta un mínimo de casi -750 nT, que parece indicar la presencia de una gran falla. En la sección entre los 35 y 41 kilómetros destacan dos grandes montañas en el Perfil de elevación que no generan un cambio considerable en los valores magnéticos, debido a la existencia de grandes cuerpos rocosos con un valor de susceptibilidad magnética bajo comparado con los materiales cercanos a ellos. Figura 14: Sección P-02; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. El perfil P-03 (Figura 15) se caracteriza por un máximo que llega a casi 400 nT y que podría correlacionarse con la elevación también destacable en la topografía, que se encuentra aproximadamente a una distancia de 53 - 57 Km y que tiene alrededor de 2250 33 m.s.m. El máximo podría deberse a cuerpos intrusivos que afectan el magnetismo de la zona. Cabe mencionar además que el perfil de anomalías magnéticas presenta pequeños máximos los cuales no se correlacionan con el perfil de elevación. Estos máximos presentan un alto número de onda, razón por la cual para el modelado directo de los perfiles se realizará un filtrado de datos, y en este caso el más apropiado son las continuaciones analíticas ascendentes. Figura 15: Sección P-03; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. La topografía del perfil P-04 (Figura 16) es compleja variando constantemente, con un hundimiento en el kilómetro 2 (290.78 m.s.m.), a su vez cuenta con dos montañas que resaltan sobre las otras elevaciones; la montaña del kilómetro 30 llega a los 1714.9 m.s.m. y la montaña del kilómetro 78 alcanza los 1511.12 m.s.m. 34 Los perfiles de la Figura 16 muestran una ligera similitud, especialmente en el tramo entre los 20 y los 45 km. A grandes rasgos se pueden observar dos grupos de anomalías de amplitud positiva, destacando las anomalías positivas que alcanzan casi 300 nT hacia el extremo norte del perfil. Figura 16: Sección P-04; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. El Perfil P-05 (Figura 17) presenta una similitud considerable a lo largo de los 80 Kilometros, especialmente en la sección entre los 5 y 21 km, a grandes rasgos se observa dos anomalias de amplitud positiva en el centro de este rango, donde las anomalias positivas alcanzan un valor aproximado de 800 nT. La topografía del perfil P-05 es muy suave aumentando prograsivamente en el centro del perfil, destacando el abrupto levantamiendo del kilómetro 45 formando una montaña (2050 m.s.m.) en el kilómetro 51. 35 Además destaca el perfil de anomalias magneticas presenta una gran cantidad de pequeñas máximos que no se correlacionan con el perfil de elevación, dichos maximos presentan un alto numero de onda, razón por la cual se realizará un filtrado por continuación analitica ascendente para el modelado directo. Figura 17: Sección P-05; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. En los primeros 19 km del perfil P-06 se presenta una topografía suave que va de los 1016.5 m.s.m. a los 1350 m.s.m., después se tiene una topografía más irregular. Del kilómetro 40 al 50 se encuentra la zona más baja llegando hasta 252.5 m.s.m. (42 km) luego la topografía sube de manera progresiva hasta los 1791 m.s.m. En el perfil P-06 (Figura 18) destaca la correlacion inversa ubicada al centro del perfil aproximadamente del kilometro 35 al 55, el máximo magnetico en esta zona podria 36 deberse a la existencia de cuerpos intrusivos con un gran valor de susceptibilidad magnética afectando el magnetismo. La correlación es menos clara en el tramo inicial del perfil (aproximadamente 15 km), asimismo en esta zona destaca un máximo de aproximadamente 300 nT seguido de un mínimo de casi -500 nT lo cual parece indicar la presencia de una gran falla. Figura 18: Sección P-06; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. La topografía del perfil P-07 (Figura 19) presenta un comportamiento suave los primeros 30 kilómetros, aumentando progresivamente hasta la montaña del kilómetro 20 (2456.38 m.s.m.) seguido de un hundimiento abrupto en el kilómetro 29 (510.96 m.s.m.), después presenta una topografía más irregular, aumentando hasta la montaña del kilómetro 62. 37 El perfil P-07 (Figura 19) no parece presentar una similitud considerable, lo que posiblemente se deba a la existencia de materiales rocosos incrustados en la región. Una posible correlación inversa se puede observar durante los primeros 12 kilómetros. En el perfil magnético resalta un mínimo de casi -600 nT en la distancia de 40 a 50 km, seguido de un alto magnético de aproximadamente 100 nT, lo que parece indicar la presencia de una falla. En la sección aproximada al kilómetro 20 se observa una gran montaña en el Perfil de elevación que no genera un cambio considerable en los valores magnéticos, posiblemente debido a la composición superficial de la misma. Figura 19: Sección P-07; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. El perfil P-08 (Figura 20) inicia con una topografía casi constante durante los primeros 20 km, variando de 1000 a 1500 m.s.m., después baja abruptamente hasta los 378.6 nT (21.4 38 km) y vuelve a subir hasta los 1800 m.s.m. Posteriormente se mantiene una topografía suave aumentando progresivamente debido a una montaña de 2226.4 m.s.m. en su punto más alto (55 km). El perfil P-08 presenta una amplitud absoluta de aproximadamente 500 nT, con varios máximos y mínimos caracterizan al perfil. Existe una ligera correlación directa entre el valor magnético y la elevación en varios puntos a lo largo del perfil. Destacando un máximo de aproximadamente 200 nT hacia el kilómetro 17 del perfil. Figura 20: Sección P-08; a) Perfil de intensidad de campo magnético (anomalías aeromagnéticas) y b) Perfil de elevación. Distancia en ambos casos en Kilómetros. Elevación en metros sobre el nivel del mar y datos aeromagnéticos en nano Teslas. 39 5.2 CONTINUACIONES ASCENDENTES DE PERFILES AEROMAGNÉTICOS Debido a la existencia de cierta correlación entre la topografía y las anomalías aeromagnéticas descrita en el capítulo previo, considerando filtrar altos números de onda que se deben comúnmente a rasgos geológicos someros y de poca amplitud, y además para tener una interpretación de los rasgos regionales que son de interés para este estudio, se han realizado varias continuaciones ascendentes en 2D de los perfiles que se interpretan. La idea de realizar varias continuaciones analíticas es para elegir una anomalía ascendida adecuada para la interpretación. Los resultados obtenidos se muestran en las Figuras 21 a 28. Se realizaron 8 continuaciones entre los 200 metros y los 3700 metros de altura y a cada 500 metros; es decir, que las anomalías aeromagnéticas fueron prolongadas en el rango entre 500 y 4000 metros sobre el nivel del terreno. En cualquier Figura entre los numerales 21 a 28 se puede observar un comportamiento semejante; esto es, se observa una fuerte atenuación de los rasgos anómalos más allá de los 1500 metros sobre el nivel del terreno, perdiéndose totalmente los rasgos de alto número de onda a los 4000 metros sobre el terreno. Por ello, se ha elegido interpretar las anomalías que se han continuado a 1500 metros. 40 Figura 21: Continuación Ascendente para el Perfil P-01; cada 500 m en el intervalo de 200 a 3700 m. Figura 22: Continuación Ascendente para el Perfil P-02; cada 500 m en el intervalo de 200 a 3700 m. 41 Figura 23: Continuación Ascendente para el Perfil P-03; cada 500 m en el intervalo de 200 a 3700 m. Figura 24: Continuación Ascendente para el Perfil P-04; cada 500 m en el intervalo de 200 a 3700 m. 42 Figura 25: Continuación Ascendente para el Perfil P-05; cada 500 m en el intervalo de 200 a 3700 m. Figura 26: Continuación Ascendente para el Perfil P-06; cada 500 m en el intervalo de 200 a 3700 m. 43 Figura 27: Continuación Ascendente para el Perfil P-07; cada 500 m en el intervalo de 200 a 3700 m. Figura 28: Continuación Ascendente para el Perfil P-08; cada 500 m en el intervalo de 200 a 3700 m. 44 5.3 PERFILES GEOLOGICOS 2D OBTENIDOS POR INTERPRETACIÓN DIRECTA El modelado directo realizado pretende comprender y caracterizar las estructuras subsuperficiales que existen a lo largo de cada perfil. Utilizar el mapa Geológico-Minero (Figura 12) permite establecer cuáles son los materiales presentes en la superficie de cada perfil. Además, existe un problema de ambigüedad ([40]), donde una situación se puede interpretar de más de una manera (es decir, existe más de una solución para justificar la señal magnética). Dicho problema solo se puede disminuir mediante el conocimiento geológico de la zona y un manejo adecuado de los datos y propiedades físicas. Los datos aeromagnéticos fueron tratados con un filtrado de continuación ascendente a 1200m con la intención de suavizar las anomalías, reduciendo los rasgos anómalos de alto número de onda que son producidos por fuentes pequeñas y someras, y dando una mayor importancia a las estructuras y cuerpos geológicos profundos sobre las fuentes superficiales. Con esto la información magnética, obtenida originalmente a 300 m sobre el terreno, queda a 1500m de elevación para el proceso de interpretación directa mediante el algoritmo de Talwani (1965) [35]. Los resultados de la interpretación directa se pueden observar en las Figuras 29, 32, 35, 38 y 41. En la gráfica superior de cada figura se presentan las anomalías aeromagnéticas observada y calculada, mientras que la gráfica inferior se esboza el modelo de estructuras y cuerpos geológicos que se proponen como fuentes de las anomalías mencionadas. La explicación de la simbología empleada en el modelado se expone en la Tabla 1. Con la finalidad de observar a mayor detalle las estructuras geológicas inferidas en los perfiles, se dividieron en dos secciones de aproximadamente 40 km de distancia. Ambas secciones se discuten después de la discusión del perfil completo. 45 Tabla 1: Simbología y valores de susceptibilidad magnética de cada material presente la zona de estudio. Símbolo Abreviatura Material Susceptibilidad Qal Aluvión Cuaternario 0.000002 Qpu Pumicita Cuaternario 0.000013 QB Basalto Cuaternario 0.00015 QBvb Brecha Volcánica Básica Cuaternario 0.0009 TplA Andesita Terciario Neógeno 0.00035 TmplB Basalto Terciario Neógeno 0.0001 TmplTR-R Toba Riolítica Terciario Neógeno 0.0005 TomTR Toba Riolitica Terciario Paleógeno 0.00045 TomR-Ig Riolita - Ignimbrita Terciario Paleógeno 0.0006 TomA-TA Andesita - Toba Andesítica Terciario Paleógeno 0.0004 TmPA Pórfido Andesitico Rocas Ígneas Intrusivas 0.00004 KsTpaTR Toba Riolítica Cretácico Superior 0.001 KsGr-Gd Granito Granodiorita Rocas Ígneas Intrusivas 0.002 46 5.2.1 MODELO DEL PERFIL P-01 El perfil P-01 (Figura 29) tiene una longitud de 81.6 Km, con estaciones interpoladas cada 200 m tanto para la elevación como para los datos magnéticos. Figura 29: P-01; Anomalías magnéticas observada y calculada a 1500 m (superior) y modelo geológico-estructural (inferior). Distancia en km, profundidad en km y datos magnéticos en nT. Con la finalidad de reducir el problema de ambigüedad de la interpretación, se emplean la información geológica de mapas geológico-mineros [4] y valores de susceptibilidad magnética extraídos de la literatura (Tabla 1). En el modelo de la Figura 29, que contiene al perfil P-01, es posible apreciar la predominancia de rocas graníticas y tobas riolíticas cretácicas (Bloque Jalisco), y rocas volcánicas del Terciario (andesitas y basaltos del Eje Neovolcánico). 47 El basamento granítico del Cretácico muestra una leve inclinación hacia el sur, con profundidades que alcanzan casi los trece km, y además exhibe una fuerte dislocación que se inclina hacia el norte, localizada también en la parte sur. Dos grandes cuerpos de tobas riolíticas cretácicas ocurren en la mitad meridional; el primero de ellos con un comportamiento similar al dislocamiento del basamento granítico. Esta estructura puede estar asociada con una falla normal profunda inclinada hacia el sur, formando una depresión con el segundo cuerpo de tobas riolíticas; dicha depresión se propone rellena de basaltos terciarios. Figura 30: P-01’; Sección meridional del perfil P-01 (distancia de 0 - 40 km). Anomalías magnéticas observada y calculada continuadas a 1500 m (superior) y modelo geológico- estructural (inferior). Distancia y profundidad en km y datos magnéticos en nT. En la parte inicial del perfil P-01’ (Figura 30) se observa un gran bloque de Toba Riolítica del Cretácico Superior de 20 km de largo con dirección casi NW-SE, con una extravasación que logra aflorar durante los primeros 9 km en la superficie de la sección 48 geológica. Después del kilómetro 12 al 15 y 20 presenta una acumulación de material volcánico, causando el alto magnético del kilómetro 13. Cabe Mencionar la existencia de un desgarre de Basalto Terciario Neógeno que se depositó sobre el cuerpo de Topa Riolítica del Cretácico Superior, con un espesor superior a un kilómetro. Se infieren 2 formaciones basálticas del Cuaternario que afloran en el tramo de 9 a 14 km. Este depósito de Basalto tiene una profundidad aproximada
Compartir