Resumen:
El semidesierto de Zacatecas ha sido un gran impulsor del desarrollo económico que se ha sostenido a nivel industrial debido a la cantidad, calidad y cotización de minerales que se encuentran dentro de las 13 imponentes formaciones litológicas del estado. La región de estudio abarca aproximadamente 3,744 km2 y se caracteriza por tener rocas de composición basáltica, riolítica, rocas intrusivas como granodioritas, monzonitas de la Era Cenozoica (Terciaria) y también rocas metamórficas caracterizadas por el tipo de metamorfismo de contacto de grado medio a bajo. El análisis de multivariables fue utilizado para construir diagramas jerárquicos (dendrogramas) y mapas de similitudes. Los principales puntos magnéticos de alta susceptibilidad magnética lo obtienen las formaciones de la Era Cuaternaria, Terciaria, Cretácica Temprana y Jurásica Tardía, lo que implica que sobre estas capas de roca y en algunas estructuras minerales profundas se obtenga el mayor contenido de minerales ferromagnéticos.
Palabras clave:formación litológicaformación litológica, metamorfismo de contacto metamorfismo de contacto, análisis multivariado análisis multivariado, diagramas jerárquicos diagramas jerárquicos, susceptibilidad magnética susceptibilidad magnética, estructuras minerales estructuras minerales.
Abstract:
The semidesert region of Zacatecas has been a major driver of economic development sustained at an industrial level by the quantity, quality and price of the minerals found in 13 stunning lithological formations in this state. The study area of approximately 3,744 km2 is characterized by rocks of basaltic composition, rhyolite, intrusive rocks such as granodiorite, monzonite of the terciary period at the beginning of the Cenozoic Era, and also low to medium grade metamorphic rocks also characterized by the type of contact metamorphism of low average degree. Multivariate
analysis was used to build hierarchical diagrams (dendrograms) and maps of
similarities. The main spots of high magnetic susceptibility were obtained to
be in formations from the quaternary, terciary, early Creataceous and late
Jurassic periods, which implies that the highest content of ferromagnetic
minerals is to be found in this layers of rock and in certain deep mineral
structures.
Keywords: lithological formation, multivariate analysis, hierarchical diagrams, magnetic susceptibility, mineral structures.
Análisis de un sistema de información geografica para la interpretación de cartas de campo magnético y geológico del semidesierto de Zacatecas
Interpretation of magnetic field and geological charts of the semidesert of Zacatecas
Recepción: 25 Junio 2015
Aprobación: 21 Marzo 2016
La región del semidesierto de Zacatecas tiene importancia a nivel internacional, debido a que allí se localizan las principales formaciones geológicas que aumentan la explotación y producción futuras de minerales de fierro y magnesio. El oro se encuentra asociado a la mineralización de especularita, hematita, calcopirita y pirita. Las amplias variedades de rocas en la región, desde minerales metálicos hasta los no metálicos, permiten al geocientífico aportar e indagar el área de estudio con métodos actuales.
Con el resultado del método aplicado, que fue un análisis multivariable, se demuestra adecuadamente que hay importantes concentra-ciones de minerales en diferentes partes de las localidades del semidesierto zacatecano, y que debido a los emplazamientos, yacimientos, estructuras y formaciones que varían de dimensión, se pueden obtener interesantes valores de hematita, magnetita, cobre y oro, los cuales son sustanciales y necesarios para la industria en general.
El área de estudio se encuentra en la porción semiárida que comprende el estado de Zacatecas y se ubica entre las coordenadas 24° 45’, 24° 15’ de latitud norte a 101° 20’, 101° 40’ de longitud oeste. Ha sido afectada por eventos geológicos durante la Era Mesozoica que consistieron en cambios tectónicos consecuentemente acompañados por amplios eventos volcánicos producidos en tiempos laramídicos (Damon et al., 1981).
El basamento consiste en rocas metamórficas y sedimentarias proterozoicas, producto de los bloques de Oaxaquia, como indican Centeno García et al. (2008). Actualmente, el área de estudio forma parte de una zona que se encuentra en constante crecimiento de actividades de exploración y explotación de minerales dada su importancia económica, como el oro, plata, cobre y zinc. Existen varios estudios geológicos y geofísicos realizados en esta región por empresas, en los que se percataron de ciertas anomalías magnéticas y permitieron delimitar algunas estructuras minerales como es el caso de Chile Colorado, Brecha Azul, Brecha, entre otros documentados por el Servicio Geológico Mexicano (SGM, 2004a y b).
El objetivo del presente artículo fue lograr un agrupamiento de las unidades de rocas encontradas en el área de estudio para emplearlo en cartografía geológica y otras áreas del conocimiento. Esta técnica puede ser utilizada para la relación de litologías, reconocimiento potencial de contenido de altos puntos magnéticos y observación de anomalías magnéticas. Otra técnica nombrada en este trabajo es el llamado gradiente vertical, sinónimo de la definición de derivada vertical por Milligan y Gunn (1997), físicamente equivalente a la medición de campo magnético simultáneamente en dos puntos verticalmente, uno encima del otro, substrayendo los datos y dividiendo los resultados entre la separación vertical espacial de los puntos de medición.
De lo nombrado anteriormente, se dará un entendimiento del panorama magnético que permita a su vez promover beneficios e información adicional que sea trascendente a la hora de aprovechar recursos naturales. Los resultados sugieren que los valores de los principales puntos magnéticos de alta susceptibilidad magnética que son las localidades El Sol y La Luna, La Felicidad y La Escondida, se obtienen de las formaciones que datan del Cretácico Temprano y Jurásico Tardío y las formaciones litológicas son parte de una fuente autóctona, como las rocas cuaternarias y terciarias.
El análisis consistió en varias fases. Las más importantes fueron dos: en la primera fueron revisadas las variables del campo magnético total (CMT) con respecto a las formaciones litológicas y la segunda fue la mención de la técnica del gradiente vertical, de la reducción al polo con respecto a las formaciones litológicas. Para el análisis de las estructuras litológicas que se encuentran dentro de la región, se utilizó un método multivariables que consistió en un análisis de conglomerados (cluster en inglés) y un agrupamiento jerárquico de datos (AJD), definido por Facchinelli et al. (2001). Algunos métodos y estudios similares han sido integrados en bases de datos sobre el noreste de México, que según Aiken et al. (1997), en colaboración con el Instituto de Geografía de la UNAM en 1990, contribuyeron a la generación del modelamiento por computadora e interpretación de la corteza sísmica con métodos gravimétricos que para las formaciones geológicas que construyen la corteza terrestre hasta cierta profundidad es determinada por el alcance de este método. Las rocas de mayor densidad aumentan la aceleración y se realizan mediciones; es decir, se miden las atracciones laterales de la fuerza gravitatoria de un lugar al otro y posteriormente se intentan relacionar con las características geométricas del cuerpo, así como de sus contrastes de densidades. Algunas de las revisiones del modelamiento por computadora para propósitos de interpretación fueron realizadas a finales de la década de los setenta del siglo XX por el ingeniero hindú Bhattacharyya (1978).
Para cubrir los objetivos propuestos en este trabajo se consultaron las cartas geológicas mineras y magnéticas del SGM, así como diversos trabajos que contienen metadatos y reportes sobre temas de carácter geológico y geofísico-magnético. Se hizo un procesamiento de 11 cartas de escala 1:50,000, de las cuales, cuatro fueron cartas geológicas-mineras, cuatro geofísicas de carácter magnético y tres de derivada vertical, debido a que la carta magnética del primer gradiente vertical correspondiente a Concepción del Oro en 2013 no había sido realizada porque para producirlas se requiere equipo especializado como magnetómetros y programas de cómputo más sofisticados, como el procesamiento de datos en un sistema informático llamado Oasis. Hay cuatro localidades del municipio de Mazapíl, Concepción del Oro G14-C62, con curvas 40 nT, Cedros G14-C61 con curvas cada 20 nT y 4 pT/m, El Cardito G14-C71 con curvas cada 10 nT y 2 pT/m y Tanquecillos G14-C72 con curvas de 20nT y 10 pT/m, las que destacamos por su importancia geológica. Para la realización de este artículo no se hizo ningún trabajo de campo, básicamente consistió en la creación y procesamiento de datos vectoriales, cuya información se adquirió de las cartas geológicas antes mencionadas.
Se capturaron capas de información geográfica del campo magnético total (CTM), gradiente vertical y litológica en formato vectorial. La sobreposición de información con la mezcla de sus atributos, mediante capas de entrada, se complementa hasta que se consiga otra capa deseada, esto fue explicado por Davis (2001). Las operaciones permiten que la información de una capa sea superpuesta sobre la información de otra donde, dependiendo de los objetos, se pueden obtener diferente información de salida. En la clasificación de los datos se aplicó un método multivariado para agrupar instancias (variables) que permite el reconocimiento de patrones.
El método utilizado consistió en un análisis de agrupación jerárquica del tipo de proceso por conglomerados acumulativos (cluster analysis) propuesto por Jain y Dubes (1988) y para la visualización de la información se utilizaron gráficos jerárquicos que en inglés son conocidos como dendograms (dendogramas) y constituyen la representación que mejor ayuda a interpretar un análisis de conglomerados. La técnica de agrupación jerárquica es apropiada para mostrar las relaciones de las propiedades físicas de las rocas. Este análisis de grupos jerárquicos es una técnica exploratoria que revela agrupaciones naturales y representa, a su vez, las relaciones con las formaciones de las rocas. El paquete informático de visualización geográfica fue de gran utilidad; correspondientemente fue seleccionado conforme a su capacidad de visualizar las relaciones espaciales entre el ambiente de datos y las características de la superficie de la región.
El análisis tuvo varias fases, las más importantes fueron dos: en la primera se analizaron variables del campo magnético total (CMT) con respecto a las formaciones litológicas y en la segunda fue la mención de la técnica del gradiente vertical, de la reducción al polo con respecto a las formaciones litológicas, según Milligan y Gun (1997). Se aclara que los valores obtenidos en tablas es la sumatoria de toda la geometría de capas geográficas creadas en las cartas litológicas.
La formación Cupido ocupa mayor extensión superficial en el área de estudio, tiene 24,102.08 ha y es seguida por la formación Cuesta del Cura, la cual tiene 12,168.25 ha. Además, se encuentran tres formaciones litológicas que contribuyen a los valores de alto gradiente magnético y altos valores de intensidad magnético de 2,100 nT, mismas que pertenecen a las formaciones Cupido, Taraises y La Caja, sus valores de desviación estándar están sobre 747.52 nT y 770.61 nT, por lo que sus datos son dispersos, lo que provoca una variación grande con respecto al bajo gradiente del campo magnético total. Consecuentemente, debajo de estos valores se encuentra la formación Zuloaga, con valor de 2,085 nT y promedio de 791 nT. Esto evidencia cuerpos ocultos mineralizados en estas tres unidades litológicas, lo que causa la aparición de anomalías magnéticas en la región y, en consecuencia, entre sus capas se asimila un papel preponderante en la preservación de los depósitos minerales, lo que a su vez causa que se encuentren algunos halos de metamorfismo originados por el contacto de las secuencias sedimentarias con el cuerpo intrusivo. Ello produce gran porcentaje de mineralización compuesta por óxidos y sulfuros de fierro encima de la otra capa secuencial de la formación Zuloaga, esta puede ser la razón por la que en la Figura 1 se exponen altos valores de susceptibilidad para rocas sedimentarias marinas. En general, el contenido de magnetita es extremadamente variable y la susceptibilidad de las rocas puede ser considerable en el solapamiento entre diferentes litologías. Las rocas metamórficas también son variables en su carácter magnético, si la presión de oxígeno parcial es relativamente baja, la magnetita es reabsorbida, el fierro y oxígeno son incorporados dentro de otra fase mineral como el grado de metamorfismo incrementa (Kearey et al., 2002).
La desviación estándar de 5.91 a 23.50 nT y promedios de -40, -49 y -11 nT de las formaciones Nazas, La Joya y Taray exhiben un bajo gradiente magnético. Esto implica que los datos de la desviación estándar no son muy altos, están cerca de los valores promedio, la formación Taray es el basamento y es de gradiente menor, también tiene la característica de ser no magnético, decrece la cantidad de minerales ferromagnéticos.

