Isotopic and Petrologic Investigation, and a Thermomechanical Model of Genesis of Large-Volume Rhyolites in Arc Environments: Karymshina Volcanic Complex, Kamchatka, Russia

The Kamchatka Peninsula of eastern Russia is currently one of the most volcanically active areas on Earth where a combination of > 8 cm/yr subduction convergence rate and thick continental crust generates large silicic magma chambers, reflected by abundant large calderas and caldera complexes. This study examines the largest center of silicic 4-0.5 Ma Karymshina Volcanic Complex, which includes the 25 × 15 km Karymshina caldera, the largest in Kamchatka. A series of rhyolitic tuff eruptions at 4 Ma were followed by the main eruption at 1.78 Ma and produced an estimated 800 km3 of rhyolitic ignimbrites followed by high-silica rhyolitic post-caldera extrusions. The postcaldera domes trace the 1.78 Ma right fracture and form a continuous compositional series with ignimbrites. We here present results of a geologic, petrologic, and isotopic study of the Karymshina eruptive complex, and present new Ar-Ar ages, and isotopic values of rocks for the oldest pre- 1.78 Ma caldera ignimbrites and intrusions, which include a diversity of compositions from basalts to rhyolites. Temporal trends in δ18O, 87Sr/86Sr, and 144Nd/143Nd indicate values comparable to neighboring volcanoes, increase in homogeneity, and temporal increase in mantle-derived Sr and Nd with increasing differentiation over the last 4 million years. Data are consistent with a batholithic scale magma chamber formed by primarily fractional crystallization of mantle derived composition and assimilation of Cretaceous and younger crust, driven by basaltic volcanism and mantle delaminations. All rocks have 35–45% quartz, plagioclase, biotite, and amphibole phenocrysts. Rhyolite-MELTS crystallization models favor shallow (2 kbar) differentiation conditions and varying quantities of assimilated amphibolite partial melt and hydrothermally-altered silicic rock. Thermomechanical modeling with a typical 0.001 km3/yr eruption rate of hydrous basalt into a 38 km Kamchatkan arc crust produces two magma bodies, one near the Moho and the other engulfing the entire section of upper crust. Rising basalts are trapped in the lower portion of an upper crustal magma body, which exists in a partially molten to solid state. Differentiation products of basalt periodically mix with the resident magma diluting its crustal isotopic signatures. At the end of the magmatism crust is thickened by 8 km. Thermomechanical modeling show that the most likely way to generate large spikes of rhyolitic magmatism is through delamination of cumulates and mantle lithosphere after many millions of years of crustal thickening. The paper also presents a chemical dataset for Pacific ashes from ODDP 882 and 883 and compares them to Karymshina ignimbrites and two other Pleistocene calderas studied by us in earlier works

Initiation of a proto‐transform fault prior to seafloor spreading

Transform faults are a fundamental tenet of plate tectonics, connecting offset extensional segments of mid‐ocean ridges in ocean basins worldwide. The current consensus is that oceanic transform faults initiate after the onset of seafloor spreading. However, this inference has been difficult to test given the lack of direct observations of transform fault formation. Here, we integrate evidence from surface faults, geodetic measurements, local seismicity, and numerical modelling of the subaerial Afar continental rift and show that a proto‐transform fault is initiating during the final stages of continental breakup. This is the first direct observation of proto‐transform fault initiation in a continental rift, and sheds unprecedented light on their formation mechanisms. We demonstrate that they can initiate during late‐stage continental rifting, earlier in the rifting cycle than previously thought. Future studies of volcanic rifted margins cannot assume that oceanic transform faults initiated after the onset of seafloor spreading

Slab detachment in laterally varying subduction zones: 3-D numerical modeling

