SEARCH · Search NASA

Results for “Surface topography”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2021

ESS-DIVE Reporting Format for Hydrologic Monitoring Data and Metadata

The ESS-DIVE Reporting Format (RF) for Hydrologic Monitoring Data and Metadata was designed for water level, temperature, electrical conductivity, specific conductance, dissolved oxygen, and pH measured in situ using sensors deployed in water bodies. It can also be applied to a broader suite of variables beyond its original scope. The RF is intended to provide clear and easy guidance to data generators aiming to upload consistent findable, accessible, interoperable, and reusable (FAIR) datasets. The RF includes instructions (pdf), a term guide (csv), templates (csv and xlsx), and a list of recommended vocabulary (csv) to equip the data generator with all of the tools needed to successfully apply the RF. There is also a crosswalk to some key data resources and standards (csv). Users should cite this dataset when using the RF but are encouraged to access the RF files on GitHub (https://github.com/ess-dive-community/essdive-hydrologic-monitoring) where the RF documentation is maintained and updated on the ESS-DIVE Community Space GitHub. Questions and suggestions for improvement can be submitted as GitHub issues.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Factors Controlling a Synthetic Aperture Radar (SAR) Derived Root-Zone Soil Moisture Product over The Seward Peninsula of Alaska

Root-zone soil moisture exerts a fundamental control on vegetation, energy balance, and the carbon cycle in Arctic ecosystems, but it is still not well understood in vast, remote, and understudied regions of discontinuous permafrost. The root-zone soil moisture product (30 m resolution) used in this analysis was retrieved from a time-series P-Band (420–440 MHz) synthetic aperture radar (SAR) backscatter observations (August 2017 & October 2017). While similar approaches have been taken to retrieve surface (0 cm to 5 cm) soil moisture from L-Band (1.2 GHz) SAR backscatter, this is one of the first known attempts at reaching the root-zone in permafrost regions. Here, we analyze secondary factors (excluding primary factors, such as precipitation) controlling summer (August) soil moisture at depths of 6 cm, 12 cm, and 20 cm over a 4500 km2 area on the Seward Peninsula of Alaska. Using a random forest model, we quantify the impact of topography, vegetation, and meteorological factors on soil moisture distributions. In developing the random forest model, we explore a variety of feature scales (30 m, 60 m, 90 m, 120 m, 180 m, and 240 m), tune hyperparameters (the structure of individual decision trees making up the ensemble including the number and depth of trees), and perform the final feature selection using cross-validated recursive feature elimination. Results suggest that root-zone soil moisture on the Seward Peninsula is primarily controlled by vegetation at 6 cm, but deeper in the soil column topography and meteorological factors, such as predominant winter wind direction and summer insolation, play a larger role. The random forest model accounts for 40% to 60% of the variation observed (R2 = 0.44 at 6 cm, R2 = 0.52 at 12 cm, R2 = 0.58 at 20 cm). These results indicate that vegetation is the dominant control on soil moisture shallow in the soil column, but the impact of vegetation does not extend to deeper layers retrieved from P-Band SAR backscatter.

Dann, Julian↗

DOE OSTI2023

Amazon drought resilience – emerging results point to new empirical needs

Will climate change turn tropical forests from a carbon sink to a source (Pan et al., 2011)? Tropical forests cycle more carbon, water, and energy than any other biome (Bonan, 2008). The Amazon is the largest tropical forest, spanning one-third of South America, and any perturbation in its function has repercussions on the global climate. Mild and severe droughts in the Amazon are predicted to double and triple in area, respectively, by 2100, while the area under wet extremes will increase (Duffy et al., 2015). Our understanding of the role of plants' hydrological environments in determining their drought response, however, is limited (Chitra-Tarak et al., 2018, 2021). An important review by Costa et al. (2023; pp. 714–733) published in this issue of New Phytologist highlights that shallow-water table (WT) forests constitute c. 50% of the Amazon and may act as hydrological refugia during droughts, yet they have been neglected in Amazon forest research to date. Emerging results suggest that the shallow-WT forests that dominate the Amazon basin may increase in productivity under moderate droughts due to relief from hypoxia, suggesting a potential buffer from drought relative to deeper WT forests (Fig. 1). During severe droughts, however, shallow-WT forests may be vulnerable to collapse due to drought-intolerant traits. Addressing the underrepresentation of relatively drought-resilient shallow-WT forests in the Amazon's forest inventories may help resolve the much-debated incongruence in recent studies of the Amazon's drought resilience. Data from forest inventories (biased toward deep-WT forests) indicate that the Amazon's carbon sink is declining, partially because of drought impacts (Phillips et al., 2009; Brienen et al., 2015), whereas basin-wide satellite-based measures of gross primary productivity, which include shallow-WT forests, indicate varied regional responses to droughts (Saleska et al., 2007; Brando et al., 2010). Costa et al.'s (2023) descriptive insights of the intricate balance of geology, topography, hydrology, vegetation, and drought on ecosystem function, as well as their conceptual predictive framework, are useful for developing new empirical research in these understudied ecosystems and improving Earth system models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Modeling the topographic influence on aboveground biomass using a coupled model of hillslope hydrology and ecosystem dynamics

Abstract. Topographic heterogeneity and lateral subsurface flow at the hillslope scale of ≤1 km may have outsized impacts on tropical forest through their impacts on water available to plants under water-stressed conditions. However, vegetation dynamics and finer-scale hydrologic processes are not concurrently represented in Earth system models. In this study, we integrate the Energy Exascale Earth System Model (E3SM) land model (ELM) that includes the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), with a three-dimensional hydrology model (ParFlow) to explicitly resolve hillslope topography and subsurface flow and perform numerical experiments to understand how hillslope-scale hydrologic processes modulate vegetation along water availability gradients at Barro Colorado Island (BCI), Panama. Our simulations show that groundwater table depth (WTD) can play a large role in governing aboveground biomass (AGB) when drought-induced tree mortality is triggered by hydraulic failure. Analyzing the simulations using random forest (RF) models, we find that the domain-wide simulated AGB and WTD can be well predicted by static topographic attributes, including surface elevation, slope, and convexity, and adding soil moisture or groundwater table depth as predictors further improves the RF models. Different model representations of mortality due to hydraulic failure can change the dominant topographic driver for the simulated AGB. Contrary to the simulations, the observed AGB in the well-drained 50 ha forest census plot within BCI cannot be well predicted by the RF models using topographic attributes and observed soil moisture as predictors, suggesting other factors such as nutrient status may have a larger influence on the observed AGB. The new coupled model may be useful for understanding the diverse impact of local heterogeneity by isolating the water availability and nutrient availability from the other external and internal factors in ecosystem modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Scaling Arctic landscape and permafrost features improves active layer depth modeling

