Subsurface State Research Articles

The role of mobile muds in basin form and fill in the tectonic and sedimentologically complex region of the southeastern margin of the Caribbean Plate, Trinidad and Tobago, West Indies

The southern portions of the Barbados Accretionary Prism and the eastern deepwater region of Trinidad and Tobago is a seafloor terrane marked by scattered mud volcanoes, mud walls, and mud canopies that rise several hundred meters above the modern seafloor. We use 3D seismic and recently drilled wells to examine this area dividing it into three tectonomorphic provinces based on the nature of the mud architecture, the structuration, basin types developed and sedimentary systems occupying the basins. Mud architectures become more massive as one moves northward toward the active tectonic prism. The mud wall in the far northern Province 1 is classified as a mud canopy system; similar to those documented by Morley et al. (2022) in the North Sabah Pagasa Wedge, Borneo. The Trinidad mud canopy exhibits sutures in the mud mass margins (akin to salt “sutures”) that separate sub-masses likely composed of different materials and exhibiting different pressures and properties. These sub-masses and associated sutures are imageable using the sweetness attributes on the near stack seismic volume. Anticlines in Provinces 2 and 3, like those that are the dominant hydrocarbon-producing structures onshore and offshore Trinidad, show mobile mud cores and are associated with extensional and compressional faulting that provide pathways for hydrocarbon migration into the surrounding lutokinetic stratigraphy. A bottom simulating reflector associated with the base of the gas hydrate zone; extensive across the study area, may be a good indicator of fluid flux activity and subsurface pressure states in mud masses. It is important to employ all the seismic data volumes; full, near, mid, and far seismic data stacks, to assess different aspects of these mud diapiric systems.

Impact of ocean in-situ observations on ECMWF sub-seasonal forecasts

We assess for the first time the impact of in-situ ocean observations on European Centre for Medium-Range Weather Forecasts (ECMWF) sub-seasonal forecasts of both ocean and atmospheric conditions. A series of coupled reforecasts have been conducted for the period 1993-2015, in which different sets of ocean observations were withdrawn in the production of the ocean initial conditions. Removal of all ocean in-situ observations in the initial conditions leads to significant degradation in the forecasts of ocean surface and subsurface mean state at lead times from week 1 to week 4. The negative impact is predominantly caused by the removal of the Argo observing system in recent decades. Changes in the mean state of atmospheric variables are comparatively small but significant in the forecasts of lower and upper atmospheric circulation over large regions. Our results highlight the value of continuous, real-time in-situ observations of the surface and subsurface ocean for coupled forecasts in the sub-seasonal range.

Progressive Induction Hardening: Measurement and Alteration of Residual Stresses

Progressive induction hardening is an in-line steel heat treatment method commonly used to surface harden powertrain components. It produces a martensitic case layer with a sharp transition zone to the base material. This rapid process will induce large residual stresses, where a compressive state in the case layer will shift to a tensile state in the transition zone. For fatigue performance, it is important to quantify the magnitude and distribution of these stresses, and moreover how they depend on material and processing parameters. In this work, x-ray diffraction in combination with a layer removal method is used for efficient and robust quantification of the subsurface stress state, which combines electropolishing with either turning or milling. Characterization is done on C45E steel samples that were progressively induction hardened using either a fast or slow (27.5 or 5 mm/s, respectively) scanning speed. The results show that although the hardening procedures will meet arbitrary requirements on surface hardness, case depth and microstructure, the subsurface tensile stress peak magnitude is doubled when using a fast scanning speed. However, the near-surface compressive residual stresses are comparable. In addition, the subsurface tensile residual stress peak is compared with the on-surface tensile stresses in the fade-out zone.

An efficient 3D finite element procedure for simulating wheel–rail cyclic contact and ratcheting

