Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vegetation dynamics 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Forest responses to simulated elevated CO 2 under alternate hypotheses of size- and age-dependent mortality

Elevated atmospheric carbon dioxide (eCO 2 ) is predicted to increase growth rates of forest trees. The extent to which increased growth translates to changes in biomass is dependent on the turnover time of the carbon, and thus tree mortality rates. Size- or age-dependent mortality combined with increased growth rates could result in either decreased carbon turnover from a speeding up of tree life cycles, or increased biomass from trees reaching larger sizes, respectively. However, most vegetation models currently lack any representation of size- or age-dependent mortality and the effect of eCO 2 on changes in biomass and carbon turnover times is thus a major source of uncertainty in predictions of future vegetation dynamics. Using a reduced-complexity form of the vegetation demographic model the Functionally Assembled Terrestrial Ecosystem Simulator to simulate an idealised tropical forest, we find increases in biomass despite reductions in carbon turnover time in both size- and age-dependent mortality scenarios in response to a hypothetical eCO 2 -driven 25% increase in woody net primary productivity (wNPP). Carbon turnover times decreased by 9.6% in size-dependent mortality scenarios due to a speeding up of tree life cycles, but also by 2.0% when mortality was age-dependent, as larger crowns led to increased light competition. Increases in aboveground biomass (AGB) were much larger when mortality was age-dependent (24.3%) compared with size-dependent (13.4%) as trees reached larger sizes before death. In simulations with a constant background mortality rate, carbon turnover time decreased by 2.1% and AGB increased by 24.0%, however, absolute values of AGB and carbon turnover were higher than in either size- or age-dependent mortality scenario. The extent to which AGB increases and carbon turnover decreases will thus depend on the mechanisms of large tree mortality: if increased size itself results in elevated mortality rates, then this could reduce by about half the increase in AGB relative to the increase in wNPP.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Three Decades of Wetland Methane Surface Flux Modeling by Earth System Models–Advances, Applications, and Challenges

Earth System Models (ESMs) simulate the exchange of mass and energy between the land surface and the atmosphere, with a key focus on modeling natural greenhouse gas feedbacks. Methane is the second most important greenhouse gas after carbon dioxide. There are growing concerns over the rapidly increasing methane concentration in the atmosphere, underscoring the need for accurate global modeling of its emissions using ESMs. Of the multitude of sources of methane globally, wetlands are the largest natural emitters for methane, leading to significant efforts targeting their representation in ESMs with a special focus on their methane emissions. In this review, we first provide a historical overview of including wetland–methane components in ESMs and how methane modeling approaches have evolved over time. Second, we discuss recent modeling advancements that show promise for improvements in methane emissions predictions, namely the coupling of surface and atmospheric modules of ESMs, the representation of microtopography and transport mechanisms, the resolution of microbial processes at different spatial–temporal scales, and the improved mapping of wetland area extent across the different wetland types. Third, we shed light on the different challenges hindering accurate estimations of wetland–methane emissions, as shown by the consistent discrepancy between bottom–up and top–down models' predictions. Finally, we emphasize that more detailed representation of biogeochemistry and dynamic hydrology while resolving the within–wetland vegetation heterogeneity should improve model predictions, especially when coupled with expanding ground–based measurement networks and high–resolution remote sensing mapping of methane–relevant variables, such as water elevation, water table depth, and methane concentration.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Characterizing 3D Vegetation Structure from Space: Mission Requirements