Tundra ecosystems in the Arctic store up to 40% of global below-ground organic carbon but are exposed to the fastest climate warming on Earth. However, accurately monitoring landscape changes in the Arctic is challenging due to the complex interactions among permafrost, micro-topography, climate, vegetation, and disturbance. This complexity results in high spatiotemporal variability in permafrost distribution and active layer depth (ALD). Moreover, these key tundra processes interact at different scales, and an observational mismatch can limit our understanding of intrinsic connections and dynamics between above and below-ground processes. Consequently, this could limit our ability to model and anticipate how ALD will respond to climate change and disturbances across tundra ecosystems. In this paper, we studied the fine-scale heterogeneity of ALD and its connections with land surface characteristics across spatial and spectral scales using a combination of ground, unoccupied aerial system, airborne, and satellite observations. We showed that airborne sensors such as AVIRIS-NG and medium-resolution satellite Earth observation systems like Sentinel-2 can capture the average ALD at the landscape scale. We found that the best observational scale for ALD modeling is heavily influenced by the vegetation and landform patterns occurring on the landscape. Landscapes characterized by small-scale permafrost features such as polygon tussock tundra require high-resolution observations to capture the intrinsic connections between permafrost and small-scale land surface and disturbance patterns. Conversely, in landscapes dominated by water tracks and shrubs, permafrost features manifest at a larger scale and our model results indicate the best performance at medium resolution (5 m), outperforming both higher (0.4 m) and lower resolution (10 m) models. This transcends our study to show that permafrost response to climate change may vary across dominant ecosystem types, driven by different above- and below-ground connections and the scales at which these connections are happening. We thus recommend tailoring observational scales based on landforms and characteristics for modeling permafrost distribution, thereby mitigating the influences of spatial-scale mismatches and improving the understanding of vegetation and permafrost changes for the Arctic region.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