Various models for simulating rail ratcheting behaviour were developed to study rolling contact fatigue (RCF) damage in rails. However, limitations remain in terms of the accuracy of wheel–rail contact modelling and computational efficiency of the cyclic loading simulation. This study developed an efficient 3D finite element (FE) procedure to simulate ratcheting in rails subjected to numerous load cycles. The procedure simulates a wheel rolling repeatedly over a rail section with updated stress–strain states, enabling automatically executed cyclic loading simulation given a predefined number of cycles. To ensure the accuracy of the contact modelling, the effect of meshing schemes on subsurface stress distribution was examined. In addition, the FE contact model with the selected meshing scheme, which balances accuracy and computational efficiency, was verified against the widely accepted CONTACT program. Subsequently, a non-linear kinematic hardening (NLKH) steel material was used in the FE model for ratcheting simulations with up to 100 wheel-loading cycles. The rail surface and subsurface stress states were replicated under partial-slip wheel–rail rolling contact conditions with traction coefficients of 0.10, 0.20 and 0.35, respectively. The ratcheting behaviour was extensively analysed in terms of plastic deformation, contact patch evolution, and ratcheting rates. The simulated plastic deformation was found to alter the contact geometry and thus contact stresses, which in turn affect further accumulation of plastic deformation and subsequent ratcheting strains. These findings highlighted the importance of considering the interplay between the rail ratcheting behaviour of the rail and evolving contact conditions for predicting ratcheting and RCF damage in rails.

Untangling the environmental and tectonic drivers of the Noto earthquake swarm in Japan.

The underlying mechanism of the ongoing seismic swarm in the Noto Peninsula, Japan, which generates earthquakes at 10 times the average regional rate, remains elusive. We capture the evolution of the subsurface stress state by monitoring changes in seismic wave velocities over an 11-year period. A sustained long-term increase in seismic velocity that is seasonally modulated drops before the earthquake swarm. We use a three-dimensional hydromechanical model to quantify environmentally driven variations in excess pore pressure, revealing its crucial role in governing the seasonal modulation with a stress sensitivity of 6 × 10-9 per pascal. The decrease in seismic velocity aligns with vertical surface uplift, suggesting potential fluid migration from a high-pore pressure zone at depth. Stress changes induced by abnormally intense snow falls contribute to initiating the swarm through subsequent perturbations to crustal pore pressure.

Enhancing high-temperature fretting wear resistance of TC21 titanium alloys by laser cladding self-lubricating composite coatings

Fretting wear is one of the main factors restricting the fastener service reliability for the high-temperature titanium alloys. Enhancing the high-temperature fretting wear resistance of titanium alloys with laser cladding coatings is an effective method. Thus, in this work, a self-lubricating composite coating was designed and fabricated on the surface of TC21 titanium alloy with laser cladding. Special attentions were made to the effect of the fractions of the hard and self-lubricating phases on high-temperature fretting resistance. In addition, the contribution of laser cladding power was also concerned. The primary phases of the self-lubricating composite coatings are α-Co, CaF2, WC, TiC and Cr23C6, and the reasons for the formation of TiC and Cr23C6 are investigated from a thermodynamic point. The prepared composite coatings are dense and uniform. The crystallographic characteristics and the volume fractions of the hard phases (WC, TiC and Cr23C6) and self-lubricating phase (CaF2) were revealed, and the effect of grain size on the coating properties was analyzed. The hardness of the coatings were all maintained around 700 HV0.2, which was significantly higher than that of the substrate. The composite coating significantly improves the friction and wear resistance at 500 ℃. The hard phase was also found to be the most resistant to friction and wear. The frictional performance of the coating gets better when the hard phase and the self-lubricating phase increase. However, the hard phase contributes more to the high-temperature friction and wear resistance than the self-lubricating phase and the wear volume is 1/40 that of substrate. For the laser cladding power, the higher power can lead to a higher fraction of CaF2 but a lower fraction of hard phase, causing a decline in COF but an increase in wear volume. The wear mechanisms of the coating and the substrate were analyzed, the subsurface stress distribution state of the wear scars were further analyzed. The current results present new insights on the anti-wear design strategy to the composite coatings consist of hard and self-lubricating phase.

The Porothermoelastic Stress Evolution of Inhomogeneous Solid due to Near-Surface Reservoir

