Results for “shallow” — Page 57

DOE OSTI2022

Donor doping of corundum (Al x Ga 1- x ) 2 O 3

We report corundum (Al x Ga 1-x ) 2 O 3 alloys have been proposed as a candidate ultrawide-bandgap oxide for a number of applications, but doping is unexplored. We examine the prospects for n-type doping with H, Si, Ge, Sn, Hf, Zr, and Ta in corundum (Al x Ga 1-x ) 2 O 3 alloys using first-principles calculations. All of the dopants are shallow donors in corundum Ga 2 O 3 . In the (Al x Ga 1-x ) 2 O 3 alloy, they transition from shallow to deep donors at Al compositions that are unique to each donor. Si and Hf remain shallow donors up to the highest Al contents in corundum (Al x Ga 1-x ) 2 O 3 alloys and are still shallow even as the (Al x Ga 1-x ) 2 O 3 bandgap exceeds 6.5 eV. Finally, we address the detrimental role of cation vacancies as compensating deep acceptors and suggest that doping in a hydrogen-rich environment under cation-rich conditions can be used to overcome this problem.

36 MATERIALS SCIENCE↗

DOE OSTI2019

The hunt for the third acceptor in CuInSe 2 and Cu(In,Ga)Se 2 absorber layers

We show the model for intrinsic defects in Cu(In,Ga)Se 2 semiconductor layers is still under debate for the full range between CuInSe 2 and CuGaSe 2 . It is commonly agreed by theory and experiment, that there are at least one shallow donor and two shallow acceptors. Spatially resolved photoluminescence on CuGaSe 2 previously revealed a third acceptor. In this study we show with the same method that the photoluminescence peak at 0.94 eV in CuInSe 2 , previously attributed to a third acceptor, is a phonon replica. However another pronounced peak at 0.9 eV is detected on polycrystalline CuInSe 2 samples grown with high copper and selenium excess. Intensity and temperature dependent photoluminescence measurements reveal that this peak originates from a DA-transition from a shallow donor (<8 meV) into a shallow acceptor A3 (135 [Formula: see text] 10) meV. The DA3 transition has three distinct phonon replicas with 28 meV spectral spacing and a Huang Rhys factor of 0.75. Complementary admittance measurements are dominated by one main step with an activation energy of 125 meV which corresponds well with the found A3 defect. The same defect is also observed in Cu(In,Ga)Se 2 samples with low gallium content. For [Ga]/([Ga] + [In])-ratios of up to 0.15 both methods show a concordant increase of the activation energy with increasing gallium content shifting the defect deeper into the bandgap. The indium vacancy ν In is discussed as a possible origin of the third acceptor level in CuInSe 2 and ν III in Cu(In,Ga)Se 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Role of carbon and hydrogen in limiting n -type doping of monoclinic ( Al x Ga 1 – x ) 2 O 3

In this study, we use hybrid density functional calculations to assess n -type doping in monoclinic ( Al x Ga 1 – x ) 2 O 3 alloys. We focus on silicon, the most promising donor dopant, and study the structural properties, formation energies, and charge-state transition levels of its various configurations. We also explore the impact of carbon and hydrogen, which are common impurities in metal-organic chemical vapor deposition (MOCVD). In Ga 2 O 3 , Si Ga is an effective shallow donor, but in Al 2 O 3 Si Al acts as a DX center with a ( + / – ) transition level in the band gap. Interstitial hydrogen acts as a shallow donor in Ga 2 O 3 but behaves as a compensating acceptor in n -type Al 2 O 3 . Interpolation indicates that Si is an effective donor in ( Al x Ga 1 – x ) 2 O 3 up to 70% Al, but it can be compensated by hydrogen already at 1% Al. We also assess the diffusivity of hydrogen and study complex formation. Si cation – H complexes have relatively low binding energies. Substitutional carbon on a cation site acts as a shallow donor in Ga 2 O 3 , but can be stable in a negative charge state in ( Al x Ga 1 – x ) 2 O 3 when x > 5 % . Substitutional carbon on an oxygen site ( C O ) always acts as an acceptor in n -type ( Al x Ga 1 – x ) 2 O 3 , but will incorporate only under relatively oxygen-poor conditions. C O – H complexes can actually incorporate more easily, explaining observations of carbon-related compensation in Ga 2 O 3 grown by MOCVD. We also investigate C cation – H complexes, finding they have high binding energies and act as compensating acceptors when x > 56 % ; otherwise the hydrogen just passivates the unintentional carbon donors. C-H complex formation explains why MOCVD-grown Ga 2 O 3 can exhibit record-low free-carrier concentrations, in spite of the unavoidable incorporation of carbon. Our study highlights that, while Si is in principle a suitable shallow donor in ( Al x Ga 1 – x ) 2 O 3 alloys up to high Al compositions, control of unintentional impurities is essential to avoid compensation.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Strong Landau-quantization effects in high-magnetic-field superconductivity of a two-dimensional multiple-band metal near the Lifshitz transition

