Search NASA⌕ Search

SEARCH · Search NASA

Results for “GPP”

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 127 records · Page 7

Contribution of Phenological and Physiological Variations on Northern Vegetation Productivity Changes over Last Three Decades

Plant phenology and maximum photosynthetic state determine spatiotemporal variability of gross primary productivity (GPP) of vegetation. Recent warming induced impacts accelerate shifts of phenology and physiological status over Northern vegetated land. Thus, understanding and quantifying these changes are very important. Here, we investigate 1) how vegetation phenology and physiological status (maximum photosynthesis) are evolved over last three decades and 2) how such components (phenology and physiological status) contribute on inter-annual variation of the GPP during the last three decades. We utilized both long-term remotely sensed (GIMMS (Global Inventory Modeling and Mapping Studies), NDVI3g (Normalized Difference Vegetation Index 3rd generation) and MODIS (Moderate Resolution Imaging Spectroradiometer)) to extract larger scale phenology metrics (growing season start, end and duration); and productivity (i.e., growing season integrated vegetation index, GSIVI) to answer these questions. For evaluation purpose, we also introduced field-measured phenology and productivity datasets (e.g., FLUXNET) and possible remotely-sensed and modeled metrics at continental and regional scales. From this investigation, we found that onset of the growing season has advanced by 1.61 days per decade and the growing season end has delayed by 0.67 days per decade over the circumpolar region. This asymmetric extension of growing season results in a longer growing-season trend (2.96 days per decade) and widespread increasing vegetation-productivity trend (2.96 GSIVI per decade) over Northern land. However, the regionally-diverged phenology shift and maximum photosynthetic state contribute differently characterized productivity, inter-annual variability and trend. We quantified that about 50 percent, 13 percent and 6.5 percent of Northern land's inter-annual variability are dominantly controlled by the onset of the growing season, the end of the growing season and the maximum photosynthetic state, respectively. Productivity characterization over the other approximately 30 percent region has been driven by these co-dominant drivers. Our study clearly shows that regionally different contribution of phenological and physiological components on characterizing vegetation production over the last three decades.

productivity↗

Combining Livestock Production Information in a Process-Based Vegetation Model to Reconstruct the History of Grassland Management

Grassland management type (grazed or mown) and intensity (intensive or extensive) play a crucial role in the greenhouse gas balance and surface energy budget of this biome, both at field scale and at large spatial scale. However, global gridded historical information on grassland management intensity is not available. Combining modelled grass-biomass productivity with statistics of the grass-biomass demand by livestock, we reconstruct gridded maps of grassland management intensity from 1901 to 2012. These maps include the minimum area of managed vs. maximum area of unmanaged grasslands and the fraction of mown vs. grazed area at a resolution of 0.5deg by 0.5deg. The grass-biomass demand is derived from a livestock dataset for 2000, extended to cover the period 19012012. The grass-biomass supply (i.e. forage grass from mown grassland and biomass grazed) is simulated by the process-based model ORCHIDEE-GM driven by historical climate change, risingCO2 concentration, and changes in nitrogen fertilization. The global area of managed grassland obtained in this study increases from 6.1 x 10(exp 6) km(exp 2) in 1901 to 12.3 x 10(exp 6) kmI(exp 2) in 2000, although the expansion pathway varies between different regions. ORCHIDEE-GM also simulated augmentation in global mean productivity and herbage-use efficiency over managed grassland during the 20th century, indicating a general intensification of grassland management at global scale but with regional differences. The gridded grassland management intensity maps are model dependent because they depend on modelled productivity. Thus specific attention was given to the evaluation of modelled productivity against a series of observations from site-level net primary productivity (NPP) measurements to two global satellite products of gross primary productivity (GPP) (MODIS-GPP and SIF data). Generally, ORCHIDEE-GM captures the spatial pattern, seasonal cycle, and inter-annual variability of grassland productivity at global scale well and thus is appropriate for global applications presented here.

Livestock dataset↗

Soil Water Content Effects on Net Ecosystem CO2 Exchange and Actual Evapotranspiration in a Mediterranean Semiarid Savanna of Central Chile

