SEARCH · Engineering Papers

Results for “surface wave”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. Finally, the patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Long‐Term Observations of Turbulence Vertical Velocity Spectra in a Convective Mixed Layer: Dependence on Land‐Surface Forcing in the U.S. Southern Great Plains

Doppler lidar vertical velocity retrievals were analyzed for the scale and structure of mixed-layer turbulence over a 7-year period, on fair-weather warm season days (442 cases) in the U.S. Southern Great Plains. Based on the hypothesized influence of land-surface forcing and updraft size on convective boundary layer clouds, a spectral analysis was performed to quantify effects of surface forcing (surface buoyancy flux, friction velocity, and evaporative fraction) on turbulence scale. Significant (order-of-magnitude) variations in spectral density were found in the energy-production subrange and mesoscale regimes. Integral scale decreased with increasing buoyancy flux (and Monin-Obukhov stability parameter), while spectral density in the energy-production subrange increased, implying a transition to buoyancy-driven cellular structures with narrower updrafts. The influence of stability parameter was limited to the neutral to convective transition, and could not explain the wide variation in spectral density in the mesoscale regime. However, high friction velocity (u* > 0.5 m s –1 ) was associated with larger integral scale (updraft width), and a greater portion of spectral density in the mesoscale regime. Lidar profiles and radar reflectivity for individual cases revealed evidence of roll and wave-like structures contributing to mesoscale variability on high friction velocity days, suggesting shear instabilities on days with wetter land surface conditions and smaller buoyancy flux. The role of friction velocity emphasizes the multivariate nature of surface influences on turbulence, and implies that variance-based turbulence closures may not be adequate for capturing effects of surface forcing on updraft size.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗

DOE OSTI2023

Three-dimensional Simulations of Massive Stars. II. Age Dependence

We present 3D full star simulations, reaching up to 90% of the total stellar radius, for three 7 M⊙ stars of different ages: zero-age main sequence (ZAMS), mid–main sequence (midMS), and terminal-age main sequence (TAMS). A comparison with several theoretical prescriptions shows that the generation spectra for all three ages are dominated by convective plumes. Two distinct overshooting layers are observed, with most plumes stopped within the layer situated directly above the convective boundary; overshooting to the second, deeper layer becomes progressively more infrequent with increasing stellar age. Internal gravity wave (IGW) propagation is significantly impacted in the midMS and TAMS models as a result of some IGWs getting trapped within their Brunt–Väisälä frequency spikes. A fundamental change in the wave structure across radius is also observed, driven by the effect of density stratification on IGW propagation causing waves to become evanescent within the radiative zone, with older stars being affected more strongly. We find that the steepness of the frequency spectrum at the surface increases from ZAMS to the older models, with older stars also showing more modes in their spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

General Analysis of Data Collected from DOE Lidar Buoy Deployments Off Virginia and New Jersey

Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinel lidar buoys for the U.S. Department of Energy’s Wind Energy Technologies. The purpose of these buoys is to collect hub-height winds and supporting meteorological and oceanographic information to facilitate the development of offshore wind energy in the U.S. In general, each buoy is deployed for a year or more at a given location in order to capture at least a full annual cycle of weather conditions. The initial deployment for one buoy was off the coast of Virginia beginning in 2015, and the other buoy was first deployed off the coast of New Jersey beginning in 2016. Over the last two years, PNNL has had an opportunity to analyze the data collected during these first two deployments. This report describes a substantial analysis of data collected by the two lidar buoys operated off the coasts of Virginia and New Jersey. The centerpiece instrument for each buoy as deployed off Virginia and New Jersey is a lidar system, which is designed to measure the horizontal wind vector from approximately 40 m to 200 m above the sea surface with vertical resolution of 40 m. Since the initial deployments, the original lidars have been replaced with more powerful Leosphere 866 v2 systems. The analyses in this report will apply to the original Vindicator systems. In addition to the wind profiles from the lidars, the buoys collect near-surface measurements of wind speed and direction, air temperature, relative humidity, barometric pressure, and solar irradiance. Oceanographic variables measured include the two-dimensional wave spectrum, water temperature and conductivity, and ocean current vectors to a depth of 90 m. An assessment of overall data recovery and a basic analysis of the data collected was provided in a previous report This report substantially extends that analysis. The various sections describe the development and application of an inertial measurement unit (IMU) data recovery scheme for the New Jersey deployment; a climatological analysis of winds at hub height and at the surface together with thermodynamic variables measured at the surface for the full deployment periods; an analysis of oceanographic observations describe sea state; the development of a refinement for NOAA’s WaveWatch III model to allow its application to near-shore areas; a basic climatology of ocean currents observed from the buoys; an analysis of observed winds in the framework of Monin-Obukhov Similarity Theory; and the development and evaluation of techniques to extract turbulence intensity and turbulence kinetic energy from the lidars. The analyses contained in this report provide a great deal of new information about offshore conditions on the U.S. East Coast. In addition, the experience gained will inform both configurations and analysis of data from future deployments of these lidar buoy systems.

