Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “proxies”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development of Heterogeneous Ice Nucleation Rate Coefficient Parameterizations From Ambient Measurements

Ice nucleating particles (INPs) are a rare subset of particles that can have an outsized impact relative to their prevalence. To simulate INP variability, models require parameterizations for the most important sources of INPs. Most parameterizations in the literature were developed from laboratory experiments that used proxies for ambient particles. Whether these laboratory models effectively represent ambient aerosol is still not well understood and parameterizations developed from field measurements may be more atmospherically relevant. Expanding on previous work in which we used an ice chamber and a single-particle mass spectrometer (SPMS) to characterize INP composition, we calculate heterogeneous ice nucleation rate coefficients for ambient particles. We find good agreement between our dust and previous measurements of airborne dust. We also use a Monte Carlo approach to assess the relative contribution of each measurement to the uncertainty and find that the biggest source is the sampling efficiency of the SPMS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Increasing Functional Diversity in a Global Land Surface Model Illustrates Uncertainties Related to Parameter Simplification

Abstract Simulations of the land surface carbon cycle typically compress functional diversity into a small set of plant functional types (PFT), with parameters defined by the average value of measurements of functional traits. In most earth system models, all wild plant life is represented by between five and 14 PFTs and a typical grid cell (≈100 × 100 km) may contain a single PFT. Model logic applied to this coarse representation of ecological functional diversity provides a reasonable proxy for the carbon cycle, but does not capture the non‐linear influence of functional traits on productivity. Here we show through simulations using the Energy Exascale Land Surface Model in 15 diverse terrestrial landscapes, that better accounting for functional diversity markedly alters predicted total carbon uptake. The shift in carbon uptake is as great as 30% and 10% in boreal and tropical regions, respectively, when compared to a single PFT parameterized with the trait means. The traits that best predict gross primary production vary based on vegetation phenology, which broadly determines where traits fall within the global distribution. Carbon uptake is more closely associated with specific leaf area for evergreen PFTs and the leaf carbon to nitrogen ratio in deciduous PFTs.

58 GEOSCIENCES↗
DOE OSTI2021

Model Evaluation and Intercomparison of Marine Warm Low Cloud Fractions With Neural Network Ensembles

Abstract Low cloud fractions (LCFs) and meteorological factors (MFs) over an oceanic region containing multiple cloud regimes are examined for three data sets: one Energy Exascale Earth System Model (E3SM) simulation with the default 72‐layer vertical grid (E3SM72), another one with 8‐times vertical resolution via the Framework for Improvement by Vertical Enhancement (E3SM 8), and one with MFs from ERA5 reanalysis and LCFs from the CERES SSF product (ERA5‐SSF). Neural networks (NNs) are trained to capture the relationship between MFs and LCF and to select the best‐performing MF subsets for predicting LCF. NN ensembles are used to (a) confirm the performance of selected MF subsets, (b) to serve as proxy models for each data set to predict LCFs for MFs from all data sets, and (c) to classify MFs into those in shared and uniquely occupied MF subspaces. Overall, E3SM72 and E3SM 8 have large fractions of MFs in shared MF subspace, but less so near the Californian and Peruvian stratocumulus decks. E3SM 8 and ERA5 have small fractions of MFs in shared MF subspace but greater than E3SM72 and ERA5, especially in the Southeast Pacific. The differences in LCFs between three pairs of data sets are decomposed into those associated with the differences in the LCF‐MF relationship and those involving different MFs. Given the same MFs, LCFs produced by E3SM 8 are greater than those produced by E3SM72 but are still different from those in ERA5‐SSF. In general, the shift in MFs dominates the difference in the LCFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Achievements and Prospects of Global Broadband Seismographic Networks After 30 Years of Continuous Geophysical Observations

Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.

58 GEOSCIENCES↗
DOE OSTI2023

Earlier Ecological Drought Detection by Involving the Interaction of Phenology and Eco-Physiological Function