Human and natural forces are rapidly modifying the global distribution and structure of terrestrial ecosystems on which all of life depends, altering the global carbon cycle, affecting our climate now and for the foreseeable future, causing steep reductions in species diversity, and endangering Earth s sustainability. To understand changes and trends in terrestrial ecosystems and their functioning as carbon sources and sinks, and to characterize the impact of their changes on climate, habitat and biodiversity, new space assets are urgently needed to produce high spatial resolution global maps of the three-dimensional (3D) structure of vegetation, its biomass above ground, the carbon stored within and the implications for atmospheric green house gas concentrations and climate. These needs were articulated in a 2007 National Research Council (NRC) report (NRC, 2007) recommending a new satellite mission, DESDynI, carrying an L-band Polarized Synthetic Aperture Radar (Pol-SAR) and a multi-beam lidar (Light RAnging And Detection) operating at 1064 nm. The objectives of this paper are to articulate the importance of these new, multi-year, 3D vegetation structure and biomass measurements, to briefly review the feasibility of radar and lidar remote sensing technology to meet these requirements, to define the data products and measurement requirements, and to consider implications of mission durations. The paper addresses these objectives by synthesizing research results and other input from a broad community of terrestrial ecology, carbon cycle, and remote sensing scientists and working groups. We conclude that: (1) current global biomass and 3-D vegetation structure information is unsuitable for both science and management and policy. The only existing global datasets of biomass are approximations based on combining land cover type and representative carbon values, instead of measurements of actual biomass. Current measurement attempts based on radar and multispectral data have low explanatory power outside low biomass areas. There is no current capability for repeatable disturbance and regrowth estimates. (2) The science and policy needs for information on vegetation 3D structure can be successfully addressed by a mission capable of producing (i) a first global inventory of forest biomass with a spatial resolution 1km or finer and unprecedented accuracy (ii) annual global disturbance maps at a spatial resolution of 1 ha with subsequent biomass accumulation rates at resolutions of 1km or finer, and (iii) transects of vertical and horizontal forest structure with 30 m along-transect measurements globally at 25 m spatial resolution, essential for habitat characterization. We also show from the literature that lidar profile samples together with wall-to53 wall L-band quad-pol-SAR imagery and ecosystem dynamics models can work together to satisfy these vegetation 3D structure and biomass measurement requirements. Finally we argue that the technology readiness levels of combined pol-SAR and lidar instruments are adequate for space flight. Remaining to be worked out, are the particulars of a lidar/pol-SAR mission design that is feasible and at a minimum satisfies the information and measurement requirement articulated herein.

Hall, Forrest G.↗
DOE OSTI2021

Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate

Abstract Plant growth and distribution in high-latitude tundra ecosystems is strongly limited by nutrient availability and is critical for quantifying centennial-scale carbon-climate interactions. However, land model representations of plant–nutrient interactions are uncertain, leading to poor comparisons with high-latitude observations. Although it has been recognized for decades in the observational community that plants continue to acquire nutrients well past when aboveground activity has ceased, most large-scale land models ignore this process. Here we address the role tundra plant nutrient acquisition during the non-growing season (NGS) has on centennial-scale vegetation growth and dynamics, with a focus on shrub expansion. We apply a well-tested mechanistic model of coupled plant, microbial, hydrological, and thermal dynamics that explicitly represents nutrient acquisition based on plant and microbial traits, thereby allowing a prognostic assessment of NGS nutrient uptake. We first show that the model accurately represents observed seasonality of NGS plant nutrient uptake in a northern Alaskan tundra site. Applying the model across the North America tundra indicates that NGS nutrient uptake is consistent with observations and ranges between ∼5% and 50% of annual uptake, with large spatial variability and dependence on plant functional type. We show that NGS plant nutrient acquisition strongly enhances modeled 21 st century tundra shrub growth and expansion rates. Our results suggest that without NGS nutrient uptake, total shrub aboveground dominance would be ∼50% lower, limited primarily by their inability to grow tall enough to maximize their inherent capacity for light competition. Evergreen shrubs would be more strongly affected because of their relatively lower capacity for nutrient remobilization and acquisition compared to deciduous shrubs. Our results highlight the importance of NGS plant and soil processes on high-latitude biogeochemistry and vegetation dynamics and motivate new observations and model structures to represent these dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Root Foraging Alters Global Patterns of Ecosystem Legacy From Climate Perturbations

