Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “watershed modeling”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Overview of the NASA Earth Action Strategies Wildland Fire Initiative

As part of NASA’s new Earth Action strategy, the Wildland Fire initiative was established, which includes both the NASA Wildland Fire Program (WFP) and the FireSense project. NASA has over 50 years of experience generating data and technology to enhance fire science and operational management. The WFP’s mission is threefold: 1) assemble communities of practice through collaborative efforts with government, academia, and the private sector; 2) co-develop knowledge and applications with relevant partners and stakeholders in the wildfire community; and 3) improve wildland fire management through the transitioning of NASA data, technology, tools, and science to stakeholder organizations. The WFP is focusing on supporting proactive fire management, including situational awareness, preparedness, and risk mitigation. This will be accomplished through selected projects that identify management challenges, relevant to partners and end users, and the NASA data that will be utilized to deliver innovative solutions to enhance the management of wildland fires. Examples include: i) investigation of evaporative stress from OpenET to help predict the risk of wildfire occurrence in watersheds; ii) incorporation of space based LiDAR for the generation of 3-dimensional forest fuel metrics, used to improve wildfire risk and behavior models; iii) integration of global, multi-platform geostationary active-fire data in near-real-time into NASA’s Fire Information for Resource Management System (FIRMS); and iv) identification of post-fire ecohydrological conditions using thermal, multispectral, synthetic aperture radar (SAR), and hyperspectral remotely-sensed data to improve flood hazard forecast models. The FireSense project is a US-focused 5-year project that will focus on delivering NASA’s unique Earth science and technological capabilities to operational agencies, striving towards enhancing fire fighting and air quality management. The project will include airborne campaigns and new technology that will likely have global implications. Initial stakeholder engagement led FireSense to focus on four use-cases focused on the characterization and measurement of: (i) pre-fire fuels conditions, (ii) active fire-dynamics; (iii) post-fire impact and threats; and iv) air quality impacts and forecasting, each-developed with identified stakeholders.

Wildland Fire program↗
DOE OSTI2024

Delineation of endorheic drainage basins in the MERIT-Plus dataset for 5 and 15 minute upscaled river networks

Abstract The MERIT-Hydro networks re-gridded by the Iterative Hydrography Upscaling (IHU) algorithm do not retain exo- or endorheic basin attributes from the original data. Here we developed methods to assign such attributes to those and any other digital river networks. The motivation is that endorheic inland drainage basins are essential for hydrologic modelling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. First, we create basin attributes that explicitly label endorheic and exorheic catchments by the criteria of direct or hidden connectivity to the ocean without changing their flow direction grid. In the second step we alter the delineation of endorheic basins by the merging algorithm that eliminates small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion and topology of the inland basins. The data was validated using the Water Balance Model by comparing volume of endorheic inland depressions with modelled water accumulation in their inland lakes.

58 GEOSCIENCES↗
DOE OSTI2023

Projected changes in extreme streamflow and inland flooding in the mid-21st century over Northeastern United States using ensemble WRF-Hydro simulations

Study region: Northeastern United States (NEUS). Study focus: We investigate the potential impacts of climate change on precipitation, streamflow, and inland flooding in the NEUS during the mid-21st century. Dynamically downscaled climate projections from three global climate models for historical (1995-2004) and future (2045-2054) periods under business-as-usual scenarios were used to force the hydrologic model WRF-Hydro at 200-meter resolution and create ensemble hydrologic simulations. Additionally, an extreme value model was developed to project the risks associated with low-frequency hydrologic events. New hydrological insights for the region: Results from four major watersheds indicate a significantly wetter regime in winter months and potential drier conditions during late spring to early summer. Discharges in fall are projected to decrease in the northern watersheds and increase toward the south. Extreme flow and water depths resulting from extreme inland flooding are projected to increase by 5-20% and > 100%, respectively. The extent of the total flooded area is likely to be 20% greater by the mid-century. These increased risks can be attributed to (i) an approximate 25% increase in decadal mean and > 40% increase in decadal extreme precipitation intensity, (ii) up to 30% lower snow availability and 5-25% higher evapotranspiration throughout the year, and (iii) a projected 5% increase in soil moisture in all seasons except summer. Furthermore, rapid snow melting in winter will likely cause an earlier peak flow in the rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advancing process-based flood frequency analysis for assessing flood hazard and population flood exposure

