The Central Andes is a key global location to study the enigmatic relation between volcanism and plutonism because it has been the site of large ignimbrite-forming eruptions during the past several million years and currently hosts the world’s largest zone of silicic partial melt in the form of the Altiplano-Puna Magma (or Mush) Body (APMB) and the Southern Puna Magma Body (SPMB). In this themed issue, results from the recently completed PLUTONS project are synthesized. This project focused an interdisciplinary study on two regions of large-scale surface uplift that have been found to represent ongoing movement of magmatic fluids in the middle to upper crust. The locations are Uturuncu in Bolivia near the center of the APMB and Lazufre on the Chile-Argentina border, on the edge of the SPMB. These studies use a suite of geological, geochemical, geophysical (seismology, gravity, surface deformation, and electromagnetic methods), petrological, and geomorphological techniques with numerical modeling to infer the subsurface distribution, quantity, and movements of magmatic fluids, as well as the past history of eruptions. Both Uturuncu and Lazufre show separate geophysical anomalies in the upper, middle, and lower crust (e.g., low seismic velocity, low resistivity, etc.) indicating multiple distinct reservoirs of magma and/or hydrothermal fluids with different physical properties. The characteristics of the geophysical anomalies differ somewhat depending on the technique used—reflecting the different sensitivity of each method to subsurface melt (or fluid) of different compositions, connectivity, and volatile content and highlight the need for integrated, multidisciplinary studies. While the PLUTONS project has led to significant progress, many unresolved issues remain and new questions have been raised.

Crustal magmatism plays a fundamental role in the generation of continental crust and the thermal and rheological evolution of mountain belts and plateaus (e.g., Hamilton and Myers, 1967; Fyfe, 1973; Brown, 2001). Magmas resulting from the chemical differentiation of mantle-derived magmas coupled with crustal partial melting are manifest as granite batholiths and large-volume volcanic eruptions (ignimbrites and their associated calderas). While there is compelling evidence that some large silicic ignimbrites are related to pluton-scale granite magmatic systems (e.g., Lipman, 1984; Chappell et al., 1987), there is controversy over whether the intrusive and extrusive processes occur simultaneously or indeed are related at all (e.g., Bachmann et al., 2007; Miller et al., 2011; Lundstrom and Glazner, 2016). For example, if large ignimbrite eruptions empty their magma chambers, an equivalent intrusion might not be formed (e.g., Glazner et al., 2004). Also, large magma chambers might grow without eruption (e.g., Jellinek and DePaolo, 2003; de Silva and Gregg, 2014; Cao et al., 2016). Intrusions may grow too slowly to form chambers of eruptible magma (Annen, 2009).

There are few examples of both plutonic and volcanic systems of similar age in the same area, although there are some, such as the Southern Rocky Mountain volcanic field (e.g., Lipman and Bachmann, 2015), the Kohistan arc, Pakistan (e.g., Jagoutz and Schmidt, 2012), and the Yerington batholith, Nevada (Dilles and Wright, 1988). Granite plutons are the integrated result of intrusion over hundreds of thousands to millions of years (e.g., Hamilton and Myers, 1967; Glazner et al., 2004; Paterson et al., 2011) and are typically studied as fossil systems exposed by erosion and tectonic exhumation. Volcanic products provide information on a large magma body at a single moment in time but typically have been eroded away when the intrusions are exposed.

The international PLUTONS project investigated the relationship between presumed ongoing intrusion within the middle to upper crust at two little-studied areas in the Central Andes and past, present, and future volcanism through multidisciplinary observations and numerical modeling. The two focus sites (Uturuncu and Lazufre, Fig. 1A) exhibit persistent uplift at a rate of 1–2 cm/yr and are unusual (e.g., compared to the compilation of calderas in Ruch et al., 2008 and 220 globally deforming volcanoes by Biggs and Pritchard, 2017) in that they are >40 km in diameter, and deformation has persisted for ten or more years. Enigmatically, the uplift centers are not directly below volcanoes with Holocene eruptions. Yet, they are near geophysically imaged zones of partial melt (Fig. 1B; Delph et al., 2017; Ward et al., 2017) and regions of Pleistocene volcanism.

Uturuncu is a Pleistocene dacitic volcano located behind the modern volcanic arc that last erupted 250 ka (Muir et al., 2015). It lies above a large region of inferred crustal melt called the APMB (Figs. 1 and 2), which is near an ignimbrite flare-up over the past ten million years that has formed the Altiplano-Puna Volcanic Complex (APVC, Fig. 1) and which is now waning (Table 1). Deformation at Uturuncu is spatially larger than at Lazufre and has a longer duration deformation history (Table 1). Before the start of the PLUTONS project, Uturuncu was the only known Pleistocene volcano in the world that had active fumaroles, was uplifting, and was seismically active (Sparks et al., 2008).

Lazufre (Pritchard and Simons, 2002) is ∼10 km from two Pleistocene–Holocene andesitic-dacitic arc volcanoes (Lastarria and Cordón del Azufre). Hundreds of older (Quaternary) volcanic vents are present within the deforming area (e.g., Froger et al., 2007) with circumferential and radial patterns (Ruch and Walter, 2010; Perkins et al., 2016a). A connection between Lazufre and Lastarria is suspected, because a smaller and shallower uplift began at Lastarria sometime between 2000 and 2003 (e.g., Froger et al., 2007; Ruch et al., 2009), plausibly related to its active hydrothermal system that has produced liquid sulfur flows (Naranjo, 1985, 1988). Deformation at Lazufre induced extension at Lastarria (e.g., Ruch et al., 2009; Nikkhoo et al., 2016).

The two study sites provide different perspectives on magmatic systems from intrusion to eruption within the Central Andes. One site is within the modern frontal volcanic arc (Lazufre), and the other is 70 km east of the modern arc in an area currently underlain by mid-crustal partial melt (Uturuncu). However, the two sites did not receive equal study. Uturuncu was the primary site with more papers in the themed issue than Lazufre. The choice of two sites mitigated against the possibility that one of the sites could cease activity during the project. We found several indications of decadal temporal variability at the two sites (Henderson and Pritchard, 2017; Henderson et al., 2017). We compare the characteristics and possible interpretations of the two sites in the Discussion section and Table 1.

The PLUTONS project began in 2009 with funding from the UK National Environmental Research Council (NERC) for four years (focused solely on Uturuncu) and the U.S. National Science Foundation (NSF) Continental Dynamics program for five years. The name PLUTONS is a shortened acronym from the names of the magmatic systems and the funding and partner agencies: Probing Lazufre and Uturuncu TOgether: Nsf (USA), Nerc (UK), Nserc (Canada), Sergeotecmin (Bolivia), Sernageomin (Chile), Sernap (Bolivia), observatorio San Calixto (Bolivia), Chilean Seismological service, Universidad de Salta (Argentina).

This themed volume summarizes the key results of the PLUTONS project with the following objectives:

  • (1) quantify the rates of intrusion and eruption at Lazufre, Uturuncu, and nearby centers, and determine whether those rates have been steady or episodic over the past several million years;

  • (2) identify and characterize active intrusions, including their depth, volume, temperature, melt content, and shape and how they relate to magma chambers; and

  • (3) constrain how the active intrusions interact with their surroundings through deformation, faulting, associated fluid flow, metamorphic reactions, and geothermal activity.