We investigate the onset of superconductivity in a magnetic field for a clean two-dimensional multiple-band superconductor in the vicinity of the Lifshitz transition when one of the bands is very shallow. Due to the small number of carriers in this band, the quasiclassical Werthamer-Helfand approximation breaks down and Landau quantization has to be taken into account. We found that the transition temperature T C2 ( H ) has giant oscillations and is resonantly enhanced at the magnetic fields corresponding to the matching of the chemical potential with the Landau levels in the shallow band. This enhancement is especially pronounced for the lowest Landau level. As a consequence, the reentrant superconducting regions in the temperature-field phase diagram emerge at low temperatures near the magnetic fields at which the shallow-band Landau levels cross the chemical potential. The specific behavior depends on the relative strength of the intraband and interband pairing interactions and the reentrance is most pronounced in the purely interband coupling scenario. The reentrant behavior is suppressed by the Zeeman spin splitting in the shallow band; the separated regions disappear already for very small spin-splitting factors. On the other hand, the reentrance is restored in the resonance cases when the spin-splitting energy exactly matches the separation between the Landau levels. As a result, the predicted behavior may be realized in the gate-tuned FeSe monolayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2018

Interplay between orbital-quantization effects and the Fulde-Ferrell-Larkin-Ovchinnikov instability in multiple-band layered superconductors

We explore superconducting instability for a clean two-band layered superconductor with deep and shallow bands in the magnetic field applied perpendicular to the layers. In the shallow band, the quasiclassical approximation is not applicable, and Landau quantization has to be accounted for exactly. The electronic spectrum of this band in the magnetic field is composed of the onedimensional Landau-level minibands. With increasing magnetic field the system experiences series of Lifshitz transitions when the chemical potential enters and exits the minibands. These transitions profoundly in uence the shape of the upper critical field at low temperatures. In addition, the Zeeman spin splitting may cause the nonuniform state with interlayer modulation of the superconducting order parameter (Fulde-Ferrell-Larkin-Ovchinnikov state). Typically, the quantization effects in the shallow band strongly promote the formation of this state. The uniform state remains favorable only in the exceptional resonance cases when the spin-splitting energy exactly matches the Landau-level spacing. Furthermore, for specific relations between electronic spectrum parameters, the alternating FFLO state may realize, in which the order parameter changes sign between the neighboring layers. For all above cases, the reentrant high-field superconducting states may emerge at low temperatures if the shallow band has significant contribution to the Cooper pairing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Drainage enhances modern soil carbon contribution but reduces old soil carbon contribution to ecosystem respiration in tundra ecosystems

Abstract Warming temperatures are likely to accelerate permafrost thaw in the Arctic, potentially leading to the release of old carbon previously stored in deep frozen soil layers. Deeper thaw depths in combination with geomorphological changes due to the loss of ice structures in permafrost, may modify soil water distribution, creating wetter or drier soil conditions. Previous studies revealed higher ecosystem respiration rates under drier conditions, and this study investigated the cause of the increased ecosystem respiration rates using radiocarbon signatures of respired CO 2 from two drying manipulation experiments: one in moist and the other in wet tundra. We demonstrate that higher contributions of CO 2 from shallow soil layers (0–15 cm; modern soil carbon) drive the increased ecosystem respiration rates, while contributions from deeper soil (below 15 cm from surface and down to the permafrost table; old soil carbon) decreased. These changes can be attributed to more aerobic conditions in shallow soil layers, but also the soil temperature increases in shallow layers but decreases in deep layers, due to the altered thermal properties of organic soils. Decreased abundance of aerenchymatous plant species following drainage in wet tundra reduced old carbon release but increased aboveground plant biomass elevated contributions of autotrophic respiration to ecosystem respiration. The results of this study suggest that drier soils following drainage may accelerate decomposition of modern soil carbon in shallow layers but slow down decomposition of old soil carbon in deep layers, which may offset some of the old soil carbon loss from thawing permafrost.