Recent studies have showcased the use of process-based hydrological models with Stochastic Storm Transposition (SST) techniques to conduct Flood Frequency Analysis (FFA). This framework, referred hereby FFA-SST, has proved to be a robust strategy to estimate peak flows of specific annual exceedance probability (e.g., 100-year peak flow) that can reflect natural and anthropogenic disturbances, including changes in land use and meteorological patterns. With the objective of advancing the FFA-SST framework, this study presents for the first time the use of an Integrated Surface-Subsurface Hydrological Model (ISSHM) to conduct FFA-SST by extending the analysis from peak flow responses to flood extent, enabling a unique view and analysis of flood hazard and population flood exposure at the basin scale. As a proof-of-concept, we used the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct FFA-SST by simulating the flood response to 5,000 annual synthetic storm events in a 2,227 $km^2$ Southeast Texas watershed. We demonstrate that ATS, without site-specific calibration, provides a robust process-based representation of peak flows, flood extent, streamflow, evapotranspiration, soil moisture content, and water storage changes. Our results and analyses, covering frequency curves up to a 500-year return period for peak flows, basin inundation fractions, and the number of people exposed to flooding, offer a unique perspective to analyze flood impacts across spatial scales. Overall, this study provides critical insights for flood risk management by extending the FFA-SST framework to include both flood hazard and population flood exposure analyses at the basin scale. Such an approach will empower stakeholders and disaster emergency agencies with a more comprehensive understanding of flood impacts across the entire basin domain, facilitating informed decision-making for flood risk assessment and management.

58 GEOSCIENCES↗
NASA NTRS1986

The effect of Thematic Mapper spectral properties on land cover mapping for hydrologic modeling

The accuracy of unsupervised land-cover classification from all seven Landsat TM bands and from six combinations of three or four bands is evaluated using images of the Clinton River Basin, a suburban watershed near Detroit. Data from aerial TMS photography, USGS topographic maps, and ground surveys are employed to determine the classification accuracy. The mapping accuracy of all seven bands is found to be significantly better (6 percent overall, 12 percent for residential areas, and 13 percent for commercial districts) than that with bands 2, 3, and 4; but almost the same accuracy is obtained by including at least one band from each major spectral region (visible, NIR, and mid-IR).

Gervin, J. C.↗
DOE OSTI2024

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗
DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

TDD LoRa and Delta Encoding in Low-Power Networks of Environmental Sensor Arrays for Temperature and Deformation Monitoring: Supporting Data

This dataset was used to assess the performance of a compression algorithm for Soil Deformation and Temperature Monitoring data. The development of this technology is presented in the article named time division duplex (TDD) "TDD LoRa and Delta Encoding in Low-Power Networks of Environmental Sensor Arrays for Temperature and Deformation Monitoring" and published in Springer's Journal on Signal Processing Systems. [Time Division Duplex (TDD); Long Range (LoRa)]Field experiment data includes: 1) soil temperature measurements, 2) snow temperature data, 3) soil temperature and deformation data. In these files, the first column contains the time (UTC), followed by the sensor's battery voltage and temperature and acceleration values (X, Y, Z) in subsequent columns. The measurements were acquired every 15 minutes. Lists included for probe, datalogger, and packaging components. The two temperature data sets were acquired at a site at the East River Watershed, Colorado, while the temperature/deformation data were collected at the Teller Road (mile 27) site near Nome, Alaska. The package includes 8 comma-separated files (.csv) and one .pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Multi-scale impacts of climate change on hydropower for long-term water-energy planning in the contiguous United States

