Search

Chalcophile element concentrations in melt inclusions and matrix glasses may be used to investigate low pressure degassing processes, as well as sulfide saturation during crustal fractionation, and mantle melting. Erupted products from Kīlauea Volcano, Hawaiʻi, record three stages of sulfide saturation (in the mantle, crust, and within lava lakes), separated by episodes of sulfide resorption (i.e., sulfide undersaturation) during ascent through the thick Hawaiian lithosphere, and during syn-eruptive degassing. The presence of residual sulfides in the mantle source throughout the melting interval accounts for the high S concentrations of Kīlauean primary melts (1387–1600 ppm). Residual sulfides retain chalcophile elements during melting, decoupling the variability of these elements in high MgO melts from that of lithophile elements. Decompression associated with magma ascent through the thick Hawaiian lithosphere drives an increase in the sulfide concentration at sulfide saturation (SCSS2-), resulting in shallow storage reservoirs (∼1–5 km depth) being supplied with sulfide-undersaturated melts. A drop in temperature, coupled with major element changes during the fractionation of olivine, causes the SCSS2- to decrease. Combined with an increase in melt S contents during fractionation, this initiates a second stage of sulfide saturation at relatively high MgO contents (∼12 wt% MgO). Syn-eruptive degassing of S drives the resorption of sulfides in contact with the carrier liquid. The covariance structure of Cu, MgO and Ni contents in melt inclusions and matrix glasses indicates that the dissolution of sulfides effectively liberates sulfide-hosted Cu and Ni back into the melt, rather than the vapour phase. The contrasting behaviour of Cu, Ni, Se and S during sulfide resorption indicates that the chalcophile element signature of the Kīlauean plume is largely controlled by silicate melt-vapour partitioning, rather than sulfide-vapour partitioning. The participation of dense sulfide liquids in shallow degassing processes may result from their direct attachment to buoyant vapour bubbles, or olivine crystals which were remobilized prior to eruption. Sulfide resorption obscures the textural and chemical record of sulfide saturation in matrix glasses, but not in melt inclusions, which are isolated from this late-stage release of chalcophile elements. The partitioning of S between the dissolving sulfide, melt and the vapour phase accounts for approximately 20% of the total S release into the atmosphere.

Intra-basinal highs within classic ‘block and basin’ style tectonic frameworks, in areas such as northern England, are underpinned by large granite bodies. This is widely believed to relate to the relative ‘rigidity’ and ‘buoyancy’ of granite in relation to accommodating basement. It has been suggested that during periods of tectonic extension, normal faulting around the peripheral regions of granite batholiths permits granite-cored blocks to resist subsidence, thus forming stable areas during periods of widespread faulting-induced subsidence. However, one-dimensional modelling indicates that relatively less dense crust is incapable of resisting subsidence in this way. Instead, when local isostasy is assumed, the occurrence of granite-cored, intra-basinal highs relates to initial isostatic compensation following granite emplacement. Differential sediment loading during extensional tectonism exaggerates this profile. An integrated two-dimensional lithospheric numerical modelling approach highlights the role of flexural rigidity in limiting the amplitude whilst increasing the wavelength of isostatic deflection. In light of these models, it is suggested that such a response leaves residual second-order stresses associated with the under-compensated buoyancy of the granite body and flexural tension. The observed basin geometries of the Carboniferous North Pennine Basin can be replicated by incorporating a density deficiency within the crust, flexural rigidity, simple shear deformation within the shallower subsurface and pure shear deformation within the deeper subsurface. In adopting this technique, the regional flexural profile in response to underlying granite bodies and large extensional faults can be reproduced and thus, to an extent, validated. It is proposed that the interaction of three factors dictate the tectonic framework within a partially granitic, brittle-ductile lithosphere and the occurrence of inter-basinal highs: 1) non-tectonic, ‘second-order’ stresses such as the flexural response of the lithosphere and residual, under-compensated buoyancy forces in relation to granite bodies; 2) extensional tectonic stress and importantly; 3) inherited basement fabric.