Kwon, Min Jung↗

DOE OSTI2018

Early burial mud diapirism and its impact on stratigraphic architecture in the Carboniferous of the Shannon Basin, County Clare, Ireland

Delta fronts are often characterized by high rates of sediment supply that result in unstable slopes and a wide variety of soft-sediment deformation, including the formation of overpressured and mobile muds that may flow plastically during early burial, potentially forming mud diapirs. The coastal cliffs of County Clare, western Ireland, expose Pennsylvanian (Namurian) delta-front deposits of the Shannon Basin at large scale and in three dimensions. These deposits include decametre-scale, internally chaotic mudstone masses that clearly impact the surrounding sedimentary strata. Evidence indicates that these were true mud (unlithified sediment) diapirs that pierced overlying strata. This study documents a well-exposed ca 20 m tall mud diapir and its impact on the surrounding mouth-bar deposits of the Tullig Cyclothem. A synsedimentary fault and associated rollover dome, evident from stratal thicknesses and the dip of the beds, define one edge of the diapir. These features are interpreted as recording the reactive rise of the mud diapir in response to extensional faulting along its margin. Above the diapir, heterolithic sandstones and siltstones contain evidence for the creation of localized accommodation, suggesting synsedimentary filling, tilting and erosion of a shallow sag basin accommodated by the progressive collapse of the diapir. Two other diapirs are investigated using three-dimensional models built from ‘structure from motion’ drone imagery. Both diapirs are interpreted to have grown predominantly through passive rise (downbuilding). Stratal relationships for all three diapirs indicate that they were uncompacted and fluid-rich mud beds that became mobilized through soft-sediment deformation during early burial (i.e. <50 m, likely <10 m depth). Each diapir locally controlled the stratigraphic architecture in the shallow subsurface and potentially influenced local palaeocurrents on the delta. The mud diapirs studied herein are distinct from deeper ‘shale diapirs’ that have been inferred from seismic sections worldwide, now largely disputed.

defects, mechanical behavior, carbon sequestration↗

DOE OSTI2018

Interaction Rates of Group-III and Group-V Impurities with Intrinsic Point Defects in Irradiated Si and Ge

A comparative study of interactions of shallow impurities with primary defects in oxygen- and carbon-lean moderately doped Si and Ge subjected to irradiation with 0.9 MeV electrons, {sup 60}Co gamma-rays, and 15 MeV protons at room temperature is presented and discussed. For the quantitative characterization of such interactions, changes in the total concentration of the original shallow group-V donor or group-III acceptor impurities in the irradiated materials are determined by Hall effect measurements over a wide temperature range. Losses of the shallow donor or acceptor states in the irradiated Si and Ge are indicative of their removal rates that can be used for estimation of production rates of primary defects interacting with the dopants. Some important factors affecting the interactions between primary defects and shallow impurities in Si and Ge are highlighted.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Observed Covariations in Boundary Layer and Cumulus Cloud Layer Processes

In this paper we examine variations in boundary-layer processes spanning the shallow to deep cumulus transition. This is accomplished by differentiating boundary layer properties on the basis of convective outcomes, ranging from shallow to deep, as observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, USA. Doppler lidar, radar, and radiosonde data are combined to determine statistical differences in boundary layer and cloud layer properties using a large sample (236) of days with a range of convective outcomes: shallow, congestus, and deep convection. In these analyses, the radar characterizes diurnal cloud depth, the lidar quantifies updraft and downdraft properties in the subcloud layer, and daily radiosonde data provides the convective inhibition (CIN). Combined, these data are used to test the hypothesis that deep convection occurs when the strength of the boundary layer turbulence (i.e., TKE) exceeds the strength of the energy barrier (i.e., CIN) at the top of the CBL. Results show that days with deep convective clouds have significantly lower vertical velocity variance and weaker updrafts within the subcloud layer. However, CIN values are also found to be significantly lower on deep convective days, allowing for these weaker updrafts to penetrate the energy barrier and reach the level of free convection (LFC). In contrast, shallow convective outcomes occur when the updrafts are strong in an absolute sense, but are weak when compared to the strength of the energy barrier. These findings support the use of the CIN/TKE framework in parameterizing convection in coarse resolution models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

