Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface hydrology”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Past fires and post-fire impacts reconstructed from a southwest Australian stalagmite

Speleothem records of past environmental change provide an important opportunity to explore fire frequency and intensity in the past, and the antecedent climatic conditions leading to fire events. Here, fire sensitive geochemical signals in a stalagmite from Yonderup Cave, a shallow cave in Western Australia, are compared to well-documented wildfire events that occurred in recent decades. The results are extended to identify wildfires during the growth interval of the speleothem (1760 CE–2005 CE). Principal component analysis of the stalagmite time-series revealed distinct peaks in a combination of phosphorus and metal (aluminium, zinc, copper and lead) concentrations in response to known fire events, which are interpreted to have come from ash. Varying responses in the geochemical signal in the pre- and post-European colonisation period are likely linked to changes in land management, fire frequency and fire intensities. Of note is what we infer to be a particularly intense fire event concluded to have occurred in 1897 ± 5 CE, which climate sensitive proxies (δ18O and bedrock-derived elements) indicate was preceded by a multi-decadal dry period that began in the late 1860s. The intensity of the deduced fire event is supported by a peak in P that is 6.3 times higher than elsewhere in the speleothem, accompanied by a peak in Zn. It was also inferred that changes to the surface-cave hydrology occurred as a result of this event, indicated by increased input of colloidal organic material onto the stalagmite that was interpreted to be caused by shallow karst bedrock fracturing from the heat-induced deformation that is observed after intense fires. These findings show the potential for speleothems to provide dated records of fire intensity and recurrence intervals. Further development could lead to a better understanding of the climate-fire relationship and the effects of land-management practices on wildfire frequency and intensity.

58 GEOSCIENCES↗
DOE OSTI2020

Evapotranspiration depletes groundwater under warming over the contiguous United States

A warmer climate increases evaporative demand. However, response to warming depends on water availability. Existing earth system models represent soil moisture but simplify groundwater connections, a primary control on soil moisture. Here we apply an integrated surface-groundwater hydrologic model to evaluate the sensitivity of shallow groundwater to warming across the majority of the US. We show that as warming shifts the balance between water supply and demand, shallow groundwater storage can buffer plant water stress; but only where shallow groundwater connections are present, and not indefinitely. As warming persists, storage can be depleted and connections lost. Similarly, in the arid western US warming does not result in significant groundwater changes because this area is already largely water limited. The direct response of shallow groundwater storage to warming demonstrates the strong and early effect that low to moderate warming may have on groundwater storage and evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global gridded crop harvested area, production, yield, and monthly physical area data circa 2015

Abstract Here we provide an update to global gridded annual and monthly crop datasets. This new dataset uses the crop categories established by the Global Agro-Ecological Zones (GAEZ) Version 3 model, which is based on the Food and Agricultural Organization of the United Nations (FAO) crop production data. We used publicly available data from the FAOSTAT database as well as GAEZ Version 4 global gridded dataset to generate circa 2015 annual crop harvested area, production, and yields by crop production system (irrigated and rainfed) for 26 crops and crop categories globally at 5-minute resolution. We additionally used available data on crop rotations, cropping intensity, and planting and harvest dates to generate monthly gridded cropland data for physical areas for the 26 crops by production system. These data are in standard georeferenced gridded format, and can be used by any global hydrology, land surface, or other earth system model that requires gridded annual or monthly crop data inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving the Estimation of the Atmospheric Water Vapor Pressure Using Interpretable Long Short-Term Memory Networks: Dataset, Python code, and trained models

