Search

Here we define and report the relationship between the maximum seismic magnitude (M) and injection volume (ΔV) through fluid-injection fault-reactivation experiments and analysis. This relationship incorporates the in situ shear modulus (G) and fault pre-stress as a fraction of the strength drop (c), expressed as M = c/(1-c) GΔV. Injection response defines a sigmoidal relation in M - Δ V space with unit gradient limbs linked by an intermediate up-step. Both laboratory observations and analysis for a rigid fault with slip limited to the zone of pressurization show trajectories of cumulative M - Δ V that evolve at a gradient of unity, are offset in order of increasing pre-stress and are capable of step changes in moment with shear reactivation at elevated critical-stresses - key features apparent in field observations. The model and confirmatory laboratory observations explain the occurrence of some triggered earthquakes at EGS sites significantly larger than expected relative to injection volumes and based on previous models., (© 2024. The Author(s).)

Decades of studies at divergent plate margins have revealed networks of magmatic sills at the crust-mantle boundary. However, a lack of direct observations of deep magma motion limits our understanding of magma inflow from the mantle into the lower crust and the mechanism of sill formation. Here, satellite geodesy reveals rift-scale deformation caused by magma inflow in the deep crust in the Afar rift (East Africa). Simultaneous inflation of four sills, laterally separated by 10s of km and at depths ranging 9-28 km, caused uplift across a ~ 100-km-wide zone, suggesting the sills are linked to a common mantle source. Our results show the supply of magma into the lower crust is temporally episodic, occurring across a network of sills. This process reflects inherent instability of melt migration through porous mantle flow and may be the fundamental process that builds the thick igneous crust beneath magmatic rifts and rifted margins globally., (© 2024. The Author(s).)

Conflicting results remain on the impacts of climate change on marine organisms, hindering our capacity to predict the future state of marine ecosystems. To account for species-specific responses and for the ambiguous relation of most metrics to fitness, we develop a meta-analytical approach based on the deviation of responses from reference values (absolute change) to complement meta-analyses of directional (relative) changes in responses. Using this approach, we evaluate responses of fish and invertebrates to warming and acidification. We find that climate drivers induce directional changes in calcification, survival, and metabolism, and significant deviations in twice as many biological responses, including physiology, reproduction, behavior, and development. Widespread deviations of responses are detected even under moderate intensity levels of warming and acidification, while directional changes are mostly limited to more severe intensity levels. Because such deviations may result in ecological shifts impacting ecosystem structures and processes, our results suggest that climate change will likely have stronger impacts than those previously predicted based on directional changes alone., (© 2024. The Author(s).)

We link changes in crustal permeability to informative features of microearthquakes (MEQs) using two field hydraulic stimulation experiments where both MEQs and permeability evolution are recorded simultaneously. The Bidirectional Long Short-Term Memory (Bi-LSTM) model effectively predicts permeability evolution and ultimate permeability increase. Our findings confirm the form of key features linking the MEQs to permeability, offering mechanistically consistent interpretations of this association. Transfer learning correctly predicts permeability evolution of one experiment from a model trained on an alternate dataset and locale, which further reinforces the innate interdependency of permeability-to-seismicity. Models representing permeability evolution on reactivated fractures in both shear and tension suggest scaling relationships in which changes in permeability ( Δ k ) are linearly related to the seismic moment ( M ) of individual MEQs as Δ k ∝ M . This scaling relation rationalizes our observation of the permeability-to-seismicity linkage, contributes to its predictive robustness and accentuates its potential in characterizing crustal permeability evolution using MEQs., (© 2024. The Author(s).)

Volcano sector collapse and pyroclastic density currents are common phenomena on active volcanoes and potentially a fatal source of tsunami waves which constitute a serious hazard for local as well as distant coastal population. Several examples in recent history, warn us on the urgent need to improve our mitigation counter-actions when tsunamis have volcanic origin. However, instrumental record of tsunami generated by mass movement along a volcano flank are still rare and not well understood yet. Small tsunamis (≤1 m) induced by pyroclastic density currents associated to violent explosions of Stromboli volcano were recorded in near-source conditions (<1.6 km). We show how tsunami waveform remains unaltered regardless of the two orders of variability in the landslide volume and dynamics. This unprecedented record is also providing the lesson to develop unconventional warning strategies necessary when the tsunamigenic source is expected to be very close (<10 minutes) to densely populated coasts and with a limited time to issue an alert based on simulation of wave propagation and inundation., (© 2024. The Author(s).)

