Subsurface Volume Research Articles

Mechanochemically responsive polymer enables shockwave visualization

Understanding the physical and chemical response of materials to impulsive deformation is crucial for applications ranging from soft robotic locomotion to space exploration to seismology. However, investigating material properties at extreme strain rates remains challenging due to temporal and spatial resolution limitations. Combining high-strain-rate testing with mechanochemistry encodes the molecular-level deformation within the material itself, thus enabling the direct quantification of the material response. Here, we demonstrate a mechanophore-functionalized block copolymer that self-reports energy dissipation mechanisms, such as bond rupture and acoustic wave dissipation, in response to high-strain-rate impacts. A microprojectile accelerated towards the polymer permanently deforms the material at a shallow depth. At intersonic velocities, the polymer reports significant subsurface energy absorption due to shockwave attenuation, a mechanism traditionally considered negligible compared to plasticity and not well explored in polymers. The acoustic wave velocity of the material is directly recovered from the mechanochemically-activated subsurface volume recorded in the material, which is validated by simulations, theory, and acoustic measurements. This integration of mechanochemistry with microballistic testing enables characterization of high-strain-rate mechanical properties and elucidates important insights applicable to nanomaterials, particle-reinforced composites, and biocompatible polymers.

A pore-scale study of fracture sealing through enzymatically-induced carbonate precipitation (EICP) method demonstrates its potential for CO2 storage management

Geological fractures are mechanical breaks in subsurface rock volumes that provide important subsurface flow pathways. However, the presence of fractures can cause unwanted challenges, such as gas leakage through fractured caprocks, which must be addressed. In this study, the dynamics of enzymatically induced carbonate precipitation in rock fractures and their subsequent influence on CO2 leakage were investigated from a pore-scale perspective for the first time. This was achieved through real-time monitoring of the injection of the solution into a rock-microfluidic flow cell using optical and scanning electron microscopy. It was revealed that the main growth dynamics occur during the first three injection cycles, with growth continuing until the fracture aperture is fully closed in the 6th cycle. Based on the flow simulation, a significant reduction of up to 25% in the CO2 conductivity of the original fracture is expected even after the first treatment cycle. Future studies are suggested to explore different resolutions, testing conditions, and to conduct 3-dimensional investigations of the growth dynamics.

Estimating differential penetration of green (532 nm) laser light over sea ice with NASA's Airborne Topographic Mapper: observations and models

Abstract. Differential penetration of green laser light into snow and ice has long been considered a possible cause of range and thus elevation bias in laser altimeters. Over snow, ice, and water, green photons can penetrate the surface and experience multiple scattering events in the subsurface volume before being scattered back to the surface and subsequently the instrument's detector, therefore biasing the range of the measurement. Newly formed sea ice adjacent to open-water leads provides an opportunity to identify differential penetration without the need for an absolute reference surface or dual-color lidar data. We use co-located, coincident high-resolution natural-color imagery and airborne lidar data to identify surface and ice types and evaluate elevation differences between those surfaces. The lidar data reveals that apparent elevations of thin ice and finger-rafted thin ice can be several tens of centimeters below the water surface of surrounding leads, but not over dry snow. These lower elevations coincide with broadening of the laser pulse, suggesting that subsurface volume scattering is causing the pulse broadening and elevation shift. To complement our analysis of pulse shapes and help interpret the physical mechanism behind the observed elevation biases, we match the waveform shapes with a model of scattering of light in snow and ice that predicts the shape of lidar waveforms reflecting from snow and ice surfaces based on the shape of the transmitted pulse, the surface roughness, and the optical scattering properties of the medium. We parameterize the scattering in our model based on the scattering length Lscat, the mean distance a photon travels between isotropic scattering events. The largest scattering lengths are found for thin ice that exhibits the largest negative elevation biases, where scattering lengths of several centimeters allow photons to build up considerable range biases over multiple scattering events, indicating that biased elevations exist in lower-level Airborne Topographic Mapper (ATM) data products. Preliminary analysis of ICESat-2 ATL10 data shows that a similar relationship between subsurface elevations (restored negative freeboard) and “pulse width” is present in ICESat-2 data over sea ice, suggesting that biased elevations caused by differential penetration likely also exist in lower-level ICESat-2 data products. The spatial correlation of observed differential penetration in ATM data with surface and ice type suggests that elevation biases could also have a seasonal component, increasing the challenge of applying a simple bias correction.