An Improved Double-Gaussian Closure for the Subgrid Vertical Velocity Probability Distribution Function

Abstract The vertical velocity probability distribution function (PDF) is analyzed throughout the depth of the lower atmosphere, including the subcloud and cloud layers, in four large-eddy simulation (LES) cases of shallow cumulus and stratocumulus. Double-Gaussian PDF closures are examined to test their ability to represent a wide range of turbulence statistics, from stratocumulus cloud layers characterized by Gaussian turbulence to shallow cumulus cloud layers displaying strongly non-Gaussian turbulence statistics. While the majority of the model closures are found to perform well in the former case, the latter presents a considerable challenge. A new model closure is suggested that accounts for high skewness and kurtosis seen in shallow cumulus cloud layers. The well-established parabolic relationship between skewness and kurtosis is examined, with results in agreement with previous studies for the subcloud layer. In cumulus cloud layers, however, a modified relationship is necessary to improve performance. The new closure significantly improves the estimation of the vertical velocity PDF for shallow cumulus cloud layers, in addition to performing well for stratocumulus. In particular, the long updraft tail representing the bulk of cloudy points is much better represented and higher-order moments diagnosed from the PDF are also greatly improved. However, some deficiencies remain owing to fundamental limitations of representing highly non-Gaussian turbulence statistics with a double-Gaussian PDF.

Fitch, A. C.↗

DOE OSTI2016

Investigation of the Residual in Column-Integrated Atmospheric Energy Balance Using Cloud Objects

Observationally based atmospheric energy balance is analyzed using Clouds and the Earth’s Radiant Energy System (CERES)-derived TOA and surface irradiance, Global Precipitation Climatology Project (GPCP)-derived precipitation, dry static and kinetic energy tendency and divergence estimated from ERA-Interim, and surface sensible heat flux from SeaFlux. The residual tends to be negative over the tropics and positive over midlatitudes. A negative residual implies that the precipitation rate is too small, divergence is too large, or radiative cooling is too large. The residual of atmospheric energy is spatially and temporally correlated with cloud objects to identify cloud types associated with the residual. Spatially, shallow cumulus, cirrostratus, and deep convective cloud-object occurrence are positively correlated with the absolute value of the residual. The temporal correlation coefficient between the number of deep convective cloud objects and individual energy components, net atmospheric irradiance, precipitation rate, and the sum of dry static and kinetic energy divergence and their tendency over the western Pacific are 0.84, 0.95, and 0.93, respectively. However, when all energy components are added, the atmospheric energy residual over the tropical Pacific is temporally correlated well with the number of shallow cumulus cloud objects over tropical Pacific. Because shallow cumulus alters not enough atmospheric energy compared to the residual, this suggests the following: 1) if retrieval errors associated with deep convective clouds are causing the column-integrated atmospheric energy residual, the errors vary among individual deep convective clouds, and 2) it is possible that the residual is associated with processes in which shallow cumulus clouds affect deep convective clouds and hence atmospheric energy budget over the tropical western Pacific.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Rayleigh Wave Propagation in the Bighorn Mountains Region, Wyoming

Short-period Rayleigh waves, Rg , provide strong constraints on the depth of shallow seismic events and are of interest for monitoring small explosions. Characterizing the seismic sources that generate Rg requires an understanding of how shallow crustal structure affects Rayleigh wave propagation. Here, in support of these efforts, this study utilizes observed waveforms from small shallow explosions recorded on temporary seismic network deployments in the Bighorn region, Wyoming. We study regional near-surface geology by measuring changes in surface-wave amplitude and polarization during propagation through basins, foothills, and mountains. We develop additional insight by carrying out surface-wave eigenfunction analyses and numerical-wave simulations, which together reproduce many characteristics seen in the observed waveforms. Our results show how sedimentary basins in the Bighorn region allow for amplified prograde-polarized higher-mode and retrograde-polarized fundamental-mode Rayleigh waves, whereas adjacent mountains only support retrograde motion. These different modes provide distinct constraints on the Earth structure and source characteristics, potentially enabling targeted inversions in future studies. Finally, our findings provide insight into Rg propagation through complex near-surface geology, improving our understanding of shallow propagation and source effects that are relevant to explosion monitoring efforts.