Climate change impacts on watersheds can potentially exacerbate water scarcity issues where water serves multiple purposes including hydropower. The long-term management of water and energy resources is still mostly approached in a siloed manner at different basins or watersheds, failing to consider the potential impacts that may concurrently affect many regions at once. There is a need for a large-scale hydropower modeling framework that can examine climate impacts across adjoining river basins and balancing authorities (BAs) and provide a periodic assessment at regional to national scales. Expanding from our prior assessment only for the United States (US) federal hydropower plants, we enhance and extend two regional hydropower models to cover over 85% of the total hydropower nameplate capacity and present the first contiguous US-wide assessment of future hydropower production under Coupled Model Intercomparison Project phase 6’s high-end Shared Socioeconomic Pathway 5-8.5 emission scenario using an uncertainty-aware multi-model ensemble approach. We present regional hydropower projections, using both BA regions and US Hydrologic Subregions (HUC4s), to consistently inform the energy and water communities for two future periods—the near-term (2020–2039) and the mid-term (2040–2059) relative to a historical baseline period (1980–2019). We find that the median projected changes in annual hydropower generation are typically positive—approximately 5% in the near-term, and 10% in the mid-term. However, since the risk of regional droughts is also projected to increase, future planning cannot overly rely on the ensemble median, as the potential of severe hydropower reductions could be overlooked. The assessment offers an ensemble of future hydropower generation projections, providing regional utilities and power system operators with consistent data to develop drought scenarios, design long duration storage and evaluate energy infrastructure reliability under intensified inter-annual and seasonal variability.

13 HYDRO ENERGY↗
NASA NTRS2014

An Efficient, Multi-Layered Crown Delineation Algorithm for Mapping Individual Tree Structure Across Multiple Ecosystems

Deriving individual tree information from discrete return, small footprint LiDAR data may improve forest above ground biomass estimates, and provide tree-level information that is important in many ecological studies. Several crown delineation algorithms have been developed to extract individual tree information from LiDAR point clouds or rasterized canopy height models (CHM), but many of these algorithms have difficulty discriminating between overlapping crowns, and also may fail to detect understory trees. Our approach uses a watershed based delineation of a CHM, which is subsequently refined using the LiDAR point cloud. Individual tree detection was validated with stem mapped field data from the Smithsonian Environmental Research Center (SERC), Maryland, and on a plot and stand level through comparisons of stem density and basal area to delineated metrics at both SERC and a study area in the Sierra Nevada, California. For individual tree detection, the algorithm correctly identified 70% of dominant trees, 58% of co-dominant trees, 35% of intermediate trees and 21% of suppressed trees at SERC. The algorithm had difficulty distinguishing between crowns of small, dense understory trees of approximately the same height. Delineated crown volume alone explained 53% and 84% of the variability in basal area at the SERC and Sierra Nevada sites, respectively. The algorithm produced crown area distributions comparable to diameter at breast height (DBH) size class distributions observed in the field in both study sites. The algorithm detected understory crowns better in the conifer-dominated Sierra Nevada site than in the closed-canopy deciduous site in Maryland. The ability for the algorithm to reproduce both accurate tree size distributions and individual crown geometries in two dissimilar and complex forests suggests great promise for applicability to a wide range of forest systems.

LiDAR↗
DOE OSTI2021

Probabilistic Evaluation of Geoscientific Hypotheses with Geophysical Data: Application to Electrical Resistivity Imaging of a Fractured Bedrock Zone

As climate changes and populations grow, groundwater sustainability is becoming increasingly important. Groundwater models, based on a conceptual understanding of the subsurface structure, are crucial tools for making sustainable management decisions. Conceptual models of the subsurface are based on knowledge of geological processes, and, frequently, observations from geophysical data. A frequent problem in groundwater model development occurs when multiple geological phenomena could explain a single subsurface observation. Uncertainty in geophysical data makes it even more difficult to discern which explanations are consistent with the geophysics. Here, we present a framework for testing geological when a geological feature is observed in geophysical data, but its physical characteristics are uncertain. The framework builds on Popper-Bayes methods developed in previous work, and is applied to study a fractured bedrock zone in a mountainous watershed in southwest Colorado. First, we outline six hypotheses based on the geological history of the watershed. Then, using the proposed Popper-Bayes approach, we demonstrate that three of the six hypotheses are inconsistent with measured electrical resistivity data, even after accounting for uncertainty. Finally, we discuss the importance of the prior model, and how this framework for handling geophysical uncertainty can be applied in other settings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Designing Iowa Agricultural Landscapes to Improve Environmental Co-Benefits of Bioenergy Production