Long-term risk assessment of subsurface carbon storage: analogues, workflows and quantification

A key aspect of transferrable oil and gas expertise to subsurface CO 2 storage (SCS) relates to risk assessment. While initial subsurface risk and volume assessments for SCS projects are similar to oil and gas prospect evaluations, full life-cycle risk assessment requires evaluation of the current knowledge of the storage complex and future potential events that may have impacts over long time frames. It is important to learn from past water, gas, and CO 2 injection and storage projects, examples of which are reviewed here. The concepts of aleatory and epistemic risk and uncertainty are discussed, and the use of standard risk matrices for evaluation of long-term, low-frequency but potentially high-impact events is challenged. Drawing on the large volume of previous work, this paper highlights key elements for development of a robust risking framework, including thorough project framing and implementation of a staged approach to project execution. The paper assesses methods for ensuring all potential hazards are captured and addresses the challenges of defining a quantitative risking scale suitable for long-term SCS projects. Example quantitative risk profiles can be used to calculate the timing and duration of ‘Peak Risk’, augment monitoring and mitigation planning for management of risk and capital exposure, and help to ensure successful project outcomes. Thematic collection: This article is part of the Geoscience workflows for CO 2 storage collection available at: https://www.lyellcollection.org/topic/collections/geoscience-workflows-for-CO2-storage

How borehole wall elasticity affects wellbore storage capacity in periodic hydraulic tests

Characterizing fluid transport and pore pressure diffusion is key for monitoring and understanding many natural (e.g., seismically active zones and volcanic systems) and engineered environments (e.g., enhanced geothermal reservoirs and CO2 underground storage). Borehole hydraulic testing allows for inferring relevant properties of the probed subsurface volume, such as, for example, its transmissivity and diffusivity, assessing the governing flow regime as well as detecting the presence of hydraulic boundaries. Periodic hydraulic tests (PHT) achieve these objectives using an oscillatory fluid injection procedure while measuring the fluid pressure response in monitoring boreholes. The relevant information on the pressure diffusion process occurring in a probed formation is retrieved from the phase shifts and amplitude ratios between the flow rate and the interval pressure. In general, the interpretation of PHT data relies on the assumption that the pressure diffusion process is largely unaffected by the deformation of the probed rock volume. To assess this assumption, we present a conceptual PHT model based on the theory of poroelasticity to investigate the role played by hydromechanical coupling (HMC) effects during PHT and to assess whether and to what extent additional mechanical information can be extracted from these tests. These effects include the influence of the borehole wall deformation on the wellbore storage coefficient Sw, which quantifies the difference between the injected flow rates and those actually entering the porous formation. We show that, in homogeneous formations, the commonly taken approach based on the uncoupled diffusion solution for radial flow reproduces the results of our HMC model. However, neglecting the deformation effect of the borehole wall on the storage capacity of the system during injection can lead to underestimations of both transmissivity and diffusivity, particularly at shorter oscillation periods. We present analytical expressions for the borehole wall deformation contribution to Sw, which depends on the shear modulus and does not change with the oscillatory period. We show that it is possible to obtain the effective Sw along with the transmissivity and diffusivity using observations at various oscillatory periods. In the presence of hydraulic boundaries, such as no-flow or constant-pressure boundaries, at a given distance from the injection borehole, Sw becomes period-dependent, introducing apparent variations in the inferred shear modulus of the formation if the associated HMC effects are not accounted for. Our results thus highlight the importance of accounting for HMC effects in PHT data interpretation, not only to reliably and accurately characterize and monitor hydraulic properties but also to fully exploit the potential of PHT to provide additional information on the mechanical behavior of the probed rock volume.