58 GEOSCIENCES↗

DOE OSTI2018

Control of Basin Water Depth On Channel Morphology and Autogenic Timescales in Deltaic Systems

River channel geometry is often controlled by upstream boundary conditions, including fluvial discharge and sediment properties. At the coast, downstream boundary conditions (e.g., tides, waves, water depth) also strongly influence channel formation and evolution. Here, we conducted a set of experiments to determine the effects of basin water depth (i.e., a downstream boundary condition) on the evolution and geometry of fluviodeltaic channels and lobes. Internal dynamics (autogenic processes) in the fluviodeltaic system drive channel avulsion through cycles of sediment storage and release. Experimental results indicate an increase in the timescale of autogenic storage and release with increasing basin water depth. Deeper basin water requires a larger volume to be filled within the delta front, thus more time to complete one autogenic storage and release cycle for a given sediment discharge. While a relationship between delta-front volume and autogenic storage and release timescales is expected, we show that autogenically generated morphological changes in the delta topset and distributary channels also exert control on timescales of storage and release. Deltas building into deeper basins develop steeper topsets, and deeper distributary channels that cause high-magnitude topset slope fluctuations, which contribute to the long autogenic timescales. Deposits in shallow basins exhibit both shallower topset slopes and shallower channels. Channel bed slopes are similar (~0.06) across all experiments, but lateral channel migration rates varied with basin depth. Deltas building into shallow basins had rapid lateral channel migration, such that channels quickly reworked the delta topset. For deep basins, channel migration rates were much slower, so the topset was reworked less often, allowing the topset to build to steeper slopes before being reworked. These experiments indicate an intimate relationship between lateral channel migration and topset aggradation. In addition, the deeper and more stable channels in deeper basins generally developed a wider range of channel widths, some of which produced elongate lobes. Finally, we conclude that the downstream boundary has a strong control on fluviodeltaic morphology, which can result in a striking difference in the autogenic timescale.

58 GEOSCIENCES↗

DOE OSTI2016

Deep Vadose Zone Respiration Contributions to Carbon Dioxide Fluxes from a Semiarid Floodplain

Core Ideas A significant fraction of CO 2 fluxes from a semiarid floodplain originates from below the root zone. Measured and calculated field CO 2 fluxes are consistent. Laboratory‐measured respiration rates are consistent with field results. Although CO 2 fluxes from soils are often assumed to originate within shallow soil horizons (<1‐m depth), relatively little is known about respiration rates at greater depths. We compared measured and calculated CO 2 fluxes at the Rifle floodplain along the Colorado River and measured CO 2 production rates of floodplain sediments to determine the relative importance of deeper vadose zone respiration. Calculations based on measured CO 2 gradients and estimated effective diffusion coefficients yielded fluxes that are generally consistent with measurements obtained at the soil surface (326 g C m −2 yr −1 ). Carbon dioxide production from the 2.0‐ to 3.5‐m depth interval was calculated to contribute 17% of the total floodplain respiration, with rates that were larger than some parts of the shallower vadose zone and underlying aquifer. Microbial respiration rates determined from laboratory incubation tests of the sediments support this conclusion. The deeper unsaturated zone typically maintains intermediate water and air saturations, lacks extreme temperatures and salinities, and is annually resupplied with organic carbon from snowmelt‐driven recharge and by water table decline. This combination of favorable conditions supports deeper unsaturated zone microbial respiration throughout the year.

Tokunaga, Tetsu K.↗

DOE OSTI2013

Groundwater monitoring program plan and conceptual site model for the Al-Tuwaitha Nuclear Research Center in Iraq.

