Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling vegetation change”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

DeepPhenoMem V1.0: deep learning modelling of canopy greenness dynamics accounting for multi-variate meteorological memory effects on vegetation phenology

Abstract. Vegetation phenology plays a key role in controlling the seasonality of ecosystem processes that modulate carbon, water and energy fluxes between the biosphere and atmosphere. Accurate modelling of vegetation phenology in the interplay of Earth's surface and the atmosphere is thus crucial to understand how the coupled system will respond to and shape climatic changes. Phenology is controlled by meteorological conditions at different timescales: on the one hand, changes in key meteorological variables (temperature, water, radiation) can have immediate effects on the vegetation development; on the other hand, phenological changes can be driven by past environmental conditions, known as memory effects. However, the processes governing meteorological memory effects on phenology are not completely understood, resulting in their limited performance of vegetation phenology represented in land surface models. A deep learning model, specifically a long short-term memory network (LSTM), has the potential to capture and model the meteorological memory effects on vegetation phenology. Here, we apply the LSTM to model the vegetation phenology using meteorological drivers and high-temporal-resolution canopy greenness observations through digital repeat photography by the PhenoCam network. We compare a multiple linear regression model, a no-memory-effect LSTM model and a full-memory-effect LSTM model to predict the whole seasonal greenness trajectory and the corresponding phenological transition dates across 50 sites and 317 site years during 2009–2018, covering deciduous broadleaf forests, evergreen needleleaf forests and grasslands. Results show that the deep learning model outperforms the multiple linear regression model, and the full-memory-effect LSTM model performs better than the no-memory-effect model for all three plant function types (median R2 of 0.878, 0.957 and 0.955 for broadleaf forests, evergreen needleleaf forests and grasslands). We also find that the full-memory-effect LSTM model is capable of predicting the seasonal dynamic variations of canopy greenness and reproducing trends in shifting phenological transition dates. We also performed a sensitivity analysis of the full-memory-effect LSTM model to assess its plausibility, revealing its coherence with established knowledge of vegetation phenology sensitivity to meteorological conditions, particularly changes in temperature. Our study highlights that (1) multi-variate meteorological memory effects play a crucial role in vegetation phenology, and (2) deep learning opens up new avenues for improving the representation of vegetation phenological processes in land surface models via a hybrid modelling approach.

Geology↗

Multispectral UAV imagery of experimental freshwater wetlands under 5 ppt saltwater intrusion, Louisiana, 2023 and 2024

Multispectral imagery was collected using an unmanned aerial vehicle (UAV) to evaluate how freshwater vegetation responds to short-term simulated saltwater intrusion events. The purpose of this data collection was to understand how plant health changes in response to acute salinity exposure, which is increasingly relevant in coastal wetland ecosystems facing sea level rise and storm surge events, such as in coastal Louisiana. Three experimental saltwater intrusions were conducted at a salinity of approximately 5 parts per thousand (ppt) for durations of 6-days, 10-days, and 17-days. UAV flights occurred both before and after each treatment. The resulting imagery was processed using Pix4DMapper software to georeference the images and generate orthomosaics. The multispectral sensor used in this study captures reflectance in five bands: blue, green, red, red-edge, and near-infrared. The uploaded data consist of georeferenced .tif orthomosaics for each spectral band, which are compatible with GIS software for vegetation analysis. This imagery can be utilized in investigations into vegetation stress, remote sensing of freshwater wetland ecosystems, and modeling of plant response to environmental changes.

EARTH SCIENCE > BIOSPHERE > ECOSYSTEMS↗

Comparing deuterium excess to large-scale precipitation recycling models in the tropics

