Search NASA⌕ Search

SEARCH · Search NASA

Results for “Catchment scale”

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

Machine Learning Classification Strategy to Improve Streamflow Estimates in Diverse River Basins in the Colorado River Basin

Streamflow in the Colorado River Basin (CRB) is significantly altered by human activities including land use/cover alterations, reservoir operation, irrigation, and water exports. Climate is also highly varied across the CRB which contains snowpack-dominated watersheds and arid, precipitation-dominated basins. Recently, machine learning methods have improved the generalizability and accuracy of streamflow models. Previous successes with LSTM modeling have primarily focused on unimpacted basins, and few studies have included human impacted systems in either regional or single-basin modeling. We demonstrate that the diverse hydrological behavior of river basins in the CRB are too difficult to model with a single, regional model. We propose a method to delineate catchments into categories based on the level of predictability, hydrological characteristics, and the level of human influence. Lastly, we model streamflow in each category with climate and anthropogenic proxy data sets and use feature importance methods to assess whether model performance improves with additional relevant data. Overall, land use cover data at a low temporal resolution was not sufficient to capture the irregular patterns of reservoir releases, demonstrating the importance of having high-resolution reservoir release data sets at a global scale. On the other hand, the classification approach reduced the complexity of the data and has the potential to improve streamflow forecasts in human-altered regions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Tectonic-Climate Interactions in Action Orogenic Belts: Quantification of Dynamic Topography with SRTM data

This project was undertaken to examine the approach to steady state in collisional mountain belts. Although the primary thrust of this grant was to look at larger collisional mountain belts, such as the Himalaya, the Tien Shan, and Southern Alps, we began by looking at smaller structures represented by growing and propagating folds. Like ranges that are evolving toward a topographic steady state, these folds undergo a series of morphologic changes as they are progressively uplifted and eroded. We wanted to document the nature of these changes and to try to discern some of the underlying controls on them. We initially focused on the Wheeler Ridge anticline in southern California. Subsequently, we progressed to looking at the topographic development and the effects of differential uplift and glaciation on the Kyrgyz Range in the northern Tien Shan. This range is unusual inasmuch as it is transformed along its length from a simple uplift with a largely preserved Mesozoic erosion surface arching across it to a highly dissected and heavily glaciated uplift in the region where uplift has been sustained at higher rates over longer intervals. In efforts to understand the distribution of erosion rates at 10(exp 3) - 10(exp 5) year time scales, cosmogenic radionuclide (CRN) concentrations have been gaining increasingly widespread usage (Brown et al., 1995; Riebe et al., 2004; Riebe et al., 2001; Vance et al., 2003). Most studies to date, however, have been conducted in slowly eroding ranges. In rapidly eroding mountains where landslides deliver most of the sediments to the rivers, we hypothesized that CRN concentrations could be highly perturbed by the stochastic processes of landsliding. Therefore, we undertook the development of a numerical model that simulated the effects of both landsliding and grain-by-grain attrition within fluvial catchments. This modeling effort has shown the effects of catchment size and erosion rate on CRN concentrations and allows a prediction of where to sample to obtain the optimal erosion rate estimates using CRN techniques. Finally, we developed computational techniques to operate on DEMs to extract useful information that would enable quantification of climate-erosion interactions. In particular, we worked on rapid techniques to define catchments of any given range of sizes, to extract channel gradients, to combine precipitation information to calculate discharge, and to utilize various stream-power models to determine the erosional energy within any given catchment within a transect. We briefly describe results from Wheeler Ridge, the Kyrgyz Range, the Nepal Himalaya, and our numerical modeling.

Burbank, Douglas W.↗
NASA NTRS2012

Inferring Land Surface Model Parameters for the Assimilation of Satellite-Based L-Band Brightness Temperature Observations into a Soil Moisture Analysis System