Atmospheric water vapor pressure is an essential meteorological control on land surface and hydrologic processes. It is not as frequently observed as other meteorologic conditions, but often inferred through the August–Roche–Magnus formula by simply assuming dew point and daily minimum temperatures are equivalent or by empirically correlating the two temperatures using an aridity correction. The performance of both methods varies considerably across different regions and during different time periods; obtaining consistently accurate estimates across space and time remains a great challenge. We applied an interpretable Long Short-Term Memory (iLSTM) network conditioned on static, location specific attributes to estimate daily vapor pressure for 83 FLUXNET sites in the United States and Canada. This data package includes all raw data of the 83 FLUXNET sites, input data for model training/validation/test, trained models and results, and python codes for the manuscript "Improving the Estimation of the Atmospheric Water Vapor Pressure Using an Interpretable Long Short-term Memory Network". Specifically, it consists of five parts. - First, "1_Daymet_data_83sites.zip" includes raw data downloaded from Daymet for the 83 sites used in the paper according to their longitude and latitude, in which vapor pressure is used. It also includes a pre-processed CSV data file combining all data from the 83 sites which is specifically used for the paper. - Second, "2_Fluxnet2015_data_83sites.zip" includes raw half hourly data of the 83 sites downloaded from FLUXNET2015 data portal, pre-processed daily data of the 83 sites, a CSV file including combined pre-processed daily data of the 83 sites, and a CSV file including the information (site ID, site name, latitude, longitude, data available period) of the 83 sites. - Third, "3_MODIS_LAI_data_83sites_raw.zip" includes raw leaf area index (LAI) data downloaded from the AppEEARs data portal. - Fourth, "4_Scripts.zip" includes all scripts related to model training and post-processing of a trained model, and a jupyter notebook showing an example for model post-processing. Two typo errors in files titled "run2get_args.py" and "postprocess.py" were corrected on March 27, 2024 to avoid confusions. - Finally, "Trained_models_and_results.zip" includes three folders and three files with suffix ".npy", and each folder corresponds to one file with suffix ".npy" with the same title. Each of the three folders include all trained models associated with one iLSTM model configuration (35 models for each configuration, details are described in the paper). Each file with suffix ".npy" includes the post-processed results of the corresponding 35 models under one iLSTM model configuration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Representing farmer irrigated crop area adaptation in a large-scale hydrological model

Abstract. Large-scale hydrological models (LHMs) are commonly used for regional and global assessment of future water shortage outcomes under climate and socioeconomic scenarios. The irrigation of croplands, which accounts for the lion's share of human water consumption, is critical in understanding these water shortage trajectories. Despite irrigation's defining role, LHM frameworks typically impose trajectories of land use that underlie irrigation demand, neglecting potential dynamic feedbacks in the form of human instigation of and subsequent adaptation to water shortages via irrigated crop area changes. We extend an LHM, MOSART-WM, with adaptive farmer agents, applying the model to the continental United States to explore water shortage outcomes that emerge from the interplay between hydrologic-driven surface water availability, reservoir management, and farmer irrigated crop area adaptation. The extended modeling framework is used to conduct a hypothetical computational experiment comparing differences between a model run with and without the incorporation of adaptive farmer agents. These comparative simulations reveal that accounting for farmer adaptation via irrigated crop area changes substantially alters modeled water shortage outcomes, with US-wide annual water shortages being reduced by as much as 42 % when comparing adaptive and non-adaptive versions of the model forced with US climatology from the period 1950–2009.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗
DOE OSTI2024

Managed aquifer recharge as a strategy to redistribute excess surface flow to baseflow in snowmelt hydrologic regimes

Water management in snowmelt hydrologic regimes, characterized by large annual fluctuations in stream flow driven by seasonal snow melt, faces the challenge of highly variable supply that often does not align with timing of demand. Climate change may exacerbate management challenges by significantly reducing snowpack or shifting snow melt earlier. Here, managed aquifer recharge (MAR) is evaluated as a potential strategy to reallocate excess early-season stream flow to time periods when less surface water is available. This strategy differs from traditional MAR, where the goal is to minimize loss to surface water. We assess how to site MAR operations such that groundwater recharge flows back to the surface water system in a lagged manor to benefit water management objectives, which we term “enhanced baseflow.” We use a regional groundwater model for the Treasure Valley aquifer located in southwestern Idaho, United States to demonstrate a generalizable approach using regional groundwater models as tools to identify favorable baseflow enhancement locations. Hypothetical MAR is simulated at 197 candidate locations, which are then evaluated for how effectively they meet potential management objectives. In addition to demonstrating the modeling and evaluation approach, we discuss lessons learned from applying a pre-existing regional groundwater model to MAR for enhanced baseflow and also describe important considerations, such as the physical and institutional availability of surface flows and specific management objectives, when assessing regional and site-specific suitability of MAR for enhanced baseflow as a potential management strategy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Estimating Future Surface Water Availability Through an Integrated Climate‐Hydrology‐Management Modeling Framework at a Basin Scale Under CMIP6 Scenarios

