SEARCH · Engineering Papers
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.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
An integrated approach was developed to assess the sustainability, origin, age, evolution, and groundwater potential of large fossil aquifers using the Mega Aquifer System (MAS) (area: 2 x 10(6) km(2)) in the Arabian Peninsula as a test site. A two-step exercise was adopted, the sustainability of the MAS was first investigated through the analysis of GRACE data and land surface model (LSM) outputs. Then, integrated geochemical, remote sensing, and field studies were conducted to address the age, origin, hydrogeological setting, and evolution of the southern sections of the MAS (Rub Al Khali Aquifer System (RAKAS)), the area identified from the analysis of GRACE and ISMs as being the most suited for sustainable development. Analysis of GRACE and ISMs revealed: (1) the MAS central and northern sections are experiencing high groundwater extraction (6.6 km(3)/yr) and depletion rates (-2.8 +/- 0.8 km(3)/yr) with minimum balancing potential through aquifer capture processes, yet sustainability could be achieved by reducing annual extraction by 2.8 km(3), and (2) the MAS southern sections (Rub Al Khali Aquifer System (RAKAS)) are experiencing low groundwater depletion rates (eastern RAKAS: -1.8 +/- 1.4 mm/yr) to steady-state conditions (western RAKAS: -0.73 +/- 1.4 mm/yr). Geochemical, remote sensing, and field investigations over the RAKAS revealed: (1) the presence of west to east trending drainage networks and large-scale groundwater flow systems consistent with a meteoric source (precipitation over Red Sea Hills); (2) increasing Cl-36 model ages along groundwater flow directions (up to 970,000 years), indicating aquifer recharge in wet Pleistocene periods; (3) progressive depletion in the O and H stable isotopic compositions of aquifers with increasing distance from the Red Sea Hills basement outcrops, indicating modest recharge during prevailing dry conditions; and (4) the presence of relatively fresh (TDS: 800-2800 mg/L) and nonradioactive (Ra-226 + Ra-228 < 0.185 Bq/kg) water in western RAKAS. Findings suggest that sustainable agricultural development is achievable at current extraction rates in western RAKAS and provide a replicable and cost-effective model.
Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.
This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.
Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.
The sedimentary fill of the Kalahari Basin, which extends across several countries in southern and central Africa, records landscape evolution processes and holds archeological evidence of early hominid occupation. Recent studies have demonstrated that the majority of the Kalahari Group sediments were deposited between the Pliocene and the recent present. However, due to limited access to natural outcrops in the flat Kalahari topography, the chronology of the sequence, mostly its earlier part, is not well constrained.
To date, only three Homo habilis specimens have been discovered that have associated craniodental and postcranial elements, providing a limited fossil record of the ontogeny and morphology of early members of the genus Homo. Recently, a nearly complete dentition, likely attributable to H. habilis, was discovered and excavated from early Pleistocene-age fluvial-lacustrine sediments of the upper Burgi Member of the Koobi Fora Formation at site F25787 in Area 13, near Ileret, Kenya. On the surface less than 15 m away, at site F25966, postcranial elements were found, which, if from the same individual as the nearby dentition, would represent the fourth associated craniodental and postcranial assemblage of this species. We developed a geochemical taphonomic history of these ca. 2 Ma hominin fossils using nondestructive X-ray based microanalytical tools (synchrotron and benchtop X-ray fluorescence chemical imaging and micro- and nano-computed tomography volumetric reconstruction), bulk analyses of sediments and paleosols at the excavation sites, and sedimentologic and stratigraphic observations. We integrate the chemical and physical taphonomic histories to test whether teeth (excavated in situ) and postcranial bones (eroded onto the outcrop surface) derive from a single individual. Minor differences in taphonomic history are attributable to the different biomineral properties of the dental and osseous components and to differences in physical damage during early post-mortem scavenging, dispersal, and burial in adjacent depositional settings. Microscale geochemical mapping enabled the temporal ordination of chemical and physical events in the specimens’ chemical taphonomic histories. Specifically, authigenic Fe- and K-bearing clays and Y, U, and Sr uptake occurred in post-burial fractures in bones and were also incorporated pervasively throughout dentin in teeth. Barite mineralization occurred along the latest fractures in both materials, and as a coating on tooth roots. The stratigraphic, taphonomic, and geochemical evidence supports the interpretation that the hominin fossils represent a single individual. Finally, successful application of these nondestructive sample characterization methods demonstrates capabilities for thorough interrogation of the taphonomic histories of other potential hominin fossil associations, enabling more robust and accurate palaeontologic constraints using those relationships.
This study reports a previously unidentified paleoenvironmental record found in sand dunes of the Atacama Desert, Chile. Long-term aeolian deflation by prevailing onshore winds has resulted in the deposition of sand on the irregular surface of a Miocene-aged anhydrite outcrop. Two deposits ~25 km apart, along the prevailing wind trajectory, were hand excavated then analyzed for vertical (and temporal) changes in physical and chemical composition. Radiocarbon ages of organic matter embedded within the deposits show that rapid accumulation of sediment began at the last glacial maximum and slowed considerably after the Pacific Ocean attained its present post-glacial level. Over this time period, grain sizes are seen to increase while accumulation rates simultaneously decrease, suggesting greater wind speeds and/or a change or decrease in sediment supply. Changes in δ 34 S values of sulfate in the sediment beginning ~10 ka indicate an increase in marine sources. Similarly, δ 2 H values from palmitic acid show a steady increase at ~10 ka, likely resulting from aridification of the region during the Holocene. Finally, due to the extreme aridity in the region, these sand dunes retain a well-preserved chemical record that reflects changes in elevation and coastal proximity after the last glacial maximum.
In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.
Pore structure is an important parameter to quantify the reservoir rock adsorption capability and diffusivity, both of which are fundamental reservoir properties to evaluate the gas production and carbon sequestration potential for coalbed methane (CBM) and shale gas reservoirs. In this study, we applied small-angle neutron scattering (SANS) to characterize the total and accessible pore structures for two coal and two shale samples. We carried out in situ SANS measurements to probe the accessible pore structure differences under argon, deuterated methane (CD 4 ), and CO 2 penetrations. The results show that the total porosity ranges between 0.25 and 5.8% for the four samples. Less than 50% of the total pores are accessible to CD 4 for the two coals, while more than 75% of the pores were found to be accessible for the two shales. This result suggests that organic matter pores tend to be disconnected compared to mineral matter pores. Argon pressurization can induce pore contraction because of the mechanical compression of the solid skeleton in both the coal and shale samples. Hydrostatic compression has a higher effect on the nanopores of coal and shale with a higher accessible porosity. Both methane and CO 2 injection can reduce the accessible nanopore volume due to a combination of mechanical compression, sorption-induced matrix swelling, and adsorbed molecule occupation. CO 2 has higher effects on sorption-induced matrix swelling and pore filling compared to methane for both the coal and shale samples. Gas densification and pore filling could occur at higher pressures and smaller pore sizes. In conclusion, the compression and adsorption could create nanopores in the San Juan coal and Marcellus shale drilled core but could have an opposite effect in the other samples, namely, the processes could damage the nanopores in the Hazleton coal and Marcellus shale outcrop.
Abstract A strong decrease of volume and density of North Pacific Subtropical Mode Water (NPSTMW) in 1999 was analyzed in a regional high‐resolution (0.1°) numerical ocean circulation model simulation. Both shoaling of the bottom and deepening of the top of the NPSTMW layer contributed equally to volume decrease. They were locally governed by different physical processes, but both seem to be associated with basin‐wide changes in wind. A westward propagating negative thermocline depth anomaly that developed in the Central Pacific when the Pacific Decadal Oscillation index changed from a positive to a negative phase in 1998 caused shoaling of the bottom of the NPSTMW layer. Deepening of the top of the NPSTMW layer was due to an increase in the near surface stratification, caused by an increase in wind‐driven lateral heat transport convergence by the Kuroshio Extension jet starting in 1997. Both processes increased the potential vorticity in the NPSTMW region, decreasing the volume of water in the NPSTMW density range that satisfied the low potential vorticity constraint that is part of the definition of “mode water.” The strong near‐surface density decrease provided preconditioning for preferential surface formation of a lighter variety of NPSTMW, further decreasing its density. It also resulted in decrease of the outcrop window in the NPSTMW density range, strongly reducing its formation rate in 1998 and 1999 despite strong surface heat loss.
In this work we investigate the effects of geological stress on fluid flow and tracer transport in natural fracture networks. We show the emergence of non-Fickian (anomalous) transport from the interplay among fracture network geometry, aperture heterogeneity, and geological stress. In this study, we extract the fracture network geometry from the geological map of an actual rock outcrop, and we simulate the geomechanical behavior of fractured rock using a hybrid finite-discrete element method. We analyze the impact of stress on the aperture distribution, fluid flow field, and tracer transport properties. Both stress magnitude and orientation have strong effects on the fracture aperture field, which in turn affects fluid flow and tracer transport through the system. We observe that stress anisotropy may cause significant shear dilation along long, curved fractures that are preferentially oriented to the stress loading. This, in turn, induces preferential flow paths and anomalous early arrival of tracers. An increase in stress magnitude enhances aperture heterogeneity by introducing more small apertures, which exacerbates late-time tailing. This effect is stronger when there is higher heterogeneity in the initial aperture field. To honor the flow field with strong preferential flow paths, we extend the Bernoulli Continuous Time Random Walk model to incorporate dual velocity correlation length scales. The proposed upscaled transport model captures anomalous transport through stressed fracture networks and agrees quantitatively with the high-fidelity numerical simulations.
Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.
Dehydration reactions within subducted oceanic crust are important for fluid-mediated element transfer within the subducting plate and potentially to the mantle wedge. The effects of metamorphic reactions and fluid flow on element recycling that occur during retrogression and exhumation of subducted continental crust from mantle depths are poorly understood. We study two metabasite pods with fresh eclogite cores and retrogressed amphibolite-facies rims and surrounding host gneiss within the Western Gneiss Region (WGR), Norway, to better understand element mobility and mass transfer during exhumation of subducted continental crust. Bulk-rock data were collected from samples taken across the pod and into the host gneiss. Phengite breakdown in eclogite and epidote recrystallization in veins and/or gneiss within pod cores contributed large ion lithophile elements and REE to retrogressed eclogite closest to the pod cores. In gneiss hosting the pods, phengite and epidote breakdown provided fluid that mediated elemental transfer and redistribution to the pod rim or tail. Compared to the studied pod in the southern WGR, the pod in the northern WGR underwent higher P-T conditions, partial melting and higher strain rates. This resulted in the infiltration of external fluid farther into the pod interior from the rim and facilitated larger mass gain in trace elements in the amphibolite tail of the pod relative to fresh eclogite in the core. The results show clear evidence for retrogression dehydration reactions driving significant fluid-mediated element redistribution as observed on the outcrop scale during exhumation following ultrahigh-pressure metamorphism, which directly impacts element signatures within the exhuming crust.
Abstract A multimethod geoelectric survey was implemented between January and March 2022 along a 220‐m long reach of the bedrock‐lined streambed of East Fork Poplar Creek in Oak Ridge, Tennessee to identify locations of surface‐water and groundwater exchange and characterize the subsurface flow paths that convey water between the stream and flood plain. A waterborne self‐potential (WaSP) survey was completed in January 2022 to measure the electric streaming‐potential field in the stream. Electric resistivity tomography (ERT) was performed in March 2022 on the flood plain adjacent to the WaSP survey reach to map the electric resistivity distribution and characterize the hydrogeology and subsurface flow paths that facilitate surface‐water and groundwater exchange in the bedrock‐lined stream. The combination of WaSP and ERT data support the qualitative interpretation that surface‐water and groundwater exchange likely occurs along fractures in outcropping bedrock and along two fault lines that intersect the limestone creek bed.
Buoyancy forcing in the subpolar North Atlantic Ocean (SPNA) is an important driver of the Atlantic Meridional Overturning Circulation (AMOC). To advance understanding of the mechanisms connecting the two processes and their relative importance in sub-basins within the SPNA, we apply the Water Mass Transformation Framework to a matched-pair of forced, ocean-sea ice simulations configured at eddying and non-eddying resolution. Within the first decade of simulation, the non-eddying simulation produces a weak AMOC, ~11 Sv at 26.5 ° N, while the AMOC in the eddying simulation is more realistic and is thus analyzed for comparison. Surface water mass transformation and boundary transport are calculated in both density and temperature–salinity coordinates in three separate deep water formation regions, the Iceland Basin, the Irminger Sea and Labrador Sea, during the first decade of both simulations. We identify strong surface freshening in the Irminger and Labrador seas in the non-eddying simulation during the AMOC decline. This freshening significantly reduces the density of the surface outcrops where surface water mass formation occurs, essentially removing this contribution to deep water formation. Concurrently, boundary transports in these two regions are warmer and saltier in the non-eddying simulation compared to the eddying simulation where the density of surface water mass formation is stable. The warming and salinification of boundary transports in the non-eddying simulation is interpreted as an obstacle to deep water formation. Labrador Sea surface water mass transformation is important, despite its small contribution to buoyancy overturning.
Significance To understand groundwater in weathered bedrock aquifers, we explore (bio)geochemical reactions using geophysical measurements. The result leads to a conceptual model for the three-dimensional distribution of geochemical alteration in the subsurface—well beyond the typical approach of inferring (bio)geochemical reactions in the subsurface from pits, boreholes, and outcrops. With this geophysical technique we discovered the prevalence of gas beneath the water table at this site. This geophysical viewing of the subsurface highlights patterns of water movement, reactions of clays, and microbial and abiotic gas-producing reactions recorded in the subsurface and thus illuminates storage, flow, and reactions of groundwater in bedrock systems. Protecting the groundwater resource beneath our feet will rely on such enhanced capacities to see into the subsurface.
Here, three-dimensional (3-D) nonlinear site response analyses are conducted using finite element models of actual soil profiles from ten nuclear power plant (NPP) sites in the United States to investigate the effects of soil properties and input motions on site amplification. The modeling approach developed in this study combines several novel elements, such as 3-D analysis (including vertical motions), nonlinear inelastic behavior of soil (strain-dependent shear modulus reduction and hysteretic damping), formulation of nonreflecting boundary conditions at the base, and generation of realistic outcrop ground motions for specific sites. All these elements of the modeling approach are first validated using actual data from five earthquakes at three downhole array stations recorded in the Kiban-Kyoshin network (KiK-net), Japan. The same approach is then used to develop site models of ten NPP sites in the United States and corresponding ground motions that are spectrally matched to the site hazard spectra. Eight sets of three-component input motions are used in the study and are categorized on the basis of presence or absence of a near-field pulse in the seed ground motions used for spectral matching. It is found that all sites retain a definite site amplification function regardless of the input motion, provided that the seed motion is spectrally matched to the site hazard spectra. The magnitude of site amplification and frequencies at which they occur depend upon soil properties, particularly the shear wave velocity profile and the constitutive relationship (strain-dependent shear modulus reduction and hysteretic damping) of soil. Amplification of spectral acceleration in the vertical direction (up-down motion) is found to be just as much as, if not more than, the amplification in the horizontal direction. Peak shear strain is found to be about 20% larger for near-field motions compared to far-field motions whereas maximum horizontal site amplification for far-field motions is found to be consistently larger than that of near-field motions, even though the differences between the two remain within the scatter resulting from individual ground motions.