The Soil Moisture and Ocean Salinity (SMOS) satellite mission provides global measurements of L-band brightness temperatures at horizontal and vertical polarization and a variety of incidence angles that are sensitive to moisture and temperature conditions in the top few centimeters of the soil. These L-band observations can therefore be assimilated into a land surface model to obtain surface and root zone soil moisture estimates. As part of the observation operator, such an assimilation system requires a radiative transfer model (RTM) that converts geophysical fields (including soil moisture and soil temperature) into modeled L-band brightness temperatures. At the global scale, the RTM parameters and the climatological soil moisture conditions are still poorly known. Using look-up tables from the literature to estimate the RTM parameters usually results in modeled L-band brightness temperatures that are strongly biased against the SMOS observations, with biases varying regionally and seasonally. Such biases must be addressed within the land data assimilation system. In this presentation, the estimation of the RTM parameters is discussed for the NASA GEOS-5 land data assimilation system, which is based on the ensemble Kalman filter (EnKF) and the Catchment land surface model. In the GEOS-5 land data assimilation system, soil moisture and brightness temperature biases are addressed in three stages. First, the global soil properties and soil hydraulic parameters that are used in the Catchment model were revised to minimize the bias in the modeled soil moisture, as verified against available in situ soil moisture measurements. Second, key parameters of the "tau-omega" RTM were calibrated prior to data assimilation using an objective function that minimizes the climatological differences between the modeled L-band brightness temperatures and the corresponding SMOS observations. Calibrated parameters include soil roughness parameters, vegetation structure parameters, and the single scattering albedo. After this climatological calibration, the modeling system can provide L-band brightness temperatures with a global mean absolute bias of less than 10K against SMOS observations, across multiple incidence angles and for horizontal and vertical polarization. Third, seasonal and regional variations in the residual biases are addressed by estimating the vegetation optical depth through state augmentation during the assimilation of the L-band brightness temperatures. This strategy, tested here with SMOS data, is part of the baseline approach for the Level 4 Surface and Root Zone Soil Moisture data product from the planned Soil Moisture Active Passive (SMAP) satellite mission.

Reichle, Rolf H.↗
NASA NTRSDate not supplied

Systematic Errors in Simulated L-Band Brightness Temperature in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF), the SMAP Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications. The EnKF-based L4_SM analysis assumes unbiased forecast errors. Consequently, the seasonally varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations. The L4_SM system is thus designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (possibly wrong) climatology. In this paper, we examine the Tb observation-minus-forecast (O-F) residuals from the L4_SM Version 7 product (computed after rescaling the Tb observations to the mean seasonal cycle of the simulated Tb). The long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. The model forecast Tb, however, still exhibits undesirable systematic errors relative to the (rescaled) Tb observations. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to improve the L4_SM algorithm.

Rolf Reichle↗
DOE OSTI2025

Ensemble‐Based Spatially Distributed CLM5 Hydrological Parameter Estimation for the Continental United States

Abstract One of the major challenges in large‐domain hydrological modeling efforts lies in the estimation of spatially distributed hydrological parameters while simultaneously accounting for their associated uncertainties. Addressing this challenge is particularly difficult in ungauged locations. With growing societal demands for large‐scale streamflow projections to inform water resource management and long‐term planning, evaluating and constraining hydrological parameter uncertainty is increasingly vital. This study introduces a hybrid regionalization approach to enhance hydrological predictions of the Community Land Model version 5 (CLM5) across the Continental United States (CONUS), with a total of 50,629 1/8° grid cells. This hybrid method combines the strengths of two existing techniques: parameter regionalization and streamflow signature regionalization. It identifies ensemble behavioral parameters for each 1/8° grid cell across the CONUS domain, tailored to three distinct streamflow signatures focused on low flows, high flows, and annual water balance. Evaluating this hybrid method for 464 CAMELS (Catchment Attributes and Meteorology for Large‐sample Studies) basins demonstrates a significant improvement in CLM5 hydrological predictions, even in challenging arid regions. In CONUS applications, the derived spatially distributed parameter sets capture both spatial continuity and variation of parameters, highlighting their heterogeneous nature within specific regions. Overall, this hybrid regionalization approach offers a promising solution to the complex task of improving hydrological modeling over large domains for important hydrological applications.

CLM5↗
NASA NTRS2023

Regional Analysis of the 2015–16 Lower Mekong River Basin Drought Using NASA Satellite Observations

Study region Lower Mekong River Basin (LMRB) Study focus Satellite remote sensing products are widely used for monitoring droughts. Using NASA satellite sensors of precipitation (Global Measurement Mission, GPM), soil moisture (Soil Moisture Active and Passive, SMAP), and terrestrial water storage (Gravity Recovery and Climate Experiment, GRACE), this study evaluates the historical drought in the LMRB during 2015–16. SMAP soil moisture was validated against in-situ soil moisture, and GPM precipitation and SMAP soil moisture were cross-validated with streamflow observations. The spatiotemporal dynamics of soil moisture were also examined in different ranges of catchment areas. In performing the analysis, we used lagged correlations between hydrological variables and the indices of the Standardized Precipitation Index (SPI) and Standardized Streamflow Index (SSI). New hydrological insights for the regions Spatio-temporal patterns of drought in 2015–16 were examined from the entire basin to small watersheds. A mismatch occurs when using GRACE data to study droughts in small watersheds (many of the small watersheds would be a fraction of the few 100 km2 spatial resolutions of GRACE pixel). In smaller watersheds, hydrological drought (SSI) was closely defined with SMAP soil moisture downscaled to 1 km rather than the meteorological drought index (SPI). By leveraging satellite-based observations across a range of spatial scales, this study highlights the utility of Earth observations in informing water resources and land management decisions at the regional scale.