17 WIND ENERGY↗

DOE OSTI2020

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗

DOE OSTI2024

Structural phase transition, s ± -wave pairing, and magnetic stripe order in bilayered superconductor La 3 Ni 2 O 7 under pressure

Motivated by the recently discovered high-T c superconductor La 3 Ni 2 O 7 , we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y 2- mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y 2- mode amplitude, the two phases may coexist in the range between 10.6 and 14 GPa, beyond which the Fmmm phase dominates. In addition, the magnetic stripe-type spin order with wavevector (π, 0) was stable at the intermediate region. Pairing is induced in the s ± -wave channel due to partial nesting between the M = (π, π) centered pockets and portions of the Fermi surface centered at the X = (π, 0) and Y = (0, π) points. This resembles results for iron-based superconductors but has a fundamental difference with iron pnictides and selenides. Moreover, our present efforts also suggest La 3 Ni 2 O 7 is qualitatively different from infinite-layer nickelates and cuprate superconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2022

The rapid destruction of toroidal magnetic surfaces

The operation of ITER will require reliable simulations in order to avoid major damage to the device from disruptions. Disruptions are the sudden breakup of magnetic surfaces across the plasma volume—a fast magnetic reconnection. This reconnection can be caused by the growth of perturbations outside of the plasma core causing an ideal perturbation to the core. This causes an increasing ratio of the maximum to the minimum separation, Δmax/Δmin, between neighboring magnetic surfaces. Magnetic reconnection becomes a dominant process when magnetic field lines can quickly interchange connections over a spatial scale ar. This occurs when Δmax/Δmin≳ar/Δd, where Δd is the scale over which non-ideal effects make magnetic field lines indistinguishable. Traditional reconnection theory is fundamentally different. It is a study of the steady-state cancelation of oppositely directed magnetic field components across a thin layer. During more than sixty years, mathematical implications of Faraday's law have been derived that clarify and constrain the physics of fast magnetic reconnection. These are reviewed because they are not commonly known but are needed to understand and to place in context how an ideal magnetic evolution can cause reconnection to quickly become a dominant process no matter how small Δd/ar may be.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Theoretical and global simulation analysis of collisional microtearing modes

Microtearing modes (MTMs) are suggested as a candidate for anomalous thermal transport in tokamak H-mode discharges. This study investigates MTMs in tokamak plasmas, employing simulations in the BOUT++ framework. It simplifies and linearizes the governing equations in detailed linear simulations. The study meticulously evaluates various conductivity models under diverse plasma conditions and collision regimes. The research thoroughly assesses different conductivity models across a range of plasma conditions and collision regimes. A unified dispersion relation that includes both MTM and Drift-Alfvén Wave (DAW) instabilities is derived, showing that DAW and MTM instabilities occur at varying distances from the rational surface. Specifically, MTMs become unstable near the rational surface but stabilize farther away, while drift-Alfvén instability appears farther from the rational surface. The study also re-derives MTM dispersion relations using Ohm's law and the vorticity equation, providing a thorough analysis of electromagnetic and electrostatic interactions in tokamaks. Global simulations demonstrate an inverse correlation between MTM growth rates and collisionality, and a direct correlation with temperature gradients. The nonalignment of the rational surface with the peak ω*e stabilizes the MTMs. Nonlinear simulations highlight electron temperature relaxation as the primary saturation mechanism for MTMs, with magnetic flutter identified as the dominant mode of electron thermal transport.

Fan, K. (ORCID:0000000227518809)↗

DOE OSTI2021

Multigap electron-phonon superconductivity in the quasi-one-dimensional pnictide K 2 Mo 3 As 3

In this work, we show, using density functional calculations, that quasi-one-dimensional K 2 Mo 3 As 3 , which is closely related to K 2 Cr 3 As 3 and has very similar superconducting properties, is not close to magnetism and has conventional s wave electron-phonon superconductivity. This superconductivity is of multigap character due to different coupling on different Fermi surface sheets. This is discussed in relation to the properties of this family of quasi-one-dimensional pnictide superconductors. The results show that this family of superconductors provides a unique opportunity for studying the interplay of spin fluctuations and electron-phonon superconductivity in transition metal pnictides, and offer a path for sorting out the different proposed superconducting scenarios in this fascinating family of pnictide superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Superconductivity at carrier density 10 17 cm – 3 in quasi-one-dimensional Li 0.9 Mo 6 O 17

