Seismic Tomography of Kamchatkan Volcanoes

—The Kamchatka Peninsula is one of the most tectonically active regions in the world, where intensive and diverse modern volcanic activity takes place. In the recent decade, substantial progress in the investigation of deep structures beneath Kamchatka has been achieved owing to numerous tomography studies based on seismological data provided by permanent stations and temporary networks deployed in some key areas. The goal of this review is summarizing and systematizing dozens of separate multiscale geophysical studies in Kamchatka and constructing an integral model of volcano-feeding systems. An important part of this review contains the description of results of various seismic studies related to the Klyuchevskoy group volcanoes, which can now be considered one of the best studied volcanic areas in the world. The results of the regional-scale seismic tomography reveal the existence of the Pacific slab window, which determines the particular activity of the Klyuchevskoy group volcanoes. Middle-scale tomography studies have found traces of an ascending hot mantle flow that passes through the slab window, reaches the bottom of the crust below Shiveluch Volcano, and then propagates laterally toward the Klyuchevskoy group. Seismic models of the entire crust in the area of the Klyuchevskoy group were used to identify different mechanisms of magmatic feeding of three most active volcanoes: Klyuchevskoy, Bezymianny, and Tolbachik. The data of local networks deployed on several volcanoes of Kamchatka were used to image the magma sources in the upper crust, which are directly responsible for the current eruption activity. The comparison of the results for the Kamchatka volcanoes with tomography models of several other volcanoes of the world allowed determining some common features and differences in feeding active magmatic systems.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume–like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a “mantle wind” toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Lateral variation in slab window viscosity inferred from global navigation satellite system (GNSS)–observed uplift due to recent mass loss at Patagonia ice fields

The geographic coincidence of the Chile Ridge slab window and the Patagonia ice fields offers a unique opportunity for assessing the effects of slab window rheology on glacial isostatic adjustment (GIA). Mass loss of these ice fields since the Little Ice Age causes rapid but variable crustal uplift, 12–24 mm/yr around the North Patagonia ice field, increasing to a maximum of 41 mm/yr around the South Patagonia ice field, as determined from newly collected or processed geodetic data. We used these observational constraints in a three-dimensional Maxwell viscoelastic finite element model of GIA response above both the subducting slab and slab window in which the upper-mantle viscosity was parameterized to be uniform with depth. We found that the viscosity of the northern part of the slab window, ~2 × 1018 Pa·s, is lower than that of the southern part by approximately an order of magnitude. We propose that this along-strike viscosity contrast is due to late Cenozoic ridge subduction beneath the northern part of the slab window, which increases asthenospheric temperature and reduces viscosity.

Supplemental Material: Lateral variation in slab window viscosity inferred from global navigation satellite system (GNSS)–observed uplift due to recent mass loss at Patagonia ice fields

Modeling method, Table S1 (newly obtained uplift rates for the ten continuous GNSS sites shown in Figure1A),Table S2 (slab window viscosity values and RMS data), supplemental figures, and the weekly solutions for GNSS data used in the study.<br>

Supplemental Material: Lateral variation in slab window viscosity inferred from global navigation satellite system (GNSS)–observed uplift due to recent mass loss at Patagonia ice fields

Modeling method, Table S1 (newly obtained uplift rates for the ten continuous GNSS sites shown in Figure1A),Table S2 (slab window viscosity values and RMS data), supplemental figures, and the weekly solutions for GNSS data used in the study.<br>

Liquiñe-Ofqui’s fast slipping intra-volcanic arc crustal faulting above the subducted Chile Ridge

The southernmost portion of the Liquiñe-Ofqui fault zone (LOFZ) lies within the proposed slab window which formed due to oblique subduction of the Chile Ridge in Patagonia. Mapping of paleo-surface ruptures, offsets, and lithological separations along the master fault allowed us to constrain geologic slip rates for the first time with dextral rates of 11.6–24.6 mm/year (Quaternary) and 3.6–18.9 mm/year (Late-Cenozoic) respectively. We had trouble mapping the LOFZ in one local because of a partially collapsed and previously undiscovered volcanic complex, Volcan Mate Grande (VMG: 1,280 m high and thus Vesuvius-sized) that grew in a caldera also offset along the LOFZ and has distinct geochemistry from adjacent stratovolcanoes. Besides the clear seismic and volcanic hazard implications, the structural connection along the main trace of the fast slipping LOFZ and geochemistry of VMG provides evidence for the slab window and insight into interplay between fast-slipping crustal intra-arc crustal faults and volcanoes.

Differential uplift of three Pliocene sea level indicator sites in southern Argentina driven by upwelling asthenosphere through the Patagonian slab window

&lt;div&gt; &lt;p&gt;Bivalve and gastropod shell beds deposited during the Early Pliocene (4.69-5.23 Ma) occur in uplifted outcrops (36 &amp;#8211; 180 m above sea level) along the east coast of Patagonian Argentina. These rock units provide a record of sea level during a geologic period when atmospheric CO2 and temperatures were higher than today. As such, reconstructing the elevation of global mean sea level (GMSL) during this time allows us to better understand how sensitive ice sheets are to increased past and future warming. However, reconstructing GMSL from local sea level indicators is hindered by effects such as mantle dynamic topography and glacial isostatic adjustment (GIA) that cause local sea level to deviate from the global mean. Here we use geodynamic modeling to better understand this complex dynamic setting and quantify the amount of uplift along this coastline.&lt;/p&gt; &lt;/div&gt;&lt;div&gt; &lt;p&gt;Despite being located on a relatively stable passive margin, significant variations in the elevation of the paleo shoreline indicators imply that the underlying convecting mantle is deforming the coastline. In particular, the subduction of the Chile Rise beginning ~18 Ma beneath Patagonia has generated a slab window underneath this region through which hot asthenosphere ascends. However, the former slab is still present deeper in the mantle, which causes a complex interplay between the downwelling slab and the upwelling asthenosphere. To quantify the effects of dynamic topography change since the Pliocene, we run 3D mantle convection simulations using the code ASPECT. We initialize our global model with a composite temperature structure derived from recent tomographic studies and a calibrated parameterization of upper mantle anelasticity. Independent estimates of pressure and temperature from thermobarometric calculations of proximal Pali-Aike xenoliths agree with the thermal structure of the tomography-based Earth model. We back-advect temperature perturbations and extract the resulting change in dynamic topography. Pairing GIA models and a suite of convection simulations in which we vary the viscosity and buoyancy structure with the observed differential paleo shoreline elevations allows us to forward model the most likely scenario for uplift along this coast.&lt;/p&gt; &lt;/div&gt;