The response of terrestrial ecosystems to climate perturbations typically persist longer than the timescale of the forcing, a phenomenon broadly referred to as legacy. Understanding the strength of legacy is critical for predicting ecosystem sensitivity to climate extremes and the extent that persistent changes in surface-atmosphere exchange might feedback onto climate. The cause of ecosystem legacy has been associated with myriad factors such as changes in aboveground biomass, however, few studies have tested how changes in the depth distribution of fine roots in response to perturbation might alter an ecosystem's recovery time. We explore this question using an Earth System Model that includes a dynamic root module where vegetation can forage for water and nutrients by altering their root profiles. The global simulations presented here show that in response to water stress events most ecosystems deepen their root profiles. In semi-arid ecosystems, this response expedites recovery (i.e., less legacy) relative to simulations without dynamics roots because access to deeper water pools after the initial event remains favorable. In wetter ecosystems, the development of deeper root profiles slows down the recovery timescale (i.e., more legacy) because the deeper root profile reduces access to nutrients and is thus unfavorable. Here, the simulations show that while the inclusion of dynamic roots might only minimally affect global patterns of Gross Primary Productivity and Evapotranspiration, the shift in root profile alters the timescale of recovery. Studies interested in the sustained response of land surfaces fluxes to climate disturbances should consider models that include dynamic root capability.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Incorporating Plant Phenology Dynamics in a Biophysical Canopy Model

The Multi-Layer Canopy Model (MLCan) is a vegetation model created to capture plant responses to environmental change. Themodel vertically resolves carbon uptake, water vapor and energy exchange at each canopy level by coupling photosynthesis, stomatal conductance and leaf energy balance. The model is forced by incoming shortwave and longwave radiation, as well as near-surface meteorological conditions. The original formulation of MLCan utilized canopy structural traits derived from observations. This project aims to incorporate a plant phenology scheme within MLCan allowing these structural traits to vary dynamically. In the plant phenology scheme implemented here, plant growth is dependent on environmental conditions such as air temperature and soil moisture. The scheme includes functionality that models plant germination, growth, and senescence. These growth stages dictate the variation in six different vegetative carbon pools: storage, leaves, stem, coarse roots, fine roots, and reproductive. The magnitudes of these carbon pools determine land surface parameters such as leaf area index, canopy height, rooting depth and root water uptake capacity. Coupling this phenology scheme with MLCan allows for a more flexible representation of the structure and function of vegetation as it responds to changing environmental conditions.

environmental changes↗
NASA NTRS1982

Ecological optimality in water-limited natural soil-vegetation systems. I - Theory and hypothesis

The solution space of an approximate statistical-dynamic model of the average annual water balance is explored with respect to the hydrologic parameters of both soil and vegetation. Within the accuracy of this model it is shown that water-limited natural vegetation systems are in stable equilibrium with their climatic and pedologic environments when the canopy density and species act to minimize average water demand stress. Theory shows a climatic limit to this equilibrium above which it is hypothesized that ecological pressure is toward maximization of biomass productivity. It is further hypothesized that natural soil-vegetation systems will develop gradually and synergistically, through vegetation-induced changes in soil structure, toward a set of hydraulic soil properties for which the minimum stress canopy density of a given species is maximum in a given climate. Using these hypotheses, only the soil effective porosity need be known to determine the optimum soil and vegetation parameters in a given climate.

Eagleson, P. S.↗
DOE OSTI2022

Multi-century dynamics of the climate and carbon cycle under both high and net negative emissions scenarios

