After ca. 250 kyr without a known eruption, in recent decades Uturuncu volcano in Bolivia has exhibited multiple signs of unrest, making the classification of this system as “active”, “dormant”, or “extinct” a complex question. Previous work identified anomalous low resistivity zones at <10 km depth with ambiguous interpretations. We investigate subsurface structure at Uturuncu with new gravity data and analysis, and compare these data with existing geophysical data sets. We collected new gravity data on the edifice in November 2018 with 1.5 km spacing, improving the resolution of existing gravity data at Uturuncu. Gradient analysis and geophysical inversion of these and older gravity data revealed a 5 km diameter, positive density anomaly beneath the summit of Uturuncu (1-3 km elevation) and a 20 km diameter arc-shaped negative density anomaly around the volcano (-3 to 4 km elevation). These structures often align with resistivity anomalies previously detected beneath Uturuncu, although the relationship is complex, with the two models highlighting different components of a common structure. Based on a joint analysis of the density and resistivity models, we interpret the positive density anomaly as a zone of sulfide deposition with connected brines, and the negative density arc as a surrounding zone of hydrothermal alteration. Based on this analysis we suggest that the unrest at Uturuncu is unlikely to be pre-eruptive. This study shows the value of joint analysis of multiple types of geophysical data in evaluating volcanic subsurface structure at a waning volcanic center.

Uturuncu volcano in southern Bolivia is a member of a distinctive class of volcanoes – systems that show unrest despite not having erupted in the Holocene. Uturuncu has not erupted in 250 kyr, but has been deforming (uplift with a moat of subsidence) for several decades, along with seismic swarms and active, sulfur-encrusted fumaroles. Our work builds on previous geophysical imaging at Uturuncu by jointly analyzing multidisciplinary datasets, focusing on imaging the shallow (<15 km depth below surface) structure of the system with geophysical and geochemical data. Whereas previous research pointed to andesite melt at depths >15 km depth, results were ambiguous as to what proportions of melts vs. brines are present at shallower depths. Identifying fluids (melt, brine, etc.) and structures at shallow depths is key for evaluating the hazard potential of the volcano and understanding the source of the unrest. We present new results from gravimetry, seismology (hypocenter relocation, seismic velocity and attenuation tomography), gas geochemistry, and InSAR observations. The results point to an extensive and active hydrothermal system extending ~20 km laterally and ~10 km vertically from Uturuncu, with possible connections at depth to the deeper magmatic system. A combined view of the new density, seismic velocity and attenuation models, and the existing resistivity model is crucial for revealing key features of the hydrothermal system: a vapour-rich conduit beneath Uturuncu (low resistivity/high attenuation column extending from 1.5 to 12.5 km depth), an extensive alteration zone surrounding Uturuncu (complex zone of annular shaped anomalies surrounding Uturuncu from 1.5 to 12.5 km depth), and a possible zone of sulfide deposition just below the western flank of Uturuncu at 1.5 km depth (high density/low resistivity/high attenuation). High fluxes of diffuse CO2 degassing at sub-magmatic temperatures and a small area directly above a low resistivity anomaly subsiding from 2014 to 2017 show that the hydrothermal system is currently active. Analyzed jointly, this multidisciplinary data set suggests that current activity within the shallow structure at Uturuncu is dominated by hydrothermal, rather than magmatic processes.

We present evidence of volcano-tectonic interactions at Sabancaya volcano that we relate to episodic magma injection and high regional fluid pore pressures. We present a surface deformation time series at Sabancaya including observations from ERS-1/2, Envisat, Sentinel-1, COSMO-SkyMed, and TerraSAR-X that spans June 1992 - February 2019. These data show deep seated inflation northwest of Sabancaya from 1992-1997 and 2013-2019, as well as creep and rupture on multiple faults. Afterslip on the Mojopampa fault following a 2013 Mw 5.9 earthquake is anomalously long-lived, continuing for at least six years. The best fit fault plane for the afterslip is right-lateral motion on an EW striking fault at 1 km depth. We also model surface deformation from two 2017 earthquakes (Mw 4.4 and Mw 5.2) on unnamed faults, for which the best fit models are NW striking normal faults at 1-2 km depth. Our best fit model for a magmatic inflation source (13 km depth, volume change of 0.04 to 0.05 km^3 yr^-1), induces positive Coulomb static stress changes on these modeled fault planes. Comparing these deformation results with evidence from satellite thermal and degassing data, field observations, and seismic records, we interpret strong pre-eruptive seismicity at Sabancaya as a consequence of magmatic intrusions destabilizing tectonic faults critically stressed by regionally high fluid pressures. High fluid pressure likely also promotes fault creep driven by static stress transfer from the inflation source. We speculate that strong seismicity near volcanoes will be most likely with high pore fluid pressures and significant, offset magmatic inflation.

