Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LaI”

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 199 records · Page 11

Analysis of Leaf Area Index and Fraction of PAR Absorbed by Vegetation Products from the Terra MODIS Sensor: 2000-2005

The analysis of two years of Collection 3 and five years of Collection 4 Terra Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) data sets is presented in this article with the goal of understanding product quality with respect to version (Collection 3 versus 4), algorithm (main versus backup), snow (snow-free versus snow on the ground), and cloud (cloud-free versus cloudy) conditions. Retrievals from the main radiative transfer algorithm increased from 55% in Collection 3 to 67% in Collection 4 due to algorithm refinements and improved inputs. Anomalously high LAI/FPAR values observed in Collection 3 product in some vegetation types were corrected in Collection 4. The problem of reflectance saturation and too few main algorithm retrievals in broadleaf forests persisted in Collection 4. The spurious seasonality in needleleaf LAI/FPAR fields was traced to fewer reliable input data and retrievals during the boreal winter period. About 97% of the snow covered pixels were processed by the backup Normalized Difference Vegetation Index-based algorithm. Similarly, a majority of retrievals under cloudy conditions were obtained from the backup algorithm. For these reasons, the users are advised to consult the quality flags accompanying the LAI and FPAR product.

Yang, Wenze↗

Potential Scientific Output from the Visible Channels of EPIC Spectroradiometer as Part of the DSCOVR L1 Mission

In addition to 4 UV channels, the EPIC spectroradiometer will provide measurements at 6 visible (and near IR) channels (443,551,680,687.75, 764 and 779.9 nm) at roughly 10 km spatial resolution. The scattering angles are near backscattering and vary between 165 and 176 degree. Two pairs {680 and 687.75 nm} and {764 and 779.9 nm} represent 02 B- and A-bands (and their references) channels; they will be used for cloud height measurements over land and ocean. The B-band channel will contribute to the more absorbing A-band measurements over bright vegetation. A pair {680 and 779.9 nm} will be used for retrieving vegetation properties. Due to the special EPIC geometry, the illuminated part of the leaves will be always observed. As a result, in additional to the traditional Leaf Area Index (LAI), these observations for the first time will provide the sunlit fraction of LAI. Since the sunlit and shaded leaves exhibit different photosynthetic response to incident radiation, these measurements will help to improve global ecological and biogeochemistry models. Finally, a pair {443 and 551 nm} will be used for atmospheric correction. As a by-product of the atmospheric correction algorithm, we also expect to get aerosol optical thickness (AOT) and surface bidirectional reflection function (BRF). The presentation will briefly overview the proposed science algorithms.

Marshak, Alexander↗

Hydroclimatic Controls over Global Variations in Phenology and Carbon Flux

The connection between phenological and hydroclimatological variations are quantified through joint analyses of global NDVI, LAI, and precipitation datasets. The global distributions of both NDVI and LAI in the warm season are strongly controlled by three quantities: mean annual precipitation, the standard deviation of annual precipitation, and Budyko's index of dryness. Upon demonstrating that these same basic (if biased) relationships are produced by a dynamic vegetation model (the dynamic vegetation and carbon storage components of the NCAR Community Land Model version 4 combined with the water and energy balance framework of the Catchment Land Surface Model of the NASA Global Modeling and Assimilation Office), we use the model to perform a sensitivity study focusing on how phenology and carbon flux might respond to climatic change. The offline (decoupled from the atmosphere) simulations show us, for example, where on the globe a given small increment in precipitation mean or variability would have the greatest impact on carbon uptake. The analysis framework allows us in addition to quantify the degree to which climatic biases in a free-running GCM are manifested as biases in simulated phenology.

Koster, Randal↗

Simulations of Seasonal and Latitudinal Variations in Leaf Inclination Angle Distribution: Implications for Remote Sensing

The leaf inclination angle distribution (LAD) is an important characteristic of vegetation canopy structure affecting light interception within the canopy. However, LADs are difficult and time consuming to measure. To examine possible global patterns of LAD and their implications in remote sensing, a model was developed to predict leaf angles within canopies. Canopies were simulated using the SAIL radiative transfer model combined with a simple photosynthesis model. This model calculated leaf inclination angles for horizontal layers of leaves within the canopy by choosing the leaf inclination angle that maximized production over a day in each layer. LADs were calculated for five latitude bands for spring and summer solar declinations. Three distinct LAD types emerged: tropical, boreal, and an intermediate temperate distribution. In tropical LAD, the upper layers have a leaf angle around 35 with the lower layers having horizontal inclination angles. While the boreal LAD has vertical leaf inclination angles throughout the canopy. The latitude bands where each LAD type occurred changed with the seasons. The different LADs affected the fraction of absorbed photosynthetically active radiation (fAPAR) and Normalized Difference Vegetation Index (NDVI) with similar relationships between fAPAR and leaf area index (LAI), but different relationships between NDVI and LAI for the different LAD types. These differences resulted in significantly different relationships between NDVI and fAPAR for each LAD type. Since leaf inclination angles affect light interception, variations in LAD also affect the estimation of leaf area based on transmittance of light or lidar returns.

