Implications of Hf Isotopes for the Evolution of the Mantle Source of Magmas Associated with the Giant El Teniente Cu-Mo Megabreccia Deposit, Central Chile

: We have determined the Hf isotopic compositions of 12 samples associated with the giant El Teniente Cu-Mo megabreccia deposit, central Chile. The samples range in age from ≥ 8.9 to 2.3 Ma and provide information about the temporal evolution of their magmatic sources from the Late Miocene to Pliocene. Together with previously published data, the new analysis indicates a temporal decrease of 10 ε Hf(t) units, from + 11.6 down to + 1.6, in the 12.7 m.y. from 15 to 2.3 Ma. These variations imply increasing incorporation of continental crust through time in the magmas that formed these rocks. The fact that the samples include mantle-derived olivine basalts and olivine lamprophyres suggests that these continental components were incorporated into their mantle source, and not by intra-crustal contamination (MASH). We attribute the increase, between the Middle Miocene and Pliocene, of crustal components in the subarc mantle source below El Teniente to be due to increased rates of subduction erosion and transport of crust into the mantle. The deposit formed above a large, long-lived, vertically zoned magma chamber that developed due to compressive deformation and persisted between ~7 to 4.6 Ma. Progressively more hydrous mantle-derived maﬁc magmas feed this chamber from below, providing heat, H 2 O, S and metals, but no unique “fertile” Cu-rich magma was involved in the formation of the deposit. As the volume of these mantle-derived magmas decreased from the Late Miocene into the Pliocene, the chamber crystallized and solidiﬁed, producing felsic plutons and large metal-rich magmatic-hydrothermal breccias that emplaced Cu and S into the older ( ≥ 8.9 Ma) maﬁc host rocks of this megabreccia deposit.


