Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subsurface 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On the transferability of residence time distributions in two 10-km long river sections with similar hydromorphic units

Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface and their residence time distributions (RTDs) in the subsurface are important for managing the water quality and ecosystem health in dynamic river corridors. However, direct simulating high-spatial resolution HEFs and RTDs can be time-consuming, especially for watershed-scale modeling. Efficient surrogate models linking RTDs to hydromorphic units (HUs) can be alternatives for simulating RTDs in large-scale models. A common concern of these surrogate models, though, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this issue, this work evaluates the HEFs and resulting RTD-HU relationships for two 10-km long river corridors along the Columbia River leveraging a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework we previously developed. Applying such a framework at the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. Finally, our comparison shows that the similarity and transferability of the RTD-HU relationship is very low for the two investigated river sections, which suggests that devising a general algorithm to estimate RTDs based solely on surface water hydrodynamics and short-distance river channel topography data, as well as HU classification, might be nearly impossible.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

A Proposal for an Integrated Geophysical Strategy to "Follow the Water" on Mars

The search for subsurface water has become a primary focus of Mars exploration. Its abundance and distribution (both as ground ice and groundwater) have important implications for understanding the geologic, hydrologic, and climatic evolution of the planet; the potential origin and continued survival of life; and the accessibility of a critical in situ resource for sustaining future human explorers. For these reasons, a principal goal of the Mars science, astrobiology, and the HEDS programs is to determine the 3-D distribution and state of subsurface H2O, at a resolution sufficient to permit reaching any desired volatile target by drilling. The three targets most often discussed are: groundwater, massive deposits of near-surface ground ice (associated with the ponded discharge of the outflow channels or the relic of a former ocean), and ice-saturated frozen ground. Based on the present best estimates of mean annual surface temperature, crustal thermal conductivity, geothermal heat flow, and groundwater freezing temperature, the mean thickness of frozen ground on Mars is expected to vary from approx. = 2.5 - 5 km at the equator to approx. = 6.5 - 13 km at the poles. However, natural variations in both crustal heat flow and thermal conductivity are likely to result in significant local departures from these predicted values. The recent discovery of "young" fluvial-like features, emanating from the slopes of local scarps, raises the possibility that liquid water may also exist episodically at shallow (approx. = 100 - 500 m) depth; however, the true nature and absolute age of these features remains highly uncertain. Although the belief that Mars is water-rich is supported by a wide variety geologic evidence, our ignorance about the heterogeneous nature and thermal evolution of the planet's crust effectively precludes geomorphic or theoretical attempts to quantitatively assess the current geographic and subsurface vertical distribution of ground ice and groundwater . For this reason, any exploration activity (such as drilling) whose success is contingent on the presence of subsurface water, must be preceded by a comprehensive high-resolution geophysical survey capable of assessing whether local reservoirs of water and ice actually exist. Terrestrial experience has demonstrated that the accurate identification of such targets is likely to require the application of multiple geophysical techniques. In this abstract we propose an integrated strategy for the geophysical exploration of Mars that we believe represents the fastest, most cost-effect, and technically capable approach to identifying the state and distribution of subsurface water. Additional information is contained in the original extended abstract.

Clifford, S. M.↗
NASA NTRS2019

Climates of Warm Earth-Like Planets III: Fractional Habitability from a Water Cycle Perspective

The habitable fraction of a planet's surface is important for the detectability of surface biosignatures. The extent and distribution of habitable areas are influenced by external parameters that control the planet's climate, atmospheric circulation, and hydrological cycle. We explore these issues using the ROCKE-3D general circulation model, focusing on terrestrial water fluxes and thus the potential for the existence of complex life on land. Habitability is examined as a function of insolation and planet rotation for an Earth-like world with zero obliquity and eccentricity orbiting the Sun. We assess fractional habitability using an aridity index that measures the net supply of water to the land. Earth-like planets become "superhabitable" (a larger habitable surface area than Earth) as insolation and day-length increase because their climates become more equable, reminiscent of past warm periods on Earth when complex life was abundant and widespread. The most slowly rotating, most highly irradiated planets, though, occupy a hydrological regime unlike any on Earth, with extremely warm, humid conditions at high latitudes but little rain and subsurface water storage. Clouds increasingly obscure the surface as insolation increases, but visibility improves for modest increases in rotation period. Thus, moderately slowly rotating rocky planets with insolation near or somewhat greater than modern Earth's appear to be promising targets for surface characterization by a future direct imaging mission.

