Geophysical Constraints Research Articles

Can the Artificial Release of Fluorinated Gases Offset Global Cooling Due to Supervolcanic Eruptions?

Large volcanic eruptions, such as the prehistoric Yellowstone eruption, induce abrupt global cooling—by some estimates at a rate of ~1 °C/year, lasting for more than a decade. An abrupt global cooling of several °C—even if only lasting a few years—would present immediate, drastic stress on biodiversity and food production. This cooling poses a global catastrophic risk to human society beyond the immediate and direct impact of eruptions. Using a simple climate model, this paper discusses the possibility of counteracting large volcanic cooling with the intentional release of greenhouse gases. Longer-lived compounds (e.g., CO2 and CH4) are unsuitable for this purpose, but selected fluorinated gases (F-gases), either individually or in combinations, could be released at gigaton scale to offset large volcanic cooling substantially. We identify candidate F-gases (e.g., C4F6 and CH3F) and derive radiative and chemical properties of ‘ideal’ compounds matching specific cooling events. Geophysical constraints on manufacturing and stockpiling due to mineral availability are considered, alongside technical and economic implications based on present-day market assumptions. The effects and uncertainty due to atmospheric chemistry related to aerosol injection, F-gases release, and solar dimming are discussed in the context of large volcanic perturbation. The caveats and future steps using more complex chemistry–climate models are discussed. Despite the speculative nature of the magnitude and composition of F-gases, our conceptual analysis has implications for testing the possibility of mitigating certain global catastrophic cooling risks (e.g., nuclear winter, asteroid impact, and glacier transition) via intentional intervention.

An individual‐based model trained on multiple data sources estimates population connectivity and facilitates aggregation of harvest management units

Abstract Management boundaries are often delimited by political and social factors, whereas animal movements are affected by ecological and geophysical constraints. Thus, understanding connectivity among distinct management units is of considerable importance, particularly for harvested species, where quotas set in ignorance of connectivity may fail to meet management goals. We constructed an individual‐based model (IBM) to better understand wild turkey movements at large scales, benefiting from multiple data sources that are often available for harvested species. We built an IBM from spring seasonal movements of wild turkey, using data from ringed, radio‐, and GPS‐marked turkeys captured in Maine, USA. Our IBM accommodated variation in individual turkey response to landscape connectivity metrics and identified emergent migratory connectivity dynamics among harvest management regions. We calculated a low degree of connectivity among wildlife management districts (WMD) which, in combination with the substantial number of boundary crossings observed, indicated a more diffuse distribution of turkeys among WMDs. Synthesis and applications: Estimates of turkeys moving between districts provided a clear delineation of where immigration was strongest, identifying which WMDs should be managed as singular population units. This approach has widespread utility for any species or system where harvest management decisions are made at finer spatial scales than the movement dynamics affecting population processes.

Tectonics of the Mw 6.8 Al Haouz earthquake (Morocco) reveals minor role of asthenospheric upwelling

A reliable identification of the fault responsible for the magnitude 6.8 Al Haouz earthquake that struck Morocco on 8 September 2023 has so far been hampered by a lack of accurate tectonic analyses. Here we provide the first updated tectonic framework of the earthquake epicentral area based on original field data. We cast our results into the context of available geomorphological, thermochronological and geophysical constraints, and discuss the earthquake characteristics within the framework of competing tectonic models either based on asthenospheric upwelling or transpressional tectonics. We found that the Al Haouz earthquake was likely generated by rupture along a north-dipping high-angle fault, linking former fault planes belonging to an orogen-scale WSW-ESE transpressional shear zone. The geological evolution and seismotectonic structure of the region are largely governed by the oblique convergence of tectonic plates. The impact of asthenospheric upwelling, if any, remains limited and may only influence the geomorphological evolution of the Western High Atlas, but cannot explain the seismotectonic and geological features observed today at the surface, which are instead effects of transpressional tectonics.

Use of tsunami observations for resolving coseismic slip of recent and historic large submarine earthquakes