The Radiation Protection Center of the Iraqi Ministry of Environment is developing a groundwater monitoring program (GMP) for the Al-Tuwaitha Nuclear Research Center located near Baghdad, Iraq. The Al-Tuwaitha Nuclear Research Center was established in about 1960 and is currently being cleaned-up and decommissioned by Iraq s Ministry of Science and Technology. This Groundwater Monitoring Program Plan (GMPP) and Conceptual Site Model (CSM) support the Radiation Protection Center by providing: A CSM describing the hydrogeologic regime and contaminant issues, recommendations for future groundwater characterization activities, and descriptions of the organizational elements of a groundwater monitoring program. The Conceptual Site Model identifies a number of potential sources of groundwater contamination at Al-Tuwaitha. The model also identifies two water-bearing zones (a shallow groundwater zone and a regional aquifer). The depth to the shallow groundwater zone varies from approximately 7 to 10 meters (m) across the facility. The shallow groundwater zone is composed of a layer of silty sand and fine sand that does not extend laterally across the entire facility. An approximately 4-m thick layer of clay underlies the shallow groundwater zone. The depth to the regional aquifer varies from approximately 14 to 17 m across the facility. The regional aquifer is composed of interfingering layers of silty sand, fine-grained sand, and medium-grained sand. Based on the limited analyses described in this report, there is no severe contamination of the groundwater at Al-Tuwaitha with radioactive constituents. However, significant data gaps exist and this plan recommends the installation of additional groundwater monitoring wells and conducting additional types of radiological and chemical analyses.

Copland, John Robin↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2023

Recent Developments in Hydrogeologic Applications for Strain Tensor Analyses

Changes in fluid pressure deform porous media and this effect occurs in a variety of hydrogeologic processes, from the change in storage during pumping or injection to fluctuations in water levels caused by barometric pressure. We have developed instruments for measuring small strains in porous media, and we have used the resulting strain data to evaluate well testing, hydraulic fracturing, manual loading at the ground surface, and ambient hydrologic processes, like rainfall and evaporation. A particularly important application is the use of strain tensor data measured at shallow depths to analyze well tests or hydraulic fractures conducted at much greater depths. An early demonstration of this technique was conducted at the North Avant Field north of Tulsa, Oklahoma, where Pennsylvanian sandstone creates a confined aquifer and oil reservoir at a depth of 530m. We have showed that the strains caused by injecting into the aquifer could be measured at a depth of 30m and used to evaluate the properties of the aquifer. We recently expanded the array of strainmeters at the North Avant Field by deploying three more instruments at shallow depth (30m) along with a deep instrument at 520m depth in the winter, 2021. The deep instrument is deployed in shale caprock slightly above the aquifer. To our knowledge, the deep strainmeter at the North Avant Field is the deepest strainmeter ever deployed and this required refining methods originally developed for shallow deployments. The instrument was lowered to depth on oil field tubing and cemented in place using techniques and materials developed for use in oil wells. Optical fiber used to communicate with the instrument was cemented in the annulus on the outside of tubing. This is significant because the techniques we used could readily be extended to greater depth, suggesting that strainmeters can be deployed over a wide range of depths for monitoring critical subsurface processes. For example, it suggests that strainmeters could be deployed through the caprock to monitor for leaks from underlying CO2 storage reservoirs. The strainmeter array at the North Avant Field has been used to characterize deformation during a series of injection tests in the spring and summer, 2021. All the new strainmeters respond to pumping, and the strainmeters we installed earlier also responded. To our knowledge, this is the first application of well testing in a deep aquifer that was monitored by an array of strainmeters—our earlier work used strainmeters at a single location. We are currently analyzing the strain data using an analytical solution, a proxy-based Bayesian inversion algorithm, and other methods. Strainmeter data has also been used to characterize periodic pumping tests by us and Riley Blais. A periodically varying pumping rate causes both hydraulic head and strain signals that vary with the same period as the pumping. The peaks and troughs of the head in monitoring wells lag behind the peaks and troughs of the head in the pumping well, and this lag time increases with distance from the pumping well. The lag time of the pressure and the distance to the monitoring well can be used in a simple analysis to estimate the hydraulic diffusivity of the aquifer. The lag time determined from strain data can be used to estimate aquifer properties using the same analysis that works for the pressure only for strain data measured at particular locations. That is because the strain field in a confining unit advances upward, laterally and then downward even though the pressure in the underlying aquifer only advances laterally, according to our recent simulations. We have field data showing that a small periodic signal superimposed on an injection rate at the North Avant Field will create a periodic strain signal at shallow strainmeters. The field data and the recent simulations suggest that including a periodic component to injection or pumping and then monitoring the resulting strain signal could be a way to monitor the subsurface.

Murdoch, Larry↗

DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