Cellulosic bioenergy feedstocks are needed to improve carbon (C) management while provisioning biomass for bioproducts and biofuel. The transition to increased cellulosic biomass production can be guided by land management plans designed to improve economic, environmental, and ecological performance. We constructed a sustainability model to compare landscape designs for biofuel production from corn (Zea mays L.) stover and switchgrass (Panicum virgatum L.) in central Iowa, USA. We used the model to compare environmental and socioeconomic outcomes associated with four landscape management strategies, with and without cellulosic biomass markets. We evaluated (1) a fuelshed area containing over 1.2 million ha (3 million acres) of corn and soybean (Glycine max (L.) Merr.) within 80 km (50 miles) of a commercial-scale cellulosic biorefinery in Nevada, Iowa, and (2) the South Fork watershed containing over 72,000 ha (178,000 acres) of these row crops within eight north central Iowa HUC-12 (hydrologic unit code) watersheds. At both landscape scales, we found that it is possible to achieve multiple environmental and socioeconomic benefits concomitantly with cellulosic biomass production by strategically collecting corn stover and converting the 10% of the lowest-profitability row crop land to perennial switchgrass. Potential benefits from landscape design include increased biodiversity, soil and water quality improvements, increased soil carbon sequestration for climate change mitigation, and reduced fertilizer use and cost. Our model results showed that increasing benefits can accrue when complementary conservation practices (e.g., reduced tillage, use of a rye cover crop) are combined and integrated throughout a fuelshed or watershed area. We conclude that ecologically based landscape designs offer valuable insights about costs and benefits of land management alternatives, with relevance for achieving stakeholder goals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

How extreme rainfall and failing dams unleashed the Derna flood disaster

On September 11, 2023, Storm Daniel unleashed unprecedented rainfall over the Wadi Derna watershed, triggering one of the most devastating floods in modern history, striking Derna, a coastal city in Libya. This study reconstructs the disaster using an integrated modeling approach that combines satellite imagery, hydrologic, hydraulic, and geotechnical simulations, machine learning, eyewitness accounts, and digital elevation data to assess the impact of cascading dam failures. Our findings reveal that the region’s dams, even if structurally sound, would have provided minimal protection against the extreme runoff. However, their failure unleashed a destructive surge wave, amplifying the disaster’s magnitude and devastation. Here, we show that the collapse of aging flood control infrastructures, compounded by inadequate risk assessment and emergency preparedness, dramatically escalated the disaster’s impact. Our findings underscore the urgent need for systematic dam safety evaluations, enhanced flood forecasting, and adaptive risk management strategies that address climate extremes and infrastructure vulnerabilities.

Hydrology↗
DOE OSTI2023

Relationship between isotope ratios in precipitation and stream water across watersheds of the National Ecological Observation Network

The timescales associated with precipitation moving through watersheds reveal processes that are critical to understanding many hydrologic systems. Measurements of environmental stable water isotope ratios (δ 2 H and δ 18 O) have been used as tracers to study hydrologic timescales by examining how long it takes for incoming precipitation tracers become stream discharge, yet limited measurements both spatially and temporally have bounded macroscale evaluations so far. Here in this observation driven study across North American biomes within the National Ecological Observation Network (NEON), we examined δ 18 O and δ 2 H stable water isotope in precipitation (δP) and stream water (δQ) at 26 co-located sites. With an average 54 precipitation samples and 139 stream water samples per site collected over 2014–2022, assessment of local meteoric water lines and local stream water lines showed geographic variation across North America. Taking the ratio of estimated seasonal amplitudes of δP and δQ to calculate young water fractions (F yw ), showed a Fyw range from 1% to 93% with most sites having F yw below 20%. Calculate d mean transit times (MTT) based on a gamma convolution model showed a MTT range from 0.10 to 13.2 years, with half of the sites having MTT estimates lower than 2 years. Significant correlations were found between the F yw and watershed area, longest flow length, and the longest flow length/slope. Significant correlations were found between MTT and site latitude, longitude, slope, clay fraction, temperature, precipitation magnitude, and precipitation frequency. The significant correlations between water timescale metrics and the environmental characteristics we report share some similarities with those reported in prior studies, demonstrating that these quantities are primarily driven by site or area specific factors. The analysis of isotope data presented here provides important constraints on isotope variation in North American biomes and the timescales of water movement through NEON study sites.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Comparing Landsat-7 ETM+ and ASTER Imageries to Estimate Daily Evapotranspiration Within a Mediterranean Vineyard Watershed