Superconductivity is a fascinating phenomenon that involves an attractive interaction by which electrons are paired and can move without electrical resistance below a critical temperature T c . A paradigm for understanding superconductive pairing, the successful “BCS” theory, developed more than 65 years ago, is challenged by certain classes of materials for which parameters (e.g. the value of T c , the density of electrons, etc.) fall outside the range for which BCS theory is applicable. The study of such materials is important because new physics may be revealed and because superconductivity has many potential applications in sensing and electronics. Superconductivity at very low carrier density is a rare and fascinating phenomenon that is of considerable importance and interest in condensed matter physics. The conditions of high density of states and effective screening of Coulomb repulsion, both key ingredients of stable Cooper pairing in the BCS theory, may not be met at low carrier density. Thus, such systems may also harbor unconventional electron pairing (e.g. non-phononic) or pairing without phase coherence as a precursor to Bose-Einstein condensation. An intriguing route toward extreme low-density (and possibly unconventional) superconductivity may be found in quasi-one-dimensional metals for which density-wave instabilities compete for the ground state and can lead to gapping of substantial portions of the Fermi surface (FS), leaving a residual (reconstructed) FS comprising a very small carrier density. Our work highlights a novel material in this unusual category: Li 0.9 Mo 6 O 17 (T c =2 K) (also known as "lithium purple bronze" or "LiPB"). We present both electrical and thermoelectric magnetotransport measurements that reveal superconductivity to occur at carrier density 2 x10 17 cm -3 , among the lowest known for any superconductor, and four orders of magnitude lower than is expected from its chemical valence and band structure. With its exceptionally large ratio of superconducting to Fermi temperature, T c /TF ≳ 0.1, LiPB is properly placed at the top of a very small group of superconductors (including cuprates, heavy-fermion and Fe-based superconductors) for which this ratio falls in the BCS-BEC crossover regime and for which pairing is believed to be spin mediated.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Localized basis set for plutonium

The implementation of optimal strictly localized atomic orbitals basis for plutonium (Pu) using norm-conserving pseudopotential density-functional theory (DFT) is presented. The basis set was applied to the α, β, γ, δ, δ', and ε phases of Pu, δ-Pu surface, and δ-PuGa alloys. The computed properties of the Pu phases and δ-Pu surface were in good agreement with both available experimental data and prior DFT calculations based on plane-wave methodologies. Results for the δ-PuGa alloys were also in good agreement with experimental data. Finally, the reliability of the basis set was further demonstrated by using ab initio molecular dynamics to model the diffusion coefficient and activation barrier for atomic diffusion in a δ-PuGa alloy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Acoustic Travel-Time Tomography for Wind Energy

New remote sensing technologies for direct observation of the atmospheric flow in the proximity of operating wind turbines are necessary for the advancement of wind turbine design, wind plant control, and the validation of high-fidelity numerical models; fundamental goals of the DOE Wind Energy Technology Office's (WETO's) A2e program. Modern remote sensing technologies for commercial and research applications rely on the backscatter of light or sound waves, which limits spatial and temporal resolution and cannot produce reliable flow estimates near solid bodies or reflective surfaces. Acoustic tomography (AT) is an innovative sensing technology based on fundamental thermodynamic properties of the atmosphere and an inverse algorithm that can reconstruct both velocity and temperature fields within a network of acoustic transducers. AT relies on the direct travel path of acoustic signals between speakers and microphones in the network and can achieve much higher resolutions than conventional backscatter technologies. For wind energy R&D applications, AT represents an extremely economical option over other remote sensing technologies, integrating mature and commonly available acoustic and signal processing hardware. This report is intended as an update on the development of the AT array on the NREL Flatirons Campus (FC) and includes a brief review of the instrumentation, the calibration of equipment, signal processing, and resulting flow field estimates. Designs for future AT systems and applications are detailed in this report, including summaries of new capabilities, engineering and theoretical challenges, and anticipated costs help to outline future directions for AT technology R\&D. The potential for application of AT to utility-scale turbines is also discussed. Investing in the development of the proposed instrument will establish AT as a transformational RS technique that represents a significant advance in measurement instrumentation, enabling new fundamental insights into the structure of the atmospheric boundary layer and supporting applications for wind energy, pollution transport, and resource characterization. AT is also an extremely economical option over other RS instruments, integrating mature and commonly available acoustic and signal processing hardware.

17 WIND ENERGY↗

DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗

DOE OSTI2025

Uncertainty-Guided Prediction Horizon of Phase-Resolved Ocean Wave Forecasting Under Data Sparsity: Experimental and Numerical Evaluation