Metamorphic fluids, faults, and shear zones are carriers of carbon from the deep Earth to shallower reservoirs. Some of these fluids are reduced and transport energy sources, like H 2 and light hydrocarbons. Mechanisms and pathways capable of transporting these deep energy sources towards shallower reservoirs remain unidentified. Here we present geological evidence of failure of mechanically strong rocks due to the accumulation of CH 4 -H 2 -rich fluids at deep forearc depths, which ultimately reached supralithostatic pore fluid pressure. These fluids originated from adjacent reduction of carbonates by H 2 -rich fluids during serpentinization at eclogite-to-blueschist-facies conditions. Thermodynamic modeling predicts that the production and accumulation of CH 4 -H 2 -rich aqueous fluids can produce fluid overpressure more easily than carbon-poor and CO 2 -rich aqueous fluids. This study provides evidence for the migration of deep Earth energy sources along tectonic discontinuities, and suggests causal relationships with brittle failure of hard rock types that may trigger seismic activity at forearc depths., (© 2024. The Author(s).)

Predicting failure in solids has broad applications including earthquake prediction which remains an unattainable goal. However, recent machine learning work shows that laboratory earthquakes can be predicted using micro-failure events and temporal evolution of fault zone elastic properties. Remarkably, these results come from purely data-driven models trained with large datasets. Such data are equivalent to centuries of fault motion rendering application to tectonic faulting unclear. In addition, the underlying physics of such predictions is poorly understood. Here, we address scalability using a novel Physics-Informed Neural Network (PINN). Our model encodes fault physics in the deep learning loss function using time-lapse ultrasonic data. PINN models outperform data-driven models and significantly improve transfer learning for small training datasets and conditions outside those used in training. Our work suggests that PINN offers a promising path for machine learning-based failure prediction and, ultimately for improving our understanding of earthquake physics and prediction., (© 2023. The Author(s).)

Lava domes exhibit highly unpredictable and hazardous behavior, which is why imaging their morphological evolution to decipher the underlying governing mechanisms remains a major challenge. Using high-resolution satellite radar imagery enhanced with deep-learning, we image the repetitive dome construction-subsidence cycles at Popocatépetl volcano (Mexico) with very high temporal and spatial resolution. We show that these cycles resemble gas-driven rise and fall of the upper magma column, where buoyant bubble-rich magma is extruded from the conduit (in ~hours-days), and successively drained back (in ~days-months) as magma degasses and crystallizes. These cycles are superimposed on a progressive decadal crater deepening, accompanied by heat and gas flux decrease, which could be partially explained by gas depletion within the magma plumbing system. Results reinforce the idea that gas retention and escape from the magma column play a key role in the short- and long-term morphological evolution of low-viscosity lava domes and their associated hazards., (© 2023. The Author(s).)

The gut is continuously invaded by diverse bacteria from the diet and the environment, yet microbiome composition is relatively stable over time for host species ranging from mammals to insects, suggesting host-specific factors may selectively maintain key species of bacteria. To investigate host specificity, we used gnotobiotic Drosophila, microbial pulse-chase protocols, and microscopy to investigate the stability of different strains of bacteria in the fly gut. We show that a host-constructed physical niche in the foregut selectively binds bacteria with strain-level specificity, stabilizing their colonization. Primary colonizers saturate the niche and exclude secondary colonizers of the same strain, but initial colonization by Lactobacillus species physically remodels the niche through production of a glycan-rich secretion to favor secondary colonization by unrelated commensals in the Acetobacter genus. Our results provide a mechanistic framework for understanding the establishment and stability of a multi-species intestinal microbiome., (© 2023. The Author(s).)

Planktonic calcifying organisms play a key role in regulating ocean carbonate chemistry and atmospheric CO 2 . Surprisingly, references to the absolute and relative contribution of these organisms to calcium carbonate production are lacking. Here we report quantification of pelagic calcium carbonate production in the North Pacific, providing new insights on the contribution of the three main planktonic calcifying groups. Our results show that coccolithophores dominate the living calcium carbonate (CaCO 3 ) standing stock, with coccolithophore calcite comprising ~90% of total CaCO 3 production, and pteropods and foraminifera playing a secondary role. We show that pelagic CaCO 3 production is higher than the sinking flux of CaCO 3 at 150 and 200 m at ocean stations ALOHA and PAPA, implying that a large portion of pelagic calcium carbonate is remineralised within the photic zone; this extensive shallow dissolution explains the apparent discrepancy between previous estimates of CaCO 3 production derived from satellite observations/biogeochemical modeling versus estimates from shallow sediment traps. We suggest future changes in the CaCO 3 cycle and its impact on atmospheric CO 2 will largely depend on how the poorly-understood processes that determine whether CaCO 3 is remineralised in the photic zone or exported to depth respond to anthropogenic warming and acidification., (© 2023. The Author(s).)