We compared the capabilities of Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imageries for mapping daily evapotranspiration (ET) within a Mediterranean vineyard watershed. We used Landsat and ASTER data simultaneously collected on four dates in 2007 and 2008, along with the simplified surface energy balance index (S-SEBI) model. We used previously ground-validated good quality ASTER estimates as reference, and we analyzed the differences with Landsat retrievals in light of the instrumental factors and methodology. Although Landsat and ASTER retrievals of S-SEBI inputs were different, estimates of daily ET from the two imageries were similar. This is ascribed to the S-SEBI spatial differencing in temperature, and opens the path for using historical Landsat time series over vineyards.

Temperature measurement↗
NASA NTRS2020

Investigation of Magmatic Activities on Early Mars Using Igneous Mineral Chemistry in Gale Crater, Mars

One objective of rover missions is exploring the geological context of the surroundings. Over the years, igneous petrology and sedimentology have been disconnected, the first investigating magmatic processes and volcanic activities, and the second seeking environmental conditions in the past and assessing the habitability of the planet. Although different, one is related to the other: igneous rocks are altered and broken down, leading to the formation of sedimentary rocks, which can in turn be used to back out the nature of their magmatic source. The Curiosity rover that landed in the 3.7 Gyr old impact crater Gale is traveling through sedimentary rocks. About fifty float rocks have been observed, and several of them with ambiguous texture and composition have been classified as igneous or sedimentary depending on studies such as Jake_M. The composition of several unambiguous igneous rocks has been analyzed [4- 6] but their heterogeneity at a larger instrumental (measurement size < 2 cm) scale prevents the measurement of a bulk composition as performed on Earth. An original approach avoiding these two last issues is to consider igneous mineral chemistry analyzed within igneous and sedimentary rocks to assess magmatic processes that could have formed them. Most Curiosity data are used to explore ancient environmental conditions, and a significant number of compositional analyses are under-explored for constraining magmatic activities. We will present how we can make use of sedimentary data for investigating igneous processes in the vicinity of Gale crater. Geological Context: We focus on the first 750 martian days, corresponding to measurements in a coherent lacustrine sedimentary unit called Bradbury, because all sedimentary rocks were sourced from the same watershed and appear to have a consistent source with minimal alteration [2-3]. Igneous detrital minerals including feldspar and pyroxene, are observed in sedimentary rocks. Monte Carlo models showed that minimal cation loss is observed based on the composition of all Bradbury rocks, implying negligible weathering [3]. Although clay minerals are detected in few rocks [7], chemical compositions of rocks can be explained by a mixture of primary igneous minerals [3]. Variation of composition within Bradbury rocks can be explained by mineral sorting and one distinct source component. While a common magmatic source is suggested, Bradbury sediments likely come from several volcanic eruptions from a single magmatic chamber [9- 10]. The occurrence of alkali minerals like sanidine and K-rich rocks throughout Bradbury supports the presence of a potassic component, likely trachytic, while plagioclase and a mafic composition suggest a basaltic component [8-9]. Instruments: Mineral chemistry can be estimated by three instruments onboard Curiosity. The CheMin instrument enables detection of mineral assemblages using X-ray diffraction (XRD). Using Rietveld refinement, each mineral is identified according to their 1D XRD pattern [11]. Note that distinction between pyroxene minerals is challenging with the CheMin instrument due to overlapping peaks on XRD patterns and low angular resolution of the instrument [12]. Then, using least square regression and optimization algorithms based on unit-cell parameters, mineral chemistry has been estimated by [11]. Plagioclase composition has been estimated using the NaAlSi3O8- CaAl2Si2O8 system and alkali feldspar is based on the NaAlSi3O8-KAlSi3O8 system (stars in Fig.1). Two mudstone samples (John Klein and Cumberland) and one sandstone sample (Windjana) were analyzed by CheMin at Yellowknife Bay and Kimberley, respectively. Figure 1. Ternary diagrams of feldspars (top) and pyroxene (bottom) quadrilateral. Stars correspond to CheMin composition and the gray patches to ChemCam composition. The colored dots are the composition of feldspar and pyroxene that crystallized during fractional crystallization at FMQ+1 of a melt extracted at distinct melting degree during the adiabatic ascent of a primitive mantle composition, without any water (left panels) and with 0.5 wt.% of water (right panels) at distinct pressure. The ChemCam instrument enables the analysis of the chemical compositions of rocks at hundreds of micrometer scale (350-550 μm) using laser induced breakdown spectroscopy (LIBS), which may provide the composition of minerals when they are larger than the beam spot (>550 μm) [13]. Within >5000 LIBS points, we performed a typical stoichiometric filtering allowing us to distinguish 56 feldspar and 10 pyroxene mineral compositions (grey patches in Fig. 1). Finally, the Alpha Particle X-ray Spectrometer (APXS) analyzes the composition of rocks with a 1.6 cm diameter spot size. Monte Carlo mass balance modeling allowed [3] to decipher a feldspar range varying between An30 and An40 (Fig. 1). Discussion: Although there could be a more complex history and other ways to form the whole compositional range of igneous minerals analyzed within the Bradbury formation, we are presenting here simple magmatic pathways commonly occurring on Earth using the thermodynamical softwares pMELTS and rhyoliteMELTS [14]. The objective is to find reasonable igneous processes that produce minerals that parallel the compositions of feldspar and pyroxene analyzed by the Curiosity rover. As commonly observed for mid-ocean ridge basalts, the adiabatic ascent of a primitive mantle composition [15] partially melting at 2 GPa has been modeled, followed by the extraction of a liquid at distinct degrees of partial melting, which undergoes fractional crystallization at an oxygen fugacity +1 log unit above the fayalite-magnetite-quartz (FMQ) buffer within the crust (0.02-0.4 GPa) with H2O = 0- 0.5 wt. %. These latter conditions correspond to those recorded within igneous clasts from the Noachian martian breccia NWA 7034 and paired and within Gale igneous rocks (colored dots in Fig. 1) [16-17]. To check the reliability of these 2-step models, we also tested fractional crystallization at similar conditions (FMQ+1; P=0.02-0.4 GPa; H2O = 0-0.5 wt. %) of starting compositions corresponding to that of magmas with distinct melting degrees obtained from isobaric experiments at 2 GPa [18]. Mineral compositions obtained from both models are similar. As shown on Fig. 1, the whole range of observed feldspar compositions cannot be reproduced by fractionation of one magma only. Indeed, while alkali feldspar and Na-plagioclase likely crystallized from fractional crystallization of a low-degree melt (here <15%), plagioclase and pyroxene can only be formed by fractional crystallization of a higher degree melt (here >19%). The corresponding liquid descent lines are broadly in agreement with compositions estimated by ChemCam corresponding to float igneous rocks (Fig. 2) [4-6]. Figure 2. Silica versus alkali content. Lines show the liquid lines of descent from magmas with distinct degrees of melting. Gray patches represents the composition of Gale igneous rocks [4-6]. Trachytic to rhyolitic magmas crystallize alkali feldspar, and andesite to dacite magmas likely form plagioclase. Therefore, at least two starting magmas at distinct melting degrees, which could easily come from a single mantle source, are necessary to explain the whole compositional range of feldspar and pyroxene analyzed within Bradbury rocks. Conclusion: Because rocks from the Bradbury formation are likely originating from the same magmatic source with minimal weathering as supported by several studies using different approaches, igneous mineral chemistry analyzed by CheMin and ChemCam allows us to back out reasonable magmatic pathways that could have crystallized them. Fractional crystallization of at least two starting magmas originating from distinct melting degrees of a single mantle source can explain the whole range of feldspar and pyroxene composition. Both alkaline and sub-alkaline liquids can be produced, with compositions corresponding to those of the igneous rocks analyzed by ChemCam within the Bradbury formation, highlighting the complexity of Mars magmatism.