Quantifying Subsurface Tar Volumes at Three Former Manufactured Gas Plants

Manufactured gas plants (MGPs) were a staple of the Industrial Revolution, operating in the US primarily on the East Coast, starting in the early 19th before phase-out in the mid-twentieth century. Their legacy includes the release of copious quantities of subsurface tar, and when located adjacent to waterways, tar in sediment due to direct discharge and leaking bulkheads. Computing the depth and extent of tar impacts is a three-dimensional problem required to facilitate appropriate remediation efforts. At three former MGPs, 306, 223, and 254 boreholes with 1452, 544, and 276 intervals were coded and %saturation attributed to saturated (Sn = 70%) and residual (S Nr = 20%) tar. These data were incorporated into the 3-D krieging software EnviroInsite to estimate the spatial extent of tar and oil impacts and their affiliated volumes. At the former Citizens, Fulton, and Metropolitan MGPs, the model estimated a most probable volume (MPV) of 7.21, 4.36, and 0.8 million gallons of tar. As a sensitivity analysis, Sn was decreased to 60%, reducing the MPV to 6.3, 3.4, and 0.6 million gallons, respectively, Regardless, it is irrefutable that a large amount of tar remains in the subsurface ∼70 years after the facilities ceased operating.

A rolling contact fatigue life prediction model for bearing steel considering its gradient structure due to cyclic hardening

This paper presents a rolling contact fatigue life prediction model for bearing steel. In the initial stage of rolling contact, the gradient structure appears in the subsurface region of bearing steel and the hardness shows a gradient distribution along the contact depth direction due to the roller compaction effect. This indicates that the fatigue resistance of bearing steel varies in the subsurface region. Besides, each subsurface material volume element is subject to different stress cycles. Based on Weibull theory, the survival possibility of subsurface volume elements is formulated. Then the rolling contact fatigue life is evaluated considering the stress state and anti-fatigue performance of each elementary volume of the subsurface material. According to the phenomenon and assumption, the model proposed for gradient material was applied in the rolling contact fatigue life prediction of the bearing steel GCr15 and validated with the fatigue experiment data in the open literature. Furthermore, the accuracy of the proposed model results was compared with the traditional empirical rolling contact fatigue life prediction models.

High-resolution geoelectrical characterization and monitoring of natural fluids emission systems to understand possible gas leakages from geological carbon storage reservoirs

Gas leakage from deep geologic storage formations to the Earth’s surface is one of the main hazards in geological carbon sequestration and storage. Permeable sediment covers together with natural pathways, such as faults and/or fracture systems, are the main factors controlling surface leakages. Therefore, the characterization of natural systems, where large amounts of natural gases are released, can be helpful for understanding the effects of potential gas leaks from carbon dioxide storage systems. In this framework, we propose a combined use of high-resolution geoelectrical investigations (i.e. resistivity tomography and self-potential surveys) for reconstructing shallow buried fracture networks in the caprock and detecting preferential gas migration pathways before it enters the atmosphere. Such methodologies appear to be among the most suitable for the research purposes because of the strong dependence of the electrical properties of water-bearing permeable rock, or unconsolidated materials, on many factors relevant to CO2 storage (i.e. porosity, fracturing, water saturation, etc.). The effectiveness of the suggested geoelectrical approach is tested in an area of natural gas degassing (mainly CH4) located in the active fault zone of the Bolle della Malvizza (Southern Apennines, Italy), which could represent a natural analogue of gas storage sites due to the significant thicknesses (hundreds of meters) of impermeable rock (caprock) that is generally required to prevent carbon dioxide stored at depth from rising to the surface. The obtained 3D geophysical model, validated by the good correlation with geochemical data acquired in the study area and the available geological information, provided a structural and physical characterization of the investigated subsurface volume. Moreover, the time variations of the observed geophysical parameters allowed the identification of possible migration pathways of fluids to the surface.