Introduction
The isotopic compositions ( 87 Sr/ 86 Sr ≥ 0.7050; ε Nd ≤ −2; ε Hf(t) ≤ +2) of some recently erupted Andean volcanic rocks indicate that they have incorporated continental crust [1][2][3][4][5]. Two different processes for the incorporation of crust into Andean magmas have been proposed ( Figure 1). One, MASH (Mixing, Assimilation, Storage and Homogenization), involves intra-crustal assimilation as mantle-derived magmas rise from their source through the crust to the surface [2]. Another, PEELED, involves a mantle source region contaminated by crustal components subducted below the mantle Figure 1. Two alternative explanations, as originally illustrated by Hickey-Vargas in 1991 [1], for the incorporation of continental crust into Andean magmas. MASHED on the left involves Mixing, Assimilation, Storage and Homogenization in a MASH zone located at either the lower crust/mantle boundary or within the crust [2]. PEELED, on the right, involves recycling of subduction eroded crust into the mantle wedge source [3][4][5]. J-K and P indicate the position of the Jurassic-Cretaceous and Pliocene plutonic belts, respectively, in the Andes of northern Chile. The recently active volcanic arc is ~250 km east of the Jurassic arc as a result in part of subduction erosion [6]. (Reprinted from Nature, volume 350, Hickey-Vargas, R. Peeled or MASHed?, page 350, Copyright 1991, with permission from Springer Nature.).
The amount of crustal components incorporated into recently erupted Andean volcanic rocks varies along strike from south-to-north, as do various geologic and tectonic features such as crustal thickness and age, subduction angle, trench depth and sediment fill, and the localized presence of buoyant features such as oceanic seamount chains and spreading ridges being subducted below the South American continental margin (Figure 2; [4]). For instance, in the Central Volcanic Zone (CVZ) of northern Chile, the crust is much thicker (≥70 km) than below the central part of the Southern Volcanic Zone (≤35 km; CSVZ; Figure 2) and the isotopic compositions of the recently erupted volcanic rocks suggests more crustal components have been incorporated in CVZ than CSVZ magmas [3][4][5]. This could reflect greater amounts of crustal assimilation by mantle-derived magmas as they rose through the thicker CVZ crust, but in this same region, the trench is devoid of sediment due to the arid conditions in northern Chile, and compared to the CSVZ, much more significant amounts of subduction erosion, at rates estimated to be 50-70 km 3 /km/m.y. since ~150 Ma (Figure 2; [5]), have truncated the continental margin [6] and caused the arc to migrate eastward ~250 km since the Jurassic (Figure 1). The subduction of tectonically eroded crust has been invoked not only to explain the crustal components in CVZ magmas [3][4][5]7,8], but also the uplift and crustal thickening observed in this part of the Andes [9].
In the same way, at the northern end of the Andean SVZ (NSVZ; Figures 2 and 3), where the MASH model was first presented [2], the crust thickens to ≥60 km compared to ≤35 km below the CSVZ, and the isotopic compositions of recently erupted volcanic rocks indicate that they contain greater amount of crustal components [3][4][5]10,11]. However, the subduction of the Juan Fernández Ridge has caused the subduction angle to flatten, increasing the rate of subduction erosion west of the NSVZ to between 115 and230 km 3 /km/m.y. over the last 10 m.y., during which time the active volcanic arc between 33-34° S migrated ~50 km to the east ( Figure 3). Stern [3][4][5] therefore proposed that crustal components were preferentially incorporated in NSVZ compared to CSVZ magmas by increased rates of subduction erosion and mantle source region contamination by subducted components, and not by intra-crustal MASH processes.
Temporal variations in the isotopic compositions of Andean igneous rocks also are not unambiguously related to either changes in crustal thickness or rates of subduction erosion. Below the NSVZ, in the region where the El Teniente deposit formed in the Late Miocene (Figures 2 and 3), Figure 1. Two alternative explanations, as originally illustrated by Hickey-Vargas in 1991 [1], for the incorporation of continental crust into Andean magmas. MASHED on the left involves Mixing, Assimilation, Storage and Homogenization in a MASH zone located at either the lower crust/mantle boundary or within the crust [2]. PEELED, on the right, involves recycling of subduction eroded crust into the mantle wedge source [3][4][5]. J-K and P indicate the position of the Jurassic-Cretaceous and Pliocene plutonic belts, respectively, in the Andes of northern Chile. The recently active volcanic arc is 250 km east of the Jurassic arc as a result in part of subduction erosion [6]. (Reprinted from Nature, volume 350, Hickey-Vargas, R. Peeled or MASHed?, page 350, Copyright 1991, with permission from Springer Nature.).  Figure 3. The figure also shows tectonic features such as the Chile Trench, which is the boundary between the Nazca and South American plates, and the depth in kilometers (100 and 150 km dashed lines) to the Benioff Zone below South America. The locus of subduction of the Juan Fernández Ridge occurs at the boundary between the Andean Flat-Slab segment, below which the subduction angle is very low as indicated by the depths to the upper boundary of the subducted slab, and the northern and central parts of the Andean Southern Volcanic Zone (NSVZ and CSVZ) of active volcanoes, below which subduction angle is steeper. The numbers in the green boxes are the estimated rates of subduction erosion in km 3 /km/m.y. in each region [5]. Subduction erosion rates are higher in northern compared to south-central Chile because of the arid conditions in the north, resulting in a trench devoid of sediment, and humid conditions in the south, resulting in a trench filled with sediment [4,5,9]. Rates are highest just south of the locus of subduction of the Juan Fernández Ridge due to the effects of the subduction of this buoyant feature [5,26,27]. (Reprinted from Geology, volume 22, Skewes, M.A. and Stern, C.R. Tectonic trigger for the formation of late Miocene Cu-rich breccia pipes in the Andes of central Chile, page 552, Copyright 1994, with fair use permission from the Geologic Society of America.).
Temporal variations in the isotopic compositions of Andean igneous rocks also are not unambiguously related to either changes in crustal thickness or rates of subduction erosion. Below the NSVZ, in the region where the El Teniente deposit formed in the Late Miocene (Figures 2 and 3), the crust has thickened and been uplifted since the Middle Miocene [12][13][14][15][16][17][18][19][20], and Sr, Nd and Pb isotopic data indicate increased proportions of continental crust in the igneous rocks emplaced through time (Figures 4 and 5; [10,[21][22][23][24]). However, as described above, the southward migration of the locus of subduction of the Juan Fernández Ridge has also caused, during this same time period, flattening of the angle of subduction, >50 km of eastward migration of the volcanic arc, and increased rates of subduction erosion (Figure 3; [3][4][5]10,[25][26][27]).
In summary, the spatial and temporal variations in geologic factors such as crustal thickness and rates of subduction erosion rates do not unambiguously favor MASH versus source region contamination in producing the spatial and temporal isotopic variations that indicate increased incorporation of continental crust in some Andean magmas. An alternative approach to distinguish these two processes is geochemical modelling focused specifically on mafic olivine-bearing basaltic rocks clearly derived from the mantle, without significant evidence for or possibility of intra-crustal contamination. For instance, to produce the northward increase of 87 Sr/ 86 Sr from 0.7039 to 0.7050 and decrease in ε Nd from >+5 to <−1 observed in the olivine basalts erupted from CSVZ compared to NSVZ volcanoes (Figures 3-5), would require intra-crustal assimilation (MASH) of >40 wt. % of Andean crust ( Figure 4; [3,10,21,22]) because of the relative high Sr ≥ 50 ppm and Nd ≥ 9 ppm of Andean basalts. For NSVZ olivine basalts with~50 wt. % SiO 2 , but with no macro or microscopic evidence of crustal xenoliths or xenocrysts, this is considered to be an unacceptably large amount of crustal assimilation, as it would increase SiO 2 to~60 wt. %. In contrast, because of the relatively low Sr = 36 ppm and Nd = 1.8 ppm of Andean subarc mantle, these variations could be produced by incorporation in the mantle source region of only 5 wt. % more subducted crust and sediment below the NSVZ (Figure 4), where current subduction erosion rates are estimated as between 115 and 230 km 3 /km/m.y. (Figure 2), than the 1 wt. % estimated to have been incorporated in the mantle below the CSVZ, where subduction erosion rates are estimated as only~35 km 3 /km/m.y. (Figure 2). Similar temporal isotopic changes to those observed between recently erupted CSVZ and NSVZ basalts occurred between the Miocene and Pliocene in igneous rocks associated with the El Teniente Cu-Mo deposit (Figures 3-5; [21][22][23]). As a contribution to the understanding of Andean magmagenesis and its relation to the formation of the giant El Teniente ore deposit, we have determined the Hf isotopes of igneous rocks ranging in age from ≥8.9 to 2.3 Ma from the area of the deposit. These samples (Table 1) include an olivine basalts (8.9 Ma), an olivine lamprophyre (3.1 Ma), and an olivine-bearing basaltic andesites (2.3 Ma), as well as other mafic and felsic plutons associated with the deposit. The data obtained provide information about the temporal evolution of their magmatic sources from the Late Miocene to Pliocene, and constrain the relative role of MASH and mantle source region contamination processes in the generation of Andean magmas. assimilation, as it would increase SiO2 to ~60 wt. %. In contrast, because of the relatively low Sr = 36 ppm and Nd = 1.8 ppm of Andean subarc mantle, these variations could be produced by incorporation in the mantle source region of only 5 wt. % more subducted crust and sediment below the NSVZ (Figure 4), where current subduction erosion rates are estimated as between 115 and230 km 3 /km/m.y. (Figure 2), than the 1 wt. % estimated to have been incorporated in the mantle below the CSVZ, where subduction erosion rates are estimated as only ~35 km 3 /km/m.y. (Figure 2).  [10]. Arrows show relative amounts of frontal-arc migration, forearc loss, and backarc shortening. Northwest-southeast trending dashed lines show offsets in the modern volcanic front that separate the Southern Volcanic Zone (SVZ) into the northern, transitional, and central segments [4]. In the active arc region, lines connect outcrop patterns marking early Miocene (pink), middle to  [10]. Arrows show relative amounts of frontal-arc migration, forearc loss, and backarc shortening. Northwest-southeast trending dashed lines show offsets in the modern volcanic front that separate the Southern Volcanic Zone (SVZ) into the northern, transitional, and central segments [4]. In the active arc region, lines connect outcrop patterns marking early Miocene (pink), middle to late Miocene (red), and the Pleistocene to Holocene volcanic centers of the SVZ (triangles) magmatic fronts. Arrows between the lines indicate inferred distance (given in circles) of frontal-arc migration from 19 to 16 Ma and from 7 to 3 Ma. In the forearc, lines between the trench and the coast show inferred early Miocene (short dashed) and middle to late Miocene (long dashed) coastlines under the assumption that the distance of frontal-arc migration equals the width of missing coast. Arrows between the lines in this region indicate distance (shown in circles) of inferred loss by subduction erosion from ca. 19 to 16 Ma. and 7 to 3 Ma. In the backarc, the length and position of arrows show the location and proportional amounts of crustal shortening over the past 20 m.y. inferred from structural profiles. Also shown are other outcrop patterns that have long been used as evidence for forearc subduction erosion along this margin [6]. The first is the northward narrowing and disappearance of the Paleozoic high pressure (P1) and low pressure (P2) paired metamorphic and granitoid (Pzg) belts along the coast. The second is the presence of Jurassic arc rocks (marked by J) along the coast north of 33 • S, but inland near the SVZ at~38 • S. K indicates Cretaceous magmatic rocks. ε Nd is +5 for active SVZ volcanoes south of 36 • S, and ≤−1 for those in the northern SVZ north of 34 • S, as a result of increased mantle source region contamination by subducted crust due to the northward increase in the rate of subduction erosion associated with the subduction of the Juan Fernández Ridge at 33 • S [3][4][5]. El Teniente occurs at the blue diamond. of 34° S, as a result of increased mantle source region contamination by subducted crust due to the northward increase in the rate of subduction erosion associated with the subduction of the Juan Fernández Ridge at 33° S [3][4][5]. El Teniente occurs at the blue diamond.   determined the Hf isotopes of igneous rocks ranging in age from ≥8.9 to 2.3 Ma from the area of the deposit. These samples (Table 1) include an olivine basalts (8.9 Ma), an olivine lamprophyre (3.1 Ma), and an olivine-bearing basaltic andesites (2.3 Ma), as well as other mafic and felsic plutons associated with the deposit. The data obtained provide information about the temporal evolution of their magmatic sources from the Late Miocene to Pliocene, and constrain the relative role of MASH and mantle source region contamination processes in the generation of Andean magmas.  [23]). Three curves illustrate the progressive changes in SiO2 and εNd for  crystal-liquid fractionation without assimilation,  intra-crustal assimilation combined with crystalliquid fractionation (MASH), and  mantle source region contamination. Crystal-liquid fractionation increases SiO2 without changing isotopic composition. MASH increases SiO2 as εNd decreases. Source region contamination decreases εNd in the mantle, but mantle partial melting (curves ) for the modified mantle still produces mafic olivine-bearing basalts and lamprophyres. Tic marks along each curve indicate the same amount of crustal contamination as in the MASH and source region contamination models illustrated in Figure 4. Although the youngest felsic plutons in the deposit are the same age (4.6 to 7.1 Ma) as rocks from the regionally defined Young Plutonic Complex of small barren plutons outside the deposit (5 to 6.6 Ma; [10]), they are isotopically similar to the older hostrocks of the deposit as well as to the Teniente Volcanic and Plutonic Complex rocks. This is due to the fact that they formed above a large magma chamber within which mantle-derived magmas were mixed, stored and homogenized (but without crustal assimilation), so that the temporal isotopic changes in these mantle-derived magmas through time are not observed in the progressively younger felsic plutons within the deposit as a result of the large volume of the magma chamber into which these mantle-derived magmas mixed. In contrast, the small barren plutons of the Young Plutonic Complex outside the deposit do reflect these temporal changes in isotopic composition because they do not form above such a large magma chamber.   [23]). Three curves illustrate the progressive changes in SiO 2 and ε Nd for 1 crystal-liquid fractionation without assimilation, 2 intra-crustal assimilation combined with crystal-liquid fractionation (MASH), and 3 mantle source region contamination. Crystal-liquid fractionation increases SiO 2 without changing isotopic composition. MASH increases SiO 2 as ε Nd decreases. Source region contamination decreases ε Nd in the mantle, but mantle partial melting (curves 4 ) for the modified mantle still produces mafic olivine-bearing basalts and lamprophyres. Tic marks along each curve indicate the same amount of crustal contamination as in the MASH and source region contamination models illustrated in Figure 4. Although the youngest felsic plutons in the deposit are the same age (4.6 to 7.1 Ma) as rocks from the regionally defined Young Plutonic Complex of small barren plutons outside the deposit (5 to 6.6 Ma; [10]), they are isotopically similar to the older host-rocks of the deposit as well as to the Teniente Volcanic and Plutonic Complex rocks. This is due to the fact that they formed above a large magma chamber within which mantle-derived magmas were mixed, stored and homogenized (but without crustal assimilation), so that the temporal isotopic changes in these mantle-derived magmas through time are not observed in the progressively younger felsic plutons within the deposit as a result of the large volume of the magma chamber into which these mantle-derived magmas mixed. In contrast, the small barren plutons of the Young Plutonic Complex outside the deposit do reflect these temporal changes in isotopic composition because they do not form above such a large magma chamber.   [21][22][23]. ε Hf(t) from [28].