Precipitation recycling is essential to sustaining regional ecosystems and water supplies, and it is impacted by land development and climate change. This is especially true in the tropics, where dense vegetation greatly influences recycling. Unfortunately, large-scale models of recycling often exhibit high uncertainty, complicating efforts to estimate recycling. Here, we examine how deuterium excess (d-excess), a stable-isotope quantity sensitive to recycling effects, acts as an observational proxy for recycling. While past studies have connected variability in d-excess to precipitation origins at local or regional scales, our study leverages >3000 precipitation isotope samples to quantitatively compare d-excess against three contemporary recycling models across the global tropics. Using rank-correlation, we find statistically significant agreement ( $\bar \tau = 0.52$ τ ¯ = 0.52 to $0.70$ 0.70 ) between tropical d-excess and recycling that is strongly mediated by seasonal precipitation, vegetation density, and scale mismatch. Our results detail the complex relationship between d-excess and precipitation recycling, suggesting avenues for further investigation.

54 ENVIRONMENTAL SCIENCES↗

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗

Simulating environmentally-sensitive tree recruitment in vegetation demographic models

Vegetation demographic models (VDMs) endeavor to predict how global forests will respond to climate change. This requires simulating which trees, if any, are able to recruit under changing environmental conditions. We present a new recruitment scheme for VDMs in which functional-type-specific recruitment rates are sensitive to light, soil moisture and the productivity of reproductive trees. We evaluate the scheme by predicting tree recruitment for four tropical tree functional types under varying meteorology and canopy structure at Barro Colorado Island, Panama. We compare predictions to those of a current VDM, quantitative observations and ecological expectations. We find that the scheme improves the magnitude and rank order of recruitment rates among functional types and captures recruitment limitations in response to variable understory light, soil moisture and precipitation regimes. Furthermore, our results indicate that adopting this framework will improve VDM capacity to predict functional-type-specific tree recruitment in response to climate change, thereby improving predictions of future forest distribution, composition and function.

59 BASIC BIOLOGICAL SCIENCES↗

Global Carbon Budget 2021

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize datasets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the first time, an approach is shown to reconcile the difference in our ELUC estimate with the one from national greenhouse gas inventories, supporting the assessment of collective countries' climate progress. For the year 2020, EFOS declined by 5.4 % relative to 2019, with fossil emissions at 9.5 ± 0.5 GtC yr−1 (9.3 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 0.9 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission of 10.2 ± 0.8 GtC yr−1 (37.4 ± 2.9 GtCO2). Also, for 2020, GATM was 5.0 ± 0.2 GtC yr−1 (2.4 ± 0.1 ppm yr−1), SOCEAN was 3.0 ± 0.4 GtC yr−1, and SLAND was 2.9 ± 1 GtC yr−1, with a BIM of −0.8 GtC yr−1. The global atmospheric CO2 concentration averaged over 2020 reached 412.45 ± 0.1 ppm. Preliminary data for 2021 suggest a rebound in EFOS relative to 2020 of +4.8 % (4.2 % to 5.4 %) globally. Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2020, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use changes emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and datasets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this dataset (Friedlingstein et al., 2020, 2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014, 2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2021 (Friedlingstein et al., 2021).

Friedlingstein, Pierre (ORCID:0000000333094739)↗

Global Carbon Budget 2022

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere in a changing climate is critical to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe and synthesize data sets and methodologies to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFOS) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly, and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) is estimated with global ocean biogeochemistry models and observation-based data products. The terrestrial CO2 sink (SLAND) is estimated with dynamic global vegetation models. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, with fossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission (including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1 (40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with a BIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low or sinks were too high). The global atmospheric CO2 concentration averaged over 2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest an increase in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %) globally and atmospheric CO2 concentration reaching 417.2 ppm, more than 50 % above pre-industrial levels (around 278 ppm). Overall, the mean and trend in the components of the global carbon budget are consistently estimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadal variability in CO2 fluxes. Comparison of estimates from multiple approaches and observations shows (1) a persistent large uncertainty in the estimate of land-use change emissions, (2) a low agreement between the different methods on the magnitude of the land CO2 flux in the northern extratropics, and (3) a discrepancy between the different methods on the strength of the ocean sink over the last decade. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set. The data presented in this work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

54 ENVIRONMENTAL SCIENCES↗

Utility of near‐surface phenology in estimating productivity and evapotranspiration across diverse ecosystems