Abstract Climate change and increasing water demand due to population growth pose serious threats to surface water availability. The biggest challenge in addressing these threats is the gap between climate science and water management practices. Local water planning often lacks the integration of climate change information, especially with regard to its impacts on surface water storage and evaporation as well as the associated uncertainties. Using Texas as an example, state and regional water planning relies on the use of reservoir “Firm Yield” (FY)—an important metric that quantifies surface water availability. However, this existing planning methodology does not account for the impacts of climate change on future inflows and on reservoir evaporation. To bridge this knowledge gap, an integrated climate‐hydrology‐management (CHM) modeling framework was developed, which is generally applicable to river basins with geographical, hydrological, and water right settings similar to those in Texas. The framework leverages the advantages of two modeling approaches—the Distributed Hydrology Soil Vegetation Model (DHSVM) and Water Availability Modeling (WAM). Additionally, the Double Bias Correction Constructed Analogues method is utilized to downscale and incorporate Coupled Model Intercomparison Project Phase 6 GCMs. Finally, the DHSVM simulated naturalized streamflow and reservoir evaporation rate are input to WAM to simulate reservoir FY. A new term—“Ratio of Firm Yield” (RFY)—is created to compare how much FY changes under different climate scenarios. The results indicate that climate change has a significant impact on surface water availability by increasing reservoir evaporation, altering the seasonal pattern of naturalized streamflow, and reducing FY.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Water balance model (WBM) v.1.0.0: a scalable gridded global hydrologic model with water-tracking functionality

Abstract. This paper describes the University of New Hampshire Water Balance Model, WBM, a process-based gridded global hydrologic model that simulates the land surface components of the global water cycle and includes water extraction for use in agriculture and domestic sectors. The WBM was first published in 1989; here, we describe the first fully open-source WBM version (v.1.0.0). Earlier descriptions of WBM methods provide the foundation for the most recent model version that is detailed here. We present an overview of the model functionality, utility, and evaluation of simulated global river discharge and irrigation water use. This new version adds a novel suite of water source tracking modules that enable the analysis of flow-path histories on water supply. A key feature of WBM v.1.0.0 is the ability to identify the partitioning of sources for each stock or flux within the model. Three different categories of tracking are available: (1) primary inputs of water to the surface of the terrestrial hydrologic cycle (liquid precipitation, snowmelt, glacier melt, and unsustainable groundwater); (2) water that has been extracted for human use and returned to the terrestrial hydrologic system; and (3) runoff originating from user-defined spatial land units. Such component tracking provides a more fully transparent model in that users can identify the underlying mechanisms generating the simulated behavior. We find that WBM v.1.0.0 simulates global river discharge and irrigation water withdrawals well, even with default parameter settings, and for the first time, we are able to show how the simulation arrives at these fluxes by using the novel tracking functions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Temperature, Humidity, and Time-Lapse Video Data from Yosemite National Park, Water Year 2024

This dataset contains time-lapse imagery and distributed measurements of air temperature, relative humidity, dew point, and soil temperature across Yosemite National Park from October 2023 to September 2024. Instruments were deployed at 10 sites in two cross-valley transects as part of the DOE Grant: Seasonal Cycles Unravel Mysteries of Missing Mountain Water organized by Jessica Lundquist (University of Washington), Rosemary Carroll (Desert Research Institute), and Ethan Gutmann (National Center for Atmospheric Research). The data are intended to support hydrologic modeling efforts to better resolve the fate of mountain water, and are published to support studies of surface climate or hydrologic processes in complex terrain. Measurements were collected with low-cost data loggers installed 2 m high on evergreen trees or buried just below the soil surface. A time-lapse camera at one site captures valley-scale seasonal snow cover variability.Dataset files are organized by site and variable (air measurements, ground measurements, or time-lapse video). Air and ground measurements are packaged in LoggerData.zip, and time-lapse imagery is compiled into a short video stored in TimelapseVideos.zip. File-level metadata contains details for each file included in the dataset. A data dictionary provides units and descriptions for column or row names in all files. The locations metadata file describes site characteristics, locations, and associated GPS methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spatial Study 2022: Surface Water Samples, Cotton Strip Degradation, and Hydrologic Sensor Data across the Yakima River Basin, Washington, USA (v3)

