Abstract
The Montecristo district, northern Chile, is one of the few places worldwide where there is a direct relationship between magnetite-(apatite) (MtAp) mineralization and iron oxide-copper–gold (IOCG) mineralization. The MtAp mineralization includes Ti-poor magnetite, fluorapatite, and actinolite and is crosscut and partially replaced by a younger IOCG mineralization that includes a second generation of actinolite and magnetite with quartz, chalcopyrite, pyrite, and molybdenite. The MtAp stage at Montecristo is interpreted as the crystallized iron-rich melts that used the pre-existing structures of the Atacama Fault System as conduits. These rocks later acted as a trap for hydrothermal IOCG mineralization. Geochronology data at Montecristo indicate that the host diorite (U–Pb zircon 153.3 ± 1.8 Ma, 2-sigma), MtAp mineralization (40Ar-39Ar in actinolite, 154 ± 2 Ma and 153 ± 4 Ma, 2-sigma), and the IOCG event (Re-Os on molybdenite, 151.8 ± 0.6 Ma, 2-sigma) are coeval within error and took place in a time span of less than 3.4 Ma. The εHfi and εNdi values of the host diorite are + 8.0 to + 9.8 and + 4.3 to + 5.4, respectively. The whole-rock 87Sr/86Sri values of the IOCG mineralization (0.70425 to 0.70442) are in the lower end of those of the MtAp mineralization (0.70426–0.70629). In contrast, εNdi values for the IOCG mineralization (+ 5.4 and + 5.7) fall between those of the MtAp rocks (+ 6.6 to + 7.2) and the host diorite, which suggests that the IOCG event was related to fluids having a more crustal Nd (εNdi < + 5.7) composition than the MtAp mineralization. This likely reflects the mixing of Nd from the MtAp protolith and a deep magmatic-hydrothermal source, very likely an unexposed intrusion equivalent to the host diorite. Sulfur isotope compositions (δ34S, + 0.3 to + 3.4‰) are consistent with a magmatic source.
Similar content being viewed by others
Introduction
Iron oxide-copper–gold (IOCG) and magnetite-(apatite) (MtAp) mineralization are a controversial group of ore deposits that share some common characteristics, such as the abundance of iron oxides and are associated with calcic-iron-alkali alteration. The largest deposit on Earth is the world-class Olympic Dam, which is thought to represent the IOCG-MtAp system as a whole (Hitzman et al. 1992). IOCG sensu stricto mineralization is widely interpreted as hydrothermal and consists of magnetite and/or hematite and chalcopyrite, with subordinate pyrite, pyrrhotite, and bornite, and accompanied by actinolite and alkali feldspar in deep environments and chlorite-sericite-feldspar at shallow depths. IOCG deposits occur from the Late Archean (e.g., Carajás district in Brazil) to the Mesozoic (e.g., Chilean and Peruvian belts) (Hitzman et al. 1992; Williams et al. 2005), and most of these deposits are replacive on and show a space–time association with felsic and mafic igneous rocks, from intrusive bodies to volcanic rocks (Williams et al. 2005; Skirrow 2022).
Two contrasting models have been proposed to explain IOCG mineralization: 1) a magmatic-hydrothermal origin related to the exsolution of fluids derived from crystallizing intrusions (Vila et al. 1996; Vivallo and Henríquez 1997, 1998; Pollard 2000; Haynes 2000; Marschik and Fontboté 2001b; Tristá and Kojima 2003; Sillitoe 2003; Williams et al. 2005; Pollard 2006; Groves et al. 2010; Tornos et al. 2010; Barton et al. 2011; Richards and Mumin 2013); or, 2) a hydrothermal fluid derived from the interaction of connate brines with evaporite deposits or fluids derived from them in convective cells (Hitzman 2000; Barton and Johnson 1996, 2000).
Magnetite-(apatite) (MtAp) mineralization, also referred to as Kiruna-type (Geijer 1931), or iron oxide apatite (IOA), are massive low-Ti (< 1 wt.% Ti) magnetite ores with variable amounts of apatite, actinolite or clinopyroxene (dominantly diopside), scapolite, and anhydrite (Badham and Morton 1976; Nyström and Henríquez 1994; Naslund et al. 2002; Valley et al. 2010; Day et al. 2016; Tornos et al. 2017; Liu et al. 2018; Heidarian et al. 2018), with many of them containing anomalously high concentrations of REEs and U associated with the apatite (Hitzman 2000; Naslund et al. 2000; Valley et al. 2010). A key diagnostic feature of MtAp systems is the presence of large crystals of apatite and actinolite, sometimes showing pegmatitic textures (Tornos et al. 2021). MtAp deposits are known to occur from the Early Proterozoic (e.g., Kiruna, Sweden) (Geijer 1931) to the Pliocene (El Laco, Chile) (Park 1961), and are principally associated with calc-alkaline to alkaline igneous rocks (Naslund et al. 2000). They also include a pervasive metasomatic alteration of the host rocks similar to that of IOCG systems, with large zones of calcic-iron-alkali alteration (Hitzman 2000) that may grade upwards into hydrolytic alteration (Hitzman et al. 1992).
Different hypotheses have been proposed for the origin of MtAp deposits that include: 1) crystallization from an Fe-rich melt that separates from a more felsic melt (Badham and Morton 1976; Henríquez and Martin 1978; Frietsch 1978; Philpotts 1981, 1982; Nyström and Henríquez 1994; Naslund et al. 2002; Lledó and Jenkins 2008; Velasco et al. 2016; Nyström et al. 2016; Hou et al. 2018; Lledó et al. 2020); the crystallization is in many locations coeval with the exsolution of large amounts of hydrothermal fluids (Tornos et al. 2016, 2017); 2) circulation of hydrothermal fluids derived from the crystallization of mafic silicate melts that lead to the complete replacement of pre-existing rocks (usually andesite) by massive magnetite (Hildebrand 1986; Ménard 1995; Rhodes et al. 1999; Sillitoe and Burrows 2002; Gandhi 2003; Corriveau et al. 2016); 3) replacement of pre-existing rocks by heated hydrothermal fluids of non-magmatic origin such as basinal brines, or non-magmatic fluids equilibrated with evaporites (Hitzman 2000; Barton and Johnson 1996, 2000; Rhodes and Oreskes 1999); 4) the crystallization of magnetite microlites from a silicate melt followed by buoyant segregation from crystallizing mafic rocks, and flotation of the magmatic magnetite-bubble pairs, deposition of massive magnetite along faults, and posterior growth of hydrothermal magnetite (Knipping et al. 2015a); and, 5) crystallization from a sulfate and iron-rich magma being the product of the melting of shallow marine sediments by intruding andesite (Bain et al. 2020, 2021).
The relationship between both styles of mineralization remains a controversial topic. Early works (Hitzman et al 1992; Sillitoe 2003) included magnetite-(apatite) as part of the same system, while others have proposed that are part of the same clan (Williams et al. 2005) or end members on a continuous spectrum of mineralization (Hitzman et al. 1992; Vivallo and Henríquez 1997; Hitzman 2000; Gandhi 2003; Corriveau et al. 2016; Simon et al. 2018; del Real et al. 2021). This interpretation is based on the shared characteristics of IOCG and MtAp deposits, such as the similar mineral assemblage, related metasomatic alteration, and the local spatial association. This relationship is reinforced by the local presence of sparse hydrothermal late-stage Cu-Fe sulfides and gold overprinting some MtAp deposits or in their vicinity (e.g., Marcona in Peru, Chen et al. 2010a; and Cerro Negro Norte in Chile, Vivallo et al. 1995, Salazar et al. 2019), or the presence of mineralization ascribed to the IOCG group near clusters of MtAp deposits such as in the Norrbotten district, Sweden (Martinsson et al. 2016; Bauer et al. 2021) or the Great Bear Magmatic Zone, Canada (Ootes et al. 2017). However, on a global scale, MtAp systems are much more abundant than IOCG systems, and only locally do they coexist, such as in Carajás in Brazil (Xavier et al. 2012; Schutesky and de Oliveira 2020), Olympic Dam in Australia (Ehrig et al. 2012), or the Coastal Cordillera of the Andes (Espinoza 1990; Sillitoe 2003).
Other studies have suggested that MtAp and IOCG systems are genetically unconnected and form different ore systems, although they are sometimes superimposed. The detailed structural study of Bauer et al. (2018, 2021) in the Norrbotten region shows that in many cases, the IOCG mineralization can be significantly younger and unrelated to the MtAp mineralization, and the massive magnetite is just a chemical trap for the later Cu-Au-rich event.
In the Coastal Cordillera of the Andes, MtAp systems are dominantly vein-like, up to 50–100 m thick, with only a few extrusive tops, despite that locally, some stratabound replacive orebodies form when the veins crosscut rocks that are favorable for fluid-rock interactions (Espinoza 1990; Henríquez et al. 1994; Travisany et al. 1995; Tornos et al. 2021). Most well-known IOCG deposits in the Coastal Cordillera of the Andes are replacive in volcanoclastic andesite (e.g., Mina Justa, Peru, Chen et al. 2010b; Punta del Cobre district, Chile, Marschik and Fontboté 2001b, Arévalo et al. 2006, Del Real et al 2021; Dominga, Chile, Arredondo et al. 2017). However, in the northern part of the Coastal Cordillera in Chile, some smaller deposits occur as subvertical MtAp and IOCG veins in the plutonic rocks underlying the Jurassic-Cretaceous mafic volcanics such as the Montecristo, Tocopilla and Gatico districts (Boric et al. 1990; Sillitoe 2003). In these plutonic-hosted veins, Espinoza (1990) and Sillitoe (2003) have proposed a vertical evolution from deep MtAp to shallower IOCG mineralization both being channelized by large trans-crustal faults that also host diorite. A similar zonation model has been expanded by Simon et al. (2018) to the giant IOCG deposits of Mantoverde and Punta del Cobre in Chile. At the dominantly stratabound deposits in the Punta del Cobre district (which includes the Punta del Cobre and Candelaria deposits, among others), there are large bodies of apatite-poor massive magnetite and only sparse veins that could be unambiguously assigned to the MtAp style mineralization (Arévalo et al. 2006; Marschik and Fontboté. 2001b; Del Real et al. 2018; J. Carriedo, pers. com., 2020).
In the Carmen-Sierra Áspera district, which is very similar to Mantoverde (located ~ 40 km to the southeast), U–Pb dating of hematite in the Carmen de Cobre deposit suggests that the IOCG mineralization is ca. 10 My younger than the MtAp rocks of the Carmen de Fierro deposit (located ~ 5 km north of Carmen de Cobre) (Gelcich et al. 2005; Verdugo-Ihl et al. 2022). In the Marcona-Mina Justa district in western Peru, the MtAp ore and the related diopside-apatite pegmatite are at least 20 My older than the associated IOCG mineralization (Chen et al. 2011; Tornos et al. 2023) confirming that there is not a continuum between both styles of mineralization.
Therefore, one of the few places for elucidating the geochemical and geochronological relationship between unambiguous MtAp and IOCG systems are the vein-type deposits such as the Montecristo and Tocopilla districts, northern Chile, where the veins include both styles of mineralization. These veins represent a larger group of veins extending from northern Chile into southern Peru (Sillitoe 2003; Tornos et al. 2021). The absence of stratabound or shear-controlled replacive mineralization minimizes the effects of fluid-rock interaction in the mineralogy and geochemistry of the ores.
In the present study, the results from a detailed mineralogical, geochemical, and geochronological investigation of selected veins in the Montecristo district reveal that although the IOCG and MtAp deposits are in close proximity, and are broadly coeval, the geochemical data suggests that perhaps they are not part of the same ore forming system.
Geological background
The Andean magmatic arc started in the Jurassic as a consequence of the subduction of the oceanic lithosphere in the westernmost margin of the South American plate after the breakup of Gondwana, and it extended along the present-day Coastal Cordillera till nowadays (Dalziel et al. 1987; Mpodozis and Ramos 1990; Charrier et al. 2007). The oblique subduction led to the formation of a large strike-slip structure, the Atacama Fault System (AFS), which has a dominant sinistral component and extends along the Coastal Cordillera from ~ 20°30’ to ~ 29°45’S latitude (Hervé 1987; Mpodozis and Ramos 1990; Scheuber and Reutter 1992; Brown et al. 1993; Cembrano et al. 2005). North of the Bolivian orocline (Capitanio et al. 2011), equivalent structures have a WNW-ESE trend.
During the Late Jurassic-Early Cretaceous, the AFS experienced alternating both transtensional and transpressional movement (Cembrano et al. 1997, 2005; Grocott and Wilson 1997). These structures controlled the emplacement of the plutonic rocks and metasomatic alteration by hydrothermal fluids in the area (Grocott and Wilson 1997).
The Coastal Iron Belt (CIB) is situated along the Coastal Cordillera between 12º and 31ºS latitude and includes most of the IOCG and MtAp deposits in northern Chile and southern Peru, among others, such as porphyry copper deposits, Cu and Au-bearing veins, and “manto-type” copper deposits (Vivallo et al. 2000; Maksaev et al. 2007). The CIB is hosted by up to a 10-km thick regionally extensive sequence of Jurassic and Early Cretaceous subaerial basaltic andesite to andesite of the La Negra Formation and equivalent lithologies in southern Peru, represented by the Rio Grande Formation (Jaillard et al. 2000). Stratigraphically equivalent is the subaqueous Punta del Cobre Formation to the south, interpreted as deposited in an intra- to back-arc setting (Marschik and Fontboté 2001a). The La Negra Formation rocks have predominantly high-K to calc-alkaline affinities (Pichowiak et al. (1990); however, rocks of tholeiitic compositions extruded during the initial stages of the magmatic arc also exist (Pichowiak et al. 1990; Lucassen and Franz 1994). The volcanic sequence is intruded by Late Jurassic to Early Cretaceous plutonic rocks with similar geochemistry (Espinoza 1990). Together, they represent the magmatism that marks the onset of the Andean arc.
The Atacama Fault System has played a major role in the location of the IOCG and MtAp deposits (Boric et al. 1990; Freraut and Cuadra 1994; Vila et al. 1996; Gelcich et al. 1998; Sillitoe 2003; Vivallo et al. 2008; Tornos et al. 2021). Most of the MtAp mineralization is located on the main AFS, but IOCG mineralization is also related to N-S and NW–SE subsidiary structures, with both sinistral strike-slip and normal dip-slip extensional movements (Dallmeyer et al. 1996; Cembrano et al. 2005).
Systematic K–Ar and 40Ar-39Ar dating of hydrothermal silicates (e.g., actinolite, K-feldspar, biotite, white mica) indicate that most IOCG and MtAp deposits in the CIB were formed during two events of Middle to Late Jurassic (~ 170–150 Ma) and Early Cretaceous (~ 130–110 Ma) age, but there is also some IOCG mineralization of likely late Cretaceous and even Paleocene age (Boric et al. 1990; Vila et al. 1996; Vivallo and Henríquez 1998; Gelcich et al. 1998, 2003; Sillitoe 2003; Vivallo et al. 2008; Arredondo et al. 2017) which suggests that the ore forming events span more than 130 Ma.
The Montecristo district
The Montecristo mining district, the focus of the present study, is located in the westernmost Middle-Late Jurassic sub-belt of the CIB (Espinoza et al. 1996; Vivallo and Henríquez 1998; Sillitoe 2003) (Fig. 1a) where the AFS separates volcanic rocks to the east from plutonic rocks to the west. The mineralized veins in this district are hosted by the Matancilla Intrusive Complex (MIC) (Naranjo and Puig 1984; Escribano et al. 2013; Álvarez et al. 2016; Fig. 1a). The Complex includes multiple lithologies, including granodiorite, diorite, quartz diorite, and monzogranite and have been dated between 178 to 154 Ma (Mavor et al. 2022). Overall, these plutonic rocks are metaluminous and have tholeiitic to calc-alkaline, to high-K affinities (Álvarez et al. 2016). In detail, the host rock to the mineralization is a medium grained, inequigranular and unoriented amphibole-rich diorite.
The Montecristo district hosts at least twelve mineralized veins striking roughly NNE-SSW that are ~ 500 to ~ 5,000 m long and ~ 3 to ~ 30 m wide (Fig. 1b; Boric et al. 1990). The IOCG mineralization includes dominant magnetite and actinolite with lesser amounts of quartz, chalcopyrite, pyrite, bornite, molybdenite, and minor titanite. The ore includes small quantities of gold, and some veins show a geochemical anomaly enriched in Ag, Pb, Zn, and V (Espinoza et al. 1996). There is a large oxidation zone capping the vein and includes atacamite, chrysocolla, antlerite, and minor stringhamite [CaCuSiO4•2(H2O)] and chenevixite [Cu2Fe3+2(AsO4)2(OH)4]. A minor “cementation zone” between the oxidized and primary zones is composed of covellite and chalcocite (Espinoza et al. 1996).
The mineralized veins in the Montecristo district are located in the Paposo segment, one of the three parts of the AFS between Antofagasta and Paposo (Brown et al. 1993). In detail, they are structurally controlled by tensional faults oblique to the main AFS and formed during a sinistral strike slip event (Tornos et al. 2021).
The best-known deposit in this district, the Montecristo vein, is a magnetite-rich vein striking N45º-50ºE and dipping 75-80ºNW, with an average width of ~ 10 m but at some places reaching a thickness of ~ 30 m. The recognized vertical extension is about 400 m. The vein was discovered in 1850; since then, ~ 10 Mt ore has been mined. Despite the economic importance of the Montecristo vein, this study focuses on the currently mined Abundancia vein and the cropping out but unmined San Juan vein due to its key geological characteristics. The Abundancia vein contains both MtAp and IOCG mineralization and has been historically an important copper source in the region. The San Juan vein is different from the other veins in the Montecristo district since it is a classical MtAp vein composed of magnetite, apatite, and actinolite with no IOCG-like overprint.
The host diorite has been irregularly altered by hydrothermal fluids to an unknown extent. The original plagioclase is in some locations replaced from the center to the edge of the grains by potassium feldspar; this alteration is followed by a later, low temperature, pervasive replacement by adularia, albite, or white mica, and local late calcite. Magmatic amphibole has been strongly chloritized (ESM Fig. 1). The immediate halo of the veins shows a texturally destructive alteration up to some cm-thick that includes albite, quartz, magnetite, chlorite, actinolite and white mica (Vivallo and Henríquez 1998). The halo of pervasive alteration seems to increase in size and intensity downwards.
The Abundancia MtAp-IOCG vein
The Abundancia vein is a N40ºE and ~ 4 m width copper-rich (1.4% Cu) vein, with subordinate amounts of gold (0.4–0.5 g/t) and molybdenum (0.03%), and small concentrations of Co and Ni, which has been mined down to 940 m below the land surface. Several hypogene mineral assemblages are recognized, which are part of the following assemblages (Fig. 2): 1) an early MtAp assemblage (MtAp I) that includes magnetite, apatite, and actinolite predominantly arranged as cm-thick vertical subparallel bands. The pale green to dark green actinolite and the centimeter-sized crystals of white apatite can have comb-like unidirectional growth textures (Fig. 3a, b & e) growing perpendicular to the edge of the veins (Fig. 3a & b); 2) A second MtAp event (MtAp II) having the same mineral assemblage but with decreasing grain size and abundant massive magnetite occurring as crosscutting veins (Fig. 3a); 3) an IOCG stage that seems to dominate in the upper parts of the Abundancia vein and clearly crosscut and replaces the MtAp assemblages (Fig. 3a & b). Minerals in this stage include medium-grained magnetite, dark green actinolite, chalcopyrite, pyrite, quartz, titanite and gold; pyrite is replaced by later chalcopyrite; 4) a copper-rich event with up to ~ 20 cm thick massive veins with chalcopyrite and quartz (Fig. 3c); and 5) a late stage event with small veinlets of calcite or laminar molybdenite along fault surfaces, which cut both the earlier MtAp and IOCG mineralization.
In the deepest levels of the system (~ 940 to ~ 915 m), there is a brecciated monzogabbro which is partially replaced by an andradite-diopside-magnetite assemblage coeval with a halo of external albitization and has a retrograde replacement of pyroxene by actinolite and veinlets of actinolite and hydroxyapatite, with minor biotite, quartz, chalcopyrite, and pyrite (Díaz et al. 2018).
The San Juan MtAp vein
The San Juan MtAp vein, located approximately 1.5 km south of the Abundancia vein, has a NE-SW strike, ~ 4 m width, and ~ 600 m of inferred length. It contains massive magnetite with subvertical bands variably enriched in coarse-grained fluorapatite and actinolite (Fig. 2) that form pegmatite-like textures similar to those found in the Abundancia vein (Fig. 3d & e). Outcrops of this vein show a strong supergene alteration with replacement of the magnetite by hematite and late cross-cutting calcite veinlets.
Analytical methods
Petrography
Samples for petrographic studies were selected from the Abundancia and San Juan veins. Polished thin sections from the mineralization and the host rock were used for transmitted and reflected light microscopy, and for cathodoluminescence (CL) imaging on apatite to reveal the distribution and metasomatic alteration. Back-scattered electron (BSE) imaging was done on two 25 mm diameter, by 6 mm tall epoxy mounts containing several rock billets after polishing and carbon coating. The billets were examined with a JEOL 7100F field emission gun scanning electron microscope (FEG-SEM) at Memorial University of Newfoundland equipped with energy-dispersive spectroscopy (EDS) at 15 kV and 50 nA. The same epoxy mounts were used for further analyses by Electron Probe Microanalysis (EPMA) and Laser Ablation Inductively Coupled Plasma Mass Spectrometry (LA-ICPMS).
Electron probe microanalysis (EPMA)
The electron probe microanalyzer (EPMA) analyses were done using a JEOL JXA-8230 SuperProbe at Memorial University of Newfoundland. Details for the EPMA analyses are summarized in the ESM Table 2A for actinolite analyses, ESM Table 2C for magnetite, and ESM Table 2E, for apatite.
The EPMA was also used to produce X-ray maps of actinolite in the same areas where the traverses were done and analyzed for Ti, Ca, Fe, Mg, and K based on previous determinations of the actinolite chemical composition. Potassium was included in order to determine its distribution for later 40Ar-39Ar dating. For each X-ray map, an accelerating voltage of 15 kV, a sample current of 200 nA, and a pixel dwell time of 100 ms were used.
EPMA traverses and spot analyses were done on representative grains of actinolite, magnetite, and apatite for major and minor elements (Fig. 4). Analyses include Ti, Mn, K, Ca, S, Na, Al, Si, Mg, Fe, F, and Cl for the actinolite. Special attention was devoted to the analysis of potassium, which as mentioned above is critical for selecting samples suitable for later 40Ar-39Ar geochronology (see below). The magnetite EPMA analyses included Cu, Ni, Mn, Ti, V, K, Ca, S, P, Na, Al, Si, Mg, Zn, Co, Fe, Cr, and Ba, and apatite was analyzed for F, Fe, Mn, As, Ca, Sr, S, Na, Si, Mg, Cl, Y, and P. EPMA traverses were done across mineral grains avoiding inclusions, cracks, and obvious alteration zones.
LA-ICPMS trace elements of actinolite
The same epoxy mounts containing actinolite from the Abundancia and the San Juan veins were used for in situ laser ablation inductively coupled plasma-mass spectrometry (LA-ICPMS) trace elements analyses (Fig. 4a, c, d, f & g). The analyses were done on either the traverse or the spots where the EPMA analyses were previously done. A total of 43 trace elements were measured by LA-ICPMS, shown in the ESM Table 2A.
The actinolite trace element analyses were done at the Micro Analysis Facility (MAF) at Memorial University of Newfoundland, using a Thermo Finnigan Element XR high-resolution (HR), double-focusing magnetic sector-inductively coupled-plasma mass spectrometer (HR-ICPMS), coupled to a GeoLas 193-nm Excimer laser system. Details for the LA-ICPMS analyses are summarized in the ESM Table 2A. The Ca concentration in the actinolite determined by EPMA was used as the internal standard for calibrating each LA-ICPMS analysis. Analyses were done using a pulse frequency of 10 Hz, an energy density of 4 J/cm2, and a 40 µm diameter spot size. For each analysis, the background was measured for around 30 s, followed by 60 s of laser ablation. NIST 610 was used as the primary standard, and USGS BCR-2G and NIST 612 were used as secondary standards and were interspersed during the trace element analyses.
40 Ar/ 39 Ar geochronology in actinolite
The study includes two 40Ar/39Ar dates of actinolite, one fine-grained actinolite associated with the MtAp I event in the Abundancia vein (sample MOC-18–09), and one coarse-grained actinolite with moderate chloritization from the San Juan vein (sample MOC-18–05). Actinolite grains were crushed using a Plattner mortar-and-pestle; individual fragments were handpicked under a binocular microscope using tweezers, avoiding altered fragments and impurities, such as calcite. Both samples were subsequently crushed again using a Plattner mortar-and-pestle and sieved to a size of 90–150 µm to obtain a concentrate of pure actinolite.
Sample separates (15 and 18 mg, ESM Table 4) were irradiated in the McMaster University nuclear reactor (Ontario, Canada) for 30 MWh, carefully avoiding Cd shielding. The age monitor was hornblende MMhb1 (McClure Mountain Complex, Colorado) with an assumed age of 523.98 Ma (Schoene and Bowring 2006). 40Ar/39Ar analyses by step-heating following the procedures described by Bosio et al. (2020) were done on the Nu Instruments Noblesse noble gas mass spectrometer at the University of Milano-Bicocca, Milan, Italy, equipped with one Faraday cup with a 1011 Ω resistor and two ion counters. Before each sample was analyzed, a blank measurement was followed by two successive measurements of atmospheric Ar delivered by a pipette system. The analytical protocol consisted of four measurement cycles to allow an in-run cross-calibration of the relative collector gains, F/IC0 and IC0/IC1, in each individual run under the (variable) mass spectrometer conditions. The four cycles were: 40Ar (F), 38Ar (IC0), 36Ar (IC1); 39Ar (F), 37Ar (IC0), 35Cl (IC1); 41C3H5 (F), 39Ar (IC0), 37Ar (IC1); 38Ar (F), 36Ar (IC0). Each measurement consisted of 25 repetitions of the four cycles to avoid artefacts due to ion counter nonlinearities (Barberini and Villa 2015). The ion counter gains are similar, but not necessarily equal, to those determined from the atmospheric pipettes (which also independently quantify the mass fractionation originating in the source). The raw data from the mass spectrometer zero-time extrapolation were processed with an in-house Excel spreadsheet by correcting for mass spectrometer background, ion counter gains, blank measurements, source fractionation, and decay of 37Ar since irradiation, as well as propagating all the associated uncertainties. The total concentrations of 39Ar, 38Ar, and 37Ar were converted to absolute concentrations of K, Cl, and Ca, respectively, to control stoichiometry; and provide the Ca/K and Cl/K ratios. Since the stoichiometry of the present samples was determined by EPMA, for the age calculation, we consider only “isochemical steps” (Villa et al. 2006) with constant Ca/K and Cl/K ratios matching the EPMA measurements. The 40 K decay constant used for the age calculation was 5.543 × 10–4 Ma−1 (Steiger and Jäger 1977), even if it is very probably slightly inaccurate (Min et al. 2000; Naumenko-Dèzes et al. 2018) was used, by convention, in order to make intercomparisons with literature data less cumbersome. Different levels of uncertainty must also be considered. The internal concordance/discordance of step ages is assessed based on in-run uncertainties only. The inter comparison of samples analyzed in one, and the same laboratory requires compounding the in-run uncertainty with that on the neutron flux gradient. Finally, the inter comparison of the present samples with those reported in the literature requires inclusion of the systematic uncertainties on the age of the irradiation monitor (if the information is provided at all) and of the uncertainty of the 40 K decay constant, where necessary.
LASS U–Pb geochronology and Hf tracer isotopes on zircon
Representative zircon crystals from the Abundancia host diorite (sample MOC-18–02) were analyzed in situ for U–Pb, and Lu–Hf isotopes. The diorite was processed using a jaw crusher, and a dish mill afterwards, retaining > 90% of the crushed material between < 500 µm and > 63 µm. Magnetic separation and heavy liquids, bromoform and methylene iodide, separated and concentrated the zircon from the other minerals. The zircon grains were handpicked under the binocular microscope and mounted in epoxy. The mounts were polished to expose the crystal centers and carbon coated. Cathodoluminescence images were obtained using a JEOL 7100F field emission gun scanning electron microscope (FEG-SEM) at Memorial University of Newfoundland equipped with energy-dispersive spectroscopy (EDS) at 15 kV and 50 nA. These images were acquired in order to select the best zircon grains and areas for the in situ analyses, avoiding inclusions of other minerals, inherited cores, or cracks.
Lu–Hf isotopes paired with U–Pb isotopes were analyzed by laser ablation split stream (LASS) at the Micro Analysis Facility at Memorial University of Newfoundland. We followed the instrument configuration, operating parameters, and data reduction methods outlined by Fisher et al. (2011; 2014a, b) with the exception that N2 was added to the Ar carrier gas to both mass spectrometers for increased sensitivity.
The Lu–Hf analyses were done using a Thermo Finnigan Neptune multi-collector inductively coupled plasma mass spectrometer (MC-ICPMS). The U–Pb analyses were done simultaneously using a Thermo Finnigan Element XR high-resolution magnetic-sector inductively coupled plasma mass spectrometer (HR-ICPMS). Both mass spectrometers were interfaced to a GeoLas Pro 193 nm Ar-F excimer laser, operating at 10 Hz, 5 J/cm2, and a 40 μm diameter spot size. After collecting a background for 30 s, the samples were ablated for 60 s. The data were reduced using Iolite (Paton et al. 2011).
The Lu–Hf analyses done on the Montecristo MOC-18–02 zircon crystals were interspersed with analyses of several zircon reference materials to assess accuracy, mass bias corrections, and external reproducibility. While all of these reference materials listed below were analyzed during the LASS analyses, not all the reference materials were used for both the Hf-Lu and U–Pb data reduction.
The zircon reference materials used in this study for the Lu–Hf LASS analyses (determined in previous studies by solution MC-ICP-MS) covered the range of (Lu + Yb)/Hf of the MOC-18–02 zircon crystals included MUNZirc-1 (B-140) (176Hf/177Hf = 0.282135 ± 7; Fisher et al. 2011), MUNZirc-4 (B-144) (176Hf/177Hf = 0.282135 ± 7; Fisher et al. 2011), Plešovice (176Hf/177Hf = 0.282482 ± 13; Sláma et al. 2008), FC-1 (176Hf/177Hf = 0.282182 ± 14; Fisher et al. 2014a, b); and R33 (176Hf/177Hf = 0.282767 ± 18; Fisher et al. 2014a, b).
The zircon reference materials used in this study for the U–Pb LASS analyses included 91,500 (1065 Ma, Wiedenbeck et al. 2004), which was used as the primary U–Pb reference material, and Plešovice (337.13 Ma, Sláma et al. 2008); FC-1 (1099 Ma, Paces and Miller 1993), R-33 (419 Ma, Black et al. 2004); Temora 2 (417 Ma, Black et al. 2004), and 02,123 (295 Ma, Ketchum et al. 2001) were used as secondary reference materials.
Initial εHf values were calculated using a present-day CHUR value of 176Hf/177Hf = 0.282785, 176Lu/177Hf = 0.0336 as reported by Bouvier et al. (2008) and the 176Lu decay constant of 1.867 × 10–11 yr–1 from Söderlund et al. (2004).
The 238U counts in the zircon grains from sample MOC-18–02 were anomalously high in the first batch of analyses (e.g., individual analyses were done on 46 grains) due to the high U concentrations in the zircon (e.g., 920–7,330 ppm). The in situ U–Pb analyses showed significant discordance caused by Pb loss, radiation damage, or a combination of both factors. A second batch of zircon grains (27 grains, one analysis per grain) also from sample MOC-18–02 were chemically abraded, generally following the method in Mattinson (2005) prior to the LASS analyses for U–Pb and Hf isotopes. The zircon grains were annealed for 48 h at 1100 °C in the air in a fused silica crucible. This was followed by partial dissolution in concentrated (i.e., 49%) HF sealed in a Parr Bomb in an oven at 190 °C for 2 h (the Parr Bomb was put into an oven at pre-set at 190 °C, and then after two hours at 190 °C the Parr Bomb was removed from the oven and allowed to cool for one hour to room temperature before opening). This procedure removes radiation-damaged, altered, or metamict zones in zircon that cannot be restored with annealing, thus virtually eliminating secondary lead loss. This method often results in, for simple grains of one age, concordant or near-concordant analyses.
Preliminary U–Pb data after annealing, however, produced similar results to the first batch of LASS analyses, likely because the uranium content in the zircon grains was still anomalously high (e.g., 560 to 5,320 ppm). Thus, the zircon accumulated radiation damage dose (i.e., alpha decay events per milligram of zircon) was calculated for each zircon using the equation from Murakami et al. (1991) to eliminate the extensive radiation-damaged zircon grains, or regions in zircon grains, in an attempt to get a concordant age. After the radiation damage calculations, zircon grains with radiation damage in stages 2 and 3 (Nasdala et al. 2004 modified the damage stages initially reported by Murakami et al. in 1991) were removed, keeping the 12 zircon grains in stage 1 that were > 90% concordant for the final age calculation and used for the εHf calculation.
Re-Os geochronology in molybdenite
Molybdenite (MoS2) has very high Re/Os during crystallization; this mineral incorporates almost no Os (Stein et al. 2001). Therefore, 187Re-187Os dating of molybdenite uses the following simplified equation: 187Osmeasured = 187Remeasured × (eλt—1), where λ is the decay constant for 187Re (1.666 × 10–11 yr–1; Smoliar et al. 1996).
Re-Os isotopic analyses were done on two molybdenite samples from level 940 (m) of the Abundancia vein. There is no evidence of supergene alteration, since the petrography shows that the molybdenite is pure and has no mineral intergrowths. Molybdenite separates were obtained using a small hand-held drill and creating a molybdenite powder under a binocular microscope.
Re and Os isotopic concentrations were determined by isotope dilution using a Thermo Triton NTIMS (Negative Thermal Ionization Mass Spectrometry) machine at AIRIE. Precisely weighed samples were loaded into a Carius tube with 8 mL inverse aqua regia and Re-Os spikes for sample dissolution and sample-spike equilibration. A mixed Re-double Os spike (185Re-188Os-190Os) permits a mass fractionation correction for Os and assessment of any common Os in the molybdenite (Markey et al. 2003). Both samples had negligible common Os. All uncertainties are reported at 2-sigma (ESM Table 7) and include the 187Re decay constant uncertainty (λ). Re and Os blanks (Re blank = 11.77 ± 0.03 pg, Os blank = 0.130 ± 0.003 pg with 187Os/188Os = 0.350 ± 0.007) do not affect the calculated ages.
Whole rock Rb–Sr and Sm–Nd radiogenic isotopes
Three whole rock samples, rich in actinolite, from the IOCG stage of the Abundancia vein were analyzed for their bulk rock Sr and Nd radiogenic isotope composition. Samples were initially crushed using a Plattner mortar-and-pestle and later reduced to powder in ethanol with an agate mortar-and-pestle. The samples might contain some minor contamination from the MtAp stage as the fine-grained nature of the veinlets makes it sometimes difficult to separate minerals from both events.
Rb–Sr and Sm–Nd bulk rock analyses were done at the Unidad de Geocronología (CAI de Ciencias de la Tierra y Arqueometría) of the Universidad Complutense de Madrid using Isotope Dilution-Thermal Ionization Mass Spectrometry (ID-TIMS) with an Isotopx Phoenix TIMS. The samples were spiked with 84Sr, 87Rb, and a mix of 149Sm-150Nd, and digested using ultra clean reagents. They were processed using chromatography, where Rb, Sr, and rare earth elements (REEs), were separated in DOWEX AG (50WÅ ~ 12 Resin, 200–400 mesh) columns. In order to isolate the Sm and Nd, REE fractions were separated in HDEHP-impregnated Teflon-powder columns. Fractionation effects were corrected using a normalization of 86Sr/88Sr = 0.1194 and 146Nd/144Nd = 0.7219. The procedural blank was 0.5 ng for Sr, and 0.1 ng for Nd. The standard materials used during the analyses were the following: Sr standard NBS-987 (0.710245 ± 0.000004; n = 58), Nd standards La Jolla (0.511850 ± 0.000004, n = 36), and JNdi-1 (0.512108 ± 0.000003, n = 33).
δ 34 S stable isotopes in sulfides
The sulfur isotopic composition was determined in fourteen samples of chalcopyrite and pyrite from the IOCG event at the Abundancia vein. In four of them, two aliquots were taken for control. The sulfides were crushed in a stainless-steel mortar and handpicked under a binocular microscope using titanium tweezers, followed by magnetic separation to remove the magnetite. Pure sulfide grains were ground to ~1 micron size in an agate mortar using a pestle, obtaining a powder of ~0.5 g. Sulfides with other minerals attached to the surface were ground and sieved and retained to 63 microns, magnetite was removed with a magnet, and the sample was processed using heavy liquids, bromoform, and methylene iodide, in order to separate the sulfides from the other minerals.
The sulfur isotope measurements were done at the Stable Isotope Laboratory of the Instituto Andaluz de Ciencias de la Tierra (CSIC-UGR, Granada, Spain). Samples were analyzed by combusting the samples with V2O5 and O2 at 1030 °C in a Carlo Elba NC1500 elemental analyzer online with a Delta Plus XL mass spectrometer (EA-IRMS). The stable isotope composition is reported as δ values per mil, calculated using the δ = (Rsample/Rstandard – 1)*1000 equation, where R = 34S/32S for δ34S. Commercial SO2 was used as the internal standard for sulfur isotopic analyses. For sulfur, five internal standards (organic and inorganic material) ranging between − 6.38 ‰ to + 23.25 ‰ (CDT, Canyon Diablo Troilite), along with the IAEA international reference materials IAEA-S1, IAEA-S2, IAEA-S3, NBS-127 and CP-1 were analyzed. This study used three internal standards of + 23.25 ‰, + 6.03, and − 6.38 ‰ (CDT). After correction of the mass spectrometer daily drift, the precision calculated from standards systematically interspersed in analytical batches was better than ± 0.2‰. The V-CDT (Vienna-Canyon Diablo Troilite) is the standard reference material for reporting sulfur isotope data.
Results
Mineral textures and chemical composition
Actinolite
Our study identifies two types of actinolite in the Montecristo system: 1) medium to coarse grained, pale to dark green actinolite I, related to the MtAp I and II events; and 2) fine grained dark green actinolite II, associated with the later IOCG stage. Both types can be distinguished based on mineral textures, associated mineral assemblage, and crosscutting relationships.
Early actinolite I in the Abundancia vein occurs as subhedral, coarse-grained (up to 0.5 cm long) crystals, often arranged in cm-thick bands hosted by, or interlayered with, the magnetite. Actinolite I also occur disseminated in the magnetite or enclosed by apatite crystals in a “poikilitic” type texture (see Fig. 4b and ESM Fig. 1) that is similar to those in cumulate rocks (Wager et al. 1960; McBirney and Noyes 1979), and chilled parts of mafic and ultramafic intrusions (Wager 1961). Later crosscutting magnetite-rich mineralization MtAp II event also contains actinolite I with the same characteristics but with a smaller (up to 200 µm) grain size. Actinolite I in the San Juan MtAp vein is intergrown with magnetite, coarse-grained (up to 7 mm) and forms massive aggregates, some moderately altered to chlorite. Actinolite appears to predate at least some of the magnetite (Fig. 4g).
Actinolite II is fine-grained (up to 200 µm) and occurs in veinlets paragenetically related to chalcopyrite, pyrite, quartz, and minor magnetite and titanite of the IOCG event. There are also abundant remnants, partially replaced, of actinolite I and apatite in the IOCG zone (see ESM Fig. 1).
The chemical composition of the actinolite is in ESM Table 2B, and a comparative oxide composition is shown in Fig. 5b. All samples analyzed classify as actinolite (Fig. 5a; Leake et al. 1997), with the Mg# [molar (Mg/(Mg + Fe)] ranging from 0.67 to 0.79, and the Si a.p.f.u. between 7.796 and 7.910. In general, the actinolite samples show the following compositional range for major elements: MgO (15.20–18.70 wt.%), CaO (12.20–12.60 wt.%), FeO (8.70–13.10 wt.%), Al2O3 (0.90–1.70 wt.%), Na2O (0.10–0.20 wt.%), MnO (0.10–0.30 wt.%), TiO2 (0.01–0.04 wt.%), SO3 (0.01–0.02 wt.%), and small amounts of K2O (0.02–0.10 wt.%). It has between 190 and 940 ppm Cl and F is below the LOD.
The EPMA analyses show that actinolite I and II from Abundancia are relatively more depleted in Fe (range of 6.80–7.10 wt.%), while actinolite I from San Juan is enriched in Fe (average of 10.14 wt.%) and Mn (average of 2,562 ppm). Actinolite II in Abundancia has low average Ti (57 ppm), Cl (188 ppm), and K (197 ppm) compared to Actinolite I, which has an average of 112 ppm, 658 ppm, and 697 ppm, respectively. This is reinforced by the Fe, Ca, Mg, K, and Ti X-ray maps (ESM Fig. 3). Alternatively, the San Juan actinolite I is enriched throughout the grains in Ti (average of 253 ppm) compared to the other samples, but this can be due to the presence of microinclusions of a Ti–rich phase. In addition, the X-ray maps show the presence of abundant disseminated inclusions of titanite (< 10 µm) (ESM Fig. 3). The actinolite compositional data obtained in this study generally agree with other data previously published for actinolite in MtAp deposits in the Coastal Cordillera (Rojas et al. 2018) and IOCG (Del Real et al. 2021).
LA-ICPMS trace element data for actinolite are presented in the ESM Table 2B and Fig. 5c and are normalized to continental crust (Rudnick and Gao 2003). The data indicate that MtAp and IOCG-related actinolite, types I and II, respectively, from Abundancia generally have similar trace element patterns. Actinolite II has slightly higher Ba, Sr, and REE contents, and lower Rb, K, and Ti than actinolite I. Actinolite I from the San Juan vein, however, shows a different trace element composition. It is depleted in Nb, La, Ce, Nd, Sm, and Y and is slightly enriched in Hf, Zr, and Ti, suggesting different chemistry (i.e., heterogeneity) between individual veins.
Magnetite
Magnetite from the MtAp stage in the Abundancia and San Juan veins have similar textures and are grouped as magnetite type I. In some cases, the magnetite grains exhibit exsolution textures, less than 100 µm in size, of titanite (ESM Fig. 1v). Magnetite I from the Abundancia vein hosts abundant acicular actinolite inclusions (Fig. 4b). At the same time, the San Juan magnetite is inclusion-free. Magnetite II from the IOCG event is relatively homogeneous and has only quartz inclusions (See ESM Fig. 1). Major, minor, and trace elements data for both types I and II of the magnetite samples are presented in ESM Table 2D and comparative oxide average composition graphs are plotted in Fig. 5d.
All the magnetite I samples have a nearly identical chemical composition and do not show noticeable systematic trends or generations. The most representative elements include very low Ti contents (100–190 ppm) as well as Cu (50–60 ppm), Ni (200–240 ppm), V (1,420–3,010 ppm), Cr (30–60 ppm), S (10–50 ppm), Si (1,570–3,510 ppm), Al (650–940 ppm), Mn (100–180 ppm), Mg (340–1,610 ppm), Ca (710–1,320 ppm), and Zn (130–160 ppm). The relatively constant contents suggest that these elements are not hosted in nanoinclusions as happens in other MtAp Coastal Cordillera deposits (e.g., Knipping et al. 2015b). Chemical compositions similar to magnetite I are reported in magnetite from other MtAp deposits in the CIB (e.g., Knipping et al. 2015b, Los Colorados; Salazar et al. 2019, Cerro Negro Norte).
Magnetite II has a quite similar composition but significantly higher Cu contents, which are constant throughout the crystals. The average composition of magnetite II includes: Ti (282 ppm), Cu (406 ppm), Ni (150 ppm), V (6,193 ppm), Cr (81 ppm), S (162 ppm), Si (871 ppm), Al (730 ppm), Mn (149 ppm), Mg (211 ppm), Ca (251 ppm), and Zn (106 ppm). The geochemistry of magnetite in IOCG deposits is variable; however, similarities exist with previously published data on IOCG deposits worldwide (Rusk et al. 2009; Zhang et al. 2017).
Apatite
Apatite is abundant in the MtAp stage, but we have not found apatite in the IOCG stage. Apatite crystals in the Abundancia MtAp mineralization are generally euhedral, coarse-grained (up to 1 cm) and lack evidence in BSE or CL of any superimposed metasomatic alteration. The apatite commonly encloses small (< 250 µm) actinolite crystals in a “poikilitic-like” texture (Fig. 4b) and shows a yellow CL emission (caused by Mn2+) typical of apatite having a magmatic origin (Dempster et al. 2003; Bouzari et al. 2016) (see ESM Fig. 1). Apatite in the San Juan vein is also euhedral and coarse-grained (up to 2 cm) (Fig. 4h). The distribution of apatite is identical to that described in the nearby Carmen de Fierro and Fresia MtAp deposits (Tornos et al. 2021) and elsewhere in MtAp pegmatites.
The EPMA analyses of apatite are provided in ESM Table 2F, and the average chemical composition for each apatite sample is plotted in Fig. 5e. Our data are similar to published apatite composition data from other MtAp deposits in the Coastal Cordillera (Treloar and Colley 1996; La Cruz et al. 2019). All samples correspond to fluorapatite (2.10–2.50 wt.% F) with a minor chlorapatite component (0.40–0.80 wt.% Cl) and include Na (160–330 ppm), Mn (180–360 ppm), Mg (40–3,010 ppm), Fe (470–3,980 ppm), As (570–1,020 ppm), Sr (170–280 ppm), Y (140–200 ppm), Si (85–5,890 ppm), and S (50–110 ppm). However, the highly variable Fe, Si, and Mg contents are most likely due to the presence of micro- or nanoinclusions of actinolite that was not visible with BSE or EPMA.
Actinolite 40 Ar/ 39 Ar dates
40Ar/39Ar data and plots for the two actinolite samples are reported in Fig. 6, and the full data set is available in the ESM Table 4. Actinolite MOC-18–09 from the Abundancia vein yielded a well-constrained plateau age of 154 ± 2 Ma (2-sigma), representing 92.3% of the 39Ar released. The plateau age is indistinguishable from the isochron age (159 ± 17 Ma, 2-sigma). The three plateau steps show a constant Ca/K ratio. In contrast, the Ar-poor low-temperature steps show lower ratios, indicating the presence of a contaminant phase with a higher amount of K or lower Ca content. The broadly constant Cl/K ratios mean that the sample is homogeneous and likely contains little if any, Cl-rich phase such as fluid inclusions. The coarse-grained actinolite MOC-18–05 from the San Juan vein yielded a plateau age of 153 ± 4 Ma (2-sigma), representing the 76.2% of the 39Ar released, and homogeneous Ca/K and Cl/K ratios in the plateau steps. The plateau age is indistinguishable at the 2-sigma level from the isochron age (160 ± 4 Ma).
A literature search of K–Ar dates of actinolite associated with the alteration halo in the Mantos deposit in the Montecristo district has an average age of 164 ± 11 Ma (JICA-MMAJ 1986), an age consistent at the 1-sigma level with the present results.
Zircon U–Pb age and Hf tracer isotopes from the host rock
The U–Pb zircon data from the diorite host rock (sample MOC-18–02) are reported in the ESM Table 5 and in Fig. 7b to d. The Lu–Hf data are presented in Fig. 7e, and all the data and calculations are shown in the ESM Table 6. The zircon grains are of relatively small size (up to 100 μm long), euhedral to subhedral, and inclusion-free. Some grains show regular growth zoning, a characteristic of igneous zircon (Fig. 7a). No inherited cores were identified in any of the zircon grains.
The Concordia diagram shown in Fig. 7b yields a Concordia age of 153.3 ± 1.8 Ma (2-sigma; MSWD = 2.6; n = 12). This date overlaps with the weighted mean 206Pb/238U age of 152.1 ± 2.6 Ma (2-sigma; MSWD = 1.3; n = 12]), and with the inverse isochron (i.e., Tera-Wasserburg) date of 153.5 ± 3.0 Ma (2-sigma; MSWD = 15; n = 12) also obtained in this study. The Concordia age agrees with previously reported U–Pb dates for the earliest crystallization stages of the Matancilla Intrusive Complex (Álvarez et al. 2016; Mavor et al. 2022 and references therein). In situ initial Hf values from the above mentioned 12 individual zircon grains, and calculated using the Concordia age, range from εHfi of + 8.0 (analysis 11a) to an εHfi of + 9.7 (analysis 37a), which suggests a dominantly juvenile source, i.e., produced in the primitive mantle, but with some incorporation of less radiogenic Hf inherited from continental crust through assimilation during emplacement (εHfi = + 10.59 to + 16.56 for Mesozoic rocks of the Canadian Cordillera, Vervoort and Blichert-Toft 1999).
Molybdenite Re-Os ages
The Re-Os ages for two molybdenite samples belonging to the IOCG stage are presented in ESM Table 7. The data show a significant age difference between the two samples: 162.4 ± 0.6 Ma (2-sigma) (MOC-18-07A) and 151.8 ± 0.6 Ma (2-sigma) (MOC-18-07G). Re concentrations are very different, with 07A yielding 23.14 ppm and 07G yielding 369.7 ppm Re. The geologically unrealistic old age for sample MOC-18-07A could reflect within-molybdenite mobilization and redistribution of Re and Os or incipient oxidation (Stein et al. 2003). In an earlier study of a Chilean IOCG deposit at Raúl-Condestable (De Haller et al. 2006), two Re-Os ages several million years older than the associated intrusion were explained by Re loss. At Raúl-Condestable, however, molybdenite is markedly early in the paragenesis, preceding chalcopyrite and pyrite. In our study, molybdenite is late in the paragenesis (Fig. 2), and follows chalcopyrite ore mineralization, making a direct analogy difficult. Given the association of our dated samples with a zone of intense slickenside, we suggest that redistribution of Re and radiogenic Os within molybdenite is the most plausible explanation.
Sr and Nd radiogenic isotopes
Whole rock Rb–Sr and Sm–Nd results are shown in Fig. 8 and Table 1. The approximate deposit formation age of 154 Ma calculated by 40Ar/39Ar geochronology was used to calculate the initial 87Sr/86Sr and εNdi values. These calculations yield εNdi values of + 5.4 to + 7.0 and 87Sr/86Sri isotopic compositions of 0.70425 to 0.70442 for the IOCG event. This suggests a dominant juvenile source with values similar to the associated igneous rocks at the district scale. Other studies were previously done on calcite at Montecristo (87Sr/86Sri = ~ 0.7058), Julia (87Sr/86Sri = ~ 0.7046 to ~ 0.7048), and the Toldo-Velarde deposit in the Gatico district (87Sr/86Sri = ~ 0.7041 to ~ 0.7043), suggesting similar 87Sr/86Sri values, and thus a similar source (Vivallo and Henríquez 1998).
The εNdi value of + 7.0 is interpreted as inherited from the MtAp stage during the replacement of the MtAp by the IOCG mineralization. One MtAp value with a relatively low εNdi value of + 5.1 and high 87Sr/86Sri ratios (0.70629) from Tornos et al. (2021) is also here interpreted as being partially re-equilibrated with the late IOCG-related fluids (see below).
δ 34 S isotopes in sulfides
δ334S values of sulfides from the IOCG event are listed in Table 2 and shown in Fig. 9. Sulfur isotopes of fifteen chalcopyrite and pyrite aliquots show a very restricted range between + 0.3 to + 3.4 per mil, with a mean of + 1.8 per mil. δ34S values in the Abundancia IOCG event agree with previously published sulfur isotopic studies in the Montecristo deposit (0 to + 5‰, Vivallo and Henríquez 1998). In addition, they are similar to other sulfur isotopic compositions determined for IOCG deposits in the Coastal Cordillera of Chile (Tocopilla, − 0.2 to + 0.6‰, unpub. data, F. Tornos; Gatico district, 0 to + 5‰, Vivallo and Henríquez 1998; Naguayán-Desesperado, − 1 to + 1‰, Vivallo and Henríquez 1998; Candelaria, 0.3 to + 3.1‰, Marschik and Fontboté 2001b; and Julia, − 4 to + 4‰, Vivallo and Henríquez 1998) or the sparse sulfides in the MtAp mineralization such as Los Colorados (− 3.2 to + 2‰; Tornos et al. 2021) (Fig. 9). These data suggest that the sulfur is dominantly of a juvenile magmatic derivation (δ34S = -3 to + 2‰; Ohmoto and Rye 1979; Ohmoto and Goldhaber 1997), with only a minor contribution of either crustally contaminated igneous rocks or a sedimentary source (Poulson et al. 1991). However, hematite-rich IOCG deposits show a wider range of δ34S values with a systematic displacement towards positive values, indicating a likely input of sulfur derived from the abiogenic reduction of seawater sulfate, or microbial reduction in a closed system but external to the magmatic-hydrothermal system, such as in Raul-Condestable (+ 1.0 to + 26.3‰ main ore stage; De Haller and Fontboté 2009), Teresa de Colmo (− 5.5 to 18.2‰; Ledlie 1998), and Mantoverde (− 6.8 to + 11.2 main ore stage, and + 26.4 to + 36.2 in later stages; Benavides et al. 2007; Rieger et al. 2010). This is consistent with previous studies by Chen (2013) suggesting the importance of external sulfur with δ34S values > + 10‰ in the ore-forming processes of hematite-rich IOCG systems.
Discussion
Timing of mineralization
The U–Pb zircon crystallization age of 153.3 ± 1.8 Ma (2-sigma) for the host diorite is consistent within error with the 40Ar/39Ar dates of the MtAp event (153 ± 4 Ma and 154 ± 2 Ma, 2-sigma), and the Re-Os molybdenite formation age of 151.8 ± 0.6 Ma (2-sigma) for the IOCG mineralization. The dates obtained in the present study indicate that the crystallization of the host diorite, the formation of the MtAp rocks, and the hydrothermal IOCG mineralization, are coeval within error and took place in a period of ~ 3.4 Ma that cannot be temporally resolved with geochronology; probably, the mineralization formed soon after the crystallization of the host diorite, especially considering that the Ar cooling ages are younger than zircon dates (Renne et al. 1998; Min et al. 2000; Kuiper et al. 2008).
The age of the host diorite indicates that it is related to one of the youngest pulses of the Matancilla Igneous Complex and suggests that mineralization took place during the waning stages of magmatic activity.
The evolution of the Montecristo system
The Montecristo vein system consists of multi-stage mineralization with an early MtAp stage (MtAp I and II) crosscut by a younger IOCG mineralization event (Events III and IV) (see Fig. 2).
The textures of the early MtAp stage are similar to those of the intrusive part of the El Laco deposit or the MtAp deposits of the Coastal Cordillera that include the diagnostic pegmatite-like textures with sometimes banded, coarse-grained fluorapatite with unidirectional growth textures intergrown with diopside or actinolite and magnetite (Naranjo et al. 2010; Tornos et al. 2017, 2021). In the case of the Abundancia vein, they also include “poikilitic” fluorapatite enclosing well-formed acicular, mostly euhedral actinolite crystals. Some of these rocks likely represent the magmatic to hydrothermal transition in MtAp systems and a pulsatile regime similar to granite-related pegmatite (Tornos et al. 2021).
The direct crystallization of immiscible melts can explain the formation of MtAp mineralization. Geological (Naslund et al. 2002; Chen et al. 2010a; Mungall et al. 2018), isotopic (Tornos et al. 2016, 2017, 2021; Troll et al. 2019; Weis et al 2021), experimental (Lledó 2005; Veksler et al. 2007; Hou et al. 2018; Mungall et al. 2018; Lledó et al. 2020, between many others), numerical modeling (Keller et al. 2022), and melt inclusion data (Jacobsen et al. 2005; Naslund et al. 2009; Velasco et al. 2016; Bain et al. 2020, 2021 Pietruszka et al. in 2022) provide overwhelming evidence, often times neglected by supporters of alternative models, that these rocks crystallized from iron-rich silica-poor immiscible melts.
Overall, the composition of magnetite and actinolite in the IOCG stage (Event III) is similar to that of the MtAp stage. The actinolite and magnetite EPMA X-ray maps and compositions (ESM Fig. 3) show only slight chemical differences, suggesting that magnetite and actinolite compositions are not useful tracers for discriminating types of mineralization, at least within these vein-like MtAp and IOCG systems. The only observed major difference is that magnetite II is distinctive in containing abundant quartz inclusions not present in magnetite I. Furthermore, magnetite II has, on average, an order of magnitude more Cu (avg. 406 ppm) than magnetite I (50–64 ppm), something consistent with the abundance of Cu during the IOCG event.
These results suggest that the Montecristo MtAp-IOCG system has a much simpler mineralogical evolution of magnetite and actinolite than that described by, for example, Knipping et al. (2015b), Rojas et al. (2018), Huang and Beaudoin (2019), Palma et al. (2020) in other MtAp and IOCG deposits of the Coastal Cordillera. Our preferred interpretation is that the monotonous compositions at Montecristo are due to the fact that the mineralization is paragenetically simple and infills structures with little, if any, interaction with the host rocks.
For example, the Los Colorados MtAp system has very different types of ore, including massive and pegmatite infilling structures, as well as large amounts of disseminated and stockwork-like magnetite ore replacing andesite (see Knipping et al 2015b; Tornos et al 2021). The stratabound IOCG deposits of the CIB also include very different types of magnetite from the replacement of different types of andesitic rocks or direct precipitation along structures (del Real et al. 2021). Magnetite, in both cases, can have abundant micro- and nanoinclusions and/or records dramatic changes in fluid compositions, fO2 and temperatures, all influencing to an unknown extent in an unbuffered system the chemical composition of the minerals.
Different sources for the IOCG and MtAp mineralization
Except for one sample by Tornos et al. (2021), the Montecristo mineralization has rather restricted 87Sr/86Sri values (0.7042–0.7047). However, variable εNdi compositions (+ 5.4 to + 7.5) that plot in the mantle array of Fig. 8 are within the field of Late Jurassic-Early Cretaceous igneous rocks of the Coastal Cordillera. Specifically, the IOCG mineralization has lower (+ 5.4 to + 7.0), more crustal values than the MtAp rocks (+ 6.7 to + 7.5). The ultimate reason why sample ABU-7 of Tornos et al. (2021) has such a high radiogenic 87Sr/86Sri value is unknown but could reflect some inheritance from an old crustal basement or Sr-enriched slab fluids. Extremely radiogenic Sr values are a constant in some MtAp deposits of the CIB.
Nd isotopes are more robust and more difficult to reset than Sr isotopes and are a more reliable tracer of the source of fluids. Initial εNd values of the MtAp mineralization in Montecristo at 154 Ma are similar to those of nearby MtAp rocks, such as the Julia and Tocopilla vein deposits, and also of Jurassic age, that are + 7.5 and + 6.6. to + 6.8, respectively (Tornos et al. 2021). The lower εNdi values of the IOCG mineralization cannot solely be derived from the same source as the MtAp mineralization or fluids in equilibrium. The εNdi values of the Cu-Au mineralization probably inherited some Nd from the replaced MtAp rocks mixed with Nd derived from a reservoir with a significantly larger contribution of continental crust.
The εNdi values for the host diorite calculated from the εHfi of the zircon following Vervoot et al. (2011), assuming that there is no inheritance and that the zircon Hf is representative of the bulk rock, are between + 4.3 and + 5.4 (see ESM Table 6). Thus, our best interpretation is that the εNdi values of the IOCG likely track mixing between the replaced MtAp mineralization and the Nd inherited from the host diorite or a non-exposed underlying intrusion of equivalent isotopic composition. The relatively low εNdi values of the IOCG mineralization at Montecristo are broadly similar to those of the MtAp deposits located between Taltal and La Serena (Chile), as well as to the εNdi values of several intrusions in the Coastal Batholith (Lucassen et al. 2006). In particular, the εNdi values are like those of the Copiapó Plutonic Complex (CPC) of Early Cretaceous age (+ 4.7 to + 5.8; Marschik et al. 2003) that is thought to be the source of the ore forming fluids of the IOCG mineralization in the Punta del Cobre district deposits (Marschik and Fontboté 2001b; Arévalo et al. 2006).
The sulfur isotope results in this study from both Events III and IV in the IOCG mineralization are consistent with the sulfur being derived from an H2S-bearing magmatic-hydrothermal fluid that transported the Cu, Mo, and Au (i.e., Rusk et al. 2004; Audétat et al. 2008, 2019; Kamenetsky and Kamenetsky 2010); derivation of reduced sulfur from the hydrothermal leaching of juvenile igneous rocks seems unlikely. These fluids also likely precipitated magnetite II, actinolite II, quartz and sulfides by cooling and interacting fluids with the older MtAp mineralized rocks.
Thus, the MtAp rocks only acted as structural and geochemical traps for the later IOCG mineralizing fluids due to the brittle and oxidized nature of the magnetite in a scheme similar to that proposed by Bauer et al (2018) in the Malmberget-Kiruna area. Such fluid-rock interactions would also create the calcic-iron-alkali alteration observed in the deeper parts of the Montecristo system.
In the Coastal Cordillera, the formation of iron-rich melts during the Middle-Late Jurassic seems to be controlled by the contamination of the mantle wedge by Sr-rich crustally derived dehydration fluids. This would promote partial melting of the mantle wedge, triggering the separation of the iron-rich melts in the most contaminated zones from a parental mafic melt (Tornos et al. 2021). After the emplacement and crystallization of the MIC, these iron-rich melts ascended along restricted extensional domains within the transcrustal faults of the AFS and when attaining neutral buoyancy crystallized magnetite I, actinolite I, and fluorapatite at ca. 800–1200° C (Fig. 10; Bain et al. 2021 and references therein) coeval with dewatering. High temperatures and the oxidized nature of these iron-rich melts—systematically enriched in anhydrite (Tornos et al. 2017, 2021)—and probably the low Cu contents manifested in the magnetite composition inhibit the formation of significant amounts of sulfides. The small amounts of available reduced sulfur allowed the precipitation at temperatures below ca. 700 °C of the sparse amounts of sulfides found in MtAp systems. The homogeneous structures in the MtAp veins and the lack of explosion breccias suggest that these veins crystallized below the two-phase surface at ca. 0.2 GPa fluid pressure.
Geologic features and Lu–Hf, Sr–Nd, and S isotope compositions are consistent with the source of the IOCG mineralization event at Montecristo being magmatic-hydrothermal fluids derived from more crustal, but likely, primitive intrusions of likely dioritic composition. Here, upwelling fluids derived from the crystallization of the diorite would react with the earlier MtAp mineralization, promoting the formation of large amounts of magnetite II, actinolite II, quartz, and sulfides. Chalcopyrite and pyrite would precipitate due to the alkalinization and oxidation of fluids due to a reaction with magnetite I, a classical mechanism for destabilizing the chloride complexes that control the transport of Cu at high temperatures (Liu and McPhail 2005). A drop in the aH2Saq due to sulfide precipitation would, in turn, destabilize the HS complexes that control gold solubility (Shenberger and Barnes 1989; Loucks and Mavrogenes 1999; Zezin et al. 2007; Pokrovski et al. 2014), leading to its precipitation. Later collapse of the hydrothermal system should lead to the formation of the late calcite veins. Again, we have not found any evidence suggesting that the IOCG mineralization precipitated above the two-phase surface; thus, these veins are also probably deep.
In this scenario, the brecciated monzogabbro found underneath the Abundancia vein and recorded by Díaz et al. (2018) may represent the upper part of a crystallizing subvolcanic cupola that would be the deep magmatic-hydrothermal root of the IOCG vein-like mineralization. Equivalent breccias related to Cu-Mo-(Au) mineralization hosted by diorite have previously been described by Tornos et al. (2010) a few km from the Montecristo district in the Tropezón IOCG deposit. There, breccia pipes supported by tourmaline or quartz grade into a replacive IOCG-like mineralization with magnetite, sulfides, and actinolite, replacing mafic plutonic rocks. In addition, a direct link between the IOCG mineralization and the nearby silicate intrusions has previously been proposed for other IOCG deposits in the Coastal Cordillera (Boric et al. 1990; Vila et al. 1996; Hopper and Correa 2000; Marschik and Fontboté 2001b; Ray and Dick 2002; Gelcich et al. 2003; Sillitoe 2003; Tornos 2011).
Erosion prevents our knowledge of the mineralization that may have originally overlain these vein-like deposits. However, despite the age difference, these veins likely represent the roots of the abundant stratabound IOCG deposits in the Coastal Cordillera, such as in the Punta del Cobre district. The mineralization replaces favorable horizons in the Late Jurassic-Early Cretaceous andesite (Marschik and Fontboté 2001b). However, it is rooted in large NNW-SSE to WNW-ESE tensional sub-vertical structures infilled with an assemblage similar to the IOCG stage at Montecristo. The late calcite veins are identical to the uppermost part of the feeder structures in the Punta del Cobre district (N. Pop, pers. com., 2015).
The data obtained in this research might lead to different interpretations. However, based on our petrographic, geochemical, mineralogical, geochronological stable and radiogenic isotope, and field data, these data suggest that these systems reveal the superposition of genetically unrelated MtAp and IOCG systems that were channelized along deep-rooted tensional structures. The existence of independent sources for MtAp and IOCG systems most likely explains why, at a global scale, they are only rarely juxtaposed. Their similar mineralogy also can guide the misidentification of these mineralized systems in other locations. However, the large size of some of these IOCG deposits is likely due to the fact that magnetite is an excellent geochemical trap, but the copper-rich event can be several million years younger than the associated magnetite-rich rock (Rotherham 1997; Bauer et al. 2018).
Conclusions
Montecristo is one of the few places worldwide with a clear superposition of an IOCG mineralization on an older MtAp system. The Montecristo vein system is hosted by a very late pulse of diorite in the late Jurassic Matancilla Intrusive Complex. It includes the early formation at ca. 154 Ma of a MtAp assemblage, including magnetite, fluorapatite, and actinolite, having relatively high εNdi values. These rocks crystallized from iron-rich melts intruding along secondary structures of the transcrustal Atacama Fault System. This early magmatic event was most likely later followed by the circulation of magmatic-hydrothermal fluids along the same secondary structures, and the partial replacement of the earlier mineralization by an IOCG assemblage that precipitated from magmatic-hydrothermal fluids derived from the crystallization and dewatering of igneous rock with more crustal contamination, and likely equivalent to the host diorite. All mineralization events took place in a relatively short period of time, less than ~ 3.4 million years.
Although the genetic relationship between MtAp and IOCG deposits remains debatable, the petrographic, geochemical, mineralogical, geochronological stable and radiogenic isotope, and field data results from the present study suggest that the IOCG and MtAp mineralization events in the Montecristo district are genetically independent of each other.
Data availability
All data are available in the Supplementary Materials.
References
Álvarez J, Jorquera R, Miralles C, Padel M, Martínez P (2016) Cartas Punta Posallaves y Sierra Vicuña Mackenna. Carta Geol Chile Ser Geol Básica 183, 184:147 (1 map)
Arévalo C, Grocott J, Martin W, Pringle M, Taylor G (2006) Structural setting of the Candelaria Fe oxide Cu-Au deposit, Chilean Andes (27°30’S). Econ Geol 101:819–841
Arredondo C, Moscoso R, Prieto X, Carrasco R, Ortega R, Vivallo W, Mateo L, Pantoja G, Ulloa M, Ercilla O, Ridelle E (2017) Depósitos minerales de la región de Coquimbo. Servicio Nac Geol Minería Carta Geol Chile Serie Recur Minerales Energéticos 35:134 (2 maps)
Audétat A, Pettke T, Heinrich CA, Bodnar R (2008) The composition of magmatic-hydrothermal fluids in barren versus mineralizaed intrusions. Econ Geol 103:1–30
Badham JPN, Morton RD (1976) Magnetite–apatite intrusions and calc-alkaline magmatism, Camsell River, N.W.T. Can J Earth Sci 13:348–354. https://doi.org/10.1139/e76-037
Bain WM, Steele-MacInnis M, Li K, Li L, Mazdab FK, Marsh EE (2020) A fundamental role of carbonate-sulfate melts in the formation of iron oxide-apatite deposits. Nature Geosci 13:751–757. https://doi.org/10.1038/s41561-020-0635-9
Bain WM, Steele-MacInnis M, Tornos F, Hanchar JM, Creaser EC, Pietruszka DK (2021) Evidence for iron-rich sulfate melt during magnetite(-apatite) mineralization at El Laco, Chile. Geology 49:1044–1048
Barberini V, Villa IM (2015) Ion counter nonlinearities. In: Proceeding 25th Goldschmidt Conference, Prague: A4058
Barton MD, Johnson DA (2000) Alternative brine sources for Fe-oxide(-Cu-Au) systems: Implications for hydrothermal alteration and metals. In: Porter TM (ed), Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PGC Publishing, Adelaide, pp 43–60
Barton MD, Kreiner DC, Jensen EP, Girardi JD (2011) Multiple igneous-related hydrothermal systems and related IOCG mineralization, near Copiapó, Chile. In: Proceedings 11th SGA Biennial Meeting, Antofagasta, pp 521–523
Barton MD, Johnson DA (1996) Evaporitic source model for igneous related Fe oxide-(REE-Cu-Au-U) mineralization. Geology 24:259–262
Bauer TE, Andersson JBH, Sarlus Z, Lund C, Kearny T (2018) Structural controls on the setting, shape, and hydrothermal alteration of Malmberget iron oxide-apatite deposit, Northern Sweden. Econ Geol 113:377–395
Bauer TE, Lynch EP, Sarlus Z, Drejing-Carroll D, Martinsson O, Metzger N, Wanhainen C (2021) Structural controls on iron oxide copper-gold mineralization and related alteration in a Paleoproterozoic supracrustal belt: insights from the Nautanen deformation zone and surroundings, Northern Sweden. Econ Geol. https://doi.org/10.5382/econgeo.4862
Benavides J, Kyser TK, Clark AH (2007) The Mantoverde iron oxide-copper-gold district, III region, Chile: the role of regionally derived, nonmagmatic fluids in chalcopyrite mineralization. Econ Geol 102:415–440
Black LP, Kamo SL, Allen CM, Davis DW, Aleinikoff JN, Valley JW, Mundil R, Campbell IH, Korsch RJ, Williams IS, Foudoulis C (2004) Improved Pb-206/U-218 microprobe geochronology by the monitoring of a trace-element-related matrix effect; SHRIMP, ID-TIMS, ELA-ICP-MS and oxygen isotope documentation for a series of zircon standards. Chem Geol 205:115–140
Boric P, Díaz F, Maksaev J (1990) Geología y yacimientos metalíferos de la Región de Antofagasta. Boletín Servicio Nacional de Geología y Minería 40:1–246
Bosio G, Malinverno E, Villa IM, Di Celma C, Gariboldi K, Gioncada A, Barberini V, Urbina M, Bianucci G (2020) Tephrochronology and chronostratigraphy of the Miocene Chilcatay and Pisco formations (East Pisco Basin, Peru). Newslet Stratigr 53:213–247
Bouvier A, Vervoort JD, Patchett PJ (2008) The Lu-Hf and Sm-Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of the terrestrial planets. Earth Planet Sci Let 273:48–57
Bouzari F, Hart CJR, Bissig T, Barker S (2016) Hydrothermal alteration revealed by apatite luminescense and chemistry: A potential indicator mineral for exploring covered porphyry copper deposits. Econ Geol 111:1397–1410
Brown M, Díaz F, Grocott J (1993) Displacement history of the Atacama Fault System, 25º00’ to 27º00’S, Northern Chile. Geol Soc Amer Bull 105:1165–1174
Capitanio FA, Faccenna C, Zlotnik S, Stegman DR (2011) Subduction dynamics and the origin of Andean orogeny and the Bolivian orocline. Nature 480:83–86. https://doi.org/10.1038/nature10596
Cembrano J, Gonzalez G, Arancibia G, Ahumada I, Olivares V, Herrera V (2005) Fault zone development and strain partitioning in an extensional strike-slip duplex: a case study from the Mesozoic Atacama fault system, Northern Chile. Tectonophysics 400:105–125
Cembrano J, Lavenu A, Arancibia G, Sanhueza A, Reynolds P (1997) Coeval transpressional and transtensional magmatic arc tectonics in the southern Andes. In: Proceedings VIII Congreso Geológico Chileno, Concepción, pp 1613–1616
Charrier R, Pinto L, Rodríguez MP (2007) Tectonostratigraphic evolution of the Andean Orogen in Chile. In: Moreno T, Gibbons W (eds) The geology of Chile, Geol Soc, London, pp 21–114
Chen H (2013) External sulphur in IOCG mineralization: Implications on definition and classification of the IOCG clan. Ore Geol Rev 51:74–78
Chen H, Kyser TK, Clark AH (2011) Contrasting fluids and reservoirs in the contiguous Marcona and Mina Justa iron oxide-Cu (-Ag-Au) deposits, south-central Peru. Miner Deposita 46:677–706
Chen H, Clark AH, Kyser TK (2010) The Marcona Magnetite Deposit, Ica, South-Central Peru: A Product of Hydrous, Iron Oxide-Rich Melts? Econ Geol 105:1441–1456
Chen H, Clark AH, Kyser TK, Ullrich TD, Baxter R, Chen Y, Moody TC (2010) Evolution of the giant Marcona-Mina Justa iron oxide-copper-gold district, South-Central Peru. Econ Geol 105:155–185
Corriveau L, Montreuil JF, Potter E (2016) Alteration facies linkages among iron oxide copper-gold, iron oxide-apatite, and affiliated deposits in the Great Bear Magmatic Zone, Northwest Territories, Canada. Econ Geol 111:2045–2072
Dallmeyer D, Brown M, Grocott J, Taylor G, Treloar P (1996) Mesozoic magmatic and tectonic events within the Andean plate boundary zone, 26º-27º30’, northern Chile: constraints from 40Ar/39Ar mineral ages. J Geol 104:19–40
Dalziel IWD, Storey BC, Garrett SW, Grunow AM, Herrod LDB, Pankhurst RJ (1987) Extensional tectonics and the fragmentation of Gondwanaland. In: Coward MP, Dewey JF, Hancock PL (eds) Continental Extensional Tectonics. Geol Soc London Spec Publ 28:433–441
Day WC, Slack JF, Ayuso R, Seeger CM (2016) Regional Geologic and Petrologic Framework for Iron Oxide ± Apatite ± Rare Earth Element and Iron Oxide Copper-Gold Deposits of the Mesoproterozoic St. Francois Mountains Terrane, Southeast Missouri, USA. Econ Geol 111:1825–1858
De Haller A, Fontboté L (2009) The Raúl-Condestable Iron Oxide Copper-Gold Deposit, Central Coast of Peru: Ore and Related Hydrothermal Alteration, Sulfur Isotopes, and Thermodynamic Constraints. Econ Geol 104:365–384
De Haller A, Corfu F, Fontboté L, Schaltegger U, Barra F, Chiaradia M, Frank M, Zúniga Alvarado J (2006) Geology, Geochronology, and Hf and Pb Isotope Data of the Raúl-Condestable Iron Oxide-Copper-Gold Deposit, Central Coast of Peru. Econ Geol 101:281–310
Del Real I, Thompson JFH, Carriedo J (2018) Lithological and structural controls on the genesis of the Candelaria-Punta del Cobre iron oxide copper gold district, northern Chile. Ore Geol Rev 102:106–153
Del Real I, Reich M, Simon AC, Deditius A, Barra F, Rodríguez-Mustafa MA, Thompson JFH, Roberts MP (2021) Formation of giant iron oxide-copper-gold deposits by superimposed, episodic hydrothermal pulses. Commun Earth Environ 2:192
Dempster TJ, Jolivet M, Tubrett MN, Braithwaite CJR (2003) Magmatic zoning in apatite: a monitor of porosity and permeability change in granites. Contrib Miner Petrol 145:568–577
Díaz RA, Villela D, Fonseca E (2018) Presencia de diópsido-andradita en Mina Montecristo y su relación con la génesis de depósitos de Fe-P y del tipo IOCG presentes la Cordillera de la Costa, Región de Antofagasta, Chile. In: Proceedings XV Congreso Geológico Chileno Concepción pp 434
Ehrig K, McPhie J, Kamenetsky VS (2012) Geology and mineralogical zonation of the Olympic Dam iron oxide Cu-U-Au-Ag deposit, South Australia. In: Hedenquist JW, Harris M, Camus F (eds) Geology and genesis of major copper deposits and districts of the world: A tribute to Richard H. Sillitoe. Soc Econ Geol, pp 237–267
Escribano J, Martínez P, Domagala J, Padel M, Espinoza M, Jorquera R, Contreras JP, De La Cruz R, Calderón M (2013) Cartas Bahía Isla Blanca y Taltal, Región de Antofagasta. Servicio Nacional de Geología y Minería, Carta Geológica de Chile, Serie Geología Básica 164–165:72 pp, 1 map
Espinoza S, Véliz H, Esquivel J, Arias J, Moraga A (1996) The cupriferous province of the Coastal Range, Northern Chile. In: Camus F, Sillitoe RH, Petersen R (eds) Andean copper deposits: new discoveries, mineralization, styles and metallogeny. Econ Geol Spec Publ 5:19–32
Espinoza S (1990) The Atacama-Coquimbo ferriferous belt, northern Chile In: Fontboté L, Amstutz GC, Cardozo M, Cedillo E, Frutos J (eds) Spec Publ SGA. Springer Verlag, Berlin, pp 353-364
Fisher CM, Hanchar JM, Dhuime B, Samson SD, Blichert-Toft J, Vervoort J, Lam R (2011) Synthetic zircon doped with hafnium and rare earth elements for use as reference material for hafnium isotopic analyses. Chem Geol 286:32–47. https://doi.org/10.1016/j.chemgeo.2011.04.013
Fisher CM, Vervoort JD, Hanchar JM (2014) Guidelines for reporting zircon Hf isotopic data by LA-MC-ICPMS and potential pitfalls in the interpretation of these data. Chem Geol 363:125–133. https://doi.org/10.1016/j.chemgeo.2013.10.019
Fisher CM, Vervoort JD, DuFrane SA (2014a) Accurate Hf isotope determinations of complex zircons using the “laser ablation split stream” method. Geochem Geophys Geosyst 15:121–139. https://doi.org/10.1002/2013GC004962
Fox KA (2000) Fe oxide (Cu-Au) mineralization and alteration at the Productora prospect, Chile. MSc Thesis, Colorado School of Mines
Freraut R, Cuadra W (1994) Mineralización de Fe, Cu y Au en la Franja Cretácica de la Costa, Regiones III y IV de Chile. In: Proceedings VII Congreso Geológico Chileno, Concepción, pp 1046–1050
Frietsch R (1978) On the magmatic origin of iron ores of the Kiruna type. Econ Geol 73:478–485. https://doi.org/10.2113/gsecongeo.73.4.478
Gandhi SS (2003) An overview of the Fe oxide-Cu-Au deposits and related deposit types. Canadian Institute of Mining Technical Paper, CD-ROM
Geijer P (1931) The iron ores of the Kiruna type. Geol Surv Sweden 367:39
Gelcich S, Espinoza C, Vivallo W (1998) Yacimientos metalíferos de las hojas Chañaral y Diego de Almagro, Región de Atacama Chile. Servicio Nac Geol Mineria Mapa Recur Minerales Chile 3:15 (1 map)
Gelcich S, Davis DW, Spooner ETC (2005) Testing the apatite-magnetite geochronometer: U-Pb and 40Ar/39Ar geochronology of plutonic rocks, massive magnetite-apatite tabular bodies and IOCG mineralization in Northern Chile. Geochim Cosmochim Acta 69:3367–3384
Gelcich S, Davis D, Spooner ETC (2003) New U-Pb ages for host rocks, mineralization and alteration of Iron Oxide (Cu-Au) deposits in the Coastal Cordillera of Northern Chile. In: Proceedings 4th South American Symposium on Isotope Geology, Brazil, pp 63–66
Grocott J, Wilson J (1997) Ascent and emplacement of granitic plutonic complexes in subduction-related extensional environments. Mineral Soc Ser 8:173–195
Groves DI, Bierlein FP, Meinert LD, Hitzman MW (2010) Iron oxide copper-Gold (IOCG) deposits through Earth history: implications for origin, lithospheric setting, and distinction from other epigenetic iron oxide deposits. Econ Geol 105:641–654
Haynes D (2000) Iron oxide Copper (-Gold) Deposits: Their position in the ore deposit spectrum and modes of origin. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PCG Publishing, Adelaide, pp 71–90
Heidarian H, Lentz DR, Alirezaei S, McFarlane CRM, Peighambari S (2018) Multiple stage ore formation in the Chadormalu iron deposit, Bafq Metallogenic Province, Central Iran: evidence from BSE imaging and apatite EPMA and LA-ICP-MS U-Pb geochronology. Minerals 8:87
Henríquez F, Dobbs F, Espinoza S, Nyström J, Travisany V, Vivallo W (1994) Origin of Chilean magnetite-apatite ore deposits. In: Proceedings VII Congreso Geológico Chileno, Concepción, pp 822–824
Henríquez F, Martin R (1978) Crystal-growth textures in magnetite flows and feeder dykes, El Laco, Chile. Can Mineral 16:581–589
Hervé M (1987) Movimiento normal de la falla Paposo, zona de falla Atacama, en el Mioceno, Chile. Rev Geol Chile 31:31–36
Hildebrand RS (1986) Kiruna-type deposits: Their origin and relationship to intermediate subvolcanic plutons in the Great Bear magmatic zone, Northwest Canada. Econ Geol 81:640–659
Hitzman MW (2000) Iron oxide-Cu-Au deposits: what, where, when and why. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PGC Publishing, Adelaide, pp 9–25
Hitzman MW, Oreskes N, Einaudi M (1992) Geological characteristics and tectonics setting of Proterozoic iron oxide (Cu-U-Au-REE) deposits. Precambrian Res 58:241–287
Hopper D, Correa A (2000) The Panulcillo and Teresa de Colmo copper deposits: two contrasting examples of Fe-Ox Cu-Au mineralisation from the Coastal Cordillera of Chile. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PCG Publishing, Adelaide, pp 177–189
Hou T, Charlier B, Holtz F, Veksler I, Zhang Z, Thomas R, Namur O (2018) Immiscible hydrous Fe–Ca–P melt and the origin of iron oxide-apatite ore deposits. Nat Comm 9:1415
Huang XW, Beaudoin G (2019) Textures and chemical compositions of magnetite from iron oxide copper-gold (IOCG) and Kiruna-type iron oxide-apatite deposits and their implications for ore genesis and magnetite classification schemes. Econ Geol 114:953–979
Jacobsen JK, Veksler IV, Tegner C, Brooks CK (2005) Immiscible iron- and silica-rich melts in basalt petrogenesis documented in the Skaergaard intrusion. Geology 33:885–888
Jaillard E, Hérail G, Monfret T, Díaz-Martínez E, Baby P, Lavenu A, Dumont JF (2000) Tectonic evolution of the Andes of Ecuador, Peru, Bolivia and northernmost Chile. In: Cordani UG, Milani EJ, Thomaz Filho A, Campos DA (eds) Tectonic evolution of South America. Sociedad Brasileira de Geologia, pp 481–559
JICA-MMAJ (1986) Informe de estudios básicos sobre la exploración de recursos minerales del área sur de Antofagasta de la República de Chile: fase I. Japan International Cooperation Agency-Metal Mining Agency of Japan
Kamenetsky V, Kamenetsky M (2010) Magmatic fluids immiscible with silicate melts: examples from inclusions in phenocrysts and glasses, and implications for magma evolution and metal transport. Geofluids 10:293–311
Keller T, Tornos F, Hanchar JM, Pietruszka DK, Soldati A, Dingwell DB, Suckale J (2022) Genetic model of the El Laco magnetite-apatite deposits by extrusion of iron-rich melt. Nat Commun 13:6114
Ketchum JWF, Jackson SE, Culshaw NG, Barr SM (2001) Depositional and tectonic setting of the Paleoproterozoic Lower Aillik Group, Makkovik Province, Canada: evolution of a passive margin-foredeep sequence based on petrochemistry and U-Pb (TIMS and LA-ICP-MS). Precambrian Res 105:331–356
Knipping JL, Bilenker LD, Simon AC, Reich M, Barra F, Deditius AP, Lundstrom C, Bindeman I, Munizaga R (2015) Giant Kiruna-type deposits form by efficient flotation of magmatic magnetite suspensions. Geology 43:591–594
Knipping JL, Bilenker LD, Simon AC, Reich M, Barra F, Deditius AP, Wälle M, Heinrich CA, Holtz F, Munizaga R (2015) Trace elements in magnetite from massive iron oxide-apatite deposits indicate a combined formation by igneous and magmatic-hydrothermal processes. Geochim Cosmochim Acta 171:15–38
Kovacic Villarroel P, Barra F, Tornos F, Morata D, Cerda A (2012) Nuevos antecedentes geológicos y geoquímicos del yacimiento tipo IOCG Casualidad, Distrito Sierra Overa, II Región de Antofagasta, Chile. In: Proceedings XIII Congreso Geológico Chileno, Antofagasta, pp 43–45
Kuiper KF, Deino A, Hilgen FJ, Krijgsman W, Renne PR, Wijbrans JR (2008) Synchronizing rock clocks of Earth History. Science 320:500–504
La Cruz NL, Simon AC, Wolf AS, Reich M, Barra F, Gagnon JE (2019) The geochemistry of apatite from the Los Colorados iron oxide–apatite deposit, Chile: implications for ore genesis. Miner Deposita 54:1143–1156
Leake BE, Woolley Leake BE, Woolley AR, Arps CES, Birch WD, Gilbert MC, Grice JD, Hawthorne FC, Kato A, Kisch HJ, Krivovichev VG, Linthout K, Laird J, Mandarino JA, Maresch WV, Nickel EH, Rock NMS, Schumacher JC, Smith DC, Stephenson NCN, Ungaretti L, Whittaker EJW, Guo Y (1997) Nomenclature of amphiboles: Report of the subcommittee on amphiboles of the International Mineralogical Association, Commission on New Minerals and Mineral Names. Can Mineral 35:219–246
Ledlie I (1998) Geology of the Teresa de Colmo Fe-Cu Breccia deposit. MSc Thesis, University of Tasmania
Li R, Chen H, Xia X, Yang Q, Li L, Xu J, Huang C, Danyushevsky LV (2017) Ore fluid evolution in the giant Marcona Fe-(Cu) deposit, Perú: Evidence from in-situ sulfur isotope and trace element geochemistry of sulfides. Ore Geol Rev 86:624–638
Li R, Chen H, Xia X, Yang Q, Danyushevsky LV, Lai C (2018) Using integrated in-situ sulfide trace element geochemistry and sulfur isotopes to trace ore-forming fluids: Example from the Mina Justa IOCG deposit (southern Perú). Ore Geol Rev 101:165–179
Liu W, McPhail DC (2005) Thermodynamic properties of copper chloride complexes and copper transport in magmatic-hydrothermal solutions. Chem Geol 221:21–39
Liu Y, Fan Y, Zhou T, Whit NC, Hong H, Zhang W (2018) In-situ LA-ICP-MS trace element analyses of magnetite from Mesozoic iron oxide apatite (IOA) deposits in the Luzong volcanic basin, eastern China. J Asian Earth Sci 166:233–246
Lledó HL (2005) Experimental studies on the origin of iron deposits; and mineralization of Sierra La Bandera, Chile. PhD Thesis, SUNY- Univ. Binghamton.
Lledó HL, Jenkins DM (2008) Experimental Investigation of the upper thermal stability of Mg-rich actinolite: Implications for Kiruna-type iron deposits. Jour Petrol 49:225–238
Lledó HL, Naslund HR, Jenkins DM (2020) Experiments on phosphate-silicate liquid immiscibility with potential links to iron oxide apatite and nelsonite deposits. Contrib Miner Petrol 175:111
Loucks RR, Mavrogenes JA (1999) Gold solubility in supercritical hydrothermal brines measured in synthetic fluid inclusions. Science 284:2159–2163
Lucassen F, Franz G (1994) Arc related Jurassic igneous and meta-igneous rocks in the Coastal Cordillera of northern Chile/Region Antofagasta. Lithos 32:237–298
Lucassen F, Kramer W, Bartsch V, Wilke H-G, Franz G, Romer RL, Dulski P (2006) Nd, Pb, and Sr isotope composition of juvenile magmatism in the Mesozoic large magmatic province of northern Chile (18–27 degrees S): indications for a uniform subarc mantle. Contrib Miner Petrol 152:571–589
Maksaev V, Townley B, Palacios C, Camus F (2007) Metallic ore deposits. In: Moreno T, Gibbons W (eds) The geology of Chile. Geol Soc London, pp 179–199
Markey RJ, Hannah JL, Morgan JW, Stein HJ (2003) A double spike for osmium analysis of highly radiogenic samples. Chem Geol 200:395–406
Marschik R, Fontboté L (2001) The Punta del Cobre formation, Punta del Cobre-Candelaria area, Northern Chile. J South Amer Earth Sci 14:401–433
Marschik R, Fontboté L (2001) The Candelaria-Punta del Cobre Iron oxide Cu-Au[-Zn-Ag) deposits, Chile. Econ Geol 96:1799–1826
Marschik R, Fontignie D, Chiaradia M, Voldet P (2003) Geochemical and Sr–Nd–Pb–O isotope composition of granitoids of the Early Cretaceous Copiapó plutonic complex (27 83 00S), Chile. J South Amer Earth Sci 16:381–398
Martinsson O, Billström K, Broman C, Weihed P, Wanhainen C (2016) Metallogeny of the Northern Norrbotten Ore Province, northern Fennoscandian Shield with emphasis on IOCG and apatite-iron ore deposits. Ore Geol Rev 78:447–492
Mattinson JM (2005) Zircon U-Pb chemical abrasion (‘“CA-TIMS”’) method: combined annealing and multi-step partial dissolution analysis for improved precision and ac- curacy of zircon ages. Chem Geol 220:47–66
Mavor SP, Singleton JS, Heuser G, Gomila R, Seymour NM, Williams S, Arancibia G (2022) Sinistral shear during Middle Jurassic emplacement of the Matancilla Plutonic Complex in northern Chile (25.4° S) as evidence of oblique plate convergence during the early Andean orogeny. J South Am Earth Sci 120:104047. https://doi.org/10.1016/j.jsames.2022.104047
McBirney AR, Noyes RM (1979) Crystallization and layering of the Skaergaard intrusion. J Petrol 20:487–554
Ménard JJ (1995) Relationship between altered pyroxene diorite and the magnetite mineralization in the Chilean Iron Belt, with emphasis on the El Algarrobo iron deposits (Atacama region, Chile). Miner Deposita 30:268–274
Min K, Mundil R, Renne PR, Ludwig KR (2000) A test for systematic errors in Ar-Ar geochronology through comparison with U-Pb analysis of a 1.1-Ga rhyolite. Geochim Cosmochim Acta 64:73–98
Mpodozis C, Ramos V (1990) The Andes of Chile and Argentina. In: Ericksen GE, Cañas Pinochet MT, Reinemund JA (eds), Geology of the Andes and its relation to hydrocarbon and mineral resources. Circum-Pacific Council Energy and Mineral Resources, Earth Sci Series 11:59–90
Mungall J, Long K, Brenan J, Smythe D, Naslund H (2018) Immiscible shoshonitic and Fe-P-oxide melts preserved in unconsolidated tephra at El Laco volcano, Chile. Geology 46:255–258
Murakami T, Chakoumakos BC, Ewing RC, Lumpkin GR, Weber WJ (1991) Alpha-decay event damage in zircon. Amer Mineral 76:1510–1532
Naranjo JA, Puig A (1984) Hojas Taltal y Chañaral. Servicio Nac Geol Minería Carta Geol Chile Ser Geol Básica 62–63:140
Naranjo JA, Henríquez F, Nyström JO (2010) Subvolcanic contact metasomatism at El Laco Volcanic Complex, Central Andes. Andean Geol 37:110–120
Nasdala L, Reiners PW, Garver JI, Kennedy AK, Stern RA, Balan E, Wirth R (2004) Incomplete retention of radiation damage in zircon from Sri Lanka. Amer Mineral 89:219–231
Naslund HR, Henriquez F, Nystrom JO, Vivallo W, Dobbs FM (2002) Magmatic iron ores and associated mineralisation; examples from the Chilean High Andes and Coastal Cordillera. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PCG Publishing, Adelaide, pp 207–226
Naslund HR, Aguirre R, Dobbs FM, Henríquez F, Nyström JO (2000) The origin, emplacement, and eruption of ore magmas. In: Proceedings IX Congreso Geológico Chileno, Puerto Varas, pp 135–139
Naslund HR, Mungall JE, Henriquez F, Nyström JO, Lledó H, Lester GW, Aguirre R (2009) Melt inclusions in silicate lavas and iron oxide tephra of the El Laco vocano, Chile. In: Proceedings XII Congreso Geológico Chileno, Santiago, pp S8–033
Naumenko-Dèzes MO, Nägler TF, Villa IM, Mezger K (2018) 40K–40Ca – 87Rb-87Sr age comparison: constraints on the 40K decay constant. Geochim Cosmochim Acta 220:235–247
Nyström JO, Henríquez F (1994) Magmatic features of iron ores of the Kiruna type in Chile and Sweden: Ore texturas and magnetite geochemistry. Econ Geol 89:820–839
Nyström JO, Henríquez F, Naranjo JA, Naslund HR (2016) Magnetite spherules in pyroclastic iron ore at El Laco, Chile. Am Mineral 101:587–595
Ohmoto H (1986) Stable isotope geochemistry of ore deposits. Rev Mineral 16:491–555
Ohmoto H, Goldhaber MB (1997) Sulfur and carbon isotopes. In: Barnes, HL (ed) Geochemistry of hydrothermal ore deposits, 3rd edn Wiley, pp 517–611
Ohmoto H, Rye RO (1979) Isotopes of sulfur and carbon. In: Barnes, HL (ed) Geochemistry of hydrothermal ore deposits, 2nd edn Wiley, pp 509–561
Ootes L, Snyder D, Davis WJ, Acosta-Góngora P, Corriveau L, Mumin AH, Gleeson SA, Samson IM, Montreuil JF, Potter E, Jackson VA (2017) A Paleoproterozoic Andean-type iron oxide copper-gold environment, the Great Bear magmatic zone, Northwest Canada. Ore Geol Rev 81:123–139
Paces JB, Miller JD (1993) Precise U-Pb ages of Duluth Complex and related mafic intrusions, northeastern Minnesota: geochronological insights into physical, petrogenetic, paleomagnetic and tectonomagmatic processes associated with the 1.1 Ga midcontinent rift system. J Geophys Res 98:13997–14013
Palma G, Barra F, Reich M, Simon AC, Romero R (2020) A review of magnetite geochemistry of Chilean iron oxide-apatite (IOA) deposits and its implications for ore-forming processes. Ore Geol Rev 126:103748
Parada M, López-Escobar L, Oliveros V, Fuentes F, Morata D, Calderón M, Aguirre L, Féraud G, Espinoza F, Moreno H (2007) Andean magmatism. In: Moreno T, Gibbons W (eds) The Geology of Chile. Geol Soc, London, pp 115–146
Park CF (1961) A magnetite “flow” in northern Chile. Econ Geol 56:431–441
Paton C, Hellstrom J, Paul B, Woodhead J, Hergt J (2011) Iolite: freeware for the visualisation and processing of mass spectrometric data. J Anal at Spectrom 26:2508–2518
Philpotts AR (1981) Liquid immiscibility in silicate melt inclusions in plagioclase phenocrysts. Bull Mineral 104:317–324
Philpotts AR (1982) Compositions of immiscible liquids in volcanic rocks. Contrib Miner Petrol 80:201–218
Pichowiak S, Buchelt M, Damm KW (1990) Magmatic activity and tectonic setting of the early stage of the Andean cycle in North Chile. Mahlburg, Kay S, Rapella CW (eds), Plutonism from Antartica to Alaska, Geol Soc Amer Spec Paper 241:127–144
Pietruszka DK, Hanchar JM, Tornos F, Wirth R, Graham NA, Severin KP, Velasco F, Steele-MacInnis M, Bain WM (2022) Magmatic immiscibility and the origin of magnetite-(apatite) iron deposits. Research Square https://doi.org/10.21203/rs.3.rs-2156064/v1
Pokrovski GS, Akinfiev NN, Borisova AY, Zotov AV, Kouzmanov K (2014) Gold speciation and transport in geological fluids: insights from experiments and physical-chemical modelling. In: Garofalo PS, Ridley JR (eds) Gold-transporting hydrothermal fluids in the Earth’s crust. Geol Soc Spec Publ, London, pp 9–70
Pollard PJ (2000) Evidence of a magmatic fluid and metal source for Fe-oxide Cu-Au mineralisation. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PGC Publishing, Adelaide, pp 27–41
Pollard PJ (2006) An intrusion-related origin for Cu-Au mineralization in iron oxide-copper gold (IOCG) provinces. Miner Deposita 41:179–187
Poulson SR, Kubilius WP, Ohmoto H (1991) Geochemical behaviour of sulfur in granitoids during intrusion of the South Mountain Batholith, Nova Scotia, Canada. Geochim Cosmochim Acta 55:3809–3830
Ramos V, Aleman A (2000) Tectonic evolution of the Andes. In: Cordani UG, Milani EJ, Thomaz Filho A, Campos DA (eds) Tectonic evolution of South America, Sociedad Brasileira de Geologia, pp 635–685
Ray GE, Dick LA (2002) The Productora prospect in north-central Chile: An example of an intrusion-related Candelaria type Fe–Cu–Au hydrothermal system. In: Porter TM (ed) Hydrothermal iron oxide copper-gold & related deposits: A global perspective. PGC Publishing, Adelaide, pp 131–151
Renne PR, Karner DB, Ludwig KR (1998) Absolute ages aren’t exactly. Science 282:1840–1841
Rhodes AL, Oreskes N (1999) Oxygen Isotope composition of magnetite deposits at El Laco, Chile: Evidence of formation from isotopically heavy fluids. In: Skinner, BJ (ed), Geology and ore deposits of the Central Andes, Soc Econ Geol Spec Publ 7:333–351
Rhodes AL, Oreskes N, Sheets S (1999) Geology and rare element (REE) geochemistry of magnetite deposits at El Laco, Chile. In: Skinner, BJ (ed), Geology and ore deposits of the Central Andes, Soc Econ Geol Spec Publ 7:299–332
Richards J, Mumin H (2013) Magmatic-hydrothermal processes within an evolving Earth: Iron oxide-copper-gold and porphiry Cu ± Mo ± Au deposits. Geology 41:767–770
Rieger AA, Marschik R, Díaz M, Hölzl S, Chiaradia M, Akker B, Spangenberg JE (2010) The hypogene iron oxide copper-gold mineralization in the Mantoverde district, Northern Chile. Econ Geol 105:1271–1299
Rogers G, Hawkesworth CJ (1989) A geochemical traverse across the North Chilean Andes: evidence for crust generation from the mantle wedge. Earth Planet Sci Let 91:271–285
Rojas P, Barra F, Deditius A, Reich M, Simon A, Roberts M, Rojo M (2018) New contributions to the understanding of Kiruna-type iron oxide-apatite deposits revealed by magnetite ore and gangue mineral geochemistry at the El Romeral deposit, Chile. Ore Geol Rev 93:413–435
Rotherham JF (1997) A measomatic origin for the iron-oxide Au-Cu Starra orebodies, Eastern Fold Belt, Mount Isa Inlier. Miner Deposita 32:205–218
Rudnick RL, Gao S (2003) Composition of the Continental Crust. In: Holland HD, Turekian KK (eds), Treatise on Geochemistry, Oxford, vol. 3 pp 1–64
Rusk BG, Reed MH, Dilles JH, Klemm LM, Heinrich CA (2004) Compositions of magmatic hydrothermal fluids determined by LA-ICP-MS of fluid inclusions from the porphyry copper–molybdenum deposit at Butte, MT. Chem Geol 210:173–199
Rusk BG, Oliver NHS, Zhang D, Brown A, Lilly R, Jungmann D (2009) Compositions of magnetite and sulfides from barren and mineralized IOCG deposits in the Eastern succession of the MT Isa inlier, Australia. In: Proceedings GSA Annual Meeting, Portland, abstract 26–10
Salazar E, Barra F, Reich F, Simon A, Leisen M, Palms G, Romero R, Rojo M (2019) Trace element geochemistry of magnetite from the Cerro Negro Norte iron oxide−apatite deposit, northern Chile. Miner Deposita 55:409–428
Scheuber E, Reutter KJ (1992) Magmatic arc tectonics in the Central Andes between 21º and 25ºS. Tectonophysics 205:127–140
Schoene B, Bowring SA (2006) U-Pb systematics of the McClure Mountain syenite: thermochronological constraints on the age of the 40Ar/39Ar standard MMhb. Contrib Miner Petrol 151:615–630
Schutesky ME, de Oliveira CG (2020) From the roots to the roof: An integrated model for the Neoarchean Carajás IOCG System, Brazil. Ore Geol Rev 127:103833
Shenberger DM, Barnes HL (1989) Solubility of gold in aqueous sulfide solutions from 150 to 350ºC. Geochim Cosmochim Acta 53:269–278
Sidder GB (1984) Ore genesis at the Monterrosas deposit in the Coastal batholith of Ica, Peru: Ph.D. thesis, Oregon State Univ
Sillitoe RH (2003) Iron oxide-copper-gold deposits: an Andean view. Miner Deposita 38:787–812
Sillitoe RH, Burrows D (2002) New field evidence bearing on the origin of the El Laco magnetite deposit, Northern Chile. Econ Geol 97:1101–1109
Simon JI, Renne PR, Mundil R (2008) Implications of pre-eruptive magmatic histories of zircons for U-Pb geochronology of silicic extrusions. Earth Planet Sci Let 266:182–194
Simon AC, Knipping J, Reich M, Barra F, Deditius AP, Bilenker L, Childress T (2018) Kiruna-Type Iron Oxide-Apatite (IOA) and Iron Oxide Copper-Gold (IOCG) deposits form by a combination of igneous and magmatic-hydrothermal processes: evidence from the Chilean Iron Belt. Soc Econ Geol Spec Publ, pp 89–114
Skirrow RG (2022) Iron oxide copper-gold (IOCG) deposits–A review (part 1): Settings, mineralogy, ore geochemistry and classification. Ore Geol Rev 140:104569
Sláma J, Košler J, Condon DJ, Crowley JL, Gerdes A, Hanchar JM, Horstwood MSA, Morris GA, Nasdala L, Norberg N, Schaltegger U, Tubrett M, Whitehouse MJ (2008) Plešovice zircon – a new natural standard for U-Pb and Hf isotopic microanalysis. Chem Geo 249:1–35. https://doi.org/10.1016/j.chemgeo.2007.11.005
Smoliar MI, Walker RJ, Morgan JW (1996) Re-Os ages of group IIA, IIIA, IVA, and IVB iron meteorites. Science 271:1099–1102
Söderlund U, Patchett PJ, Vervoort JD, Isachsen CE (2004) The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions. Earth Planet Sci Lett 219:311–324
Steiger RH, Jäger E (1977) Subcommission on geochronology: convention on the use of decay constants in geo- and cosmochronology. Earth Planet Sci Let 36:359–362
Stein HJ, Markey RJ, Morgan JW, Hannah JL, Scherstén A (2001) The remarkable Re-Os chronometer in molybdenite: how and why it works. Terra Nova 13:479–486
Stein H, Scherstén A, Hannah J, Markey R (2003) Sub-grain scale decoupling of Re and 187Os and assessment of laser ablation ICP-MS spot dating in molybdenite. Geochim Cosmochim Acta 67:3673–3686
Tornos F, Velasco F, Barra F, Morata D (2010) The Tropezón Cu-Mo-(Au) deposit deposit, Northern Chile: the missing link between IOCG and porphyry coppper systems? Miner Deposita 45:313–321
Tornos F, Velasco F, Hanchar JM (2016) Iron-rich melts, magmatic magnetite, and superheated hydrothermal systems: The El Laco deposit, Chile. Geology 44:427–430
Tornos F, Velasco F, Hanchar JM (2017) The magmatic to magmatic-hydrothermal evolution of the El Laco deposit (Chile) and its implications for the genesis of magnetite-apatite deposits. Econ Geol 112:1595–1628
Tornos F, Hanchar JM, Munizaga R, Velasco F, Galindo C (2021) The role of the subducting slab and melt crystallization in the formation of magnetite-(apatite) systems, Coastal Cordillera of Chile. Miner Deposita 56:253–278
Tornos F (2011) Magnetite-Apatite and IOCG deposits formed by magmatic-hydrothermal evolution of complex calc-alkaline melts. In: Proceedings 11th SGA Biennial Meeting, Antofagasta, pp 26–28
Tornos F, Hanchar JM, Colquehuanca F, Chavez E (2023) Origin of the magnetite-(apatite) and IOCG deposits of the coastal Cordillera of the Andes. In. Proceedings ProExplo 2023, Lima
Travisany V, Henríquez F, Nyström JO (1995) Magnetite lava flows in the Pleito-Melón district of the Chilean Iron Belt. Econ Geol 90:438–444
Treloar PJ, Colley H (1996) Variations in F and C1 contents in apatites from magnetite-apatite ores in northern Chile, and their ore-genetic implications. Min Mag 60:285–301
Tristá D, Kojima S (2003) Mineral paragenesis and fluid inclusions of some pluton-hosted vein-type copper deposits in the Coastal Cordillera, Northern Chile. US Geol Surv Resour Geol 53:21–28
Troll VR, Weis FA, Jonsson E, Andersson UB, Majidi SA, Högdahl K, Harris C, Millet MA, Chinnasamy SS, Kooijman E, Nilsson KP (2019) Global Fe-O isotope correlation reveals magmatic origin of Kiruna-type apatite-iron-oxide ores. Nat Comm 10:1712
Valley PM, Fisher CM, Hanchar JM, Lam R, Tubrett M (2010) Hafnium isotopes in zircon: A tracer of fluis-rock interation during magnetite-apatite (“Kiruna-type”) mineralization. Chem Geol 275:208–220
Veksler IV, Dorfman AM, Borisov AA, Wirth R, Dingwell DB (2007) Liquid Immiscibility and the Evolution of Basaltic Magma. J Petrol 48:2187–2210
Velasco F, Tornos F, Hanchar JM (2016) Immiscible iron-and silica-rich melts and magnetite geochemistry at the El Laco volcano (northern Chile): Evidence for a magmatic origin for the magnetite deposits. Ore Geol Rev 79:346–366
Verdugo-Ihl MR, Ciobanu CL, Courtney-Davies L, Cook NJ, Slattery A, Ehrig K, Tornos F, Hanchar JM (2022) U-Pb geochronology and mineralogy of hematite from Mantoverde and Carmen de Cobre, northern Chile: constraints on Andean IOCG minerlization. Econ Geol 117(4):943–960
Vervoort JD, Plank T, Prytulak J (2011) The Hf–Nd isotopic composition of marine sediments. Geochim Cosmochim Acta 75:5903–5926
Vervoort JD (2010) Hf analyses in zircon by LA-MC-ICPMS: promise and pitfalls. Geol Soc Amer Abstr Programs 42:286–289
Vervoort JD, Blichert-Toft J (1999) Evolution of the depleted mantle: Hf isotope evidence from juvenile rocks through time. Geochim Cosmochim Acta 6:533–556. https://doi.org/10.1016/S0016-7037(98)00274-9
Vila T, Lindsay N, Zamora R (1996) Geology of the Manto Verde copper deposit, northern Chile: A specularite-rich, hydrothermal-tectonic breccia related to the Atacama Fault Zone. In: Camus F, Sillitoe RH, Petersen R (eds), Andean Copper deposits: New discoveries, mineralization styles and metallogeny. Soc Econ Geol Spec Publ 5:157–169
Villa IM, Ruggieri G, Puxeddu M, Bertini G (2006) Geochronology and isotope transport systematics in a subsurface granite from the Larderello-Travale geothermal system (Italy). J Volcanol Geotherm Res 152:20–50
Vivallo W, Henríquez F (1998) Génesis común de los yacimientos estratoligados y vetiformes de cobre del Jurásico Medio a Superior en la Cordillera de la Costa, Región de Antofagasta, Chile. Rev Geol Chile 25:199–228
Vivallo W, Henríquez F, Espinoza S (1995) Los depósitos de hierro del tipo magnetita-apatita: geoquímica de las rocas volcánicas asociadas y potencialidad de la mena de hierro como fuente de mineralización de oro. Rev Geol Chile 22:159–175
Vivallo W, Díaz A, Jorquera R (2008) Yacimientos metalíferos de la Región de Atacama. Servicio Nac Geol Minería Carta Geol Chile Serie Recur Minerales Energéticos 27:72 (1 map)
Vivallo W, Henríquez F (1997) Relación genética entre los yacimientos estratoligados de Cu ("Tipo Manto"), de Cu-Fe-Au y de hierro del tipo Kiruna. In: Proceedings VIII Congreso Geológico Chileno, Antofagasta, pp 1189–1193
Vivallo W, Díaz A, Gelcich S, Lledó H (2000) Estilos y tipos de mineralización del Jurásico y Cretácico Inferior en la Cordillera de la Costa de la región de Copiapo, Chile. In: Proceedings IX Congreso Geológico Chileno, Santiago, pp 179–182
Wager LR (1961) A note on the origin of ophitic texture in the chilled olivine gabbro of the Skaergaard intrusion. Geol Mag 98:353–366
Wager LR, Brown GM, Wadsworth WJ (1960) Types of igneous cumulates. J Petrol 1:73–85
Weis F, Troll VR, Jonsson E, Högdahl K, Harris C, Sun W, Nilsson KP, Dahrén B (2021) Absence of hydrothermal oxygen variations in host rocks supports magmatic origin of the giant Grängesberg iron oxide-apatite (IOA) deposit, central Sweden. Intern J Earth Sci. https://doi.org/10.1007/s00531-021-02122-9
Wiedenbeck M, Hanchar JM, Peck WH, Sylvester P, Valley J, Whitehouse M, Kronz A, Morishita Y, Nasdala L, Fiebig J, Franchi I, Girard JP, Greenwood RC, Hinton R, Kita N, Mason PRD, Norman M, Ogasawara M, Piccoli PM, Rhede D, Satoh H, Schulz-Dobrick B, Skar O, Spicuzza MJ, Terada K, Tindle A, Togashi S, Vennemann T, Xie Q, Zheng YF (2004) Further characterisation of the 91500 zircon crystal. Geostand Geoanalytical Res 28:9–39
Williams P, Barton MD, Johnson DA, Fontboté L, Haller AD, Mark G, Oliver NHS, Marschik R (2005) Iron oxide copper-gold deposits: Geology, space-time distribution, and possible modes of origin. In: Hedenquist JW, Thompson JFH, Goldfarb RJ, Richards JP (eds) Econ Geol - One hundredth anniversary Volume Soc Econ Geol, Littleton, pp 371–406
Xavier PR, Monteiro LVS, Moreto CPN, Silva Pestilho AL, Coelho de Melo GH, Da Silva MAD, Aires B, Ribeiro C, Freitas e Silva FH (2012) The iron oxide-copper-gold systems of the Carajás mineral province, Brazil. Econ Geol Spec Publ 16:433–454
Zezin DY, Migdisov AA, Williams-Jones AE (2007) The solubility of gold in hydrogen sulfide gas: An experimental study. Geochim Cosmochim Acta 71:3070–3081
Zhang D, Oliver NHS, Rusk BG, Chang Z, Colman R (2017) Trace element geochemistry of magnetite from Ernest Henry IOCG deposit, Australia. 14th SGA Bienn Meet 3:975–978
Acknowledgements
We acknowledge Justo Esquivel (Montecristo district), and Guillermo López [SERNAGEOMIN (The National Geology and Mining Service of Chile)] for their assistance during field work. Numerous people are thanked for assistance in the labs and discussions, including Markus Wälle, Rebecca Lam, Wanda Aylward, Waldo Vivallo, and Alejandro Díaz. We acknowledge Tobias Bauer, Wyatt Bain and Bernd Lehmann for their thoughtful manuscript review.
Funding
Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature. This project was financially supported by an NSERC discovery grant to J.M. Hanchar (RGPIN/004649-2015), and a RTI2018-099157-A-I00 (MCI/AEI/FEDER, UE) research grant to F. Tornos.
Author information
Authors and Affiliations
Contributions
LM, FT, and JMH conceived and developed the project. LM collected the samples. LM, IV, HS, and AD collected and curated the data. LM synthesized and visualized the data and wrote the original manuscript. All authors contributed to the review and approval of the final manuscript prior to submission.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Editorial handling: M. Steele-Macinnis
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mateo, L., Tornos, F., Hanchar, J.M. et al. The Montecristo mining district, northern Chile: the relationship between vein-like magnetite-(apatite) and iron oxide-copper–gold deposits. Miner Deposita 58, 1023–1049 (2023). https://doi.org/10.1007/s00126-023-01172-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00126-023-01172-0