Preliminary Analysis of the Cenozoic Section from Core Drilled in the Rock Valley Fault Zone, Nevada National Security Site
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
Sierra Nevada forests transpire a significant amount of California's water resources, sparking interest in applying forest management to improve California's water supply. Determining the source water of evapotranspiration enables forest managers to make informed decisions. To this end, a significant interest in critical zone science is to develop new methods to work across time scales to predict subsurface water storage and use. Here, forest vegetation accessed young water and switched sources depending on availability, suggesting that forest drought vulnerability may depend on the range of water sources available (rain, snowmelt and deeply stored water). This finding also suggests that changes in transpiration rates may have immediate effects on water sources in close proximity to vegetation, and delayed effects on storage and runoff. New δ 18 O, δ 2 H and 3 H data were used to track precipitation, runoff, evapotranspiration and storage through the critical zone seasonally, including seasons where evapotranspiration and snowmelt were in phase (winter snowmelt) and out of phase (seasonally dry summer). The main source of this headwater catchment's runoff is derived from its meadow saturated zone water, which was dominated by snowmelt. Water that originated as snowmelt contributed to transpiration, unless other sources, such as recent rain, became available. In cases where xylem δ 18 O and δ 2 H signatures matched those of deeper saturated zone water, 3 H data showed that xylem water was distinctly younger than the deep saturated zone water. During 2016, which experienced relatively normal snowpack in winter and seasonally dry summer conditions, mean summer saturated zone water and vegetation water were similar in δ 18 O, –12.4 ± 0.04 ‰ and – 12.5 ± 0.3 ‰, respectively, but were distinctly different in 3 H, 5.5 ± 0.2 pCi/L and 13.7 ± 1.1 pCi/L, respectively. While δ 18 O shows that vegetation and meadow saturated zone water have similar origins, 3 H shows they have dissimilar ages.
Nitrogen (N) is an essential nutrient in soil that regulates plant growth, terrestrial sequestration of atmospheric carbon dioxide, and persistence of organic compounds. However, major knowledge gaps remain about how climate change may impact N accumulation and persistence, especially in deep soil and saprock (friable weakly weathered bedrock). Our objective was to understand how climate impacts the accumulation of N in soil and saprock and how climate impacts soil N distribution, persistence, and mechanisms of N persistence. We investigated N concentration in bulk soil and density fractions. We estimated N persistence along a bio-climatic sequence—sites range from a low-elevation oak savannah, mid-elevation pine-oak/mixed-conifer forest, to a high-elevation subalpine forest—in the southern Sierra Nevada in California. A combination of radiocarbon and elemental composition measurements along with a first-order kinetic model was used. The N concentration in the bulk soil and density fractions declined with depth, and there was a relatively greater mineral-associated heavy fraction (HF) N in deeper samples. The cooler/wetter mixed conifer site held 37% of profile N in saprock, which was greater than that of the entire soil profile at the drier/hotter oak savannah. The majority of N in soil, which was in the HF, was not influenced by climate proxies tested. However, both unprotected and occluded fractions of N were strongly influenced by climate. Soil N mean residence time (MRT) showed that drier/hotter climates have a shorter MRT, compared to mid-elevation sites with cooler/wetter climates. The effect of climate on deep saprock N storage might be indirect, primarily through climatic influence on the thickness of saprock. Overall, our findings suggest the mineral-associated HF N pool will not be vulnerable to changes in climate and will continue to contribute to the persistent soil N pool. The amount of topsoil and subsoil unprotected and occluded N can be explained by gross primary productivity and mean annual precipitation indicating that changes in climate can influence N partitioning. N stored in deep soil and saprock may be less vulnerable to climate than N stored in drier/hotter climates with less deeply stored N. It is critical to dig deeper to understand terrestrial ecosystems’ response to climate.
In this study, tungsten (W) occurrence and speciation was investigated in sediments collected from Fallon, Nevada where previous studies have linked elevated W levels in human body fluids to an unusual cluster of childhood leukemia cases. The speciation of sedimentary W was determined by μ-XRF mapping and μ-XANES. The W content of the analyzed surface sediments ranged between 81 and 25,908 mg/kg, which is significantly higher than the W content in deeper sediments which ranged from 37 to 373 mg/kg at 30 cm depth. The μ-XANES findings reveal that approximately 20–50% of the total W in the shallow sediment occurs in the metallic form (W 0 ); the rest occurs in the oxide form (W VI O 3 ). Because W 0 does not occur naturally, its elevated concentrations in surface sediments point toward a possible local anthropogenic origin. The oxidation of metallic W 0 with meteoric waters likely leads to the formation of W VI O 3 . The chief water-soluble W species was identified as WO 4 2– by chromatographic separation and speciation modeling. These results led us to postulate that W 0 particles from a currently unknown but local source(s) is (are) deposited onto the soils and/or surface sediments. The W 0 in interaction with meteoric water is oxidized to W VI O 3 , and as these sediment-water interactions progress, WO 4 2– is formed in the water at pH ~7. Under pH < 7, and sufficient W concentrations, tungstate tends to polymerize, and polymerized species are less likely to adsorb onto sediments. Polymerized species have lower affinity than monomers, which leads to enhanced mobility of W.
This dataset presents a 30-year high resolution meteorological dataset obtained using the WRF model (Advanced version Research WRF version 4.4). We used WRF and European Centre for Medium-Range Weather Forecasts Reanalysis v5 as initial and boundary conditions to generate gridded meteorological variables. A large number of surface weather stations was used for model validation. A multi-physics analysis was first developed to identify a good physics suite extended from 6 November 00 UTC to 10 November 23 UTC, 2018, which included the Camp Fire in northern California. Based on the best physics suite, the downscaling dataset extends from 1 December to 28 February, 1990–2021 and the horizontal domain has 1.5 km grid spacing covering the entire states of California and Nevada in the United States. Comparisons between hourly surface observations and WRF simulations of air temperature, relative humidity and wind speeds show mean absolute errors on the order of (1.6-2.0 C), (10 %) and 1.2–1.5 m s -1 , respectively.
In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.
This study characterises near-surface wind conditions at the Nevada Solar One (NSO) concentrated solar power plant using two years of measurements from a 15-m mast located 30 m west of the solar arrays. The surrounding terrain is fairly flat within 3 km, but nearby solar photovoltaic (PV) and concentrated solar power (CSP) installations introduce strong directional wind inhomogeneity, which motivates a sector-based analysis that distinguishes undisturbed inflow from waked flow. The analysis uses data from three sonic anemometers to assess whether turbulence in the surface layer can be represented as Gaussian and whether standard surface-layer scaling remains valid below 10 m height. The roughness length and displacement height are estimated for twelve directional sectors under near-neutral wind conditions, and all three velocity components show departures from Gaussian behaviour. Spectral analysis indicates that the turbulence contains low-frequency fluctuations with amplitudes higher than those predicted by standard models such as the Kaimal model or the Simiu and Scanlan model. The study introduces a new turbulence spectral model based on the eddy-surface-layer framework to better represent near-surface wind turbulence. The analysis further quantifies the influence of the surrounding power plants on mean wind profiles, one-point spectra, and vertical coherence for all velocity components. The solar collector arrays located 30 m east of the mast significantly modify the turbulence structure, whereas PV installations situated a few hundred metres away have only a limited influence on the turbulence measured at the mast. These findings highlight the limitations of conventional turbulence models for near-surface wind over solar infrastructure and provide improved parameterisations to support the structural design. The implications for the structural response of parabolic troughs and PV modules to such turbulence remain uncertain and require further investigation.
As historic drought conditions become more common in western North America, Late Quaternary hydroclimate records become vital for putting present anthropogenic conditions into a longer-term context. Here, we establish a high-resolution record of drought for the eastern Sierra Nevada (California) using lacustrine carbonates from well-dated sediment cores. We used oxygen and carbon stable-isotope ratios, combined with high-resolution scanning X-ray fluorescence counts of calcium (Ca) and titanium (Ti), to reconstruct the drought record over the last 4600 years in June Lake. We found elevated δ 18 O and δ 13 C carbonate isotope values coinciding with peaks in both total inorganic carbon and Ca/Ti, suggesting enhanced carbonate precipitation in response to evaporative concentration of lake water. At least six intervals of prolonged (centennial-scale) carbonate deposition were identified, including three pulses during the Late Holocene Dry Period (LHDP; ~3500–2000 cal yr BP), the Medieval Climate Anomaly (~1200–800 cal yr BP), and the Current Warm Period, which began around 100 cal yr BP. This record highlights the complexities of the LHDP, an interval that was more variable at June Lake than has been previously described in regional records.
Abstract Sub‐grid topographic heterogeneity has large impacts on surface energy balance and land‐atmosphere interactions. However, the impacts of representing sub‐grid topographic effects in land surface models (LSMs) on surface energy balance and boundary conditions remain unclear. This study analyzed and evaluated the impacts of sub‐grid topographic representations on surface energy balance, turbulent heat flux, and scalar (co‐)variances in the Energy Exascale Earth System Model (E3SM) land model (ELM). Three sub‐grid topographic representations in ELM were compared: (a) the default sub‐grid structure (D), (b) the recently developed sub‐grid topographic structure (T), and (c) high spatial resolution (1KM). Additionally, two different solar radiation schemes in ELM were compared: (a) the default plane‐parallel radiative transfer scheme (PP) and (b) the parameterization scheme (TOP) that accounts for sub‐grid topographic effects on solar radiation. A series of offline simulations with the three grid discretization structures (D, T, and 1KM) and two schemes of solar radiation (TOP and PP) were carried out using the Sierra Nevada, California. 1KM simulations with TOP well capture the spatial heterogeneity of surface fluxes compared to Moderate Resolution Imaging Spectroradiometer remote sensing data. There are significant differences between TOP and PP in the 1‐km simulated surface energy balance, but the differences in mean values and standard deviations become small when aggregated to the grid scale (i.e., 0.5°). The T configuration better mimics the 1KM simulations with TOP than the D configuration and better captures the sub‐grid topographic effects on surface energy balance and boundary conditions. These results underline the importance of representing sub‐grid topographic heterogeneities in LSMs and motivate future research to understand the sub‐grid topographic effects on land‐atmosphere interactions over mountainous areas.
Abstract Tectonic activity, geothermal fluids, and microseismic events (MSEs) tend to occur in similar locations as a result of spatiotemporal changes in the subsurface stress state. To quantify this association, we analyze data from a dense seismic array deployed at the San Emidio geothermal field, Nevada for 1 week in December 2016 to coincide with a 19.45‐hr shutdown of all injection and production pumping operations. 123 MSEs were detected, of which 101 occurred during the shutdown. The spatial association of the MSEs with the production wells suggests a causal relationship between the production cessation and the MSEs. Here we performed a detailed analysis to investigate reservoir material properties, distribution of seismically activated faults, and local stress state. We determined the hypocenters, magnitudes, and focal mechanisms for the MSEs, P‐wave tomographic velocity model, and local stress tensor. The results show that most MSEs occurred near the production wells. Magnitudes fall between −2.2 and 0.0 with larger events located closer to the production wells. Most MSEs occurred within a westward‐dipping normal fault zone in the reservoir associated with anomalously low P‐wave velocity values. The focal mechanism and stress inversion results show predominantly normal faulting with the maximum horizontal stress oriented north‐south. We suggest that the MSEs during shutdown were triggered on pre‐existing, small‐scale, critically stressed fault patches in the reservoir as the pore pressure increased around the production wells when the production pumping ceased. We interpret the larger MSE magnitudes closer to the production wells as a result of higher pore pressure increase.
Abstract Prior research confirmed the substantial bias from using precipitation‐based intensity‐duration‐frequency curves (PREC‐IDF) in design flood estimates and proposed next‐generation IDF curves (NG‐IDF) that represent both rainfall and snow processes in runoff generation. This study improves the NG‐IDF technology for a snow‐dominated test basin in the Sierra Nevada. A well‐validated physics‐based hydrologic model, the Distributed Hydrology Soil Vegetation Model (DHSVM), is used to continuously simulate snowmelt and streamflow that are used as benchmark data sets to systematically assess the NG‐IDF technology. We find that, for the studied small snow‐dominated basin, the use of standard rainfall hyetographs in the NG‐IDF technology leads to substantial underestimation of design floods. Thus, we propose probabilistic hyetographs that can represent unique patterns of events with different underlying mechanisms. For the test basin where flooding events are generated entirely by snowmelt, we develop a hyetograph that characterizes snowmelt temporal patterns, which greatly improves the performance of NG‐IDF technology in design flood estimates. In contrast to the standard rainfall hyetographs characterized by a symmetrically peaked, bell‐shaped curve, the snowmelt hyetograph displays a more rapid rise (i.e., greater intensity) and a distinct diurnal pattern influenced by solar energy input. The results also show that the uncertainty of hyetography plays an important role in design flood estimation and can have important implications for future flood projections.
Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.
Large uncertainty remains in the spatial distribution of deep soil organic carbon (OC) storage and how climate controls belowground OC. This research aims to quantify OC stocks, characterize soil OC age and chemical composition, and evaluate climatic impacts on OC storage from the soil surface through the deep critical zone to bedrock. These objectives were carried out at four sites along a bio-climosequence in the Sierra Nevada, California. On average, 74% of OC was stored below the A horizon, and up to 30% of OC was stored in saprock (friable weakly weathered bedrock). Radiocarbon, spectroscopic, and isotopic analyses revealed the coexistence of very old organic matter (OM) (mean radiocarbon age = 20 300 years) with relatively recent OM (mean radiocarbon age = 4800 years) and highly decomposed organic compounds with relatively less decomposed material in deep soil and saprock. This co-mingling of OM suggests OC is prone to both active cycling and long-term protection from degradation. In addition to having direct effects on OC cycling, climate indirectly controls deep OC storage through its impact on the degree of regolith weathering (e.g. thickening). Although deep OC concentrations are low relative to soil, thick saprock represents a large, previously unrealized OC pool.
Three-dimensionally preserved Ediacaran fossils occur globally within sandstone beds. Sandy siliciclastic deposits of the Ediacaran Wood Canyon Formation (WCF) in the Montgomery Mountains, Nevada, contain two fossil morphologies with similar shapes and sizes: one exhibits mm-scale ridges and a distinct lower boundary and the other is devoid of these diagnostic features. We interpret these as taphomorphs of erniettomorphs, soft-bodied organisms with uncertain taxonomic affinities. We explore the cast-and-mould preservation of both taphomorphs by petrography, Raman spectroscopy, X-ray fluorescence microprobe and X-ray diffraction. All fossils and the surrounding sedimentary matrix contain quartz grains, iron-rich chlorite and muscovite. The ridged fossils contain about 70% larger quartz grains compared to the ridgeless taphomorph, indicating a lower abundance of clay minerals in the ridged fossil. Chlorite and muscovite likely originated from smectite and kaolinite precursors that underwent lower greenschist facies metamorphism. Kaolinite and smectite are inferred to have been abundant in sediments around the ridged fossil, which enabled the preservation of a continuous, distinct, clay- and kerogen-rich bottom boundary. The prevalence of quartz in the ridged fossils of the WCF and in erniettomorphs from other localities also suggests a role for this mineral in three-dimensional preservation of erniettomorphs in sandstone and siltstone deposits.
Abstract In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid-flow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fluid-flow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k -means clustering (NMF k ), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.
In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.
This submission includes a high-precision seismic event catalog estimated from seismic data collected at San Emidio, Nevada from April to May 2022. The catalog lists the precise time, location, and magnitude of microseismic events recorded during this period. Both the seismic data and microseismic event catalog were produced as part of the Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses (WHOLESCALE) project. Attached here are the microseismic event catalog, a link to the raw seismic data, and a detailed description of methods used to create the catalog.