The PLUTONS project builds on previous regional geochemical, geological, and geophysical studies in the Central Andes. Several previous publications summarize the magmatic evolution of this region in the past ∼25 m.y. (e.g., Hawkesworth et al., 1982; de Silva, 1989; Coira et al., 1993; Ort et al., 1996; Lindsay et al., 2001; Schmitt et al., 2001; de Silva et al., 2006; Kay and Coira, 2009; Folkes et al., 2013; Burns et al., 2015; Freymuth et al., 2015) and the geophysical characteristics of the subsurface (e.g., Schmitz et al., 1997; Chmielowski et al., 1999; Schilling et al., 2006; Beck et al., 2015; Pritchard and Gregg, 2016). We note a few key points here for context. In the Central Andes, volcanic activity has taken place from the volcanic arc to hundreds of km east in the backarc, including large-volume ignimbrites (Fig. 1C), as well as other centers of mostly silicic effusive and explosive magmatism with subsidiary mafic centers (e.g., de Silva and Francis, 1991; Kay and Coira, 2009). The timing, location, and type of activity have varied along strike, but two important regions have been identified (e.g., Kay and Coira, 2009; Beck et al., 2015):

  • 1. Between ∼21° and 24° S, an ignimbrite province (>15,000 km3 erupted) between 11 and 1 Ma called the APVC (de Silva, 1989) surrounds Uturuncu (Fig. 1C). The onset of volcanism in the APVC has been correlated with marked steepening of the subducting Nazca plate from nearly flat-slab subduction prior to 16 Ma to subduction with a ∼30° dip today (e.g., Barazangi and Isacks, 1976; Allmendinger et al., 1997). The onset of steepening is thought to have induced decompression and dehydration melting in the mantle and may have facilitated delamination of the base of the continental lithosphere (e.g., Allmendinger et al., 1997; Kay and Coira, 2009). Below the APVC is a large volume (∼200 km in diameter) of inferred partial melt, called the Altiplano-Puna Magma (or Mush) Body (APMB), which has been identified by seismic (e.g., Chmielowski et al., 1999; Zandt et al., 2003), gravity (Prezzi et al., 2009), and electromagnetic methods (e.g., Brasse et al., 2002; Schilling et al., 2006).

  • 2. A region between 24° and 28° S that includes the large-volume backarc Cerro Galán ignimbrites (6–2 Ma), numerous, smaller centers, and the volcanic arc that includes Lazufre (Fig. 1C). The zone of partial melt below Cerro Galán, extending to Lazufre and other areas has been called the Southern Puna Magma Body (SPMB, Fig. 1A) (e.g., Bianchi et al., 2013. Several smaller features include the Lazufre Magma Body (LMB) and Cerro Galán Magma Body (CGMB) (Delph et al., 2017; Ward et al., 2017), but the seismic velocity reduction in the mid-crust (Fig. 1B) and the inferred amount of partial melt are less than in the APMB.

One proposed explanation (e.g., Beck et al., 2015) for the older age of ignimbrite volcanism and larger inferred amount of crustal melt in the APMB compared to the SPMB is that slab-steepening and subsequent lithospheric delamination that are driving volcanism are older in the APMB than in the SPMB. Additional factors are N-S variations in crustal thickening, tectonic variations in the upper plate, and oblique subduction of the Juan de Fuca Ridge (e.g., Kay and Coira, 2009; Freymuth et al., 2015; Brandmeier and Wörner, 2016; Naranjo et al., 2018a, 2018b).

Figure 3 shows a map of the geophysical networks deployed by the PLUTONS project, and in the Appendix, we summarize new geochemical and geophysical data and how they were analyzed. Results of these studies and numerical modeling are described in the subsequent sections.

Below we summarize key results of the PLUTONS project and other recent projects in the region in the context of previous work. We describe the results from individual studies of Uturuncu and Lazufre and synthesize them into estimates of the amount of partial melt as a function of depth. Because the regional altitudes near Uturuncu and Lazufre are 5 and 4.3 km above sea level, respectively, we label all depths as either beneath these locally defined surfaces or below sea level (bsl).

What Did We Learn About Uturuncu and the APMB?

Seismic Velocities as a Function of Depth

The thickness of the APMB seismic anomaly using newly available data and analyses through ambient noise tomography and receiver functions is on average 9–11 km (Ward et al., 2014b, 2017; McFarlin et al., 2018), as opposed to ∼1–2 km estimated in previous work (e.g., Leidig and Zandt, 2003). The dimensions of the APMB can be defined by the 2.9 km/s Vs velocity contour, because velocities below this value likely require fluids (such as melt or brines) and are difficult to explain by anisotropy or temperature variations alone (Ward et al., 2014b). Below Uturuncu, the Vs = 2.9 km/s contour is as shallow as 4 km bsl with a minimum velocity and a presumed highest melt fraction at 10–15 km bsl (Fig. 2C). PLUTONS seismic data show lateral and vertical variations in the low-velocity zone not imaged in previous work (e.g., Yuan et al., 2000), with lower Vs (<2.0 km/s) as well as varying thickness (5–20 km) and depth (9–19 km) below the surface (Ward et al., 2014b, 2017; McFarlin et al., 2018).

A different view of seismic velocities is provided by earthquake tomography. Kukarina et al. (2017) image a vertical structure 20–30 km wide below Uturuncu; this structure extends from ∼6 km below the surface into the lower crust with low Vs (Fig. 4E) and Vp/Vs >1.9 (Fig. 5A). The bottom is not well resolved because of vertical smearing of the tomography. The APMB is largely invisible to Vs earthquake tomography because the local seismicity lies mostly above it (Jay et al., 2012; Kukarina et al., 2017), while the subducting slab seismicity cannot resolve uniform layers of low velocity that extend across the entire network (Kukarina et al., 2017). The ambient noise tomography and receiver functions analysis of Ward et al. (2014b) could not resolve this vertical feature very well because it is smaller than the spatial resolution of the long-wavelength surface waves used in that study (30–100 km). However, a vertical anomaly below Uturuncu is just visible at depths greater than 20 km bsl (Fig. 4D) in the analysis of Ward et al. (2017). The Vp earthquake tomography image (Fig. 4F) is quite different from the Vs image (Fig. 4E); there is a regional low Vp at the depth of the APMB, but there is a Vp high directly below Uturuncu at 10 km depth.

Receiver functions show topography on the top of the APMB and a velocity change in the lower crust that varies laterally. Specifically, there is a change in the depth of a sharp velocity contrast from 60 km regionally (possibly the Moho) to 40 km below Uturuncu (Fig. 5) (McFarlin et al., 2018). The 40 km reflector is not the Moho since the seismic velocities below it are not mantle velocities (Ward et al., 2017), but the reflector could indicate the top of a localized mafic lower crust, possibly formed as a plexus of mafic and/or ultramafic sills or cumulates due to basaltic magma intrusion (e.g., Müntener and Ulmer, 2006; Nandedkar et al., 2014). A high Vp anomaly in the lower crust below Mount St. Helens has been explained by ultramafic cumulates (Kiser et al., 2016) and is similar in spatial scale to the observations at Uturuncu. A Vp anomaly is not detected in the lower crust of Figure 4F but would be hard to resolve (Kukarina et al., 2017).

Subsurface Electrical Resistivity

Previous work revealed a large low-resistivity subsurface zone around Uturuncu (e.g., Schilling et al., 2006), but the improved data and methods used by Comeau et al. (2015, 2016) reveal several resistive anomalies in the upper crust that may be connected with the APMB and surface volcanic features (Figs. 3B, 3C, 4B, 6B). The largest region with the lowest resistivity (C2 in Figs. 3B, 3C, and 4) is inferred to be the top of the APMB and is 5 km shallower in the 3D inversion (13–14 km bsl; Figs. 3C and 6B) than the 2D inversion (18–19 km bsl; Fig. 4B). The depth from the 3D inversion (∼13–14 km bsl) is identified by the inflection point in 1D resistivity profiles (Fig. 2B) based on synthetic resolution studies (Comeau et al., 2016). Thus the top of the APMB differs by almost 10 km from that of the joint ambient noise tomography and receiver function solution (4–5 km bsl; Ward et al., 2014b). There is a similar difference between the depth to the minimum Vs and lowest resistivity (Fig. 2) (Comeau et al., 2015). This difference is attributed to the different sensitivities of these methods to partial melt of different compositions (discussed in the next section and Fig. 2; Table 2; Comeau et al., 2016; Pritchard and Gregg, 2016).

Gravity and Subsurface Density

A low isostatic and Bouguer gravity anomaly over the APMB (e.g., Götze et al., 1994; Prezzi et al., 2009) and a dense grid of observations used to generate new models have produced a range of plausible subsurface density structures (del Potro et al., 2013). These models show low-density structures connecting the APMB to the surface (dotted lines in Fig. 4B). These structures overlap the MT anomalies.

“Fossil Dacite Magma Chamber” and Hydrothermal System

The volcanic evolution of Uturuncu was summarized by Sparks et al. (2008), Michelfelder et al. (2014), and Muir et al. (2015), with additional geochemical, geochronological, and petrologic data from Michelfelder et al. (2013) and Muir et al. (2014a, 2014b). Volcanic rocks from the edifice of Uturuncu are dominantly crystal-rich, high-K, calc-alkaline, dacite lava flows and domes with phenocrysts of plagioclase > orthopyroxene > biotite > quartz. Groundmass minerals and glass make up 70%–50% of the modal volume. Eruptive activity at Uturuncu was episodic with six eruptive periods separated by >50 k.y. hiatuses with major edifice construction occurring between 595 and 505 ka (Muir et al., 2015).

Petrological analysis (experiments and melt inclusions) suggests that the dacite lavas that erupted at Uturuncu >250 ka were stored 2–4.5 km below the surface (at or above sea level, Muir et al., 2014a, 2014b). There are geophysical anomalies (Table 2: low resistivity, earthquakes, and low Vs; Figs. 3C and 4; e.g., Jay et al., 2012) at approximately the same depth region today. While a shallow dacite magma chamber fed Uturuncu eruptions, it is unlikely that magma is present at shallow depth today. In order to explain the observed resistivities, 80%–100% dacite partial melt would be required (Fig. 6C) because dacitic melt is resistive (e.g., Laumonier et al., 2015; Comeau et al., 2015, 2016). This amount of melt at these depths is implausible—being ruled out by the nonzero values of Vs (e.g., Jay et al., 2012; Kukarina et al., 2017). Instead, the low resistivity of C4 (Fig. 6B) can be explained by saline fluids of magmatic and/or meteoric origin (Comeau et al., 2015, 2016). The resistivity of brines at 300–500 °C (Fig. 2) is 0.1–0.01 Ωm (e.g., Unsworth and Rondenay, 2013), which is consistent with the observed bulk resistivity for plausible values of porosity and fluid interconnectivities (Fig. 6D). For example, if we assume salinity to be 30,000 ppm (3%) and temperature 200–400 °C, the brine resistivity is calculated to be ∼0.03 Ωm (Ussher et al., 2000). For plausible upper-crustal interconnected porosity of 5%, we calculate the bulk resistivity observed in the MT-derived model (5 Ωm) using a modification of Archie’s Law (Fig. 6D).

The proposed NaCl content is realistic given the composition of Uturuncu dacites, phase relations in the system NaCl-H2O, and the NaCl contents (10%–50%) measured at other volcanoes (e.g., Aguilera et al., 2016). The dacite magma prior to eruption, as shown by melt inclusions (Muir et al., 2015), contained ≤3.9 wt% H2O and 0.23 wt% chlorine. If we convert this composition into a bulk equivalent salinity (NaCleq, i.e., assuming all Cl is complexed by cations such as K, Na, and Ca), we calculate 9 wt% or 2.7 mol%. At a solidus temperature of ∼600 °C, such a fluid will be supercritical at pressure above 90 MPa but condense to a mixture of high-salinity brine and low-salinity vapor at lower pressures (Driesner and Heinrich, 2007). The corresponding depth depends on the density of the surrounding rock and the extent to which the fluids are interconnected, i.e., lithostatic or hydrostatic pressure. In the former case, with a rock density of 2600 kg/m3 (Fig. 2A), 90 MPa corresponds to 3.5 km depth. In the latter case, taking fluid density as 1300 kg/m3, the depth is 7 km. Thus as magmatic fluids ascend and separate in the depth range 3.5–7 km, brines with compositions consistent with the low-resistivity anomaly are formed (Fig. 6B). Numerical modeling of fluid release and brine condensation above solidifying plutons (Weis, 2015) indicates that the brine can accumulate in cylindrical bodies similar in scale to those imaged. For a pluton at 5 km depth exsolving fluids with 10 wt% NaCleq, Weis’ (2015) calculations show brine accumulation at 1–5 km depth as shown in Figure 6D.

Inferred Partial Melt

Converting geophysical observables to estimates of the amount of partial melt requires assumptions about the nature of crustal anisotropy (e.g., Ward et al., 2014b), the distribution of melt at the grain scale (i.e., dihedral angle; see Forsyth, 1993) and in veins, dikes, and sills too small to be resolved by geophysical inversions, and the level of melt interconnectivity (e.g., Unsworth and Rondenay, 2013). For example, the calculations shown in Figure 6 assume that melt is interconnected (or not) at the microscale along grain boundaries, but larger-scale melt-rich regions (composed of layers, dikes, sills, veins, etc. spanning meters to km) are considered plausible (e.g., Schmitz et al., 1997; Jackson et al., 2003; Schilling et al., 2006; Solano et al., 2012, 2014; Sparks and Cashman, 2017). They cannot be geophysically resolved as individual features but might be seen in seismic anisotropy if there was a dominant orientation (e.g., Jaxybulatov et al., 2014).

The melt percentage of the APMB is at least 10%–25% based on combining the seismic and MT observations with laboratory-derived relations between Vs, resistivity, attenuation, and melt percentage using appropriate pressure, temperature, and compositional inferences from volcanic samples (e.g., Fig. 6; Schilling et al., 2006; Ward et al., 2014b; Comeau et al., 2015, 2016; Farrell et al., 2017). Melt volumes up to 80%–90% are permitted by the gravity and MT data (e.g., del Potro et al., 2013; Comeau et al., 2015). But bulk values of partial melt above 35%–50% would produce a Vs “shadow zone” in silicic melts (e.g., Chu et al., 2010)—i.e., no S waves would be recorded at the surface above this zone, which has not been observed (e.g., Schilling et al., 2006).

Temporal Variations in Ground Deformation

A globally anomalous pattern of surface deformation is centered approximately on Uturuncu with a central uplift and surrounding ring of subsidence with a diameter of 150 km that some have called a “sombrero” (e.g., Fialko and Pearse, 2012; Henderson and Pritchard, 2013). The sombrero is superimposed on an ∼1-km-high topographic anomaly across the APVC with a wavelength of ∼230 km (Perkins et al., 2016b). The deformation seems to have been ongoing from at least 1965–2010 with a maximum uplift rate of 1.2 ± 0.1 cm/yr and subsidence of 0.3 ± 0.03 cm/yr coincident with the interferometric synthetic aperture radar (InSAR) observed uplift and subsidence (Gottsmann et al., 2017). But, the rate is variable—in the 1990s, vertical rates are greater than 0.7 ± 2 cm/yr at GPS station UTUR; rates decrease to 0–0.1 cm/yr in 2010 and then increase to 0.9 cm/yr in 2015 (Henderson and Pritchard, 2017). There is no net geomorphic deformation of river profiles or lake shorelines across the current uplift field at Uturuncu, which implies that the uplift must have been of less than 100 years duration at the average uplift rates from 1992 to 2011 (Perkins et al., 2016a). Any similar episodes of uplift must have been balanced by subsidence to create no net deformation over the past 16,000 years (Perkins et al., 2016a). A model for the distribution of lineaments around Uturuncu is presented by Walter and Motagh (2014) as evidence of subsidence, although the structural interpretation of these features has been questioned by the original authors as well as by Perkins et al. (2016a). Given the wide spatial footprint of the deformation and coupled uplift and/or subsidence, the deformation is likely caused by processes in the middle to lower crust (regional crustal thickness is 60–70 km). Possible processes include an active diapir (e.g., Fialko and Pearse, 2012; del Potro et al., 2013) and the transfer of magma and magmatic fluids from deeper to shallower levels within a large-scale mushy magma system (Fialko and Pearse, 2012; Henderson and Pritchard, 2013, 2017; Gottsmann et al., 2017; Sparks and Cashman, 2017). A mid-crustal source is consistent with the lack of gravity change associated with the deformation (Gottsmann et al., 2017). We discuss the merits and limits of these models in the section on “Conceptual Models for Deformation at Uturuncu and Lazufre.”

Synthesis of Uturuncu Observations into a Depth Profile

Here we explore the possibility of making a self-consistent interpretation of all data sets, considering that each one is sensitive to different properties, following the approach of Comeau et al. (2016). Our approach is non-unique and can only be considered semiquantitative because we have not created a unified numerical model that is consistent with all data sets. We put forward this model (Fig. 2; Table 2) with some falsifiable predictions to summarize our current understanding. The key interpretations are at three depths: (1) shallow depths less than 5 km below the surface; (2) between 9 and 18 km depths; and (3) 18 km and below.

  • (1) At depths 2–4.5 km below surface, the low resistivity (Fig. 4C) and low Vp and Vs (Fig. 4C; Jay et al., 2012; Kukarina et al., 2017) suggest brines in the hydrothermal system (Fig. 5D). As depth increases, the permeability decreases, and so the geophysical expression of these brines decreases.

  • (2) At ∼9 km below the surface, the top of the APMB is defined by Vs <2.9 km/s (presumed to be at the dacite solidus). The partially molten dacite reduces the Vs, but it is relatively dry (3.2 ± 0.7 wt% H2O; Muir et al., 2014b) with high viscosity. Consequently, the melt phase may not be well connected, and the melt thus does not reduce the resistivity (Fig. 6C). The high Vp at this depth (Fig. 4F; Kukarina et al., 2017) may reflect a mafic contribution.

  • (3) There is an inflection point in the resistivity depth profile (Fig. 2B; Comeau et al., 2016) ∼18–19 km below the surface where the resistivity decreases significantly. This inflection point is interpreted to be the top of the partially molten andesite. This inferred andesitic melt extends possibly to 25–39 km below the surface (Table 2). The low resistivity is explained by the high water content of the andesite (∼9 wt% H2O; Laumonier et al., 2017) and the melt phase being well connected (Fig. 6E).

Our conceptual model implies approximate temperatures at the bottom of the seismogenic zone, at the top of the APMB as defined seismically by Ward et al. (2014b), and at the top of the APMB as defined using magnetotellurics by Comeau et al. (2016). To estimate the present-day surface heat flow, we assume that heat loss is purely by thermal conduction, that the system is in steady state, use the inferred temperatures in Figure 2D and Table 2, along with the equation for heat flow = k*dT/dz (where k is thermal conductivity and dT is the temperature difference between depths dz), to find a heat flow of 150–200 mW/m2. These values can be tested with future heat-flow observations. We emphasize that the depths shown in Figure 2 and Table 2 correspond to depths below Uturuncu, and we know that these vary spatially at other locations within the APMB.

What Did We Learn About Lazufre?

Seismic and Electrical Resistivity as a Function of Depth

Prior to 2002, the Lazufre magmatic system was not known, and only limited geophysical data had been collected. Yuan et al. (2000) identified a region of potential partial melt from one seismometer near Lazufre. Later work using regional Vp earthquake tomography identified a zone of partial melt in the Puna backarc called the SPMB (e.g., Bianchi et al., 2013). Ambient noise tomography that included PLUTONS seismometers (Fig. 1B; Ward et al., 2013; Delph et al., 2017; Ward et al., 2017) showed a low Vs zone below Lazufre ∼50 km in diameter from ∼10–20 km bsl that could be an isolated body or related to the SPMB (e.g., Fig. 1B; Ward et al., 2013; Delph et al., 2017; Ward et al., 2017). Shallower low-velocity zones have been found in the Lastarria-Lazufre magmatic system using seismic networks focused on Lastarria (Fig. 7; Table 3; Spica et al., 2015) using ambient noise tomography. Two of these areas within the upper 6 km beneath Lastarria have minimum Vs 1.25–1.8 km/s (implying >30% melt, Fig. 6A) and are thought to be related to a small magmatic system feeding a hydrothermal reservoir (Fig. 7B; Spica et al., 2015). The locations of these low seismic velocity zones also correspond to low resistivity as measured by MT that could be either brines or partial melts (implying melt >40% for andesite; Fig. 7D) in the upper 7 km below Lastarria (Díaz et al., 2015; Fig. 7C) and below 20 km at Lazufre (Budach et al., 2013).

Temporal Variations in Ground Deformation

Ground uplift at Lazufre occurs on two spatial scales: (1) the broad, ∼40-km-diameter deformation from Lazufre that began in 1998, and (2) the smaller pattern at Lastarria volcano that spans a few kilometers and began ∼2003 (e.g., Pritchard and Simons, 2002; Froger et al., 2007; Ruch et al., 2009; Fig. 7A). The Lazufre deformation was ongoing in 2016, although the rate of deformation has varied in time from ∼1 cm/yr from 1998 to 2003, 3 cm/yr from 2003 to 2010, and ∼1.5 cm/yr from 2010 to 2016 (e.g., Ruch et al., 2009; Henderson and Pritchard, 2013; Pearse and Lundgren, 2013; Remy et al., 2014; Henderson et al., 2017). Whether deformation has continued at Lastarria after 2008 is not currently known. Analysis of topography suggests that this is not the first episode of uplift at Lazufre because the deformation pattern is coincident with a 500-m-high topographic swell (Froger et al., 2007), implying a source characterized by inelastic deformation. To assess the age of this swell, Perkins et al. (2016a) measured the orientation of lava flows of different ages in the area and concluded that lava flows older than ca. 0.4 Ma were not affected by the presence of the swell—thus, the swell must be younger than 0.4 Ma. The swell formation is thought to be older than ∼16,000 years, since the lake shorelines near Lastarria are not tilted (Perkins et al., 2016a). Thus the current uplift is very short-lived, or multiple uplift episodes are counterbalanced by subsidence to produce no net deformation. The implications of this are further explored in the Discussion section.

Summary of Lastarria-Lazufre Depth Profile

We have summarized all of the available observations for the Lastarria-Lazufre magmatic system in Table 3. Given currently available data sets, Vs is lower in the shallow crust beneath Lastarria than at Uturuncu; so we expect the heat flow to be higher (assuming thermal equilibrium and the composition is similar), and indeed the fumarole temperatures are higher (e.g., Table 1). If the low-resistivity and Vs region in the upper 6 km is dacite with >30% partial melt (Figs. 6A and 7D), the solidus curve in Figure 2 suggests a temperature of 660 °C—or hotter still, if the composition is andesite. This implies a heat flow of >700 mW/m2. The temperature and heat flow could be lower if the seismic and MT anomalies are due at least in part to hydrothermal fluids and not partially molten rock.

The Plutonic-Volcanic Connection at Uturuncu and Lazufre

We now return to the question of the plutonic-volcanic connection that motivated the PLUTONS project. First of all, how do we identify a pluton or batholith in the subsurface? We do not have a simple or straightforward discriminator. For example, solidified batholiths can exhibit a high resistivity if they are composed of unaltered granite with low permeability (e.g., Comeau et al., 2015) or a low resistivity if hydrothermal circulation produced alteration and clay minerals. Low density in the subsurface can indicate solidified plutons, but the interpretation is non-unique (e.g., del Potro et al., 2013; Lipman and Bachmann, 2015). For a region that is currently partially molten, the geophysical anomaly could indicate country rock that has partially melted (i.e., a migmatite) or alternatively near-solidus igneous rock that will eventually become a batholith (e.g., Ward et al., 2014b). For the APMB, both interpretations are possible at different depths as part of a mature trans-crustal magmatic system. Based on the geophysical evidence for the presence of partial melt and the history of large ignimbrite eruptions in the region, many workers have concluded there is strong evidence for a large composite upper-crustal batholith beneath the APVC, but there are varying estimates of the volume and shape of this body or bodies (e.g., Zandt, et al., 2003; de Silva and Gosnold, 2007; Ward et al., 2014b; Kern et al., 2016). The PLUTONS data resolve some of these varying estimates.

When interpreting the geophysical observations in terms of the crustal-scale plutonic structure, we consider two end members developed at different trans-crustal magmatic systems around the world: (1) discrete depths of magma accumulation surrounded by country rock, including partially molten migmatite (Fig. 8C) or (2) a continuous composite batholith that extends through the entire crustal column (Figs. 8A and 8B). At Uturuncu, an interpretation based on a single geophysical measurement might (erroneously) appear to favor one end member over the other. At one extreme, the MT zones of low resistivity (e.g., C2 and C3, Fig. 9) are clearly separated by ∼10 km, possibly suggesting separate zones of current or previous partial melt, such as in Figure 8C. At the other extreme, the regional sources of the ignimbrites (4–8 km below the surface, e.g., de Silva et al., 2006) could represent the top of the APMB. Thus, this entire region, from ∼4 km to the bottom of the APMB ∼30 km below the surface, is interpreted to be igneous rock that forms a batholith (Burns et al., 2015; Ward et al., 2014b; similar to Fig. 8A).

Of course, these end members are both simplifications—for example, the vertically continuous batholith shown in Figure 8A was inferred from gravity data, and this type of data analysis has reduced resolution and may not image vertically elongate plutons separated by country rock (Lipman and Bachmann, 2015, and references therein). We favor an interpretation with elements of both end members at Uturuncu (Figs. 2 and 5), where a large batholith has been formed from an ensemble of plutons (e.g., Fig. 8B; Ward et al., 2014b, 2017), but with lateral heterogeneity in the upper 18 km (at least) resolved by the gravity and MT (Figs. 4B and 4C), but not by seismology (although there are hints of upper-crustal variability in Figs. 4E and 4F; Jay et al., 2012). We revise the Ward et al. (2014b) interpretation of a single continuous batholith from ∼2–30 km below the surface to instead include (in the upper 18 km at least) igneous material surrounded by unmelted and partially melted country rock at scales that are seen in the field (Fig. 8B; Jagoutz and Schmidt, 2012; Fig. 4 of Paterson et al., 2011), inferred from geochemistry (e.g., Sparks et al., 2008; Michelfelder et al., 2013, 2014) and geochronology (e.g., Schmitt et al., 2003; Kern et al., 2016), but are not well resolved by the available geophysical data. Even after the ∼10 m.y. of existence and extensive magmatic processing of the APVC, the APMB has not completely obliterated all traces of country rock. The consequence of our interpretation is the intruded:extruded ratio is likely below the 35:1 extreme from Ward et al. (2014b, their figure 7); although see Tierney et al. (2016) for a 75:1 petrological estimate in the region. Because we do not yet have a good discriminant between what is country rock versus igneous rock that is part of the batholith between 2 and 30 km depth, we cannot estimate this ratio.

Focusing on just the upper 10 km, we can relate erupted products to geophysical anomalies of potential plutonic origin. At Uturuncu, there are several EW zones of high and low resistivity that all extend from ∼3 km above sea level to 3 km bsl and are ∼15 km long and ∼4.5 km wide (Fig. 9). If the anomalies signify solidified or partially solidified dacite, their total volume >1000 km3, compared to the Uturuncu erupted volume of 80 km3 (Sparks et al., 2008), provides an intrusion:extrusion ratio of >12:1, at the low end of the values proposed by Ward et al. (2014b). We interpret the pattern of high- and low-resistivity zones around Uturuncu as a propylitic hydrothermal alteration (e.g., Sillitoe, 2010) of varying permeability, temperature, and clay composition with tectonically controlled upflow of low pH magmatic fluids and regional inflow of neutral pH groundwater. A possible analogy of such an active zone of hydrothermal alteration imaged with MT is described by Gíslason et al. (2015), although the Uturuncu pattern is more complex. We do not think that the EW low-resistivity features are dikes as proposed by Comeau et al. (2016). The high-resistivity zones could be a shallow pluton related to intrusion of dacite magma, whereas the low-resistivity zones are attributed to saline fluids and hydrothermal alteration. Note that zones more distant from the summit of Uturuncu appear smoother because in those regions there are fewer MT stations.

For the Lastarria-Lazufre magmatic system, the end member of Figure 8C depicting discrete magma bodies is likely more appropriate than the continuous pluton (Figs. 8A and 8B). Geophysical inversions indicate isolated pockets of partially molten rocks at different depths that are separated presumably by country rock or solidified plutonic rock (Figs. 7B and 7C). The volume of the seismic anomaly from 3 to 6 km is 81 km3, compared to the Lastarria erupted volume of ∼10 km3, also suggesting an intruded:extruded ratio (8:1) at the low end of the range of Ward et al. (2014b). For Lazufre, Naranjo et al. (2018b) estimated a lava volume of 35 km3 that has erupted since 0.4 Ma with another 85 km3 of lava erupted between 3.5 and 0.5 Ma, compared to the partially molten volume estimate of Remy et al. (2014) >400 km3 (or 2769 km3 for the best-fitting model)—providing a range of intruded:extruded ratio from 3:1 to 23:1.

The difference between Lazufre (more similar to Fig. 8C) and Uturuncu (more similar to Figs. 8A and 8B) could reflect these two different systems being at different stages in their evolution as mentioned above—Uturuncu and the APMB are mature, and activity is waning; while the SPMB and Lazufre are less mature, and activity could be increasing in the area, such that the zones of partial melt beneath Cerro Galán, Lazufre, and the southern terminus of the active arc (Fig. 1B) might eventually coalesce into an APMB-like magma body. Alternatively, perhaps the amount of melting will never be as extreme in the SPMB as the APMB because it is tectonically and compositionally distinct such that the crust is less fertile for melting, or the heat source is smaller (e.g., Beck et al., 2015). Certainly there is no evidence at Lazufre for a major magmatic flare-up such as the APVC ignimbrite province within which Uturuncu developed, although there is a long history of explosive volcanism and tectonomagmatic interaction (Naranjo et al., 2018a).

Future Eruptions at Uturuncu and Lazufre?

Considering that neither Uturuncu nor Lazufre are unambiguously associated with a Holocene volcano, the eruptive threat from these deformation centers is unclear. Uturuncu was the first area of active ground deformation observed at a Pleistocene or older volcano, but now a handful of these are known, including Sillajhuay in the Central Andes (Pritchard et al., 2014a) and elsewhere (e.g., Africa: Biggs et al., 2013; or New Zealand: Hamling et al., 2016). These features have been called “zombie volcanoes” by a headline writer of a popular science news Web site because they are “alive” even though they appear “dead” (Pritchard et al., 2014b).

We argue that Sillajhuay and Uturuncu should be added to the list of potentially active volcanoes, both because of their current activity and the fact that volcanoes in the Central Andes can have >105 years between eruptions (e.g., Clavero et al., 2004; Klemetti and Grunder, 2008). Sparks et al. (2008) noted that given the long time interval since the last eruption at Uturuncu and the evidence for ongoing magma intrusion from the ground deformation, there was the potential for a large eruption, perhaps even a caldera-forming eruption. But there is still no geophysical evidence for the shallow accumulation of eruptible magma at Uturuncu or Sillajhuay that could lead to eruption in the near future. In contrast, there are currently potentially small volumes of magma near the surface in the Lazufre-Lastarria system (e.g., Tamburello et al., 2014; Spica et al., 2015; Lopez et al., 2018). However, these shallow accumulations may not be relevant because there is increasing evidence that prior to some caldera-forming eruptions, large volumes of silicic magma are transferred into shallow crustal chambers from much deeper in only decades to centuries (e.g., Druitt et al., 2012; Allan et al., 2013). In any case, we recommend continued vigilance at these potentially active volcanoes, to better understand their background levels of activity and characterize temporal changes.

Conceptual Models for Deformation at Uturuncu and Lazufre

How does the ground deformation at Uturuncu and Lazufre relate to the model of batholith formation and eruption potential discussed above? Previous work inferred that a magma intrusion was the most likely cause of deformation because of the deep sources at Lazufre (8–14 km below the surface; e.g., Pearse and Lundgren, 2013; Remy et al., 2014) and Uturuncu (6–19.5 km below the surface; e.g., Hickey et al., 2013, and Gottsmann et al., 2017). The hydrothermal systems are considered to be shallower (e.g., Pritchard and Simons, 2002; Sparks et al., 2008). If the volume changes are assumed to be due to magma intrusion, then the present-day rate (∼10–2 km3/yr or 1 m3/s at Uturuncu) is much higher than the long-term eruptive average (∼6 × 10–5 km3/yr at Uturuncu; Muir et al., 2015). Similarly, the rate of inferred magmatic intrusion at Lazufre (∼7 × 10–3 km3/yr; Henderson et al., 2017) is more than two orders of magnitude higher than the long-term average since 4 Ma (∼3 × 10–5 km3/yr; Naranjo et al., 2018b). The putative magma intrusion must therefore be short-lived and is consistent with a model of batholith growth by numerous increments of magma addition over millions of years (e.g., Coleman et al., 2004; Glazner et al., 2004).

However, as described above, the lack of geomorphic uplift over the past 16,000 years at Uturuncu and Lazufre causes us to consider alternatives to the model of magma intrusion as the cause of the surface deformation. To produce effectively zero net deformation over millennia, uplift from magma injection and/or melting of surrounding country rock (i.e., anatexis that involves a significant volume change; e.g., Fialko et al., 2001) would have to be balanced by equal subsidence from cooling and crystallization of the new magmatic intrusion and surrounding country rock including processes such as “reverse diapirism,” where mafic cumulates in a magma reservoir cause sinking (e.g., Roman and Jaupart, 2016). These subsidence mechanisms may be occurring at calderas that have been subsiding for decades without eruption, such as Cerro Blanco in the Central Andes (Henderson and Pritchard, 2013) and elsewhere (e.g., Parker et al., 2014). However, to get the uplift and subsidence to be nearly equal requires a narrow range of parameter space as demonstrated by viscoelastic and thermoelastic modeling including crystallization in the case of the Socorro magma body, New Mexico (Pearse and Fialko, 2010). This narrow range of parameters, while possible, may not be very likely.

An alternative conceptual model is that supercritical fluids from the partially molten zones at depths >10 km at Lazufre and Uturuncu could become temporarily trapped, causing pressurization of a subsurface reservoir and ground uplift (e.g., Gottsmann et al., 2017). Such a model has been proposed for Soufriére Hills Volcano, Montserrat (Christopher et al., 2015; Sparks and Cashman, 2017), although these papers do not explore the implications of the model for ground deformation. At Uturuncu, this pressurization would occur in a column above a bulging andesitic APMB (with a base at ∼18 km below the surface; Fig. 4; Table 2) and would occur at the same time that lateral depressurization occurred in the APMB—the result would be the distinct pattern of uplift surrounded by a moat of subsidence observed (Gottsmann et al., 2017). At Lazufre, the gases might be only partially trapped, since there is a significant flux of gas out of Lastarria (Tamburello et al., 2014; Lopez et al., 2018).

The impact of the gas pressurization and/or depressurization cycles on ground deformation are likely to be of equal magnitude due to mass conservation, creating no net geomorphic deformation. In fact, such an equal and opposite pattern of uplift and subsidence was seen in the Central Andes at Cerro Overo between 1992 and 2011 (Henderson and Pritchard, 2013) where deformation changed sign over a short timespan (years or less), which would be more consistent with gas movement than bulk thermodynamic, thermoelastic, and viscoelastic processes. Compared to Uturuncu, the Cerro Overo deformation has a smaller maximum amplitude (∼4 cm) and rate (∼0.5 cm/yr), temporal and spatial scale (<10 years for a pulse of uplift or subsidence spanning 20 km diameter), perhaps explained by its location on the periphery of the APMB (with lower melt and gas content), while Uturuncu is near the center (Fig. 1A).

Whether the current uplift pattern at Lazufre will be counterbalanced by subsidence in the future is not clear—while there is geomorphic evidence for a 500 m topographic swell associated with the Lazufre uplift and radially oriented volcanic vents formed ca. 0.4 Ma, net uplift over the past ∼16,000 years is negligible (figure 2 in Perkins et al., 2016a). Is the current Lazufre uplift an intrusion that will cause net uplift consistent with the long-lived topographic swell, or is it a transient feature that will produce no net deformation over timescales of centuries to millennia? We favor the latter interpretation for two reasons: (1) the topographic swell does not seem to have been growing during the recent millennia, and (2) although the current uplift and the swell are in the same location, the current uplift has a radius that is 10–15 km smaller than the swell (Perkins et al., 2016a), suggesting that different processes or source depths could be responsible. The difference between current uplift and the regional topographic anomaly is also true at Uturuncu—the ∼70 km regional uplift is much smaller than the ∼230-km-diameter, 1-km-amplitude regional topographic anomaly associated with the APMB (Perkins et al., 2016b). While it is possible that repeated uplift events such as those at Uturuncu and Lazufre could build the observed larger topographic features if the centroid of deformation moved around, the lack of any deflected geomorphic features in the past 16,000 years (Perkins et al., 2016a) leads us to suggest that the process or processes contributing to the creation of the 0.5–1 km topographic features have not been active over the past millennia at these two sites.

Although the millennial-scale deformation has only been studied at a few volcanic systems, there seem to be alternating episodes of uplift and subsidence that produce little net deformation in areas of current active deformation, such as Yellowstone and Socorro (e.g., Pierce et al., 2002; Finnegan and Pritchard, 2009). On the other hand, we know that such episodes of net deformation do occur with timescales of centuries to millennia at restless calderas (e.g., Campi Flegrei, Cinque et al., 1985; Laguna del Maule, Chile, Singer et al., 2015) and resurgent calderas (Kaizuka et al., 1989; Chen et al., 1995; Kennedy et al., 2012; de Silva et al., 2015). Based on recent crustal volcanism in the Central Andes, we suspect that episodes of magma intrusion in the middle crust (that could be related to batholith formation as described above) and surface uplift are still occurring in the region, but that these intrusions have not yet been observed in real time. We therefore propose that intrusions that cause long-term uplift are less common than uplift from pressurization by magmatic fluids and that gas and/or volatile processes are the cause of the five uplifting and/or subsiding areas in the Central Andes between 1992 and 2016 (e.g., Pritchard and Simons, 2004; Henderson and Pritchard, 2013; Pritchard et al., 2014a). Of course, some process is needed to cause the proposed change in volatile movements, and one such mechanism is mush destabilization as proposed by Sparks and Cashman (2017) driven by intrusion or inherent instabilities in the system. Further numerical modeling is needed to determine if these processes would still cause net geomorphic uplift or be nearly volume neutral.

Following our proposed conceptual model, the observed ground deformation at Uturuncu and Lazufre and most other deformation sites in the Central Andes (e.g., Putana, Cerro Overo, and Sillajhuay) is not fundamentally part of the process of batholith creation but is a consequence of large bodies of partially molten rock in the mid-crust releasing volatiles. These volatile phases are important for the generation of ore bodies (e.g., Hedenquist and Lowenstern, 1994; Cloos, 2001; Einaudi et al., 2003; Ingebritsen et al., 2010; Blundy et al., 2015; Weis, 2015). If correct, this conceptual model suggests that deformation could occur over any large zones of partial melt that are exsolving volatiles where the volatiles do not have a direct path to the surface and can undergo periodic trapping and release. This hypothesis would predict that if we looked over a longer time period, other areas above zones of partial melt in the Central Andes would show uplift and subsidence. We suggest that any of the topographic swells of similar scale to Lazufre seen in the long-wavelength topography of Figure 1 of Perkins et al. (2016a), including Putana, Llullaillaco, Cerro Galán, Cordón de Puntas Negras, El Queva, Antofalla, etc., could host partial melt and ground deformation.

In fact, there are several signs of activity in these areas—bursts of uplift and seismic activity have been seen at Putana (Pritchard et al., 2014a), uplift and subsidence at Cordón de Puntas Negras (at Cerro Overo; Henderson and Pritchard, 2013), as well as seismic swarms at Cerro Galán (Mulcahy et al., 2014). One potential trigger for these bursts could be shaking from subduction zone earthquakes (as proposed for the 1998 onset at Lazufre by Pritchard and Simons, 2002), but not all earthquakes impact changes in ground deformation (e.g., Henderson and Pritchard, 2017); therefore, the link is still speculative.

Regional and Global Context for Uturuncu and Lazufre

Ground Deformation

The seven large deformation patterns in the Central Andes (Uturuncu, Lazufre, Hualca Hualca, Cerro Overo, Cerro Blanco, Putana, and Sillajhuay) have rather simple ellipsoidal shapes whose centers do not move significantly over years to decades, unlike the patterns in some other areas of silicic magmatism, where the pattern, loci, and sign of deformation are more variable (e.g., Yellowstone: Wicks et al., 2006; and Cordón Caulle: Jay et al., 2014). We interpret the simpler deformation patterns in the Central Andes to reflect less spatial heterogeneity in the magma storage systems at least over the time period of observation—the ground deformation at Cordón Caulle and Yellowstone is caused by shallower magmatic sources (∼5–10 km below the surface) that provide thermal energy to shallow hydrothermal systems that could be more extensive in those areas with more ground water than the arid Central Andes. For comparison, Yellowstone Caldera is thought to have a magma reservoir with a volume of ∼46,000 km3 partial melt (Huang et al., 2015), while Ward et al. (2014b) estimated that the APMB contains a partially molten region of ∼500,000 km3. Among the locations undergoing deformation in the Central Andes and globally, the Uturuncu deformation is unique in having a moat of subsidence, which is attributed to the large partially molten zone beneath the volcano (Gottsmann et al., 2017). Although the Putana uplift is also near the APMB, it is much shallower, probably related to the hydrothermal system—which raises the question of whether the Uturuncu hydrothermal system also undergoes deformation. Deformation caused by the Uturuncu hydrothermal system has not yet been observed (although see Lau et al., 2017).

Surface Temperature and Fumaroles

Movements of magmatic fluids and gases have been proposed as the cause of ground deformation at Uturuncu and Lastarria-Lazufre, and both areas have fumarolic fields. Uturuncu has one of the largest satellite-detected hotspots and fumarolic fields in the Central Andes in terms of surface area (Fig. 10; Jay et al., 2013) and is the only hotspot at a Pleistocene volcano in the region. At Lazufre, there are no active fumaroles yet detected near the center of the Lazufre deformation (e.g., Pritchard and Simons, 2004; Jay et al., 2013), but the nearby Lastarria volcano has a well-studied fumarolic field (e.g., Naranjo and Cornejo, 1992, 2010; Aguilera et al., 2012; 2016; Lopez et al., 2018) and satellite-detected thermal anomaly (Jay et al., 2013). It has been suggested that the Lastarria fumarole field serves as a pressure valve for Lazufre (Froger et al., 2007).

The Uturuncu and Lastarria satellite detected hotspots, and in fact all volcanic hotspots in the Central Andes and their associated thermal features (fumaroles, hot springs, geysers, etc.), are restricted to areas with a spatial extent of a few km2 (Jay et al., 2013). In comparison, satellite hotspots and thermal features at Cordón Caulle and Yellowstone have areas greater than 10 km2 (e.g., Vaughan et al., 2012; Jay et al., 2013). As with the deformation patterns, the thermal characteristics of Cordón Caulle and Yellowstone suggest those systems are more active at shallower depths than Uturuncu and Lazufre, with significant lateral heterogeneity.

Samples from the Lastarria hydrothermal system show evidence for mixing of (1) saline fluids superheated, by a deep, magmatic source that was ultimately derived from mantle melting (3He/4He ratio of 4–6 times the average atmospheric composition) with (2) cooler, shallower, less saline groundwater (Aguilera et al., 2012, 2016; Lopez et al., 2018). The emission rate of SO2 was measured by mini–differential optical absorption spectroscopy (DOAS) by Tamburello et al. (2014) in 2012, and in 2014 using several techniques by Lopez et al. (2018), who measured gas composition with the MultiGAS instrument to calculate the total emission rate of all volatiles of ∼12,400–13,500 t/d. The high flux of SO2 and the CO2/SO2 ratio <2 are both consistent with shallow magmatic degassing. Temporal variations in several gas ratios (increase in H2O/CO2 and maybe HCl/St along with a decrease in CO2/St and H2O/St) between the first campaign measurements in 2006 and the most recent in 2014 suggest a drying out of the hydrothermal system and an increase in degassing of magma soluble gases with time. Lopez et al. (2018) interpret this as evidence for magma ascent between 2006 and 2012–2014 or isobaric crystallization-induced degassing of a stalled magma body.

In contrast, there are no published reports on the gas discharge from Uturuncu. One argument against the role of gas in causing deformation is that the gas release at Uturuncu is apparently low based on the surface appearance and limited measurements (see Appendix); but, perhaps the gas discharge is diffuse (as suggested by the large satellite thermal anomaly, Fig. 10), or the gas emission is episodic (as suggested for Montserrat, with several thousand years of discharge occurring in a few years; Christopher et al., 2015).

Shallow Earthquakes

While ground deformation and thermal anomalies can be measured from space at Uturuncu and Lazufre and compared, it is more difficult to compare the seismic activity of these volcanoes. Uturuncu and Lazufre are both hundreds of kilometers from the nearest permanent seismic networks, and there are no teleseismic events (Mw greater than ∼5) with reliable depths recorded near these volcanoes (e.g., Devlin et al., 2012; Pritchard et al., 2014a). Local networks near both Lastarria and Uturuncu from the PLUTONS and related projects recorded hundreds of small events (Mw <3) and reveal a wealth of information about the processes occurring including volcanic tremor and low-frequency earthquakes (Lastarria only), earthquake swarms, non-double-couple earthquakes, and triggered seismicity (e.g., Jay et al., 2012; Alvizuri and Tape, 2016). Although there are fewer data available on shallow earthquakes at other volcanoes in the Central Andes, based on available reconnaissance-level observations, it seems that the seismic activity at Uturuncu and Lastarria is among the highest in the region (Pritchard et al., 2014a). In fact, to our knowledge, Uturuncu is the only seismically active, Pleistocene volcano in the Central Andes, although seismic activity has been detected near the Pampa Lirima geothermal project in the Central Andes (Arcos et al., 2011) and could be related to the Pliocene Sillajhuay volcano (Pritchard et al., 2014a).

Remaining Questions

There are some questions that the PLUTONS project was not able to answer and new questions that have been raised by the research. Below, we describe these questions and suggest new observations and models that are needed to provide answers.

Temporal Variations in Ground Deformation and Seismicity

A single GPS station at Uturuncu spanning 2010–2015 shows a variable rate of uplift (Henderson and Pritchard, 2017). Similarly, ground deformation at the Lastarria-Lazufre magmatic system is seen to vary in time over the past two decades (e.g., Henderson et al., 2017). Future InSAR observations that span a large enough region to sample the entire deformation field (e.g., from the Sentinel-1a and 1-b satellites) and more GPS data (including reoccupation of campaign GPS sites installed by the PLUTONS project) are needed to better understand the temporal and possibly spatial variations in deformation that are used to inform models of the processes causing deformation. For example, our model of gas pressurization and/or depressurization suggests that subsidence of approximately equal magnitude to the uplift is needed to explain the lack of net geomorphic uplift at both Uturuncu and Lazufre and could occur in the coming decades.

There are no permanent seismometers near Uturuncu or seismic networks at Lazufre (although the Southern Andean Volcano Observatory [OVDAS] has two seismometers at Lastarria, on the Chilean side of the border), and so the characteristics and temporal variations in the earthquake activity at these volcanoes are not known. Jay et al. (2012) suggest some temporal variations at Uturuncu, but the data are very limited. In particular, it would be useful to understand if there are changes in the number, size, and type of earthquakes associated with the observed changes in ground deformation described in the last paragraph, or to determine if magma might be moving closer to the surface. We recommend installation of small, permanent, seismic networks near both Uturuncu and Lazufre.

Seismic and MT Anisotropy

Partially molten zones contain significant spatial heterogeneity that can be detected as anisotropy by seismic and electromagnetic (EM) wave propagation. To unravel the melt geometry and determine the amount of partial melt, observations of seismic anisotropy (e.g., Jaxybulatov et al., 2014; Spica et al., 2017), along with lab experiments and field observations, are required (e.g., Hammond and Kendall, 2016), and seismic anisotropy can also be used to understand the orientation of faults and fractures (e.g., Baird et al., 2015). Previous work in the APVC has suggested seismic anisotropy in the crust (Leidig and Zandt, 2003), particularly in the upper 7 km (Maher, 2017), as well as in the mantle (Sims et al., 2015), but further work is required, including studies at Lazufre. Magnetotelluric anisotropy has been studied in terms of understanding the need for a 3D inversion (e.g., Comeau et al., 2015, 2016; Díaz et al., 2015) but may also provide independent constraints on the amount of partial melt (e.g., as suggested for the mantle; Pommier et al., 2015). Even with the detailed observations from PLUTONS, the observations were logistically and financially constrained and so could be improved by better networks of geophysical observations (denser networks and/or spanning a larger spatial scale). For example, additional MT stations to establish a grid instead of profiles at both Uturuncu and Lazufre would help to better understand the 3D subsurface characteristics.

Expanded Data Collection at Lastarria-Lazufre

Since this area was the secondary target for the PLUTONS project, the study is only a reconnaissance effort. In particular, a denser network of stations spanning the international border is needed to better illuminate the earthquake activity and image the subsurface below the Lazufre uplift swell. A complete geochemical study of the Azufre, Cordón del Azufre, and Bayo complexes, including Los Colorados ca. 9 Ma ignimbrite and lavas within the caldera (Naranjo et al., 2018b), will help better determine the age of this complex structure, whether the volcanic vents are related, and the characteristics of the zone of melt generation (e.g., Froger et al., 2007; Perkins et al., 2016a).

Degassing Emissions and Composition

No constraints on gas composition or flux are available at Uturuncu (see Appendix), and only temporally discrete observations of these parameters are available at Lastarria. Considering our proposed model for ground deformation involves gas pressurization and/or depressurization, measuring the discharge of gases and their temporal variations would be important tests of this model. Specifically, measurements of diffuse degassing should be made at both Uturuncu and Lazufre and more broadly in the Central Andes, including deforming centers such as Cerro Overo (whose cycles of uplift and subsidence could also be related to degassing) and Cerro Blanco as well as non-deforming volcanoes.

Eruptive History of Volcanic Systems in the Central Andes

Uturuncu is one of the few volcanic centers in the region where the lava flows have been dated and their volumes have been measured to calculate the temporal variations in erupted volume. Do other centers above large zones of partial melt, such as the APMB and SPMB, show the same temporal variations, reflecting perturbations to these large zones of partial melt, or do they reflect smaller, individual characteristics (e.g., erosion or flank instabilities of the edifice promoting eruption)?

Joint Geophysical and Petrological Inversion

As demonstrated in Figure 2, by combining different geophysical inferences of the subsurface, we can understand the limitations of each individual data set to develop a unified model of subsurface characteristics. For example, a zone of low seismic velocity between 9 and 18 km below the surface of Uturuncu could have a relatively high resistivity if partial melt in that region is resistive dacite (e.g., Comeau et al., 2016). The next step forward is to combine various geophysical data sets including seismic, MT, and gravity data in a formal joint inversion that accounts for the different sensitivities of each technique (e.g., Moorkamp et al., 2011). But, for these joint inversions to succeed, we need experimentally determined relations between subsurface properties, such as composition, geometry, and the amount of partial melt and observations of seismic velocity, density, and resistivity. Laboratory measurements of rock properties have been extremely valuable—for example, they have shown that basaltic melts have lower resistivity than sometimes assumed (e.g., Pommier and Garnero, 2014) and that dacitic melt actually has a relatively high resistivity (Laumonier et al., 2015). However, laboratory experiments have not explored all relevant portions of parameter space and should consider the effects of three-phase systems on resistivity, Vp, Vs, including the role of CO2 and other non-aqueous volatiles. Further, accurate assessment of the amount of melt must consider anisotropy (mentioned above) as well as the effects of smoothing on geophysical inversions that can underestimate melt percentages by more than 20% (e.g., Paulatto et al., 2012).

Geomorphology and Numerical Modeling

Our interpretation that the uplift at Uturuncu and Lazufre is related to cycles of gas pressurization and/or depressurization could be challenged by new observations or models. Specifically, new observations of net geomorphic uplift over the past several thousand years as found at Laguna del Maule, Chile (Singer et al., 2015) would challenge the gas pressurization and/or depressurization model. Because we currently do not have evidence of net uplift over these timescales, we think that a diapir model is no longer favored for Uturuncu (e.g., Perkins et al., 2016a), but future numerical models of diapir ascent may show decadal variations in uplift rate (e.g., related to spatial variations in melt fraction) that could produce no net deformation over millennia. In particular, we need self-consistent models that include not just movement of partial melt and pressurization and/or depressurization of reservoirs but those that model the deformation effects of thermal and thermodynamic consequences of fluid movements in a viscoelastic material—for example, anatexis, heating and then cooling of surrounding country rock and crystallization. These models would include the temperature and heat-flow evolution with time as activity in the shallow dacite chamber waxes and wanes; such activity can be compared with petrological and heat-flow measurements.

Temperature and Rheology

To constrain the numerical models described in the last paragraph, we need better constraints on the subsurface temperatures and rheology at both Uturuncu and Lazufre. Observations of heat flow in both locations are valuable and could be used to constrain different inferences about the subsurface including whether the physical cause of the low seismic velocity zones are caused by partial melt or something else (e.g., comparing 150 and 200 mW/m2 lines in Fig. 2D). An independent constraint on the subsurface temperature could be provided by using Vp and Vs to determine the depth of the alpha-beta quartz transition or the Curie temperature, as done elsewhere (e.g., Hacker et al., 2014; Saibi et al., 2015). In terms of subsurface rheology, the viscosity of the partially molten zone can be estimated from the resistivity derived from the MT data (Comeau et al., 2016). However, the bulk effective viscosity on timescales of years to millennia is also needed for modeling the geodetic and geomorphic observations. This viscosity could be estimated by measuring the lithospheric response to unloading of large paleolakes, such as Lake Minchin, which is just to the north of the APMB (e.g., Bills et al., 1994). Finally, seismic attenuation (e.g., Farrell et al., 2017) could be used to estimate viscosity (e.g., Priestley and McKenzie, 2013).

Relationship between Uturuncu and the APVC Ignimbrites

The relationship between magmatism at Uturuncu and the ignimbrite-dominated activity of the APVC is not fully understood. The spatial and temporal connection between Uturuncu and the APVC is obvious, and the calc-alkaline high-K dacites and minor andesitic compositions erupted at Uturuncu share many characteristics with those of the APVC—that they are ultimately derived from APMB is not in question (Muir et al., 2014a, 2014b; de Silva et al., 2006). However, Michelfelder et al. (2014) showed that Uturuncu is an outlier compared to other Pleistocene arc front and other composite volcanoes in the CVZ. Bulk-rock isotope data and the Sr isotopes from plagioclase phenocryst cores from Uturuncu are more radiogenic compared to regional arc dacites and andesites but concordant with those from the APVC (particularly the eastern APVC ignimbrites, such as Panizos), although Uturuncu dacites do not share the peraluminous character of these ignimbrites. Model entrapment pressures based on plagioclase-hosted melt inclusions for Uturuncu (Muir et al., 2014a, 2014b) and those from quartz-hosted melt inclusions from APVC dacites (Schmitt et al., 2001; Grocke et al., 2018) show broad concordance with the ignimbrite magmas extending to slightly greater pressures. There are some important differences. Uturuncu’s anhydrous assemblage of plagioclase, orthopyroxene, and biotite is a subset of the APVC magmas that typically have plagioclase, quartz, amphibole, and biotite with some units containing sanidine and others that contain orthopyroxene. These differences might be reconciled with local differences in the source and local storage regions, but these remain to be investigated, and the differences are not well understood. As also noted by Muir et al. (2014a, 2014b), the broad similarity of Uturuncu and APVC dacites brings to our attention that we need to better understand how “eruptions of vastly differing magnitude and energy can be sourced from magmas with similar compositions and similar dissolved volatile contents stored at similar depths in the upper crust.”

Regional and Global Studies of Magmatic Systems

The PLUTONS project is one of a handful that has used a suite of detailed, coincident geochemical and geophysical observations to study the entire crustal expression of a magmatic system (e.g., see summary in Pritchard and Gregg, 2016). To better understand the significance of the results at Uturuncu and Lazufre, similar measurements are needed at other locations in the Andes and globally, to determine for example, how unique are the amounts and extent of partial melt? For example, are the low-velocity zones imaged in the lower crust extending from the APMB to perhaps the Moho (Kukarina et al., 2017) and the lower crustal reflector beneath Uturuncu (McFarlin et al., 2018) long-lived or transient features? Do they reflect a partially molten pathway that allows magma to ascend from the mantle to the mid-crust through porous flow as proposed for Hawaii (e.g., Pritchard et al., 2007), feeding the APMB? If so, has it always been located below Uturuncu, or does it move, depending on heat and melt supply from the mantle? Is the vertical anomaly from earthquake tomography below Uturuncu unique? Is this the only “feeder” to the APMB? Do other large partially molten zones (such as the SPMB) have a similar feature? Is this a long-lived feature, or does it move around over millions of years? Is there an APMB- or SPMB-like feature elsewhere in the northern Andes (Cotopxi or Nevado de Ruiz), southern Andes (Laguna del Maule or Cordón Caulle), or elsewhere in the world?

Note that if the APMB is mostly silicic and andesitic, then there must have been at least an equal volume of deeper lower-crustal complementary mafic and ultramafic cumulates—while some of this material may still be evident below the lower-crustal reflector observed by receiver functions (Fig. 5; McFarlin et al., 2018), some may have been lost by the delamination that is necessary to provide the heat to maintain a partially molten APMB for millions of years.

Some have posited that the APMB has been partially molten for the entire duration of the ignimbrite flare-up, although Babeyko et al. (2002) argue that to keep temperatures high enough in the mid-crust requires convection upwellings from the lower crust to transport heat from the mantle to the mid-crust. Maybe the low-velocity region imaged by Kukarina et al. (2017) shows the upwelling limb of the convection pattern in the lower crust that was suggested by Babeyko et al. (2002)? Partial melts moving upward would be partially compensated by the subsidence of dense refractory products of partial meting and fractional crystallization as envisaged in models of compaction and larger-scale crustal dynamics (e.g., Solano et al., 2012; Roman and Jaupart, 2016) that may be imaged by receiver functions in the lower crust (Fig. 5; McFarlin et al., 2018) at the bottom of the vertical low-velocity zone (Kukarina et al., 2017).

Alternatively, the uppermost APMB may have grown with time as reflected by the spatial, temporal, and volumetric pattern of ignimbrite eruptions through time (e.g., de Silva and Gosnold, 2007; Kern et al., 2016). Linking the record of the flare-up to the deeper APMB remains a challenge. An integrated numerical model of the evolution of the APMB over tens of millions of years constrained with PLUTONS data showing the newly imaged structures in the lower crust that could address the questions raised here has yet to be developed, and the fundamental ambiguity of the relevance of the geophysical observations to the >10 m.y. evolution of the APMB remains.

By combining multiple geophysical and geochemical methods, we have constrained the crustal-scale distribution of partial melt at Lazufre and Uturuncu. In the lower crust below Uturuncu, there is a 20–30-km-wide vertical low Vs zone (Kukarina et al., 2017) that possibly extends from ∼40 km deep to the mid-crust, where it overlaps with the 200-km-diameter APMB. Beneath 40 km and the vertical low Vs zone is a zone of mafic cumulates in the mid- to lower crust that is visible on the receiver functions (McFarlin et al., 2018). We interpret the 7–10 km difference in depths between the top of the APMB from seismic tomography and/or receiver functions (Ward et al., 2014b) and from magnetotellurics to indicate that the APMB is compositionally zoned, with an upper region of dacitic melt (9–18 km below the surface) and andesitic melts below (18–30 km below the surface) (Fig. 2; Table 2). There are several zones of low resistivity between the APMB and the surface (Comeau et al., 2015, 2016), including one at ∼5 km depth below Uturuncu itself; this zone is likely related to the current hydrothermal system developed in the region of the dacite magma chambers that supplied eruptions of Uturuncu before 250 ka (Muir et al., 2014a, 2014b, 2015). The presence of such a system and the likelihood that it contains brines may have implications for ore formation in this area because brines are known to strongly sequester and concentrate metals such as Cu and Mo (Candela and Holland, 1984; Simon et al., 2006; Frank et al., 2011; Zajacz et al., 2011; Tattitch et al., 2015; Tattitch and Blundy, 2017).

The PLUTONS project demonstrates the need for multiple geophysical and petrological techniques in order to constrain a crustal cross section such as that presented in Figure 2. All geophysical methods are non-unique (e.g., low seismic velocities can be caused by changing porosity, partial melt, composition, volatile content, temperature, attenuation, anisotropy, etc.), and some methods have “blind spots” (e.g., dacite melt has relatively high resistivity or MT cannot easily detect melt that is not well interconnected). Furthermore, the different geophysical techniques have different capabilities to resolve structures at different spatial resolutions—for example, long-wavelength seismic waves (Fig. 4D) averages over subsurface features that can be resolved by higher-resolution methods (Figs. 4C, 4E, and 4F), but these high-resolution methods may, in turn, miss the larger-scale structures (Fig. 4E). All geophysical techniques provide a snapshot of conditions in the subsurface at the present and raise questions about their temporal evolution (e.g., Ward et al., 2017). We do not know how the geometry of the APMB or SPMB has changed with time, except by inferences from the record of surface volcanic eruptions (e.g., Kay and Coira, 2009; Kern et al., 2016). Further, we do not know if structures imaged in the lower crust are long-lived or transient, but we have concluded that the “fossil dacite chamber” at Uturuncu once contained eruptible dacitic magma, although we do not think that there is partial melt there now—only an active hydrothermal system.

The Central Andes is unique among volcanic arcs in having very large zones of crustal partial melt (APMB and SPMB, total >0.5 million km3) at >10 km depth. These zones are a product of the unusually thick crust and a history of subduction (e.g., Kay and Coira, 2010; Beck et al., 2015). Uturuncu and Lazufre are not the only deforming or seismically active volcanoes possibly related to these zones of partial melt, but they are globally unique in terms of the size of the deforming regions, the decadal timescales over which they have been deforming, and the ambiguous relationship between the deformation and seismic activity and Holocene eruptive centers. However, our observations of ground deformation in these areas span only a few decades, and considering the repose interval of >105 years, our interpretation of the cause of deformation requires a longer period of observation. This longer time period can be achieved through geomorphic observations of net uplift and subsidence or through the ergodic hypothesis that we can substitute space for time—which allows us to observe hundreds of other magmatic systems in different stages of their eruptive cycle to better understand the complete cycle. By combining these approaches, the PLUTONS project motivated a new interpretation for the cause of deformation. Because of the temporal variations in the deformation rate recorded by geodesy (e.g., Henderson and Pritchard, 2017; Henderson et al., 2017) and the lack of long-term uplift from geomorphology (Perkins et al., 2016a), we do not favor en masse diapiric movements. Instead, we suggest that the current uplift is a transient caused by exsolved volatiles that are being temporarily trapped at depths of up to 10–15 km and that will eventually be released with concomitant subsidence, causing no net deformation (Gottsmann et al., 2017). This model is consistent with other magmatic systems in the Central Andes such as Cerro Overo, which has been undergoing alternating periods of uplift and subsidence, or Cerro Blanco subsiding for decades (e.g., Henderson and Pritchard, 2013). This model can be tested by longer-term geodetic observations, time-series observations of total gas emissions, and viscoelastic numerical models of deformation that incorporate self-consistent compositional and thermodynamic characteristics and temporal evolution.

While Uturuncu and Lazufre share some characteristics with caldera systems such as Yellowstone and Cordón Caulle, there are significant differences in their manifestations of deformation, surface thermal features, and seismicity, possibly related to their hydrothermal systems, influence of faults, and the depth of partially molten rock. There might be partial melt at a depth of 3 km beneath Lastarria, but this needs to be tested with heat-flow observations. Furthermore, this magma chamber is relatively small (∼81 km3) and is on average <40% partial melt; so it is not considered eruptible (Spica et al., 2015). At Uturuncu, partial melt of ∼10% may be present at 9 km depth (Ward et al., 2014b), but eruptible quantities do not exist in large volumes (e.g., Pritchard and Gregg, 2016). At this time, there is no indication that the large volumes of partial melt observed beneath Uturuncu or Lazufre will erupt in the future. However, further monitoring of these volcanoes, including the installation of permanent seismic arrays, is suggested to better understand the changing dynamics.

Funding for this project was provided by the U.S. National Science Foundation through grants EAR-0908850 to N.J.F., EAR-0908324 to S.L. de S., EAR-0901148 to Todd Feeley, EAR-0909254 to University of Alaska, Fairbanks, and EAR-0908281 to Cornell with a subaward to the University of Alberta, the UK Natural Environmental Research Council (grant NE/G01843X/1), a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery grant awarded to M.J.U., and a Benjamin Meaker Visiting Professorship to M.E.P from the Institute for Advanced Study, University of Bristol. Total funding from all sources was about $4.5 million. We thank project collaborators who contributed to the project without direct PLUTONS funding. This project would not have been possible without the help of the Observatorio San Calixto, Bolivian Servicio de Areas Protegidas (especially the staff at the Reserva Eduardo Avaroa), the residents of Quetena Chico and Quetena Grande, the Servicio Nacional de Geología y Técnico de Minas, SERNAGEOMIN Chile, and our outstanding field drivers. We dedicate this paper and themed issue of Geosphere to the memory of our colleague and PLUTONS co-principal investigator, Todd Feeley (Michelfelder et al., 2016). We thank Robert Trumbull and an anonymous reviewer for critical reviews that improved the manuscript. Data from the PLUTONS project are archived for future use—seismic waveforms are at the IRIS Data Management Center, raw GNSS data are at UNAVCO, and access to other data sets are described in the respective papers.



1. Gravity: 160 new survey (or static) benchmarks for a regional map (del Potro et al., 2013) and 25 sites that were reoccupied twice to assess temporal variations in gravity (called dynamic gravimetric data; Gottsmann et al., 2017).

2. Ground Deformation: Interferometric synthetic aperture radar time series from 1992 to 2011, two functioning continuous GPS, 2010–present (plus an additional site that was vandalized after one year; Henderson and Pritchard, 2017). GPS reoccupation in 2010 of leveling benchmarks from 1965 (Gottsmann et al., 2018).

3. Seismometry: 28 broadband seismometers from 2010 to 2012 (e.g., Farrell et al., 2017) have been used for: (1) joint receiver function and ambient noise tomography imaging using a subset of PLUTONS stations with others (Ward et al., 2014b); (2) earthquake tomography and earthquake locations (Kukarina et al., 2017); (3) P-wave attenuation and time delays (Farrell et al., 2017); (4) detailed studies of receiver functions using all PLUTONS data (McFarlin et al., 2018); (5) full moment tensors for 63 shallow events beneath the volcano (Alvizuri and Tape, 2016); and (6) a one-dimensional velocity model based on receiver functions, Rayleigh-wave horizontal-vertical amplitude ratios, and Rayleigh-wave phase velocities from earthquakes and ambient noise (Shen et al., 2017). A related project used 15 short- and intermediate-period seismometers at Uturuncu from 2009 to 2010, and Jay et al. (2012) describes earthquakes during that time and provides a 3D ambient noise tomographic map of Vs velocities in the upper 5 km.

4. Magnetotellurics: 180 broadband stations occupied for about one day per station during field campaigns in 2011, 2012, and 2013. Inversion for the subsurface resistivity structure was done using all stations in a 2D inversion (Comeau et al., 2015), as well as 73 (Comeau et al., 2015) and 96 stations (Comeau et al., 2016) in 3D inversions.

5. Geomorphology: Deflections of river longitudinal profiles and shorelines of lakes around Uturuncu were assessed with GPS and a 35 × 35 km 0.5 m/pixel photogrammetry-derived digital elevation model—selected shoreline features were dated using optically stimulated luminescence (Perkins et al., 2016a). Elevation of basement exposures and a long-wavelength uplift over the APMB have been analyzed to determine the extent of isostatic compensation and potential magmatic contributions to the modern surface topography (Perkins et al., 2016b).

6. Geochronology and Isotope Geochemistry:Kern et al. (2016) used in situ secondary ionization mass spectrometry (SIMS) analyses of U-Pb spot ages and trace elements in zircon to assess the magmatic evolution of 16 ignimbrites from the APVC for regional context. The volcanic output from Uturuncu has been constrained by 26 new 39Ar/40Ar plateau ages for 23 lavas and domes (Muir et al., 2015). Whole-rock analysis of 121 Uturuncu lava flow and dome samples for major- and trace-element concentrations were determined by X-ray fluorescence spectrometry (XRF). A subset of 60 samples was analyzed for trace-element concentrations including rare-earth elements by inductively coupled plasma mass spectrometry (ICP MS), Nd, Pb, and Sr isotopes by thermal ionization mass spectrometry (TIMS), and δ18O values by laser fluorination (Michelfelder, 2015). Michelfelder (2015) provide major-, trace-element, and in situ Sr-isotope data on Uturuncu plagioclase phenocrysts combined with the whole-rock Sr-isotope ratios.

7. Petrology and Geobarometry:Grocke et al. (2018) analyzed samples from the 3.64 Ma Tara ignimbrite (from the 5.6–1.8 Ma Guacha caldera immediately south of Uturuncu) for whole-rock major and trace elements and isotopic compositions. Major- and trace-element concentrations of minerals, matrix glasses, and melt inclusions by electron microprobe, as well as the H2O and CO2 concentrations in quartz- and sanidine-hosted melt inclusions using secondary ionization mass spectrometry (SIMS), laser ablation (LA) ICP MS, and Raman spectroscopy have been analyzed. Muir et al. (2015) analyzed Uturuncu lavas for major-element concentrations of minerals, matrix glasses, and melt inclusions by electron microprobe, as well as the H2O and CO2 concentrations in plagioclase-hosted melt inclusions from lava flows by SIMS. Petrological phase equilibria experiments were run at a range of pressure and temperature conditions and compared with natural phase assemblages, modes, and mineral and glass compositions to understand the storage conditions of the dacitic magmas (Muir et al., 2014a, 2014b).

8. Gas Emissions: Measurements of gas flux and composition were attempted during the PLUTONS project but are not published elsewhere. We mention the available measurements here to motivate future work: Georg Zellmer (2011, personal commun.) took MultiGAS measurements in 2009 and noted diffuse degassing, with hundreds of tiny vents covering a large area close to the saddle of the mountain. He collected data during only a few brief periods that provide a rough estimate of the gas composition: CO2 >6000 ppm (off scale); H2S peaking at just above 80 ppm; SO2 <170 ppm (though in most cases less than H2S). Relative humidity went up to 100% (at 10 °C), confirming significant water in this otherwise dry region. In 2012, Taryn Lopez used a FLYSPEC (scanning ultaviolet spectrometer) at the highest seismic station on Uturuncu (1–2 km from the degassing area), but SO2 was below the detection limit, implying the emission rate is probably less than ∼10–15 t/d. The shallow hydrothermal system at Uturuncu could remove the highly water-soluble SO2 from the gas phase and is consistent with the relatively high proportion of H2S to SO2. On the other hand, there is abundant sulfur on the surface at Uturuncu and a former sulfur mine a few km from the summit, suggesting higher deposition of sulfur in the past.


1. Seismology: Eight broadband seismometers deployed by PLUTONS from 2011 to 2013—a three-month subset of these data were combined with a two-month data set from 17 seismometers installed by Chilean and German researchers in 2008 (Spica et al., 2015) for ambient noise tomography mapping of the seismic shear-wave velocity variations beneath Lastarria. Analysis of these data for local earthquakes is ongoing (Heather McFarlin, 2018, personal commun.).

2. Magnetotellurics: The University of Alberta group collected 38 broadband MT stations in an EW profile through Lazufre during a field campaign in 2013. An additional 30 broadband MT sites from Díaz et al. (2015) and 11 broadband and long-period MT stations from Budach et al. (2013) focused around Lastarria are also included in the inversion (Comeau, 2015).

3. Geochronology and Field Mapping: During the past five years, the Chilean Geological Survey (Sernageomin) has created 1:100,000 maps along the volcanic arc between 25° and 27°S, including Lastarria and Lazufre, with 110 new 39Ar/40Ar ages to go along with the ∼300 K-Ar ages from the 1990s (Naranjo and Cornejo, 1992; Naranjo et al., 2013).

4. Geomorphology:Perkins et al. (2016a) describe shorelines at a late Pleistocene lake near Lastarria volcano that could indicate the longevity of deformation. These authors also use the orientation of multiple dated lava flows to date the age of uplift of a topographic swell co-located with Lazufre, and they map the orientation of volcanic vents around the swell.

5. Ground Deformation: There are two functioning continuous GPS in the Lazufre deformation zone, with data between 2010 and 2014 at one station, and 2010–2015 at the other, plus a reference site 150 km from the deformation with data from 2011 to 2015 (Henderson et al., 2017). An additional reference site was vandalized, and useful data were not collected. There are ten campaign GPS sites installed by PLUTONS and occupied in 2012–2014, and seven other sites are also available from 2006, 2007, and 2008 from another project (Remy et al., 2014). An InSAR time series of ground deformation has been created from 1995 to 2011, and additional interferograms constrain deformation through 2016 (Henderson et al., 2017).

6. Petrology and Geochemistry: Lava flows were sampled at Lastarria (69 samples) and Cordón del Azufre (21), and the major- and trace-element geochemistry was studied by XRF, while a subset was analyzed for trace-element concentrations including rare-earth elements by ICP MS. Nd-, Pb-, and Sr-isotope ratios were measured by TIMS for 15 whole-rock samples (Wilder, 2015).

7. Gas Emissions: At Lastarria, fumarolic gases were sampled and analyzed by Aguilera et al. (2012) and Lopez et al. (2018), and vapor condensates and fumarolic solid samples (including fluid inclusions) were analyzed by Aguilera et al. (2016).

Science Editor: Raymond M. Russo
Gold Open Access: This paper is published under the terms of the CC-BY-NC license.