The occurrence of rhyolite melts in the mantle has been predicted by high pressure-high temperature experiments but never observed in nature. Here we report natural quartz-bearing rhyolitic melt inclusions and interstitial glass within peridotite xenoliths. The oxygen isotope composition of quartz crystals shows the unequivocal continental crustal derivation of these melts, which approximate the minimum composition in the quartz-albite-orthoclase system. Thermodynamic modelling suggests rhyolite was originated from partial melting of near-anhydrous garnet-bearing metapelites at temperatures ~1000 °C and interacted with peridotite at pressure ~1 GPa. Reaction of rhyolite with olivine converted lherzolite rocks into orthopyroxene-domains and orthopyroxene + plagioclase veins. The recognition of rhyolitic melts in the mantle provides direct evidence for element cycling through earth's reservoirs, accommodated by dehydration and melting of crustal material, brought into the mantle by subduction, chemically modifying the mantle source, and ultimately returning to surface by arc magmatism., (© 2022. The Author(s).)

Surface faulting earthquakes are known to cluster in time from historical and palaeoseismic studies, but the mechanism(s) responsible for clustering, such as fault interaction, strain-storage, and evolving dynamic topography, are poorly quantified, and hence not well understood. We present a quantified replication of observed earthquake clustering in central Italy. Six active normal faults are studied using 36 Cl cosmogenic dating, revealing out-of-phase periods of high or low surface slip-rate on neighboring structures that we interpret as earthquake clusters and anticlusters. Our calculations link stress transfer caused by slip averaged over clusters and anti-clusters on coupled fault/shear-zone structures to viscous flow laws. We show that (1) differential stress fluctuates during fault/shear-zone interactions, and (2) these fluctuations are of sufficient magnitude to produce changes in strain-rate on viscous shear zones that explain slip-rate changes on their overlying brittle faults. These results suggest that fault/shear-zone interactions are a plausible explanation for clustering, opening the path towards process-led seismic hazard assessments., (© 2022. The Author(s).)

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters., (© 2022. The Author(s).)

Analysis of seismicity can illuminate active fault zone structures but also deformation within large volumes of the seismogenic zone. For the M w 6.5 2016-2017 Central Italy seismic sequence, seismicity not only localizes along the major structures hosting the mainshocks (on-fault seismicity), but also occurs within volumes of Triassic Evaporites, TE, composed of alternated anhydrites and dolostones. These volumes of distributed microseismicity show a different frequency-magnitude distribution than on-fault seismicity. We interpret that, during the sequence, shear strain-rate increase, and fluid overpressure promoted widespread ductile deformation within TE that light-up with distributed microseismicity. This interpretation is supported by field and laboratory observations showing that TE background ductile deformation is complex and dominated by distributed failure and folding of the anhydrites associated with boudinage hydro-fracturing and faulting of dolostones. Our results indicate that ductile crustal deformation can cause distributed microseismicity, which obeys to different scaling laws than on-fault seismicity occurring on structures characterized by elasto-frictional stick-slip behaviour., (© 2022. The Author(s).)

The Magnitude-Frequency-Distribution (MFD) of earthquakes is typically modeled with the (tapered) Gutenberg-Richter relation. The main parameter of this relation, the b-value, controls the relative rate of small and large earthquakes. Resolving spatiotemporal variations of the b-value is critical to understanding the earthquake occurrence process and improving earthquake forecasting. However, this variation is not well understood. Here we present remarkable MFD variability during the complex 2016/17 central Italy sequence using a high-resolution earthquake catalog. Isolating seismically active volumes ('clusters') reveals that the MFD differed in nearby clusters, varied or remained constant in time depending on the cluster, and increased in b-value in the cluster where the largest earthquake eventually occurred. These findings suggest that the fault system's heterogeneity and complexity influence the MFD. Our findings raise the question "b-value of what?": interpreting and using MFD variability needs a spatiotemporal scale that is physically meaningful, like the one proposed here., (© 2022. The Author(s).)