Biosphere-atmosphere water and carbon fluxes depend on ecosystem structure, and their magnitudes and seasonal behavior are driven by environmental and biological factors. We studied the seasonal behavior of net ecosystem CO2 exchange (NEE), Gross Primary Productivity (GPP), Ecosystem Respiration (RE), and actual evapotranspiration (ETa) obtained by eddy covariance measurements during two years in a Mediterranean Acacia savanna ecosystem (Acacia caven) in Central Chile. The annual carbon balance was −53 g C/sq. m in 2011 and −111 g C/sq. m in 2012, showing that the ecosystem acts as a net sink of CO2, notwithstanding water limitations on photosynthesis observed in this particularly dry period. Total annual ETa was of 128 mm in 2011 and 139 mm in 2012. Both NEE and ETa exhibited strong seasonality with peak values recorded in the winter season (July to September), as a result of ecosystem phenology, soil water content and rainfall occurrence. Consequently, the maximum carbon assimilation rate occurred in wintertime. Results show that soil water content is a major driver of GPP and RE, defining their seasonal patterns and the annual carbon assimilation capacity of the ecosystem, and also modulating the effect that solar radiation and air temperature have on NEE components at shorter time scales.

Meza, Francisco J.↗

Identifying and Forecasting Potential Biophysical Risk Areas within a Tropical Mangrove Ecosystem Using Multi-Sensor Data

Mangroves are one of the most productive ecosystems known for provisioning of various ecosystem goods and services. They help in sequestering large amounts of carbon, protecting coastline against erosion, and reducing impacts of natural disasters such as hurricanes. Bhitarkanika Wildlife Sanctuary in Odisha harbors the second largest mangrove ecosystem in India. This study used Terra, Landsat and Sentinel-1 satellite data for spatio-temporal monitoring of mangrove forest within Bhitarkanika Wildlife Sanctuary between 2000 and 2016. Three biophysical parameters were used to assess mangrove ecosystem health: leaf chlorophyll (CHL), Leaf Area Index (LAI), and Gross Primary Productivity (GPP). A long-term analysis of meteorological data such as precipitation and temperature was performed to determine an association between these parameters and mangrove biophysical characteristics. The correlation between meteorological parameters and mangrove biophysical characteristics enabled forecasting of mangrove health and productivity for year 2050 by incorporating IPCC projected climate data. A historical analysis of land cover maps was also performed using Landsat 5 and 8 data to determine changes in mangrove area estimates in years 1995, 2004 and 2017. There was a decrease in dense mangrove extent with an increase in open mangroves and agricultural area. Despite conservation efforts, the current extent of dense mangrove is projected to decrease up to 10% by the year 2050. All three biophysical characteristics including GPP, LAI and CHL, are projected to experience a net decrease of 7.7%, 20.83% and 25.96% respectively by 2050 compared to the mean annual value in 2016. This study will help the Forest Department, Government of Odisha in managing and taking appropriate decisions for conserving and sustaining the remaining mangrove forest under the changing climate and developmental activities.

Shanti Shrestha↗

Modeling Phenological Controls on Carbon Dynamics in Dryland Sagebrush Ecosystems

Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbonfluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle,we must also be able to understand and model its effects on dryland vegetation. Dynamic Global VegetationModels (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduceseasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to bothtotal productivity and competitive interactions with other plants, errors in how process-based models representphenology hinder our ability to predict climate change impacts. This may be particularly problematic in drylandecosystems where many species have developed a complex phenology in response to seasonal variability in bothmoisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure ofexisting phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) andgross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized modelparameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, wemodified the model to include a semi-deciduous phenology type and added flexibility to the representation ofgrass phenology. While optimizing parameters reduced model bias, the largest gains in model performance wereassociated with the development of our new representation of phenology. This modified model was able to bettercapture seasonal patterns of both leaf area index (R2=0.75) and gross primary productivity (R2=0.84), thoughits ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new modelalso resulted in a more realistic outcome of modeled competition between grass and shrubs. These findingsdemonstrate the importance of improving how DGVMs represent phenology in order to accurately forecastclimate change impacts in dryland ecosystems.

Ecosystem model Phenology Parameter optimization E↗