El Teniente Cu-Mo Deposit
The giant El Teniente Cu-Mo deposit, located in central Chile (34 • 05 S, 70 • 21 W; Figures 2 and 3) is one of the largest such deposits in the world [28][29][30][31][32][33][34][35][36][37][38]. It is the southernmost in a belt of three giant Late Miocene Cu-Mo deposits which also includes Los Pelambres and Río Blanco-Los Bronces [39,40]. El Teniente has been described as a porphyry deposit formed around the central Teniente Dacite Porphyry dike ( Figure 6; [31]), or alternatively as a magmatic-hydrothermal breccia deposit formed either in association with the intrusion of dacite to diorite dikes [37] or by the generation of multiple large mineralized breccias, including the central Braden breccia pipe, by exsolution of high-temperature mineral-rich magmatic-hydrothermal fluids from the roof of an underlying large, long-lived and vertically stratified magma chamber recharged from below by mantle-derived magmas [23,33,34,38]. The two other Miocene deposits in this belt are also characterized by the presence of multiple large mineralized magmatic-hydrothermal breccias [38][39][40][41][42].
Pb, S, O and H isotopic studies of the igneous rocks and ore minerals in the deposit indicate that El Teniente is clearly an orthomagmatic deposit, with Cu, Mo and S being derived from the associated igneous rocks [33,34,39,[43][44][45]. Igneous rocks spatially related with the deposit range in age from mid Tertiary to Pliocene. Over this time period, the volume of magmatic rocks decreased progressively [22,23,38,39], and in the Pliocene the locus of Andean arc magmatism at the latitude of El Teniente (34 • S), near the current northern end of the Andean Southern Volcanic Zone (NSVZ; Figures 2 and 3), migrated~50 km to the east as the subduction angle flattened [25]. This, we suggest, was due to the southward migration of the locus of subduction of the Juan Fernández Ridge [21,[25][26][27]38,39].
The igneous rock in the vicinity of El Teniente include continental volcanic rocks, up to 3300 m thick, of the Oligocene to Early Miocene (≥15 Ma) Coya-Machali or Abanico Formation ( Figure 3; [15,46,47]), erupted above thin (<30 km) crust, or in a transtensional intra-arc basin [48]. These formed by relatively high degrees of partial melting of subarc mantle modified to a small degree by the influx from below of slab-derived components ( Figure 5; [10,17,46]). This is indicated by their low La/Yb ratios (2 to 6), as well as their low initial 87 Sr/ 86 Sr (0.7033-0.7039) and high ε Nd (+6.5 to +4; Figures 4 and 5). Although the Coya-Machali Formation volcanic rocks do not crop out either within or in the immediate vicinity of the El Teniente deposit (Figures 3 and 6), these rocks occur both to the west and the east of the deposit, and they almost certainly also occur at depth below the deposit.
Minerals 2019, 9, x FOR PEER REVIEW 8 of 23 Figure 6. Geological map of level Teniente 5 (2284 m above sea level) in the mine [33,34]. The Dacitic Porphyries north of the Sewell Tonalite, are mapped as a distal portion of this pluton, although they are younger [35] and have an independent origin [37]. The spatial extent of biotite breccias is projected onto this level from where they have been mapped between levels Teniente 4 and 8. Extrusive rocks of the Miocene Farellones Formation, locally referred to as the Teniente Volcanic Complex [10], are the oldest rocks exposed in the immediate area surrounding the deposit (Figure 3). The Farellones Formation, a sequence of >2500 m of lavas, volcanoclastic rocks, dikes, sills and stocks of basaltic to rhyolitic composition, was erupted after an Early Miocene (19-16 Ma episode of crustal deformation, thickening and uplift [10,14]. The Teniente Volcanic Complex near the deposit has been correlated with the upper part of this formation and dated between 14.4 and 6.5 Ma [10,32]. No extrusive rocks with ages less than 6.5 Ma have been found within the El Teniente deposit.  [33,34]. The Dacitic Porphyries north of the Sewell Tonalite, are mapped as a distal portion of this pluton, although they are younger [35] and have an independent origin [37]. The spatial extent of biotite breccias is projected onto this level from where they have been mapped between levels Teniente 4 and 8. Extrusive rocks of the Miocene Farellones Formation, locally referred to as the Teniente Volcanic Complex [10], are the oldest rocks exposed in the immediate area surrounding the deposit (Figure 3).
The Farellones Formation, a sequence of >2500 m of lavas, volcanoclastic rocks, dikes, sills and stocks of basaltic to rhyolitic composition, was erupted after an Early Miocene  Ma episode of crustal deformation, thickening and uplift [10,14]. The Teniente Volcanic Complex near the deposit has been correlated with the upper part of this formation and dated between 14.4 and 6.5 Ma [10,32]. No extrusive rocks with ages less than 6.5 Ma have been found within the El Teniente deposit.
The Teniente Volcanic Complex consists of tholeiitic to calc-alkaline extrusive rocks, which plot in the medium to high-K group of convergent plate boundary arc magmas [10]. Mafic rocks of the Teniente Volcanic Complex generally have higher La/Yb (4 to 9) compared with rocks of the older Coya-Machali Formation, and mafic, intermediate and silicic rocks also have higher initial 87 Sr/ 86 Sr (0.7039-0.7041; Figure 4) and lower initial ε Nd (+2.7 to +3.6; Figures 4 and 5). These differences are interpreted to represent a change from magma genesis below relatively thin continental crust during the mid-Tertiary, when the Coya-Machali Formation was generated, to conditions of thickened continental crust when the Teniente Volcanic Complex formed in the Middle to Late Miocene.

Previous Hf Isotopic Studies
Hf isotopes have previously been determined in zircons from felsic plutons both within [28] and in the vicinity [50] of the El Teniente deposit (Figure 7). The Hf isotopic composition of zircons in six samples of felsic plutons within the deposit, that range in age from 7.1 to 4.8 Ma, have overall high initial average of ε Hf(t) = +7.4 ± 1.2 and range from +6 to +10 (Table 1), which rules out significant crustal contamination [28]. There was little change in these compositions during this time, but the older four plutons (7.1 to 6.1 Ma) average ε Hf(t) = +7.8, while the younger Teniente Dacite Porphyry (5.3 Ma) has ε Hf(t) = +7.2 and the youngest latite dike (4.8 Ma) has ε Hf(t) = +6.9 (Table 1), so some small but significant temporal change, with decreasing ε Hf(t) through time, is apparent.
Hf isotopic compositions were also determined in four barren Middle and Late Miocene plutonic complexes in the vicinity of El Teniente [50]. These can be divided into an older group (~15 Ma) and a younger group (10 to 12 Ma), both groups being older than the felsic plutons within the El Teniente deposit (Figure 7). The older group is characterized by quartz monzogabbroic compositions with 54.0 to 57.5 wt. % SiO 2 , while the younger group are quartz monzodioritic with higher SiO 2 from 61.6 to 66.8 wt. %. ε Hf(t) for the older group ranges from +9.8 to +13.7 and averages +11.6 ( Figure 7), while ε Hf(t) for the younger group has a lower range of +7.6 to +8.5 and a lower average of +8.0. These differences are also associated with increasing ε Os from +32 to +160 for the older compared to the younger group. Both the Hf and Os isotopic data suggests increasing incorporation of crustal components with time, and the younger plutons, which have higher SiO 2 , are noteworthy for the presence of abundant mafic enclaves indicating intra-crustal (MASH) assimilation on a macroscopic scale. These are absent from the older plutonic group which have more primitive isotopic compositions and lower SiO 2 .  Table 4; red diamonds from [28]) and Middle Miocene barren plutons in the same general area (green squares from [50]). The size of the symbols is greater than the error in the 176 Hf/ 177 Hf values. The mantle isotopic evolution line is drawn from the average of the older group of Middle Miocene barren plutons through the Late Miocene Teniente Mafic Complex (CMET) olivine basalts and gabbros, and then to the Pliocene olivine lamprophyres and olivine-bearing basaltic andesites, ignoring the younger group of barren plutons because these have macroscopic evidence for intra-crustal assimilation [50]. The amount of subducted crust plus sediment added to the mantle to generate this line is indicated in Figure 9.
Hf and Nd isotopic compositions have also been determined for five Late Miocene to Pliocene (8.3 to ~3 Ma) plutons from the Eastern Principal Cordillera ~40 km east of El Teniente (Figure 8; [24]). These have lower more radiogenic εHf(t) (+4 to −4) and εNd (+3 to 0) than the felsic plutons in El Teniente. Crustal-contamination (MASH) processes involving assimilation of Paleozoic and Mesozoic basement rocks are evidently involved in the genesis of the igneous rocks from the Eastern Cordillera (Figure 8), both because some contain zircon crystal with older (Paleozoic) inherited cores and others have wide ranges (+1 to −4) of zircon εHf(t) values [24]. However, it has also been suggested that the low εNd values of the Pliocene olivine lamprophyres (+1.2 to +0.7) and olivine-bearing basaltic andesites (−1.1; Figures 4 and 5; Table 1; [21][22][23]) emplaced in the Western Cordillera in the vicinity of El Teniente reflect a westward propagation of the Eastern Cordillera magmatic isotopic signatures, caused by westward thrusting and material transport of the Eastern Cordillera basement under the Western Cordillera in the Pliocene (Figure 8; [24]). We disagree with this interpretation as it fails to take into account the fact that the Pliocene high-Mg olivine lamprophyres, as well as the olivinebearing basaltic andesite lavas, are likely to be mantle-derived, without significant crustal contamination. For high-Mg olivine lamprophyres in general, an origin by hydrous (6-16 wt. % H2O) melting of mantle has been suggested [51][52][53][54][55]. As discussed in more detail below, we therefore interpret the temporal change in both εHf(t) and εNd to lower more radiogenic values for Pliocene olivine-bearing mafic rocks associated with El Teniente in the Western Cordillera as reflecting a progressive change in mantle composition due to subduction of increased proportions of continental crust, and not to intra-crustal (MASH) contamination processes.  Table 4; red diamonds from [28]) and Middle Miocene barren plutons in the same general area (green squares from [50]). The size of the symbols is greater than the error in the 176 Hf/ 177 Hf values. The mantle isotopic evolution line is drawn from the average of the older group of Middle Miocene barren plutons through the Late Miocene Teniente Mafic Complex (CMET) olivine basalts and gabbros, and then to the Pliocene olivine lamprophyres and olivine-bearing basaltic andesites, ignoring the younger group of barren plutons because these have macroscopic evidence for intra-crustal assimilation [50]. The amount of subducted crust plus sediment added to the mantle to generate this line is indicated in Figure 9.
Hf and Nd isotopic compositions have also been determined for five Late Miocene to Pliocene (8.3 to~3 Ma) plutons from the Eastern Principal Cordillera~40 km east of El Teniente (Figure 8; [24]). These have lower more radiogenic ε Hf(t) (+4 to −4) and ε Nd (+3 to 0) than the felsic plutons in El Teniente. Crustal-contamination (MASH) processes involving assimilation of Paleozoic and Mesozoic basement rocks are evidently involved in the genesis of the igneous rocks from the Eastern Cordillera (Figure 8), both because some contain zircon crystal with older (Paleozoic) inherited cores and others have wide ranges (+1 to −4) of zircon ε Hf(t) values [24]. However, it has also been suggested that the low ε Nd values of the Pliocene olivine lamprophyres (+1.2 to +0.7) and olivine-bearing basaltic andesites (−1.1; Figures 4 and 5; Table 1; [21][22][23]) emplaced in the Western Cordillera in the vicinity of El Teniente reflect a westward propagation of the Eastern Cordillera magmatic isotopic signatures, caused by westward thrusting and material transport of the Eastern Cordillera basement under the Western Cordillera in the Pliocene (Figure 8; [24]). We disagree with this interpretation as it fails to take into account the fact that the Pliocene high-Mg olivine lamprophyres, as well as the olivine-bearing basaltic andesite lavas, are likely to be mantle-derived, without significant crustal contamination. For high-Mg olivine lamprophyres in general, an origin by hydrous (6-16 wt. % H 2 O) melting of mantle has been suggested [51][52][53][54][55]. As discussed in more detail below, we therefore interpret the temporal change in both ε Hf(t) and ε Nd to lower more radiogenic values for Pliocene olivine-bearing mafic rocks associated with El Teniente in the Western Cordillera as reflecting a progressive change in mantle composition due to subduction of increased proportions of continental crust, and not to intra-crustal (MASH) contamination processes. This also produced increasing crustal thickness, which led to uplift of the entire belt. Deformation amounts suggest that the Cuyania block descended immediately below the eastern part of the magmatic arc, changing composition of the base of the crust and thus changing the crustal magmatic source (MASH zone) below the Eastern Principle Cordillera (EPC). However, below the Western Principle Cordillera (WPC), mafic magmas continued to rise from the mantle without intra-crustal assimilation in a MASH zone. Between 7 and4 Ma, these mantle-derived mafic magmas did not reach the surface, but were mixed, stored and homogenized in a large magma chamber (blue) below El Teniente. Within this chamber, the felsic plutons of the deposit developed without crustal assimilation and above it the mineralized breccias in the deposit were emplaced (see Figure 10). As subduction angle flattened and magma supply from the mantle decreased this chamber solidified. After 4 Ma, high-MgO olivine lamprophyres and olivine-bearing basaltic andesites, derived from a mantle hydrated and isotopically modified by increased proportions of subducted crust and sediment (Figure 7), rose through the now crystallized magma chamber, without any intra-crustal assimilation, to form dikes and lavas. Finally, εHf(t) values ranging from +6.9 to +9.6 were determined for 8 plutonic rocks between the ages of 28.1 and 11.5 Ma from the region north of El Teniente between 32-33° S [47]. The oldest  [18][19][20]24] of crustal evolution of the Andes of central Chile and Argentina at 7 Ma (A) and 4 Ma (B). Evolution is characterized by increased shortening rates along the eastern flank of the orogenic belt after 7 Ma, where the Chilenia-Choiyoi block overrides the Cuyania terrane, producing uplift of the Frontal Cordillera. This also produced increasing crustal thickness, which led to uplift of the entire belt. Deformation amounts suggest that the Cuyania block descended immediately below the eastern part of the magmatic arc, changing composition of the base of the crust and thus changing the crustal magmatic source (MASH zone) below the Eastern Principle Cordillera (EPC). However, below the Western Principle Cordillera (WPC), mafic magmas continued to rise from the mantle without intra-crustal assimilation in a MASH zone. Between 7 and4 Ma, these mantle-derived mafic magmas did not reach the surface, but were mixed, stored and homogenized in a large magma chamber (blue) below El Teniente. Within this chamber, the felsic plutons of the deposit developed without crustal assimilation and above it the mineralized breccias in the deposit were emplaced (see Figure 10). As subduction angle flattened and magma supply from the mantle decreased this chamber solidified. After 4 Ma, high-MgO olivine lamprophyres and olivine-bearing basaltic andesites, derived from a mantle hydrated and isotopically modified by increased proportions of subducted crust and sediment (Figure 7), rose through the now crystallized magma chamber, without any intra-crustal assimilation, to form dikes and lavas. Finally, ε Hf(t) values ranging from +6.9 to +9.6 were determined for 8 plutonic rocks between the ages of 28.1 and 11.5 Ma from the region north of El Teniente between 32-33 • S [47]. The oldest sample, a dacite sill (sample Ab-143; 28.1 Ma), which has ( 87 Sr/ 86 Sr)i = 0.70661 and is likely to have assimilated crustal material, has lower ε Hf(t) of +7.8 than younger (22.2-17.7 Ma) gabbro, microgabbro and granodiorite samples which average ε Hf(t) of +8.9. All the still younger plutons have lower ε Hf and the youngest (11.5 Ma; sample Z-132) has an ε Hf(t) of +7.0. This suggest decreasing ε Hf(t) with time beginning after~17 Ma, and possibly before this time. Since these samples are all north of the current locus of subduction of the Juan Fernández Ridge, in an area where ridge subduction, decreasing subduction angle, eastward arc migration and increased rates of subduction erosion all happened earlier than at the latitude of El Teniente, we do not consider them appropriate for evaluating the temporal evolution of the mantle source of El Teniente magmatic rocks.

