Abstract: The Cu- and Au-rich Tatatila-Las Minas IOCG skarn deposits in Veracruz (central-east Mexico) are circumscribed to the earliest stages of the Trans-Mexican Volcanic Belt (TMVB) and stand for a metallogenic province directly linked to its tectonomagmatic dynamics. This is the first well-documented case for such metallogenic province. These deposits were formed as skarns between rocks of the Mesozoic carbonate series and Miocene intermediate to acid hypabyssal rocks. New U-Pb zircon and40Ar/39Ar ages provide evidence for four epochs of magmatic activity in the area: (1) early Permian (Artinskian), in association with the Paleozoic basement, (2) late Oligocene to early Miocene suite of pre-TMVB intrusive rocks, (3) middle to late Miocene suite of early TMVB-related intrusive rocks, and (4) Pliocene intrusive and extrusive rocks of the TMVB, possibly associated with the Los Humeros post-caldera stage. The obtained ages range between 24.60 ± 1.10 and 19.04 ± 0.69 Ma for stage 2, and between 16.34 ± 0.20 and 13.92 ± 0.22 Ma for stage 3. Stage 2 corresponds to a magmatic stage unheard of in the area, until this study. Only stage 3 rocks are associated with the IOCG skarn mineralization, with retrograde stages dated at 12.44 ± 0.09 (chromian muscovite, phyllic association) and 12.18 ± 0.21 Ma (zircon, potassic association). Therefore, the ages of stage-3 intrusive rocks are interpreted to date the formation of the prograde skarn associations (mostly ~15.4 to <14 Ma). The petrogenetic affinity of stage-2 and stage-3 rocks is about the same-the main difference has to do with higher Y and Yb contents in stage-3 rocks (although no affinity with within-plate granites was found), which is suggestive of an interaction of their parental magmas with alkaline magmas that most likely belong to the conterminous and contemporaneous Eastern Mexico Alkaline Province. Petrological indicators (elemental and isotopic) in Cenozoic rocks consistently point to intermediate to acid, metaluminous, I- and S-type rocks that were emplaced in a subduction-related continental arc, within the medium- to high-potassium calc-alkaline series, with high-silica adakitic signatures due associated to deep-sourced magmas that underwent crustal contamination to some degree. The various possible sources for the magmas with adakitic signature in this context can be narrowed down to two of them that are not mutually exclusive: adakitic derived from subducted slab melting and melting-assimilation-storage-homogenization (MASH)-derived adakites. Both sources are, in principle, capable of generating magmas that would eventually produce magmatic-hydrothermal mineralizing systems with an associated variety of ore deposit types, including IOCG. Also, both possible sources for adakites are compatible with the renewed steepening of the subducted slab after a period of flat subduction, for the earliest stage in the evolution of the TMVB.
Keywords: IOCGIOCG,adakitesadakites,MioceneMiocene,Trans-Mexican Volcanic BeltTrans-Mexican Volcanic Belt,skarnskarn,magmatic-hydrothermalmagmatic-hydrothermal,iron oxidesiron oxides.
Resumen: Los skarns IOCG ricos en Cu y Au de Tatatila-Las Minas en Veracruz (centro-oriente de México) están circunscritos a los estadíos más tempranos de la Faja Volcánica Transmexicana (FVTM) e indican directamente la existencia de una provincia metalogenética vinculada a su dinámica tectonomagmática. Este es el primer caso bien documentado para dicha provincia metalogenética. Estos depósitos se formaron como skarns entre rocas de la secuencia carbonatada del Mesozoico y rocas hipabisales indermedias a ácidas del Mioceno. Los nuevos fechamientos U-Pb en zircón y40Ar/39Ar evidencian la existencia de cuatro épocas de actividad magmática en el área: (1) en el Pérmico temprano (Artinskiano), en asociación con el basamento paleozoico de las secuencias del Mesozoico, (2) un conjunto de intrusivos pre-FVTM entre del Oligoceno tardío y el Mioceno temprano, (3) un conjunto de intrusivos del Mioceno medio y tardío asociados a la FVTM, y (4) rocas intrusivas extrusivas del Plioceno de la FVTM, posiblemente asociadas a los depósitos del estadio post-caldera de Los Humeros. Las edades obtenidas varían entre 24.60 ± 1.10 y 19.04 ± 0.69 Ma para el estadío 2, y entre 16.34 ± 0.20 y 13.92 ± 0.22 Ma para el estadío 3. El estadío 2 corresponde a una etapa magmática hasta el presente estudio desconocida en el área. Sólo las rocas del estadío 3 están asociadas a las mineralizaciones de skarn IOCG, cuyas etapas retrógradas han sido fechadas en 12.44 ± 0.09 (moscovita crómica, asociación fílica) y 12.18 ± 0.21 Ma (zircón, asociación potásica). Por tanto, las edades de las rocas intrusivas del estadío 3 se interpretan como parte de las asociaciones de skarn prógrado (mayormente, de ~15.4 a <14 Ma). La afinidad petrogenética de las rocas correspondientes a los estadíos 2 y 3 es prácticamente la misma-su principal diferencia estriba los contenidos más altos de Y e Yb en rocas del estadío 3 (aunque no se encontró afinidad alguna con granitos de intraplaca), lo cual sugiere la interacción de sus magmas primigenios con magmas alcalinos que posiblemente pertenecieron a la contigua y contemporánea Provincia Alcalina Oriental Mexicana. Los indicadores petrogenéticos (elementales e isotópicos) en las rocas del Cenozoico apuntan consistentemente a rocas intermedias a ácidas, metalumínicas, de tipo I y S, emplazadas en un arco continental debido a subducción y pertenecen a las series calci-alcalinas de potasio medio a alto, con (mayormente) firmas de adakitas altas en sílice debidas a un origen profundo de magmas que experimentaron cierto grado de contaminación cortical. La diversidad de posibles orígenes para las fimas adakíticas en este contexto pueden reducirse a sólo dos de ellas, que no son mutuamente exclusivas: adakitas derivadas de la fusión de la placa subducida y adakitas derivadas de procesos tipo fusión-asimilación-almacenamiento-homogeneización (MASH, por sus siglas en inglés). Ambas fuentes, en principio, poseen la capacidad de generar magmas que eventualmente pudieran producir sistemas mineralizantes magmático-hidrotermales con una cierta variedad de tipos de depósitos minerales asociados, incluyendo depósitos IOCG. Además, ambas posibles fuentes de adakitas son compatibles con la reverticalización de la placa subducida tras un periodo de subducción plana para el estadío más temprano en la evolución de la FVTM.
Palabras clave: IOCG, adakitas, Mioceno, Faja Volcánica Transmexicana, skarn, magmático-hidrotermal, óxidos de hierro.
Articles
The Miocene Tatatila-Las Minas IOCG skarn deposits (Veracruz) as a result of adakitic magmatism in the Trans-Mexican Volcanic Belt
Los depósitos de tipo skarn IOCG miocénicos de Tatatila-Las Minas (Veracruz) como resultado del magmatismo adakítico de la Faja Volcánica Trans-Mexicana
Sociedad Geológica Mexicana A.C.
Received: 28 October 2020
Accepted: 10 May 2020
Funding source: CONACyT
Contract number: 155662
Recent assessment has shown that the metallogenic potential of the mid-Miocene to Holocene Trans-Mexican Volcanic Belt (TMVB) and the potential of Miocene to Holocene ore deposits in Mexico are greater than previously believed (Camprubí, 2009,2013;Clark and Fitch, 2009;Poliquin, 2009;Jansenet al., 2017;Camprubíet al., 2019,2020;Fuentes-Guzmánet al., 2020). The metallogeny of Miocene to Holocene epochs in Mexico is, in fact, distributed across several regions, namely (1) the southernmost part of the Sierra Madre Occidental, in association with its last ignimbritic flare-up, (2) the Trans-Mexican Volcanic Belt (TMVB), (3) the southern part of the Eastern Mexico Alkaline Province (EMAP) and northern Chiapas, (4) the easternmost part of the Sierra Madre del Sur (in Oaxaca), and (5) the Gulf of California. As (a) the easternmost ending of the TMVB coincides with the N-S geographic distribution of the EMAP, (b) the metallogeny of the TMVB is still poorly understood, and (c) there is a wide variety of types of ore deposits across the EMAP-including, skarns, metalliferous porphyries, epithermal deposits, IOCG deposits and carbonatites-, the identification of whether an ore deposit in such a region is geologically associated with the TMVB or the EMAP is not a straightforward task.
The Tatatila-Las Minas district in Veracruz State is located precisely in the region in which the TMVB and the EMAP overlap geographically, in the Palma Sola area. The ore deposits in the Tatatila-Las Minas have a magmatic-hydrothermal origin and are essentially Cu-Au iron oxide skarns, part of the IOCG “clan”, and epithermal deposits (Camprubí, 2013). Therefore, in order to investigate the origin of these deposits, the first necessary step would be to elucidate their genetic affinity with either magmatic province.Camprubí (2013) deduced a plausible age of ~11 Ma and some affinity with alkaline magmatism for the deposits in the Tatatila-Las Minas district, based onNegendanket al. (1985) andFerrariet al. (2005a), which linked the Palma Sola massif with the EMAP. However, the middle Miocene to Recent alkaline and calc-alkaline volcanism of the Palma Sola area was ascribed to the TMVB, and to the subduction along the Pacific trench, as inBeschet al. (1988),Gómez-Tuenaet al. (2003), andOrozco-Esquivelet al. (2007). The relevance of the EMAP, besides its petrotectonic affinity, as a major metallogenic province was already stressed byCamprubí (2009,2013). However, the age of magmatism with which these ore deposits were plausibly associated corresponds well to the middle and late Miocene arc at the beginning of the TMVB (~19 to 10 Ma;Gómez-Tuenaet al., 2005,2007).
In summary, we may use as a starting hypothesis the fact that neither the EMAP nor the TMVB are implausible magmatic provinces to have produced the parental magmatism to the Tatatila-Las Minas deposits. The implications for regional mineral exploration that may arise from either possibility are very different, nonetheless. In this paper, we analyze the petrologic affinity of the hypabyssal intrusive bodies with which the formation of the IOCG deposits of the Tatatila-Las Minas district is associated. This will enable a discrimination between the ascription of these deposits to the metallogeny of the TMVB or the EMAP. The proximal-to-source character of these magmatic-hydrothermal deposits (i.e., iron skarns) allows to soundly elucidate the linkage between the magmatism and the hydrothermal activity that generated the deposits. In addition, this paper contributes to a long-standing program that aims to the geochronological characterization of Mexican mineral deposits and the geologic events with which they are genetically associated (Camprubíet al., 2015,2016a,2016b,2017,2018,2019,2020;Farfán-Panamáet al., 2015;Martínez-Reyeset al., 2015;González-Jiménezet al., 2017a,2017b;Enríquezet al., 2018;Fuentes Guzmánet al., 2020) to better constrain the metallogenic evolution of Mexico, as documented byCamprubí (2009,2013,2017).
The Tatatila-Las Minas mining district is located in the central-eastern part of the state of Veracruz (Figure 1) within the Palma Sola massif. It is characterized by the intrusion of Neogene stocks. Stock compositions are described to vary between gabbro and granodiorite, with dominantly monzodioritic to dioritic compositions, and intruded middle Jurassic, red beds and lower Cretaceous carbonate rocks. The latter rocks are part of the continental to marine sequences of the Sierra Madre Oriental that were deformed during the orogenic pulses of the Mexican Fold-and-Thrust Belt between the late Cretaceous and the Paleocene (Centeno-García, 2017;Fitz-Díazet al., 2018; and references therein). The Middle Jurassic red bed sequence in the area correlates with the Cahuasas Formation, and is overlain by carbonates and lutites of the Pimienta (Tithonian-Barriasian) and Orizaba (Albian-Cenomanian) formations. The host carbonate series in the study area consists essentially of platform carbonates that correspond to the Orizaba Formation (Ortuño-Arzateet al., 2005). The Lower Cretaceous sequence unconformably overlies Permo-Triassic schists intruded by granitic rocks. The latter can be mistaken for Neogene intrusive bodies with similar compositions, as the thick vegetation cover commonly hinders their visualization and the identification of the lithologic contacts; both groups of intrusive rocks come in contact by faulting in the northernmost termination of the mineralized area (Figure 1). The Mesozoic sedimentation was controlled by the horst-and-graben configuration that resulted from the opening of the Gulf of Mexico during the breakup of Pangea (Martini and Ortega-Gutiérrez, 2018), thus developing simultaneously shallow platforms and relatively deep open-sea facies, hence the Córdoba platform (Ortuño-Arzateet al., 2005) on which the upper Jurassic and Lower Cretaceous sedimentary units developed.
Figure 1
Geological map of the Tatatila-Las Minas mining district, east of the Palma Sola massif. Adapted fromServicio Geológico Mexicano (2007,2010). Purple circles denote the location of samples on which this study is based, with indication of the obtained ages.
Neogene intrusive bodies generated typical skarn associations, with prograde mineralization by contact metamorphism (Ca silicate-rich) that was followed by retrograde IOCG-type hydrothermal stages of mineralization (Figure 2). Such intrusive bodies made up a NE-SW striking ~20 km long and ~10 km wide intrusive ensemble whose composition varies from gabbro to granodiorite, with dominantly monzodioritic to dioritic compositions (see below) with phaneritic textures. These rocks typically contain hornblende, biotite, pyroxenes, apatite and zircon this two as accessory mineral (Figure 3). Some andesite dykes, up to 30 m long and ~2 m thick crosscut the intrusive ensemble and predate the mineralization. A sequence of andesitic, basaltic and dacitic hypabyssal, this with porphyritic texture include plagioclase phenocrysts and volcanic rocks postdates the mineralization and the emplacement of the associated intrusive rocks, and comprises a variety of deposits, including volcanic conglomerates, tuffs, ash-fall and pyroclastic deposits. Such rocks are interpreted as distal Pliocene deposits associated with the post-caldera deposits of Los Humeros caldera (Carrasco-Nuñezet al., 2018;Dorantes-Castro, 2016;Sarabia-Jacinto, 2017). Ages for the Palma Sola area to the east of the Tatatila-Las Minas area were 14.6 ± 0.3 (U-Pb, zircon) and 11 ± 0.87 Ma (K-Ar, biotite), were reported byPoliquin (2009) andMurillo-Muñetón and Torres-Vargas (1987), respectively. These correspond to the ensemble of hypabyssal and volcanic rocks that allowedCamprubí (2013) to deduce a tentative age of ~11 Ma for these ore deposits, which is also constrained by the formation of capping volcanic rocks between 9 and 6.6 Ma. Contact metamorphism and mineralization of the fresh carbonate rocks can be observed in the conspicuous formation of marble in a 300 to 400 m wide zone that shows an outward decreasing degree of recrystallization. Skarn associations are distributed in the classic zonation from endoskarn to exoskarn. Endoskarns consist of grossular-andradite, clinopyroxene, and quartz in prograde associations, and magnetite, chalcopyrite, bornite, and native gold in retrograde associations (Figure 2). Exoskarns consist of wollastonite, clinopyroxene, potassium feldspar, quartz, epidote, and chromian muscovite (“fuchsite”;Figure 2).
Figure 2
Selected aspects of the IOCG skarn mineralization at the Tatatila-Las Minas deposits showing both prograde (garnet and tourmaline) and retrograde (actinolite and fuchsite) associations. (A) Hand specimen showing a garnet-rich prograde association followed by an actinolite- and fuchsite-rich retrograde association in the Santa Cruz mine. (B) Photomicrography of a garnet and tourmaline prograde association followed by an actinolite and fuchsite retrograde association; transmitted light, crossed polars; same sample as in A. Fuchsite separates from A and B were dated by argon geochronometry in this study. (C) Hand specimen of prograde patchy to partially banded magnetite ore; El Dorado mine. (D) Hand specimen of banded exoskarn magnetite- and chalcopyrite-rich retrograde ore, with martitized magnetite; El Dorado mine. Key: Amp = amphibole-group minerals (actinolite), Cal = calcite, Ccp = chalcopyrite, Ep = epidote, Fuch = chromian muscovite or “fuchsite”, Grt = garnet-group minerals (grossular-andradite), Hm = hematite, Mag = magnetite, Mc = malachite, Py = pyrite, Tur = tourmaline.
Figure 3
Photomicrographs of representative hypabyssal bodies, unaffected by hydrothermal alteration, associated with IOCG skarn mineralization in the Tatatila-Las Minas district. (A) Quartz-monzodiorite showing euhedral apatite crystals within plagioclase phenocrysts, Santa Cruz mine; transmitted light, crossed polars. (B) Quartz-monzodiorite showing myrmekitic intergrowths, surrounded by hornblende and plagioclase phenocrysts, La Virgen mine; transmitted light, crossed polars. (C) Quartz-monzodiorite showing hornblende phenocrysts, Santa Cruz mine; transmitted light, crossed polars. (D) Quartz-monzodiorite showing euhedral zircon crystals within potassium feldspar, Santa Cruz mine; transmitted light, crossed polars. (E) Monzodiorite showing hornblende intergrown with magnetite, Carbonera mine; plane-polarized transmitted light. (F) Monzodiorite showing euhedral hornblende, biotite and apatite crystals within a plagioclase-potassium feldspar assemblage, Carbonera mine; plane-polarized transmitted light. (G) Monzogranite showing rock-forming biotite crystals, same sample as in F, Rancho La Virgen; transmitted light, crossed polars. (H) Monzogranite showing late biotite crystals intergrown with magnetite, Rancho La Virgen; plane-polarized transmitted light. Key: Ap = apatite, Bt = biotite, Fp = potassium feldspar, Hb = hornblende, Mt = magnetite, Pl = plagioclase, Qz = quartz, Zr = zircon.
Mining activity in the study area can be dated back to pre-colonial epochs, when the native population of Chiconquiaco obtained gold that was mainly destined to fulfill the contributions imposed upon them by their Aztec overlords. Formal mining by the Spaniards can be dated back to at least 1680, when the exploitation of large high-grade gold and silver bonazas has been documented (Castro-Moraet al., 1994). Mining and exploration have remained intermittently active in the area ever since (Viniegra, 1965;Castro-Moraet al., 1994;Servicio Geológico Mexicano, 2007). By 1996, the exploration endeavors carried out by International Northair, in association with Battle Mountain Gold Co., allowed location of relevant Au-Cu-Fe resources in a broad area. In 2006, Bell Resources Corp. took over the property in the Las Minas area and subsequently assigned the mining rights to Chesapeake Gold Corp.
The formation of this Au-Cu-Fe rich area is generally acknowledged to belong to an IOCG model with overimposed late epithermal veins (Servicio Geológico Mexicano, 2007;Camprubí, 2009,2013;Dorantes-Castro, 2016;Castro-Moraet al., 2016;Sarabia-Jacinto, 2017). The metal grades in the deposit range between 1 and 39.3 ppm Au, between 4.11 and 127 ppm Ag, and between 0.64 and 11.7% Cu; inferred reserves are 719000 Oz Au equiv., and indicated reserves are 304000 Oz Au equiv. (Castro-Moraet al., 2016).
Representative samples from the Neogene intrusive ensemble were collected in the Tatatila-Las Minas mineralized area (47 samples; purple circles inFigure 1) in order to characterize the petrologic affinity and age of skarn-generating intrusive bodies, as well as the age of hydrothermal activity itself. The ages of intrusions are considered as representative of the age of prograde mineralization in IOCG skarns, and hydrothermal assemblages correspond to retrograde stages of these deposits. The representativeness of such samples with regard to the formation (or postdating) of mineralized bodies was determined on the basis of their distribution, their possible association with mineralized bodies, and the types of rocks thereby represented, after thorough cartography and sampling. All the analyzed samples were examined by means of petrographic studies in order to ensure that no alteration would cause any disturbances to the geochemical or geochronological analyses.
Elemental analyses were carried out on 15 g aliquots from samples at a 200 mesh. The two dated samples from retrograde hydrothermal associations are chromian muscovite, which correspond to high-temperature phyllic assemblages from the Las Minas area, and zircon within pervasive potassic alteration assemblages from the Tatatila area.
Multi-elemental geochemical analyses of host rocks were carried out by means of X-ray fluorescence (XRF) with a Rigaku Primus II equipment available at the Laboratorio Nacional de Geoquímica y Mineralogía (LANGEM) in accordance with the procedure described byLozano-Santa Cruzet al. (1995); results are presented inTable 1. Trace and rare-earth elements (REE) were analyzed by means of inductively coupled plasma quadrupole mass spectrometery (Q-ICP-MS) with a Termo ICap Qc equipment, coupled to a collision/reaction cell (He, N2, NH3 and O2) in order to minimize spectral interference, the procedure described byMoriet al. (2007), at the Laboratorio de Estudios Isotópicos (LEI) of the Centro de Geociencias (CGeo-UNAM). The obtained data are presented inTable 2. For Sr, Nd and Pb isotopic analyses, a Thermo Fisher Neptune Plus mass spectrometer available at the CGeo-UNAM. Sample preparation and measurement procedures for Sr-Nd-Pb isotopic analyses are described inGómez-Tuenaet al. (2003) for LDEO.87Sr/86Sr ratios obtained in both labs were normalized to86Sr/88Sr = 0.1194 and corrected to a NBS-987 standard ratio of87Sr/86Sr = 0.710230, and143Nd/144Nd ratios were normalized to146Nd/144Nd = 0.72190 and corrected to a La Jolla standard value of143Nd/144Nd=0.511860. At LDEO, Sr and Nd were measured by dynamic multicollection, with each analysis consisting of ~120 isotopic ratios. Sr ratios were measured using tungsten filaments and a TaCl4 activator solution (Birck, 1986). Nd isotopes were measured as NdO+. During five separate analysis intervals the measured values of the NBS-987 standard were87Sr/86Sr = 0.710245 ± 0.000016 (2σ, n = 4); 0.710271 ± 0.000014 (2σ, n = 6); 0.710274 ± 0.000016 (2σ, n=18); 0.710310 ± 0.000013 (2σ, n = 5); 0.710261 ± 0.000012 (2σ, n = 10). The measured143Nd/144Nd ratio of the La Jolla standard at LDEO was 0.511836 ± 0.000013 (2σ, n = 15), as ofTodtet al. (1996), according to the procedure described by Moriet al. (2007), obtained data are presented inTable 3.
Major elements in host intrusive rocks to the Tatatila-Las Minas IOCG deposits. All values in wt.%. Asterisks (*) correspond to analyses inDorantes-Castro (2016).
Key: LOI = loss on ignition.
Trace elements in host intrusive rocks to the Tatatila-Las Minas IOCG deposits. All values in ppm unless otherwide noted. Asterisks (*) correspond to analyses inDorantes-Castro (2016).
Sr, Nd and Pb isotopic values of selected samples from intrusive rocks associated with IOCG skarn mineralization in the Tatatila-Las Minas area.
The two dated samples from retrograde hydrothermal associations are chromian muscovite that corresponds to high-temperature phyllic assemblages from the Las Minas area, and zircon within pervasive potassic alteration assemblages from the Tatatila area, the40Ar/39Ar analyses of samples from intrusive rocks were carried out at the Noble Gas Laboratory, Pacific Centre for Isotopic and Geochemical Research, University of British Columbia (Vancouver, British Columbia, Canada). The mineral separates were step-heated at incrementally higher powers in the defocused beam of a 10W CO2 laser (New Wave Research MIR 10) until fused. The gas evolved from each step was analyzed by a VG5400 mass spectrometer equipped with an ion-counting electron multiplier. All measurements were corrected for total system blank, mass spectrometer sensitivity, mass discrimination, radioactive decay during and subsequent to irradiation, as well as interfering Ar from atmospheric contamination and the irradiation interferences of Ca, Cl and K. The plateau and correlation ages were calculated using the Isoplot 3.09 software (Ludwig, 2003). Errors are quoted at the 2-sigma (95% confidence) level and are propagated from all sources except mass spectrometer sensitivity and age of the flux monitor. The full results and spectra are reported inAppendices 1 and2 and summarized inFigure 4.
Figure 4
Outlines of40Ar/39Ar age spectra (plateau ages) of intrusive host rocks to the IOCG skarn deposits in the Tatatila-Las Minas district, Veracruz.
The40Ar/39Ar analysis were performed at the Geochronology Laboratory of the Departmento de Geología, Centro de Investigación Científica y Educación Superior de Ensenada (CICESE, Mexico). The argon isotope experiments were conducted on a few flakes of fuchsite, hornblende, K-feldspar and biotite. The mineral grains were heated with a Coherent Ar-ion Innova 370 laser. The extraction system is on line with a VG5400 mass spectrometer. The sample and irradiation monitors, were irradiated in the Uenriched research reactor of University of McMaster in Hamilton, Canada, at position 5C. To block thermal neutrons, the capsule was covered with a cadmium liner during irradiation of chromian muscovite (“fuchsite”;Figure 5A and5B) from the skarn gangue association in IOCG mantos, Santa Cruz mine (sample SC-1). The mineral grains were heated with a Coherent Ar ion Innova 370 laser. The extraction system is on line with a VG5400 mass spectrometer. The sample and irradiation monitors were irradiated in the U-enriched research reactor of University of McMaster in Hamilton, Canada, at position 5C. To block thermal neutrons, the capsule was covered with a cadmium liner during irradiation. To determine the neutron flux variations, aliquots of the irradiation monitor FCT-2 sanidine (28.201 ± 0.046 Ma;Kuiperet al., 2008) were irradiated alongside sample SC-1. Upon irradiation the monitors were fused in one step while the fuchsite sample was step-heated. The argon isotopes were corrected for blank, mass discrimination, radioactive decay of37Ar and39Ar, and atmospheric contamination. For the Ca neutron interference reactions, the factors given byMasliwec (1984) were used. The decay constants recommended bySteiger and Jäger (1977) were applied in the data processing. The equations reported byYorket al. (2004) were used in all the straight line fitting routines of the argon data reduction.40Ar/39Ar data are presented inAppendices 1 and2, which includes the results of the individual steps, and the integrated, plateau and isochron ages, and their synthetic version inFigure 5. The analytical precision is reported as standard deviation (2σ). The error in the integrated, plateau and isochron ages includes the scatter in the irradiation monitors. With the exception of the first fraction, a well-defined straight line, with mean squared weighted deviations (MSWD) of 0.55 for n = 6, indicates an isochron age of 12.49 ± 0.09 Ma.
Figure 5
40Ar/39Ar age spectra (plateau and isochron ages) of a chromian muscovite (“fuchsite”) from the magmatic-hydrothermal retrograde assemblage of the IOCG skarn deposit in the Santa Cruz mine, Tatatila-Las Minas district, Veracruz.
Zircon crystals were separated by means of panning from samples selected for U-Pb dating that are representative of various sets of rocks in the area: Au-Ag mineralized vein from Tatatila (sample TMG-5), and granodiorite to granite samples from the Santa Cruz (samples TMG-24 and SC-2), Carboneras (CR-5), Escalona (ES-3), Cinco Señores (5S-1), Boquillas (BQ-1), and Rancho Virgen (RV-2) areas. The sizes of the collected zircon crystals range between 20 and 90 μm in length. The U-Pb zircon analyses were performed with a quadrupole Thermo-X series ICP-MS with an Excimer (193 nm) laser ablation system by Resonetics, at the Isotopic Studies Laboratory (LEI), CGeo-UNAM, and following the procedure described bySolariet al. (2010). The data reduction was performed with the aid of the UPb.age in-house software (Solari and Tanner, 2011) and plotted with the Isoplot 3.0 software (Ludwig, 2003). See further technical aspects inGonzález-Leónet al. (2017). U-Pb ages are displayed inFigures 6 and7,Table 4 andAppendix 3.
Figure 6
Tera-Wasserburg U-Pb concordia diagrams and plots of weighted averages of individual206Pb/238U ages of analyzed zircons, and pre-ablation SEM-CL images of zircons from a granodiorite intrusive from the Santa Cruz Mine (A), and from a potassic alteration assemblage that was pervasively developed on a granite-granodiorite intrusion in the village of Tatatila (B), from the Tatatila-Las Minas district, Veracruz. Solid-line ellipses, with blue square centers, are data used for age calculations; gray-line ellipses are data excluded from age calculations due to different degrees of Pb-loss and/or zircon inheritance. All U-Pb data are plotted with 2-sigma errors and all calculated weighted mean ages are also listed at the 2-sigma level. Original U(Th)-Pb data can be found for inspection inTable 5.
Figure 7
Tera-Wasserburg U-Pb concordia diagrams for zircons from various intrusive bodies in the Tatatila-Las Minas area. (A) Post-mineralization dyke. (B to D) Syn-mineralization hypabyssal bodies whose age can be attributed to the prograde skarn associations. (E) Granitic intrusive that corresponds to the Permo-Triassic basement. Solid-line ellipses, with black square centers, are data used for age calculations; gray-line ellipses are data excluded from age calculations due to different degrees of Pb-loss and/or zircon inheritance. All U-Pb data are plotted with 2-sigma errors. Original U(Th)-Pb data can be found for inspection inAppendix 3.
U-Th-Pb analytical data for LA-ICPMS spot analyses on zircon grains for granitic units in Tatatila de Las Minas, Veracruz, Mexico.
n = 30 Mean206Pb/238U Age = 12.18 ± 0.21 (2 sigma, MSWD = 1.5; n = 26)# U and Th concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.† Isotopic ratios are corrected relative to 91500 standard zircon for mass bias and down-hole fractionation (91500 with an age ~1065 Ma; Wiedenbeck et al., 1995). Isotopic 207Pb/206Pb ratios, ages and errors are calculated following Paton et al. (2010).* All errors in isotopic ratios are in percentage whereas ages are reported in absolute and given at the 2-sigma level. The weighted mean 206Pb/238U age is also reported in absolute values at the 2-sigma level. The uncertenties have been propagated following the methodology discussed by Paton et al. (2010).* *Rho is the error correlation value for the isotopic ratios 206Pb/238U and 207Pb/235U calculated by dividing these two percentage errors. The Rho value is required for plotting concordia diagrams.* **Percentage discordance values are obtained using the following equation (100*[(edad 207Pb/235U)-(edad 206Pb/238U)]/edad 207Pb/235U) proposed by Ludwig (2001). Positive and negative values indicate normal and inverse discordance, respectively. Individual zircon ages in bold were used to calculate the weighted mean 206Pb/238U age and MSWD (Mean Square of Weigthed Deviates) using the computacional program Isoplot (Ludwig , 2003).
Summary of geochronometric data obtained for host intrusive rocks and IOCG mineralization at the Tatatila-Las Minas area.
The U-Pb ages of zircon crystals from granite, granodiorite, quartz-monzonite and monzodiorite are displayed inFigures 6 and7, inTable 4, andAppendices 3 and4. The sample from Carboneras (CR-5) yielded a U-Pb concordia lower intercept at 15.05 ± 0.94 Ma (MSWD = 2.5, n = 19;Figure 7C). Two samples from the Santa Cruz mine were dated; sample TMG-24 yielded a U-Pb concordant age at 15.27 ± 0.36 Ma (MSWD = 2 n = 14;Figure 6A), and sample SC-2b a weighted mean U-Pb age at 14.33 ± 0.38 Ma (MSWD = 2.6, n = 9;Figure 7B). The sample from Cinco Señores (5S-1) yielded a U-Pb weighted mean age at 15.09 ± 0.48 Ma (MSWD = 4.0, n = 8;Figure 7D).40Ar/39Ar determinations in host intrusive samples as granodiorite, granite, monzodiorite and quartz-monzonite yielded two groups of ages: (A) late Oligocene to early Miocene, between 22.12 ± 0.74 and 19.04 ± 0.69 Ma for a pre-mineralization suite of intrusive bodies, and (B) middle to late Miocene, between 16.34 ± 0.20 and 13.92 ± 0.22 Ma for a syn-mineralization suite of intrusive bodies, all reported ages correspond to plateau ages.
The samples for40Ar/39Ar different minerals such as biotite, hornblende, K-feldspar and fuchsite, were separated from each sample for analysis.
The40Ar/39Ar determination in hydrothermal chromian muscovite (“fuchsite”) of the Santa Cruz mine yielded a plateau age of 12.49 ± 0.09 Ma (isochron age at 12.39 ± 0.1 Ma;Figure 5). The sample (TMG-5) from a potassic alteration assemblage that was pervasively developed on a granite-granodiorite intrusion in the village of Tatatila (thus corresponding to hydrothermal associations) yielded a U-Pb age of 12.18 ± 0.21 Ma, (2σ, MSWD = 1.5; n = 26;Figure 6B).
The sample ES-3 from Escalona corresponds to a dyke that crosscuts the IOCG mineralization and yielded a U-Pb weighted mean age at 4.11 ± 0.11 Ma (MSWD = 0.53, n = 8;Figure 7A). A sample from Boquillas (BQ-1a, BQ-1b) yielded a U-Pb weighted mean age at 286 ± 2 Ma (MSWD = 1.02, n = 6; Artinskian, early Permian;Figure 7E).
The intrusive rocks associated with the formation of IOCG deposits in the Tatatila-Las Minas area span compositions between those of sub-alkaline gabbros and granodiorites, and mostly concentrate in the granite, diorite and monzodiorite fields (Figure 8A). The geochemical affinity of the rocks is essentially metaluminous (Figure 8B), calc-alkaline (Figure 8C), and they plot within the fields of volcanic-arc granites (VAG) (Figure 9A) and I- and S-type granites (Figure 9B). Some samples have adakitic signatures (Figure 9D), mostly of the high-silica type (Figure 9E), thus indicating that their compositional variation is controlled mainly by partial melting (Figure 9C). Light rare-earth and large-ion lithophile elements (LREE and LILE) are slightly enriched in such rocks (Figure 10) with respect to heavy rare-earth and high field strength elements (HREE and HFSE), as is characteristic for rocks associated with subduction, and conform with the results obtained byDorantes-Castro (2016). Radiogenic isotope data range as follows:87Sr/86Sr between 0.7040 and 0.7059, ɛSr between -11.4 and 19.9,143Nd/144Nd between 0.5123 and 0.5128, ɛNd between -6.6 and 3.2, and206Pb/204Pb between 18.65 and 18.75 (Table 3;Figure 11). The distribution of such data is in accordance with that determined byGómez-Tuenaet al. (2003) for rocks from the Trans-Mexican Volcanic Belt.
Figure 8
Petrological discrimination diagrams from major elements in intrusions associated with IOCG skarn mineralization in the Tatatila-Las Minas district, Veracruz. (A) Silica vs. alkaline element bivariant diagram, adapted fromCoxet al. (1979). (B) Alumina saturation diagram, adapted fromFrostet al. (2001), with compositions of skarns fromMeinert (1995). (C) AFM diagram, adapted fromIrvine and Baragar (1971).
Figure 9
Petrological discrimination diagrams from trace elements in intrusive rocks associated with IOCG skarn mineralization in the Tatatila-Las Minas district, Veracruz. (A) Y+Nbvs. Rb, Yvs. Nb, Ta+Ybvs. Rb, and Ybvs. Ta diagrams for discriminating tectonic settings, adapted fromPearceet al. (1984). (B) Discrimination diagram for different granite sources, adapted fromWhalenet al. (1987). (C) Discrimination diagram for the generation of magmas by fractional crystallizationvs. variable degree of partial melting, adapted fromThirlwallet al. (1994). (D) Discrimination diagram for adakitic affinity, adapted fromMartin (1986) with chondrite-normalized valuesSun and McDonough (1989). (E) Discrimination diagrams for high-silica (HSA) and low-silica adakites (LSA), adapted from Martin and Moyen (2002, 2003) and Martinet al. (2005). Key: HSA = high-silica adakites (>60% SiO2), LSA = low-silica adakites (<60% SiO2), ORG = ocean ridge granites, VAG = volcanic arc granites, syn-COLG = syn-collision granites, WPG = within plate granites.
Figure 10
Spider diagrams of REE (A) and trace element contents (B) normalized to chondrite (Sun and McDonough, 1989).
The ages (Figures 4 and12;Table 5) of magmatic and hydrothermal episodes the Tatatila-Las Minas deposits range between 16.34 and 13.92 Ma for the associated intrusive bodies (all of them observed as direct contributors to prograde skarn formation), and between 12.49 and 12.18 Ma for hydrothermal minerals (retrograde skarn stages). It is important to emphasize that the analyzed rocks are not merely terms of an intrusive suite that included IOCG skarn generators, but IOCG skarn generators themselves, as the sampling strategy was directed to rocks spatially associated with such mineralization-whether prograde or retrograde. The discussion to follow relies on this fact. The maximum time gap between prograde and retrograde skarn associations thus determined spans ~1.5 My, which is similar to that defined for other skarn deposits (i.e.,Camprubíet al., 2015). A late dyke that crosscuts the mineralization, in association with capping volcanic rocks of the Trans-Mexican Volcanic Belt, was dated at 4.11 Ma. The early Permian age obtained for intrusive rocks in the Las Minas area (286 ± 2 Ma) is likely to correspond to the Carboniferous-Permian arc (Ortega-Obregónet al., 2013;Kirschet al., 2012), known as the Teziutlán massif, that constitutes the basement in the region and was dated at 269-252 Ma (K-Ar;López-Infanzón, 1991) and at 281-268 Ma (40Ar/39Ar;Iriondoet al., 2003).
Figure 11
Isotope variation diagrams for the Tatatila-Las Minas Miocene intrusive bodies associated with IOCG skarn mineralization. (A) Sr-Nd isotopes variation diagram. (B) Pb-Nd isotopes variation diagram. (C) Pb isotopes variation diagram. (D) ɛNd vs. ɛSr diagram that illustrates possible end-member sources for magmas, afterDePaolo and Wasserburg (1979a,1979b). Key = DMM = depleted MORB-mantle, EMI = enriched mantle I, EMII = enriched mantle II, HIMU = mantle component, MORB = 5°-15° NE Pacific Rise mid-ocean ridge basalts, NHRL = northern hemisphere reference line, TMVB = current volcanic front of the Trans-Mexican Volcanic Belt. See sources for all reference values in Gómez-Tuenaet al. (2003), which is also the source of values represented as green dots in diagrams A to C that correspond to volcanic rocks from the Palma Sola area in the eastern TMVB. The magmatic fractionation and sediment recycling trends in the zoomed view of A are simplified afterHoffman and White (1982).
Figure 12
Summary of the U-Pb and39Ar/40Ar ages obtained in this study for the intrusive rocks and IOCG skarn mineralization at the Tatatila-Las Minas area, Veracruz.
A consistent range of ages between 24.60 and 19.04 Ma (late Oligocene to early Miocene;Figures 4 and12;Table 5) has been additionally obtained, which corresponds to intrusive rocks that predate the syn-mineralization suite. Such ages also predate the earliest stage of magmatism that is associated with the Trans-Mexican Volcanic Belt (Gómez-Tuenaet al., 2005,2007) and are similar to those characteristic of the final stage of magmatic activity of the Sierra Madre Occidental (Ferrariet al., 2005b,2007).
The multielemental and isotopic geochemical determinations of IOCG skarn-related intrusive rocks at Tatatila-Las Minas are sound and congruent indicators of mostly intermediate to acid (Figure 8A), metaluminous (Figure 8B), and I- and S-type rocks (Figure 9B) that were emplaced in a subduction-related continental arc (Figure 9A), and high La/Yb ratios could also be obtained through high pressures in basaltic melt (Figure 9C; McPhersonet al., 2006), since the late Oligoce to Miocene. In addition, these rocks are part of the medium- to high-potassium (not shown) calc-alkaline series, with adakitic signatures and a compelling isotopic affinity with the Trans-Mexican Volcanic Belt (TMVB). A sound adakitic affinity of most analyzed samples in the study area is determined by a general geochemical behavior (Tables 1 to3;Figure 9D,E) that meets most of the characteristics of such petrological association (Table 6). If anything, Y and Yb contents appear to be significantly higher than in adakitic (Tables 3 and6), a characteristic that will be addressed later on. Despite the possible occurrence of alkaline magmatism in the Palma Sola region in association with the Eastern Mexico Alkaline Province (EMAP;Demand and Robin, 1975;Negendanket al., 1985;Ferrariet al., 2005a), the formation of IOCG skarn deposits in the Tatatila-Las Minas district can be solely attributed to the TMVB, as no adakitic affinity has been consistently reported for the magmatism associated with the EMAP (see references inCamprubí, 2013). However, some ages of alkaline rocks in Palma Sola are much younger than syn-mineralization ages, with no associated mineralization. Then, the adakitic signatures found in the Palma Sola region are more likely to correspond to the volcanism of the TMVB rather than that of the EMAP. This is the first instance in which adakites are directly associated with the formation of any ore deposits in the TMVB-in this case, IOCG skarn deposits.
Comparative table between the general geochemical composition of adakites (as ofMoriet al., 2007;Richards and Kerrich, 2007) and of intrusive rocks at the Tatatila-Las Minas area.
However, anomalously high Y and Yb contents (with respect to typical adakitic signatures) similar to those found in the Tatatila-Las Minas host rocks have been explained in adakites as to reflect some degree of interaction with alkaline or ultrapotassic rocks (Luet al., 2013;Liuet al., 2017)-hence the high-potassium character of many of the studied rocks (?)-or due to crustal contamination (Zhanget al., 2017). Therefore, despite the likely dominant affinity of these rocks with the TMVB, some degree of interaction between their parental TMVB magmas and EMAP magmas cannot be ruled out at this stage of research. As a matter of fact, magmas with either affinity coexisted in the region, as evidenced by the formation of the Tatatila-Las Minas deposits (Negendanket al., 1985;Ferrariet al., 2005a; see alsoFigure 7 inCamprubí, 2009). Also, the occurrence of A-type granites (alkaline) is hinted at in some of the analyzed samples despite mostly belonging to I- and S-types (Figure 9B), but no affinity with within-plate granites was found (Figure 9A).
In addition, the data in this paper stand for the idea of a metallogeny of the TMVB in its own right, as established byCamprubí (2013). The ages of Miocene IOCG skarn-related magmatism in the Tatatila-Las Minas area (16.34 to 13.92 Ma) fit well within the ~19 to 10 Ma bracket defined byGómez-Tuenaet al. (2005,2007) for the early stages of the TMVB, particularly in its eastern region, in which the adakitic signature of volcanism is conspicuous. Such continental magmatism display geochemical signatures that strongly evoke those of adakites, with the inherent likeliness that it may be associated with melting of the flattened subducted slab (Gómez-Tuenaet al., 2005,2007;Moriet al., 2007). Adakite is the common term that refers to magmas produced by melting of subducted oceanic crust under high pressures and in the presence of water (due to dehydration of the subducted slab). However, other processes for magma generation are possible in the generation of magmas with adakitic geochemical signatures (Defantet al., 2002;Richards and Kerrich, 2007;Rodríguezet al., 2007;Richards, 2011;Maet al., 2015;Ribeiroet al., 2016;Denget al., 2017;Keevilet al., 2019). The adakitic signatures at a regional scale in the TMVB are the very high Sr/Y ratios, depletion in Y and HREE, and Sr, Nd and Pb isotopic compositions that approximate to those of mid-ocean ridge basalts in the East Pacific Rise (Gómez-Tuenaet al., 2005,2007;Moriet al., 2007). Nonetheless, adakitic affinities do not necessarily imply that these magmas are derived from the melting of the subducted slab alone, and other geological mechanisms are also plausible for their inception or as relevant contributors to adakitic signatures, as discussed below.
The linkage between adakitic magmas and the variety of tectonomagmatic settings that the generation of such magmas entails is suggestive of a significant potential for the formation of associated ore deposits (González-Partidaet al., 2003a,2003b;Chiaradiaet al., 2004;Sunet al., 2011;Denget al., 2017;Keevilet al., 2019). Although the association between “adakites” and ore deposits normally refers to the classic definition of adakite magmas, the generation of such magma through melting of a subducted slab has been questioned (Richards and Kerrich, 2007;Richards, 2011). In the case of Tatatila-Las Minas, however, the intrusive rocks of adakitic-affinity associated with IOCG skarn mineralization have dominantly high-silica compositions (Figure 9E). This denotes that melting of basalt from the subducted slab would have effectively occurred, with subsequent reaction of the resulting melts with peridotites during their ascent through the mantle wedge (Defant and Drummond, 1990;Drummond and Defant, 1990;Martinet al., 2005). Also, the distribution of Nd and Sr isotopic compositions in the Tatatila-Las Minas intrusions point to magma fractionation as per their distribution (Figure 11A). ɛNd values in the analyzed rocks (between -6.6 and 3.2;Table 3) point to contributions of both relatively isotopically enriched and depleted magma sources for Nd, and represent mantle derived melts that were contaminated by continental crust lithologies, especially when correlated with ɛSr values (Figure 11D). As already highlighted byGómez-Tuenaet al. (2003), Pb isotopic compositions lie between those expected for subducted sediments and MORB (Figure 11B,11C andTable 3), thus requiring an isotopically depleted source.
An association between adakites and the formation of IOCG skarn deposits was earlier established in Mexico for the late Cretaceous-early Paleocene Mezcala deposits in the Sierra Madre del Sur (Camprubí and González-Partida, 2017, and references therein). The formation of adakites in that locality has been linked to early stages of a subduction-related continental arc (González-Partidaet al., 2003b), a feature that is explained by the switch from subduction-related oceanic arcs to continental arcs in southern Mexico during the Late Cretaceous (Camprubí, 2013,2017). Besides the particular case of Mezcala, in these and the Tatatila-Las Minas deposits the formation of associated adakitic magmas can be explained by slab rollback or flattening subduction as younger portions of the subducted slab were being consumed (Morán-Zentenoet al., 1999;Ferrari and Rosas-Elguera, 1999;Gutscheret al., 2000;Gómez-Tuenaet al. 2003;Keppie and Morán-Zenteno, 2005). Also, in both regions similar associations of different magmatic-hydrothermal types of deposits (i.e., IOCG, sulfide skarns, metalliferous porphyries, epithermal deposits;Camprubí, 2013,2017) were produced. Such flattening of the subducted slab has been extensively documented along the entire Western Cordillera of North America and the Andes and explains the historical distribution of metallogenic provinces within them (Camprubí, 2017, and references within).
However, magmatic processes such as assimilation and fractional crystallization (AFC) or those occurring in melting-assimilation-storage-homogenization (MASH) zones in “normal” continental arc magmas may also account for adakitic compositions of intrusions in association with the subsequent formation of magmatic-hydrothermal ore deposits (Richards and Kerrich, 2007;Richards, 2011;Gatzoubaroset al., 2014;Lohmeieret al., 2019). fact, these processes can generate andesitic to dacitic differentiates with HREE-depleted normalized REE patterns, and high La/Yb and Sr/Y ratios (Feeley and Davison, 1994;Kayet al., 1999;Klepeiset al., 2003;Richards, 2011). However, AFC processes can be virtually ruled out as important contributors to the adakitic signal because Eu anomalies in this case are weak (Figure 10; seeChenet al., 2014). The absence of Eu anomalies would support the model byRichards (2011), as high water contents in typical adakitic rocks are characteristic of MASH zones. MASH interactions may involve partial melts of lower crustal rocks that may imprint high La/Yb and Sr/Y. Such signature is derived from high pressure fractionation in MASH zones with amphibole and garnet, which would produce high La/Yb ratios, and from the suppression of plagioclase fractionation due to high water content in the magmas, thus resulting in high Sr/Y ratios (see references inRichards, 2011). In lowf S2 and highf O2 conditions underneath “normal” continental arcs, MASH processes may induce the formation of IOCG deposits in intra-arc settings (Richards and Mumin, 2013), thus producing an alternative scenario for the association between adakite-like and IOCG deposits. With regard to slab flattening underneath a continental arc due to steep subduction,Richards and Mumin (2013) argued about scarce to nil associated magmatic activity or the migration of magmatism toward back-arc settings. Interestingly, slab flattening would cause the dehydration of the slab and the subsequent hydration of the lithosphere, which would be too cold to melt. However, once the slab re-steepened, the temperature of the hydrated lithosphere would rise in contact with the asthenosphere, generating the partial melting of sub-continental mantle and subsequent vigorous volcanic flare-ups, thus reactivating the formation of magmatic-hydrothermal ore deposits-among them, IOCG deposits (seeFigure 1 inRichards and Mumin, 2013). The formation of adakites in such specific settings and in association with magmatic-hydrothermal ore deposits has not been reported. However, the involvement of MASH-type processes in metallogeny has actually been invoked in the formation of continental-arc related magmatic-hydrothermal ore deposits nonetheless (Sunet al., 2011). The possibility of magma generation by MASH-type processes that followed re-steepening of the subducted slab with which the formation of IOCG deposits would be linked is particularly significant for the Tatatila-Las Minas case. Indeed, the formation of these deposits occurred during the late Miocene, once the subducted slab underneath the Trans-Mexican Volcanic Belt, in fact, re-steepened (seeFigure 13 inGómez-Tuenaet al., 2003).
In summary, the most likely settings for the formation of parental adakitic magmas to the IOCG skarn deposits at Tatatila-Las Minas would be (1) a “normal adakitic” slab-melt setting with some crustal contamination, or (2) MASH-related adakitic compositions. However, these settings do not necessarily have to be considered as mutually exclusive in the generation of adakites with associated magmatic-hydrothermal ore deposits (Chenet al., 2014;Sunet al., 2018). To our reckoning, these settings cannot be effectively discriminated given the current wealth of data from the Tatatila-Las Minas district. In addition, it is possible that TMVB calc-alkaline and EMAP alkaline magmas underwent some kind of interaction that produced the intrusive bodies with which the studied IOCG skarn deposits are associated. Interestingly, despite the common tectonomagmatic affinity of all the Cenozoic magmatic rocks, the only samples that show high Y and Yb contents are those whose ages correspond entirely to the initial stages of the TMVB (not those older than 19 Ma). This, again, stands for different magmatic processes-albeit slightly-between TMVB and pre-TMVB rocks.
The iron oxide-Cu-Au deposits at the Tatatila-Las Minas district (central Veracruz) are skarn-related deposits that belong to the IOCG family, and associated Au-rich epithermal deposits also occur in the area. U-Pb and40Ar/39Ar dating of these IOCG skarns yielded early to middle Miocene ages for prograde (16.34 to 13.92 Ma for the associated intrusive bodies) and retrograde (12.44 to 12.18 Ma for hydrothermal minerals) associations. Such ages and the geochemical affinity of host intrusive rocks (calc-alkaline to adakitic) that are directly involved in the formation of IOCG skarns match well with those previously established for the early stages of evolution of the Trans-Mexican Volcanic Belt (TMVB). A set of pre-TMVB Cenozoic rocks has been also dated between ~24.6 and 19 Ma.
The multi-elemental and isotopic geochemical study of IOCG skarn-related intrusive rocks determined that these are intermediate to acid, metaluminous, I- and S-type, medium- to high-potassium, typical calc-alkaline to adakitic rocks that are compatible with those expected for a continental volcanic arc such as the TMVB. Therefore, the studied deposits are likely to be ascribed to the metallogeny of the TMVB, which can be rightfully spoken of as an actual metallogenic province. Such a fact broadens the economic expectations of a province that has traditionally been overlooked by mineral exploration.
The prominent adakitic signal as found in the IOCG skarn-generating intrusive rocks has been regionally attributed to adakitic melts associated with flat subduction and the subsequent resteepening of the subducted slab-with independent evidence for crustal contamination. The results in this paper concur with such an interpretation. The general geochemical characteristics of these rocks, however, do not rule out the possibility that melting-assimilation-storage-homogenization (MASH) processes were involved in the generation of parental magmas. There are also hints that these magmas interacted with alkaline melts, which would likely be associated with the nearly contemporaneous EMAP. Only TMVB rocks display Y and Yb contents that would suggest such interaction-all other petrologic indicators suggest common characteristics for TMVB and pre-TMVB Cenozoic rocks. In both a adakitic and MASH scenarios, the most plausible stage at which the formation of IOCG skarn-associated magmas occurred would be once the flattened subducted slab re-steepened, thus allowing melting of either (or both) slab material or the hydrated lower lithosphere.
SC-2b biotite
J = 0.00381413 ± 0.00000572 Volume 39ArK = 1.073 x E-13 cm3 NPT Integrated Date = 14.45 ± 0.15 Ma Plateau age = 14.46 ± 0.15 Ma (2s, including J-error of .2%) MSWD = 1.3, probability=0.24 Includes 98.8% of the 39Ar steps 1 through 12 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 9 points Age = 14.32 ± 0.17 Ma Initial 40Ar/36Ar =285 ± 20 MSWD = 0.89 Probability = 0.51
SC-2a biotite
J = 0.00383360 ± 0.00000575 Volume 39ArK = 1.146 x E-13 cm3 NPT Integrated Date = 15.44 ± 0.16 Ma Plateau age = 15.43 ± 0.16 Ma (2s, including J-error of .2%) MSWD = 1.4, probability=0.16 Includes 99.58% of the 39Ar steps 3 through 13 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 11 points Age = 15.30 ± 0.19 Ma Initial 40Ar/36Ar =287 ± 17 MSWD = 1.4 Probability = 0.19
SC-2-b1 biotite
J = 0.00385940 ± 0.00000579 Volume 39ArK = 1.775 x E-13 cm3 NPT Integrated Date = 16.34 ± 0.20 Ma Plateau age = 16.34 ± 0.20 Ma (2s, including J-error of .2%) MSWD = 0.89, probability=0.55 Includes 99.97% of the 39Ar steps 2 through 13 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 13 points Age = 16.01 ± 0.22 Ma Initial 40Ar/36Ar =308 ± 11 MSWD = 0.76 Probability = 0.68
RV-2 feldspar
J = 0.00391100 ± 0.00000587 Volume 39ArK = 1.264 x E-13 cm3 NPT Integrated Date = 22.10 ± 0.45 Ma Plateau age = 20.67 ± 0.57 Ma (2s, including J-error of .2%) MSWD = 0.46, probability=0.90 Includes 86.6% of the 39Ar steps 2 through 11 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 11 points Age = 21.5 ± 1.1 Ma Initial 40Ar/36Ar =288.5 ± 8.2 MSWD = 0.37 Probability = 0.95
CR-1 biotite
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.369 x E-13 cm3 NPT Integrated Date = 25.10 ± 1.07 Ma Plateau age = 24.6 ± 1.1 Ma (2s, including J-error of .2%) MSWD = 1.02, probability=0.39 Includes 89.2% of the 39Ar steps 4 through 8 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 8 points Age = 21.0 ± 2.5 Ma Initial 40Ar/36Ar =335 ± 16 MSWD = 2.0 Probability = 0.06
CR-1 hornblende
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.079 x E-13 cm3 NPT Integrated Date = 26.06 ± 2.30 Ma Plateau age = 23.4 ± 2.5 Ma (2s, including J-error of .2%) MSWD = 1.19, probability=0.31 Includes 80.7% of the 39Ar steps 4 through 7 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 7 points Age = 22.1 ± 2.8 Ma Initial 40Ar/36Ar =332 ± 23 MSWD = 0.76 Probability = 0.958
CR-1 feldspar
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.565 x E-13 cm3 NPT Integrated Date = 20.10 ± 0.60 Ma Plateau age = 19.04 ± 0.69 Ma (2s, including J-error of .2%) MSWD = 0.91, probability=0.50 Includes 95.1% of the 39Ar steps 1 through 9 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 9 points Age = 19.01 ± 0.98 Ma Initial 40Ar/36Ar =296.0 ± 7.7 MSWD = 0.98 Probability = 0.44
BQ-1c feldspar
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.272 x E-13 cm3 NPT Integrated Date = 24.15 ± 0.67 Ma Plateau age = 22.12 ± 0.74 Ma (2s, including J-error of .2%) MSWD = 0.80, probability=0.52 Includes 73.2% of the 39Ar steps 5 through 9 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 5 points Age = 21.4 ± 2.9 Ma Initial 40Ar/36Ar =304 ± 35 MSWD = 1.03 Probability = 0.36
BQ-1b biotite
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.483 x E-13 cm3 NPT Integrated Date = 13.92 ± 0.22 Ma Plateau age = 13.92 ± 0.22 Ma (2s, including J-error of .2%) MSWD = 1.14, probability=0.33 Includes 100% of the 39Ar steps 1 through 10 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 10 points Age = 14.39 ± 0.26 Ma Initial 40Ar/36Ar =292.0 ± 7.5 MSWD = 0.96 Probability = 0.46
BQ-1b Biotite
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.264 x E-13 cm3 NPT Integrated Date = 14.13 ± 0.33 Ma Plateau age = 14.60 ± 0.34 Ma (2s, including J-error of .2%) MSWD = 0.099, probability=0.96 Includes 73% of the 39Ar steps 4 through 7 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 4 points Age = 14.33 ± 0.87 Ma Initial 40Ar/36Ar =307 ± 130 MSWD = 0.14 Probability = 0.87
FSC-1 Fuchsita
Integrated results
J = 0.002419 ± 0.000010 Plateau age tp = 12.49 ± 0.09 Ma Weighted mean of fractions e to i, representing 92.99% of39Ar released in 5 consecutive fractions, MSWD = 0.18 Isochron age tc = 12.45 ± 0.11 Ma; (40Ar/36Ar)i = 300 ± 15, MSWD = 0.2 for n = 5 (e to i)
Plateau age spectra and normal isochron diagrams in host intrusive bodies to the IOCG skam deposits at the Tatatila-Las Minas district, Veracruz.
data from zircons in the intrusive rocks associated with IOCG skarns in the Tatatila-Las Minas area.
# U and Th concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.† Isotopic ratios are corrected relative to 91500 standard zircon for mass bias and down-hole fractionation (91500 with an age ~1065 Ma; Wiedenbecket al, 1995). Isotopic207Pb/206Pb ratios, ages and errors are calculated following Patonet al. (2010).* All errors in isotopic ratios are in percentage whereas ages are reported in absolute and given at the 2-sigma level. The weighted mean206Pb/238U age is also reported in absolute values at the 2-sigma level. The uncertenties have been propagated following the methodology discussed by Patonet al. (2010).* *Rho is the error correlation value for the isotopic ratios206Pb/238U and207Pb/235U calculated by dividing these two percentage errors. The Rho value is required for plotting concordia diagrams.
trace element data for LA-ICPMS spot analyses on zircon grains for intrusive units in Tatatila de Las Minas, Veracruz, Mexico.
Element concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.
This paper constitutes a part of the dissertations of E.F.G. and G.H.A., who acknowledge the support of CONACyT through PhD and MSc grants, respectively. The Instituto de Geología UNAM is acknowledged for authorizing E.F.G. to carry on her PhD research along with her academic duties. Funding for this work was provided by CONACyT through the research grants 155662 to A.C. and “GEMex: Cooperación México-Europa para la investigación de sistemas geotérmicos mejorados y sistemas geotérmicos supercalientes” (within the 4.1 and 8.2 research sections:Determinación de propiedades petrológicas, de alteración hidrotermal, microtermométricas, geoquímicas, de isótopos estables y geocronológicos de afloramientos basamentales de áreas aledañas a Los Humeros y Acoculco, Pue.” to GEOMINCO S.A. de C.V. Additional funding was provided by the Instituto de Geología UNAM and the Centro de Geociencias UNAM through personal allocations. The radiogenic isotope determinations were carried out at the Centro de Geociencias UNAM with the assistance of Ofelia Pérez Arvizu and Carlos Ortega Obregón. The thin sections were elaborated by Juan Tomás Vázquez Ramírez of the Centro de Geociencias UNAM. FRX determinations were carried out at the Laboratorio Nacional de Geoquímica y Mineralogía-Instituto de Geología UNAM with the assistance of Rufino Lozano Santacruz. The separation of zircon crystals was carried out with the assistance of Teodoro Hernández Treviño of the Instituto de Geofísica UNAM. Assistance during field work was provided by Jesús Castro and Dunia Figueroa. Also,Figure 1 was drawn with the assistance of Rodrigo Delgado Sánchez. The authors are also grateful to Carl Nelson, Lisard Torró and Joaquín Proenza, the guest editors of the present special issue, and to three anonymous referees, whose comments helped to significantly improve this manuscript.
*Corresponding author: (A. Camprubí)camprubi@comunidad.unam.mx
Figure 1
Geological map of the Tatatila-Las Minas mining district, east of the Palma Sola massif. Adapted fromServicio Geológico Mexicano (2007,2010). Purple circles denote the location of samples on which this study is based, with indication of the obtained ages.
Figure 2
Selected aspects of the IOCG skarn mineralization at the Tatatila-Las Minas deposits showing both prograde (garnet and tourmaline) and retrograde (actinolite and fuchsite) associations. (A) Hand specimen showing a garnet-rich prograde association followed by an actinolite- and fuchsite-rich retrograde association in the Santa Cruz mine. (B) Photomicrography of a garnet and tourmaline prograde association followed by an actinolite and fuchsite retrograde association; transmitted light, crossed polars; same sample as in A. Fuchsite separates from A and B were dated by argon geochronometry in this study. (C) Hand specimen of prograde patchy to partially banded magnetite ore; El Dorado mine. (D) Hand specimen of banded exoskarn magnetite- and chalcopyrite-rich retrograde ore, with martitized magnetite; El Dorado mine. Key: Amp = amphibole-group minerals (actinolite), Cal = calcite, Ccp = chalcopyrite, Ep = epidote, Fuch = chromian muscovite or “fuchsite”, Grt = garnet-group minerals (grossular-andradite), Hm = hematite, Mag = magnetite, Mc = malachite, Py = pyrite, Tur = tourmaline.
Figure 3
Photomicrographs of representative hypabyssal bodies, unaffected by hydrothermal alteration, associated with IOCG skarn mineralization in the Tatatila-Las Minas district. (A) Quartz-monzodiorite showing euhedral apatite crystals within plagioclase phenocrysts, Santa Cruz mine; transmitted light, crossed polars. (B) Quartz-monzodiorite showing myrmekitic intergrowths, surrounded by hornblende and plagioclase phenocrysts, La Virgen mine; transmitted light, crossed polars. (C) Quartz-monzodiorite showing hornblende phenocrysts, Santa Cruz mine; transmitted light, crossed polars. (D) Quartz-monzodiorite showing euhedral zircon crystals within potassium feldspar, Santa Cruz mine; transmitted light, crossed polars. (E) Monzodiorite showing hornblende intergrown with magnetite, Carbonera mine; plane-polarized transmitted light. (F) Monzodiorite showing euhedral hornblende, biotite and apatite crystals within a plagioclase-potassium feldspar assemblage, Carbonera mine; plane-polarized transmitted light. (G) Monzogranite showing rock-forming biotite crystals, same sample as in F, Rancho La Virgen; transmitted light, crossed polars. (H) Monzogranite showing late biotite crystals intergrown with magnetite, Rancho La Virgen; plane-polarized transmitted light. Key: Ap = apatite, Bt = biotite, Fp = potassium feldspar, Hb = hornblende, Mt = magnetite, Pl = plagioclase, Qz = quartz, Zr = zircon.
Major elements in host intrusive rocks to the Tatatila-Las Minas IOCG deposits. All values in wt.%. Asterisks (*) correspond to analyses inDorantes-Castro (2016).
Key: LOI = loss on ignition.
Trace elements in host intrusive rocks to the Tatatila-Las Minas IOCG deposits. All values in ppm unless otherwide noted. Asterisks (*) correspond to analyses inDorantes-Castro (2016).
Sr, Nd and Pb isotopic values of selected samples from intrusive rocks associated with IOCG skarn mineralization in the Tatatila-Las Minas area.
Figure 4
Outlines of40Ar/39Ar age spectra (plateau ages) of intrusive host rocks to the IOCG skarn deposits in the Tatatila-Las Minas district, Veracruz.
Figure 5
40Ar/39Ar age spectra (plateau and isochron ages) of a chromian muscovite (“fuchsite”) from the magmatic-hydrothermal retrograde assemblage of the IOCG skarn deposit in the Santa Cruz mine, Tatatila-Las Minas district, Veracruz.
Figure 6
Tera-Wasserburg U-Pb concordia diagrams and plots of weighted averages of individual206Pb/238U ages of analyzed zircons, and pre-ablation SEM-CL images of zircons from a granodiorite intrusive from the Santa Cruz Mine (A), and from a potassic alteration assemblage that was pervasively developed on a granite-granodiorite intrusion in the village of Tatatila (B), from the Tatatila-Las Minas district, Veracruz. Solid-line ellipses, with blue square centers, are data used for age calculations; gray-line ellipses are data excluded from age calculations due to different degrees of Pb-loss and/or zircon inheritance. All U-Pb data are plotted with 2-sigma errors and all calculated weighted mean ages are also listed at the 2-sigma level. Original U(Th)-Pb data can be found for inspection inTable 5.
Figure 7
Tera-Wasserburg U-Pb concordia diagrams for zircons from various intrusive bodies in the Tatatila-Las Minas area. (A) Post-mineralization dyke. (B to D) Syn-mineralization hypabyssal bodies whose age can be attributed to the prograde skarn associations. (E) Granitic intrusive that corresponds to the Permo-Triassic basement. Solid-line ellipses, with black square centers, are data used for age calculations; gray-line ellipses are data excluded from age calculations due to different degrees of Pb-loss and/or zircon inheritance. All U-Pb data are plotted with 2-sigma errors. Original U(Th)-Pb data can be found for inspection inAppendix 3.
U-Th-Pb analytical data for LA-ICPMS spot analyses on zircon grains for granitic units in Tatatila de Las Minas, Veracruz, Mexico.
n = 30 Mean206Pb/238U Age = 12.18 ± 0.21 (2 sigma, MSWD = 1.5; n = 26)# U and Th concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.† Isotopic ratios are corrected relative to 91500 standard zircon for mass bias and down-hole fractionation (91500 with an age ~1065 Ma; Wiedenbeck et al., 1995). Isotopic 207Pb/206Pb ratios, ages and errors are calculated following Paton et al. (2010).* All errors in isotopic ratios are in percentage whereas ages are reported in absolute and given at the 2-sigma level. The weighted mean 206Pb/238U age is also reported in absolute values at the 2-sigma level. The uncertenties have been propagated following the methodology discussed by Paton et al. (2010).* *Rho is the error correlation value for the isotopic ratios 206Pb/238U and 207Pb/235U calculated by dividing these two percentage errors. The Rho value is required for plotting concordia diagrams.* **Percentage discordance values are obtained using the following equation (100*[(edad 207Pb/235U)-(edad 206Pb/238U)]/edad 207Pb/235U) proposed by Ludwig (2001). Positive and negative values indicate normal and inverse discordance, respectively. Individual zircon ages in bold were used to calculate the weighted mean 206Pb/238U age and MSWD (Mean Square of Weigthed Deviates) using the computacional program Isoplot (Ludwig , 2003).
Summary of geochronometric data obtained for host intrusive rocks and IOCG mineralization at the Tatatila-Las Minas area.
Figure 8
Petrological discrimination diagrams from major elements in intrusions associated with IOCG skarn mineralization in the Tatatila-Las Minas district, Veracruz. (A) Silica vs. alkaline element bivariant diagram, adapted fromCoxet al. (1979). (B) Alumina saturation diagram, adapted fromFrostet al. (2001), with compositions of skarns fromMeinert (1995). (C) AFM diagram, adapted fromIrvine and Baragar (1971).
Figure 9
Petrological discrimination diagrams from trace elements in intrusive rocks associated with IOCG skarn mineralization in the Tatatila-Las Minas district, Veracruz. (A) Y+Nbvs. Rb, Yvs. Nb, Ta+Ybvs. Rb, and Ybvs. Ta diagrams for discriminating tectonic settings, adapted fromPearceet al. (1984). (B) Discrimination diagram for different granite sources, adapted fromWhalenet al. (1987). (C) Discrimination diagram for the generation of magmas by fractional crystallizationvs. variable degree of partial melting, adapted fromThirlwallet al. (1994). (D) Discrimination diagram for adakitic affinity, adapted fromMartin (1986) with chondrite-normalized valuesSun and McDonough (1989). (E) Discrimination diagrams for high-silica (HSA) and low-silica adakites (LSA), adapted from Martin and Moyen (2002, 2003) and Martinet al. (2005). Key: HSA = high-silica adakites (>60% SiO2), LSA = low-silica adakites (<60% SiO2), ORG = ocean ridge granites, VAG = volcanic arc granites, syn-COLG = syn-collision granites, WPG = within plate granites.
Figure 10
Spider diagrams of REE (A) and trace element contents (B) normalized to chondrite (Sun and McDonough, 1989).
Figure 11
Isotope variation diagrams for the Tatatila-Las Minas Miocene intrusive bodies associated with IOCG skarn mineralization. (A) Sr-Nd isotopes variation diagram. (B) Pb-Nd isotopes variation diagram. (C) Pb isotopes variation diagram. (D) ɛNd vs. ɛSr diagram that illustrates possible end-member sources for magmas, afterDePaolo and Wasserburg (1979a,1979b). Key = DMM = depleted MORB-mantle, EMI = enriched mantle I, EMII = enriched mantle II, HIMU = mantle component, MORB = 5°-15° NE Pacific Rise mid-ocean ridge basalts, NHRL = northern hemisphere reference line, TMVB = current volcanic front of the Trans-Mexican Volcanic Belt. See sources for all reference values in Gómez-Tuenaet al. (2003), which is also the source of values represented as green dots in diagrams A to C that correspond to volcanic rocks from the Palma Sola area in the eastern TMVB. The magmatic fractionation and sediment recycling trends in the zoomed view of A are simplified afterHoffman and White (1982).
Figure 12
Summary of the U-Pb and39Ar/40Ar ages obtained in this study for the intrusive rocks and IOCG skarn mineralization at the Tatatila-Las Minas area, Veracruz.
Comparative table between the general geochemical composition of adakites (as ofMoriet al., 2007;Richards and Kerrich, 2007) and of intrusive rocks at the Tatatila-Las Minas area.
SC-2b biotite
J = 0.00381413 ± 0.00000572 Volume 39ArK = 1.073 x E-13 cm3 NPT Integrated Date = 14.45 ± 0.15 Ma Plateau age = 14.46 ± 0.15 Ma (2s, including J-error of .2%) MSWD = 1.3, probability=0.24 Includes 98.8% of the 39Ar steps 1 through 12 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 9 points Age = 14.32 ± 0.17 Ma Initial 40Ar/36Ar =285 ± 20 MSWD = 0.89 Probability = 0.51
SC-2a biotite
J = 0.00383360 ± 0.00000575 Volume 39ArK = 1.146 x E-13 cm3 NPT Integrated Date = 15.44 ± 0.16 Ma Plateau age = 15.43 ± 0.16 Ma (2s, including J-error of .2%) MSWD = 1.4, probability=0.16 Includes 99.58% of the 39Ar steps 3 through 13 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 11 points Age = 15.30 ± 0.19 Ma Initial 40Ar/36Ar =287 ± 17 MSWD = 1.4 Probability = 0.19
SC-2-b1 biotite
J = 0.00385940 ± 0.00000579 Volume 39ArK = 1.775 x E-13 cm3 NPT Integrated Date = 16.34 ± 0.20 Ma Plateau age = 16.34 ± 0.20 Ma (2s, including J-error of .2%) MSWD = 0.89, probability=0.55 Includes 99.97% of the 39Ar steps 2 through 13 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 13 points Age = 16.01 ± 0.22 Ma Initial 40Ar/36Ar =308 ± 11 MSWD = 0.76 Probability = 0.68
RV-2 feldspar
J = 0.00391100 ± 0.00000587 Volume 39ArK = 1.264 x E-13 cm3 NPT Integrated Date = 22.10 ± 0.45 Ma Plateau age = 20.67 ± 0.57 Ma (2s, including J-error of .2%) MSWD = 0.46, probability=0.90 Includes 86.6% of the 39Ar steps 2 through 11 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 11 points Age = 21.5 ± 1.1 Ma Initial 40Ar/36Ar =288.5 ± 8.2 MSWD = 0.37 Probability = 0.95
CR-1 biotite
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.369 x E-13 cm3 NPT Integrated Date = 25.10 ± 1.07 Ma Plateau age = 24.6 ± 1.1 Ma (2s, including J-error of .2%) MSWD = 1.02, probability=0.39 Includes 89.2% of the 39Ar steps 4 through 8 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 8 points Age = 21.0 ± 2.5 Ma Initial 40Ar/36Ar =335 ± 16 MSWD = 2.0 Probability = 0.06
CR-1 hornblende
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.079 x E-13 cm3 NPT Integrated Date = 26.06 ± 2.30 Ma Plateau age = 23.4 ± 2.5 Ma (2s, including J-error of .2%) MSWD = 1.19, probability=0.31 Includes 80.7% of the 39Ar steps 4 through 7 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 7 points Age = 22.1 ± 2.8 Ma Initial 40Ar/36Ar =332 ± 23 MSWD = 0.76 Probability = 0.958
CR-1 feldspar
J = 0.00393680 ± 0.00000591 Volume 39ArK = 0.565 x E-13 cm3 NPT Integrated Date = 20.10 ± 0.60 Ma Plateau age = 19.04 ± 0.69 Ma (2s, including J-error of .2%) MSWD = 0.91, probability=0.50 Includes 95.1% of the 39Ar steps 1 through 9 Inverse isochron (correlation age) results: Model 1 Solution (±95%-conf.) on 9 points Age = 19.01 ± 0.98 Ma Initial 40Ar/36Ar =296.0 ± 7.7 MSWD = 0.98 Probability = 0.44
BQ-1c feldspar
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.272 x E-13 cm3 NPT Integrated Date = 24.15 ± 0.67 Ma Plateau age = 22.12 ± 0.74 Ma (2s, including J-error of .2%) MSWD = 0.80, probability=0.52 Includes 73.2% of the 39Ar steps 5 through 9 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 5 points Age = 21.4 ± 2.9 Ma Initial 40Ar/36Ar =304 ± 35 MSWD = 1.03 Probability = 0.36
BQ-1b biotite
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.483 x E-13 cm3 NPT Integrated Date = 13.92 ± 0.22 Ma Plateau age = 13.92 ± 0.22 Ma (2s, including J-error of .2%) MSWD = 1.14, probability=0.33 Includes 100% of the 39Ar steps 1 through 10 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 10 points Age = 14.39 ± 0.26 Ma Initial 40Ar/36Ar =292.0 ± 7.5 MSWD = 0.96 Probability = 0.46
BQ-1b Biotite
J = 0.00401420 ± 0.00000602 Volume 39ArK = 0.264 x E-13 cm3 NPT Integrated Date = 14.13 ± 0.33 Ma Plateau age = 14.60 ± 0.34 Ma (2s, including J-error of .2%) MSWD = 0.099, probability=0.96 Includes 73% of the 39Ar steps 4 through 7 Inverse isochron (correlation age) results, plateau steps: Model 1 Solution (±95%-conf.) on 4 points Age = 14.33 ± 0.87 Ma Initial 40Ar/36Ar =307 ± 130 MSWD = 0.14 Probability = 0.87
FSC-1 Fuchsita
Integrated results
J = 0.002419 ± 0.000010 Plateau age tp = 12.49 ± 0.09 Ma Weighted mean of fractions e to i, representing 92.99% of39Ar released in 5 consecutive fractions, MSWD = 0.18 Isochron age tc = 12.45 ± 0.11 Ma; (40Ar/36Ar)i = 300 ± 15, MSWD = 0.2 for n = 5 (e to i)
Plateau age spectra and normal isochron diagrams in host intrusive bodies to the IOCG skam deposits at the Tatatila-Las Minas district, Veracruz.
data from zircons in the intrusive rocks associated with IOCG skarns in the Tatatila-Las Minas area.
# U and Th concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.† Isotopic ratios are corrected relative to 91500 standard zircon for mass bias and down-hole fractionation (91500 with an age ~1065 Ma; Wiedenbecket al, 1995). Isotopic207Pb/206Pb ratios, ages and errors are calculated following Patonet al. (2010).* All errors in isotopic ratios are in percentage whereas ages are reported in absolute and given at the 2-sigma level. The weighted mean206Pb/238U age is also reported in absolute values at the 2-sigma level. The uncertenties have been propagated following the methodology discussed by Patonet al. (2010).* *Rho is the error correlation value for the isotopic ratios206Pb/238U and207Pb/235U calculated by dividing these two percentage errors. The Rho value is required for plotting concordia diagrams.
trace element data for LA-ICPMS spot analyses on zircon grains for intrusive units in Tatatila de Las Minas, Veracruz, Mexico.
Element concentrations (ppm) are calculated relative to analyses of trace-element glass standard NIST 610.