The segregation and accumulation of felsic melts, from crystallising crustal magma reservoirs, is essential for the chemical evolution of the crust and is a phenomenon preceding some of the largest eruptions on Earth. The physical properties of residual melt and magma and the time over which the conditions remain appropriate for melt extraction are important factors controlling the efficiency of melt extraction and the distribution of melt in magma reservoirs. Here we focus on the initial water content (H2Oi) of magma as it affects both the physical properties of the residual melt and the timescales at which conditions remain appropriate for melt extraction during progressive magma crystallisation. We use rhyolite-MELTS simulations to evaluate the physical evolution of crystallising granodioritic (or dacitic) hydrous magma (i.e. ≥1 wt.% H 2 O i ) at shallow depth at 200 MPa. To constrain the solidification timescales of reservoirs containing magmas with initially different water content, we perform 2.5D thermal modelling. We combine these results with calculations of melt extraction velocity by compaction and hindered settling to identify the optimal conditions at which melt segregation occurs. These calculations suggest that hydrous felsic magmas that attain water saturation after 40 wt.% crystallisation (rheological locking point) are best suited for melt extraction. Once water-saturation is achieved, the rate of release of latent heat of crystallisation and with it the time magma spends within a given temperature interval increases while the viscosity of the residual liquid and crystal–liquid density contrast remain favourable for melt segregation. We first test our findings on the Takidani pluton (Japan) because it shows evidence of residual melt segregation from crystallising magma, and is associated with caldera-forming eruptions. We finally generalise our results to crustal magma reservoirs containing hydrous felsic magmas. Our results suggest that if segregation starts at rheological locking (i.e. crystallinity of 40–50 wt.%) upper crustal reservoirs of ≥100 km3 granodioritic (i.e. dacitic) magma with >2 wt.% H 2 O i can produce large melt-rich caps at the top of partially crystallised magma reservoirs in few hundreds to few thousands of years. The formation of separate melt lenses becomes more likely when segregation of melt starts at crystallinities >0.6. Our results suggest that H 2 O i plays an important role in modulating the distribution of eruptible melt in upper crustal reservoirs. Reservoirs of felsic and water-poor magma ( H 2 O i ) tend to be associated with the formation of isolated pockets of crystal-poor and eruptible magma, which could account for the often-observed geochemical heterogeneity of the products of large caldera-forming eruptions in the Snake River Plane. The limited dimensions of these eruptible magma pockets make their detection by geophysical methods challenging.

To meet the Paris Agreement target of limiting global warming to 2ºC or below it is widely accepted that Carbon Capture and Storage (CCS) will have to be deployed at scale. The influence of residual trapping on CO2 well injectivity and its response over time has a major impact on the injection efficiency and storage capacity of CO2 storage sites. For the first time, experiments have been undertaken over six cycles of water and supercritical CO2 injection using a state of the art high flow rig recreating in-situ conditions of near wellbore injection into analogue storage reservoir rocks. The results show that differential pressure continuously increases over multiple injection cycles. Our interpretation is that multiple cycles of injection results in a reduced effective permeability due to increased residual trapping acting as a barrier to flow resulting in reduced injectivity rates. This is supported by numerical modelling and field observations that show CO2 injectivity and its response over time will be affected by multiple cycles of injection. These results suggest that loss of injectivity must be incorporated into the injection strategy and that careful management of cyclic injection will create the opportunity to increase residual trapping.

International audience; Sediment flux signals from source to sink in foreland basins preserve a record of tectonics, sea level and climate through erosion and sedimentation. However, longitudinal sediment transport often occurs in foreland basins, thus removing part of the orogenic material flux from foreland basin records. Here we use mass balance calculation and stratigraphic simulations of sediment fluxes for the Taiwan orogen to provide an order of magnitude estimate of how much orogenic material may bypass a foreland basin. Our results indicate a significant, potentially more than 50%, mismatch between sediment volume currently preserved in the basin and the amount of material presumably eroded from the orogen since the onset of collision in Taiwan. This suggests either a significant overestimation of average erosion rates over the period concerned with orogenic development of Taiwan, or it supports previous paleogeographic work suggesting that longitudinal sediment transport in the paleo-Taiwan Strait served as a major bypass conduit of importance for the establishment of a steady state orogen. We identify candidate submarine topography in the South China Sea that may preserve Taiwan's missing erosional mass.

In the Southern Steep Belt, Italian Central Alps, relicts of the pre-Alpine continental crust are preserved. Between Valtellina and Val Camonica, a poly-metamorphic rock association occurs, which belongs to the Austroalpine units and includes two classically subdivided units: the Languard-Campo nappe (LCN) and the Tonale Series (TS). The outcropping rocks are low to medium grade muscovite, biotite and minor staurolite-bearing gneisses and micaschists, which include interlayered garnet- and biotite-bearing amphibolites, marbles, quartzites and pegmatites, as well as sillimanite-bearing gneisses and micaschists. Permian intrusives (granitoids, diorites and minor gabbros) emplaced in the metamorphic rocks. We performed a detailed structural, petrological and geochronological analysis focusing on the two main lithotypes, namely, staurolite-bearing micaschists and sillimanite- bearing paragneisses, to reconstruct the Variscan and Permian-Triassic history of this crustal section. The reconstruction of the tectono-metamorphic evolution allows for the distinction between two different tectono-metamorphic units during the early pre-Alpine evolution (D1) and predates the Permian intrusives, which comprise rocks from both TS and LCN. In the staurolite- bearing micaschists, D1 developed under amphibolite facies conditions (P=0.7- 1.1 GPa, T=580-660◦C), while in the sillimanite-bearing paragneisses formed under granulite facies conditions (P=0.6-1.0 GPa, T>780◦C). The two tectono-metamorphic units coupled together during the second pre-Alpine stage (D2) under granulite-amphibolite facies conditions at a lower pressure (P=0.4-0.6 GPa, T=620-750◦C) forming a single tectono-metamorphic unit (Languard- Tonale Tectono-Metamorphic Unit), which comprised the previously distinguished LCN and TS. Geochronological analyses on zircon rims indicate ages ranging between 250 and 275 Ma for D2, contemporaneous with the emplacement of Permian intrusives. This event developed under a high thermal state, which is compatible with an extensional tectonic setting that occurred during the exhumation of the Languard-Tonale Tectono-Metamorphic Unit. The extensional regime is interpreted as being responsible for the thinning of the Adriatic continental lithosphere during the Permian, which may be related to an early rifting phase of Pangea.