Anisotropic and heterogeneous solids, comprising a minimum of two or more elements with different properties, appear pervasively in rock materials including the pore structure and mineral composition of quartz and granite and usually have extensive applications in construction aggregates, dimension stone, geotechnical engineering, and petrology/geology. The elastic stress fields of inhomogeneous materials or composites inevitably change due to the presence of heterogeneous microstructure under applied external conditions, and hence the total mechanical and physical properties are affected by the temperature and pore pressure fluctuations beneath the surface. In this paper, the thermoporoelastic (TPE) approach on the basis of eigenstrain concept is introduced to predict the stress fields in fluid-saturated porous geologic materials like hydrocarbon reservoirs or aquifers by accounting for the coupling between thermal, poroelastic, and mechanical effects. It is an extension of the micromechanical theory that also incorporates thermal effects like pore pressure changes and temperature alternations. In addition, the TPE approach provides an important multiphysical modeling tool for understanding subsurface fluid-rock interactions and stress states in applications like unconventional oil/gas and geothermal energy.

Enhanced Microseismicity During Production Pumping Cessation at the San Emidio Geothermal Field (Nevada, USA) in December 2016

AbstractTectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.

Reconstruction of subsurface ocean state variables using Convolutional Neural Networks with combined satellite and in situ data

Subsurface ocean measurements are extremely sparse and irregularly distributed, narrowing our ability to describe deep ocean processes and thus also limiting our understanding of the role of ocean and marine ecosystems in the Earth system. To overcome these observational limitations, neural networks combining remotely-sensed surface measurements and in situ vertical profiles are increasingly being used to retrieve high-quality three-dimensional estimates of the ocean state. This study proposes a convolutional neural network (CNN) architecture for the reconstruction of vertical profiles of temperature and salinity starting from surface observation-based data. The model is trained on satellite and in situ data collected between 2005 and 2020 in the Atlantic Ocean. Rather than using spatially gridded in situ observations, we use directly measured vertical profiles. Different combinations of surface variables are analyzed and compared in order to determine the most effective inputs for the CNN. Furthermore, the relative importance of each of these variables in the vertical reconstruction is assessed using Shapley values, originally developed in the framework of cooperative game theory. The model performance is shown to be superior to current state-of-the-art methods and the same approach can easily be extended to other basins or to the global ocean.

Nanometric Cutting Mechanism of Cerium–Lanthanum Alloy

Cerium–lanthanum alloy is widely used in the green energy industry, and the nanoscale smooth surface of this material is in demand. Nanometric cutting is an effective approach to achieving the ultra-precision machining surface. Molecular dynamics (MD) simulation is usually used to reveal the atomic-scale details of the material removal mechanism in nanometric cutting. In this study, the effects of cutting speed and undeformed chip thickness (UCT) on cutting force and subsurface deformation of the cerium–lanthanum alloy during nanometric cutting are analyzed through MD simulation. The results illustrate that the dislocations, stacking faults, and phase transitions occur in the subsurface during cutting. The dislocations are mainly Shockley partial dislocation, and the increase of temperature and pressure during the cutting process leads to the phase transformation of γ-Ce (FCC) into β-Ce (HCP) and δ-Ce (BCC). β-Ce is mainly distributed in the stacking fault area, while δ-Ce is distributed in the boundary area between the dislocation atoms and γ-Ce atoms. The cutting speed and UCT affect the distribution of subsurface damage. A thicker deformed layer including dislocations, stacking faults and phase-transformation atoms on the machined surface is generated with the increase in the cutting speed and UCT. Simultaneously, the cutting speed and UCT significantly affect the cutting force, material removal rate, and generated subsurface state. The fluctuations in the cutting force are related to the generation and disappearance of dislocations. This research first studied the nanometric cutting mechanism of the cerium–lanthanum ally, providing a theoretical basis for the development of ultra-precision machining techniques of these materials.

Earthquake focal mechanisms with distributed acoustic sensing