Roles of the Indo-Pacific subsurface Kelvin waves and volume transport in prolonging the triple-dip 2020–2023 La Niña

The rare triple-dip 2020–2023 La Niña event has resulted in a series of extreme climate events across the globe. Here, we reveal the role of tropical Indo-Pacific oceanic interactions in driving the first triple-dip La Niña of the twenty-first century. Specifically, we found that the eastern Indian Ocean subsurface warming anomalies were associated with the re-intensification of the subsequent La Niña event. The subsurface warming anomaly signals were propagated eastward by equatorial and coastal subsurface Kelvin waves from the eastern Indian Ocean to the western Pacific Ocean through the Indo-Pacific oceanic pathway, which contributes to the accumulation of heat content and deepens the thermocline in the western tropical Pacific. The westward Indonesian Throughflow (ITF) transported more heat during multi-year La Niña events from the western Pacific Ocean to the eastern Indian Ocean than during single-year events, resulting in the injection of more warm water into the eastern Indian Ocean. The combination of subsurface Kelvin wave propagation and increased ITF volume transport in the Indo-Pacific region acted to prolong the heat content in the western Pacific during the decay phase of La Niña, ultimately leading to the rare triple-dip 2020–2023 La Niña event.

Prediction of subsurface settlement induced by shield tunnelling in sandy cobble stratum

This study focuses on the analytical prediction of subsurface settlement induced by shield tunnelling in sandy cobble stratum considering the volumetric deformation modes of the soil above the tunnel crown. A series of numerical analyses is performed to examine the effects of cover depth ratio (C/D), tunnel volume loss rate (ηt) and volumetric block proportion (VBP) on the characteristics of subsurface settlement trough and soil volume loss. Considering the ground loss variation with depth, three modes are deduced from the volumetric deformation responses of the soil above the tunnel crown. Then, analytical solutions to predict subsurface settlement for each mode are presented using stochastic medium theory. The influences of C/D, ηt and VBP on the key parameters (i.e. B and N) in the analytical expressions are discussed to determine the fitting formulae of B and N. Finally, the proposed analytical solutions are validated by the comparisons with the results of model test and numerical simulation. Results show that the fitting formulae provide a convenient and reliable way to evaluate the key parameters. Besides, the analytical solutions are reasonable and available in predicting the subsurface settlement induced by shield tunnelling in sandy cobble stratum.

Surface microrelief induced by tillage management alters the pathway and composition of dissolved organic matter exports from soils to runoff during rainfall

Rainfall-runoff process mobilizes considerable dissolved organic matter (DOM) from soils to aquatic systems via surface and sub-surface flow pathways. Microrelief induced by tillage management can alter this flow partitioning and thus likely affect the associated pathway and composition of DOM exports during rainfall. This study conducted rainfall simulation experiments, combined with three-dimensional fluorescence spectra analysis and Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) technology, to explore the effects of different surface microreliefs treatments on the quantity and composition of DOM exports at a plot scale. Four typical microrelief treatments (i.e., contour tillage (CP), longitudinal tillage (LP), artificial digging tillage (AP), and flat tillage (CK, as control)) were selected. Results showed that ratios of surface to sub-surface flow volume under four treatments were in order of LP (193:1)>CK (73:1)>AP (4.5:1)>CP (0.5:1). DOM concentrations in sub-surface flow (13.37∼33.50 mg L-1) were 7∼18 times of that in surface flow (0.03∼4.56 mg L-1). The total export fluxes of DOM were 134, 139, 563, and 1214 mg m-2 at LP, CK, AP, and CP treatments, respectively, with proportions of 8%, 17%, 82%, and 98% via sub-surface flow. Compared to surface flow, DOM molecular composition in sub-surface flow showed a significant feature of higher oxygen to carbon ratio, higher molecular weight, and lower hydrogen to carbon ratio. The findings indicated that microrelief with higher surface storage capacity tends to favor a large flux of DOM export, primarily via sub-surface flow, which might significantly affect the DOM cycling in the receiving aquatic ecosystems.

