The Volcanic Evolution of Cerro Uturuncu : A High-K , Composite Volcano in the Back-Arc of the Central Andes of SW Bolivia

Cerro Uturuncu, southwest Bolivia, is a high-K, calc-alkaline, composite volcano constructed upon extremely thick continental crust approximately 125 km behind the arc-front of the Andean Central Volcanic Zone (CVZ). Eruptive activity occurred between 890 271 ka in a single phase of volcanism lasting ~620,000 years. The edifice consists of a central cone and several flank vents where dacitic and andesitic lava flows and domes erupted. Volumes of individual eruptive units range from 0.1 to ~10 km3; the composite volume of Uturuncu is ~89 km3. In this paper, we present new field, petrographic, and geochemical data in an effort to understand the volcanic and magmatic evolution of Uturuncu. Lava flows and domes have a restricted range in whole rock compositions ranging from 61 wt% 67 wt% SiO2; magmatic inclusions contained within these units have a larger range from 53 wt% 64 wt% SiO2. Typical phenocryst assemblages are plagioclase > orthopyroxene > biotite >> quartz and Fe-Ti oxides. Pb isotope ratios are characteristic of the southern CVZ by containing high 207Pb/204Pb and 206Pb/204Pb and moderate to high 208Pb/204Pb. Sr and Nd isotope ratios indicate that Uturuncu magmas were modified by high 87Sr/86Sr and low 143Nd/144Nd felsic basement lithology during magma migration and differentiation. In all eruptive units, there is petrographic and geochemical evidence for magma mixing and mingling. In this regard, magma mixing and mingling is considered to be responsible for the small range in lava flow and dome compositions throughout the eruptive history of the center.