Earthquake focal mechanisms provide critical in-situ insights about the subsurface faulting geometry and stress state. For frequent small earthquakes (magnitude< 3.5), their focal mechanisms are routinely determined using first-arrival polarities picked on the vertical component of seismometers. Nevertheless, their quality is usually limited by the azimuthal coverage of the local seismic network. The emerging distributed acoustic sensing (DAS) technology, which can convert pre-existing telecommunication cables into arrays of strain/strain-rate meters, can potentially fill the azimuthal gap and enhance constraints on the nodal plane orientation through its long sensing range and dense spatial sampling. However, determining first-arrival polarities on DAS is challenging due to its single-component sensing and low signal-to-noise ratio for direct body waves. Here, we present a data-driven method that measures P-wave polarities on a DAS array based on cross-correlations between earthquake pairs. We validate the inferred polarities using the regional network catalog on two DAS arrays, deployed in California and each comprising ~ 5000 channels. We demonstrate that a joint focal mechanism inversion combining conventional and DAS polarity picks improves the accuracy and reduces the uncertainty in the focal plane orientation. Our results highlight the significant potential of integrating DAS with conventional networks for investigating high-resolution earthquake source mechanisms.

Comprehensive excited state carrier dynamics of 2D selenium: One-photon and multi-photon absorption regimes

Semiconductors play a critical role in optoelectronic applications, and recent research has identified group-VI 2D semiconductors as promising materials for this purpose. Here, we report the comprehensive excited state carrier dynamics of bilayer, two-dimensional (2D) selenium (Se) in one-photon and multi-photon absorption regimes using transient reflection (TR) spectroscopy. Carrier lifetime obtained from TR measurement is used to theoretically predict the photo-responsivity for 2D Se photo-detectors operating in the one-photon-absorption regime. We also calculate a giant two-photon absorption cross section of 2.9×105 GM at 750 nm hinting possible application of 2D Se as a sub-bandgap photo-detector. The carrier recombination process is dominated by surface and sub-surface defect states in one- and multi-photon absorption regimes, respectively, resulting nearly one order increased carrier lifetime in a three-photon-absorption regime (1700 ps) compared to a one-photon-absorption regime (103 ps). Femtosecond Z-scan measurement shows saturation behavior for above bandgap excitation, further indicating the possibility of 2D Se as a saturable absorber material for passive Q-switching. Our study provides comprehensive insight into the excited state carrier dynamics of bilayer 2D Se and highlights its potential as a versatile material for various linear and non-linear optoelectronic applications.

Rate-induced tipping to metastable Zombie fires

Zombie firesin peatlands disappear from the surface, smoulder underground during the winter, and ‘come back to life’ in the spring. They can release hundreds of megatonnes of carbon into the atmosphere per year and are believed to be caused by surface wildfires. Here, we propose rate-induced tipping (R-tipping) to a subsurface hot metastable state in bioactive peat soils as a main cause of Zombie fires. Our hypothesis is based on a conceptual soil-carbon model subjected to realistic changes in weather and climate patterns, including global warming scenarios and summer heatwaves. Mathematically speaking, R-tipping to the hot metastable state is a genuine nonautonomous instability, due to crossing an elusive quasi-threshold, in a multiple-timescale dynamical system. To explain this instability, we provide a framework combining a special compactification technique with concepts from geometric singular perturbation theory. This framework allows us to reduce an R-tipping problem due to crossing a quasi-threshold to a heteroclinic orbit problem in a singular limit. We identify generic cases of tracking–tipping transitions via: (i) unfolding of a codimension-twoheteroclinic folded-saddle-node type-I singularityfor global warming and (ii) analysis of a codimension-onesaddle-to-saddle hetroclinic orbitfor summer heatwaves, in turn revealing new types of excitability quasi-thresholds.

Headcut Migration in Earthen Embankment Induced by Varying Sub-surface and Seepage State under Overflow

Headcut Migration in Earthen Embankment Induced by Varying Sub-surface and Seepage State under Overflow

Hydrological forecasting at impact scale: the integrated ParFlow hydrological model at 0.6 km for climate resilient water resource management over Germany