Abstract Agroecosystems, which include row crops, pasture, and grass and shrub grazing lands, are sensitive to changes in management, weather, and genetics. To better understand how these systems are responding to changes, we need to improve monitoring and modeling carbon and water dynamics. Vegetation Indices (VIs) are commonly used to estimate gross primary productivity (GPP) and evapotranspiration (ET), but these empirical relationships are often location and crop specific. There is a need to evaluate if VIs can be effective and, more general, predictors of ecosystem processes through time and across different agroecosystems. Near‐surface photographic (red‐green‐blue) images from PhenoCam can be used to calculate the VI green chromatic coordinate (G CC ) and offer a pathway to improve understanding of field‐scale relationships between VIs and GPP and ET. We synthesized observations spanning 76 site‐years across 15 agroecosystem sites with PhenoCam G CC and GPP or ET estimates from eddy covariance (EC) to quantify interannual variability (IAV) in the relationship between GPP and ET and G CC across. We uncovered a high degree of variability in the strength and slopes of the G CC ∼ GPP and ET relationships (R 2 = 0.1 ‐ 0.9) within and across production systems. Overall, G CC is a better predictor of GPP than ET (R 2 = 0.64 and 0.54, respectively), performing best in croplands (R 2 = 0.91). Shrub‐dominated systems exhibit the lowest predictive power of G CC for GPP and ET but have less IAV in slope. We propose that PhenoCam estimates of G CC could provide an alternative approach for predictions of ecosystem processes.

Environmental Sciences & Ecology↗

Attributing human mortality from fire PM 2.5 to climate change

Climate change intensifies fire smoke, emitting hazardous air pollutants that impact human health. However, the global influence of climate change on fire-induced health impacts remains unquantified. Here, in this study, we used three well-tested fire-vegetation models in combination with a chemical transport model and health risk assessment framework to attribute global human mortality from fire fine particulate matter (PM 2.5 ) emissions to climate change. Of the 46401 (1960s) –to 98748 (2010s) annual fire PM 2.5 mortalities, 669 (1.2%, 1960s) –to 12566 (12.8%, 2010s) were attributed to climate change. The most substantial influence of climate change on fire mortality occurred in South America, Australia, and Europe, coinciding with decreased relative humidity, and in boreal forests with increased air temperature. Increasing relative humidity lowered fire mortality in other regions, like South Asia. Our study highlights the role of climate change in fire mortality, aiding public health authorities in spatial targeting adaptation measures for sensitive fire-prone areas.

54 ENVIRONMENTAL SCIENCES↗

How deep should we go to understand roots at the top of the world?

Harsh environmental conditions and the short summers of northern, high-latitude biomes impose unique constraints on the plants that live in the arctic tundra and the boreal forest. To escape the harsh aboveground environment, plants in these habitats often allocate a large portion of their biomass belowground to facilitate nutrient acquisition. In turn, the proximity of living plant roots to vast stores of sequestered soil carbon in these biomes means that shifts in rooting depth distribution and the size of the root–soil interface could significantly contribute to ongoing climate change. Indeed, plant ‘priming’ of rhizosphere decomposition via root exudation, particularly from shallowly distributed roots, can lead to losses of carbon from tundra soils. While we have a hard-won understanding of the distribution of plant communities across the arctic tundra and the boreal forest from direct field observations scaled to the landscape level using climate-informed mapping techniques (i.e. the Circumpolar Arctic Vegetation Map (CAVM); Walker et al. , 2005), these vegetation maps are only the tip of the iceberg (Iversen et al. , 2015). Root form and function remain hidden beneath the land surface. In an article recently published in New Phytologist , Blume-Werry et al. (2023, 10.1111/nph.18998) asked whether rooting depth distribution, and ensuing carbon emissions, could be inferred from commonly used vegetation mapping classifications across the pan-Arctic. An important question to guide our understanding, mapping, and prediction of belowground characteristics and ecosystem feedbacks at the top of the world. Unfortunately, they found that the answer was ‘not quite’. While rooting depth distribution varied demonstrably, in turn causing substantial changes in modeled carbon emissions via rhizosphere ‘priming’, variation across rooting depth profiles did not correspond with vegetation mapping classes. If we are unable to predict belowground rooting depth distributions across large spatial scales by leveraging aboveground vegetation community distributions, how then should belowground researchers proceed?