Modeling Phenological Controls on Carbon Dynamics in Dryland Sagebrush Ecosystems

Dryland ecosystems play an important role in determining how precipitation anomalies affect terrestrial carbon fluxes at regional to global scales. Thus, to understand how climate change may affect the global carbon cycle, we must also be able to understand and model its effects on dryland vegetation. Dynamic Global Vegetation Models (DGVMs) are an important tool for modeling ecosystem dynamics, but they often struggle to reproduce seasonal patterns of plant productivity. Because the phenological niche of many plant species is linked to both total productivity and competitive interactions with other plants, errors in how process-based models represent phenology hinder our ability to predict climate change impacts. This may be particularly problematic in dryland ecosystems where many species have developed a complex phenology in response to seasonal variability in both moisture and temperature. Here, we examine how uncertainty in key parameters as well as the structure of existing phenology routines affect the ability of a DGVM to match seasonal patterns of leaf area index (LAI) and gross primary productivity (GPP) across a temperature and precipitation gradient. First, we optimized model parameters using a combination of site-level eddy covariance data and remotely-sensed LAI data. Second, we modified the model to include a semi-deciduous phenology type and added flexibility to the representation of grass phenology. While optimizing parameters reduced model bias, the largest gains in model performance were associated with the development of our new representation of phenology. This modified model was able to better capture seasonal patterns of both leaf area index (R(exp 2) = 0.75) and gross primary productivity (R(exp 2) = 0.84), though its ability to estimate total annual GPP depended on using eddy covariance data for optimization. The new model also resulted in a more realistic outcome of modeled competition between grass and shrubs. These findings demonstrate the importance of improving how DGVMs represent phenology in order to accurately forecast climate change impacts in dryland ecosystems.

Ecosystem model↗

Larch (Larix dahurica Turcz) growth response to climate change in the Siberian permafrost zone

Larch-dominant communities are the most extensive high-latitude forests in Eurasia and are experiencing the strongest impacts from warming temperatures. We analyzed larch (Larix dahurica Turcz) growth index (GI) response to climate change. The studied larch-dominant communities are located within the permafrost zone of Northern Siberia at the northern tree limit (ca. N 67° 38′, E 99° 07′). Methods included dendrochronology, analysis of climate variables, root zone moisture content, and satellitederived gross (GPP) and net (NPP) primary productivity. It was found that larch response to warming included a period of increased annual growth increment (GI) (from the 1970s to ca. 1995) with a follow on GI decline. Increase in GI correlated with summer air temperature, whereas an observed decrease in GI was caused by water stress (vapor pressure deficit and drought increase). Water stress impact on larch growth in permafrost was not observed before the onset of warming (ca. 1970). Water limitation was also indicated by GI dependence on soil moisture stored during the previous year. Water stress was especially pronounced for stands growing on rocky soils with low water-holding capacity. GPP of larch communities showed an increasing trend, whereas NPP stagnated. A similar pattern of GI response to climate warming has also been observed for Larix sibirica Ledeb, Pinus sibirica Du Tour, and Abies sibirica Ledeb in the forests of southern Siberia. Thus, warming in northern Siberia permafrost zone resulted in an initial increase in larch growth from the 1970s to the mid-1990s. After that time, larch growth increment has decreased. Since ca. 1990, water stress at the beginning of the vegetative period became, along with air temperature, a main factor affecting larch growth within the permafrost zone.

Larix dahurica↗

Impacts of Degradation on Water, Energy, and Carbon Cycling of the Amazon Tropical Forests

Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high‐resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED‐2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED‐2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.

Marcos Longo↗

A Multi-scaled Analysis of Forest Structure using Individual-Based Modeling in a Costa Rican Rainforest