In the context of the repeated droughts that have affected central Europe over the last years (2018–2020, 2022), climate-resilient management of water resources, based on timely information about the current state of the terrestrial water cycle and forecasts of its evolution, has gained an increasing importance. To achieve this, we propose a new setup for simulations of the terrestrial water cycle using the integrated hydrological model ParFlow/CLM at high spatial and temporal resolution (i.e., 0.611 km, hourly time step) over Germany and the neighboring regions. We show that this setup can be used as a basis for a monitoring and forecasting system that aims to provide stakeholders from many sectors, but especially agriculture, with diagnostics and indicators highlighting different aspects of subsurface water states and fluxes, such as subsurface water storage, seepage water, capillary rise, or fraction of plant available water for different (root-)depths. The validation of the new simulation setup with observation-based data monthly over the period 2011–2020 yields good results for all major components of the terrestrial water cycle analyzed here, i.e., volumetric soil moisture, evapotranspiration, water table depth, and river discharge. As this setup relies on a standardized grid definition and recent globally available static fields and parameters (e.g., topography, soil hydraulic properties, land cover), the workflow could easily be transferred to many regions of the Earth, including sparsely gauged regions, since ParFlow/CLM does not require calibration.

The subsurface thermal state of Svalbard and implications for geothermal potential

Svalbard is a High Arctic Archipelago at 74–81°N and 15–35 °E under the sovereignty of Norway. All settlements in Svalbard, including the capital of Longyearbyen (population 2400), currently have isolated energy systems with coal or diesel as the main energy source. Geothermal energy is considered as a possible alternative for electricity production, as a heat source in district heating systems or harnessed for heating and cooling using geothermal heat pump installations. In this contribution we present the until now fragmented data sets relevant to characterize and assess the geothermal potential of Svalbard. Data sets include petroleum and deep research boreholes drilled onshore Svalbard, 14 of which have recorded subsurface temperature data at depths below 200 m. Geothermal gradients on Spitsbergen vary from 24 °C/km in the west to 55 °C/km in the south-east, with an average of 33 °C/km. Four deep research boreholes were fully cored and analyzed for thermal conductivity. These analyses were complemented by thermal conductivity calculated from wireline logs in selected boreholes and four measurements on outcrop samples. 1D heat flow modelling on five boreholes calibrated with the measured thermal conductivities offers insights into heat transfer through the heterogeneous sedimentary succession. Offshore petroleum boreholes in the south-western Barents Sea and marine heat flow stations around Svalbard provide a regional framework for discussing spatial variation in heat flow onshore Svalbard, with emphasis on the effects of erosion and deposition on the thermal regime. We conclude that Svalbard's geology is well suited for geothermal exploration and potential production, though challenges related to permafrost, the presence of natural gas, heterogeneous reservoir quality and strongly lateral varying heat flow need to be adequately addressed prior to geothermal energy production. Specifically for Longyearbyen, high geothermal gradients of 40–43 °C/km in the nearest borehole (DH4) suggest promising sub-surface thermal conditions for further exploration of deep geothermal potential near the settlement.

Unidirectional Seebeck effect from Rashba spin-orbit coupling on the subsurface of Ge(111)

A new nonlinear magnetothermal effect, namely unidirectional Seebeck effect(USE), has been recently reported in magnetic/nonmagnetic topological insulator (TI) heterostructure and is ascribed to the asymmetry magnon scattering. Here, we show a new mechanism to generate USE in the Rashba two dimensional electron gas (2DEG), in which the magnetism or magnetic order is completely absent. It's found that the USE has a quantum origin from the spin-momentum locking generated from Rashba spin-orbit coupling. We find that the USE exhibits a sine dependence on the orientation of the magnetic field with respect to the direction of the temperature gradient and is dominant from the states above the Lifshitz point. The USE in the semiconductor Ge(111) subsurface states has been theoretically and systematically investigated.

Planktic foraminifera iodine/calcium ratios from plankton tows

