Search NASA⌕ Search

SEARCH · Search NASA

Results for “BEDROCK POROSITY”

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.

25 records · Page 2

Geophysical borehole logging data of wells ER-GLS1, ER-GUM1, ER-PLM7, and ER-PLM8 at the East River Watershed, Colorado

The purpose of this data set is to provide a means of characterizing the bedrock variability both within a single borehole location and across the East River Watershed. The East River is part of the Watershed Function Scientific Focus Area (WFSFA) located in the Upper Colorado River Basin, United States, and all boreholes are located within Mancos shale bedrock. The data set contains geophysical borehole logging data acquired in October 2018 in boreholes ER-GUM1 (East River- Gothic Upper Montane Borehole #1) and ER-GLS1 (East River- Gothic Lower Subalpine Borehole #1), and acquired in October 2019 in boreholes ER-PLM7 (East River- Pumphouse Lower Montane Borehole #7) and ER-PLM8 (East River- Pumphouse Lower Montane Borehole #8). It contains depth, natural gamma ray emission, spectral gamma ray emission(Th, K, U), magnetic susceptibility, fluid temperature and electrical conductivity, electromagnetic (EM) induction derived rock electrical resistivity, heat pulse and spinner flow measurements to estimate groundwater flow, and Nuclear Magnetic Resonance (NMR) derived water content and hydraulic conductivity. Data for each borehole are provided as comma delimited text files. For each borehole accurate location information is provided (based on RTK GPS (Real Time Kinematic Global Positioning System) surveying). Borehole logging data were acquired using Mount Sopris borehole tools, and Vista Clara Javelin for NMR data.

54 ENVIRONMENTAL SCIENCES↗

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗

Instabilities and finger formation in replacement fronts driven by an oversaturated solution: FINGER FORMATION IN MINERAL REPLACEMENT

We consider a simple model of infiltration-driven mineral replacement, in which the chemical coupling between precipitation and dissolution leads to the appearance of a reaction front advancing into the system. Such fronts are usually accompanied by a local increase of porosity. We analyze the linear stability of the replacement front to establish whether such a localized porosity increase can lead to global instability and pattern formation in these systems. We find that for a wide range of control parameters such fronts are unstable. However, both short- and long-wavelength perturbations are stabilized, whereas in a purely dissolutional instability only short wavelengths are stable. We analyze the morphologies of the dissolution patterns emerging in the later stages of the evolution of the system, when the dynamics are beyond the linear regime. Implications of these results for the natural systems are discussed, particularly in the context of karst formation in terra rossa-covered carbonate bedrock.

58 GEOSCIENCES↗

Relating the depth of the water table to the depth of weathering

Weathering of bedrock creates and occludes permeability, affecting subsurface water flow. Often, weathering intensifies above the water table. On the other hand, weathering can also commence below the water table. To explore relationships between weathering and the water table, a simplified weathering model for an eroding hillslope was formulated that takes into account both saturated and unsaturated subsurface water flow (but does not fully account for changes in dissolved gas chemistry). The phreatic line was calculated using solutions to mathematical treatments for both zones. In the model, the infiltration rate at the hill surface sets both the original and the eventual steady-state position of the water table with respect to the weathering reaction front. Depending upon parameters, the weathering front can locate either above or below the water table at steady state. Erosion also affects the water table position by changing porosity and permeability even when other hydrological conditions (e.g. hydraulic conductivity of parent material, infiltration rate at the surface) do not change. The total porosity in a hill (water storage capacity) was found to increase with infiltration rate (all else held constant). Furthermore, this effect was diminished by increasing the erosion rate. We also show examples of how the infiltration rate affects the position of the water table and how infiltration rate affects weathering advance.

58 GEOSCIENCES↗

Optimization of Deep Borehole Systems for HLW Disposal

This is the final report on a project to update and improve the conceptual design of deep boreholes for high level nuclear waste disposal. The effort was concentrated on application to intact US legacy LWR fuel assemblies, but conducted in a way in which straightforward extension to other waste forms, host rock types and countries was preserved. The reference fuel design version consists of a vertical borehole drilled into granitic bedrock, with the uppermost kilometer serving as a caprock zone containing a diverse and redundant series of plugs. There follows a one to two kilometer waste canister emplacement zone having a hole diameter of approximately 40-50 cm. Individual holes are spaced 200-300 m apart to form a repository field. The choice of verticality and the use of a graphite based mud as filler between the waste canisters and the borehole wall liner was strongly influenced by the expectation that retrievability would continue to be emphasized in US and worldwide repository regulatory criteria. An advanced version was scoped out using zinc alloy cast in place to fill void space inside a disposal canister and its encapsulated fuel assembly. This excludes water and greatly improves both crush resistance and thermal conductivity. However the simpler option of using a sand fill was found adequate and is recommended for near-term use. Thermal-hydraulic modeling of the low permeability and porosity host rock and its small (≤ 1%) saline water content showed that vertical convection induced by the waste’s decay heat should not transport nuclides from the emplacement zone up to the biosphere atop the caprock. First order economic analysis indicated that borehole repositories should be cost-competitive with shallower mined repositories. It is concluded that proceeding with plans to drill a demonstration borehole to confirm expectations, and to carry out priority experiments, such as retention and replenishment of in-hole water is in order.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Coupled land surface–subsurface hydrogeophysical inverse modeling to estimate soil organic carbon content and explore associated hydrological and thermal dynamics in the Arctic tundra

Quantitative characterization of soil organic carbon (OC) content is essential due to its significant impacts on surface–subsurface hydrological–thermal processes and microbial decomposition of OC, which both in turn are important for predicting carbon–climate feedbacks. While such quantification is particularly important in the vulnerable organic-rich Arctic region, it is challenging to achieve due to the general limitations of conventional core sampling and analysis methods, and to the extremely dynamic nature of hydrological–thermal processes associated with annual freeze–thaw events. In this study, we develop and test an inversion scheme that can flexibly use single or multiple datasets – including soil liquid water content, temperature and electrical resistivity tomography (ERT) data – to estimate the vertical distribution of OC content. Our approach relies on the fact that OC content strongly influences soil hydrological–thermal parameters and, therefore, indirectly controls the spatiotemporal dynamics of soil liquid water content, temperature and their correlated electrical resistivity. We employ the Community Land Model to simulate nonisothermal surface–subsurface hydrological dynamics from the bedrock to the top of canopy, with consideration of land surface processes (e.g., solar radiation balance, evapotranspiration, snow accumulation and melting) and ice–liquid water phase transitions. For inversion, we combine a deterministic and an adaptive Markov chain Monte Carlo (MCMC) optimization algorithm to estimate a posteriori distributions of desired model parameters. For hydrological–thermal-to-geophysical variable transformation, the simulated subsurface temperature, liquid water content and ice content are explicitly linked to soil electrical resistivity via petrophysical and geophysical models. We validate the developed scheme using different numerical experiments and evaluate the influence of measurement errors and benefit of joint inversion on the estimation of OC and other parameters. We also quantify the propagation of uncertainty from the estimated parameters to prediction of hydrological–thermal responses. We find that, compared to inversion of single dataset (temperature, liquid water content or apparent resistivity), joint inversion of these datasets significantly reduces parameter uncertainty. We find that the joint inversion approach is able to estimate OC and sand content within the shallow active layer (top 0.3 m of soil) with high reliability. Due to the small variations of temperature and moisture within the shallow permafrost (here at about 0.6 m depth), the approach is unable to estimate OC with confidence. However, if the soil porosity is functionally related to the OC and mineral content, which is often observed in organic-rich Arctic soil, the uncertainty of OC estimate at this depth remarkably decreases. Our study documents the value of the new surface–subsurface, deterministic–stochastic inversion approach, as well as the benefit of including multiple types of data to estimate OC and associated hydrological–thermal dynamics.

58 GEOSCIENCES↗

Rates of subcritical cracking and long-term rock erosion

Bedrock fracture is a key element of rock erosion and subsequent surface processes. We test the hypothesis that rock's susceptibility to subcritical cracking, a specific type of fracturing, significantly drives and limits rock erosion. We measured 10 Be-derived erosion rates, compressive strength, and crack characteristics on 20 outcrops of different rock units (quartzite, granite, and two metasandstones) in the northern Blue Ridge Mountains of Virginia (USA). We also measured the subcritical cracking index (n), Charles’s law velocity constant (A), and fracture toughness (K IC ) of samples from four of the same outcrops, representative of each rock type. Erosion rates range from 1.16 ± 0.67 to 32.3 ± 7.8 m/m.y. We find strong correlations—across the four rock units—between average erosion rates and the three subcritical cracking parameters (R 2 > 0.85, p < 0.05), but not compressive strength (R 2 = 0.6; p > 0.1). We also find a correlative relationship between n and outcrop fracture length (R 2 = 0.91; p < 0.05). The latter correlation is consistent with that of published model predictions, further indicating a mechanistic link between subcritical cracking and rock erosion. We infer that subcritical cracking parameters closely tie to erosion rates, because subcritical cracking is the dominant process of mechanical weathering, leading to positive feedbacks relating subcritical cracking rates, crack length, porosity, and water accessibility. Finally, these data are the first that directly test and support the hypothesis that subcritical cracking can set the pace of long-term rock erosion.

58 GEOSCIENCES↗