An Integrated Geophysical and Unmanned Aerial Systems Surveys for Multi-Sensory, Multi-Scale and Multi-Resolution Cave Detection: The Gravaglione Site (Canale di Pirro Polje, Apulia)

Gravaglione represents one of the main swallow holes of the Canale di Pirro, low Murge, Apulia region, Italy. Here, after an intense rainstorm, a huge volume of rainwater accumulates at the surface. The drainage dynamics suggest that the Gravaglione could be part of a large, and potentially unknown, karst system. To verify this hypothesis and to acquire useful information on the possible karst environment features, an integrated aerial and geophysical multiscale and multimethod approach was applied. In particular, aerial photogrammetry, ground penetrating radar measurements and electrical resistivity tomography surveys were hence conducted and integrated to potentially detect the caves, define the subsurface volume possibly affected by karst systems and to verify the existence of links between the surficial morphology and the subsoil structure. The results provided interesting insights regarding the presence of a complex karst system extending up to 200 m b.g.l. and with a marked 3D nature. Overall, the Gravaglione case study demonstrates the geophysical approach validity and poses the basis for the development of an expeditive and low-cost high-resolution strategy for detecting and characterizing karst caves.

Reconnaissance Survey for Potential Energy Storage and Carbon Dioxide Storage Resources of Petroleum Reservoirs in Western Europe

Energy producers and utilities use oil and gas reservoirs for gas storage to meet peak seasonal demand or to supplement intermittent energy production. These reservoirs are also suitable for the long-term storage of carbon dioxide (CO2), a greenhouse gas. This study reports on a reconnaissance analysis of the potential magnitude of storage resources in 9424 known oil and gas reservoirs from 24 countries within highly industrialized western Europe. To standardize the storage resources of the oil and gas reservoirs, their volumetric capacity is expressed in terms of metric tons (mass) of CO2. Estimates of recoverable oil and gas at the surface are converted to subsurface volumes and then converted to the equivalent mass of CO2 at reservoir conditions. The results indicate 36.7 gigatons of CO2 could be stored, with oil reservoirs accounting for 32% of that total and natural gas reservoirs comprising the remaining 68%. About four-fifths of the reservoir storage resource is offshore, with about three-fourths of that offshore resource at water depths of 200 m or less. Most countries do not have the reservoir storage resources to store 15 years of CO2 at 2017 emission levels. With few exceptions the bulk of the storage is offshore for countries that do have at least 15 years of storage. The expansion of natural gas storage for strategic purposes in abandoned onshore gas reservoirs is not expected to seriously impact CO2 storage. The contribution of this analysis is the description of the spatial distribution of potential storage and physical accessibility.

Bayesian Monte Carlo Inversion of InSAR Time Series Deformation Induced by Wastewater Injection: A Case Study in West Texas

AbstractWastewater disposal can induce detectable surface uplift, which may cause ground instability and threaten infrastructure. The distributions of local hydro‐geomechanical parameters, especially Young's modulus and hydraulic conductivity, play an essential role in these geohazards. To constrain these parameters, we have inverted spatio‐temporal deformation measured by Interferometric Synthetic Aperture Radar (InSAR) and injection information using a Bayesian Monte Carlo approach with a poroelastic finite element model. Sentinel‐1A/B imagery from 2014 to 2020 is processed to track the spatio‐temporal deformation in Winkler county, West Texas, USA. The posterior distribution of subsurface effective volumes reveals under‐reported volumes in the well near the deformation center. In addition, the inversion results provide better constraints for the parameters than those solely obtained based on the cumulative spatial deformation or temporal development of the deformation center.

Mapping existing wellbore locations to compare technical risks between onshore and offshore CCS activities in Texas