Abstract. Future climate projections from Earth system models (ESMs) typically focus on the timescale of this century. We use a set of five ESMs and one Earth system model of intermediate complexity (EMIC) to explore the dynamics of the Earth's climate and carbon cycles under contrasting emissions trajectories beyond this century to the year 2300. The trajectories include a very-high-emissions, unmitigated fossil-fuel-driven scenario, as well as a mitigation scenario that diverges from the first scenario after 2040 and features an “overshoot”, followed by a decrease in atmospheric CO2 concentrations by means of large net negative CO2 emissions. In both scenarios and for all models considered here, the terrestrial system switches from being a net sink to either a neutral state or a net source of carbon, though for different reasons and centered in different geographic regions, depending on both the model and the scenario. The ocean carbon system remains a sink, albeit weakened by carbon cycle feedbacks, in all models under the high-emissions scenario and switches from sink to source in the overshoot scenario. The global mean temperature anomaly is generally proportional to cumulative carbon emissions, with a deviation from proportionality in the overshoot scenario that is governed by the zero emissions commitment. Additionally, 23rd century warming continues after the cessation of carbon emissions in several models in the high-emissions scenario and in one model in the overshoot scenario. While ocean carbon cycle responses qualitatively agree in both globally integrated and zonal mean dynamics in both scenarios, the land models qualitatively disagree in zonal mean dynamics, in the relative roles of vegetation and soil in driving C fluxes, in the response of the sink to CO2, and in the timing of the sink–source transition, particularly in the high-emissions scenario. The lack of agreement among land models on the mechanisms and geographic patterns of carbon cycle feedbacks, alongside the potential for lagged physical climate dynamics to cause warming long after CO2 concentrations have stabilized, points to the possibility of surprises in the climate system beyond the 21st century time horizon, even under relatively mitigated global warming scenarios, which should be taken into consideration when setting global climate policy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Global fire modelling and control attributions based on the ensemble machine learning and satellite observations

Contemporary fire dynamics is one of the most complex and least understood land surface phenomena. Global fire controls related to climate, vegetation, and anthropogenic activity are usually intertwined, and difficult to disentangle in a quantitative way. Here, we leveraged an ensemble of five machine learning (ML) models and multiple satellite-based observations to conduct global fire modeling for three fire metrics (burned area, fire number, and fire size), and quantified driving mechanisms underlying annual fire changes in a spatially resolved manner for the period 2003–2019. Ensemble learning is a meta-approach that combines multiple ML predictions to improve accuracy, robustness, and generalization performance. We found that the optimized ensemble ML well reproduced annual dynamics of global burned area (R 2 = 0.90, P < 0.001), total fire numbers (R 2 = 0.86, P < 0.001), and averaged fire size (R 2 = 0.70, P < 0.001). Additionally, the ensemble ML captured key spatial patterns of multi-year mean magnitudes, annual variabilities, anomalies, and trends for different fire metrics. Our ML-based fire attributions further highlighted the dominant role of enhanced anthropogenic activity in reducing global burned area (–1.9 Mha/yr, P < 0.01), followed by climate control (–1.3 Mha/yr, P < 0.01) and insignificant positive vegetation control (0.4 Mha/yr, P = 0.60). Spatially, climate dominated a much larger burned area (53.7%) than human (23.4%) or vegetation control (22.9%); however, the counteracting effects from regional wetting and drying trends weakened the net climate impacts on global burned area. The fire number and fire size exhibited similar spatial control patterns with burned area; globally, however, fire number tended to be more affected by climate while fire size more influenced by human activities. Overall, our study confirmed the feasibility and efficiency of ensemble ML in global fire modeling and subsequent control attributions, providing a better understanding of contemporary fire regimes and contributing to robust fire projections in a changing environment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land-surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snow pack evolution now support a range of drought indicators that better reflect coupled water, energy and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Drought↗
DOE OSTI2022

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Studying Rate of Fire Spread Using Landscape Statistics and FIRETEC