DOE OSTI2023

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Arctic tundra shrubification: a review of mechanisms and impacts on ecosystem carbon balance

Abstract Vegetation composition shifts, and in particular, shrub expansion across the Arctic tundra are some of the most important and widely observed responses of high-latitude ecosystems to rapid climate warming. These changes in vegetation potentially alter ecosystem carbon balances by affecting a complex set of soil–plant–atmosphere interactions. In this review, we synthesize the literature on (a) observed shrub expansion, (b) key climatic and environmental controls and mechanisms that affect shrub expansion, (c) impacts of shrub expansion on ecosystem carbon balance, and (d) research gaps and future directions to improve process representations in land models. A broad range of evidence, including in-situ observations, warming experiments, and remotely sensed vegetation indices have shown increases in growth and abundance of woody plants, particularly tall deciduous shrubs, and advancing shrublines across the circumpolar Arctic. This recent shrub expansion is affected by several interacting factors including climate warming, accelerated nutrient cycling, changing disturbance regimes, and local variation in topography and hydrology. Under warmer conditions, tall deciduous shrubs can be more competitive than other plant functional types in tundra ecosystems because of their taller maximum canopy heights and often dense canopy structure. Competitive abilities of tall deciduous shrubs vs herbaceous plants are also controlled by variation in traits that affect carbon and nutrient investments and retention strategies in leaves, stems, and roots. Overall, shrub expansion may affect tundra carbon balances by enhancing ecosystem carbon uptake and altering ecosystem respiration, and through complex feedback mechanisms that affect snowpack dynamics, permafrost degradation, surface energy balance, and litter inputs. Observed and projected tall deciduous shrub expansion and the subsequent effects on surface energy and carbon balances may alter feedbacks to the climate system. Land models, including those integrated in Earth System Models, need to account for differences in plant traits that control competitive interactions to accurately predict decadal- to centennial-scale tundra vegetation and carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

The First Detection of an Earthquake From a Balloon Using Its Acoustic Signature

Extreme temperature and pressure conditions on the surface of Venus present formidable technological challenges against performing ground-based seismology. Efficient coupling between the Venusian atmosphere and the solid planet theoretically allows the study of seismically generated acoustic waves using balloons in the upper atmosphere, where conditions are far more clement. However, earthquake detection from a balloon has never been demonstrated. We present the first detection of an earthquake from a balloon-borne microbarometer near Ridgecrest, CA in July 2019 and include a detailed analysis of the dependence of seismic infrasound, as measured from a balloon on earthquake source parameters, topography, and crustal and atmospheric structure. Our comprehensive analysis of seismo-acoustic phenomenology demonstrates that seismic activity is detectable from a high-altitude platform on Earth, and that Rayleigh wave-induced infrasound can be used to constrain subsurface velocities, paving the way for the detection and characterization of such signals on Venus.

58 GEOSCIENCES↗

DOE OSTI2025

The Anisotropic Adsorption of De Novo Allosteric Two‐Component Protein Fibers on Mica Surfaces

Protein adsorption at solid–liquid interfaces underlies many biomedical and materials applications, yet the mechanisms governing adsorption of proteins and their assemblies remain poorly understood. Here we investigate de novo–designed proteins that self-assemble into three fibrillar morphologies — small (S), large (L), and helical (H)—on muscovite mica. S-fibers are metastable, forming first but diminishing as L- and H-fibers develop and deposit. Adsorption of both fibers and monomers depends on fiber morphology and solvent environment, which is modulated by the substrate. The anisotropic surface features of the fiber types—long axis of S- and L-fibers and helical grooves of H-fibers—correlate with specific crystallographic directions on mica (001). S- and L-fibers align along the unique lattice axis, while the substrate-facing groove of H-fibers preferentially aligns along the remaining symmetry-related directions. Increasing potassium chloride (KCl) concentrations to molar levels alters adsorption, decreasing monomer coverage relative to fibers. These observations, interpreted through the effect of interfacial interface, indicate that ion-mediated modulation of hydration layers governs electrostatic interactions and alignment. This study reveals how coupling between protein topography, substrate crystallography, and hydration structure dictates selective adsorption and orientation of protein assemblies, offering insights for rational biomolecular material design.