Seasonal↗

Spin-Up and Tuning of the Global Carbon Cycle Model Inside the GISS ModelE2 GCM

Planetary carbon cycle involves multiple phenomena, acting at variety of temporal and spacial scales. The typical times range from minutes for leaf stomata physiology to centuries for passive soil carbon pools and deep ocean layers. So, finding a satisfactory equilibrium state becomes a challenging and computationally expensive task. Here we present the spin-up processes for different configurations of the GISS Carbon Cycle model from the model forced with MODIS observed Leaf Area Index (LAI) and prescribed ocean to the prognostic LAI and to the model fully coupled to the dynamic ocean and ocean biology. We investigate the time it takes the model to reach the equilibrium and discuss the ways to speed up this process. NASA Goddard Institute for Space Studies General Circulation Model (GISS ModelE2) is currently equipped with all major algorithms necessary for the simulation of the Global Carbon Cycle. The terrestrial part is presented by Ent Terrestrial Biosphere Model (Ent TBM), which includes leaf biophysics, prognostic phenology and soil biogeochemistry module (based on Carnegie-Ames-Stanford model). The ocean part is based on the NASA Ocean Biogeochemistry Model (NOBM). The transport of atmospheric CO2 is performed by the atmospheric part of ModelE2, which employs quadratic upstream algorithm for this purpose.

Climate model↗

Detection of Chlorophyll and Leaf Area Index Dynamics from Sub-weekly Hyperspectral Imagery

Temporally rich hyperspectral time-series can provide unique time critical information on within-field variations in vegetation health and distribution needed by farmers to effectively optimize crop production. In this study, a dense time series of images were acquired from the Earth Observing-1 (EO-1) Hyperion sensor over an intensive farming area in the center of Saudi Arabia. After correction for atmospheric effects, optimal links between carefully selected explanatory hyperspectral vegetation indices and target vegetation characteristics were established using a machine learning approach. A dataset of in-situ measured leaf chlorophyll (Chll) and leaf area index (LAI), collected during five intensive field campaigns over a variety of crop types, were used to train the rule-based predictive models. The ability of the narrow-band hyperspectral reflectance information to robustly assess and discriminate dynamics in foliar biochemistry and biomass through empirical relationships were investigated. This also involved evaluations of the generalization and reproducibility of the predictions beyond the conditions of the training dataset. The very high temporal resolution of the satellite retrievals constituted a specifically intriguing feature that facilitated detection of total canopy Chl and LAI dynamics down to sub-weekly intervals. The study advocates the benefits associated with the availability of optimum spectral and temporal resolution spaceborne observations for agricultural management purposes.

Houborg, Rasmus↗

GEONEX: Land Monitoring From a New Generation of Geostationary Satellite Sensors

The latest generation of geostationary satellites carry sensors such as ABI (Advanced Baseline Imager on GOES-16) and the AHI (Advanced Himawari Imager on Himawari) that closely mimic the spatial and spectral characteristics of Earth Observing System flagship MODIS for monitoring land surface conditions. More importantly they provide observations at 5-15 minute intervals. Such high frequency data offer exciting possibilities for producing robust estimates of land surface conditions by overcoming cloud cover, enabling studies of diurnally varying local-to-regional biosphere-atmosphere interactions, and operational decision-making in agriculture, forestry and disaster management. But the data come with challenges that need special attention. For instance, geostationary data feature changing sun angle at constant view for each pixel, which is reciprocal to sun-synchronous observations, and thus require careful adaptation of EOS algorithms. Our goal is to produce a set of land surface products from geostationary sensors by leveraging NASA's investments in EOS algorithms and in the data/compute facility NEX. The land surface variables of interest include atmospherically corrected surface reflectances, snow cover, vegetation indices and leaf area index (LAI)/fraction of photosynthetically absorbed radiation (FPAR), as well as land surface temperature and fires. In order to get ready to produce operational products over the US from GOES-16 starting 2018, we have utilized 18 months of data from Himawari AHI over Australia to test the production pipeline and the performance of various algorithms for our initial tests. The end-to-end processing pipeline consists of a suite of modules to (a) perform calibration and automatic georeference correction of the AHI L1b data, (b) adopt the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm to produce surface spectral reflectances along with compositing schemes and QA, and (c) modify relevant EOS retrieval algorithms (e.g., LAI and FPAR, GPP, etc.) for subsequent science product generation. Initial evaluation of Himawari AHI products against standard MODIS products indicate general agreement, suggesting that data from geostationary sensors can augment low earth orbit (LEO) satellite observations.