Meteorological drought indices like the Standardized Precipitation Evaporation Index (SPEI) are frequently used to diagnose “ecological drought,” despite the fact that they were not explicitly designed for this purpose. More recently developed indices like the Evaporative Stress Index (ESI), which is based on the degree of coupling between actual to potential evapotranspiration, may better capture dynamic plant response to moisture limitations. However, the skill of these indices at describing plant water stress is rarely evaluated at sub-seasonal timescales over which drought evolves. Moreover, it remains unclear how variability in phenological timing impacts and complicates early drought detection. Here, we compared the ability of ESI and SPEI to reflect the dynamics of ecological drought in forests and grasslands, based on anomalies of Gross Primary Productivity (GPP), surface conductance (G s , a proxy for stomatal conductance), soil moisture, and vapor pressure deficit. ESI performed better than SPEI in capturing the dynamics of GPP and G s , but still missed early ecological drought signals due to biases linked to earlier onset of spring leaf development. Thus, we developed a modified variant of the ESI (ESI LAI ) that accounts for the complicating effects of phenological shifts in leaf area index (LAI). The ESI LAI detected drought onset up to 7–10 weeks earlier than SPEI and ESI. Additionally, drought onset dates determined from ESI LAI are close to (±2 weeks) the dates determined from LAI-corrected anomalies of G s , and GPP, as well as the onset dates of soil water deficit and atmospheric aridity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Uncertainty in Reconstructing Paleo‐Elevation of the Antarctic Ice Sheet From Temperature‐Sensitive Ice Core Records

Paleotemperature reconstructions from ice cores are mixed signals of changes in climate and ice‐surface elevation. A common, temperature‐based paleoaltimetry method suggests these signals can be disentangled by comparing two proxy locations with similar climates. The difference between the records is assumed to be due to elevation, which is estimated by scaling the temperature difference by a lapse rate. We investigate the uncertainty associated with this approach using a case study of the Antarctic Ice Sheet during the Last Glacial Maximum. From an ensemble of climate simulations, we extract modeled temperatures at locations of real ice cores. We find uncertainty on the order of hundreds of meters that results from spatial heterogeneity in non‐adiabatic temperature change, which itself stems in part from elevation‐induced atmospheric circulation change. Our findings suggest that caution is needed when interpreting temperature‐based paleoaltimetry results for ice sheets.

58 GEOSCIENCES↗
DOE OSTI2022

Using Multiscale Ethane/Methane Observations to Attribute Coal Mine Vent Emissions in the San Juan Basin From 2013 to 2021

Abstract Source attribution of natural gas emissions from fossil fuels in New Mexico's San Juan Basin (SJB) is challenging due to source heterogeneity and emissions transience. We demonstrate that ethane (C 2 H 6 ) to methane (CH 4 ) mixing ratios can identify and separate sources over different scales using various measurement techniques. We report simultaneous CH 4 and C 2 H 6 observations near a coal mine vent and oil and gas (O&G) emission sources using ground‐based in situ measurements in 2020/2021. During these campaigns, we observed a stable coal vent C 2 H 6 :CH 4 ratio of 1.28% ± 0.11%, discernibly different than nearby O&G source ratios ranging from 0.9% to 16.8%. We analyze airborne observations of the SJB taken in 2014/2015 that exhibit similar coal vent ratios and further show the region's heterogeneity. We identify episodic O&G sources, including a gas plant source detected in 2014/2015 that is absent in our 2020/2021 data. We examine total column observations of C 2 H 6 and CH 4 made in 2013 with a solar spectrometer and find a C 2 H 6 :CH 4 ratio of 1.3% ± 0.4% for the coal vent. The stable and unique coal vent ratio relative to other O&G sources in the region is used to demonstrate that consistent attribution is possible using various measurement methods at multiple scales across many years. Finally, we demonstrate that using C 2 H 6 as a proxy for fossil CH 4 inversions can inform detailed basin‐scale inversions, provided we understand source specific changes in the C 2 H 6 :CH 4 ratio like we report in the SJB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Microtopography Matters: Belowground CH 4 Cycling Regulated by Differing Microbial Processes in Peatland Hummocks and Lawns

Water table depth and vegetation are key controls of methane (CH 4 ) emissions from peatlands. Microtopography integrates these factors into features called microforms. Microforms often differ in CH 4 emissions, but microform-dependent patterns of belowground CH 4 cycling remain less clearly resolved. Here, to investigate the impact of microtopography on belowground CH 4 cycling, we characterized depth profiles of the community composition and activity of CH 4 -cycling microbes using 16S rRNA amplicon sequencing, incubations, and measurements of porewater CH 4 concentration and isotopic composition from hummocks and lawns at Sallie's Fen in NH, USA. Geochemical proxies of methanogenesis and methanotrophy indicated that microforms differ in dominant microbial CH 4 cycling processes. Hummocks, where water table depth is lower, had higher porewater redox potential (Eh) and higher porewater δ 13 C-CH 4 values in the upper 30 cm than lawns, where water table depth is closer to the peat surface. Porewater δ 13 C-CH 4 and δD-CH 3 D values were highest at the surface of hummocks where the ratio of methanotrophs to methanogens was also greatest. These results suggest that belowground CH 4 cycling in hummocks is more strongly regulated by methanotrophy, while in lawns methanogenesis is more dominant. We also investigated controls of porewater CH 4 chemistry. The ratio of the relative abundance of methanotrophs to methanogens was the strongest predictor of porewater CH 4 concentration and δ 13 C-CH 4 , while vegetation composition had minimal influence. As microbial community composition was strongly influenced by redox conditions but not vegetation, we conclude that water table depth is a stronger control of belowground CH 4 cycling across microforms than vegetation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Characterizing Performance of Freshwater Wetland Methane Models Across Time Scales at FLUXNET-CH 4 Sites Using Wavelet Analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