Accurate short-term wave forecasting is critical for the safe and efficient operation of marine structures that rely on real-time, phase-resolved ocean wave information for control and monitoring purposes (e.g., digital twins). These systems often depend on environmental sensors (e.g., waverider buoys, wave-sensing LIDAR). Challenges arise when upstream sensor data are missing, sparse, or phase-shifted due to drift. This study investigates the performance of two machine learning models, time-series dense encoder (TiDE) and long short-term memory (LSTM), for forecasting phase-resolved ocean surface elevations under varying degrees of data degradation. We introduce the τ-trimming algorithm, which adapts the prediction horizon based on uncertainty thresholds derived from historical forecasts. Numerical wave tank (NWT) and wave basin experiments are used to benchmark model performance under short- and long-term data masking, spatially coarse sensor grids, and upstream phase shifts. Results show under a 50% probability of upstream data loss, the τ-trimmed TiDE model achieves a 46% reduction in error at the most upstream target, compared to 22% for LSTM. Furthermore, phase misalignment in upstream data introduces a near-linear increase in forecast error. Under moderate model settings, a ±3 s misalignment increases the mean absolute error by approximately 0.5 m, while the same error is accumulated at ±4 s using the more conservative approach. These findings inform the design of resilient, uncertainty-aware wave forecasting systems suited for realistic offshore sensing environments.

42 ENGINEERING↗

DOE OSTI2025

TEAMER Technical Support on Biofouling and Corrosion for E-Wave Technologies LLC Novel (Abstract)

National Technology & Engineering Solutions of Sandia (NTESS) in collaboration with the Pacific Northwest National Laboratory (PNNL), and with guidance from E-Wave, will assess corrosion and biofouling challenges specific to E-Wave’s innovative WEC designed to retrofit into and power existing aquaculture farm infrastructure from ocean waves. The teams will examine coatings via a static raw seawater test and determine the efficacy of coatings for E-Wave’s applications. Samples will include panels of materials provided by E-Wave and coated with three different paint treatments. Fouling will be determined using a staining technique followed by photographic analysis to evaluate fouling intensity on the surface. Materials will also be weighed before and after submersion to facilitate weight change analysis resulting from fouling organisms aggregating on the materials.

16 TIDAL AND WAVE POWER↗

DOE OSTI2022

Utah FORGE - Development of a Reservoir Seismic Velocity Model and Seismic Resolution Study

This is data from and a final report on the development of a 3D velocity model for the larger FORGE area and on the seismic resolution in the stimulated fracture volume at the bottom of well 16A-32. The velocity model was developed using RMS velocities of the seismic reflection survey and seismic velocity logs from borehole measurements as an input model. To improve the accuracy of the model in the shallow subsurface, travel times phase arrivals of the direct propagating P-waves were determined from the seismic reflection data, using PhaseNet, a deep-neural-network-based seismic arrival time picking method. The travel times were subsequently inverted using the input velocity model. The seismic resolution study used borehole and surface seismic sensors as well as the seismicity observed during the April 2022 stimulation experiment to estimate the seismic resolution in the activated fracture reservoir. The data contain a 3D P- and S-wave velocity model for the larger FORGE area.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Seismic Elastic Double-Beam Characterization of Faults and Fractures for CO₂ Storage Site Selection

Site characterization for underground injection and storage of gigatonne-scale CO₂ requires reliable and cost-effective methods to detect and characterize faults and fractures and to assess their stress state and fault activation potential. This is critical, as wastewater injection and disposal have been shown to activate faults and induce earthquakes, and CO₂ leakage remains a key concern for long-term storage. In this project, we developed seismic methods to detect and characterize large-scale sedimentary and crystalline basement faults and associated small-scale fractures below conventional seismic imaging resolution using multicomponent (9C) surface seismic data. Machine learning was used to automatically interpret large-scale faults, providing key information for estimating the maximum magnitude of potential induced earthquakes. High-fidelity imaging was achieved by exploiting redundancy across multiple elastic wave modes, where independent images from different modes and frequencies cross-validate each other. We also used our nonlinear signal comparison (NLSC) method for ground roll removal, improving data quality in complex near-surface conditions. The methods were validated using field data acquired in central Montana. Results show that basement faults extend into the sedimentary section and that small-scale fractures are widespread above the basement. The inferred stress orientation is consistent with regional stress data, and the estimated maximum induced earthquake magnitude is small (Mw ~2.3). The developed workflow provides a practical approach for fault and fracture characterization and for assessing induced seismicity and leakage risk. It is directly applicable to CO₂ storage site selection and to other subsurface systems.

02 PETROLEUM↗