Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BEDROCK”

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.

At least 37 records · Page 2

Deep Roots Supply Reactivity and Enhance Silicate Weathering in the Bedrock Vadose Zone

In upland environments, roots commonly extend deep below soil into partially saturated bedrock. This Bedrock Vadose Zone (BVZ) has been shown to store and circulate water, host organic carbon respiration and serve as a critical source of rock-derived nutrients. However, the extent to which deep roots influence chemical weathering rates remains poorly understood. Here, we report 4 years of depth-resolved major ion chemistry over a 16-m thick BVZ hosting a deep rhizosphere in a catchment subject to a Mediterranean climate. These data allow development and validation of a reactive transport model (RTM), revealing that the timescales of water storage and drainage in the BVZ are sufficient to facilitate substantial chemical weathering of the shale bedrock. However, observed solute concentrations are only reproduced by the RTM when we explicitly include measured rates of CO 2(g) production meters below soil driven by the deeply rooted forest. By combining direct observations and a process-based RTM we conclude that the carbon respiration promoted by deep roots significantly enhances chemical weathering rates in the BVZ, constituting 43% ± 3% of total solute flux from the base of the BVZ to the water table.

Osorio‐Leon, Ivan D. [University of Illinois Urban↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗

Detection of bedrock topography beneath a thin cover of alluvium using thermal remote sensing

Modeling of the annual heat flow within a thin alluvium veneer on a granitic bedrock substrate in desert environments, such as found in the southwestern United States, predicts that at certain times of the year the depth to bedrock has a measurable effect on the surface temperature if the alluvium cover is less than 2 m thick. Changes in the thickness of the alluvial cover caused by bedrock topography will produce contrasts in the surface temperature. If temperature contrasts as small as 0.1 C can be resolved, a linear topographic feature having several metres of relief buried by 1.5 m of alluvium may be visible in thermal imagery acquired during January or August in the southwestern U.S. under optimal conditions. Thermal remote sensing may provide a means for delineating some buried faults, fluvial channels, and other features of interest on buried, granitic pediment surfaces.

Nash, D. B.↗

An evaluation of probable bedrock exposure in the Sinus Meridiani region of the Martian highlands

Estimates of the amount of probable bedrock exposures visible in the highest-resolution Viking imaging data for a selected region of the Martian surface are presented. The study region includes portions of the classical low-albedo regions of Sinus Meridiani and Sinus Sabaeus and the high-albedo region of Arabia, and is completely within the cratered highlands of Mars. More than 400 images were used to estimate the percentage of the area of each frame having locally steep slopes. It is seen that 42 percent of the high-resolution frames have less than one percent probable bedrock exposure while three percent of the frames have 15 percent probable bedrock exposure.

Zimbelman, J. R.↗

Effects of Bedrock Landsliding on Cosmogenically Determined Erosion Rates

The successful quantification of long-term erosion rates underpins our understanding of landscape. formation, the topographic evolution of mountain ranges, and the mass balance within active orogens. The measurement of in situ-produced cosmogenic radionuclides (CRNs) in fluvial and alluvial sediments is perhaps the method with the greatest ability to provide such long-term erosion rates. In active orogens, however, deep-seated bedrock landsliding is an important erosional process, the effect of which on CRN-derived erosion rates is largely unquantified. We present a numerical simulation of cosmogenic nuclide production and distribution in landslide-dominated catchments to address the effect of bedrock landsliding on cosmogenic erosion rates in actively eroding landscapes. Results of the simulation indicate that the temporal stability of erosion rates determined from CRN concentrations in sediment decreases with increased ratios of landsliding to sediment detachment rates within a given catchment area, and that larger catchment areas must be sampled with increased frequency of landsliding in order to accurately evaluate long-term erosion rates. In addition, results of this simulation suggest that sediment sampling for CRNs is the appropriate method for determining long-term erosion rates in regions dominated by mass-wasting processes, while bedrock surface sampling for CRNs is generally an ineffective means of determining long-term erosion rates. Response times of CRN concentrations to changes in erosion rate indicate that climatically driven cycles of erosion may be detected relatively quickly after such changes occur, but that complete equilibration of CRN concentrations to new erosional conditions may take tens of thousands of years. Simulation results of CRN erosion rates are compared with a new, rich dataset of CRN concentrations from the Nepalese Himalaya, supporting conclusions drawn from the simulation.

Niemi, Nathan↗

Variations in Bedrock and Vegetation Cover Modulate Subsurface Water Flow Dynamics of a Mountainous Hillslope

Abstract Predicting the hydrological response of watersheds to climate disturbances requires a detailed understanding of the processes connecting hillslopes and streams. Using a network of soil moisture and temperature sensors, electrical resistivity tomography monitoring, and a weather station we assess the above and below‐ground processes driving the hydrological response of a hillslope during snowmelt and summer monsoon. The transect covers bedrock and vegetation gradients, with a steep upper part characterized by shallow bedrock, and gentle lower part underlain by colluvium. The main vegetation cover is conifers on the upper, and grass and veratrum on the lower part. Combined with a simplified hydrological model, we show that the thin soil layer of the steep slope acts as a preferential flow path, leading to mostly shallow lateral flow, interrupted by vertical flow, mostly at tree locations, and likely facilitated by flow along fractures and roots. Vertical flow and upstream‐driven groundwater dynamics are prevailing at the colluvium, presenting a very different hydrological behavior compared to the upper part. These results show that subsurface structure and features have a strong control on the hydrological response of a hillslope and that those can create considerably varying hydrological dynamics across small spatial scales.

54 ENVIRONMENTAL SCIENCES↗

Predicting sedimentary bedrock subsurface weathering fronts and weathering rates: Dataset

Although bedrock weathering strongly influences water quality and global carbon and nitrogen budgets, the weathering depths and rates within subsurface are not well understood nor predictable. Determination of both porewater chemistry and subsurface water flow are needed in order to develop more complete understanding and obtain weathering rates. In a long-term field study, we applied a multiphase approach along a mountainous watershed hillslope transect underlain by marine shale. Here we report three findings. First, the deepest extent of the water table determines the weathering front, and the range of annually water table oscillations determines the thickness of the weathering zone. Below the lowest water table, permanently water-saturated bedrock remains reducing, preventing deeper pyrite oxidation. Secondly, carbonate minerals and potentially rock organic matter share the same weathering front depth with pyrite, contrary to models where weathering fronts are stratified. Thirdly, the measurements-based weathering rates from subsurface shale are high, amounting to base cation exports of about 70 kmolc ha−1 y−1, yet consistent with weathering of marine shale. Finally, by integrating geochemical and hydrological data we present a new conceptual model that can be applied in other settings to predict weathering and water quality responses to climate change.

54 ENVIRONMENTAL SCIENCES↗

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗

Quantifying Dynamic Water Storage in Unsaturated Bedrock with Borehole Nuclear Magnetic Resonance

Quantifying the volume of water that is stored in the subsurface is critical to studies of water availability to ecosystems, slope stability, and water-rock interactions. In a variety of settings, water is stored in fractured and weathered bedrock as rock moisture. However, few techniques are available to measure rock moisture in unsaturated rock, making direct estimates of water storage dynamics difficult to obtain. Here, we use borehole nuclear magnetic resonance (NMR) at two sites in seasonally dry California to quantify dynamic rock moisture storage. Additionally, we show strong agreement between NMR estimates of dynamic storage and estimates derived from neutron logging and mass balance techniques. The depths of dynamic storage are up to 9 m and likely reflect the depth extent of root water uptake. To our knowledge, these data are the first to quantify the volume and depths of dynamic water storage in the bedrock vadose zone via NMR.

58 GEOSCIENCES↗

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

Self‐Potential Tomography Preconditioned by Particle Swarm Optimization—Application to Monitoring Hyporheic Exchange in a Bedrock River

Abstract A self‐potential (SP) data‐inversion algorithm was developed and tested on an analytical model of electrical‐potential profile data attributed to single and multiple polarized electrical sources. The developed algorithm was then validated by an application to SP‐monitoring field data measured on the floodplain of East Fork Poplar Creek, Oak Ridge, Tennessee, to image electrical sources in areas conducive to preferential flow into the flood plain from the bedrock‐lined riverbed. The algorithm combined stochastic source‐localization by particle‐swarm‐optimization (PSO) of electrical sources characterized by simplified geometries with source tomography by regularized weighted least‐squares minimization of a quadratic objective function. Prior information was incorporated by preconditioning the tomography algorithm by PSO results. Variable percentages of random noise were added to analytical‐model data to evaluate the algorithm performance. Results indicated that true parameters of single‐source models were inverted and approximated with small residual error, whereas inversion of analytical‐model data representing multiple electrical sources accurately approximated the locations of the sources but miscalculated some parameters because of the non‐uniqueness of the inverse‐model solution. Source tomography applied to analytical model data during testing produced a spatially continuous parameter field that identified the locations of point‐scale synthetic dipole sources of electrical current flow with varying degrees of accuracy depending on the prior information incorporated into the tomography. When applied to SP‐monitoring field data, the algorithm imaged electrical sources within a known fault that intersects the bedrock riverbed and flood plain of East Fork Poplar Creek and depicted dynamic electrical conditions attributed to hyporheic exchange.

54 ENVIRONMENTAL SCIENCES↗

River thorium concentrations can record bedrock fracture processes including some triggered by distant seismic events

Fractures are integral to the hydrology and geochemistry of watersheds, but our understanding of fracture dynamics is very limited because of the challenge of monitoring the subsurface. Here we provide evidence that long-term, high-frequency measurements of the river concentration of the ultra-trace element thorium (Th) can provide a signature of bedrock fracture processes spanning neighboring watersheds in Colorado. River Th concentrations show abrupt (subdaily) excursions and biexponential decay with approximately 1-day and 1-week time constants, concentration patterns that are distinct from all other solutes except beryllium and arsenic. The patterns are uncorrelated with daily precipitation records or seasonal trends in atmospheric deposition. Groundwater Th analyses are consistent with bedrock release and dilution upon mixing with river water. Most Th excursions have no seismic signatures that are detectable 50 km from the site, suggesting the Th concentrations can reveal aseismic fracture or fault events. We find, however, a weak statistical correlation between Th and seismic motion caused by distant earthquakes, possibly the first chemical signature of dynamic earthquake triggering, a phenomenon previously identified only through geophysical methods.

58 GEOSCIENCES↗

Predicting sedimentary bedrock subsurface weathering fronts and weathering rates

Although bedrock weathering strongly influences water quality and global carbon and nitrogen budgets, the weathering depths and rates within subsurface are not well understood nor predictable. Determination of both porewater chemistry and subsurface water flow are needed in order to develop more complete understanding and obtain weathering rates. In a long-term field study, we applied a multiphase approach along a mountainous watershed hillslope transect underlain by marine shale. Here we report three findings. First, the deepest extent of the water table determines the weathering front, and the range of annually water table oscillations determines the thickness of the weathering zone. Below the lowest water table, permanently water-saturated bedrock remains reducing, preventing deeper pyrite oxidation. Secondly, carbonate minerals and potentially rock organic matter share the same weathering front depth with pyrite, contrary to models where weathering fronts are stratified. Thirdly, the measurements-based weathering rates from subsurface shale are high, amounting to base cation exports of about 70 kmol c ha -1 y -1 , yet consistent with weathering of marine shale. Finally, by integrating geochemical and hydrological data we present a new conceptual model that can be applied in other settings to predict weathering and water quality responses to climate change.

58 GEOSCIENCES↗

Evolution of Titan's coupled ocean-atmosphere system and interaction of ocean with bedrock

A recent model for the surface state of Titan proposes a liquid ethane-methane-molecular nitrogen layer of order one kilometer thick which because of stratospheric methane photolysis has become increasingly ethane-rich with time. The interaction of such an ocean with the underlying 'bedrock' of Titan (assumed to be water-ice or ammonia hydrate) and with the primarily nitrogen atmosphere is explored. It is concluded that although modest exchange of oceanic hydrocarbons with enclathrated methane in the bedrock can in principle occur, it is unlikely for reasonable regolith depths. The surprisingly high solubility of water-ice in liquid methane implies that topographic features on Titan of order 100 meters in height can be eroded away on a time scale of one-billion years. The large solubility difference of N2 in methane versus ethane implies that the ocean composition is a strong determinant of atmospheric pressure; a simple radiative model of the Titan atmosphere is employed to demonstrate that significant surface pressure and temperature changes can occur as the oceanic composition evolves with time.

Lunine, J. I.↗

Lithologic controls on AIRSAR signatures of bedrock and alluvium, at Lunar Crater, Nevada

Radar backscatter intensity as measured by calibrated synthetic aperture radar (SAR) systems is primarily controlled by three factors: local incidence angle, wavelength-scale roughness, and dielectric permittivity of surface materials. In order to make adequate use of radar observations for geological investigations of surface type, the relationships between lithology and the above characteristics must be adequately understood. In arid terrains weathering signatures (e.g. fracturing, debris grain size and shape, slope characteristics) are controlled to some extent by lithologic characteristics of the parent bedrock. These textural features of outcrops and their associated debris control radar backscatter to varying degrees. The quad-polarization JPL AIRSAR system allows sampling of textures at three distinct wavelength scales: C-band (5.66 cm), L-band (23.98 cm), and P-band (68.13 cm). This paper presents a discussion of AIRSAR data using recent field observations of weathered felsic and basaltic volcanic rock units exposed in the southern part of the Lunar Crater Volcanic Field, in the Pancake Range of central Nevada. The focus is on the relationship of radar backscatter at multiple wavelengths to weathering style and parent bedrock lithology.

Rivard, Benoit↗

Bedrock Denudation on Titan: Estimates of Vertical Extent and Lateral Debris Dispersion

Methane rainfall and runoff, along with aeolian activity, have dominated the sculpting of Titan s landscape. A knowledge of the vertical extent of bedrock erosion and the lateral extent of the resulting sediment is useful for several purposes [1]. For instance, what is the magnitude and expression of modification of constructional landforms (e.g., mountains)? Does highland denudation and the filling of basins with sediment cause adjustments (uplift and subsidence) in the crustal ice shell? Here we report preliminary findings of putative eroded craters and the results of landform evolution modeling (Fig. 1) that suggest that approx. 250 m of net bedrock erosion has at least locally taken place and approx.1 km of maximum local erosion.

Moore, Jeffrey↗

Castagna, M., M. W. Becker, and A. Bellin (2011), Joint estimation of transmissivity and storativity in a bedrock fracture, Water Resour. Res., 47, W09504, doi:10.1029/2010WR009262.

The spatial variability of transmissivity (T) and storativity (S) in a sandstone bedding plane fracture was investigated through field and numerical hydraulic tomography experiments. Cross-hole pumping tests were conducted in a single subhorizontal bedrock fracture, thereby constraining the problem to two dimensions. A Bayesian method implementing pilot points was used to invert the drawdown, under the hypothesis that T and S are independent random space functions, where integral scale, mean, and variance of the random fields were parameters of inversion in addition to pilot point values. An ensemble of 20 Monte Carlo realizations was created for each model scenario. Inversion of the measured cross-hole drawdown resulted in spatially variable T and S fields. To determine the accuracy of the method, a suite of synthetic inversions modeled after the field site were also conducted. The synthetic experiments demonstrate that the assumed spatial structure of S influences the accuracy of both the T and S inversions. When the S field is spatially heterogeneous and is assumed to be heterogeneous in the inversion, a consistent representation of the T and S fields is obtained. When the S field is spatially heterogeneous but is assumed to be homogeneous in the inversion, a consistent representation of the T field is still obtained, but the S field is entirely misrepresented; if the S field is actually homogeneous but assumed to be heterogeneous, a poor representation of the T and S fields is obtained along with high uncertainty in the inversion statistics of confidence.

54 ENVIRONMENTAL SCIENCES↗