Many rifts develop in response to multiphase extension with numerical and physical models suggesting that reactivation of first-phase normal faults and rift-related variations in bulk crustal rheology control the evolution and final geometry of subsequent rifts. However, many natural multiphase rifts are deeply buried and thus poorly exposed in the field and poorly imaged in seismic reflection data, making it difficult to test these models. Here we integrate recent 3D seismic reflection and borehole data across the entire East Shetland Basin, northern North Sea, to constrain the long-term, regional development of this multiphase rift. We document the following key stages of basin development: (i) pre-Triassic to earliest Triassic development of multiple sub-basins controlled by widely distributed, NNW- to NE-trending, east- and west-dipping faults; (ii) Triassic activity on a single major, NE-trending, west-dipping fault located near the basins western margin, and formation a large half-graben; and (iii) Jurassic development of a large, E-dipping, N- to NE-trending half-graben near the eastern margin of the basin, which was associated with rift narrowing and strain focusing in the Viking Graben. In contrast to previous studies, which argue for two discrete periods of rifting during the Permian-Triassic and Late Jurassic-Early Cretaceous, we find that rifting in the East Shetland Basin was protracted from pre-Triassic to Cretaceous. We find that, during the Jurassic, most pre-Jurassic normal faults were buried and in some cases cross-cut by newly formed faults, with only a few being reactivated. Previously developed faults thus had only a limited control on the evolution and geometry of the later rift. We instead argue that strain migration and rift narrowing was linked to the evolving thermal state of the lithosphere, an interpretation supporting the predictions of lithosphere scale numerical models. Our study indicates that additional regional studies of natural rifts are required to test and refine the predictions of physical and numerical models, more specifically, our study suggests models not explicitly recognising or including thermal or rheological effects might over emphasise the role of discrete pre-existing rift structures such as normal faults.

Carbonate rocks are important archives of past ocean conditions as well as hosts of economic resources such as hydrocarbons, water, and minerals. Geologists typically perform compositional analysis of grain, matrix, cement and pore types in order to interpret depositional environments, diagenetic modification, and reservoir quality of carbonate strata. Such information can be obtained primarily from petrographic analysis, a task that is costly, labor-intensive, and requires in-depth knowledge of carbonate petrology and micropaleontology. Recent studies have leveraged machine learning-based image analysis, including Deep Convolutional Neural Networks (DCNN), to automate description, classification and interpretation of thin sections, subsurface core images and seismic facies, which would accelerate data acquisition and reproducibility for these tasks. In carbonate rocks, this approach has been applied primarily to recognize carbonate lithofacies, and no attempt has been made to individually identify and quantify various types of carbonate grains, matrix, and cement. In this study, the applicability and performance of DCNN-based object detection and image classification approaches are assessed with respect to carbonate compositional analysis. The training data comprised of more than 13,000 individually labelled objects from nearly 4000 carbonate petrographic images. The dataset is grouped into six and nine different classes for the image classification and object detection tasks, respectively. Even with a small and relatively imbalanced training set, the DCNN was able to achieve an F1 score of 92% for image classification and 84% mean precision for object detection by combining one-cycle policy, class weight, and label mixup-smoothing methods. This study highlights the inefficiency of image classification as an approach to replicating human description and classification of carbonate petrography. By contrast, DCNN-based object detection appears capable of approaching human speed and accuracy in the area of carbonate petrography because it is able to individually locate and identify different carbonate components with greater cost-efficiency, speed, and reproducibility than conventional (human) petrographic analysis.

We develop methods to incorporate paleogeographical constraints into numerical models of mantle convection. Through the solution of the convection equations, the models honor geophysical and geological data near the surface while predicting mantle flow and structure at depth and associated surface deformation. The methods consist of four constraints determined a priori from a plate history model: (1) plate velocities, (2) thermal structure of the lithosphere, (3) thermal structure of slabs in the upper mantle, and (4) velocity of slabs in the upper mantle. These constraints are implemented as temporally- and spatially-dependent conditions that are blended with the solution of the convection equations at each time step. We construct Earth-like regional models with oceanic and continental lithosphere, trench migration, oblique subduction, and asymmetric subduction to test the robustness of the methods by computing the temperature, velocity, and buoyancy flux of the lithosphere and slab. Full sphere convection models demonstrate how the methods can determine the flow associated with specific tectonic environments (e.g., back-arc basins, intraoceanic subduction zones) to address geological questions and compare with independent data, both at present-day and in the geological past (e.g., seismology, residual topography, stratigraphy). Using global models with paleogeographical constraints we demonstrate (1) subduction initiation at the Izu–Bonin–Mariana convergent margin and flat slab subduction beneath North America, (2) enhanced correlation of model slabs and fast anomalies in seismic tomography beneath North and South America, and (3) comparable amplitude of dynamic and residual topography in addition to improved spatial correlation of dynamic and residual topography lows.