geostationary↗

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↗

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↗

Determination of Chlorophyll Content in Selected Grass Communities of Krkonoše Mts. Tundra Based on Laboratory Spectroscopy and Aerial Hyperspectral Data

The study focuses on the determination of chlorophyll content in four prevailing grasses in the relict arctic-alpine tundra located in the Krkonoše Mountains National Park, Czech Republic. We compared two methods for determination of leaf chlorophyll content (LCC) – spectrophotometric determination in the laboratory, and the LCC assessed by fluorescence portable chlorophyll meter CCM-300. Relationships were established between the LCCs and vegetation indices calculated from leaf spectra acquired with contact probe coupled with an ASD FieldSpec4 Wide-Res spectroradiometer. Canopy chlorophyll contents (CCC) were computed from the LCCs and green leaf area index (LAI), and modelled based on the field spectra measured by the spectroradiometer and the hyperspectral images acquired by Headwall Nano-Hyperspec® mounted on the DJI Matrice 600 Pro drone. The calculations are performed on datasets acquired in June, July and August 2020 together and separately for species and months. In general, the correlations based on June datasets work the best at both levels: median R 2 for all indices was 0.52 for all species together at leaf level and median R 2 = 0.47 at the canopy level (vegetation indices computed from field spectra). Canopy chlorophyll content map was created based on the results of stepwise multiple linear regression. The R 2 was 0.42 when using four wavelengths from the red and red edge spectral region. We attribute the weak model performance to a combination of several factors: leaf structure may bias LCC from laboratory measurements, effects of LAI variability on CCC, and the sampling design, probably not covering the whole phenology equally for all studied species.

L. Cervena↗

Evaluating Differences Among Crop Models in Simulating Soybean in-Season Growth

Crop models are useful tools for simulating agricultural systems that require continued model development and testing to increase their robustness and improve how they describe our current understanding of processes. Coordinated and “blind” evaluation of multiple models using same protocols and experimental datasets provides unique opportunities to further improve models and enhance their reliability. For soybean [Glycine max (L.) Merr.], there has been limited coordinated multi-model evaluations for the simulation of in-season plant growth dynamics. We evaluated ten dynamic soybean crop models for their simulation of in-season plant growth using data from five experiments conducted in Argentina, Brazil, France, and USA. We evaluated models after a Blind (using only phenology data) and a Full calibration (with in-season and end-of-season variables). Calibration reduced model uncertainty by reducing standard bias for the simulation of in-season variables (biomass, leaf, pod, and stem weights, and leaf area index, LAI). However, we found that most models had difficulty in reproducing leaf growth dynamics, with normalized root mean squared error (nRMSE) of 56% for leaf weight and 43% for LAI (across locations and models after Full calibration). Models with different levels of complexity and experience were capable of simulating final seed yield at maturity with reasonable accuracy (nRMSE of 8–31% after Full calibration). However, the nRMSE for pod weight (of 17–64% after Full calibration) was two-fold larger than that of seed yield. Moreover, the models differed in how they simulated timing from sowing to beginning seed growth (47–93 days) and effective seed filling period (18–54 days), owing to model structural differences in defining the reproductive developmental stages. Overall, we identified the following processes that can benefit from further model improvement: leaf expansion and senescence, reproductive phenology, and partitioning to reproductive growth. Simulation of pod wall tissue and individual seed cohorts is another aspect that many models currently lack. Model improvement can benefit from high-temporal resolution experimental datasets that concurrently account for phenology, plant growth, and partitioning. Further, we recommend collecting reproductive phenology in the field consistent with actual dry matter allocation to organs in the models and collecting multiple observations of seed and pod weight to aid model improvement for simulation of seed growth and yield formation.

Agricultural Model Intercomparison and Improvement↗

An Observation-Driven Framework for Modeling Post-Fire Hydrologic Response: Evaluation for Two Central California Case Studies