Tsunamis generated by seafloor displacements accompanying large submarine earthquakes provide sensitivity to absolute slip position and distribution for offshore faulting analogous to that of geodetic observations for landward faulting. Tsunami recordings at deep‐water and near‐shore ocean bottom pressure sensors and tide gauges, along with runup and inundation measurements, can now be reliably modeled using detailed bathymetric structures and robust numerical codes. As a result, tsunami observations now play an important role in quantifying coseismic slip distributions for large submarine earthquakes in subduction zones and other tectonic environments. Applications of joint modeling or inversion of seismic, geodetic and tsunami observations for recent major earthquakes are described, highlighting the specific contributions of the tsunami observations to source model resolution. Tsunami observations provide unique information on the up‐dip extent of earthquake coseismic slip on subduction zone megathrust faults and occurrence of near‐trench slip, which are usually not well constrained by seismic and land‐based geodetic signals. Tsunami signals also help to detect offshore slow slip that is not evident in seismic or land‐based geodetic data and to balance geophysical constraints on ruptures that extend from on‐shore to off‐shore. Tsunami runup measurements and stratigraphic deposits further provide unique constraints on large earthquake ruptures that occurred prior to modern geophysical instrumentation.

3-D spatial distribution of concealed ore-forming granitoid intrusion and structures determined by the CSAMT survey of Wuxu Sb-Zn-polymetallic ore district, South China

Wuxu ore district is an important lead-zinc-antimony polymetallic ore field within the Danchi metallogenic belt of South China. Although previous geochemical studies have indicated that the deposit type is medium to low-temperature magmatic-hydrothermal, the existing metallogenic models remain controversial due to the thick covering layer and lack of exposure of granitic bodies. To address this issue, this study utilised four CSAMT parallel profiles across the Wuxu anticline to perform two-dimensional inversion, yielding high-resolution subsurface resistivity models. By comparing the resistivity anomalies along the regional strike shown by the profiles and integrating geochemical and ore deposit structural analyses, a detailed geological interpretation was conducted. A significant high-resistivity body beneath the core region of the Wuxu anticline is interpreted to be the Yanshanian period granitic intrusion. There are westward- and eastward-dipping high-conductivity belts existing within the western and eastern wings and appear to sandwich the inferred high-resistivity code of the Wuxu anticline, respectively. They may represent the deep faults extending underground and are the result of decoupling of lithological interface within the anticline and thereby rock fracturing. The high-conductivity feature may be the result of filling with high-conductivity saline fluids, and possible hydrothermal alteration. The regional compression stress during the Caledonian period and the early Yanshanian period promoted the formation of the anticline and rock fracturing of both wings of the anticline. After that, the regional extension background of the late Yanshanian period made the fault zones reactivated, becoming the site and channels for magmatic intrusion, upward transporting of hydrothermal fluid and alteration, causing lead-zinc-antimony polymetallic ore formation. The exposed ore deposits above may also corroborate this view point. In summary, our resistivity model visually depicts the spatial distribution of the granitic intrusion controlling ore formation and major structures such as anticlines, fault zones, and hydrothermal alteration zones in the 3-D space, providing new geophysical constraints for solving the controversy over the metallogenic model for Wuxu ore district.