Payre, V.↗
DOE OSTI2021

Hidden Processes During Seasonal Isolation of a High-Altitude Watershed

Biogeochemical processes capable of altering global carbon systems occur frequently in Earth’s Critical Zone–the area spanning from vegetation canopy to saturated bedrock–yet many of these phenomena are difficult to detect. Observation of these processes is limited by the seasonal inaccessibility of remote ecosystems, such as those in mountainous, snow- and ice-dominated areas. This isolation leads to a distinct gap in biogeochemical knowledge that ultimately affects the accuracy and confidence with which these ecosystems can be computationally modeled for the purpose of projecting change under different climate scenarios. To examine a high-altitude, headwater ecosystem’s role in methanogenesis, sulfate reduction, and groundwater-surface water exchange, water samples were continuously collected from the river and hyporheic zones (HZ) during winter isolation in the East River (ER), CO watershed. Measurements of continuously collected ER surface water revealed up to 50 μM levels of dissolved methane in July through September, while samples from 12 cm deep in the hyporheic zone at the same location showed a spring to early summer peak in methane with a strong biogenic signature (<65 μM, δ13C-CH4, -60.76‰) before declining. Continuously collected δ 18 O-H 2 O and δ 2 H-H 2 O isotopes from the water column exhibited similar patterns to discrete measurements, while samples 12 cm deep in the hyporheic zone experienced distinct fluctuations in δ 18 O-H 2 O, alluding to significant groundwater interactions. Continuously collected microbial communities in the river in the late fall and early winter revealed diverse populations that reflect the taxonomic composition of ecologically similar river systems, including taxa indicative of methane cycling in this system. These measurements captured several biogeochemical components of the high-altitude watershed in response to seasonality, strengthening our understanding of these systems during the winter months.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Hysteresis Patterns of Watershed Nitrogen Retention and Loss Over the Past 50 years in United States Hydrological Basins