Ocean sediments consist mainly of calcium carbonate and organic matter (phytoplankton debris). Once subducted, some carbon is removed from the slab and returns to the atmosphere as CO 2 in arc magmas. Its isotopic signature is thought to reflect the bulk fraction of inorganic (carbonate) and organic (graphitic) carbon in the sedimentary source. Here we challenge this assumption by experimentally investigating model sediments composed of 13 C-CaCO 3  +  12 C-graphite interacting with water at pressure, temperature and redox conditions of an average slab-mantle interface beneath arcs. We show that oxidative dissolution of graphite is the main process controlling the production of CO 2 , and its isotopic composition reflects the CO 2 /CaCO 3 rather than the bulk graphite/CaCO 3 (i.e., organic/inorganic carbon) fraction. We provide a mathematical model to relate the arc CO 2 isotopic signature with the fluid-rock ratios and the redox state in force in its subarc source., (© 2022. The Author(s).)

During explosive eruption of low viscosity magmas, pyroclasts are cooled predominantly by forced convection. Depending on the cooling efficiency relative to other timescales, a spectrum of deposits can be formed. Deposition of hot clasts, above their glass transition temperature, can form spatter mounds, ramparts and clastogenic lava flows. Clasts may also be deposited cold, producing tephra cones and blankets. Thus, the deposit and pyroclast type can provide information about eruption dynamics and magma properties. Here we examine pyroclasts from Tseax volcano, British Columbia, Canada. These newly identified inflated pyroclasts, are fluidal in form, have undergone post-depositional expansion, and are found juxtaposed with scoria. Detailed field, chemical and textural observations, coupled with high temperature rheometry and thermal modelling, reveal that abrupt transitions in eruptive behaviour - from lava fountaining to low-energy bubble bursts - created these pyroclastic deposits. These findings should help identify transitions in eruptive behaviour at other mafic volcanoes worldwide., (© 2022. The Author(s).)

Assessing the timing of great megathrust earthquakes is together crucial for seismic hazard analysis and deemed impossible. Geodetic instrumentation of subduction zones has revealed unexpected deformation patterns at subduction segments adjacent to those that hosted recent mega-earthquakes: coastal sites move landward with faster velocities than before the earthquake. Here, we show observations from the largest and best-monitored megathrust earthquakes, and from a scaled analog model, to reveal that these events create coseismic and postseismic deformation patterns typical of a complete gear-like rotation about a vertical axis, hereafter called twisting. We find that such twisting alters the interseismic velocity field of adjacent subduction segments depending on the time since the last earthquake. Early interactions accelerate while late interactions decelerate local kinematics. This finding opens the possibility of using megathrust earthquakes, the characteristics of the twisting pattern, and the ensuing geodetic velocity changes, as a proxy for estimating the timing of the seismic cycle at unruptured segments along the margin., (© 2022. The Author(s).)

Fractionation effects related to evaporation and condensation had a major impact on the current elemental and isotopic composition of the Solar System. Although isotopic fractionation of moderately volatile elements has been observed in tektites due to impact heating, the exact nature of the processes taking place during hypervelocity impacts remains poorly understood. By studying Fe in microtektites, here we show that impact events do not simply lead to melting, melt expulsion and evaporation, but involve a convoluted sequence of processes including condensation, variable degrees of mixing between isotopically distinct reservoirs and ablative evaporation during atmospheric re-entry. Hypervelocity impacts can as such not only generate isotopically heavy, but also isotopically light ejecta, with δ 56/54 Fe spanning over nearly 5‰ and likely even larger variations for more volatile elements. The mechanisms demonstrated here for terrestrial impact ejecta modify our understanding of the effects of impact processing on the isotopic evolution of planetary crusts., (© 2021. The Author(s).)