Mekong↗
NASA NTRS2014

A Spatially Coherent Global Soil Moisture Product with Improved Temporal Resolution

Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.

spatial↗
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↗
NASA NTRS2023

Investigating Hydroclimatic Impacts of the 168–158 BCE Volcanic Quartet and Their Relevance to the Nile River Basin and Egyptian History

The Ptolemaic era (305–30 BCE) is an important period of Ancient Egyptian history known for its material and scientific advances, but also intermittent political and social unrest in the form of (sometimes widespread) revolts against the Ptolemaic elites. While the role of environmental pressures has long been overlooked in this period of Egyptian history, ice-core-based volcanic histories have identified the period as experiencing multiple notable eruptions, and a repeated temporal association between explosive volcanism and revolt has recently been noted. Here we analyze the global and regional (Nile River basin) hydroclimatic response to a unique historical sequence of four large and closely timed volcanic eruptions (first a tropical one, followed by three extratropical northern hemispheric events) between 168 and 158 BCE, a particularly troubled period in Ptolemaic history for which we now provide a more detailed hydroclimatic context. The NASA (National Aeronautics and Space Administration) GISS (Goddard Institute for Space Studies) ModelE2.1 Earth system model simulates a strong radiative response with a radiative forcing (top of atmosphere) of −7.5 W m−2 (following the first eruption) and −2.5 W m−2 (after each of the three remaining eruptions) at a global scale. Associated with this, we observe a global surface cooling of the order of 1.5 ∘C following the first (tropical) eruption, with the following three extratropical eruptions extending the cooling period for more than 15 years. Consequently, this series of eruptions is observed to constrain the northward migration of the inter-tropical convergence zone (ITCZ) during the Northern Hemisphere summer monsoon season, and major monsoon zones (African, South Asian, and East Asian) were seen to experience a suppression of rainfall of >1 mm d−1 during the monsoon (JJAS) season averaged for 2 years after each eruption. A substantial suppression of the Indian and North African summer monsoon (over the Nile River headwater region) was seen to strongly affect the modeled river flow in the catchment and discharge at river mouth. River mass flow over the basin was observed to decrease by 29 % and 38 % relative to an unperturbed (non-volcanic) annual mean flow in the first and second year, respectively, after the first (i.e., tropical) eruption. A moderate decrease ranging between 5 % and 18 % was observed after the third and fourth (extratropical) eruptions. These results indicate, in sum, that the first eruption likely produced a strong hydroclimate response, with the following extratropical eruptions prolonging this. These results also support the recently hypothesized association between ice-core-based signals of explosive volcanism and hydroclimatic variability during the Ptolemaic era, including the suppression of the agriculturally critical Nile summer flooding.

hydroclimate impacts↗
NASA NTRS2019

Application of GRACE to the Estimation of Groundwater Storage Change in a Data-poor Region: A Case Study of Ngadda Catchment in the Lake Chad Basin

The present study is to explore the feasibility of GRACE-based estimation of a groundwater storage change in a data-poor region using a case study of the Ngadda catchment in the Lake Chad Basin. Although the Ngadda catchment has only one set of in situ time series data of groundwater from 2006 to 2009 and a limited number of groundwater measurements in 2005 and 2009, GRACE-based groundwater storage change can be evaluated against thein situ groundwater measurements combined with specific yield data. The cross-correlation analysis in the Ngadda catchment shows that maximum rainfall reached in July and August, whereas both the maximum total water storage anomaly and the maximum groundwater storage anomaly occurred 2months later. Whereas the mean annual amplitude of total water storage anomaly is about 17cm from both the average total water storage anomaly from three mascon products and the one from three spherical harmonic products, the mean annual amplitude of soil moisture storage anomaly is substantially varied from5.58cm for CLM to about 14cm for NOAH and Mosaic. The goodness-of-fit tests show that CLM soil moisture produces the closest estimation of groundwater storage anomaly to thein situ groundwater measurements. The present study shows that GRACE-based estimation for groundwater storage anomaly can be a cost-effective and alternative tool to observe how groundwater changes in a basin scale under the limitation of modelling and in situ data availability.