Patterns of watershed nitrogen (N) retention and loss are shaped by how watershed biogeochemical processes retain, biogeochemically transform, and lose incoming atmospheric deposition of N. Loss patterns represented by concentration, discharge, and their associated stream exports are important indicators of integrated watershed N retention behaviors. We examined continental United States (CONUS) scale N deposition (e.g., wet and dry atmospheric deposition), vegetation trends, and stream trends as potential indicators of watershed N-saturation and retention conditions, and how watershed N retention and losses vary over space and time. By synthesizing changes and modalities in watershed nitrogen loss patterns based on stream data from 2200 U.S. watersheds over a 50 years record, our work revealed two patterns of watershed N-retention and loss. One was a hysteresis pattern that reflects the integrated influence of hydrology, atmospheric inputs, land-use, stream temperature, elevation, and vegetation. The other pattern was a one-way transition to a new state. We found that regions with increasing atmospheric deposition and increasing vegetation health/biomass patterns have the highest N-retention capacity, become increasingly N-saturated over time, and are associated with the strongest declines in stream N exports—a pattern, that is, consistent across all land cover categories. We provide a conceptual model, validated at an unprecedented scale across the CONUS that links instream nitrogen signals to upstream mechanistic landscape processes. Our work can aid in the future interpretation of in-stream concentrations of DOC and DIN as indicators of watershed N-retention status and integrators of watershed hydrobiogeochemical processes.

58 GEOSCIENCES↗