58 GEOSCIENCES↗
DOE OSTI2023

Accounting for Changes in Radiation Improves the Ability of SIF to Track Water Stress-Induced Losses in Summer GPP in a Temperate Deciduous Forest

Global observations of solar-induced chlorophyll fluorescence (SIF) are available from multiple satellite platforms, and SIF is increasingly used as a proxy for photosynthetic activity and ecosystem productivity. Because the relationship between SIF and gross primary productivity (GPP) depends on a variety of factors including ecosystem type and environmental conditions, it is necessary to study SIF observations across various spatiotemporal scales and ecosystems. To explore how SIF signals relate to productivity over a temperate deciduous forest, we deployed a PhotoSpec spectrometer system at the University of Michigan Biological Station AmeriFlux site (US-UMB) in the northern Lower Peninsula of Michigan during the 2018 and 2019 growing seasons. We found that SIF correlated with GPP across diurnal and seasonal cycles (R 2 = 0.61 and 0.64 for 90-min- and daily-averaged data), but that SIF signals were more strongly related to downwelling radiation than GPP (R 2 = 0.91 for daily-averaged data). The dependence of SIF on radiation obscured the impact of intraseasonal drought in the SIF timeseries, but drought stress was apparent as a decrease in relative SIF, which exhibited a stronger correlation with GPP (R 2 = 0.56) than other remotely sensed data over the drought period. These results highlight the potential of SIF for detecting stress-induced losses in forest productivity. Additionally, we found that the red:far-red SIF ratio did not exhibit a response to water stress-induced losses in productivity, but was largely driven by seasonal and interannual changes in canopy structure, as well as by synoptic changes in downwelling radiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving E3SM Land Model Photosynthesis Parameterization via Satellite SIF, Machine Learning, and Surrogate Modeling

The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Zonal Patterns in Late Quaternary Tropical South American Precipitation

Speleothem oxygen isotope records (δ 18 O) of tropical South American rainfall in the late Quaternary show a zonal “South American Precipitation Dipole” (SAPD). The dipole is characterized by opposing east-west precipitation anomalies compared to the present—wetter in the east and drier in the west at the mid-Holocene (~7 ka), and drier in the east and wetter in the west at the Last Glacial Maximum (~21 ka). However, the SAPD remains enigmatic because it is expressed differently in western versus eastern δ18O records and isotope-enabled climate model simulations usually misrepresent the magnitude and/or spatial pattern of δ 18 O change. Here, we address the SAPD enigma in two parts. First, we re-interpret the δ 18 O data to account for upwind rainout effects that are known to be pervasive in tropical South America, but are not always considered in Quaternary paleoclimate studies. Our revised interpretation reconciles the δ 18 O data with cave infiltration and other proxy records, and indicates that the centroid of tropical South American rainfall has migrated zonally over time. Second, using an energy balance model of tropical atmospheric circulation, we hypothesize that zonal migration of the precipitation centroid can be explained by regional energy budget shifts, such as changing Saharan albedo associated with the African Humid Period, that have not been modeled in previous SAPD studies. This hypothesis of a migrating precipitation centroid presents a new framework for interpreting δ18O records from tropical South America and may help explain the zonal rainfall anomalies that predate the late Quaternary.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Detecting Permafrost Active Layer Thickness Change From Nonlinear Baseflow Recession

We report permafrost underlies about one fifth of the global land area and affects ground stability, freshwater runoff, soil chemistry, and surface-atmosphere gas exchange. The depth of thawed ground overlying permafrost (active layer thickness) has broadly increased across the Arctic in recent decades, coincident with a period of increased streamflow, especially the lowest flows (baseflow). Mechanistic links between active layer thickness and baseflow have recently been explored using linear reservoir theory, but most watersheds behave as nonlinear reservoirs. We derive theoretical nonlinear relationships between long-term average saturated soil thickness $\bar{η}$(proxy for active layer thickness) and long-term average baseflow. When applied to 38 years of daily streamflow data for the Kuparuk River basin on the North Slope of Alaska, the theory predicts $\bar{η}$ increased 0.17 ± 0.22[2σ] cm a -1 between 1983 and 2020 (6.4 ± cm total). The rate of increase nearly doubled to 0.29 ± 0.31 cm a -1 between 1990 and 2020, during which time local field measurements from Circumpolar Active Layer Monitoring sites indicate the active layer increased 0.31 ± 0.22 cm a -1 . The predicted rate of increase more than doubled again between 2002 and 2020, outpacing a near doubling of observed active layer thickening, consistent with trends in terrestrial water storage inferred from Gravity Recovery and Climate Experiment satellite gravimetry and Modern-Era Retrospective Analysis for Research and Applications climate reanalysis. Overall, hydrologic change is accelerating in the Kuparuk River basin, and we provide a theoretical framework for estimating basin-scale changes in active layer water storage from streamflow measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Thunderstorms With the Greatest Lightning Densities on Earth