59 BASIC BIOLOGICAL SCIENCES↗

Consistent time allocation fraction to vegetation green-up versus senescence across northern ecosystems despite recent climate change

Extended growing season lengths under climatic warming suggest increased time for plant growth. However, research has focused on climatic impacts to the timing or duration of distinct phenological events. Comparatively little is known about impacts to the relative time allocation to distinct phenological events, for example, the proportion of time dedicated to leaf growth versus senescence. We use multiple satellite and ground-based observations to show that, despite recent climate change during 2001 to 2020, the ratio of time allocated to vegetation green-up over senescence has remained stable [1.27 (± 0.92)] across more than 83% of northern ecosystems. This stability is independent of changes in growing season lengths and is caused by widespread positive relationships among vegetation phenological events; longer vegetation green-up results in longer vegetation senescence. These empirical observations were also partly reproduced by 13 dynamic global vegetation models. Our work demonstrates an intrinsic biotic control to vegetation phenology that could explain the timing of vegetation senescence under climate change.

54 ENVIRONMENTAL SCIENCES↗

A Continental-Scale Estimate of Soil Organic Carbon Change at NEON Sites and Their Environmental and Edaphic Controls

Current carbon cycle models focus on the effects of climate and land-use change on primary productivity and microbial-mineral dependent carbon turnover in the topsoil, while less attention has been paid to vertical soil processes and soil-dependent response to land-use change along the profile. In this study, a spatial-temporal analysis was used to estimate soil organic carbon (SOC) change in topsoil/A horizon and subsoil/B horizon at National Ecological Observatory Network (NEON) sites, USA over 30 years. To separate the effects of land-use, environmental, and edaphic factors on SOC change, space-for-time substitution was used in combination with the Continuous Change Detection and Classification algorithm and Structural Equation Modeling. Results showed that (a) under natural vegetation, Spodosols and Inceptisols found in the eastern NEON sites had substantial topsoil SOC accumulation (+0.4 to +1.2 Mg C ha –1 year –1 ), while Inceptisols and Andisols in the west had a comparable magnitude of topsoil SOC loss (–0.5 to –1.8 Mg C ha –1 year –1 ); (b) Mollisols and Alfisols in the Central Plains sites were susceptible to significant SOC loss under farming and grazing; (c) Runoff/erosion and leaching potential, vertical translocation, and mineral sorption were the most important factors controlling SOC variation across the NEON sites. Our work could be used to parameterize ecosystem models simulating SOC change.

54 ENVIRONMENTAL SCIENCES↗

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗

IAM-FIRE: a Climate Emulator–Based Framework to Project Wildfire Impacts and Risks for Integrated Assessment Models

Most Integrated Assessment Models (IAMs) underrepresent dynamic feedbacks from climate-driven disturbances such as wildfires, potentially overestimating the permanence of land-based carbon sinks. In particular, representing the impacts of forest fires is becoming increasingly important, as these are expected to intensify in the coming years. We introduce IAM-FIRE (Integrated Assessment Model – Fire Impacts & Risks Emulator), a novel framework that enables the projection of wildfire burned area (BA) and carbon emissions (CE) directly from IAM outputs. IAM-FIRE combines a spatial climate emulator, land-use downscaling, vegetation productivity modelling, and an empirical fire model to generate global annual wildfire impacts for arbitrary socioeconomic and emissions scenarios at 0.5° resolution for the period 2020–2100. Calibrated against GFEDv5 observations and using inputs from the Global Change Analysis Model (GCAM), we report projections BA and CE derived from IAM-FIRE for four scenarios: SSP1-2.6, SSP2-4.5, SSP3-6.6 and SSP5-7.6. The model reproduces historical global trends for total BA, including the observed global decline since the early 2000s, and for forest BA. Projected fire trajectories differ strongly among scenarios: total BA range from declines under SSP1-2.6 (-3.36 Mha yr-1) to increases under SSP3-6.6 (+1.6 Mha yr-1). Corresponding total CE show a similar divergence ranging from -15 to +10.6 TgC yr-1. Socioeconomic development exerts a dominant suppressing effect on wildfire impacts while climate change and CO2-driven increases in vegetation productivity amplify fire risk, particularly under high-emissions pathways. Compared with CMIP6 fire-enabled Earth System Models, IAM-FIRE exhibits greater sensitivity to radiative forcing and a stronger role for human-driven fire suppression, highlighting substantial structural uncertainties in future fire projections. By providing a computationally efficient and internally consistent approach to represent wildfire impacts within IAMs, IAM-FIRE enables systematic exploration of fire–climate–land feedbacks and supports improved assessments of mitigation permanence and climate risks in future integrated scenarios.