The mild activity of basaltic volcanoes is punctuated by violent explosive eruptions that occur without obvious precursors. Modelling the source processes of these sudden blasts is challenging. Here, we use two decades of ground deformation (tilt) records from Stromboli volcano to shed light, with unprecedented detail, on the short-term (minute-scale) conduit processes that drive such violent volcanic eruptions. We find that explosive eruptions, with source parameters spanning seven orders of magnitude, all share a common pre-blast ground inflation trend. We explain this exponential inflation using a model in which pressure build-up is caused by the rapid expansion of volatile-rich magma rising from depth into a shallow (<400 m) resident magma conduit. We show that the duration and amplitude of this inflation trend scales with the eruption magnitude, indicating that the explosive dynamics obey the same (scale-invariant) conduit process. This scale-invariance of pre-explosion ground deformation may usher in a new era of short-term eruption forecasting.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Geological sources of H 2 and abiotic CH 4 have had a critical role in the evolution of our planet and the development of life and sustainability of the deep subsurface biosphere. Yet the origins of these sources are largely unconstrained. Hydration of mantle rocks, or serpentinization, is widely recognized to produce H 2 and favour the abiotic genesis of CH 4 in shallow settings. However, deeper sources of H 2 and abiotic CH 4 are missing from current models, which mainly invoke more oxidized fluids at convergent margins. Here we combine data from exhumed subduction zone high-pressure rocks and thermodynamic modelling to show that deep serpentinization (40-80 km) generates significant amounts of H 2 and abiotic CH 4 , as well as H 2 S and NH 3 . Our results suggest that subduction, worldwide, hosts large sources of deep H 2 and abiotic CH 4 , potentially providing energy to the overlying subsurface biosphere in the forearc regions of convergent margins.

Faulting and earthquakes occur extensively along the flanks of the East African Rift System, including an offshore branch in the western Indian Ocean, resulting in remobilization of sediment in the form of landslides. To date, constraints on the occurrence of submarine landslides at margin scale are lacking, leaving unanswered a link between rifting and slope instability. Here, we show the first overview of landslide deposits in the post-Eocene stratigraphy of the Tanzania margin and we present the discovery of one of the biggest landslides on Earth: the Mafia mega-slide. The emplacement of multiple landslides, including the Mafia mega-slide, during the early-mid Miocene is coeval with cratonic rifting in Tanzania, indicating that plateau uplift and rifting in East Africa triggered large and potentially tsunamigenic landslides likely through earthquake activity and enhanced sediment supply. This study is a first step to evaluate the risk associated with submarine landslides in the region.

Flank instability and sector collapses, which pose major threats, are common on volcanic islands. On 22 Dec 2018, a sector collapse event occurred at Anak Krakatau volcano in the Sunda Strait, triggering a deadly tsunami. Here we use multiparametric ground-based and space-borne data to show that prior to its collapse, the volcano exhibited an elevated state of activity, including precursory thermal anomalies, an increase in the island's surface area, and a gradual seaward motion of its southwestern flank on a dipping décollement. Two minutes after a small earthquake, seismic signals characterize the collapse of the volcano's flank at 13:55 UTC. This sector collapse decapitated the cone-shaped edifice and triggered a tsunami that caused 430 fatalities. We discuss the nature of the precursor processes underpinning the collapse that culminated in a complex hazard cascade with important implications for the early detection of potential flank instability at other volcanoes.

Continental rift systems form by propagation of isolated rift segments that interact, and eventually evolve into continuous zones of deformation. This process impacts many aspects of rifting including rift morphology at breakup, and eventual ocean-ridge segmentation. Yet, rift segment growth and interaction remain enigmatic. Here we present geological data from the poorly documented Ririba rift (South Ethiopia) that reveals how two major sectors of the East African rift, the Kenyan and Ethiopian rifts, interact. We show that the Ririba rift formed from the southward propagation of the Ethiopian rift during the Pliocene but this propagation was short-lived and aborted close to the Pliocene-Pleistocene boundary. Seismicity data support the abandonment of laterally offset, overlapping tips of the Ethiopian and Kenyan rifts. Integration with new numerical models indicates that rift abandonment resulted from progressive focusing of the tectonic and magmatic activity into an oblique, throughgoing rift zone of near pure extension directly connecting the rift sectors.

Oxygen-isotope thermometry played a critical role in the rise of modern geochemistry and remains extensively used in (bio-)geoscience. Its theoretical foundations rest on the assumption that 18 O/ 16 O partitioning among water and carbonate minerals primarily reflects thermodynamic equilibrium. However, after decades of research, there is no consensus on the true equilibrium 18 O/ 16 O fractionation between calcite and water ( 18 α cc/w ). Here, we constrain the equilibrium relations linking temperature, 18 α cc/w , and clumped isotopes (Δ 47 ) based on the composition of extremely slow-growing calcites from Devils Hole and Laghetto Basso (Corchia Cave). Equilibrium 18 α cc/w values are systematically ~1.5‰ greater than those in biogenic and synthetic calcite traditionally considered to approach oxygen-isotope equilibrium. We further demonstrate that subtle disequilibria also affect Δ 47 in biogenic calcite. These observations provide evidence that most Earth-surface calcites fail to achieve isotopic equilibrium, highlighting the need to improve our quantitative understanding of non-equilibrium isotope fractionation effects instead of relying on phenomenological calibrations.