Wildfires are multi-system and multi-scale processes that are influenced by the land and atmosphere. Here, we study the influence of land, in terms of topography and vegetation, on uphill and downhill fire spread using data of the state of California as a testbed. The California landscape is first quantified by reclassifying its dominant vegetation types into grass, trees, and shrubs. Next, statistics of slope, vegetation height, and live fuel moisture percentage for each dominant vegetation type are gathered to create ensembles of idealized computational fluid dynamics (CFD) model simulations representative of the landscape. Rate of spread values for each CFD simulation are recorded to observe the differences in perturbing slope and fuel moisture for each vegetation type. We have found that all vegetation types prefer slopes of less than 30 degrees. Trees in the California data are 10 to 40 meters high, whereas shrub heights are less than 3 meters and grass heights are less than 1 meter. The live fuel moisture percentage of both trees and shrubs fluctuates around approximately 100% over one year (2019-2020). Using these parameters as inputs into the CFD model, fire spread rate is found to be fastest in the grass ensemble (maximum 12 m s -1 ), followed by shrubs (maximum 10 m s -1 ) and lastly by trees (maximum 7 m s -1 ). Increasing slope increases the spread rate of the uphill fire but the downhill spread decreases for slopes greater than 15 degrees. Increasing fuel moisture decreases the spread rate both up and down the hill. This study highlights the impacts that topography and vegetation have on fire spread.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snowpack evolution now support a range of drought indicators that better reflect coupled water, energy, and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity, and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Christa D Peters-lidard↗
NASA NTRSDate not supplied

Advances in Land Surface Models and Indicators for Drought Monitoring and Prediction

Millions of people across the globe are affected by droughts every year, and recent droughts have highlighted the considerable agricultural impacts and economic costs of these events. Monitoring the state of droughts depends on integrating multiple indicators that each capture particular aspects of hydrologic impact and various types and phases of drought. As the capabilities of land surface models and remote sensing have improved, important physical processes such as dynamic, interactive vegetation phenology, groundwater, and snowpack evolution now support a range of drought indicators that better reflect coupled water, energy, and carbon cycle processes. In this work, we discuss these advances, including newer classes of indicators that can be applied to improve the characterization of drought onset, severity, and duration. We utilize a new model-based drought reconstruction to illustrate the role of dynamic phenology and groundwater in drought assessment. Further, through case studies on flash droughts, snow droughts, and drought recovery, we illustrate the potential advantages of advanced model physics and observational capabilities, especially from remote sensing, in characterizing droughts.

Christa D Peters-lidard↗
DOE OSTI2024

Simulations of flaming combustion and flaming-to-smoldering transition in wildland fire spread at flame scale

Our objective in the present study is to provide basic insights into the coupling between external-gas and solid biomass vegetation processes that control the dynamics of flame spread in wildland fire problems. Here we focus on a modeling approach that resolves processes occurring at vegetation and flame scales, i.e., the formation of flammable vapors due to the thermal degradation of the solid biomass, the subsequent combustion in ambient air, the thermal feedback to the biomass through radiative and convective heat transfer, and the possible transition from flaming combustion (taking place outside of the solid biomass) to smoldering combustion (taking place inside the solid biomass). The capability uses a multiphase combustion framework and treats external-gas processes through a Large Eddy Simulation solver and solid biomass processes through a discrete particle model. The discrete particle model adopts a one-dimensional porous medium formulation, includes descriptions of drying, thermal pyrolysis, oxidative pyrolysis, and char oxidation, as well as a description of the external-gas-to-solid-biomass diffusion of oxygen mass; the discrete particle model thereby provides a treatment of in-depth oxidative processes and allows the simulation of smoldering combustion. The modeling capability is applied to the simulation of fire spread across a surrogate biomass vegetation bed corresponding to a discrete array of cylindrical-shaped, vertically-oriented, pine wood sticks, characterized by a monomodal size distribution, in horizontal flat terrain and under wind-aided conditions. The numerical results demonstrate that the model can simulate successful flaming-to-smoldering transition followed by complete biomass consumption.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

Controls and relationships of soil organic carbon abundance and persistence vary across pedo–climatic regions