Rouhette, Theo↗

Semi-automated Design of Artificial Intelligence Earth Science Models

Prediction and observation of water cycles involve not only patterns isolated in space and time, but rather modeling complex spatio-temporal relationships across multiple sources of data and domains. For instance, Evapotranspiration (ET) and Leaf Area Indexes (LAI) are two critical components in DOE’s Energy Exascale Earth System Model (E3SM). Accurate assessments of ET and LAI are critical for understanding hydrological processes, deforestation, crop yield, and irrigation impacts. However, current ET estimates for global simulations are available at very coarse spatial resolution. They are usually derived from satellite data based on broad plant functional types (PFT), which fail to capture the fine-scale variations due to change in vegetation type across the globe. Within this context and in light of the data-model integration challenges highlighted in the EESSD Strategic Plan, the new era of AI model development for geosciences calls for data-driven methods that provide domain scientists with estimations of parameters such as PFT and LAI in an efficient, interpretable, and easy-to-operate manner.

Ambrozio Dias, Philipe↗

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Shrub Expansion Simulations at Trail Valley Creek Tundra site using E3SM Land Model (ELM) Arctic-focused Version

The warming of the Arctic is causing substantial compositional, structural, and functional changes in tundra vegetation including shrub and densification in parts of the Arctic. Assessing the impact of these changes in vegetation composition on the Arctic’s carbon and energy budgets is important to constrain projected local and global surface-atmosphere exchanges. We conduct a sensitivity analysis of the projected surface energy fluxes, soil carbon pools, and carbon dioxide fluxes (net ecosystem exchange, gross primary production, and ecosystem respiration) between present day and 2100 to different shrub expansion rates and air temperature increases under future emission scenarios (intermediate – RCP4.5, and high – RCP8.5) using the Arctic-focused version of the Energy Exascale Earth System Model (E3SM) Land Model (ELM). We focus on Trail Valley Creek (TVC), a mineral upland tundra site located in the western Canadian Arctic, which is experiencing tall shrub densification and expansion. In this study, we run TVC under two different warming scenarios RCP4.5 and RCP 8.5 and simulate different shrubification rates projected until year 2100. In this repository, we include all the forcing, input, parameters, and output data corresponding to all the simulations performed. flmd.csv includes a detailed description of the datasets files.

54 ENVIRONMENTAL SCIENCES↗

Future Hurricanes Will Increase Palm Abundance and Decrease Aboveground Biomass in a Tropical Forest

Hurricanes are expected to intensify throughout the 21st century, yet the impact of frequent major hurricanes on tropical ecosystems remains unknown. To investigate tropical forest damage and recovery under different hurricane regimes, we generate a suite of scenarios based on Coupled Model Intercomparison Project Phase 6 climate projections and increased hurricane recurrence and intensity for the Luquillo Experimental Forest, Puerto Rico. We then use the Ecosystem Demography model to predict changes in carbon stocks, forest structure and composition. In this work, our results indicate that frequent hurricane disturbances in the future would decrease the overall aboveground biomass, decrease the dominance of late-successional species, but increase the dominance of palm species. Warmer climates with increased CO 2 would have little effect on the functional-type composition but increase the aboveground biomass. However, the predicted climate and CO 2 fertilization effects would not compensate for the biomass loss due to more frequent severe-hurricane disturbances.

54 ENVIRONMENTAL SCIENCES↗