This dataset supports a broader study examining the drivers of spatial variability in sediment respiration rates in the Yakima River Basin. The dataset provides data and photos generated from sample collection during the same one-week period at 48 sites within multiple rivers throughout the Yakima River Basin in Washington, USA. The contents include surface water geochemistry data; river substrate grain size photos; stream depth data; manual chamber open channel respiration data; and field metadata (including qualitative information on instream and river corridor characteristics). Grain size photos can be used to improve estimates of channel substrate D50 data. The dataset also includes tensile strength and photos from cotton strip field degradation experiments; five-week sensor time series temperature, dissolved oxygen, pressure, pH, specific conductance, chlorophyll A, and turbidity data; plots of the sensor data; and R scripts used to generate the plots. Samples collected during this study were labeled as “Second Spatial Study” or “SSS.” A subset of data from the SSS samples were published in the contiguous United States (CONUS)-Scale Model-Sample (CM) study data package available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1923689 that presents data from across the CONUS. SSS data published in the CM data package were not included in this data package. They include dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC), total nitrogen (TN), grain size, aerobic sediment respiration, dissolved oxygen (DO), and temperature. Parent IDs and Site IDs are consistent between the SSS and CM data packages, and they can be mapped directly so data across packages can be used together. Additionally, sensor data from a similar 2021 spatial study can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1892052 and 2021 sample data can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/1898914. The 2021 spatial study had some sites in common with this 2022 spatial study. This dataset is comprised of three photo folders and one main data folder with six subfolders. The photo folders contain photographs and videos of cotton strip retrieval and sediment quadrats. The main data folder consists of (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) total suspended solids (TSS) data and cotton strip tensile strength data and averages; (5) field protocol; (6) readme; (7) methods codes; (8) international generic sample number (IGSN) mapping file; (9) sensor installation methods summary; (10) stream depth and averages; and (11) Ultrameter data and averages. The Sonar subfolder consists of Sonar time-series depth data and a processing script. The BarotrollAtm, DepthHOBO, MantaRiver, miniDOT, and miniDOTManualChamber subfolders contain time-series data, plots, and summary files. All files are .csv, .pdf, .txt, .R, .Rmd, .jpg, .jpeg, .AVI, .mp4, or .mov. The data package was originally published in April 2023. It was updated in August 2023 (v2; modified files) and September 2024 (v3; modified files). See the change history section in the readme for details. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

SAIL Aerosol Vertical Profiles (SAIL-AVP) Field Campaign Report

Aerosols are central to understanding the water cycle within mountainous regions, but a complete understanding of this role cannot be provided without vertically resolved observations, particularly for aerosol-cloud interactions, since we simultaneously need to know aerosol information and meteorology near cloud bases. The goal of SAIL Aerosol Vertical Profiles (SAIL-AVP) was to perform tethered balloon system (TBS) measurements during different seasons to elucidate process-level understanding of the aerosols and associated meteorological conditions within complex mountainous terrain. Measurements were conducted during the deployment of the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility’s second Mobile Facility (AMF2) as part of the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign in the East River Watershed of the Upper Colorado River Basin (UCRB) in southwestern Colorado. In mid-latitude continental interior mountain ranges, and especially at SAIL, aerosol vertical profiles and their relationships with meteorology are, almost without exception, inferred from a combination of ground-based and remote-sensing data. These flights provided many more dimensions to aerosol-cloud interactions that are unavailable to those inferential approaches: the TBS collected in situ information in the vertical, optical, chemical, and biological dimensions for aerosol-cloud interactions. Furthermore, the collocation of these flights with longstanding collaborative resources in the region, including the ongoing surface and subsurface hydrologic observations from the DOE’s Watershed Function Science Focus Area (SFA), produced a unique set of atmospheric observations that are complemented by existing land-surface and subsurface (e.g., groundwater) observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Temperature, Humidity, and Time-Lapse Video Data from the East River Watershed, Water Year 2024

A new version of this dataset is available at doi:10.15485/3001338 and is the first citation in the 'Related References' section. It is expands on this dataset by appending another water year of data collection and additional logger sites.This dataset contains time-lapse imagery and distributed measurements of air temperature, relative humidity, dew point, and soil temperature across the East River basin from 3 October 2023 to 12 August 2024. Instruments were deployed at 14 sites as part of the DOE Grant: Seasonal Cycles Unravel Mysteries of Missing Mountain Water organized by Jessica Lundquist (University of Washington), Rosemary Carroll (Desert Research Institute), and Ethan Gutmann (National Center for Atmospheric Research). The data are published to support studies of surface climate or hydrologic processes in complex terrain. Measurements were collected with low-cost data loggers installed 2 m high on evergreen trees or buried just below the soil surface. Time-lapse cameras were deployed at three sites. Imagery from sites AP BONUS and AP5 provides insight into large-scale seasonal snow cover variability. Imagery from site EL2 shows smaller-scale snow patterns across a nearby meadow.Dataset files are organized by site and variable (air measurements, ground measurements, or time-lapse video). Air and ground measurements are packaged in LoggerData.zip, and time-lapse imagery is compiled into short videos stored in TimelapseVideos.zip. File-level metadata contains details for each file included in the dataset. A data dictionary provides units and descriptions for column or row names in all files. The locations metadata file describes site characteristics, locations, and associated GPS methods.Dataset update 2025-03-03: Resolved header and datetime formatting inconsistencies within LoggerData.zip files KP1_Air, KP3_Air, KP3_Ground, KP6_Ground, AP3_Ground, AP4_Air, and AP6_Air.Dataset update 2025-11-18: Modified abstract and related references sections to include new version of dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temperature, Humidity, and Time-Lapse Video Data from the East River Watershed, Water Years 2024 and 2025