Considerable ambiguity remains over the extent and nature of millennial/centennial-scale climate instability during the Last Interglacial (LIG). Here we analyse marine and terrestrial proxies from a deep-sea sediment sequence on the Portuguese Margin and combine results with an intensively dated Italian speleothem record and climate-model experiments. The strongest expression of climate variability occurred during the transitions into and out of the LIG. Our records also document a series of multi-centennial intra-interglacial arid events in southern Europe, coherent with cold water-mass expansions in the North Atlantic. The spatial and temporal fingerprints of these changes indicate a reorganization of ocean surface circulation, consistent with low-intensity disruptions of the Atlantic meridional overturning circulation (AMOC). The amplitude of this LIG variability is greater than that observed in Holocene records. Episodic Greenland ice melt and runoff as a result of excess warmth may have contributed to AMOC weakening and increased climate instability throughout the LIG.

Myrmeleontiformia are an ancient group of lacewing insects characterized by predatory larvae with unusual morphologies and behaviours. Mostly soil dwellers with a soft cuticle, their larvae fossilize only as amber inclusions, and thus their fossil record is remarkably sparse. Here, we document a disparate assemblage of myrmeleontiform larvae from the mid-Cretaceous amber (99 Ma) of Myanmar, evidence of a considerable diversification. Our cladistic analysis integrating extant and extinct taxa resolves the fossils as both stem- and crown-groups. Similarities between extinct and extant species permit inferences of larval ethology of the fossil species through statistical correlation analyses with high support, implying that morphological disparity matched behavioural diversity. An improved understanding of the evolutionary history of antlions and relatives supports the conclusion that hunting strategies, such as camouflage and fossoriality, were acquired early within the lineage.

Dinosaurs diversified in two steps during the Triassic. They originated about 245 Ma, during the recovery from the Permian-Triassic mass extinction, and then remained insignificant until they exploded in diversity and ecological importance during the Late Triassic. Hitherto, this Late Triassic explosion was poorly constrained and poorly dated. Here we provide evidence that it followed the Carnian Pluvial Episode (CPE), dated to 234-232 Ma, a time when climates switched from arid to humid and back to arid again. Our evidence comes from a combined analysis of skeletal evidence and footprint occurrences, and especially from the exquisitely dated ichnofaunas of the Italian Dolomites. These provide evidence of tetrapod faunal compositions through the Carnian and Norian, and show that dinosaur footprints appear exactly at the time of the CPE. We argue then that dinosaurs diversified explosively in the mid Carnian, at a time of major climate and floral change and the extinction of key herbivores, which the dinosaurs opportunistically replaced.

Estimates of dissolved CO 2 in subduction-zone fluids are based on thermodynamic models, relying on a very sparse experimental data base. Here, we present experimental data at 1-3 GPa, 800 °C, and ∆FMQ ≈ -0.5 for the volatiles and solute contents of graphite-saturated fluids in the systems COH, SiO 2 -COH ( + quartz/coesite) and MgO-SiO 2 -COH ( + forsterite and enstatite). The CO 2 content of fluids interacting with silicates exceeds the amounts measured in the pure COH system by up to 30 mol%, as a consequence of a decrease in water activity probably associated with the formation of organic complexes containing Si-O-C and Si-O-Mg bonds. The interaction of deep aqueous fluids with silicates is a novel mechanism for controlling the composition of subduction COH fluids, promoting the deep CO 2 transfer from the slab-mantle interface to the overlying mantle wedge, in particular where fluids are stable over melts.Current estimates of dissolved CO 2 in subduction-zone fluids based on thermodynamic models rely on a very sparse experimental data base. Here, the authors show that experimental graphite-saturated COH fluids interacting with silicates at 1-3 GPa and 800 °C display unpredictably high CO 2 contents.