Understanding the three-dimensional (3-D) dynamics of subduction-collision systems is a longstanding challenge in geodynamics. We investigate the impact of slab detachment in collision systems that are subjected to along-trench variations. High-resolution thermomechanical numerical models, encompassing experimentally derived flow laws and a pseudo free surface, are employed to unravel lithospheric and topographic evolutions. First, we consider coeval subduction of adjacent continental and oceanic lithospheres (SCO). This configuration yields to two-stage slab detachment during collision, topographic buildup and extrusion, variable along-trench convergence rates, and associated trench deformation. The second setting considers a convergent margin, which is laterally limited by a transform boundary (STB). Such collisional system is affected by a single slab detachment, little trench deformation, and moderately confined upper plate topography. The effect of initial thermal slab age on SCO and STB models are explored. Similarities with natural analogs along the Arabia-Eurasia collision are discussed

A Secondary Zone of Uplift Due to Megathrust Earthquakes

The 1960 M9.5 Valdivia and 1964 M9.2 Alaska earthquakes caused a decimeters high secondary zone of uplift a few hundred kilometers landward of the trench. We analyze GPS data from the 2010 M8.8 Maule and 2011 M9.0 Tohoku-Oki earthquakes to reveal the persistent existence of a secondary zone of uplift due to great earthquakes at the megathrust interface. This uplift varies in magnitude and location, but consistently occurs at a few hundred kilometers landward from the trench and is likely mainly coseismic in nature. This secondary zone of uplift is systematically predicted by our 2D visco-elasto-plastic seismothermo-mechanical numerical simulations, which model both geodynamic and seismic cycle timescales. Through testing hypotheses in both simple and realistic setups, we propose that a superposition of two physical mechanisms could be responsible for this phenomenon. First, a wavelength is introduced through elastic buckling of a visco-elastically layered fore-arc that is horizontally compressed in the interseismic period. The consequent secondary zone of interseismic subsidence is elastically rebound during the earthquake into a secondary zone of relative uplift. Second, absolute and broader uplift is ensured through a mass conservationdriven return flow following accelerated slab penetration due to the megathrust earthquake. The dip and width of the seismogenic zone and resulting (deep) coseismic slip seem to have the largest affect on location and amplitude of the secondary zone of uplift. These results imply that subduction and mantle flow do not occur at constant rates, but are rather modulated by earthquakes. This suggests a link between deep mantle and shallow surface displacements even at time scales of minutes

Devolatilization of extrasolar planetesimals by 60Fe and 26Al heating

Whilst the formation of Solar system planets is constrained by meteoritic evidence, the geophysical history of low-mass exoplanets is much less clear. The bulk composition and climate states of rocky exoplanets may vary significantly based on the composition and properties of the planetesimals they form. An important factor influenced by planetesimal composition is water content, where the desiccation of accreting planetesimals impacts the final water content of the resultant planets. While the inner planets of the Solar system are comparatively water-poor, recent observational evidence from exoplanet bulk densities and planetary formation models suggests that rocky exoplanets engulfed by substantial layers of high-pressure ices or massive steam atmospheres could be widespread. Here, we quantify variations in planetesimal desiccation due to potential fractionation of the two short-lived radioisotopes 26Al and 60Fe relevant for internal heating on planetary formation time-scales. We focus on how order of magnitude variations in 60Fe can affect the water content of planetesimals, and how this may alter the formation of extrasolar ocean worlds. We find that heating by 26Al is the dominant cause of planetesimal heating in any Solar system analogue scenario, thus validating previous works focussing only on this radioisotope. However, 60Fe can become the primary heating source in the case of high levels of supernova enrichment in massive star-forming regions. These diverging scenarios can affect the formation pathways, bulk volatile budget, and climate diversity of low-mass exoplanets

Speed of thermal adaptation of terrestrial vegetation alters Earth’s long-term climate