Consideration of scale is essential when examining structural relationships in forests. In this study, we present a parameterization of the FORMIND individual-based forest model for old growth Atlantic lowland rainforest in La Selva, Costa Rica. Results show that the simulated forest reproduces the structural complexity of Costa Rican rainforest within 2.3% of aboveground biomass values, based on comparisons with CARBONO inventory plot data. The Costa Rica FORMIND simulation was then used to investigate the relationship between canopy height and aboveground biomass (AGB), leaf area index (LAI) and gross primary productivity (GPP) at different spatial scales (20 × 20 m, 60 × 60 m, 100mx100m). The relationship between aboveground biomass and height is of particular importance toward the calibration of various remote sensing products including lidar and radar, whereas the LAI and GPP relationships are understudied in this context. We found that the relationship between all three variables and height varies considerably: the relationship is stronger at finer scales and weaker at coarser resolution. However, in all three comparisons, RMSE also decreased as scales coarsened, with the largest difference shown between 100 m and 10 m resolutions in relating AGB to Lorey's height (R^2 decreased by 0.3; RMSE decreased by 114.5 Mg/ha). This suggests that a trade-off between accuracy and precision exists, and further highlights the importance of spatial scale in determining the relatability of forest structure variables.

A.H.Armstrong↗

Skillful Seasonal Forecasts of Land Carbon Uptake in Northern Mid- and High Latitudes

Here we present a first look at the Gross Primary Production (GPP) forecast skill levels achievable with a state-of-the-art subseasonal-to-seasonal (S2S) forecast system. Using NASA’s retrospective S2S ensemble forecast in conjunction with a terrestrial biosphere model, and using an independent, remote sensing-based dataset for validation, we demonstrate an ability to accurately forecast spring-summer carbon uptake at multi-month leads. Averaged across mid-and high latitudes of the Northern Hemisphere land, the GPP forecast initialized on January 1 produces statistically significant skill through summer. The skill achieved, however, is spatially variable, with some regions appearing to extract skill from accurate forecasts of snowpack removal and others extracting skill from the initialization of carbon and nitrogen states. Our results reveal some heretofore unexplored facets of climate predictability and provide a look at what might be possible with future S2S forecast systems that are fully integrated with biogeochemical cycles.

Eunjee Lee↗

Climate drives modeled forest carbon cycling resistance and resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond to (resistance) and recover from (resilience) rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the Ecosystem Demography-2.2 (ED-2.2) model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900 (approximating the clear-cut conditions at the time), and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes such as gross primary production (GPP), and interacted with one another. We then calculated resistance and resilience, two dimensions of ecosystem stability. Modeled GPP exhibited a two-fold decrease in mean resistance across disturbance severities of 45%, 65%, and 85% mortality; conversely, resilience increased by a factor of two with increasing disturbance severity. This pattern held for net primary production and net ecosystem production, indicating a trade-off in which greater initial declines were followed by faster recovery. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