AbstractCarbon dioxide capture and geologic storage (CCS; geologic sequestration) is a promising technology for reducing anthropogenic greenhouse gas emissions to the atmosphere from industrial point sources. Aspects of CCS have been investigated for over two decades, and many large‐ and small‐scale geologic storage field demonstration projects are now underway globally. Interest in offshore CCS has been increasing in recent years (e.g., European Union, Australia, Japan, and the United States). Deep geologic storage in offshore settings is analogous to onshore CCS activities in many respects (i.e., geologic and geotechnical aspects), but is distinct from previously explored seabed sediment CO2 storage) or deep marine dissolution). Given the large subsurface geologic storage volumes available in offshore settings, much discussion of offshore CCS is focused on the benefits and risks of such activity compared to onshore settings. Similar to onshore settings, existing (legacy) wells likely present the most direct migration pathway and largest risk of noncontainment in offshore settings. As part of current studies to evaluate geologic storage options in offshore settings along the Texas coast and greater Gulf of Mexico (GoM), mapping of the geographic distribution and ages of wells in a region containing coastal counties and extending 30 miles offshore Texas indicates that both well spatial density and well age decrease moving from onshore to offshore. Results suggest reduced risk of leakage owing to more rigorous and documented well completion and abandonment practices for these generally younger wells (although many are decades old). A result of decreased well density is that larger areas are available for leasing for CCS projects that avoid legacy wells altogether (> 1 mile from any existing well). The one‐mile designation is used as an arbitrary convention, and while it is recognized that this is smaller than a typical area of review (AoR) for permitting, each site will have a different AoR radius for consideration. The combination of large subsurface storage volumes under control of a single landowner and reduced risks from legacy wells makes offshore CCS attractive in the GoM. © 2023 Society of Chemical Industry and John Wiley & Sons, Ltd.

Self-ion irradiation effects on nanoindentation-induced plasticity of crystalline iron: A joint experimental and computational study

In this paper, experimental work is supported by multi-scale numerical modeling to investigate nanomechanical response of pristine and ion irradiated with Fe2＋ ions with energy 5 MeV high purity iron specimens by nanoindentation and Electron Backscatter Diffraction. The appearance of a sudden displacement burst that is observed during the loading process in the load–displacement curves is connected with increased shear stress in a small subsurface volume due to dislocation slip activation and mobilization of pre-existing dislocations by irradiation. The molecular dynamics (MD) and 3D-discrete dislocation dynamics (3D-DDD) simulations are applied to model geometrically necessary dislocations (GNDs) nucleation mechanisms at early stages of nanoindentation test; providing an insight to the mechanical response of the material and its plastic instability and are in a qualitative agreement with GNDs density mapping images. Finally, we noted that dislocations and defects nucleated are responsible the material hardness increase, as observed in recorded load–displacement curves and pop-ins analysis.

Study of GIS-based groundwater potential zones for agricultural sustainability in the arid region

Abstract The cluster-wise area of shallow and deep aquifer zones is used to estimate the potential of groundwater. The potential of the shallow aquifer zone is estimated at 4.61 MCM (million cubic meters) and for the deep aquifer zone at 17,509.03 MCM, while the total groundwater potential for both aquifer zones is estimated at 17,513.64 MCM. The Geographical Information System (GIS) was employed efficiently to estimate the subsurface volume of the lithological rock layers using cost-effective and time-saving techniques, while the Rockwork software integrated with GIS was successfully used to visualize the subsurface lithology and stratigraphy of the aquifer zones. The estimated potential of groundwater can be uncovered by using the alternative solar pumping system to improve the agricultural system in the study area, thereby reducing the migration rate, reducing poverty, and improving the socio-economic conditions of livelihood. In the future, too, it will be essential to design water quality studies to ensure the proper use of groundwater.

Local structure and magnetic properties of a nanocrystalline Mn-rich Cantor alloy thin film down to the atomic scale