In a warming climate, wildfires are becoming increasingly common, especially in semi-arid environments. Wildfires can disrupt forest ecosystems and induce changes to the land surface. Collectively, these impacts can alter the hydrologic response of a catchment following a fire, resulting in increased potential for surface runoff, reduced evapotranspiration, and, ultimately, a higher risk for flash flooding and mass wasting. The timescale of post-fire recovery of hydrological processes to return to pre-fire conditions is not well established due to the lack of ground measurements. Accurate characterization of the impacts of fire on hydrologic response is also challenging to simulate, given the complex interplay of various processes. Here, we present a generalized framework to quantify the impacts of wildfire on runoff generation. We consider the disturbances in the vegetation and soil as the two main factors contributing to post-fire floods. Using an ensemble modeling structure to account for parameter uncertainty, remotely sensed leaf area index (LAI) is assimilated into a land surface model (LSM) to simulate vegetation disturbance, and the maximum land surface saturation LSM parameter is decreased to parameterize the soil disturbance following observed fires. We consider the impacts of fire-induced changes to LAI and soil saturation on hydrologic states like runoff and evapotranspiration for two case studies. These case studies demonstrate the general applicability of hydrophobicity formulation to serve as a guideline for exploring the range of hydrologic responses post-fire.

Fire Hydrology↗

Maize and Sorghum Establishment and Yield following Pre-Emergence Waterlogging

In 2020, early-season extreme precipitation events occurred following the planting of Sorghum bicolor (L.) Moench and Zea mays L. in central Illinois that caused ponding. Following the first rainfall event 50m transects were established to assess the waterlogging effects on seedling emergence and crop yields. Soil moisture, emergence, stem and tiller count, LAI, and yield were measured at various points in the season along these transects.

Emergence↗

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, model inputs, and analysis scripts associated with a manuscript on stream intermittency controls across spatial scales in Pacific Northwest watersheds

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 "Hydroclimatic Memory and Watershed Template Shape Stream Intermittency: Multi-scale Attribution Using Process-based Simulation and Explainable ML" by Niroula et al. (2026), submitted to Water Resources Research (WRR). The study investigates the dominant controls on stream intermittency across local, reach, and watershed scales using a coupled process-based simulation and explainable machine-learning framework. Long-term daily simulations from the Advanced Terrestrial Simulator (ATS) were used to generate wetness states and ponded-depth responses over river-corridor cells. These ATS outputs were then aggregated across scales and used to train XGBoost (eXtreme Gradient Boosting) models. SHAP (SHapley Additive exPlanations) was applied to quantify the relative importance of hydroclimatic forcings, watershed template attributes, and antecedent-memory effects in shaping intermittency behavior. The analysis is carried out for three contrasting Pacific Northwest watersheds: Oak Creek (OCW), American River Watershed (ARW), and H.J. Andrews (HJA). Across these testbeds, the package contains ATS-ready watershed inputs, ATS run configuration and selected output files, model-evaluation data products, intermittency-analysis datasets, machine-learning target-feature tables, SHAP outputs, and notebooks used to organize, analyze, and visualize results. At a high level, the package documents a workflow in which ATS provides the physically based simulation backbone and explainable machine learning is used as a post-processing attribution tool. The contents are intended to support interpretation of the manuscript figures and results, provide context for how intermittency metrics were generated at multiple scales, and preserve the key artifacts needed to understand and reuse the analysis workflow. The package contains a high-level directory summary file (`summary.txt`) and four main content folders (1) `evaluation_plots` contains evaluation figures and supporting evaluation datasets; (2) `intermittency_plots` contains intermittency-focused analysis notebook and prepared datasets; (3) `ml-training-and-shap_values_plots` contains ML training inputs, SHAP outputs, and figure-generation notebooks; and (4) `watershed_mesh_and_ats_input` contains ATS model setup materials, forcing inputs, geometry, and selected run files. More specifically, the `evaluation_plots` folder contains the notebook used for ATS evaluation plotting and site-specific evaluation datasets. These include evapotranspiration and water-balance products for three watersheds, as well as an Oak Creek field-measurement discharge file. The `intermittency_plots` folder contains the notebook used for intermittency analysis and the prepared datasets used to analyze intermittent and non-intermittent wetness behavior across the study watersheds. The `ml-training-and-shap_values_plots` folder contains notebooks and outputs for the machine-learning and explainability workflow. This includes the main XGBoost and SHAP notebook(s), a beeswarm plotting notebook, target-feature tables for machine-learning training, SHAP summary tables, and per-sample SHAP value archives. The `watershed_mesh_and_ats_input` folder contains ATS-related watershed inputs and supporting materials. This includes mesh and shape products, ATS-readable LAI and meteorological forcing inputs, selected ATS spinup and transient-run files, and a watershed workflow example notebook. Subdirectories are organized by watershed where applicable.All files are .cpg (codepage files), .csv (comma-separated values), .dbf (database files), .exo (Exodus mesh format), .h5 (HDF5 format), .ipynb (Jupyter notebooks), .pkl (Python pickle), .prj (projection files), .sh (shell scripts), .shp (shapefile geometry), .shx (shapefile index), .txt (text files), or .xml (markup data).

Advanced Terrestrial Simulator↗

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↗