Para la clasificación de los elementos en conglomerados se utilizaron diagramas jerárquicos que partiendo de un solo grupo van separados de los demás elementos y los cuales, cada vez se van haciendo más pequeños. A partir de estas unidades iniciales se van formando grupos de forma ascendente; las líneas superior e inferior con valores

ascendentes indican que cada elemento entre más se aleje de cero disminuye la similitud con respecto a sus componentes minerales; las observaciones más distantes son 4 y 12 las únicas mayores distancias y la 2, que no tiene grupo. La observación más cercana de elementos con menor distancia (más cercanas a cero) es 11, 8, 9, 13, 7, 10, 6, 15, 14. En el dendograma de las variables campo magnético total (CMT) y formaciones geológicas, las formaciones La Peña (11) y Taraises (8) forman un grupo en el primer nivel; las formaciones Indidura (13), La Caja (7), Cupido (10) y Zuloaga (6) forman el segundo grupo en el primer nivel. Las formaciones Ahuichila (15) y Caracol (14) se encuentran formando el tercer agrupamiento.
El primer grupo de rocas representado por las formaciones La Peña, Taraises, Cuesta del Cura, Indidura, La Caja, Cupido y Zuloaga, está asociado a los cambios progresivos del magma intermedio que se solidifica como cuerpo intrusivo dentro de estas capas de roca encajonante; el magma enriquecido interactúa con las diferentes unidades litológicas, que es la posible razón de anomalías fuertes y positivas magnéticas.
Las formaciones Caracol y Ahuichila están controladas por una clara asociación, sus contactos superiores se encuentran adyacentes, caracterizados por basaltos y sedimentos lacustres a consecuencia de efectos tectónicos; son formaciones afectadas por intensos plegamientos y probablemente comparten componentes minerales en estos depósitos continentales.
Si se toma en cuenta la sobreposición de capas de atributos de la derivada vertical, la combinación de los valores de las formaciones que están distribuidas en la región y los datos arrojados del campo magnético total, se obtiene que las formaciones Cuesta del Cura, La Peña y Cupido tienen los valores más altos de 282 pT/m, 280 pT/m y 279 pT/m, y sus valores mínimos, de -204 pT/m,-200 pT/m y -200 pT/m, implican que Cuesta del Cura es una de las unidades litológicas que conforma núcleos significativos de áreas mineralizantes granodioríticas y que por alguno de los diversos efectos como la presión parcial de alto oxígeno pudieron haber producido minerales accesorios como ilmenita, piroxenos, esfena y biotita, alanita, pirita Clark (1997); en la Figura 3 se reflejan sus diferentes capas y agrupaciones. Los valores máximo y mínimo se muestran en negritas en la Tabla 1.
Las observaciones que se percibieron en la formación Parras, composición de lutitas fisibles, consisten en que sus dominios magnéticos carecen de sus valores de intensidad, por lo que las capas litológicas no son geológicamente magnéticas debido a que las partículas de roca en su mayoría están estructuradas por minerales detríticos como arcillas de 0.0039 mm, los más finos en las rocas. Es improbable que entre sus minerales accesorios se encuentren minerales ferromagnéticos debido a su probable alteración del contenido de carbonato de sodio.
Como en el dendrograma anterior, las líneas superior e inferior con valores ascendentes indican que cada elemento entre más se aleje de cero disminuye la similitud con respecto a sus componentes minerales. Las observaciones más distantes son los grupos de formación Indidura (13), que no tiene grupo, skarn, hornfels (2), mármol (1), formación Caracol (14) y filita-arenisca (4). Las observaciones más cercanas, que son elementos con menor distancia de cero corresponden a 11, 8, 9, 15, 6, 7. En el dendrograma de las variables de gradiente vertical y formaciones geológicas se observan cuatro grupos, de los más significativos son 2: la formación La Peña (11), Taraises (8), forman el primer grupo correspondiente al primer nivel. De las formaciones Ahuichila (15) y Zuloaga (6) se forma el segundo agrupamiento, las rocas skarn, hornfels (2) y mármol (1) forman el tercer grupo y la formación Caracol (14) y filita-arenisca (4) forman el cuarto.