GRACE↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

Influence of Global Climate on Freshwater Changes in Africa’s Largest Endorheic Basin Using Multi-Scaled Indicators

The poor investments in gauge measurements for hydro-climatic research in Africa has necessitated the need to investigate how decision makers can leverage on sophisticated spaceborne measurements to improve knowledge on surface water hydrology that can feed directly into water accounting processes and risk assessment from extreme droughts and its impacts. To demonstrate such potential, a suite of satellite earth observations (Sentinel-2, altimetry, Landsat, GRACE, and TRMM) and model data are combined with the standardized precipitation evapotranspiration index to assess the impacts of global climate on freshwater dynamics over the LCB (Lake Chad basin), Africa’s largest endorheic basin. As shown in the results of this study, the significant relationship of climate modes (AMO; r = 0.68 and 0.59; and AMM; r = 0.2 and 0.47) with drought patterns in the LCB highlights the evidence of global climate influence in the region. The significant declines in drought extents and their intensities (2004 - 2015) over LCB coincide with the rise in surface water extent of the Lake Chad during the same period. Change detection analysis of open water features in the southern pool of Lake Chad during the 2015 - 2019 period shows that on the average, only 28.4% of inundated areas within the vicinity of the Lake persisted during the period. While the association of terrestrial water storage (TWS) with model-derived surface water storage (SWS) is strongest (r = 0.89) in the catchments that provide the most nourishment to the Lake Chad, the relationship of rainfall (2002 - 2017) with TWS (r = 0.85), model TWS (r = 0.87) and SWS (r = 0.88) confirm that the LCB’s hydrology is predominantly climate-driven. This notion is further reinforced as the predicted SWS over the LCB using a support vector machine regression scheme was found to be strongly correlated (r = 0.95 at = 0.05) with observed SWS.

Sentinel-2↗
NASA NTRS2010

Mud Volcanoes - Analogs to Martian Cones and Domes (by the Thousands!)

Mud volcanoes are mounds formed by low temperature slurries of gas, liquid, sediments and rock that erupt to the surface from depths of meters to kilometers. They are common on Earth, with estimates of thousands onshore and tens of thousands offshore. Mud volcanoes occur in basins with rapidly-deposited accumulations of fine-grained sediments. Such settings are ideal for concentration and preservation of organic materials, and mud volcanoes typically occur in sedimentary basins that are rich in organic biosignatures. Domes and cones, cited as possible mud volcanoes by previous authors, are common on the northern plains of Mars. Our analysis of selected regions in southern Acidalia Planitia has revealed over 18,000 such features, and we estimate that more than 40,000 occur across the area. These domes and cones strongly resemble terrestrial mud volcanoes in size, shape, morphology, associated flow structures and geologic setting. Geologic and mineralogic arguments rule out alternative formation mechanisms involving lava, ice and impacts. We are studying terrestrial mud volcanoes from onshore and submarine locations. The largest concentration of onshore features is in Azerbaijan, near the western edge of the Caspian Sea. These features are typically hundreds of meters to several kilometers in diameter, and tens to hundreds of meters in height. Satellite images show spatial densities of 20 to 40 eruptive centers per 1000 square km. Many of the features remain active, and fresh mud flows as long as several kilometers are common. A large field of submarine mud volcanoes is located in the Gulf of Cadiz, off the Atlantic coasts of Morocco and Spain. High-resolution sonar bathymetry reveals numerous km-scale mud volcanoes, hundreds of meters in height. Seismic profiles demonstrate that the mud erupts from depths of several hundred meters. These submarine mud volcanoes are the closest morphologic analogs yet found to the features in Acidalia Planitia. We are also conducting laboratory analyses of surface samples collected from mud volcanoes in Azerbaijan, Taiwan and Japan. X-ray diffraction, visible / near infrared reflectance spectroscopy and Raman spectroscopy show that the samples are dominated by mixed-layer smectite clays, along with quartz, calcite and pyrite. Thin section analysis by optical and scanning electron microscopy confirms the mineral identifications. These samples also contain chemical and morphological biosignatures, including common microfossils, with evidence of partial replacement by pyrite. The bulk samples contain approximately 1 wt% total organic carbon and 0.4 mg / gm volatile hydrocarbons. The thousands of features in Acidalia Planitia cited as analogous to terrestrial mud volcanoes clearly represent an important element in the sedimentary record of Mars. Their location, in the distal depocenter for massive Hesperian-age floods, suggests that they contain fine-grained sediments from a large catchment area in the martian highlands. We have proposed these features as a new class of exploration target that can provide access to minimally-altered material from significant depth. By analogy to terrestrial mud volcanoes, these features may also be excellent sites for the sampling martian organics and subsurface microbial life, if such exist or ever existed.