Planktic foraminifera test iodine to calcium ratios represent an emerging proxy method to assess subsurface seawater oxygenation states. Several core-top studies show lower planktic foraminifera I/Ca in locations with oxygen depleted subsurface waters compared to well oxygenated environments. The reasoning behind this trend is that only the oxidized species of iodine, iodate, is incorporated in foraminiferal calcite. The I/Ca of foraminiferal calcite is thought to reflect iodate contents in seawater. To test this hypothesis, we compare planktic foraminifera I/Ca ratios, obtained from plankton tows, with published and new seawater iodate concentrations from 1) the Eastern North Pacific with extensive oxygen depletion, 2) the Benguela Current System with moderately depleted oxygen concentrations, and 3) the well oxygenated North and South Atlantic. We find the lowest I/Ca ratios (0.07 µmol/mol) in planktic foraminifera retrieved from the Eastern North Pacific, and higher values for samples (up to 0.72 µmol/mol) obtained from the Benguela Current System and North and South Atlantic. The I/Ca ratios of plankton tow foraminifera from environments with well oxygenated subsurface waters, however, are an order of magnitude lower compared to core-tops from similarly well-oxygenated regions. This would suggest that planktic foraminifera gain iodine post-mortem, either when sinking through the water column, or during burial.

Geoelectrical Sounding to Identify Sub-surface and Groundwater State at Village Banauli, Singrauli District, Madhya Pradesh, India

Electrical resistivity (Geoelectrical) methods are well-known and common techniques for investigating the groundwater potential zone. These methods are economically viable and have the highest resolving power compared with other geophysical methods. A total of fifteen Vertical electrical soundings were conducted in the village of Banauli, located in Singrauli district in Madhya Pradesh, India. Vertical electrical sounding was carried out using Schlumberger electrode configuration with the maximum current electrode (AB) spacing of 200 m and potential electrode (MN) spacing of 10 m. For interpretation of measured resistivity, the Partial curve matching technique was used to calculate the layer parameters (resistivity and thickness) and further depict the depth section of the profile. In this study, the maximum five-layer model is obtained, and most curves are of HAK types. The first layer has a mean resistivity value of 12.41 Ωm and a mean thickness of 0.94 m. The second layer has mean resistivity of 7.93 Ωm and a mean thickness of 4.79 m. The third layer has a mean thickness value of 10.55 m and a mean resistivity value of 16.54 Ωm. The fourth layer has a mean resistivity value of 20.17 Ωm and a mean thickness of 9.20 m, and finally, the fifth layer, the bedrock, has a higher mean resistivity value of 59.92 Ωm. Thus, the obtained results may be used for identifying the drilling site for the groundwater potential zone.

Explicit Analytic Solutions for the Subsurface Stress Field in Single Plane Contacts of Elastically Similar Truncated Cylinders or Wedges

As has been pointed out recently, a possible solution strategy to the wear–fatigue dilemma in fretting, operating on the level of contact mechanics and profile geometries, can be the introduction of “soft” sharp edges to the contact profiles, for example, by truncating an originally smooth profile. In that regard, analysis of possible mechanical failure of a structure, due to the contact interaction, requires the knowledge of the full subsurface stress state resulting from the contact loading. In the present manuscript, a closed-form exact solution for the subsurface stress state is given for the frictional contact of elastically similar truncated cylinders or wedges, within the framework of the half-plane approximation and a local-global Amontons–Coulomb friction law. Moreover, a fast and robust semi-analytical method, based on the appropriate superposition of solutions for parabolic contact, is proposed for the determination of the subsurface stress fields in frictional plane contacts with more complex profile geometries, and compared with the exact solution. Based on the analytical solution, periodic tangential loading of a truncated cylinder is considered in detail, and important scalar characteristics of the stress state, like the von-Mises equivalent stress, maximum shear stress, and the largest principal stress, are determined. Positive (i.e., tensile) principal stresses only exist in the vicinity of the contact edge, away from the pressure singularity at the edge of the profile, and away from the maxima of the von-Mises equivalent stress, or the maximum shear stress. Therefore, the fretting contact should not be prone to fatigue crack initiation.