El resultado estadístico del primer grupo fue menos concluyente, mientras que el agrupamiento 2 confirma la explicación mineralógenico-magmática, quedan algunas dudas acerca de los grupos 3 y 4. Las posibles explicaciones son las siguientes:
— Las formaciones La Peña y Taraises, en el grupo 1, al igual que el grupo 2 de las formaciones Ahuichila y Zuloaga, están asociadas con otro factor mineralogénico; por ejemplo la presencia de rocas con sulfuros del mismo cuerpo intrusivo que los atraviesa.
— En las rocas skarn, hornfels y mármol están controladas por el proceso metamórfico regional y si se analiza la variabilidad somera del skarn es principalmente controlada por algunas anomalías magnéticas fuertes y positivas alrededor de los intrusivos riolíticos y granodioríticos en las localidades de Noche Buena, Concepción del Oro y Terminal de Providencia.
Los resultados derivados del presente estudio muestran los puntos interesantes de susceptibilidad magnética de uno de los semidesiertos más importantes de México y Latinoamérica. Estos

valores de susceptibilidad magnética se encuentran apegados a los datos expuestos obtenidos de las cartas litológicas y magnéticas de Cedros, Concepción del Oro, El Cardito y Tanquecillos expuestos en el SGM (1997, 1998, 2004, 2005, 2009, 2010) y representan una seria comprensibilidad debido a las implicaciones que el sector minero tiene sobre los recursos naturales.
Si se revisan los estudios y la aportación de datos de las cartas geológicas del SGM (2004), el valor alto magnético coincide en el área de Concepción del Oro, el área Noche Buena y Santa Rosa; tres de las más productivas de México. Se explica debido al mineral granodiorítico asociado al intrusivo que genera un mineral llamado skarn, en donde las principales estructuras formadas son chimeneas de dimensiones de hasta 500 m2, en donde se destacan las localidades El Capulín, Sol y Luna, La Escondida, La Laja, La felicidad, entre otros. Dentro de las rocas Mesozoicas, las formaciones Taraises, La Caja, Cupido y La Peña, como se explicó anteriormente, son formaciones donde sus estructuras minerales se encuentran profundas y ocultas.
Las diferencias de estos valores con respecto a los estudios de otras zonas del país no son similares. De acuerdo con un estudio de López Loera y Tristán González (2013), en el altiplano central de la parte centro oriental en el norte de la sierra San Miguelito, los valores son bajos con respecto a sus anomalías magnéticas de bloques para patrones de fallamientos normales, asimismo, para las anomalías de características negativas.
Mediante uno de los métodos más elegantes: la agrupación jerárquica del tipo de proceso por conglomerados acumulativos (cluster analysis) se puede aceptar el cumplimiento de agrupación para cada una de las formaciones litológicas, a su vez expresaría que estos patrones de agrupamiento operan no solo para formaciones litológicas y capas de roca, sino también para una buena parte de fenómenos aglomerativos, susceptibles de ser clasificados jerárquicamente como valores anómalos geoquímicos, entre otros. Es por ello que se recomienda reportar los resultados en estas representaciones gráficas para enriquecer el análisis con este modelo estadístico.
Cabe destacar que es escasa la literatura que aborda las características magnéticas en el estado de Zacatecas se exhiben más bien una serie de reportes de diferentes localidades y municipalidades en la que el autor de este trabajo cree que esta versión actual es una contribución interesante que resulta de una apertura de proceso de datos. Finalmente, el área del semidesierto y el país cuenta con varias contingencias de estudio, que actualmente ha colocado a México como líder de la producción, importación y exportación de oro.
Se utilizaron 11 mapas magnéticos y geológicos de la red SGM para conocer los agrupamientos jerárquicos, de los cuales esta investigación se benefició para el análisis de la región del semidesierto de Zacatecas, en especial algunas localidades del municipio de Mazapíl. Los resultados principales se obtuvieron de la información de los mapas del campo magnético total para el análisis multivariado de las diferentes capas geológicas y por medio de agrupación e interpretación se identificaron varias relaciones.
Las rocas de pórfido riolítico y tipo skarn y hornfels tienen similitudes en sus valores de intensidad del campo magnético total. Debido a su parentesco en remanencia y susceptibilidad magnética, ambos grupos de rocas están asociados con los procesos magmáticos. Las formaciones La Peña, Taraises, Cuesta del Cura, Indidura, Caja, Cupido y Zuloaga constituyen grupos y sus capas de rocas tienen altos valores de magnetización. La magnetización de estas rocas sedimentarias marinas es producida por el magma enriquecido que interactúa con las diferentes unidades litológicas, que es la posible razón de anomalías magnéticas fuertes y positivas; además, estos datos significan que debido a su contenido magnético semejante, las capas se encuentran solapadas posiblemente con intrusivos de la región asociados a un proceso de acreción en la porción de la mesa central. Este tipo de fenómeno sucede cuando cierta cantidad de sedimentos son empujados contra los continentes y sufren metamorfismo de alta presión, esto se puede asociar con la deformación del Cretácico Tardío y con la sucesión de cabalgamientos y depósitos originados por la Orogenia Laramide.
Referente al gradiente magnético y formaciones litológicas, la variabilidad de las formaciones La Peña y Taraises, ambas asociadas, es controlada por rocas padre de rocas magmáticas. El paquete informático SIG muestra datos de relaciones espaciales entre anomalías positivas y negativas en relación con minerales ferromagnéticos.
Por último, el agrupamiento de la edad del campo magnético total es similar al grupo de los datos de la derivada de reducción en el sentido en que fueron agrupados en relación con sus diferentes valores magnéticos, son más suaves, por lo que el gradiente es menos abstracto. Los valores obtenidos muestran que los principales puntos magnéticos de alta susceptibilidad (área de Noche Buena, Concepción del Oro y Santa Rosa) pertenecen a las formaciones de la Era Cretácica Temprana y Jurásica Tardía, que probablemente resultaron de una fuente autóctona, al igual que las rocas de la Era Cuaternaria y Terciaria.
Los sistemas de información geográfica son usados por miles de geocientíficos, técnicos y continúan expandiéndose. Esto ha provocado que el uso de esta tecnología haya sido asimilada por universidades, gobierno, industria e instituciones en las que actualmente el SIG (paquete informático) provee de acceso y conexión a internet para poder acceder a nuevos subconjuntos de datos y opciones de operatividad vía satélite, en los que varias disciplinas en conjunto tratan de entender los patrones, procesos y algunas consecuencias que tienen las formaciones litológicas dentro de la geología y la tierra, algunas zonas potenciales de mineral hacen mejorar y entender el proceso que interviene en ellas en donde se incorporan más métodos en lo que se refiere a estudios de focos de potencialidad magnética, así como la determinación de la cuantificación de los recursos naturales.
El autor desea agradecer los valiosos comentarios y contribuciones de Francisco Benita Maldonado y Simón Eduardo Carranco Lozada, los cuales sirvieron para mejorar considerablemente la calidad de este documento. También los comentarios y sugerencias de dos dictaminadores anónimos, así como a Mónica Ávalos Valladares su sugerencia sobre el traslado de parte de los textos al pie de figuras a párrafo, que mejoró la estructura del trabajo.
hcmcasas@gmail.com