Kalyn Dorheim↗

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (−30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (−19%) and ecosystem (−7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Félicien Meunier↗

The Response of the Amazon Ecosystem to the Photosynthetically Active Radiation Fields: Integrating Impacts of Biomass Burning Aerosol and Clouds in the NASA GEOS Earth System Model

The Amazon experiences fires every year, and the resulting biomass burning aerosols, together with cloud particles, influence the penetration of sunlight through the atmosphere, increasing the ratio of diffuse to direct photosynthetically active radiation (PAR) reaching the vegetation canopy and thereby potentially increasing ecosystem productivity. In this study, we use the NASA Goddard Earth Observing System (GEOS) model with coupled aerosol, cloud, radiation, and ecosystem modules to investigate the impact of Amazon biomass burning aerosols on ecosystem productivity, as well as the role of the Amazon’s clouds in tempering this impact. The study focuses on a seven-year period (2010-2016) during which the Amazon experienced a variety of dynamic environments (e.g., La Niña, normal years, and El Niño). The direct radiative impact of biomass burning aerosols on ecosystem productivity—called here the aerosol diffuse radiation fertilization effect —is found to increase Amazonian Gross Primary Production (GPP) by 2.6% via a 3.8% increase in diffuse PAR (DFPAR) despite a 5.4% decrease in direct PAR (DRPAR) on multiyear average during burning seasons. On a monthly basis, this increase in GPP can be as large as 9.9% (occurring in August 2010). Consequently, the net primary production (NPP) in Amazon is increased by 1.5%, or ~92 Tg C a-1– equivalent to ~37% of the average carbon lost due to Amazon fires over the seven years considered. Clouds, however, strongly regulate the effectiveness of the aerosol diffuse radiation fertilization effect. The efficiency of this fertilization effect is the highest in cloud-free conditions and linearly decreases with increasing cloud amount until the cloud fraction reaches ~0.8, at which point the aerosol-influenced light changes from being a stimulator to an inhibitor of plant growth. Nevertheless, interannual changes in the overall strength of the aerosol diffuse radiation fertilization effect are primarily controlled by the large interannual changes in biomass burning aerosols rather than by changes in cloudiness during the studied period.

biomass burning aerosols on ecosystem productivity↗

Climate drives modeled forest carbon cycling resistance and resilience in the Upper Great Lakes Region, USA

Forests dominate the global terrestrial carbon budget, but their ability to continue doing so in the face of a changing climate is uncertain. A key uncertainty is how forests will respond to (resistance) and recover from (resilience) rising levels of disturbance of varying intensities. This knowledge gap can optimally be addressed by integrating manipulative field experiments with ecophysiological modeling. We used the Ecosystem Demography-2.2 (ED-2.2)model to project carbon fluxes for a northern temperate deciduous forest subjected to a real-world disturbance severity manipulation experiment. ED-2.2 was run for 150 years, starting from near bare ground in 1900(approximating the clear-cut conditions at the time),and subjected to three disturbance treatments under an ensemble of climate conditions. Both disturbance severity and climate strongly affected carbon fluxes such as gross primary production (GPP),and interacted with one another. We then calculated resistance and resilience, two dimensions of ecosystem stability. Modeled GPP exhibited a two-fold decrease in mean resistance across disturbance severities of 45%, 65%, and 85% mortality; conversely, resilience increased by a factor of two with increasing disturbance severity. This pattern held for net primary production and net ecosystem production, indicating a trade-off in which greater initial declines were followed by faster recovery. Notably, however, heterotrophic respiration responded more slowly to disturbance, and its highly variable response was affected by different drivers. This work provides insight into how future conditions might affect the functional stability of mature forests in this region under ongoing climate change and changing disturbance regimes.

Kalyn Dorheim↗

Development of Carbon Flux Model Using ABI Data Over the Conterminous US

The satellite-driven carbon flux estimation has been playing important role to estimate continental-scale carbon budget. One of the biggest recent advances in the satellite-driven carbon flux modeling is utilization of high-frequent geostationary satellites to estimate diurnal cycle in carbon fluxes. The satellite based carbon flux estimation used the polar orbiting satellite sensors (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)), which allow us to observe target regions only once during the day. The new generation of geostationary satellite sensors provide frequent observations, often less than every 10 minutes. Here, we use GOES Advanced Baseline Imager (ABI) data to estimate hourly NEE over the conterminous US. We used the Terrestrial Observation Prediction System (TOPS) model for estimating hourly NEE. TOPS is a diagnostic ecosystem process model that simulates the fluxes of carbon and water through vegetation in response to climate variability. For the climate input, we developed hourly climate data using the same algorithm with NASA Earth Exchange Gridded Daily Meteorology (NEX-GDM) datasets based on machine learning techniques. The hourly climate data includes precipitation, maximum temperature, minimum temperature, dew point temperature, and solar radiation were derived from the Geostationary observations. The spatial patterns of ecosystem parameters used in TOPS are optimized using satellite Solar Induced Fluorescence (SIF) data. The high frequency GPP estimations from geostationary satellite sensors make it comparable to the instantaneous SIF data than daily GPP. We also used Ameriflux data for optimization of model parameters and the validation of the output. The derived data addresses the diurnal dynamics of carbon cycling at large scales and should help in reducing the uncertainties in carbon budget studies.

geostationary satellite↗

Data and Scripts associated with a manuscript on ecosystem responses to wildfires in the Columbia River Basin

This data package is associated with the publication “Ecosystem leaf area, gross primary production, and evapotranspiration responses to wildfire in the Columbia River Basin” submitted to Biogeosciences (Shi et al., 2024; doi: 10.22541/au.171053013.30286044/v1). In this research, data products, leaf area index (LAI), gross primary production (GPP), and evapotranspiration (ET), from the Moderate Resolution Imaging Spectroradiometer (MODIS) are used to quantify the resistance and resilience of different ecosystem types in the Columbia River Basin (CRB). A machine learning algorithm, random forest (RF), was used to examine the impacts of precipitation, vapor pressure deficit (VPD), and burn severity from Monitoring Trends in Burn Severity (MTBS) on ecosystem resilience. The data package includes the processed MODIS data products, precipitation, VPD, and burn severity in 138 fire regions in CRB and the input files for RF model training. This data package includes six folders. The MODIS products are included in three MODIS_* folders with shell scripts for data clipping and *ncl files for data processing: (1) “/MODIS_LAI_CRB”; (2) “/MODIS_GPP_CRB”; and (3) “/MODIS_ET_CRB”. All the processed data for each fire event are NetCDF formatted. The MTBS burn severity data and the shell and *ncl scripts used for data processing are in the folder named (4) “MTBS_fire”. The ERA meteorological fields and the data processing scritps are in (5) “ERA_Var_CR”. All the scripts for figure development are in the format of *ncl and in the folder (6) “paper_scripts”. See the file ending in “flmd.csv” for a list of all files contained in this data package and descriptions for each. Tabular column headers and units are described in the data dictionary file ending in “dd.csv”.

54 ENVIRONMENTAL SCIENCES↗

Data and scripts associated with a manuscript analyzing ELM-FATES parameter sensitivity under pre-fire and postfire scenarios using machine learning

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “Fire Severity-Dependent Shifts in Vegetation Parameter Sensitivity: A Pre- and Post-Fire Analysis Using ELM-FATES and Explainable AI” submitted to Journal of Advances in Modeling Earth Systems (Zahura et al. 2026). The study examines vegetation physiological parameters controlling pre-fire and post-fire vegetation dynamics. To support this analysis, 73 vegetation parameters in Functionally Assembled Terrestrial Ecosystem Simulator (FATES) (Fisher et al., 2018) , which is coupled with E3SM (Energy Exascale Earth System Model) land model (ELM, ELM-FATES), were perturbed using a Sobol sequence to generate 1,024 ensemble members for two plant functional types: needleleaf evergreen extratropical trees (NEET) and C3 grass. Simulations were conducted for the pre-fire period (2016) and post-fire period (2018–2023). Burn severity was represented by modifying the Nesterov index in FATES to 75,000, 150,000, and 300,000 for low, moderate, and high severity, respectively. A no-fire scenario was also included. Simulations were performed for 16 grid cells in the American River Watershed across different burn severities and plant functional types. XGBoost (eXtreme Gradient Boosting) models were trained using the parameter ensembles and ELM-FATES-simulated outputs, including leaf area index (LAI), gross primary productivity (GPP), aboveground biomass, vegetation evaporation, transpiration, and soil evaporation. Models were trained separately for each year and burn severity, followed by SHAP (SHapley Additive exPlanations) analysis to identify changes in dominant parameters after fire disturbance. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. The data package contains the ELM-FATES simulation data. The scripts and data related to the analysis will be added later. The inputs and outputs from ELM-FATES are inside the “FATES” folder. “FATES_domain_surface” contains the domain and surface netcdfs that were used to run ELM-FATES in the study area. “FATES_parameters” contains the 1024 ensembles that were generated using Sobol sequence. “FATES_outputs” folder contains ELM-FATES simulated variables. All files are .csv and .nc (NetCDF).

Aboveground biomass↗

Hydroclimatic Controls on the Means and Variability of Vegetation Phenology and Carbon Uptake

Long-term, global offline (land-only) simulations with a dynamic vegetation phenology model are used to examine the control of hydroclimate over vegetation-related quantities. First, with a control simulation, the model is shown to capture successfully (though with some bias) key observed relationships between hydroclimate and the spatial and temporal variations of phenological expression. In subsequent simulations, the model shows that: (i) the global spatial variation of seasonal phenological maxima is controlled mostly by hydroclimate, irrespective of distributions in vegetation type, (ii) the occurrence of high interannual moisture-related phenological variability in grassland areas is determined by hydroclimate rather than by the specific properties of grassland, and (iii) hydroclimatic means and variability have a corresponding impact on the spatial and temporal distributions of gross primary productivity (GPP).

GPP↗