Introduction
Volcanism in the Central Volcanic Zone (CVZ) of the Andes is a prime example of active continental arc magmatism.Young (e.g., <6 Ma) volcanic rocks in the CVZ are dominantly andesitic to dacitic in composition with geochemical and petrographic features indicative of extensive crustal contamination during magma differentiation (e.g., [1]- [3]).Previous studies of CVZ volcanic rocks have greatly contributed to our understanding of the regional setting of the magmas [4] [5]; and references therein in addition to ignimbrite forming eruptions [6]- [9]; among others and volcanism along the modern arc-front [3] [10].However, comprehensive studies of back-arc composite volcanoes are scarce, due in part, to their rarity and remote locations.As a result, studies of these systems are typically reconnaissance in nature and associated with regional studies [4] [5] [11]- [13].
The lack of comprehensive studies on individual back-arc composite volcanoes in the CVZ prompted a detailed field, petrologic and geochemical study of the back-arc center Cerro Uturuncu (6010 m; 22˚15'S, Figure 1).The objectives of this study are to discuss the volcanic evolution of an individual back-arc composite volcano and broadly assess petrogenetic processes affecting magma composition in the back-arc of the CVZ.This study of Uturuncu results from a reconnaissance petrologic and geophysical investigation by Sparks et al. [16].These authors concluded Uturuncu exhibits stages of volcanism separated by periods of quiescence common at other intermediate composite centers in the CVZ.Currently, the system is in a period of dormancy as the most recent eruption occurred at ~270 ka.Relatively long periods of quiescence are common at composite volcanoes as they been observed at other intermediate composition systems in the CVZ and elsewhere such as at Soufrière Hills Volcano, Montserrat [3] [17] [18].In this paper, we build on the reconnaissance investigation of Sparks et al. [16] by furthering discussion of the volcanic history, field relationships, and petrography of volcanic rocks at Uturuncu.We also present a broad geochemical model to account for the restricted compositional diversity of magmas erupted at the volcano.More detailed discussion of the petrology and petrogenesis of Uturuncu magmas will ensue in a forth-coming paper.

Figure 1.
Map showing location of Andean Central Volcanic Zone (CVZ).Shaded area shows the region where crustal thickness exceeds 60 km [14] [15]; stippled region illustrates distribution of Quaternary volcanic rocks.Modified from Feeley and Hacker [2].

Regional and Tectonic Setting
The Andes are considered the classic example of a modern Cordilleran-type orogen formed by long-term subduction of oceanic lithosphere beneath continental lithosphere.The central Andes represents the type locality for this process owing to the great width of the orogen, immense crustal thickness (up to 70 -80 km) [19] and high elevations (~4 -6 km above sea level) [19] that occur over a vast area (Figure 1).The region composes one of the youngest and largest active silicic volcanic provinces with recent caldera formation.The CVZ contains ~20 calderas that have erupted ignimbrites (large ash flow sheets) less than 10 Ma and over 1100 late Cenozoic volcanic edifices [6] [7] [20].
The southern central Andes between 21˚S -24˚S is divided into three NW to SE trending geologic provinces based on structural similarities [19].From east to west, these are the Eastern Cordillera (Cordillera Oriental), the Altiplano, and the Western Cordillera (Cordillera Occidental).The Eastern Cordillera comprises a classic thinskinned fold and thrust belt.The Altiplano is a broad plateau where undeformed late Miocene and younger ignimbrites overlie variably deformed mid-Miocene and older volcanic and sedimentary rocks.The Western Cordillera consists of a westward dipping monocline and the modern volcanic arc-front [19] [21].Depths to the subducting slab at these latitudes range from 130 km at the arc-front to ~250 km along the eastern boundary of the Altiplano [22] [23].
Late Cenozoic to modern volcanic rocks in the CVZ is divided into three groups based on similarities in composition and eruptive style [21] [24].First, between 21˚S -24˚S, large-volume, regionally extensive ignimbrites erupted from large caldera complexes on the Altiplano and Western Cordillera [6] [7] [25].These rocks are calc-alkaline, homogeneous dacites to rhyolites.Second, basaltic andesitic to dacitic lava flows erupted from 23 Ma to the present.At 21˚S -22˚S, the largest volumes of these magmas erupted from composite volcanoes forming the peaks of the Western Cordillera.Baker and Francis [12] estimated the Western Cordillera contains ~3000 km 3 of these lavas between 21˚S and 22˚S latitude.Lavas associated with composite volcanoes extend for ~200 km eastward onto the Altiplano, although volumes decrease sharply to <800 km 3 .Uturuncu is associated with this group.Third, volumetrically minor alkali basalts erupted from small monogenetic centers to the east of the arc-front on the Altiplano north of 21˚S [13] [26].

Methods
New whole rock samples of Uturuncu volcanic rocks were analyzed for major and trace element compositions at Washington State University, Pullman.Major and trace element compositions for 121 samples were acquired by X-ray fluorescence spectrometry (XRF) on a ThermoARL Advant'XP+ automated sequential wavelength spectrometer.Methods and errors for XRF analyses are described in Johnson et al. [27].Additional trace element compositions, including the rare earth elements, for 45 samples, were acquired by inductively coupled plasma mass spectrometry (ICP-MS) on an Agilent 7700 ICP-MS.Methods and errors for trace element analyses are presented in detail in Jarvis [28].Whole rock Nd and Sr isotopic ratios, for 30 samples, were acquired by thermal ionization mass spectrometry (TIMS) on a VG Sector 54, analyzed by five Faraday collectors in dynamic mode at New Mexico State University, Las Cruces.Pb isotopes were analyzed on a Thermo Finnigan Neptune multi-collector ICP-MS equipped with nine Faraday collectors and an ion counter at the University of California at Santa Barbara.Michelfelder et al. [29] describe in detail the methods for radiogenic isotopic analyses.Geologic mapping was performed using commercially purchased ASTER satellite imagery and field mapping.

Geology and Eruptive History of Cerro Uturuncu
In this section, we summarize the volcanic and petrographic features of the lava flows and domes at Cerro Uturuncu.Cerro Uturuncu is classified as a composite volcano [20] resulting from volcanism over ~620,000 years [16].Construction of the edifice of Cerro Uturuncu took place the back-arc of the CVZ approximately 125 km behind the modern arc front in the Altiplano region (Figure 1).The volcano is built upon a series of 5 Ma to 1 Ma ignimbrites erupted from large calderas complexes in the area [8] [9] [30] [31].The eruptive history of Uturuncu is dominantly effusive with no evidence of explosive activity.Eruptive activity at Uturuncu predates the last glacial episode (~11,000 ky [26]) in the area as evident by glacial valleys incising lava flows and the reformed southern and northern flanks mantled in places by moraines (Figure 2).Glacial incision has allowed for access to the potentially oldest lava flows erupted from Uturuncu.Active fumaroles producing hydrothermally altered volcanic rocks are present in the saddle between the two peaks of Uturuncu.Hydrothermal alteration, mining of altered rocks for sulfur and erosion of the stratigraphically youngest volcanic rocks has removed the interpreted location of the central vent.The presence of fumaroles suggests current activity beneath Cerro Uturuncu.Pritchard and Simmons [32] identified a 70 km diameter regional uplift centered on Uturuncu.This uplift, suggested to be the result of diapiric rise of new magma into a shallow chamber beneath the center [33], could lead to new activity at Uturuncu in the future.
Relatively few radiometric ages are available for Uturuncu lavas and domes.The available ages suggest two stages of volcanism between 890 -271 ka [16] and a period of quiescence from 540 -470 ka.The stages described by Sparks et al. [16] are compositionally and petrographically similar, though erupted volumes of vol-

Morphology and Volume of Cerro Utrunucu
Considered dormant, Cerro Uturuncu has not been the focus of intense study until recently.Fernandez et al. [34], Kussmaul et al. [11], Baker and Francis [12] and Hildreth and Moorebath [35] included Uturuncu in regional studies but did not investigate the volcanic history.Reconnaissance geologic mapping performed by Sparks et al. [16] identified ~50 lava flows and domes with volumes ranging from 0.1 -10 km 3 .Based on new mapping (Figure 2), as well as geochemical and stratigraphic criteria, the eruptive history of Uturuncu represents a single stage of activity spanning ~620,000 years with an additional 45 lava flows and 10 domes identified as part of this study.Newly identified lava flows and domes observed during mapping and verified geochemically are stratigraphically between the 540 ka and the 470 ka lava flows.Volumes of these newly identified lava flows and domes are similar to those reported by Sparks et al. [16] with a total estimated volume of ~89 km 3 .The suggested repose interval based on radiometric age dates and the number of erupted units is between 6000 and 8000 years.
The predominant rock types at Cerro Uturuncu are medium-grey, blocky to platy, orthopyroxene, biotite andesites and dacites.Dacites of relatively uniform compositions were the dominant magmas erupted, although significant volumes of andesite lava flows were erupted throughout the life span of the volcano.Flow fronts of lava flows range in thickness from <5 m to over 200 m thick.Widths of flows vary with slope and are wider than thick for gentler slopes (Figure 2 & Figure 3).Some lava flows extend for more than 10 km and traced back to the central vent (Figure 2), while smaller volume flows rarely extend more than 2 -3 km from the vent.Many flows have internal flow folds and are autobrecciated at the terminus showing several meters of oxidation (Figure 3).
Domes occur throughout the volcano and suggest the activity was not restricted to a central vent, though lava flows cannot be trace to these flank vents.Piles of glassy, prismatically jointed blocks making up the exterior walls identified domes.Rare exposures of the interiors of domes are vesiculated ranging from 8% -23% of the total volume (Figure 3).Vesicles ranged in size from millimeter to centimeter scale.

Petrography
Mineralogically and petrographically, Uturuncu lava flows and domes are porphyritic to seriate.Domes are occasionally hiatal.Some lava flows have ophitic and poikilitic textures.Glomerocrysts containing orthopyroxene (OPX) and plagioclase common with occasional olivine observed in some lava flows.Modal compositions of Uturuncu volcanic rocks are variable in crystal content, but similar in mineral assemblage (Figure 4).Crystal contents range from 30% -53% total volume for all lava flows and domes with domes containing the highest crystal volumes.Phenocryst phases include plagioclase (An 42-94 ) > OPX (En 45-83 ) > biotite >> quartz, Fe-Ti oxides, and trace clinopyroxene (CPX), hornblende and olivine in micro-inclusions (Figure 4).
Plagioclase and OPX phenocrysts exhibit multiple textures.In the same eruptive unit, plagioclase phenocrysts textures include a combination of sieving, dissolution surfaces, growth zones, and clean crystals (Figure 5).Sieving occurs in both cores and rims of phenocrysts and dissolution surfaces commonly separate clear or sieved cores from zoned rims (Figure 5(B) & Figure 5(C)).Large plagioclase phenocrysts (>0.5 mm) are interpreted as xenocrysts because they contain cores riddled with abundant, irregular shaped glass inclusions and euhedral to subeuhedral rims with oscillatory chemical zoning from core to rim (An core = 46 -92; An rims = 44 -88).Nearly all OPX crystals show some zoning.Chemically unzoned OPX are only associated with glomerocrysts and in groundmass crystals (Figure 5(A) & Figure 5(D)).OPX crystals are Ca-poor at (Wo 06 or less) and represent two zoned populations [16].The first contains Mg-poor cores (En 55 ) and reverse zoned rims (En 65-75 ), and the second population contains Mg-rich cores (En 65-80 ) and normally zoned rims (En 50-705 ).Reaction rims of Fe-Ti oxides and OPX typically surrounds biotite phenocrysts over 0.5 mm.Quartz and olivine phenocrysts, and microphenocrysts typically exhibit resorption textures suggesting they are most likely xenocrysts or antecrysts.
Groundmass of the lava flows and domes are hyalopilitic with mineral phases similar to phenocryst phases.Additional phases in the groundmass include glass, trace zircon, apatite, and occasional olivine.Groundmass plagioclase (An 45-85 ) crystals do not exhibit sieving seen in phenocrysts.Groundmass plagioclase and OPX (En 46-75 ) crystals occasionally show zoning and dissolution surfaces.Vesicles comprise approximately 8% -23% of the mode in domes interiors and 0% -12 % in lava flows.

Magmatic Inclusions
A common occurrence in Uturuncu volcanic rocks is the presence of magmatic inclusions.Individual lava flows vary greatly in the volume of magmatic inclusions ranging from ~0% -4%.Three general populations of inclusions exist in Uturuncu lava flows and domes.Mineralogically and petrographically inclusions exhibit similar modes to the lavas and domes.With the exception of inclusions defining population 3, inclusions typically lack quartz and biotite and contain a higher volume of CPX and olivine compared to the host lava flows and domes.
The first population of inclusions is ellipsoidal with a porphyritic to hiatal texture and contains a thick, vesicle-free, glassy rim with vesiculated interiors.These inclusions are characterized by euhedral to subhedral parageneitcally early phenocrysts of OPX and plagioclase in varying proportions.Plagioclase phenocrysts range in The second population of inclusions are non-vesiculated basaltic andesite clots, containing plagioclase > OPX > Fe-Ti oxides > olivine > amphibole ± clinopyroxene.Hiatal to intergranular crystal sizes dominate this population of inclusions.The second population of inclusions lack glassy rims in contact with the host lava flow or dome as described above.The contacts of the inclusions are angular and defined by crystal boundaries.Crystal sizes range from microlite-size to crystals ~1.5 mm along the long axis.Inclusion size varies from micro-inclusions (~3 mm) to ~10 cm in diameter.Hand specimen size samples of this inclusion population was only observed in two lava flows, although these inclusions may have been overlooked due to their similarity in color to the host lavas when weathered.Two additional lava flows and one dome contain micro-inclusions of this population.
Non-vesiculated micro-inclusions are distinct from the glomerocrysts observed lava flows and domes due to their crystal size and presence of rare amphibole.OPX phenocrysts are similar in size to the plagioclase though Vesiculated and non-vesiculated magmatic inclusions (population 1 and 2) are common features observed in CVZ volcanic rocks as well as at other intermediate composition volcanic systems around the world [21] [36]- [38].Magmatic inclusions, with or without vesiculation, represent blobs of mafic magma quenched in cooler, more silicic magma [36] due to the processes of magma mixing and mingling.This is our interpretation of inclusion populations 1 and 2 observed in Uturuncu lavas and domes.
The third population of inclusions observed does not exhibit the same mineralogy as the host lavas and domes.This population is dominantly hypidiomorphic to allotriomorphic granular with trace to rare glass (0% -0.5%).Plagioclase and quartz dominate the mode of these inclusions with trace OPX and biotite.Plagioclase crystals are clear with no growth zones or dissolution surfaces observed.Quartz crystals show resorption along the corners of the crystals.Margins of these inclusions are angular and defined by the crystal boundaries.Textures observed in this population of inclusions suggest the same process does not produce these inclusions as populations 1 and 2. It is likely that inclusions in population 3 are xenoliths.
Sparks et al. [16] described two xenolith populations have a similar petrographic description to inclusions in population 3 of this study.These authors interpreted observed quartz-rich xenoliths to represent crustal partial melting of rocks trapped in a MASH zone [35].It is our interpretation that xenoliths observed in this study represent basement rocks interacting with Uturuncu magmas during migration and differentiation.These xenoliths may represent the residue of a granodiorite crustal contaminant beneath Uturuncu.
A fourth inclusion population described by Sparks et al. [16] and was not observed in this study.These authors describe a noritic composition xenolith population found in multiple lava flows sampled in the study.These xenoliths contain mineral assemblages and similar mineral compositions to those observed in the lava flows and domes [16].The adcumulate textures observed in these xenoliths suggested to these authors that these xenoliths represent adcumulates of early-formed crystals that segregated from hotter, less evolved residual magma related to the Uturuncu magmas [16].

Major and Trace Elements
Figures 6-8 illustrate major-and trace-element compositions of lava flows, domes and inclusions.Table 1 and Table 2 present representative analyses of these data.These compositions define a high-K, calc-alkaline suite (Figure 6).SiO 2 ranges from 61 wt% -67 wt% for the lava flows and domes, 53 wt% -64 wt% for the magmatic inclusions (population 1 & 2), and 78 wt% -79 wt% for the xenoliths (inclusion population 3).For the suite as a whole, contents of CaO, FeO * , MgO, TiO 2 , MnO, Yb, Sr and Cr decrease, and concentrations of Na 2 O, K 2 O, P 2 O 5 , La, Zr, Ba and Rb increase with increasing SiO 2 concentrations (Figure 7 & Figure 8).Compared to other centers in the CVZ [29], higher concentrations of are observed in LIL-element concentrations for Uturuncu volcanic rocks.
Consistent with petrography, xenolith samples do not exhibit the same geochemical trends observed in the lavas, domes, or inclusions (Figure 7 & Figure 8) suggesting these xenoliths are trapped residue in a MASH zone or represent a crustal contaminant.Compared to other volcanic rock compositions at Uturuncu, compositions of the xenoliths are significantly lower for elements not compatible with quartz or plagioclase (Figure 7 & Figure 8).
The compositional variation and trends observed versus SiO 2 highlight the petrologic processes affecting the evolution of Uturuncu magmas.The linear trends observed in all major element compositions and the divergent linear trends in MgO, La, Sr, Y and Yb suggest multiple components in magma evolution.These linear trends reflect magma mixing as a dominant process in the generation of the limited compositional diversity observed of lavas and domes at Uturuncu.Magmatic inclusion compositions are an extension of this trend at lower SiO 2 contents.9 illustrate radiogenic isotopic ratios for 87 Sr/ 86 Sr, 143 Nd/ 144 Nd, and Pb isotopes of lavas and domes, inclusions and xenoliths.The range in 87 Sr/ 86 Sr and 143 Nd/ 144 Nd isotopic ratios for Uturuncu rocks is small compared to the range observed across the CVZ (Figure 9(B)) [4] [37].In comparison with other individual CVZ composite cones, Uturuncu rocks have higher 87 Sr/ 86 Sr ratios and lower 143 Nd/ 144 Nd ratios and are more similar in comparison to ignimbrites erupted from Cerro Panizos on the eastern Altiplano or the 5 Ma to 1 Ma Altiplano ignimbrites Uturuncu overlies (Figure 9(B)) [30] [39].When compared to previously report isotopic analyses in reconnaissance or regional studies including volcanic rocks from Uturuncu, this study reports more isotopic variation [11] [12].The 87 Sr/ 86 Sr isotopic compositions for the inclusions overlap the lowest ratios measured in the host lavas and domes (Figure 9

Petrogenesis of Cerro Uturuncu Magmas
In this section, we describe a simple working model to explain the limited compositional diversity observed in Uturuncu volcanic rocks.The intent of the proposed model is not a rigorous description of the petrogenesis of the magmas nor is the variables and compositions selected for assimilation-fractional crystallization and mag- mamixing models intended to quantify the magmatic system beneath Uturuncu.The intent is to constrain the processes affecting compositional diversity of the rocks and serve as a point of reference in future studies of Uturuncu volcanic rocks.[40] from the southern CVZ.We selected this composition because of petrographic similarities to a sillimanite-gneiss xenolith described by Sparks et al. [16] from Uturuncu lavas.The upper crustal contaminant used in the models is Paleozoic granite from Sierra de Moreno described by Lucassen et al. [41].These granites are petrographically similar to the xenolith observed in this study.We contaminate two primitive basalt compositions with lower and upper crustal contaminants: one from Davidson and de Silva [13] for basalts from the central Altiplano, and a second from Hernando et al. [42] for pre-caldera basalts from the Payún Matrú volcanic field in western Argentina's Southern Volcanic Zone.Table 3 presents concentrations and isotopic ratios of the contaminants and basalts.
Model curves labeled AFC (D = 1.7) and AFC (D = 1.5) were constructed to simulate the effects of differentiation on two potential parental magmas compositions contaminated by a highly evolved crustal contaminant (Figure 10).Bulk distribution coefficients of 1.5 and 1.7 are high for andesites, but the decrease of Sr content with increasing SiO 2 content and in increase of Rb/Sr ratios with increasing 87 Sr/ 86 Sr ratios suggest Sr was compatible during differentiation.The curve AFC (D = 0.1) was constructed to simulate the effects of differentiation in the lower crust compared to curves for the upper crust.A bulk distribution coefficient of 0.1 is considered appropriate for this composition based on regional studies that have shown an increase in Sr/Y ratio with a decrease in 87 Sr/ 86 Sr ratios reflects incompatibility of Sr in garnet-stable/plagioclase-unstable-hybridized crust [29] [44] [45].
87 Sr/ 86 Sr and 143 Nd/ 144 Nd ratios for Uturuncu lavas and domes create an array consistent with assimilation of old radiogenic basement rocks with relatively non-radiogenic Nd and radiogenic Sr. High Rb/Sr ratios suggest these basement rocks were also felsic in composition as opposed to noritic as could be interpreted from the presence of noritic composition xenoliths described by Sparks et al. [16].High silica, quartz-feldspathic xeno-  3.

Volcanic Hazards
Cerro Uturuncu has historically been inactive and glacial erosion suggests that the volcano is in a period of dormancy.However, new geophysical data has shown a concentric uplift centered on Uturuncu for the past 20 years with the potential to initiate future eruptions [32] [33].According to the spatial distribution of previously erupted lava flows, it is difficult to predict the likeliness of where a lava flow will occur during a future potential eruption.In addition, the injection of new magma beneath the central vent could cause instability of the highly altered hydrothermal rocks composing the summit potentially leading to partial collapse of the vent and debris avalanches towards the north.This event has precedence in other CVZ composite cones [10] [21] [26], but has not been described in volcanic history of Cerro Uturuncu.These events at other intermediate centers have also been associated with pyroclastic flows and volcanic blasts [46]- [48].The only village within the possibly affected area is Quetena Chico approximately 30 km to the northeast potentially affecting the community with ash fall in the event of one of these more explosive eruptions.

Conclusions
In this study we combine field, petrographic and geochemical data to describe the volcanic evolution of Cerro Uturuncu, at 22˚15'S, in the CVZ.The volcano has erupted uniform, crystal-rich orthopyroxene, biotite dacite and andesite lavas for ~620,000 years.Hydrothermal activity and glaciation resulted in the removal of the main vent suggesting a long period of dormancy since the last eruption.Field relationships between the lava flows and domes suggest eruptions are effusive with an average repose interval between 6000 and 8000 years.Sparks et al. [16] identified 50 lava flows; we have identified an additional 45 lava flows and 10 domes with the total number of erupted units to 105.Erupted volumes ranged between 0.1 km 3 and ~10 km 3 per eruption with a total volume erupted of ~89 km 3 .Construction of the edifice took place in one stage in the volcanic history.Typical phenocryst assemblage is plagioclase > OPX > biotite >> Fe-Ti oxides and quartz.Micro-inclusions and glomerocrysts contain trace olivine.Phenocrysts of plagioclase and OPX show multiple textures consistent with thermal disequilibrium.Lava flows and domes were erupted with a limited compositional range (61 wt% -67 wt% SiO 2 ) of magma for ~620,000 years.Magmatic inclusions and xenoliths observed in nearly all lava flows and domes record the magmatic system beneath Uturuncu.Repeated injections of contaminated parental magma produce the compositional diversity observed in these magmatic inclusions.Xenocrysts/antecrysts of olivine and quartz in both the magmatic inclusion and the lava flows and domes, and the presence of quartz-rich xenoliths suggests crustal contamination is an important process in magma genesis.
Lavas and domes at Uturuncu show abundant evidence for thermal and chemical disequilibrium through the presence of sieved and zoned phenocryst phases and magmatic inclusions.Magma mixing controlled the limited compositional range observed in Uturuncu volcanic rocks.It is possible to explain the geochemical trends of the andesite and dacite volcanic rocks by differentiation in old, Sr radiogenic, Nd non-radiogenic felsic basement rock.The variation in isotopic composition suggests large amounts of assimilation fractional crystallization during differentiation.

Figure 2 .
Figure 2. Simplified geologic map of Cerro Uturuncu based on field mapping and satellite imagery interpretation.Modified from Sparks et al. [16].
with younger age.

Figure 3 .
Figure 3. Representative views of Cerro Uturuncu geology.(A) View southeast towards the edifice of Uturuncu; (B) Typical flow folding in lava flows; (C) Hydrothermally altered remains of the edifice of Uturuncu produced by fumaroles; (D) Typical prismatically jointed block from exterior wall of a collapsed dome; (E) Andesite inclusions in an Uturuncu lava flow; (F) Typical flow front of an Uturuncu lava flow overlying APVC ignimbrite; (G) Oxidized autobrecciation found at the terminus of lava flows; (H) Pressure ridge found commonly found in lava flows over 20 m thick.

Figure 4 .
Figure 4. Modal percent phenocrysts versus SiO 2 for representative Uturuncu lava and domes.sizefrom 0.25 mm to 1.5 mm along the long axis.Large plagioclase phenocrysts are interpreted as xenocrysts or antecrysts because they contain cores riddled with abundant, irregular shaped glass inclusions and euhedral rims and oscillatory zoning patterns similar to those observed in the host lava flows and domes.The groundmasses of these inclusions are hyalopilitic and are composed of microlite-sized (<0.25 mm) plagioclase crystals and Fe-Ti oxides in a glassy matrix.These inclusions are the most common type observed, and range in size from microinclusions (<3 mm) to ~5 cm in diameter although occasional larger inclusions were observed.The second population of inclusions are non-vesiculated basaltic andesite clots, containing plagioclase > OPX > Fe-Ti oxides > olivine > amphibole ± clinopyroxene.Hiatal to intergranular crystal sizes dominate this population of inclusions.The second population of inclusions lack glassy rims in contact with the host lava flow or dome as described above.The contacts of the inclusions are angular and defined by crystal boundaries.Crystal sizes range from microlite-size to crystals ~1.5 mm along the long axis.Inclusion size varies from micro-inclusions (~3 mm) to ~10 cm in diameter.Hand specimen size samples of this inclusion population was only observed in two lava flows, although these inclusions may have been overlooked due to their similarity in color to the host lavas when weathered.Two additional lava flows and one dome contain micro-inclusions of this population.Non-vesiculated micro-inclusions are distinct from the glomerocrysts observed lava flows and domes due to their crystal size and presence of rare amphibole.OPX phenocrysts are similar in size to the plagioclase though

Figure 5 .
Figure 5. Back-scattered electron micrographs of Uturuncu volcanic rocks.(A) Plagioclase phenocryst exhibiting resoption surfaces and growth zones; (B) Glomerocryst composed of zoned OPX phenocrysts around a single sieved and zoned plagioclase phenocryst; (C) Dacite lava flow with a sieved, unzoned plagioclase phenocryst; and (D) a glomerocryst composed of unzoned OPX phenocrysts around zoned and unzoned plagioclase phenocrysts.chemical zoning is not present.Some plagioclase phenocrysts show sieving and dissolution surfaces and all crystals over 0.5 mm contain growths zones.A reaction rim of OPX typically surrounds olivine phenocrysts.The groundmasses of the non-vesiculated inclusions are intergranular to hyalopilitic and are composed of microlite-sized acicular lath of plagioclase and blocky OPX.Vesiculated and non-vesiculated magmatic inclusions (population 1 and 2) are common features observed in CVZ volcanic rocks as well as at other intermediate composition volcanic systems around the world[21] [36]-[38].Magmatic inclusions, with or without vesiculation, represent blobs of mafic magma quenched in cooler, more silicic magma[36] due to the processes of magma mixing and mingling.This is our interpretation of inclusion populations 1 and 2 observed in Uturuncu lavas and domes.The third population of inclusions observed does not exhibit the same mineralogy as the host lavas and domes.This population is dominantly hypidiomorphic to allotriomorphic granular with trace to rare glass (0% -0.5%).Plagioclase and quartz dominate the mode of these inclusions with trace OPX and biotite.Plagioclase crystals are clear with no growth zones or dissolution surfaces observed.Quartz crystals show resorption along the corners of the crystals.Margins of these inclusions are angular and defined by the crystal boundaries.Textures observed in this population of inclusions suggest the same process does not produce these inclusions as populations 1 and 2. It is likely that inclusions in population 3 are xenoliths.Sparks et al.[16] described two xenolith populations have a similar petrographic description to inclusions in population 3 of this study.These authors interpreted observed quartz-rich xenoliths to represent crustal partial melting of rocks trapped in a MASH zone[35].It is our interpretation that xenoliths observed in this study represent basement rocks interacting with Uturuncu magmas during migration and differentiation.These xenoliths may represent the residue of a granodiorite crustal contaminant beneath Uturuncu.A fourth inclusion population described by Sparks et al.[16] and was not observed in this study.These authors describe a noritic composition xenolith population found in multiple lava flows sampled in the study.These xenoliths contain mineral assemblages and similar mineral compositions to those observed in the lava flows and domes[16].The adcumulate textures observed in these xenoliths suggested to these authors that these (A) & Figure 9(B)) suggesting extensive crustal contamination of the most mafic magmas observed in this study.

Figure
Figure Trace-element concentrations versus SiO 2 for Uturuncu lavas and domes, inclusions and xenoliths.

Figure 9 .
Figure 9. (A) Nd isotope and Sr isotope ratios for Uturuncu lavas and domes; (B)143 Nd/ 144 Nd ratios versus 87 Sr/ 86 Sr ratios for Uturuncu volcanic rocks.The lavas and domes have a large range of 87 Sr/ 86 Sr ratios for a narrow range of 143 Nd/ 144 Nd ratios (Figure9(B)).Assimilation of crustal rocks with high 87 Sr/ 86 Sr ratios but similar 143 Nd/ 144 Nd ratios explains the observed trend.Figure10illustrates six scenarios using different composition contaminants and different partition coefficients dependent on the crustal composition.The lower crustal contaminant is a garnet-sillimanite-gneiss xenolith described by McLeod et al.[40] from the southern CVZ.We selected this composition because of petrographic similarities to a sillimanite-gneiss xenolith described by Sparks et al.[16] from Uturuncu lavas.The upper crustal contaminant used in the models is Paleozoic granite from Sierra de Moreno described by Lucassen et al.[41].These granites are petrographically similar to the xenolith observed in this study.We contaminate two primitive basalt compositions with lower and upper crustal contaminants: one from Davidson and de Silva[13] for basalts from the central Altiplano, and a second from Hernando et al.[42] for pre-caldera basalts from the Payún Matrú volcanic field in western Argentina's Southern Volcanic Zone.Table3presents concentrations and isotopic ratios of the contaminants and basalts.Model curves labeled AFC (D = 1.7) and AFC (D = 1.5) were constructed to simulate the effects of differentiation on two potential parental magmas compositions contaminated by a highly evolved crustal contaminant (Figure10).Bulk distribution coefficients of 1.5 and 1.7 are high for andesites, but the decrease of Sr content with increasing SiO 2 content and in increase of Rb/Sr ratios with increasing 87 Sr/ 86 Sr ratios suggest Sr was compatible during differentiation.The curve AFC (D = 0.1) was constructed to simulate the effects of differentiation in the lower crust compared to curves for the upper crust.A bulk distribution coefficient of 0.1 is considered appropriate for this composition based on regional studies that have shown an increase in Sr/Y ratio with a decrease in 87 Sr/ 86 Sr ratios reflects incompatibility of Sr in garnet-stable/plagioclase-unstable-hybridized crust[29] [44][45].87Sr/ 86 Sr and 143 Nd/ 144 Nd ratios for Uturuncu lavas and domes create an array consistent with assimilation of old radiogenic basement rocks with relatively non-radiogenic Nd and radiogenic Sr. High Rb/Sr ratios suggest these basement rocks were also felsic in composition as opposed to noritic as could be interpreted from the presence of noritic composition xenoliths described by Sparks et al.[16].High silica, quartz-feldspathic xeno-

Figure 10
illustrates six scenarios using different composition contaminants and different partition coefficients dependent on the crustal composition.The lower crustal contaminant is a garnet-sillimanite-gneiss xenolith described by McLeod et al.

Figure 10 .
Figure 10.(A) Nd Isotope and Sr isotope constraints on bulk mixing and assimilation fractional crystallization (AFC) [43] models for Uturuncu lava flows and domes.(B) 87 Sr/ 86 Sr ratios versus Rb/Sr ratios for Uturuncu lava flows, domes and inclusions with AFC models of primary melts discussed in text and in Table3.

Table 1 and
Table 2 present, and Figure

Table 1 .
Representative whole rock major-, trace-element concentrations and isotopic ratios analyses of Uturuncu lava flows and domes.

Table 2 .
Representative whole rock major-, trace-element concentrations and isotopic ratios analyses of Uturuncu magmatic inclusions and xenoliths.

Table 3 .
Trace element concentrations and isotopic ratios of basement rocks and primary basalts used in geochemical modeling.
[16]s, observed in this study, and similar to xenoliths interpreted by Sparks et al.[16]to be remnants of a crustal hot zone could represent this contaminant.