The huge atomic heterogeneity of high-entropy materials along with a possibility to unravel the behavior of individual components at the atomic scale suggests a great promise in designing new compositionally complex systems with the desired multi-functionality. Herein, we apply multi-edge X-ray absorption spectroscopy (extended X-ray absorption fine structure (EXAFS), X-ray absorption near edge structure (XANES), and X-ray magnetic circular dichroism (XMCD)) to probe the structural, electronic, and magnetic properties of all individual constituents in the single-phase face-centered cubic (fcc)-structured nanocrystalline thin film of Cr20Mn26Fe18Co19Ni17 (at.%) high-entropy alloy on the local scale. The local crystallographic ordering and component-dependent lattice displacements were explored within the reverse Monte Carlo approach applied to EXAFS spectra collected at the K absorption edges of several constituents at room temperature. A homogeneous short-range fcc atomic environment around the absorbers of each type with very similar statistically averaged interatomic distances (2.54–2.55 Å) to their nearest-neighbors and enlarged structural relaxations of Cr atoms were revealed. XANES and XMCD spectra collected at the L2,3 absorption edges of all principal components at low temperature from the oxidized and in situ cleaned surfaces were used to probe the oxidation states, the changes in the electronic structure, and magnetic behavior of all constituents at the surface and in the sub-surface volume of the film. The spin and orbital magnetic moments of Fe, Co, and Ni components were quantitatively evaluated. The presence of magnetic phase transitions and the co-existence of different magnetic phases were uncovered by conventional magnetometry in a broad temperature range.

Morphometric constraints on the formation of new terrestrial analogs for planetary pits

The origin of geological depressions abounding on Mars and other planetary bodies remains poorly understood, partially due to the limited variability in the geological settings of existing terrestrial analogs. Here, we present a new terrestrial analog that is located at the northwestern margin of the Levantine volcanic field of Harrat Ash−Shaam along the Dead Sea Transform. The analog site consists of tens of geological depressions (locally named “juba”) that morphologically resemble Martian bowl-shaped pits and occur within a Pleistocene basaltic plateau that overlies Meso-Cenozoic carbonates. To constrain plausible formation mechanisms for the juba depressions, we carried out detailed field mapping and morphometric analyses using a 0.25 m/pixel digital terrain model (DTM) derived from airborne light detection and ranging (LiDAR) survey covering 34 km2 of the study area, and centimeter-scale, ground-based LiDAR scans of selected juba depressions. We show that variable magnitudes of slope asymmetry between north- and south-facing walls within the juba depressions, along with different degrees of sediment infilling, provide effective proxies for the relative geomorphic maturity of these landforms, and in turn indicate asynchronous formation of the juba depressions after the Pleistocene emplacement of the Harrat Ash−Shaam basalts in the study area. Our findings preclude formation of the juba depressions by phreatomagmatic explosions and instead point toward collapse into missing subsurface volume. In a broader context, we propose that the morphometric analyses developed herein to distinguish between plausible juba formation mechanisms in the Harrat Ash−Shaam volcanic field can be extended to better constrain the formation mechanisms of similar pit features on Mars and other planetary bodies.

Volume Contraction in Shallow Sediments: Discrete Element Simulation

Displacements induced by mineral dissolution and subsurface volume contraction affect overlying soils. In this study, we examine the consequences of mass loss or volume contraction at shallow depths using a discrete element method. The goal of the study is to identify particle-scale and global effects as a function of the relative depth of a dissolving inclusion, initial soil density, and granular interlocking. There are successive arch formation and collapse events, and a porosity front propagates upwards as grains slide down to refill the space. Grains around and within the refilled cavity are loosely packed and have small contact forces that are sufficient to avert the buckling of granular arches that form around the dissolving zone. Denser packings and interlocking combine to exacerbate rotational frustration and lead to more pronounced force chains along granular arches, looser fill, and reduced surface settlement. In fact, surface settlement vanishes, and the sediment hides the localized dissolution when deep inclusions z/D ≥ 5 dissolve within dense sediments. While scaling relations limit the extrapolation of these numerical results to tunneling and mining applications, macroscale trends observed in the field resemble results gathered in this study.