AFM↗

DOE OSTI2021

Operando Topography and Mechanical Property Mapping of CO 2 Reduction Gas-Diffusion Electrodes Operating at High Current Densities

Electrochemical atomic force microscopy (EC-AFM) enables measurement of electrode topography and mechanical properties during electrochemical reactions. However, for aqueous-based reactions that make gas products, such as CO 2 reduction and water splitting into CO/H 2 , current densities below 1 mA cm −2 have been necessary to prevent formation of bubbles at the electrode; such bubbles can stick to the AFM probe and prevent further AFM imaging. Here, we demonstrate a novel cell design with a gas-diffusion electrode (GDE) to exhaust the gas products, thereby enabling high current density EC-AFM measurements at 1, 10, and 100 mA cm −2 that are not disturbed by bubble formation at the electrode surface. These experiments revealed a stable morphological structure of Cu catalysts deposited on GDEs during high current density operation. Systematic spatially resolved maps of deformation and adhesion showed no signs of a gas-liquid interface between catalyst particles of the GDE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Impact of Coastal Marsh Eco‐Geomorphologic Change on Saltwater Intrusion Under Future Sea Level Rise

Abstract Coastal saltwater intrusion (SWI) is one key factor that affects the hydrology, ecology, and biogeochemistry of coastal ecosystems. Future climate change, especially intensified sea level rise (SLR), is expected to trigger SWI to encroach on coastal freshwater aquifers more intensively. Numerous studies have investigated decadal/century scale SWI under SLR by assuming a static coastal landscape topography. However, coastal landscapes are highly dynamic in response to SLR, and the impact of coastal landscape evolution on SWI has received very little attention. Thus, this study used a coastal marsh landscape as an example and investigated how coastal marsh evolution affects future SWI with a physically‐based coastal hydro‐eco‐geomorphologic model, Advanced Terrestrial Simulator. Our numerical experiments showed that it is very likely that the marsh elevation increases with future SLR due to sediment deposition, and a depression zone is formed due to different marsh accretion rates between the ocean boundary and the inland. We found that marsh accretion may significantly reduce the surface saltwater inflow at the ocean boundary, and the evolved topographic depression zone may prolong the residence time of surface ponded saltwater, affecting subsurface salinity distribution differently. We also predicted that marshlands may become more sensitive to upland freshwater supply under future SLR, compared with previous predictions without marsh evolution. This study demonstrates the importance of coastal evolution to coastal freshwater‐saltwater interaction. The eco‐geomorphologic effect may not be ignored when evaluating coastal SWI under SLR at decadal or century scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Mechanisms influencing and prediction of tool influence function spots during hemispherical sub-aperture tool polishing on fused silica

Sub-aperture tool polishing of precision optics requires a detailed understanding of the local material removal [tool influence function (TIF)] at the contact spot between the workpiece and tool to achieve high removal determinism and hence precision of the optic relative to the desired/design surface figure. In this study, the mechanisms influencing and the quantitative prediction of the removal rate and shape of TIF spots during polishing of fused silica glass with cerium oxide slurry using a rotating hemispherical pad-foam tool for a wide variety of process conditions (including tool properties, kinematics, and applied displacements) are investigated. The TIF volumetric removal rate can be estimated utilizing the average relative velocity and contact area using a simple analytical model. In addition, stability of the volumetric removal rate for fixed process conditions is shown to be greatly dependent on the pad preparation and amount of tool use (affecting both pad topography and slurry buildup), whose general behavior shows an increase in removal rate followed by stabilization with polishing time. The determination of the TIF removal shape is more complex. An extended version of the Preston removal model is developed to explain a comprehensive set of measured TIF removal shapes to within ∼ <#comment/> 22 % <#comment/> . This model incorporates a number of phenomena impacting the TIF removal shape including: (a) temporal and spatial dependent relative velocity between the workpiece and tool; (b) an elastic mechanics based, as well as hydrodynamic, pressure distribution; (c) a spatially dependent friction coefficient possibly caused by both reduced slurry replenishment in low velocity regions and pad slurry islands (100 µm scale) and porosity (millimeter scale); and (d) a shear-based removal mechanism on the periphery of the contact spot.