This dataset contains time-lapse imagery and distributed measurements of air temperature, relative humidity, dew point, and soil temperature across the East River basin from 3 October 2023 to 8 August 2025. Instruments were deployed at 19 sites as part of the DOE Grant: Seasonal Cycles Unravel Mysteries of Missing Mountain Water organized by Jessica Lundquist (University of Washington), Rosemary Carroll (Desert Research Institute), and Ethan Gutmann (National Center for Atmospheric Research). The data are published to support studies of surface climate or hydrologic processes in complex terrain. Measurements were collected with low-cost data loggers installed 2 m high on evergreen trees or buried just below the soil surface. Time-lapse cameras were deployed at three sites. Imagery from sites AP BONUS and AP5 (Avery Picnic) provides insight into large-scale seasonal snow cover variability. Imagery from site EL2 (Emerald Lake) shows smaller-scale snow patterns across a nearby meadow. Dataset files are organized by site and variable (air measurements, ground measurements, or time-lapse video). Air and ground measurements are packaged in LoggerData.zip, and time-lapse imagery is compiled into short videos stored in TimelapseVideos.zip. File-level metadata contains details for each file included in the dataset. A data dictionary provides units and descriptions for column or row names in all files. The locations metadata file describes site characteristics, locations, and associated GPS methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Feedbacks Between Surface Deformation and Permafrost Degradation in Ice Wedge Polygons, Arctic Coastal Plain, Alaska

In the past three decades, an abrupt, pan-Arctic acceleration of ice wedge melting has transformed tundra landscapes, spurring the formation of hummock-like features known as high-centered polygons (HCPs). This rapid geomorphic transition profoundly alters regional hydrology and influences surface emissions of CO 2 and CH 4 . In Arctic Alaska, most recent instances of ice wedge degradation have arrested within 15–20 years of inception, stabilizing HCP microtopography. However, feedbacks between ground surface deformation and permafrost stability are incompletely understood, limiting our capacity to predict trajectories of landscape evolution in a still warmer future. In this work, we use field data from a site near Prudhoe Bay, Alaska, to develop a modeling-based framework for assessing the strength of positive (i.e., exacerbating) feedbacks on ice wedge degradation, focusing on the importance of heterogeneity in surface drainage and microtopographic conditions. Our simulations suggest that, when troughs are narrow, positive feedbacks on ice wedge melting (associated with thermokarst pool formation) are relatively weak. Positive feedbacks are markedly stronger beneath wide troughs, such as those that form above older, larger ice wedges. Seasonal thaw abruptly accelerates once a talik begins to form beneath wide and deep thermokarst pools. Once a talik initiates, winter severity and snowpack thickness increase in importance as predictors of thaw intensity in summer. Our results indicate that meter-scale heterogeneity in polygonal microtopography potentially exerts strong, nonlinear controls on thermokarst trajectories. These findings are useful for predicting future thermokarst dynamics and for interpreting the results from coarser-resolution land surface models operating at greater spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site

Numerical simulations are essential tools for understanding the complex hydrologic response of Arctic regions to a warming climate. However, strong coupling among thermal and hydrological processes on the surface and in the subsurface and the significant role that subtle variations in surface topography have in regulating flow direction and surface storage lead to significant uncertainties. Careful model evaluation against field observations is thus important to build confidence.We evaluate the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations from polygonal tundra at the Barrow Environmental Observatory. ATS couples a multiphase, 3D representation of subsurface thermal hydrology with representations of overland nonisothermal flows, snow processes, and surface energy balance. We simulated thermal hydrology of a 3D ice-wedge polygon with geometry that is abstracted but broadly consistent with the surface microtopography at our study site. The simulations were forced by meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of parameters appearing in the soil evaporation model, the 3-year simulations agreed reasonably well with snow depth, summer water table elevations in the polygon center, and high-frequency soil temperature measurements at several depths in the trough, rim, and center of the polygon. Upscaled evaporation is in good agreement with flux tower observations. The simulations were found to be sensitive to parameters in the bare soil evaporation model, snowpack, and the lateral saturated hydraulic conductivity.Timing of fall freeze-up was found to be sensitive to initial snow density, illustrating the importance of including snow aging effects. The study provides new support for an emerging class of integrated surface/subsurface permafrost simulators.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