The presence of carbonates in inclusions in diamonds coming from depths exceeding 670 km are obvious evidence that carbonates exist in the Earth's lower mantle. However, their range of stability, crystal structures and the thermodynamic conditions of the decarbonation processes remain poorly constrained. Here we investigate the behaviour of pure iron carbonate at pressures over 100 GPa and temperatures over 2,500 K using single-crystal X-ray diffraction and Mössbauer spectroscopy in laser-heated diamond anvil cells. On heating to temperatures of the Earth's geotherm at pressures to ∼50 GPa FeCO 3 partially dissociates to form various iron oxides. At higher pressures FeCO 3 forms two new structures-tetrairon(III) orthocarbonate Fe 4 3+ C 3 O 12 , and diiron(II) diiron(III) tetracarbonate Fe 2 2+ Fe 2 3+ C 4 O 13 , both phases containing CO 4 tetrahedra. Fe 4 C 4 O 13 is stable at conditions along the entire geotherm to depths of at least 2,500 km, thus demonstrating that self-oxidation-reduction reactions can preserve carbonates in the Earth's lower mantle.

Molten I-type cosmic spherules formed by heating, oxidation and melting of extraterrestrial Fe,Ni metal alloys. The entire oxygen in these spherules sources from the atmosphere. Therefore, I-type cosmic spherules are suitable tracers for the isotopic composition of the upper atmosphere at altitudes between 80 and 115 km. Here we present data on I-type cosmic spherules collected in Antarctica. Their composition is compared with the composition of tropospheric O 2 . Our data suggest that the Earth's atmospheric O 2 is isotopically homogenous up to the thermosphere. This makes fossil I-type micrometeorites ideal proxies for ancient atmospheric CO 2 levels.

Alteration of ultramafic rocks plays a major role in the production of hydrocarbons and organic compounds via abiotic processes on Earth and beyond and contributes to the redistribution of C between solid and fluid reservoirs over geological cycles. Abiotic methanogenesis in ultramafic rocks is well documented at shallow conditions, whereas natural evidence at greater depths is scarce. Here we provide evidence for intense high-pressure abiotic methanogenesis by reduction of subducted ophicarbonates. Protracted (≥0.5-1 Ma), probably episodic infiltration of reduced fluids in the ophicarbonates and methanogenesis occurred from at least ∼40 km depth to ∼15-20 km depth. Textural, petrological and isotopic data indicate that methane reached saturation triggering the precipitation of graphitic C accompanied by dissolution of the precursor antigorite. Continuous infiltration of external reducing fluids caused additional methane production by interaction with the newly formed graphite. Alteration of high-pressure carbonate-bearing ultramafic rocks may represent an important source of abiotic methane, with strong implications for the mobility of deep C reservoirs.

Melting of the mantle during continental breakup leads to magmatic intrusion and volcanism, yet our understanding of the location and dominant mechanisms of melt generation in rifting environments is impeded by a paucity of direct observations of mantle melting. It is unclear when during the rifting process the segmented nature of magma supply typical of seafloor spreading initiates. Here, we use Rayleigh-wave tomography to construct a high-resolution absolute three-dimensional shear-wave velocity model of the upper 250 km beneath the Afar triple junction, imaging the mantle response during progressive continental breakup. Our model suggests melt production is highest and melting depths deepest early during continental breakup. Elevated melt production during continental rifting is likely due to localized thinning and melt focusing when the rift is narrow. In addition, we interpret segmented zones of melt supply beneath the rift, suggesting that buoyancy-driven active upwelling of the mantle initiates early during continental rifting.

Constraining the timescales of pre-eruptive magmatic processes in active volcanic systems is paramount to understand magma chamber dynamics and the triggers for volcanic eruptions. Temporal information of magmatic processes is locked within the chemical zoning profiles of crystals but can be accessed by means of elemental diffusion chronometry. Mineral compositional zoning testifies to the occurrence of substantial temperature differences within magma chambers, which often bias the estimated timescales in the case of multi-stage zoned minerals. Here we propose a new Non-Isothermal Diffusion Incremental Step model to take into account the non-isothermal nature of pre-eruptive processes, deconstructing the main core-rim diffusion profiles of multi-zoned crystals into different isothermal steps. The Non-Isothermal Diffusion Incremental Step model represents a significant improvement in the reconstruction of crystal lifetime histories. Unravelling stepwise timescales at contrasting temperatures provides a novel approach to constraining pre-eruptive magmatic processes and greatly increases our understanding of magma chamber dynamics.