The most intense thunderstorms on Earth were surveyed using the comprehensive meteorological instrumentation on the Tropical Rainfall Measuring Mission (TRMM) satellite. Expansive land-based Mesoscale Convective Systems (MCSs) were consistently identified among the Earth's most intense thunderstorms, with their organization into many convective cells spanning a large areal extent permitting exceptional overall flash rates for these storms. In this study, we identify a new class of extreme thunderstorm. Lightning-dense thunderstorms are relatively compact convective storms whose concentrated lightning activity hinders our ability to accurately measure their flash rates. The top storms have a flash rate of one flash spanning many seconds, as there is insufficient separation to distinguish one flash from another. While any particularly active convective cell could be capable of producing high lightning densities, we find that thunderstorms with the greatest lightning densities on Earth are found in maritime thunderstorms that have not been appreciated in prior work due to the inaccurate flash rate measurements. These storms that are mostly found throughout the Gulf of Mexico and east of South Africa (among other coastal and oceanic regions) have measured TRMM proxies for convective intensity that rival the top MCS thunderstorms, but their horizontal and vertical dimensions are small by comparison. Thus, the necessary microphysical elements for electrification processes are more highly concentrated, enabling the observed extreme lightning densities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying Seasonal and Diurnal Cycles of Solar‐Induced Fluorescence With a Novel Hyperspectral Imager

Abstract Solar‐induced fluorescence (SIF) is a proxy of ecosystem photosynthesis that often scales linearly with gross primary productivity (GPP) at the canopy scale. However, the mechanistic relationship between GPP and SIF is still uncertain, especially at smaller temporal and spatial scales. We deployed a ultra‐hyperspectral imager over two grassland sites in California throughout a soil moisture dry down. The imager has high spatial resolution that limits mixed pixels, enabling differentiation between plants and leaves within one scene. We find that imager SIF correlates well with diurnal changes in leaf‐level physiology and gross primary productivity under well‐watered conditions. These relationships deteriorate throughout the dry down event. Our results demonstrate an advancement in SIF imaging with new possibilities in remotely sensing plant canopies from the leaf to the ecosystem. These data can be used to resolve outstanding questions regarding SIF's meaning and usefulness in terrestrial ecosystem monitoring.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2024

Relating Hydro–Mechanical and Elastodynamic Properties of Dynamically Stressed Tensile–Fractured Rock in Relation to Applied Normal Stress, Fracture Aperture, and Contact Area

We exploit nonlinear elastodynamic properties of fractured rock to probe the micro-scale mechanics of fractures and understand the relation between fluid transport and fracture aperture under dynamic stressing. Experiments were conducted on rough, tensile-fractured Westerly granite subject to triaxial stresses. We measure fracture permeability for steady-state fluid flow with deionized water. Pore pressure oscillations are applied at amplitudes ranging from 0.2 to 1 MPa at 1 Hz frequency. During dynamic stressing we transmit ultrasonic signals through the fracture using an array of piezoelectric transducers (PZTs) to monitor evolution of interface properties. We examine the influence of fracture aperture and contact area by conducting measurements at effective normal stresses of 10–20 MPa. Additionally, the evolution of contact area with stress is characterized using pressure sensitive film. These experiments are conducted separately with the same fracture and map contact area at stresses from 9 to 21 MPa. The measurements are a proxy for “true” contact area for the fracture surface and we relate them to elastic properties using the calculated PZT sensor footprints via numerical modeling of Fresnel zones. We compare the elastodynamic response of the fracture using the stress-induced changes in ultrasonic wave velocities for transmitter-receiver pairs to image spatial variations in contact properties. We show that nonlinear elasticity and permeability enhancement decrease with increasing normal stress. Additionally, post-oscillation wave velocity and permeability exhibit quick recoveries toward pre-oscillation values. Estimates of fracture contact area (global and local) demonstrate that the elastodynamic and permeability responses are dominated by fracture topology.

58 GEOSCIENCES↗