One of the largest uncertainties in the terrestrial carbon cycle is the timing and magnitude of soil organic carbon (SOC) response to climate and vegetation change. This uncertainty prevents models from adequately capturing SOC dynamics and challenges the assessment of management and climate change effects on soils. Reducing these uncertainties requires simultaneous investigation of factors controlling the amount (SOC abundance) and duration (SOC persistence) of stored C. We present a global synthesis of SOC and radiocarbon profiles (n Profile = 597) to assess the timescales of SOC storage. We use a combination of statistical and depth–resolved compartment models to explore key factors controlling the relationships between SOC abundance and persistence across pedo–climatic regions and with soil depth. This allows us to better understand (i) how SOC abundance and persistence covary across pedo–climatic regions and (ii) how the depth dependence of SOC dynamics relates to climatic and mineralogical controls on SOC abundance and persistence. We show that SOC abundance and persistence are differently related; the controls on these relationships differ substantially between major pedo–climatic regions and soil depth. For example, large amounts of persistent SOC can reflect climatic constraints on soils (e.g., in tundra/polar regions) or mineral absorption, reflected in slower decomposition and vertical transport rates. In contrast, lower SOC abundance can be found with lower SOC persistence (e.g., in highly weathered tropical soils) or higher SOC persistence (e.g., in drier and less productive regions). We relate variable patterns of SOC abundance and persistence to differences in the processes constraining plant C input, microbial decomposition, vertical C transport and mineral SOC stabilization potential. This process–oriented grouping of SOC abundance and persistence provides a valuable benchmark for global C models, highlighting that pedo–climatic boundary conditions are crucial for predicting the effects of climate change and soil management on future C abundance and persistence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Groundwater and Surface Water Flow (GSFLOW) model files for the East River, Colorado

The data package contains model input files and executables for the East River, Colorado (750 km2) located in the headwaters of the Upper Colorado River Basin. The code applied is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution to represent topographic complexity and a daily timestep accounts for energy and water partitioning between the snowpack, vegetation, soil zone and bedrock. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model accounts for nine stratigraphic units and extends 400 m below land surface. Using this modeling framework we explore historical water budgets (water year 1987-2022) and the influence of seasonal warming and associated mechanisms driving groundwater declines and streamflow loss. A Readme_062424.txt file provides instructions on how to download all files. Input and output files are provided for the historical simulation representative of water years 1987 to 2022 (i.e. baseline) and the all-year +4C warming scenario. Instructions are provided to run the seasonal warming scenarios with warming applied only to the autumn, winter, spring or summer months. Modeled output used in figures for Carroll et al., 2024 are also provided with metadata describing where these data were obtained. This dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.Updated on 06-24-2024: The dataset was updated to include model outputs (updated base.zip and 4C.zip files), figure source files (Figure_Source_Files.zip), an updated Readme file (Readme_062424.txt), climate input files (climate_input_files.zip), and the adoption of ESS-DIVE File Level Metadata and CSV reporting formats (inclusion of flmd.csv and dd.csv files).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial patterns of snow distribution in the sub-Arctic

Abstract. The spatial distribution of snow plays a vital role in sub-Arctic and Arctic climate, hydrology, and ecology due to its fundamental influence on the water balance, thermal regimes, vegetation, and carbon flux. However, the spatial distribution of snow is not well understood, and therefore, it is not well modeled, which can lead to substantial uncertainties in snow cover representations. To capture key hydro-ecological controls on snow spatial distribution, we carried out intensive field studies over multiple years for two small (2017–2019; ∼ 2.5 km2) sub-Arctic study sites located on the Seward Peninsula of Alaska. Using an intensive suite of field observations (> 22 000 data points), we developed simple models of the spatial distribution of snow water equivalent (SWE) using factors such as topographic characteristics, vegetation characteristics based on greenness (normalized different vegetation index, NDVI), and a simple metric for approximating winds. The most successful model was random forest, using both study sites and all years, which was able to accurately capture the complexity and variability of snow characteristics across the sites. Approximately 86 % of the SWE distribution could be accounted for, on average, by the random forest model at the study sites. Factors that impacted year-to-year snow distribution included NDVI, elevation, and a metric to represent coarse microtopography (topographic position index, TPI), while slope, wind, and fine microtopography factors were less important. The characterization of the SWE spatial distribution patterns will be used to validate and improve snow distribution modeling in the Department of Energy's Earth system model and for improved understanding of hydrology, topography, and vegetation dynamics in the sub-Arctic and Arctic regions of the globe.

54 ENVIRONMENTAL SCIENCES↗