Del Genio, Anthony D.↗
NASA NTRS2019

Climates of Warm Earth-like Planets. III. Fractional Habitability from a Water Cycle Perspective

The habitable fraction of a planet's surface is important for the detectability of surface biosignatures. The extent and distribution of habitable areas are influenced by external parameters that control the planet's climate, atmospheric circulation, and hydrological cycle. We explore these issues using the ROCKE-3D general circulation model, focusing on terrestrial water fluxes and thus the potential for the existence of complex life on land. Habitability is examined as a function of insolation and planet rotation for an Earth-like world with zero obliquity and eccentricity orbiting the Sun. We assess fractional habitability using an aridity index that measures the net supply of water to the land. Earth-like planets become "superhabitable" (a larger habitable surface area than Earth) as insolation and day-length increase because their climates become more equable, reminiscent of past warm periods on Earth when complex life was abundant and widespread. The most slowly rotating, most highly irradiated planets, though, occupy a hydrological regime unlike any on Earth, with extremely warm, humid conditions at high latitudes but little rain and subsurface water storage. Clouds increasingly obscure the surface as insolation increases, but visibility improves for modest increases in rotation period. Thus, moderately slowly rotating rocky planets with insolation near or somewhat greater than modern Earth's appear to be promising targets for surface characterization by a future direct imaging mission.

Exoplanet atmospheric variability↗
DOE OSTI2025

Controls From Above and Below: Snow, Soil, and Steepness Drive Diverging Trends of Subsurface Water and Streamflow Dynamics

ABSTRACT The importance of subsurface water dynamics, such as water storage and flow partitioning, is well recognised. Yet, our understanding of their drivers and links to streamflow generation has remained elusive, especially in small headwater streams that are often data‐limited but crucial for downstream water quantity and quality. Large‐scale analyses have focused on streamflow characteristics across rivers with varying drainage areas, often overlooking the subsurface water dynamics that shape streamflow behaviour. Here we ask the question: What are the climate and landscape characteristics that regulate subsurface dynamic storage, flow path partitioning, and dynamics of streamflow generation in headwater streams? To answer this question, we used streamflow data and a widely‐used hydrological model (HBV) for 15 headwater catchments across the contiguous United States. Results show that climate characteristics such as aridity and precipitation phase (snow or rain) and land attributes such as topography and soil texture are key drivers of streamflow generation dynamics. In particular, steeper slopes generally promoted more streamflow, regardless of aridity. Streams in flat, rainy sites (< 30% precipitation as snow) with finer soils exhibited flashier regimes than those in snowy sites (> 30% precipitation as snow) or sites with coarse soils and deeper flow paths. In snowy sites, less weathered, thinner soils promoted shallower flow paths such that discharge was more sensitive to changes in storage, but snow dampened streamflow flashiness overall. Results here indicate that land characteristics such as steepness and soil texture modify subsurface water storage and shallow and deep flow partitioning, ultimately regulating streamflow response to climate forcing. As climate change increases uncertainty in water availability, understanding the interacting climate and landscape features that regulate streamflow will be essential to predict hydrological shifts in headwater catchments and improve water resources management.