Earth’s long-term climate is driven by the cycling of carbon between geologic reservoirs and the atmosphere-ocean system. Our understanding of carbon-climate regulation remains incomplete, with large discrepancies remaining between biogeochemical model predictions and the geologic record. Here, we evaluate the importance of the continuous biological climate adaptation of vegetation as a regulation mechanism in the geologic carbon cycle since the establishment of forest ecosystems. Using a model, we show that the vegetation’s speed of adaptation to temperature changes through eco-evolutionary processes can strongly influence global rates of organic carbon burial and silicate weathering. Considering a limited thermal adaptation capacity of the vegetation results in a closer balance of reconstructed carbon fluxes into and out of the atmosphere-ocean system, which is a prerequisite to maintain habitable conditions on Earth’s surface on a multimillion-year timescale. We conclude that the long-term carbon-climate system is more sensitive to biological dynamics than previously expected, which may help to explain large shifts in Phanerozoic climate

A Secondary Zone of Uplift Due to Megathrust Earthquakes

The 1960 M9.5 Valdivia and 1964 M9.2 Alaska earthquakes caused a decimeters high secondary zone of uplift a few hundred kilometers landward of the trench. We analyze GPS data from the 2010 M8.8 Maule and 2011 M9.0 Tohoku-Oki earthquakes to reveal the persistent existence of a secondary zone of uplift due to great earthquakes at the megathrust interface. This uplift varies in magnitude and location, but consistently occurs at a few hundred kilometers landward from the trench and is likely mainly coseismic in nature. This secondary zone of uplift is systematically predicted by our 2D visco-elasto-plastic seismothermo-mechanical numerical simulations, which model both geodynamic and seismic cycle timescales. Through testing hypotheses in both simple and realistic setups, we propose that a superposition of two physical mechanisms could be responsible for this phenomenon. First, a wavelength is introduced through elastic buckling of a visco-elastically layered fore-arc that is horizontally compressed in the interseismic period. The consequent secondary zone of interseismic subsidence is elastically rebound during the earthquake into a secondary zone of relative uplift. Second, absolute and broader uplift is ensured through a mass conservationdriven return flow following accelerated slab penetration due to the megathrust earthquake. The dip and width of the seismogenic zone and resulting (deep) coseismic slip seem to have the largest affect on location and amplitude of the secondary zone of uplift. These results imply that subduction and mantle flow do not occur at constant rates, but are rather modulated by earthquakes. This suggests a link between deep mantle and shallow surface displacements even at time scales of minutes

Seismic behaviour of mountain belts controlled by plate convergence rate

The relative contribution of tectonic and kinematic processes to seismic behaviour of mountain belts is still controversial. To understand the partitioning between these processes we developed a model that simulates both tectonic and seismic processes in a continental collision setting. These 2D seismo-thermo-mechanical (STM) models obtain a Gutenberg–Richter frequency–magnitude distribution due to spontaneous events occurring throughout the orogen. Our simulations suggest that both the corresponding slope (b value) and maximum earthquake magnitude () correlate linearly with plate convergence rate. By analyzing 1D rheological profiles and isotherm depths we demonstrate that plate convergence rate controls the brittle strength through a rheological feedback with temperature and strain rate. Faster convergence leads to cooler temperatures and also results in more larger seismogenic domains, thereby increasing both and the relative number of large earthquakes (decreasing b value). This mechanism also predicts a more seismogenic lower crust, which is confirmed by a transition from uni- to bi-modal hypocentre depth distributions in our models. This transition and a linear relation between convergence rate and b value and is supported by our comparison of earthquakes recorded across the Alps, Apennines, Zagros and Himalaya. These results imply that deformation in the Alps occurs in a more ductile manner compared to the Himalayas, thereby reducing its seismic hazard. Furthermore, a second set of experiments with higher temperature and different orogenic architecture shows the same linear relation with convergence rate, suggesting that large-scale tectonic structure plays a subordinate role. We thus propose that plate convergence rate, which also controls the average differential stress of the orogen and its linear relation to the b value, is the first-order parameter controlling seismic hazard of mountain belts