Rapid heating (

Ultrahigh-temperature (UHT) metamorphism is the most thermally extreme form of regional crustal metamorphism, with temperatures exceeding 900°C. The duration and heat source of UHT metamorphism are critical constraints on the tectonic evolution of orogenic systems. Here, we report the discovery of a sapphirine-bearing granulite from the east-central Himalaya which preserves UHT evidence. The reconstructed pressure-temperature-time path indicates that the temperature increased by almost 200°C within ~2 Ma which is consistent with rapid asthenospheric heat input. Numerical simulations illustrate potential mechanisms for such heating: juxtaposition of the deep crustal accretionary system with the upwelling asthenospheric mantle through newly developed apertures generated by slab break-off and/or associated vertical tearing of down-going Indian lithosphere. Spatial-temporal consistencies among the UHT metamorphic phases, postcollisional magmatism, geophysical constraints, and crustal deformation indicate that slab break-off or tearing controls broad swaths of Himalayan tectonics. The consequent upwelling asthenosphere may have been a significant heat source for the Miocene Himalaya and for similar ancient collisional orogenic systems.

Pre-stack seismic inversion based on model-constrained generative adversarial network

Pre-stack seismic inversion usually uses various traditional algorithms to estimate elastic parameters such as P-wave velocity, S-wave velocity, and density. It is hard to derive accurate elastic parameters due to their non-uniqueness and high dimensionality between elastic parameters and seismic data, the calculation of elastic parameters is inaccurate. Convolutional Neural Networks (CNNs) have high-dimensional feature space mapping capabilities, which are utilized to establish mapping relationships between seismic data and elasticity parameters. However, their effectiveness is greatly affected by label data, and at the same time, due to the lack of enough label data, resulting in a low degree of fitting between prediction results and real data. In addition, conventional seismic inversion methods based on CNNs lack physical model constraints, resulting in low accuracy and poor interpretability of prediction results. We propose a Cycle-consistent Generative Adversarial Network based on a geophysical mechanism (SeisInv-CycleGAN). Deterministic inversion results and labeled data are combined into hybrid geophysical data as a training set of SeisInv-CycleGAN with geophysical constraints. At the same time, the residual (seismic loss) between the seismic data synthesized by forward modeling and the actual data is used as part of the loss function. The SeisInv-CycleGAN does not require building an initial model, and it can achieve higher accuracy in prediction results with a small amount of labeled data.

Hydrogen and silicon are the preferred light elements in Earth’s core

Hydrogen is an important light element in the Earth’s core for its high cosmochemical abundance and strong affinity to iron under core-formation conditions. Thus, constraining the core composition requires knowledge on the distribution of hydrogen between the liquid outer core and solid inner core. Here we investigate the chemical equilibrium of hydrogen at the inner-core boundary by calculating the chemical potential of hydrogen in solid and liquid iron-hydrogen alloys, respectively, using first-principles molecular dynamic simulations and neural network methods. We find that hydrogen partitions preferentially into the outer core and provides a major contribution to the density jump across the inner-core boundary. Combining geophysical constraints, mineral physics data, and chemical equilibrium, we evaluated light element abundances in the outer and inner cores simultaneously. Our results suggest hydrogen and silicon are the preferred light elements in the core, implying a relatively reduced environment during the Earth’s accretion and core-formation processes.

Geochemical and geochronologic evidence for a contiguous northeastern Wyoming Province

The extent and nature of the Wyoming Province, an Archean craton in southwestern Laurentia, are poorly understood due to limited exposure between spatially isolated basement-cored uplifts. This lack of exposure has led to debate about whether the northeastern Wyoming Province is underlain by contiguous Archean crust or Proterozoic rocks and suture zone associated with the Trans-Hudson orogeny. To assess these models we analyzed samples recovered from drill cores in the buried northeastern Wyoming Province that straddle the proposed Proterozoic suture. Whole-rock geochemical and Nd-Pb-O isotopic data suggest the rocks formed by partial melting of > 3.5 Ga hydrated mafic to tonalitic sources similar to elsewhere in the northern Wyoming Province. Zircon U-Pb dates record Mesoarchean magmatism and metamorphism (3.0–2.8 Ga). Published whole-rock Rb-Sr, hornblende and biotite K-Ar, and apatite fission track dates suggest these rocks have not been heated above 300 °C since 2.5–2.1 Ga. Geophysical potential field data are consistent across northeastern Wyoming contrasting with major discontinuities associated with documented Proterozoic orogens on all other margins of the Wyoming Province. Hence, geochemical, isotopic, and geochronologic data along with geophysical imaging can be most simply interpreted in terms of continuous Archean crust in the northeastern Wyoming Province. A geophysically-imaged reflector east of the Bighorn Mountains may juxtapose Archean terranes with similar ages and sources, similar to the boundary between the Montana metasedimentary terrane and Beartooth-Bighorn magmatic zone in the northwestern Wyoming Province. This work emphasizes the value of integrating geologic and geophysical constraints to constrain Archean provinces and their tectonic evolution.

Effects of crustal rheology and fault strength on the formation of the Qaidam Basin: Results from 2-D mechanical modeling

The Cenozoic Qaidam Basin on the northern Tibetan Plateau is bounded by the foreland thrust belts of the East Kunlun Mountains to the south and Qilian Mountains to the north. In this basin, the maximum thickness of the accumulated sedimentary strata is ∼17 km near the geographic center of the basin. Such a basin architecture challenges the classic foreland-basin deformation model, wherein substantial sediments should be localized near the basin margins owing to thrust loads on a bending plate. Based on the available geological and geophysical constraints, we built a two-dimensional viscoelastoplastic finite element model to investigate whether the unique pattern of the Qaidam Basin was mechanically related to the rheological structure of the crust and frictional strength of the basin-bounding thrust fault. By testing reasonable bounds on the viscosity of the ductile crust and effective frictional coefficient of the thrust faults, we modeled a basin framework that fits well with that of Qaidam Basin. The model parameters include a layered rheological structure, with the viscosity of the ductile crust of the basin being two orders of magnitude greater than that of the adjacent mountain belts and an effective frictional coefficient of >0.2. Furthermore, the initial width of the basin significantly affected its deformation style. Therefore, the layered rheological structure of the crust and frictional strength of the basin-bounding thrust faults played important roles in the formation of the Qaidam Basin. Our findings suggest that the inelastic behavior of the crust affects foreland deformation, highlighting that the classic foreland basin model must be revised for some tectonic settings.

Geologic Provinces Beneath the Greenland Ice Sheet Constrained by Geophysical Data Synthesis

AbstractPresent understanding of Greenland's subglacial geology is derived mostly from interpolation of geologic mapping of its ice‐free margins and unconstrained by geophysical data. Here we refine the extent of its geologic provinces by synthesizing geophysical constraints on subglacial geology from seismic, gravity, magnetic and topographic data. North of 72°N, no province clearly extends across the whole island, leaving three distinct subglacial regions yet to be reconciled with margin geology. Geophysically coherent anomalies and apparent province boundaries are adjacent to the onset of faster ice flow at both Petermann Glacier and the Northeast Greenland Ice Stream. Separately, based on their subaerial expression, dozens of unusually long, straight and sub‐parallel subglacial valleys cross Greenland's interior and are not yet resolved by current syntheses of its subglacial topography.

New Numerically Derived Scaling Relationships for Impact Basins on Mars

AbstractMost impact basins are believed to have formed during the early epochs of planetary evolution. The planet's gravity, internal structure, and thermal regime have the strongest control over their formation. Because of this, we can use the geophysical constraints on Mars' interior composition, structure, and geophysical evolution derived from the InSight mission to better understand the formation of impact basins on the planet. To achieve this, we performed numerical simulations of large impacts using the iSALE shock physics code. We investigated the effects of temperature and crustal thickness variations on impact basin size and morphology. Our scaling relationships indicate that: (a) basins formed in a warmer crust have larger final diameters in comparison to basins formed in a colder crust, a difference that is further accentuated as basin size gets bigger; and (b) the largest impact basins on Mars were created by impactors ranging from 35 to 680 km in diameter, up to ∼32% larger than estimates based on classical scaling. Our results expand the current understanding of the extent of early and large impact bombardment on Mars and provide a more comprehensive knowledge of impact basin formation on planetary surfaces.

A Machine Learning-Based Thermobarometer for Magmatic Liquids

Abstract Experimentally calibrated models to recover pressures and temperatures of magmas are widely used in igneous petrology. However, large errors, especially in barometry, limit the capacity of these models to resolve the architecture of crustal igneous systems. Here, we apply machine learning to a large experimental database to calibrate new regression models that recover P–T of magmas based on melt composition plus associated phase assemblage. The method is applicable to compositions from basalt to rhyolite, pressures from 0.2 to 15 kbar, and temperatures of 675°C to 1400°C. Testing and optimisation of the model with a filter that removes estimates with standard deviation above the 50th percentile show that pressures can be recovered with root-mean-square-error (RMSE) of 1.1 to 1.3 kbar and errors on temperature estimates of 21°C. Our findings demonstrate that, given constraints on the coexisting mineral assemblage, melt chemistry is a reliable recorder of magmatic variables. This is a consequence of the relatively low thermodynamic variance of natural magma compositions despite their relatively large number of constituent oxide components. We apply our model to two contrasting cases with well-constrained geophysical information: Mount St. Helens volcano (USA), and Askja caldera in Iceland. Dacite whole-rocks from Mount St Helens erupted 1980 to 1986, inferred to represent liquids extracted from cpx–hbl–opx–plag–mt–ilm mush, yield melt extraction source pressures of 5.1 to 6.7 kbar in excellent agreement with geophysical constraints. Melt inclusions and matrix glasses record lower pressures (0.7–3.8 kbar), consistent with magma crystallisation within the upper reaches of the imaged geophysical anomaly and during ascent. Magma reservoir depth estimates for historical eruptions from Askja match the location of seismic wave speed anomalies. Vp/Vs anomalies at 5 to 10 km depth correspond to hot (~990°C) rhyolite source regions, while basaltic magmas (~1120°C) were stored at 7 to 17 km depth under the caldera. These examples illustrate how our model can link petrology and geophysics to better constrain the architecture of volcanic feeding systems. Our model (MagMaTaB) is accessible through a user-friendly web application (https://igdrasil.shinyapps.io/MagmaTaBv4/).

Geophysical constraints on continental rejuvenation in central China: Implications for outward growth of the Tibetan Plateau

Continental rejuvenation results from the tectonic reactivation of crustal structures and lithospheric reworking by mantle flow. Geochemical observations and field mapping have traditionally provided the primary evidence for the secular evolution of crustal composition and tectonic processes during continental rejuvenation. Nonetheless, the impact of continental rejuvenation on the observed present-day strain rate and orogenic-scale lithospheric structure has not been well constrained. The pre-existing E-W−trending Central China Orogenic Belt has been overprinted by the N-S−trending Central Longitudinal Seismic Belt and constitutes the intracontinental West Qinling Syntaxis in central China, where the tectonic setting changes eastward from contraction to extension. Combining updated global positioning system data and high-resolution crustal seismic tomography, we reveal a modern continental rejuvenation process within the West Qinling Syntaxis in central China. The northward extrusion of the Tibetan Plateau’s weak lithospheric layer (middle-lower crust and lithospheric mantle) of southwestern China relative to the rigid Sichuan Basin/Ordos Block of the eastern West Qinling Syntaxis results in regional dextral shearing that shapes the Central Longitudinal Seismic Belt and defines the eastern Tibetan Plateau margin. The pre-existing E-W−trending Central China Orogenic Belt has been preserved above the brittle-ductile transition zone, and the northward movement of the deep lithospheric layer drives the deformation of the upper crust in the West Qinling Syntaxis. Our results, along with previous studies, suggest the presence of an intracontinental lithospheric interchange structure in central China. The continental rejuvenation of the West Qinling Syntaxis results from a combination of fault reactivation in the upper crust (Stage I, Eocene−Oligocene) and reworking of the deep lithosphere (Stage II, middle−late Miocene) related to the plateau-wide shift in stress accommodation ultimately driven by the redistribution of mass outward from the central Tibetan Plateau. At present, the transition zone between the high- and low-velocity anomalies along the Central Longitudinal Seismic Belt not only shapes the landscape boundary but controls the size and recurrence interval of earthquakes within the West Qinling Syntaxis in central China.

Gravity Investigation to Characterize Enceladus's Ocean and Interior

A key objective for the future exploration of the icy moon Enceladus is the characterization of the habitable conditions in its internal ocean. Radio science instrumentation on board a spacecraft in a low-altitude orbit about Enceladus would enable gravity measurements that are fundamental to providing constraints on its internal structure. We present here the concept of operations and expected results of the gravity investigation for a New Frontiers–class mission. Numerical simulations are carried out to determine the gravity field in spherical harmonics to degree and order 30 and the Love number k 2 and its phase. By combining Enceladus’s shape measured by Cassini and the geophysical constraints obtained through the processing of the simulated radio science data, a Bayesian inference network is used for the interior model inversion. Our results indicate that the gravity investigation would enable tight constraints on core radius and density, ocean depth and density, and ice shell rigidity. By assuming a high core rigidity and a preliminary modeling of dissipation in the ice shell, our interior model inversion also yields information on the ice shell viscosity. Further data on the hydrosphere properties might be gathered through optical navigation data by accurately measuring Enceladus’s orientation model.

Middle Miocene renewed subduction of arabian-eurasian plates: Implications for convergent tectonic mechanisms, tectonically-induced fluid overpressures, and sediment recycling dynamics

The complex interplay between thrust wedge deformation, tectonically-induced fluid overpressures, and the complex phases of wedge-top sedimentary depositional systems in active plate convergent margins presents a formidable challenge. The offshore Makran accretionary wedge (MAW) is the world's largest, formed by the convergence of Arabian plate underneath Eurasian plate in the Early Cretaceous, followed by Middle Miocene renewed subduction. In this research, geological and geophysical constraints derived from seismic and borehole stratigraphy, seismic attributes, depth structural maps, isopach maps and 3D structural geological models reveal that Middle Miocene renewed subduction controls thrust wedge deformation. Additionally, it influences deep fluid overpressure and the deformational pattern of wedge-top sedimentary depositional systems. The seismic data analysis suggested that the basement pre-kinematic Himalayan turbidities depositional system (HTDs) have been tectonically incorporated into the accretionary wedge by Middle Miocene renewed subduction, which gave rise to complex folding and N-dipping imbricate thrusting features. These N-dipping imbricate thrust faults initiated sediment underplating and served as major structural pathways for deep fluid overpressure migration, resulting in the formation of fluid escape pipes. The fluid escape pipes exhibit conical and bifurcating geometries with inclination angles of approximately 73°, 90°, and 100°, signifying the upward migration of deep overpressured fluids. Furthermore, the Middle Miocene to Pliocene syn-kinematic piggyback and growth depositional system (PGDs) show progressive thickening (Max ∼3600 m) towards onshore MAW. It lies unconformable above the pre-kinematic HTDs and carries valuable sedimentation records, growth stratal patterns, onlap terminations, and truncations against growing structures, revealed the timing and spatial deformation of thrust faults in 3D space. In contrast, the Pleistocene to Recent post-kinematic progradational depositional system (PDs) exhibits clinoforms marked by top-lapping, on-lapping, and down-lapping reflection, which are primarily controlled by sediment recycling dynamics and sea level change since Pleistocene. It shows progressive thickening (Max ∼1800m) towards the Makran subduction zone and lies unconformably above the syn-kinematic PGDs. Abridging the analyses, the tectonic reconstruction and deformation mechanisms model shows four phases, e.g. (1) initiation of subduction and accretionary wedge formation in Eocene, (2) N-dipping imbricate thrust faults and deep fluid overpressure triggered by Middle Miocene renewed subduction, (3) development of piggyback basins and growth stratal geometries in Middle Miocene to Pliocene, (4) thrust deformation ceases and development of progradational depositional system.

Geophysical Constraints to the Geological Evolution and Genesis of Rare Earth Element–Thorium–Uranium Mineralization in Pegmatites at Alces Lake, SK, Canada

This investigation establishes an integrated method for rare earth elements (REE) exploration through a very promising and advanced exploration prospect in the Alces Lake area (SK, Canada) by assessing the integrated analysis of several multisource geophysical datasets. The resulting outcome provides important lithostructural information to the well-exposed, mineralized middle-to-lower crust at Alces Lake, comprising deep-seated poly-phase folds, ductile shear zones, and brittle faults. Geophysical–geological models of the Alces Lake property were constructed at different scales. The area of interest is located within the Beaverlodge Domain, about 28 km north of the Athabasca Basin’s northern margin. It contains some of the highest-grade rare earth elements (REE) in the world with the REE hosted predominantly in monazites within quartzo-feldspathic granitic to biotite–garnet–monazite–zircon-rich restite-bearing/cumulate mush melt pegmatites of anatectic origin (abyssal). Geophysical magnetic, gravity, and radiometric data were used together with Shuttle Radar Topography Mission (SRTM) images to facilitate the processing, modeling, and interpretation. Consequently, major structures were identified at different scales; however, the emphasis was given to studying those at the district/camp scale. The REE zones discovered to date occur within a large district-scale refolded synformal anticline. The eastern limb of this folded structure comprises a 30–40 km long, NW-trending shear zone/fault corridor with deep-seated structural crustal roots that may have served as the major pathway for ascending fluids/melts and facilitated the emplacement of mineralization. Thus, shear zones, faults, and folds in combination with lithological contacts/rheological contrasts appear to control residual/cumulate pegmatite emplacement and monazite deposition. Anomalies obtained from the airborne equivalent thorium survey data prove to be the most useful for REE pegmatite exploration. The results herein provide new interpretation and modeling perspectives leading to a better understanding of the distribution and lithostructural controls of REE on the property, and to new guidelines for future exploration programs at Alces Lake and elsewhere in northern Saskatchewan.

New constraints on the shear wave velocity structure of the Ivrea geophysical body from seismic ambient noise tomography (Ivrea-Verbano Zone, Alps)

SUMMARY We performed seismic ambient noise tomography to investigate the shallow crustal structure around the Ivrea geophysical body (IGB) in the Ivrea-Verbano Zone (IVZ). We achieved higher resolution with respect to previous tomographic works covering the Western Alps, by processing seismic data collected by both permanent and temporary seismic networks (61 broad-band seismic stations in total). This included IvreaArray, a temporary, passive seismic experiment designed to investigate the IVZ crustal structure. Starting from continuous seismic ambient noise recordings, we measured and inverted the dispersion of the group velocity of surface Rayleigh waves (fundamental mode) in the period range 4–25 s. We obtained a new, 3-D vS model of the IVZ crust via the stochastic neighbourhood algorithm (NA), with the highest resolution between 3 to 40 km depth. The fast and shallow shear wave velocity anomaly associated with the IGB presents velocities of 3.6 km s−1 directly at the surface, in remarkable agreement with the location of the exposed lower-to-middle crustal and mantle outcrops. This suggests a continuity between the surface geological observations and the subsurface geophysical anomalies. The fast IGB structure reaches vS of 4 km s−1 at 20–25 km depth, at the boundary between the European and Adriatic tectonic plates, and in correspondence with the earlier identified Moho jump in the same area. The interpretation of a very shallow reaching IGB is further supported by the comparison of our new results with recent geophysical investigations, based on receiver functions and gravity anomaly data. By combining the new geophysical constraints and the geological observations at the surface, we provide a new structural interpretation of the IGB, which features lower crustal and mantle rocks at upper crustal depths. The comparison of the obtained vS values with the physical properties from laboratory analysis of local rock samples suggests that the bulk of the IGB consists of a combination of mantle peridotite, ultramafic and lower crustal rocks, bound in a heterogeneous structure. These new findings, based on vS tomography, corroborate the recent interpretation for which the Balmuccia peridotite outcrops are continuously linked to the IGB structure beneath. The new outcomes contribute to a multidisciplinary framework for the interpretation of the forthcoming results of the scientific drilling project DIVE. DIVE aims at probing the lower continental crust and its transition to the mantle, with two ongoing and one future boreholes (down to 4 km depth) in the IVZ area, providing new, complementary information on rock structure and composition across scales. In this framework, we constrain the upper crustal IGB geometries and lithology based on new evidence for vS, connecting prior crustal knowledge to recent active seismic investigations.

Geophysical Constraints on Decarbonized Systems—Building Spatio-Temporal Uncertainties into Future Electricity Grid Planning

Geophysical Constraints on Decarbonized Systems—Building Spatio-Temporal Uncertainties into Future Electricity Grid Planning

Cosmogenic Nuclide Tracking of Sediment Recycling From a Frontal Siwalik Range in the Northwestern Himalaya

AbstractThe Himalayan orogen exports millions of tons of sediment annually to the Indo‐Gangetic foreland basin, derived from both hinterland and foreland fold‐thrust belts (FTB). Although erosion rates in the hinterland are well‐constrained, erosion rates in the foreland FTB and, by extension, the sediment flux have remained poorly constrained. Here, we quantified erosion rates and sediment flux from the Mohand Range in the northwestern Himalaya by modeling and measuring the cosmogenic radionuclide (CRN) 10Be and 26Al concentrations in modern fluvial sediments. Our model uses local geological and geophysical constraints and accounts for CRN inheritance and sediment recycling, which enables us to determine the relative contributions of the hinterland and foreland FTB sources to the CRN budget of the proximal foreland deposits. Our model predictions closely match measured concentrations for a crustal shortening rate across the Mohand Range of 8.0 ± 0.5 mm yr−1 (i.e., approximately 50% of the total shortening across the Himalaya at this longitude) since Ma. This shortening implies a spatial gradient in erosion rates ranging from 0.42 ± 0.03 to 4.92 ± 0.34 mm yr−1, controlled by the geometry of the underlying structure. This erosion pattern corresponds to an annual sediment recycling of ∼2.0 megatons from the Mohand Range to the downstream Yamuna foreland. Converted to sediment fluxes per unit width along the Himalaya, the foreland FTB accounts for ∼21% ± 5% of the total flux entering the foreland. Because these sediments have lower 10Be concentrations than hinterland‐derived sediment, they would lead to ∼14% overestimation of 10Be‐derived erosion rates, based on Yamuna sediments in the proximal foreland.