Kerins, Devon [Department of Civil and Environment↗
DOE OSTI2025

Evaluating post-fire watershed response to varying burn severity and precipitation regimes using fully-distributed and integrated hydrologic models

Wildfires can cause significant changes in vegetation and soil, which may lead to increased surface runoff and soil erosion, thereby affecting water cycling within ecosystems. This study uses the Advanced Terrestrial Simulator (ATS), an integrated and fully distributed hydrologic model at the watershed scale, to examine post-fire hydrologic responses in selected watersheds with varying burn severities in the Pacific Northwest region of the United States. The model integrates surface overland flow, subsurface flow, and canopy biophysical processes. We have developed a new fire module in ATS to account for changes in soil hydraulic properties caused by fire in the topsoil layer. Modeling results indicate that, in the year following a high-severity burn, watershed-averaged evapotranspiration decreases by 25%. Additionally, post-fire peak flows increase by 18-29% in watersheds burned with medium to high severity due to changes in soil properties. Conversely, a low-severity burn results in less than a 1% increase in post-fire peak flow. Furthermore, a high-severity fire causes a 38% reduction in the infiltration rate within the affected watershed during the first post-fire wet season. Hypothetical numerical experiments with varying precipitation regimes after a high-severity fire show that post-fire peak flows can increase by 1-29% due to fire-induced changes in soil hydraulic properties. This study highlights the importance of using fully distributed hydrologic models to quantify disturbance-feedback loops, which are essential for understanding the complexities brought about by spatial heterogeneity in post-fire landscapes.

ATS↗
NASA NTRSDate not supplied

Phyllosilicate Formation on Early Mars Via Open-System Acid Alteration of Basaltic Glass

Smectites are widespread on Mars, but neutral/alkaline conditions favorable for smectite formation have surprisingly not led to abundant carbonates. Smectite formation on Mars could also occur in acidic environments unfavorable for carbonate formation. Acidic smectite formation has been demonstrated in batch (closed hydrologic) systems, however, the mechanisms and octahedral composition of smectite forming in acidic flow-through (open hydrologic) systems are not fully understood. Stapafell basaltic glass was hydrothermally altered (190◦C) at two flow rates corresponding to low and high water to rock ratio (W/R) and initial pH (pH_0) values of 2, 3, 4 and 6, and a batch low W/R experiment was conducted at pH_0 2. Kaolinite, montmorillonite and chlorite formed at pH_0 2 at low W/R; no phyllosilicates formed at pH_0 2 at high W/R; and lizardite formed at pH_0 ≥3 at both W/R ratios. Lizardite, kaolinite, and montmorillonite in these experiments formed by precipitation from solution and chlorite likely formed through alteration of montmorillonite and/or basalt. Saponite formed at pH_0 2 in batch conditions by alteration of basaltic glass. Comparison of experimental data with martian phyllosilicate assemblages indicated that smectite formation on Mars likely occurred under water-limited conditions. Al-rich smectite could form in open-system low W/R subsurface environments under a narrow range of pH (pH <3) while saponite could form in closed low W/R systems under acidic to alkaline conditions. Combined open and closed hydrological regimes could be responsible for development of clay mineral stratigraphies observed on Mars. The acidic conditions required for formation of Al-rich smectite were unfavorable for carbonate precipitation, but carbonate precipitation could occur together with Fe/Mg-smectite in closed systems at pH >4. The lack of carbonates occurring together with Fe/Mg-smectite could be cause by low pCO2 in subsurface closed environments and/or by overall lack of carbonate precipitation despite significant oversaturation.

S J Ralston↗
NASA NTRS2000

Detection and Characterization of Martian Volatile-Rich Reservoirs: The Netlander Approach

Geological and theoretical modeling do indicate that, most probably, a significant part of the volatiles present in the past is presently stocked within the Martian subsurface as ground ice, and as clay minerals (water constitution). The detection of liquid water is of prime interest and should have deep implications in the understanding of the Martian hydrological cycle and also in exobiology. In the frame of the 2005 joint CNES-NASA mission to Mars, a set of 4 NETLANDERs developed by an European consortium is expected to be launched between 2005 and 2007. The geophysical package of each lander will include a geo-radar (GPR experiment), a magnetometer (MAGNET experiment), a seismometer (SEIS experiment) and a meteorological package (ATMIS experiment). The NETLANDER mission offers a unique opportunity to explore simultaneously the subsurface as well as deeper layers of the planetary interior on 4 different landing sites. The complementary contributions of all these geophysical soundings onboard the NETLANDER stations are presented.

Banerdt, B.↗
NASA NTRS2007

Salinity Remote Sensing and the Study of the Global Water Cycle

The SMOS and AquariusISAC-D satellite missions will begin a new era to map the global sea surface salinity (SSS) field and its variability from space within the next twothree years. They will provide critical data needed to study the interactions between the ocean circulation, global water cycle and climate. Key scientific issues to address are (1) mapping large expanses of the ocean where conventional SSS data do not yet exist, (2) understanding the seasonal and interannual SSS variations and the link to precipitation, evaporation and sea-ice patterns, (3) links between SSS and variations in the oceanic overturning circulation, (4) air-sea coupling processes in the tropics that influence El Nino, and (4) closing the marine freshwater budget. There is a growing body of oceanographic evidence in the form of salinity trends that portend significant changes in the hydrologic cycle. Over the past several decades, highlatitude oceans have become fresher while the subtropical oceans have become saltier. This change is slowly spreading into the subsurface ocean layers and may be affecting the strength of the ocean's therrnohaline overturning circulation. Salinity is directly linked to the ocean dynamics through the density distribution, and provides an important signature of the global water cycle. The distribution and variation of oceanic salinity is therefore attracting increasing scientific attention due to the relationship to the global water cycle and its influence on circulation, mixing, and climate processes. The oceans dominate the water cycle by providing 86% of global surface evaporation (E) and receiving 78% of global precipitation (P). Regional differences in E-P, land runoff, and the melting or freezing of ice affect the salinity of surface water. Direct observations of E-P over the ocean have large uncertainty, with discrepancies between the various state-of-the-art precipitation analyses of a factor of two or more in many regions. Quantifying the climatic influence of the oceanic water cycle requires more accurately resolving the net air-sea water flux. Measuring global SSS trends on seasonal to interannual timescales by satellite is fundamental to this problem because the SSS trends represent detectable time-integrated signals of the variable marine hydrological cycle. Satellite measurements, coupled with an array of in situ observations, will provide global synoptic SSS fields for the first time history. These data will provide a strong constraint on climate models and data assimilation efforts, which must properly represent the freshwater budget in terms of E-P, ocean advection and surface layer mixing in order to accurately simulate the true ocean state. The SSS fields will allow us to quantify the covariability between the SSS and the strong seasonal E-P cycle in the tropics and high latitudes. Field measurement campaigns to exploit satellite and in situ measurements to close the seasonal E-P cycle over an ocean region are being considered. Lastly the satellite systems will monitor and trace the large long-lived SSS anomalies from year to year that have the potential to influence El Nino and the large scale ocean circulation.

Lagerloef, G. S. E.↗
NASA NTRS1996

Passive Microwave Remote Sensing of Soil Moisture

Microwave remote sensing provides a unique capability for direct observation of soil moisture. Remote measurements from space afford the possibility of obtaining frequent, global sampling of soil moisture over a large fraction of the Earth's land surface. Microwave measurements have the benefit of being largely unaffected by cloud cover and variable surface solar illumination, but accurate soil moisture estimates are limited to regions that have either bare soil or low to moderate amounts of vegetation cover. A particular advantage of passive microwave sensors is that in the absence of significant vegetation cover soil moisture is the dominant effect on the received signal. The spatial resolutions of passive Microwave soil moisture sensors currently considered for space operation are in the range 10-20 km. The most useful frequency range for soil moisture sensing is 1-5 GHz. System design considerations include optimum choice of frequencies, polarizations, and scanning configurations, based on trade-offs between requirements for high vegetation penetration capability, freedom from electromagnetic interference, manageable antenna size and complexity, and the requirement that a sufficient number of information channels be available to correct for perturbing geophysical effects. This paper outlines the basic principles of the passive microwave technique for soil moisture sensing, and reviews briefly the status of current retrieval methods. Particularly promising are methods for optimally assimilating passive microwave data into hydrologic models. Further studies are needed to investigate the effects on microwave observations of within-footprint spatial heterogeneity of vegetation cover and subsurface soil characteristics, and to assess the limitations imposed by heterogeneity on the retrievability of large-scale soil moisture information from remote observations.

Njoku, Eni G.↗
NASA NTRS1995

The circum-Chryse region as a possible example of a hydrologic cycle on Mars: Geologic observations and theoretical evaluation

The transection and superposition relationships among channels, chaos, surface materials units, and other features in the circum-Chryse region of Mars were used to evaluate relative age relationships and evolution of flood events. Channels and chaos in contact (with one another) were treated as single discrete flood-carved systems. Some outflow channel systems form networks and are inferred to have been created by multiple flood events. Within some outflow channel networks, several separate individual channel systems can be traced to a specific chaos which acted as flood-source area to that specific flood channel. Individual flood-carved systems were related to widespread materials units or other surface features that served as stratigraphic horizons. Chryse outflow channels are inferred to have formed over most of the perceivable history of Mars. Outflow channels are inferred to become younger with increasing proximity to the Chryse basin. In addition, outflow channels closer to the basin show a greater diversity in age. The relationship of subsequent outflow channel sources to the sources of earlier floods is inferred to disfavor episodic flooding due to the progressive tapping of a juvenile near-surface water supply. Instead, we propose the circum-Chryse region as a candidate site of past hydrological recycling. The discharge rates necessary to carve the circum-Chryse outflow channels would have inevitably formed temporary standing bodies of H2O on the Martian surface where the flood-waters stagnated and pooled (the Chryse basin is topographically enclosed). These observations and inferences have led us to formulate and evaluate two hypotheses: (1) large amounts of the sublimated H2O off the Chryse basin flood lakes precipitated (snowed) onto the flood-source highlands and this H2O was incorporated into the near surface, recharging the H2O sources, making possible subsequent deluges; and (2) ponded flood-water in Chryse basin drained back down an anti basinward dipping subsurface layer accessed along the southern edge of the lake, recharging the flood-source aquifers. H2O not redeposited in the flood-source region was largely lost to the hydrologic cycle. This loss progressively lowered the vitality of the cycle, probably by now killing it. Our numerical evaluations indicate that of the two hypotheses we formulated, the groundwater seep cycle seems by far the more viable. Optimally, approx. 3/4 of the original mass of an ice-covered cylindrical lake (albedo 0.5, 1 km deep, 100-km radius, draining along its rim for one quarter of its circumference into substrata with a permeability of 3000 darcies) can be modeled to have moved underground (on timescales of the order of 10(exp 3) years) before the competing mechanisms of sublimation and freeze down choked off further water removal. Once underground, this water can travel distances equal to the separation between Chryse basin and flood-source sites in geologically short (approx. 10(exp 6) year-scale) times. Conversely, we calculate that optimally only approx. 40% of the H2O carried from Chryse can condense at the highlands, and most of the precipitate would either collect at the base of the highlands/lowlands scarp or sublimate at rates greater than it would accumulate over the flood-source sites. Further observations from forthcoming missions may permit the determination of which mechanisms may have operated to recycle the Chryse flood-waters.

Moore, Jeffrey M.↗
NASA NTRS1995

The Circum-Chryse Region as a Possible Example of a Hydrologic Cycle on Mars: Geologic Observations and Theoretical Evaluation

The transection and superposition relationships among channels, chaos, surface materials units, and other features in the circum-Chryse region of Mars were used to evaluate relative age relationships and evolution of flood events. Channels and chaos in contact (with one another) were treated as single discrete flood-carved systems. Some outflow channel systems form networks and are inferred to have been created by multiple flood events. Within some outflow channel networks, several separate individual channel systems can be traced to a specific chaos which acted as flood-source area to that specific flood channel. Individual flood-carved systems were related to widespread materials units or other surface features that served as stratigraphic horizons. Chryse outflow channels are inferred to have formed over most of the perceivable history of Mars. Outflow channels are inferred to become younger with increasing proximity to the Chryse basin. In addition, outflow channels closer to the basin show a greater diversity in age. The relationship of subsequent outflow channel sources to the sources of earlier floods is inferred to disfavor episodic flooding due to the progressive tapping of a juvenile near-surface water supply. Instead, we propose the circum-Chryse region as a candidate site of past hydrological recycling. The discharge rates necessary to carve the circum-Chryse outflow channels would have inevitably formed temporary standing bodies of H2O on the Martian surface where the flood-waters stagnated and pooled (the Chryse basin is topographically enclosed). These observations and inferences have led us to formulate and evaluate two hypotheses: (1) large amounts of the sublimated H2O off the Chryse basin flood lakes precipitated (snowed) onto the flood-source highlands and this H2O was incorporated into the near surface, recharging the H2O sources, making possible subsequent deluges; and (2) ponded flood-water in Chryse basin drained back down an anti basinward dipping subsurface layer accessed along the southern edge of the lake, recharging the flood-source aquifers. H2O not redeposited in the flood-source region was largely lost to the hydrologic cycle. This loss progressively lowered the vitality of the cycle, probably by now killing it. Our numerical evaluations indicate that of the two hypotheses we formulated, the groundwater seep cycle seems by far the more viable. Optimally, approximately 3/4 of the original mass of an ice-covered cylindrical lake (albedo 0.5, 1 km deep, 100-km radius, draining along its rim for one quarter of its circumference into substrata with a permeability of 3000 darcies) can be modeled to have moved underground (on timescales of the order of 10(exp 3) years) before the competing mechanisms of sublimation and freeze down choked off further water removal. Once underground, this water can travel distances equal to the separation between Chryse basin and flood-source sites in geologically short (approximately 10(exp 6) year-scale) times. Conversely, we calculate that optimally only approximately 40% of the H2O carried from Chryse can condense at the highlands, and most of the precipitate would either collect at the base of the highlands/lowlands scrap or sublimate at rates greater than it would accumulate over the flood-source sites. Further observations from forthcoming missions may permit the determination of which mechanisms may have operated to recycle the Chryse flood-waters.

Moore, Jeffrey M.↗
NASA NTRS2000

Short-Term Retrospective Land Data Assimilation Schemes

Subsurface moisture and temperature and snow/ice stores exhibit persistence on various time scales that has important implications for the extended prediction of climatic and hydrologic extremes. Hence, to improve their specification of the land surface, many numerical weather prediction (NWP) centers have incorporated complex land surface schemes in their forecast models. However, because land storages are integrated states, errors in NWP forcing accumulates in these stores, which leads to incorrect surface water and energy partitioning. This has motivated the development of Land Data Assimilation Schemes (LDAS) that can be used to constrain NWP surface storages. An LDAS is an uncoupled land surface scheme that is forced primarily by observations, and is therefore less affected by NWP forcing biases. The implementation of an LDAS also provides the opportunity to correct the model's trajectory using remotely-sensed observations of soil temperature, soil moisture, and snow using data assimilation methods. The inclusion of data assimilation in LDAS will greatly increase its predictive capacity, as well as provide high-quality land surface assimilated data.

Houser, P. R.↗
NASA NTRS2004

Eastern Sahara Geology from Orbital Radar: Potential Analog to Mars

Much of the surface of Mars has been intensely reworked by aeolian processes and key evidence about the history of the Martian environment seems to be hidden beneath a widespread layer of debris (paleo lakes and rivers, faults, impact craters). In the same way, the recent geological and hydrological history of the eastern Sahara is still mainly hidden under large regions of wind-blown sand which represent a possible terrestrial analog to Mars. The subsurface geology there is generally invisible to optical remote sensing techniques, but radar images obtained from the Shuttle Imaging Radar (SIR) missions were able to penetrate the superficial sand layer to reveal parts of paleohydrological networks in southern Egypt.

Farr, T. G.↗
NASA NTRS2005

Surface Clutter Removal in Airborne Radar Sounding Data from the Dry Valleys, Antarctica

We have collected roughly 1,000 line-km of airborne radar sounding data over glaciers, rock/ice glaciers, permafrost, subsurface ice bodies, ice-covered saline lakes, and glacial deposits in Taylor and Beacon Valley. These data are being analyzed in order to develop techniques for discriminating between subsurface and off-nadir echoes and for detecting and characterizing subsurface interfaces. The identification of features on Mars exhibiting morphologies consistent with ice/rock mixtures, near-surface ice bodies and near-surface liquid water, and the importance of such features to the search for water on Mars, highlights the need for appropriate terrestrial analogs and analysis techniques in order to prepare for radar sounder missions to Mars. Climatic, hydrological, and geological conditions in the Dry Valleys of Antarctica are analogous in many ways to those on Mars. A crucial first step in the data analysis process is the discrimination of echo sources in the radar data. The goal is to identify all returns from the surface of off-nadir topography in order to positively identify subsurface echoes. This process will also be critical for radar data that will be collected in areas of Mars exhibiting significant topography, so that subsurface echoes are identified unambiguously. The positive detection and characterization of subsurface (including sub-ice) water is a primary goal of NASA's Mars exploration program. Our data over the Dry Valleys provides an opportunity to implement techniques we are developing to accomplish these goals.

Holt, J. W.↗
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2024

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