Suratwala, T. (ORCID:0000000190861039)↗

DOE OSTI2022

A novel hardmask-to-substrate pattern transfer method for creating 3D, multi-level, hierarchical, high aspect-ratio structures for applications in microfluidics and cooling technologies

Abstract This letter solves a major hurdle that mars photolithography-based fabrication of micro-mesoscale structures in silicon. Conventional photolithography is usually performed on smooth, flat wafer surfaces to lay a 2D design and subsequently etch it to create single-level features. It is, however, unable to process non-flat surfaces or already etched wafers and create more than one level in the structure. In this study, we have described a novel cleanroom-based process flow that allows for easy creation of such multi-level, hierarchical 3D structures in a substrate. This is achieved by introducing an ultra-thin sacrificial silicon dioxide hardmask layer on the substrate which is first 3D patterned via multiple rounds of lithography. This 3D pattern is then scaled vertically by a factor of 200–300 and transferred to the substrate underneath via a single shot deep etching step. The proposed method is also easily characterizable—using features of different topographies and dimensions, the etch rates and selectivities were quantified; this characterization information was later used while fabricating specific target structures. Furthermore, this study comprehensively compares the novel pattern transfer technique to already existing methods of creating multi-level structures, like grayscale lithography and chip stacking. The proposed process was found to be cheaper, faster, and easier to standardize compared to other methods—this made the overall process more reliable and repeatable. We hope it will encourage more research into hybrid structures that hold the key to dramatic performance improvements in several micro-mesoscale devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Exploratory Efforts to Constrain Geologic Material Properties from Remote Sensing Data: Joint Study (FY2020 Final Report)

Identification and characterization of underground events from surface or remote data requires a thorough understanding of the rock material properties. However, material properties usually come from borehole data, which is expensive and not always available. A potential alternative is to use topographic characteristics to approximate the strength, but this has never been done before quantitatively. Here we present the results from the first steps towards this goal. We have found that there are strong correlations between compressive and tensile strengths and slopes, but these correlations vary depending on data analysis details. Rugosity may be better correlated to strength than slope values. More comprehensive analyses are needed to fully understand the best method of predicting strength from topography for this area. We also found that misalignment of multiple GIS datasets can have a large influence on the ability to make interpretations. Lastly, these results will require further study in a variety of climatic conditions before being applicable to other sites.

58 GEOSCIENCES↗

DOE OSTI2021

Exploratory Efforts to Constrain Geologic Material Properties from Remote Sensing Data: A Joint Study

42 ENGINEERING↗

DOE OSTI2023

Watersheds and stream networks viewed longitudinally: Example insights from novel spatial portrayals of watershed characteristics

Abstract Longitudinal depictions of watershed structure and characteristics, including topography, stream networks, wetlands, ground water levels, and land use, can provide watershed knowledge and understanding unavailable from standard plan view maps. Three case studies provide examples of knowledge gained by applying longitudinal views of stream networks, watershed hydrologic behavior, and land use distributions. Longitudinal views of mountain stream networks show extreme variability in the slope‐area relationships of low Strahler order streams, large discontinuities in drainage area (large parts of drainage area space are absent in networks), and large variations in network curvature. Longitudinal views of a groundwater‐dominated headwater watershed increase the inference available from limited groundwater observations and clearly reveal how groundwater connections affect the permanence of surface water features and the distribution of vadose zone storage in the landscape. Plotting land uses longitudinally illuminates and allows a quantitative analysis of how land uses are distributed relative to topographic position. Viewing watersheds and stream networks longitudinally can provide new insights into watershed forms and processes and motivate new questions and research.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