Samples and Methods
Twelve samples were selected for Hf-isotopic analysis (Tables 1-3). These range in age from ≥8.9 to 2.3 Ma and include (1) four samples from the Teniente Mafic Complex, including one olivine (Fo74) basalt, one gabbro and two other olivine-free basalts ( Table 2); (2) five syn-mineralization felsic intrusions that range in age from 7.1 to 4.8 Ma (Table 3), including four samples previously analyzed for Hf isotopes in zircon [28]; and (3) three post-mineralization samples that range in age from 3.8 to 2.3 Ma, including one high-MgO (7.9 wt. %; Ni~190 and Cr~390 ppm) olivine (Fo88) lamprophyre [22] and one olivine-bearing (Fo64) basaltic andesite ( Table 2). All the sample have been previously analyzed for major and trace element contents and Sr and Nd isotopic compositions [21][22][23].
Ground sample weighing~200 mg was decomposed in 15 mL screw-top Savillex Teflon ® vials with a 1:1 mixture of 27 N HF and 16 N HNO 3 at~150 • C for 5 days. The closed beakers were repeatedly placed in an ultrasonic bath during the course of digestion to enhance digestion efficiency. Each individual sample was carefully checked to make sure no visible solid residue was left in the Teflon beakers, and a clear solution was obtained for each run. The sample solutions were evaporated to incipient dryness and redissolved first by 2 mL 4 N HCl and then by 1 N HCl. The dried samples were then taken into solution by 1 mL 2 N HCl and 8 mL 6 N HCl−0.2 N HF prior to chemical separation. Hafnium was separated from the matrix elements by two columns packed with Ln-spec resin (50−100 µm) from Eichrom Industries. Major elements and REE were eluted by 3 N HCl and 6 N HCl, and then Ti and Zr were eluted using a 2 N HCl−0.1 N HF mixture. Hafnium was collected from the column with 2 N HCl and 0.2 N HF, and was then evaporated to incipient dryness and taken up by 2 N HCl−0.2 N HF. The chromatographic procedures were repeated again on the same column for further purification. Hafnium isotopes were measured on a Nu Plasma II multi-collector ICP-MS connected to an Aridus II desolvator and an enhanced sensitivity system with Ni 325-294 skimmer cone and Ni 319-646 sampler cone. The measured 176 Hf/ 177 Hf values were normalized for mass fractionation based on 179 Hf/ 177 Hf = 0.7325. Long-term mean value for the Hf standard JMC475 in the laboratory is 176 Hf/ 177 Hf = 0.282164 ± 8 (n = 55). The 176 Lu/ 177 Hf ratio was calculated from the Lu and Hf contents determined by ICP-MS techniques by Activations Laboratories (Ancaster, ON, Canada).