Allen, Carlton C.↗
NASA NTRS2013

The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling

The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.

Maggioni, V.↗
DOE OSTI2026

Model scripts associated with “Revisiting controls on hyporheic respiration with knowledge-guided machine learning at continental scale”

NOTE: The manuscript associated with this data package is currently in review. The data/scripts may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final scripts and additional metadata. This data package is associated with the publication “Revisiting controls on hyporheic respiration with knowledge-guided machine learning at continental scale” submitted to Environmental Science & Technology (Zheng et al. 2026). The project combines mechanistic process modeling with knowledge-guided machine learning (KGML) to evaluate how organic matter chemistry, microbial biomass, and physical substrate accessibility regulate realized respiration rates across river corridors. All data used in this paper have been previously published and can be accessed at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1729719 (Goldman et al., 2020). This data package contains 3 R-markdown (Rmd) preprocessing scripts for the previously published data and subsequent modelling workflows. The full workflow with input and output data can be found in the associated GitHub repository at https://github.com/jianqiuz/KGML-WHONDRS.

Biogeochemistry↗
DOE OSTI2025

Water Stable Isotopes in Precipitation, Rivers, and Groundwater Across an Elevation Gradient in the Sierra Nevada Mountains (USA) Reflect Source Elevation

Understanding watershed processes is critical to predict the impacts of climate change and forest management on water resources. However, collecting hydrological data in mountainous terrain is challenging. Precipitation, river water, and groundwater H and O stable isotope data can provide insights into processes occurring at the mountain range scale. Water δ 2 H and δ 18 O values in precipitation vary with terrain elevation; thus, the resulting isotopic lapse rates of precipitation, groundwater, and river water have the potential to elucidate watershed processes and source elevations of major rivers. We analysed H and O stable isotope data of precipitation, groundwater, and river water over the course of one Water Year (Oct 2016—Oct 2017) in the Sierra Nevada mountains of California, USA. We calculated elevation-dependent isotopic lapse rates of these waters to estimate the source elevation of major rivers draining the west flank of the Sierra Nevada mountains. We also investigated the Cosumnes River's watershed in more detail to determine how river flow may be more fully partitioned. Here, we found that H and O stable isotopes in precipitation are temporally variable, but isotopic lapse rates are generally consistent with prior studies. However, groundwater samples across an elevation gradient provide a more consistent and accessible isotopic lapse rate to predict river water source elevations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Examination of L-Band Brightness Temperature Forecasts in the SMAP Level-4 Soil Moisture Analysis