Uturuncu volcano is situated in the Bolivian Andes, directly above the world’s largest crustal body of silicic partial melt, the Altiplano-Puna Magma Body (APMB). Uturuncu last erupted 250,000 years ago, yet is seismically active and lies at the centre of a 70 km diameter uplifted region. Here, we analyse seismicity from 2009 to 2012. Our earthquake locations, using a newly developed velocity model, delineate the top and bottom of the APMB, reveal individual faults, and reconcile differences in depth distribution between previous studies. Spatial clustering analysis of these earthquakes reveals the orientations of the faults, which match stress orientations from seismic anisotropy. Earthquake b-values derived from moment magnitudes (1.4) differ significantly from those using local magnitude measurements (0.8). We suggest that, if possible, moment magnitudes should always be used for accurate b-value analysis. We interpret b-values > 1 in terms of fluid-enhanced seismicity. Shallow seismicity local to Uturuncu yields b-values > 1.1 with some temporal variation, suggesting fluid migration along pre-existing faults in a shallow hydrothermal system, likely driven by advection from the APMB. Intriguingly, events deeper than the APMB also yield large b-values (1.4), mapping the ascent into the lower crust of fluids originating from a subducting slab. Cumulatively, these results provide a picture of an active magmatic system, where fluids are exchanged across the more ductile APMB, feeding a shallow, fault-controlled hydrothermal system. Such pathways of fluid ascent may influence our understanding of arc volcanism, control future volcanic eruptions and promote the accumulation of shallow hydrothermal ore deposits.

Mapping fluid accumulation in the crust is pertinent for numerous applications including volcanic hazard assessment, geothermal energy generation and mineral exploration. Here, we use seismic attenuation tomography to map the distribution of fluids in the crust below Uturuncu volcano, Bolivia. Seismic P-wave and S-wave attenuation, as well as their ratio (QP/QS), constrain where the crust is partially and fully fluid-saturated. Seismic anisotropy observations further constrain the mechanism by which the fluids accumulate, predominantly along aligned faults and fractures in this case. Furthermore, subsurface pressure-temperature profiles and conductivity data allow us to identify the most likely fluid composition. We identify shallow regions of both dry and H2O/brine-saturated crust, as well as a deeper supercritical H2O/brine column directly beneath Uturuncu. Our observations provide a greater understanding of Uturuncu’s transcrustal hydrothermal system, and act as an example of how such methods could be applied to map crustal fluid pathways and hydrothermal/geothermal systems elsewhere.

The Vavilov Ice Cap destabilized in 2013. It reached its highest annual ice loss rate of 4.48 km3/yr between 2015 and 2016, a rate that is more than half of the entire combined ice loss from all the other ice caps in the Russian Arctic. To understand the mechanics of how the surge took place and what will happen in the future, we investigate surface elevation and glacier velocities using the Cryosphere And Remote Sensing Toolkit (CARST), an open-access python toolbox designed for processing temporal changes of high-resolution remote sensing data. We use optical satellite images from WorldView, Landsat, and Sentinel-2 and their derived elevation and velocity products to track the history of the surge between 2010 and the present. We propose that the surge initiated when the ice front overrode weak marine sediments in 2013, leading to a reduction of frontal friction. Velocity time series show that the glacier reached a maximum speed of 25 m/day (9 km/yr) in late 2015, when a piedmont-like ice lobe stretched more than 10 km into the Kara Sea. However, in spring 2017 the glacier slowed down to 7-9 m/day with the development of a new channel inside the piedmont lobe, with new shear margins visible from optical imagery. The channelized flow pushes through the grounded portion of the piedmont glacier, and suggests a further reorganization of resistive forces. The unprecedented evolution of the surge at Vavilov Ice Cap shows a strong connection to the status of the terminus, which might pose general concern for the stability of marine- or lake-terminating glaciers regardless of their locations on Earth.