Slow earthquakes represent an important conundrum in earthquake physics. While regular earthquakes are catastrophic events with rupture velocities governed by elastic wave speed, the processes that underlie slow fault slip phenomena, including recent discoveries of tremor, slow-slip and low-frequency earthquakes, are less understood. Theoretical models and sparse laboratory observations have provided insights, but the physics of slow fault rupture remain enigmatic. Here we report on laboratory observations that illuminate the mechanics of slow-slip phenomena. We show that a spectrum of slow-slip behaviours arises near the threshold between stable and unstable failure, and is governed by frictional dynamics via the interplay of fault frictional properties, effective normal stress and the elastic stiffness of the surrounding material. This generalizable frictional mechanism may act in concert with other hypothesized processes that damp dynamic ruptures, and is consistent with the broad range of geologic environments where slow earthquakes are observed.

Ice streams drain large portions of ice sheets and play a fundamental role in governing their response to atmospheric and oceanic forcing, with implications for sea-level change. The mechanisms that generate ice stream flow remain elusive. Basal sliding and/or bed deformation have been hypothesized, but ice stream beds are largely inaccessible. Here we present a comprehensive, multi-scale study of the internal structure of mega-scale glacial lineations (MSGLs) formed at the bed of a palaeo ice stream. Analyses were undertaken at macro- and microscales, using multiple techniques including X-ray tomography, thin sections and ground penetrating radar (GPR) acquisitions. Results reveal homogeneity in stratigraphy, kinematics, granulometry and petrography. The consistency of the physical and geological properties demonstrates a continuously accreting, shallow-deforming, bed and invariant basal conditions. This implies that ice stream basal motion on soft sediment beds during MSGL formation is accommodated by plastic deformation, facilitated by continuous sediment supply and an inefficient drainage system.

Effusive eruptions are explained as the mechanism by which volcanoes restore the equilibrium perturbed by magma rising in a chamber deep in the crust. Seismic, ground deformation and topographic measurements are compared with effusion rate during the 2007 Stromboli eruption, drawing an eruptive scenario that shifts our attention from the interior of the crust to the surface. The eruption is modelled as a gravity-driven drainage of magma stored in the volcanic edifice with a minor contribution of magma supplied at a steady rate from a deep reservoir. Here we show that the discharge rate can be predicted by the contraction of the volcano edifice and that the very-long-period seismicity migrates downwards, tracking the residual volume of magma in the shallow reservoir. Gravity-driven magma discharge dynamics explain the initially high discharge rates observed during eruptive crises and greatly influence our ability to predict the evolution of effusive eruptions.

The Late Cretaceous 'greenhouse' world witnessed a transition from one of the warmest climates of the past 140 million years to cooler conditions, yet still without significant continental ice. Low-latitude sea surface temperature (SST) records are a vital piece of evidence required to unravel the cause of Late Cretaceous cooling, but high-quality data remain illusive. Here, using an organic geochemical palaeothermometer (TEX86), we present a record of SSTs for the Campanian-Maastrichtian interval (~83-66 Ma) from hemipelagic sediments deposited on the western North Atlantic shelf. Our record reveals that the North Atlantic at 35 °N was relatively warm in the earliest Campanian, with maximum SSTs of ~35 °C, but experienced significant cooling (~7 °C) after this to <~28 °C during the Maastrichtian. The overall stratigraphic trend is remarkably similar to records of high-latitude SSTs and bottom-water temperatures, suggesting that the cooling pattern was global rather than regional and, therefore, driven predominantly by declining atmospheric pCO2 levels.

The discovery of a natural quasicrystal, icosahedrite (Al63Cu24Fe13), accompanied by khatyrkite (CuAl2) and cupalite (CuAl) in the CV3 carbonaceous chondrite Khatyrka has posed a mystery as to what extraterrestrial processes led to the formation and preservation of these metal alloys. Here we present a range of evidence, including the discovery of high-pressure phases never observed before in a CV3 chondrite, indicating that an impact shock generated a heterogeneous distribution of pressures and temperatures in which some portions reached at least 5 GPa and 1,200 °C. The conditions were sufficient to melt Al-Cu-bearing minerals, which then rapidly solidified into icosahedrite and other Al-Cu metal phases. The meteorite also contains heretofore unobserved phases of iron-nickel and iron sulphide with substantial amounts of Al and Cu. The presence of these phases in Khatyrka provides further proof that the Al-Cu alloys are natural products of unusual processes that occurred in the early solar system.