The NASA Soil Moisture Active Passive (SMAP) mission [4] has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. By assimilating the Tb observations into the NASA Catchment land surface model [5] using a spatially distributed ensemble Kalman filter, the NASA Global Modeling and Assimilation Office generates the SMAP Level-4 Soil Moisture (L4_SM) product, which provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with ~2.5-day latency for use in research and applications [6]. The L4_SM product also includes estimates of soil temperature, land surface fluxes, and assimilation diagnostics such as the model forecast and observed Tb values [7]. The output from the L4_SM system is routinely monitored by the L4_SM team. Such monitoring provides valuable information; instances of unusually large Tb observation-minus-forecast (O-F) residuals can indicate events during which soil moisture conditions are poorly described in the land modeling system [7,8]. For example, repeated occurrences of very large Tb O-F values in central Australia were traced back to deficiencies in the gauge-based precipitation product used in the land modeling system through L4_SM Version 5; this discovery prompted the use, beginning in Version 6, of satellite- and gauge-based precipitation observations outside of North America [9]. Recently, a statistical analysis of the Tb O-F residuals revealed systematic errors in the “tau-omega” L-band radiative transfer model that converts the land-model simulated soil moisture and temperature into the forecast Tb prior to the L4_SM analysis. The L4_SM analysis is built on the ensemble Kalman filter and assumes unbiased forecast errors. The presence of systematic Tb forecast errors could thus adversely impact the quality of the analyzed soil moisture. In this paper, we examine the Tb O-F residuals of the latest Version 7 L4_SM data (Science Version ID Vv7030 and Vv7032) [10,11,12]. In the L4_SM Version 7 algorithm, key parameters of the L-band radiative transfer model, including soil roughness, scattering albedo, and a (seasonally varying) climatology of vegetation opacity, are obtained from the SMAP Level-2 Radiometer retrieval product (Version 5) [1]. Additionally, we also use ground measurements of surface soil moisture and soil temperature from the SMAP core validation sites as in situ reference of soil conditions [2,3]. As part of the L4_SM system calibration, the seasonally-varying bias between the model forecast Tb and the observed values is removed prior to the assimilation of the SMAP Tb observations [7]. That is, the L4_SM system is designed to only correct errors in synoptic-scale and interannual variations from the long-term mean seasonal cycle while maintaining the model’s (potentially erroneous) climatology. The Tb O-F residuals examined here are thus computed using SMAP Tb observations after they are rescaled to the mean seasonal cycle of Tb from the modeling system without SMAP data assimilation. Consequently, the long-term average of the Tb O-F residuals has a global mean of only 0.13 K and locally small values, ranging from -1 to 3 K. Despite the small time-average values of the Tb O-F residuals, the model forecast Tb was nevertheless found to exhibit undesirable systematic errors. At some locations, the time-average Tb O-F values strongly depend on surface soil moisture (SM). At the Yanco SMAP core validation site, for example, the Tb O-F residuals typically range from 5 to 15 K under dry soil moisture conditions (SM < 0.15 m3 m-3) yet are predominantly negative under wet soil moisture conditions (SM > 0.25 m3 m-3), with values ranging from 0 to -40 K. This results in soil moisture analysis increments that persistently make the soil drier under dry SM conditions and persistently make the soil wetter under wet SM conditions, suggesting an error in the dynamic range of the simulated Tb, soil moisture or soil temperature. In this paper, we describe the higher-order systematic Tb forecast errors in more detail, examine their impact on the L4_SM product quality, and explore potential avenues to further improve the L4_SM algorithm.

Rolf Reichle↗
DOE OSTI2024

Laboratory time series moisture manipulative experiment from sediment across the contiguous US: time series aerobic respiration and geochemistry (v2)

This dataset supports a broader study examining the effects of wetting and drying on hyporheic zone respiration across the contiguous United States (CONUS). The dataset provides data generated from a laboratory moisture manipulation experiment. The contents include time series aerobic respiration and moisture; dissolved oxygen; sediment geochemistry data; and field metadata (including qualitative information on instream and river corridor characteristics). Samples were collected as part of the WHONDRS CONUS-Scale Model-Sample Study (CM). This study was designed following ICON (integrated, coordinated, open, and networked) principles to facilitate a model-experiment (ModEx) iteration approach, leveraging crowdsourced sampling across the CONUS. The data package associated with the CM study is available at https://data.ess-dive.lbl.gov/view/doi:10.15485/1923689. CM sampling began in April 2022 and ended in October 2023. This study uses subsamples from a subset of CM samples collected between June 2022 and June 2023. The original field samples were labeled as CM_###. Subsequent subsamples for this study were labeled as EC_###. The labels from the field samples and the EC subsamples can be mapped directly based on the digits following the prefix and underscore (i.e., EC_001 is a subsample from CM_001). See the critical details section below for more details on sample naming. This data package was originally published in August 2024. It was updated in February 2026 (v2; new and modified files). See the change history section in the readme for more details. For details on how to navigate this data package, see this infographic from the River Corridor SFA https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data and one main data folder containing (1) file-level metadata; (2) data dictionary; (3) field metadata; (4) readme; (5) field protocol; and a (6) a subfolder with sediment sample data from the incubation experiment. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC); (2) total nitrogen (TN); (3) adenosine triphosphate (ATP); (4) percent carbon and nitrogen; (5) effect size; (6) iron (II); (7) gravimetric moisture; (8) respiration rates and raw dissolved oxygen values; (9) specific conductance; (10) pH; (11) temperature; (12) a summary containing median values of each data type for each treatment (wet and dry); (13) methods codes; (14) FTICR-MS methods; and (15) a subfolder of 9.4 Tesla FTICR-MS data. This folder contains three subfolders, one containing the sediment .xml data files, one containing the sediment CoreMS output files, the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .ref, or .xml.

54 ENVIRONMENTAL SCIENCES↗