Results
The ε Hf(t) determined on the bulk-rock samples of the 7.1 to 4.8 Ma felsic plutons (Table 4) are similar within analytical error to those determined previously on zircons for the same plutons [28]. The good match of isotopic compositions indicates that Hf isotopes were not fractionated by any potential undigested phases. The ε Hf(t) for all twelve samples ranges from an average of +9.1 for the four samples of ≥8.9 Ma Teniente Mafic Complex basalts and gabbros down to +1.6 for the 2.3 Ma olivine-bearing basaltic andesite lava flow (Table 4), which is the last manifestation of magmatic activity in the region of El Teniente prior to the eastward migration of the volcanic arc. This implies a drop of 7.5 ε Hf(t) units over this 6.6 m.y. interval. Although Muñoz et al. [28] concluded that there was no change in ε Hf(t) among the felsic plutons emplaced between 7.1 and 4.8 Ma, the new data for samples formed over the larger~6.6 m.y. interval from ≥8.9 to 2.3 Ma indicates that there was a clear progressive decrease in ε Hf(t) with decreasing age (Figure 7).  Muñoz et al. [28] suggested that the full ε Hf(t) range of +6 to +10 of the 7.1 to 4.8 Ma felsic plutons is "identical to preceding Cenozoic barren magmatic activity in Central Chile". However, although their average ε Hf(t) value of +7.4 (Table 1) is only somewhat lower than the younger group (10-12 Ma) of barren plutons, for which macroscopic and Os isotopic evidence of upper crustal assimilation has been documented, it is significantly lower than the older group (~15 Ma; Figure 7). When the ε Hf(t) of the older~15 Ma group of barren plutons in the region, which average +11.6 [50], is taken into account, it is clear that ε Hf(t) values have not remained constant through time (Figure 7). In fact, there has been an overall decrease of 10 ε Hf(t) units for igneous rocks emplaced in the vicinity of El Teniente in the 12.7 m.y. period from the Middle Miocene to the Pliocene.

Andean Magmagenesis
The observed 10 ε Hf(t) unit decrease for igneous rocks emplaced in the vicinity of El Teniente in the 12.7 m.y. period from the Middle Miocene to the Pliocene is consistent with decreasing ε Nd(t) from +5 to −1.1, increasing 87 Sr/ 86 Sr(i) from 0.70376 to 0.70485, and increasing 206 Pb/ 204 Pb from 18.55 to 18.68 determined for these same samples (Figures 4 and 5; [21][22][23]). These variations imply increasing incorporation of continental crust through time in the magmas that formed these rocks.
The fact that the samples include mantle-derived olivine basalts and high-Mg olivine lamprophyres suggests that these continental components were incorporated into their mantle source, and not by intra-crustal contamination (MASH). This suggestion is supported by two other lines of evidence. One is that for all the rocks of variable SiO 2 content within the igneous associations of specific ages identified in the area of El Teniente, the isotopic compositions do not vary as SiO 2 increases ( Figure 5). The isotopic compositions of the different associations vary through time, but not as a function of SiO 2 . This implies that crystal fractionation processes are not combined with intra-crustal assimilation in the generation of the range of rocks from basalts through rhyolites formed during each of these successive magmatic events. Second, models of intra-crustal assimilation (MASH) for Sr and Nd isotopes (Figure 4; [21][22][23]) indicate that relatively large proportions (>40 wt. %) of crust must be assimilated to produce the Sr and Nd isotopic changes observed between mid-Tertiary basalts and Pliocene olivine lamprophyres and olivine-bearing basaltic andesites. This is considered to be an unacceptably large amount of crust to be assimilated in olivine-bearing mafic magmas.
In contrast, models of mantle source region contamination by subducted crustal components indicates that the Sr and Nd isotopic evolution of the subarc mantle source of the mafic magmas erupted in the vicinity of El Teniente between the mid Tertiary and Pliocene can be produced by the incorporation of increased, but still relatively small proportions of subducted crust and sediment (Figures 4 and 5; [3,10,[21][22][23]). Mid Tertiary Coya-Machali group mafic volcanic rocks require the addition to a primitive subarc mantle of only 1 wt. % subducted crust and sediment, while Pliocene olivine lamprophyres and basaltic andesites require the addition of between 4-6 wt. % subducted crust to produce an appropriately isotopically modified mantle source.
The difference between these intra-crustal assimilation and mantle source region contamination models reflects the relatively high Sr (450 ppm) and Nd (9 ppm) contents of a basalt relative to the crust being assimilated, thus requiring a relatively large proportion of crust to produce the appropriate isotopic leverage. In contrast the low Sr (36 ppm) and Nd (1.8 ppm) contents of the mantle allows a much smaller proportion of subducted crust and sediment to produce the same isotopic leverage (Figures 4 and 5; [3][4][5]10,[21][22][23]). For the Hf isotopic system, the same logic applies. The olivine basalts and gabbros of the Teniente Mafic Complex contain an average of 2 ppm Hf (Table 4), while average upper continental crust contains 5.3 ppm and bulk continental crust 3.7 ppm Hf [56]. In contrast, primitive mantle contains only 0.28 ppm Hf [57], so the proportion of crust needed to modify the Hf isotopic composition of the mantle is small compared to that needed to be assimilated by a basalt.
Since the mantle source region contamination model better explains the Sr and Nd isotopic evolution through time of the same olivine-bearing mafic rocks that we have analyzed for Hf isotopes, we also interpret the temporal change in their Hf isotopic compositions to reflect the isotopic evolution of their mantle source affected by increased incorporation of subducted crust and sediment. To quantitatively evaluate the amount of source region contamination to produce these Hf isotopic variations is complicated for two reasons. The first is the inherent complexities involved in subduction zone magmatism, such as the amounts of pelagic, terrigenous and tectonically eroded continental crust being subducted, and the uncertainty in the processes of transfer of these components into the overlying mantle wedge by either dehydration and volatile transport, melting or diapiric uprise of these low density materials. The second is that information is not available for the Hf contents and isotopic compositions of all these different components relevant to magmagenesis in the vicinity of El Teniente, in particular for the coastal and offshore Paleozoic and early Mesozoic crustal rocks being tectonically eroded and subducted.
In the central Mexican magmatic belt, west of where subduction erosion rates are estimated as 79-88 km 3 /km/m.y., and where Hf concentration and isotopic data are available for both pelagic and terrigenous sediment, as well as coastal continental rocks, Straub et al. [58] concluded that tectonically eroded continental crust, consisting of coastal and off-shore granodiorite plutons, was the dominant control on the isotopic composition of the calc-alkaline arc magmas, and that neither intra-crustal assimilation nor subducted trench sediment played a significant role in the genesis of these magmas. For the El Teniente rocks, we calculated a model assuming for Hf contents the average estimated values for primitive mantle and the average upper crust for the subducted crust and sediment, as well as a low crustal 176 Hf/ 177 Hf isotopic ratio similar to that determined by Straub et al. [58] for Mexican basement biotite gneisses (Table 5). This model can generate results ( Figure 9) consistent with those previously determined for the Sr and Nd isotopic systems for which more published values for the materials involved in Andean magmagenesis are available [10,[21][22][23]. These models all indicate that the addition of 1 wt. % subducted crust and sediment can produce an isotopically modified mantle appropriate for the generation of the mid Tertiary Coya-Machali basalts. Addition of progressively greater amounts of 2, 4 and 6 wt. % subducted crust and sediment is required to further modify the mantle to produce the Late Miocene Teniente Mafic Complex olivine basalts and gabbros, the Pliocene olivine lamprophyres and younger Pliocene basaltic andesites, respectively. Table 5. Parameters used for the mantle modified by subducted crust + sediment model [22]. olivine bearing mantle-derived lamprophyres [22]. Decreased degrees of mantle partial melting may also explain the somewhat higher Hf contents of the small volume of Pliocene olivine lamprophyres (>3 ppm) and basaltic andesite lavas (5 ppm) compared to the older Teniente Mafic Complex olivine basalts and gabbros (2 ppm; Table 4). Lamprophyres in general have been attributed to small degrees, relative to basalts, of partial melting of hydrated (6 to 16 wt. % H2O) mantle [51][52][53][54][55], H2O being one of the subducted components likely to increase in significance in the mantle wedge as mantle volume decreases.

Element and Isotopes
. Figure 9. 176 Hf/ 177 Hf versus 143 Nd/ 144 Nd for samples from El Teniente (Table 4; blue circles) compared to isotopic values of a primitive mantle (Table 5) modified by the addition of various proportions (red diamonds) of subducted crust and sediment.
In summary, the Hf isotopic data are consistent with Sr, Nd and Pb isotopic data, and suggest the incorporation of small but progressively increasing proportions of continental crust in the mantle source of the mafic magmas erupted in the vicinity of El Teniente between 15 and 2.3 Ma. Trench sediment recycling has been confirmed by the detection of cosmogenic 10 Be in many, including Andean arc, lavas [59][60][61], and comprehensive trace-element and isotopic studies have also provided evidence for the incorporation of fore-arc eroded crust in the subarc mantle source of both Andean [3][4][5]7,8,10,11] and arc volcanoes along other convergent plate boundaries [58,[62][63][64]. The model developed for the El Teniente area suggests that increasing rates of subduction erosion through time has resulted in variable, but nevertheless relatively small (1 to ≤6 wt. %), amounts of subducted crust having significantly affected the isotopic composition of the mantle source of the mafic magmas through time. Models for magmagenesis in other arcs, for instance the central Mexican volcanic belt, have suggested much larger amounts of recycled crust involved in the generation of both mafic and intermediate arc magmas; from >60% in conjunction with diapiric uprise and melting of subducted crust [58] to 100% resulting in the generation of andesites by bulk melting of subducted crust "relaminated" below the mantle wedge [63].

El Teniente Deposit
The giant El Teniente Cu-Mo deposit originally contained ≥100 million metric tonnes (Mt) of Cu [33,34]. Pb, S and Os isotopic data [43][44][45]65] indicate that Cu, Mo and S were derived from the associated igneous rocks. O and H isotopic and fluid inclusion data suggest that these metals were transported and deposited by high-temperature magmatic fluids [38,39,66], and a significant proportion of the Cu ore is contained with magmatic-hydrothermal breccia pipes [33][34][35][36][37]. Important  (Table 4; blue circles) compared to isotopic values of a primitive mantle (Table 5) modified by the addition of various proportions (red diamonds) of subducted crust and sediment.
What can be said unequivocally in support of this model of the isotopic evolution of the subarc mantle underlying El Teniente is that the rate of subduction erosion and amount of Paleozoic and early Mesozoic crust being subducted did increase between the mid Tertiary and Pliocene as the locus of subduction of the Juan Fernández Ridge migrated southward. Prior to impingement of the ridge on the northern end of the SVZ, subduction erosion rates may have been similar to those estimated presently for the central SVZ of 35 km 3 /km/m.y. (Figure 2; [5]). Over the last~10 m.y., as the ridge approached it current position at the northern end of the SVZ (Figures 2 and 3; [26,27]), rates of subduction erosion west of El Teniente increased to possibly as high as between 115 and 230 km 3 /km/m.y. [5]. This increased rate of subduction erosion must be considered as the most probable cause for the increasingly crustal Hf, Sr, Nd and Pb isotopic signatures in the igneous rocks associated with the El Teniente deposit through time.
A second important factor that changed over the same time period, was that the angle of subduction decreased, thus contributing, along with subduction erosion, to the~50 km eastward migration of the Andean arc in the Pliocene (Figure 3; [10,25]). Decreasing subduction angle reduced the volume of the subarc mantle wedge and resulted in a progressive decrease in the volume of magmas being produced prior to arc migration. Even if the amount of crustal components transferred from the subducted slab into the mantle wedge remained the same through time, which it did not since subduction erosion rates increased, the decreasing volume of the wedge would result in these components becoming progressively more significant within the wedge. Stern [3][4][5] and Kay et al. [10] have shown that the decreasing volume of magma erupted in the vicinity of El Teniente through time from the mid Tertiary to Pliocene is consistent with decreasing percent of mantle partial melting as the angle of subduction decreased. This resulted in increasing Sr contents, from 450 to 700-900 ppm in association with increasing 87 Sr/ 86 Sr ratios in the sequence of mid Tertiary basalts to Pliocene olivine bearing mantle-derived lamprophyres [22]. Decreased degrees of mantle partial melting may also explain the somewhat higher Hf contents of the small volume of Pliocene olivine lamprophyres (>3 ppm) and basaltic andesite lavas (5 ppm) compared to the older Teniente Mafic Complex olivine basalts and gabbros (2 ppm; Table 4). Lamprophyres in general have been attributed to small degrees, relative to basalts, of partial melting of hydrated (6 to 16 wt. % H 2 O) mantle [51][52][53][54][55], H 2 O being one of the subducted components likely to increase in significance in the mantle wedge as mantle volume decreases.
In summary, the Hf isotopic data are consistent with Sr, Nd and Pb isotopic data, and suggest the incorporation of small but progressively increasing proportions of continental crust in the mantle source of the mafic magmas erupted in the vicinity of El Teniente between 15 and 2.3 Ma. Trench sediment recycling has been confirmed by the detection of cosmogenic 10 Be in many, including Andean arc, lavas [59][60][61], and comprehensive trace-element and isotopic studies have also provided evidence for the incorporation of fore-arc eroded crust in the subarc mantle source of both Andean [3][4][5]7,8,10,11] and arc volcanoes along other convergent plate boundaries [58,[62][63][64]. The model developed for the El Teniente area suggests that increasing rates of subduction erosion through time has resulted in variable, but nevertheless relatively small (1 to ≤6 wt. %), amounts of subducted crust having significantly affected the isotopic composition of the mantle source of the mafic magmas through time. Models for magmagenesis in other arcs, for instance the central Mexican volcanic belt, have suggested much larger amounts of recycled crust involved in the generation of both mafic and intermediate arc magmas; from >60% in conjunction with diapiric uprise and melting of subducted crust [58] to 100% resulting in the generation of andesites by bulk melting of subducted crust "relaminated" below the mantle wedge [63].

El Teniente Deposit
The giant El Teniente Cu-Mo deposit originally contained ≥100 million metric tonnes (Mt) of Cu [33,34]. Pb, S and Os isotopic data [43][44][45]65] indicate that Cu, Mo and S were derived from the associated igneous rocks. O and H isotopic and fluid inclusion data suggest that these metals were transported and deposited by high-temperature magmatic fluids [38,39,66], and a significant proportion of the Cu ore is contained with magmatic-hydrothermal breccia pipes [33][34][35][36][37]. Important questions concerning the genesis of El Teniente, and any other giant orthomagmatic ore deposit, are: (1) Was the large amount of metal in the deposit derived from a large amount of magma or a unique metal-rich magma?; and (2) How were the metals transferred from the magma into the host rocks?
Muñoz et al. [28] conclude that the high initial ε Hf(t) values (total average of +7.4) of the syn-mineralization felsic plutons "rules out involvement of any significant crustal contamination in the genesis of the El Teniente magmas". We agree, at least with regard to intra-crustal assimilation (MASH). However, they then go on to suggest that, based "entirely on the near constant Hf isotopic composition shown by the Cenozoic igneous of the region", they also "envision that Cenozoic magmatism in central Chile originates from a MASH reservoir". We strongly disagree, given the clear evidence that the ε Hf(t) values of Cenozoic magma, in particular mantle-derived mafic magmas erupted in the region, were not at all constant, but progressively decreased significantly with time.
We suggest instead that the near constancy of the Hf isotopic composition of the 7.1 to 4.8 Ma felsic plutons in the El Teniente deposit is a result of mixing, storage and homogenization (but not crustal assimilation) of mantle-derived magmas in a large vertically zoned upper crustal magma chamber ( Figure 10) which developed as a result of compressive deformation after 7 Ma (Figure 8). These mafic mantle-derived magmas provided heat, water, S and metals into the base of this magma chamber [67], which became vertically stratified with more hydrous and more silica-rich magmas closer to the surface. As the input of mantle-derived magmas decreased during the Late Miocene, and uplift and erosion brought the roof of the chamber closer to the surface [12,13], the upper part of the chamber crystallized and differentiated into the magmas that formed the felsic plutons. Although the Hf isotopic composition of the mantle-derived magmas may have been evolving over this same time period (Figure 7), the large volume of the magma chamber into which these magmas mixed limited the variations observed within the progressively younger felsic plutons, although some small decrease in their ε Hf(t) values is apparent. Once this large magma chamber had fully crystallized, mantle-derived olivine lamprophyre dikes were able to penetrate up to the surface, and these clearly exhibit lower ε Hf(t) resulting from increased source region contamination of the mantle by subducted components. Small barren plutons of the Young Plutonic Complex (5.5 to 6 Ma; Figure 5; [10]), which occur outside the area of the deposit, also have more radiogenic characteristics than the felsic porphyries of similar age within the deposit, presumably because they have not formed above such a large, long-lived magma chamber and their isotopic compositions reflect more directly that of the mantle-derived magmas that they were derived from. This is consistent with the suggestion that the difference between barren plutons and the mineralized deposit reflects the size and longevity of the magma chambers that lead to their formation. mantle derived magmas were added continuously to the base of the chamber and subsequently concentrated near the roof of the chamber by volatile transport. Clearly, because of increased source region contamination by subducted materials the mantle wedge did also became progressively more hydrated, ultimately leading to the generation of lamprophyric magmas, which require 6-16 wt. % H2O, instead of olivine basalts, which require only 2 wt. % H2O [52][53][54][55]. Increasing H2O content of the mantle-derived magmas through the time that they were mixed into the overlying magmas chamber may have resulted in more oxidizing conditions within the chamber, enhancing the process of volatile transfer of S and metals towards the roof of the chamber ( Figure 10; [49,51,71]). Figure 10. Model for the multistage development of the El Teniente deposit [23,33,34,39]. The main features of the model include (1) a large, long-lived (≥3 Myr) open-system magma chamber, crystallizing at ≥4 km depth, fed from below by mantle-derived mafic magmas and exsolving mineralrich high-temperature magmatic fluids through its roof to produce the large breccia pipes that are prominent features this deposit; (2) decreasing magma supply in the Late Miocene and Pliocene as subduction angle decreases, leading to crystallization and solidification of this chamber; (3) progressive uplift and erosion that enhances this crystallization and solidification process and results in telescoping of different types of breccia and igneous rocks [12,13]; (4) progressive igneous differentiation of the magma chamber associated with crystallization and volatile loss, generating felsic porphyries that intrude previously mineralized rocks above the chamber. No coeval volcanism is known to have occurred during mineralization, but once the chamber solidified, postmineralization mantle-derived olivine lamprophyre dikes were emplaced. Older studies of El Teniente suggested that Cu was transferred into the host rocks in association with the intrusion of the Teniente Dacite Porphyry, but new age dates indicate that much of the mineralization precedes this event, and that where this essentially barren porphyry intrudes outside the deposit, it creates only alteration but not mineralization. Instead we have suggested that the Cu and S were transferred into the host rocks of the deposit from the roof of the large magma chamber by exsolution of high-temperature metal-rich fluids generating multiple large magmatichydrothermal mineralized breccia complexes such as the central Braden Pipe, and El Teniente is thus better classified as a megabreccia deposit (Figure 10; [23,33,34,38,39]. Vry et al. [37] proposed that Figure 10. Model for the multistage development of the El Teniente deposit [23,33,34,39]. The main features of the model include (1) a large, long-lived (≥3 Myr) open-system magma chamber, crystallizing at ≥4 km depth, fed from below by mantle-derived mafic magmas and exsolving mineral-rich high-temperature magmatic fluids through its roof to produce the large breccia pipes that are prominent features this deposit; (2) decreasing magma supply in the Late Miocene and Pliocene as subduction angle decreases, leading to crystallization and solidification of this chamber; (3) progressive uplift and erosion that enhances this crystallization and solidification process and results in telescoping of different types of breccia and igneous rocks [12,13]; (4) progressive igneous differentiation of the magma chamber associated with crystallization and volatile loss, generating felsic porphyries that intrude previously mineralized rocks above the chamber. No coeval volcanism is known to have occurred during mineralization, but once the chamber solidified, post-mineralization mantle-derived olivine lamprophyre dikes were emplaced. Muñoz et al. [28] further suggest that, as a result of crustal thickening, dehydration melting in the deep crustal hot zone [68] within which the MASH reservoir was located "probably influenced the fertility of the magmas by increasing the melt component derived from this process relative to the component derived from primary" mantle melts. They do not detail what specific characteristics make a magma more "fertile", but they note that dehydration melting leads to the development of water undersaturated magmas, which are less, not more likely to be able to exsolve mineral-rich magmatic fluids.
In contrast, we conclude that the mantle-derived olivine basalts emplaced in the Teniente Mafic Complex prior to the development of the large magma chamber below the deposit, and the mantle-derived olivine lamprophyres which cut through the deposit after this magma chamber fully solidified, have changed their isotopic compositions due to increased contamination of the mantle source region with subducted components, but we do not observe [22,23] any significant change in the Cu contents between the pre-mineralization olivine basalts (~100 ppm) and the post-mineralization olivine lamprophyres (~100 ppm). Therefore, we do not believe that source region contamination generated a uniquely "fertile" magma with respect to Cu content, although it did increase the H 2 O content of the mantle. We have suggested instead that the large amount of Cu in the deposit reflects a large amount of magma in the chamber above which the deposit developed [38,39,69,70]. To produce the ≥100 × 10 6 tonnes of Cu in the deposit requires a batholith-sized (≥600 km 3 ) amount of magma with 100 ppm Cu. What made the deposit giant was that over the~3 m.y. period that a magma chamber existed below the developing deposit, Cu, S heat and water from mantle derived magmas were added continuously to the base of the chamber and subsequently concentrated near the roof of the chamber by volatile transport. Clearly, because of increased source region contamination by subducted materials the mantle wedge did also became progressively more hydrated, ultimately leading to the generation of lamprophyric magmas, which require 6-16 wt. % H 2 O, instead of olivine basalts, which require only 2 wt. % H 2 O [52][53][54][55]. Increasing H 2 O content of the mantle-derived magmas through the time that they were mixed into the overlying magmas chamber may have resulted in more oxidizing conditions within the chamber, enhancing the process of volatile transfer of S and metals towards the roof of the chamber ( Figure 10; [49,51,71]).
Older studies of El Teniente suggested that Cu was transferred into the host rocks in association with the intrusion of the Teniente Dacite Porphyry, but new age dates indicate that much of the mineralization precedes this event, and that where this essentially barren porphyry intrudes outside the deposit, it creates only alteration but not mineralization. Instead we have suggested that the Cu and S were transferred into the host rocks of the deposit from the roof of the large magma chamber by exsolution of high-temperature metal-rich fluids generating multiple large magmatic-hydrothermal mineralized breccia complexes such as the central Braden Pipe, and El Teniente is thus better classified as a megabreccia deposit ( Figure 10; [23,33,34,38,39]. Vry et al. [37] proposed that such breccias formed in association with the intrusion of the small dacite and diorite dikes to the northeast of the Braden pipe ( Figure 6). However, we suggest instead that the breccias formed first and the porphyries intruded into the zone of weakness created by the breccias, since some of these porphyry bodies, such as Porphyry A ( Figure 6) contain clasts of the breccias they intrude, and others of these small porphyries have no associated breccias surrounding them.

Conclusions
The Hf isotopic data are consistent with the Sr, Nd and Pb data, and imply the incorporation of small, but progressively increasing proportions of continental crust into the mantle source of the mafic magmas erupted in the vicinity of El Teniente between ≥15 and 2.3 Ma.
The felsic plutons in the deposit formed by mixing, storage and homogenization, combined with crystal-liquid differentiation, but not crustal assimilation (not MASH), in a large, long-lived, vertically zoned magma chamber that developed due to compressive deformation and persisted between the period~7 to 4.6 Ma.
Progressively more hydrous and radiogenic mantle-derived mafic magmas feed this chamber from below, providing heat, H 2 O, S and metals, but no unique "fertile" Cu-rich magma was involved in the formation of the deposit.
As the volume of these mantle-derived magmas decreased from the Late Miocene into the Pliocene, the chamber crystallized and solidified, producing both felsic plutons and large magmatic-hydrothermal breccias that emplaced Cu and S into the older (≥8.9 Ma) mafic Teniente Mafic Complex host rocks of the deposit. The deposit is, therefore, best classified as a megabreccia deposit.