Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Land surface”

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 1,189 records · Page 66

Portland Urban Development: Quantifying and Visualizing Urban Heat with Compounding Vulnerabilities to Support Community Depaving Initiatives

Urban heat is a pressing concern in Portland, Oregon as climate change induced heat waves increase. Cities experience higher temperatures due to the urban heat island effect (UHI), and environmental injustice and disenfranchisement in minority communities expose low-income and Black, Indigenous, and People of Color (BIPOC) residents to more extreme and debilitating heat events. Our team identified Portland’s communities on the frontlines of urban heat impacts by overlapping environmental and social vulnerabilities using NASA Earth observations. We partnered with Depave, a Portland-based nonprofit that works alongside communities to replace pavement with greenspace in historically disenfranchised areas. Using Landsat 8 Thermal Infrared Sensor (TIRS) imagery, we mapped Land Surface Temperature (LST) and developed a heat-specific Social Vulnerability Index (SVI) through a Principal Component Analysis (PCA) to identify Portland’s communities with the highest potential heat vulnerability. Then, we calculated the temperature change of depaving in six case studies to quantify Depave's efforts in heat mitigation and environmental justice. Our analysis demonstrated that, throughout Portland, there are frontline communities experiencing high potential social vulnerability to extreme temperatures due to environmental injustices and over-pavement. Finally, Depave’s impact on urban heat is observable and quantifiable using remote-sensing data and tools, with an average of 1ºF LST decrease across the six case studies. We illustrated the significance of local urban heat mitigation efforts and propose next steps for conducting inclusive and intentional research that highlights the lived experiences and resilience of frontline communities.

Environmental justice↗

Storm and Annual Time Scale Hydrological Data for the Russian River Watershed 1996-2022

These files contains observed and simulated hydrological data (discharge, precipitation) that have been aggregated to the storm-event scale and to the annual time scale for the Russian River Watershed (RRW), California. Observed data are obtained from 12 different USGS hydrological stations located throughout the watershed, and discharge is simulated using the GR5H hourly hydrological model. Precipitation data are obtained from the NASA NLDAS NOAH Community Land Surface Model. All simulated and observed discharge are aggregated to the storm and annual time scales. We investigated changes in watershed hydrological conditions in the Russian River Watershed, a Mediterranean, drought prone, wildfire-adapted ecosystem, following eleven wildfires that occurred from 2017-2020. We ask two research questions: 1) How do interacting wildfire events, drought, and atmospheric rivers impact hydrological conditions of the watershed and in particular streamflow?, and 2) What percentage of new and/or cumulative wildfire disturbance is required to initiate hydrological change? We hypothesize that sub-watersheds of the RRW Mediterranean ecoregion have not burned beyond an intrinsic, and still unknown, threshold required to initiate change. Using a series of paired burned/unburned catchments nested within the larger watershed, we examined temporal and spatial patterns of pre-and-post wildfire water hydrology using a rainfall runoff hydrological model compared with data.

54 ENVIRONMENTAL SCIENCES↗

Residential segregation and outdoor urban moist heat stress disparities in the United States

The combined impact of urbanization-induced warming and drying on large-scale heat stress disparities remains unknown, with multicity studies using satellite-derived land surface temperature as a proxy for these disparities. Here, using high-resolution urban-resolving numerical model simulations for 2014–2018, we find pervasive disparities in all-sky average maximum summertime air temperature and moist heat stress metrics across US cities, with higher outdoor heat stress exposure in poorer and primarily non-white census tracts. Ninety-four percent of the US urban population (228 million) live in cities where heat stress burdens the poor, with heat stress inequities between white and non-white populations strongly associated with residential segregation. Similarly, historically redlined neighborhoods show higher heat stress than their non-redlined counterparts, demonstrating how historical segregation relates to present-day environmental inequalities. Our results provide quantitative estimates of physiologically relevant heat stress disparities at the US national scale and highlight potential biases when using satellites as a proxy for these.

54 ENVIRONMENTAL SCIENCES↗

Historical global Earth System Model simulations with E3SM's dynamic root module

This data package contains global-scale Earth System Model simulations run with DOE's E3SM Model. The package includes one set of simulations run in the default configuration and a second set of simulations run with the dynamic root module enabled. These two datasets can be used to assess the impact that dynamic roots have on Land Surface fluxes such as carbon and water. The outputs are provided here on a 0.5 x 0.5 global grid at monthly resolution. The simulations were run with atmospheric forcing from the Global Soil Wetness Project. We provide specifically here Gross Primary Production, Transpiration, relative root fraction per soil layer as well as Air Temperature, Precipitation and all the parameterizations used for the simulations. Our purpose for generating and analyzing these simulations was to assess the role that root dynamics and foraging for water has in the recovery timescale of ecosystems following climate perturbation. All data provided are in the "netcdf" format using standard CF-1 (Climate and Forecast) Convention. The data are broken up into 20 year chunks to facilitate easier read-in. The data are all machine readable using various software platforms including Matlab, NCO, Panoply and GrADS.

54 ENVIRONMENTAL SCIENCES↗

National Climate Assessment - Land Data Assimilation System (NCA-LDAS) Data and Services at NASA GES DISC

The National Climate Assessment-Land Data Assimilation System (NCA-LDAS) is an Integrated Terrestrial Water Analysis, and is one of NASAs contributions to the NCA of the United States. The NCA-LDAS has undergone extensive development, including multi-variate assimilation of remotely-sensed water states and anomalies as well as evaluation and verification studies, led by the Goddard Space Flight Centers Hydrological Sciences Laboratory (HSL). The resulting NCA-LDAS data have recently been released to the general public and include those from the Noah land-surface model (LSM) version 3.3 (Noah-3.3) and the Catchment LSM version Fortuna-2.5 (CLSM-F2.5). Standard LSM output variables including soil moistures temperatures, surface fluxes, snow cover depth, groundwater, and runoff are provided, as well as streamflow using a river routing system. The NCA-LDAS data are archived at and distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). The data can be accessed via HTTP, OPeNDAP, Mirador search and download, and NASA Earth data Search. To further facilitate access and use, the NCA-LDAS data are integrated into the NASA Giovanni, for quick visualization and analysis, and into the Data Rods system, for retrieval of time series of long time periods. The temporal and spatial resolutions of the NCA-LDAS data are, respectively, daily-averages and 0.125x0.125 degree, covering North America (25N 53N; 125W 67W) and the period January 1979 to December 2015. The data files are in self-describing, machine-independent, CF-compliant netCDF-4 format.

data assimiliation↗

Improving the Fire Weather Index System for Peatlands Using Peat-Specific Hydrological Input Data

The Canadian Fire Weather Index (FWI) system, even though originally developed and calibrated for an upland Jack pine forest, is used globally to estimate fire danger for any fire environment. However, for some environments, such as peatlands, the applicability of the FWI in its current form, is often questioned. In this study, we replaced the original moisture codes of the FWI with hydrological estimates resulting from the assimilation of satellite-based L-band passive microwave observations into a peatland-specific land surface model. In a conservative approach that maintains the integrity of the original FWI structure, the distributions of the hydrological estimates were first matched to those of the corresponding original moisture codes before replacement. The resulting adapted FWI, hereafter called FWIpeat, was evaluated using satellite-based information on fire presence over boreal peatlands from 2010 through 2018. Adapting the FWI with model- and satellite-based hydrological information was found to be beneficial in estimating fire danger, especially when replacing the deeper moisture codes of the FWI. For late-season fires, further adaptations of the fine fuel moisture code show even more improvement due to the fact that late-season fires are more hydrologically driven. The proposed FWIpeat should enable improved monitoring of fire risk in boreal peatlands.

Canadian Fire Weather Index↗

CHESS 2025: Waveform LiDAR data from NEON AOP surveys

This dataset provides Level 1 (L1) full-waveform light detection and ranging (LiDAR) data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. Waveform LiDAR data can provide more detailed information about objects on the ground than discrete point clouds typically do, and they are often used for granular target segmentation and characterization of subcanopy vegetation. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary waveform LiDAR data delivered by NEON and are provided per flightline in compressed Pulsewaves format, an open-source binary file standard. A Pulsewaves object comprises a two files: a pulse (.pls) file, which stores the geographic origin, outgoing vector, and metadata for every laser pulse emitted by the scanner, and a wave file (.wvs), which stores the sequential amplitude samples of the outgoing pulse and the returning signals. The files are published here in their compressed forms (.plz, .wvz). All waveform data were processed following the theoretical workflow described in the NEON L0-to-L1 Waveform LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022a); however, the Pulsewaves output format differs from a legacy format described in that document. Waveform amplitude samples are recorded at 1 nanosecond intervals. All coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Waveform data for the UPTA survey area were collected without incident and the published records are complete. However, both the ALMO and CRBU collections experienced issues that resulted in incomplete data for those areas. On collection day 2018-06-16 a hardware failure caused the waveform digitizer to lose data from the eastern edge of the ALMO site (Figure 22). The waveform data for flightlines 2–20 could not be extracted from the digitizer, and the data proved unrecoverable. As a result, a portion of the site does not have coverage with waveform data. Although no hardware failure was observed during collection over the CRBU area, final waveform files generated by vendor software contained only ~25% of the expected number of return pulses. After discovery, NEON initiated troubleshooting with the vendor. The root cause of the data ablation had not been identified at the time of publication. Additional data will be published in an update to this package if further recovery proves successful. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

NASA SPoRT Land Information System Products for Soil Moisture Analysis

The NASA Short-term Prediction Research and Transition (SPoRT) Program has been producing a near real-time instance of NASA’s Land Information System (hereafter known as SPoRT-LIS) since 2010, which contains output of soil moisture and temperature at layered depths. The unique configuration of the SPoRT-LIS enables decision-making on operational timescales since it assimilates near real-time observations such as Green Vegetation Fraction from the Visible Infrared Imaging Radiometer Suite and Quantitative Precipitation Estimates from the Multi-Radar Multi-Sensor System. The SPoRT-LIS was developed initially to provide land surface initialization variables for local numerical weather prediction models, such as the Weather Research and Forecasting (WRF) model. The earliest documented uses of the SPoRT-LIS as a tool for local drought analysis came from the National Weather Service Office in Huntsville, AL in 2011, in which data were used to provide recommendations for drought category changes to the U.S. Drought Monitor. Through engagement with SPoRT collaborative partners, use of the SPoRT-LIS has gradually expanded in recent years as a component of drought analysis and feedback to the U.S. Drought Monitor. This effort was aided by feedback from end users who expressed needs for specific soil moisture variables, which led to increased applicability for analysis by others in the drought community. This presentation will provide information on the SPoRT-LIS for drought analysis and the collaborative process that has led to changes in product output to meet user needs, with focus in the Southern Appalachian region. However, applications for hydrology and fire weather will also be discussed.

Land surface modeling↗

CHESS 2025: Discrete-return LiDAR point clouds from NEON AOP surveys

This dataset provides Level 1 (L1) discrete-return light detection and ranging (LiDAR) point cloud data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary unclassified discrete-return LiDAR data delivered by NEON and are provided per flightline as LASzip (LAZ) 1.4 Format 6 files. Data were processed following the workflow described in the NEON L0-to-L1 Discrete Return LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022). Each record in the unclassified point clouds represents a geolocated laser target/return recorded by the LiDAR system, with values for X, Y, Z position and return intensity. All point coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Flight metadata describing flightline boundaries and positional uncertainty by point are also included. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Data and scripts associated with “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA"

This data package is associated with the publication “Allometric scaling of hyporheic respiration across basins in the Pacific Northwest USA” submitted to JGR-Biogeosciences (Regier et al. 2025).This study used reach-scale modeled estimates of hyporheic aerobic respiration made by the River Corridor Model (Fang et al. 2020) and watershed characteristics across the Willamette and Yakima River basins to explore potential allometric scaling (i.e., power-law relationships between size and function) of cumulative hyporheic respiration across catchment-to-basin scales. Scaling was explored quantitatively via the R2, slope, and y-intercept of relationships between cumulative hyporheic respiration and watershed area, divided into hyporheic exchange flux (HEF) quantiles. We also explored relationships between allometric scaling and other watershed characteristics through linear regression, spatial patterns, and mutual information analyses. Our results also suggest variability of hyporheic respiration allometry for middle exchange flux quantiles, and in relation to land-cover. Our findings provide initial evidence that allometric scaling may be useful for predicting hyporheic biogeochemical dynamics across watersheds from reach to basin scales. This data package is associated with the GitHub repository found at https://github.com/peterregier/rc_wrb_yrb_scaling. The data package is organized into several key directories. The “data” folder contains multiple CSV files, including landscape heterogeneity, scaling analysis, and watershed boundary data. The “figures” folder has all figure files in both PDF and PNG formats. Core analysis scripts and figure generation scripts are in the “scripts” directory, systematically numbered for sequential execution. The root directory includes essential project files; please see the file ending in “flmd.csv” for a list and description of all files contained in this data package and the file ending in “dd.csv” for data dictionaries used to describe tabular column headers.

54 ENVIRONMENTAL SCIENCES↗

Leaf 13 C data constrain the uncertainty of the carbon dynamics of temperate forest ecosystems

Stable carbon isotope discrimination occurred in plant biophysical and biogeochemical processes can help understand plant physiology and soil biogeochemistry with respect to carbon cycling. Here, we incorporated the photosynthetic carbon isotope discrimination into a process-based land surface model (iTEM) to test if stable carbon isotope composition (δ 13 C) can impose additional constraint on model parameters. Sequential data assimilation was implemented at six eddy covariance flux tower sites using carbon flux observations including gross primary productivity (GPP) and net ecosystem exchange (NEE) with and without considering foliar δ 13 C (δ 13 C f ) measurement constraints, respectively. Our model-data fusion showed that δ 13 C f can provide useful constraint on photosynthetic (V cmax25 , the maximum rate of carboxylation at 25°C) and stomatal (g1, the slope of stomatal function) parameters as well as the posterior carbon fluxes. When including δ 13 C f measurement, g1 spatially varies among the six sites and is significantly correlated with annual precipitation. We incorporated the statistical relationship between g1 and annual precipitation into iTEM, which is then used to quantify the regional carbon dynamic in temperate forest ecosystems of the Northern Hemisphere. Compared with the simulation only conditioned on carbon flux observations, regional carbon flux estimations performed slightly better against the FLUXCOM products and the uncertainties of modeled carbon fluxes were reduced by 27%. Our study demonstrated that δ 13 C f data constrains carbon flux uncertainties across space.

54 ENVIRONMENTAL SCIENCES↗

Assessing the Impact of Soil Layer Depth Specification on the Observability of Modeled Soil Moisture and Brightness Temperature

The utility of hydrologic land surface models (LSMs) can be enhanced by using information from observational platforms, but mismatches between the two are common.This study assesses the degree to which model agreement with observations (observability) is affected by two mechanisms in particular: 1) physical incongruities between the support volumes being characterized and 2) inadequateor inconsistent parameterizations of physical processes. The Noah and Noah-MP LSMs by default characterize surface soil moisture (SSM) in the top 10 cm of the soil column. This depth may be reasonable when comparing against soil moisture from in situ probes centered at 5 cm,but it is notably different from the 5 cm (or less) sensing depth of NASA’s Soil Moisture Active Passive (SMAP) satellite mission. These depth inconsistencies are examined by using thinner model layers in the Noah and Noah-MP LSMs and comparing resultant simulations to in situ and SMAP soil moisture. In addition, a forward radiative transfer model to simulate microwave brightness temperatures (Tbs) is used to facilitate direct comparisons of LSM-based and SMAP-based L-band Tb retrievals. Observability is quantified using Kolmogorov-Smirnov distance values, calculated from empirical cumulative distribution functions of SSM and Tb time series. Experiment results depend on the particular subspace being analyzed (SSM or Tb). This study concludes that therole of increasedsoil layer discretizationon LSM observability is secondary to the influence of component parameterizations, the effects of which dominate systematic differences with observations

SMAP↗

Analyzing Machine Learning Predictions of Passive Microwave Brightness Temperature Spectral Difference Over Snow-Covered Terrain in High Mountain Asia

Snow is an important component of the terrestrial freshwater budget in high mountainAsia (HMA) and contributes to the runoff in Himalayan rivers through snowmelt. Despitethe importance of snow in HMA, considerable spatiotemporal uncertainty exists across the different estimates of snow water equivalent for this region. In order to better estimate snow water equivalent, radiative transfer models are often used in conjunction with microwave brightness temperature measurements. In this study, the efficacy of support vector machines (SVMs), a machine learning technique, to predict passive microwave brightness temperature spectral difference (1Tb) as a function of geophysical variables (snow water equivalent, snow depth, snow temperature, and snow density) is explored through a sensitivity analysis. The use of machine learning (as opposed to radiative transfer models) is a relatively new and novel approach for improving snow water equivalent estimates. The Noah-MP land surface model within the NASALand Information System framework is used to simulate the hydrologic cycle over HMA and model geophysical variables that are then used for SVM training. The SVMsserve as a nonlinear map between the geophysical space (modeled in Noah-MP) andthe observation space (1Tb as measured by the radiometer). Advanced MicrowaveScanning Radiometer-Earth Observing System measured passive microwave brightness temperatures over snow-covered locations in the HMA region are used as training data during the SVM training phase. Sensitivity of well-trained SVMs to each Noah-MP modeled state variable is assessed by computing normalized sensitivity coefficients. Sensitivity analysis results generally conform with the known first-order physics. Input states that increase volume scattering of microwave radiation, such as snow density and snow water equivalent, exhibit a plurality of positive normalized sensitivity coefficients. In general, snow temperature was the most sensitive input to the SVM predictions. The sensitivity of each state is location and time dependent. The signs of normalized sensitivity coefficients that indicate physical irrationality are ascribed to significant cross-correlation between Noah-MP simulated states and decreased SVM prediction capability at specific locations due to insufficient training data. SVM prediction pitfalls do exist that serve to highlight the limitations of this particular machine learning algorithm.

high mountain Asia↗

Evaluation of a Regional Crop Model Implementation for Sub-National Yield Assessments in Kenya

CONTEXT: Cropping system models can be used to both assess regional food security and to monitor and predict agricultural drought. Agriculture in Kenya is extremely important to both the economy and food security of the country. OBJECTIVE: This study evaluated a regional implementation of a widely used crop model, the Decision Support System for Agrotechnology Transfer (DSSAT), within a coupled modeling framework, the Regional Hydrologic Extremes Assessment System (RHEAS), over Kenya. The goal of this study was to assess the ability of RHEAS to simulate the annual variability of maize yields at the county level and evaluate the uncertainty inherent in the model and inputs. METHODS: The RHEAS system implements a stochastic ensemble approach to account for field scale variabilities in crop management practices and underlying soil and weather conditions. Satellite-derived datasets were used to evaluate the land surface component of the system and seasonally disaggregated yield for 5 years was used to assess the performance of the cropping system model. RESULTS AND CONCLUSIONS: The median correlation between RHEAS and satellite-derived soil moisture and evapotranspiration estimates were 0.78, and 0.51, respectively, indicating that the model is able to capture the key drivers of the hydrological budget. Overall, RHEAS simulated yearly yield variations with a median correlation of 0.7 with reported yields, with the best performance in the short rains season. However, across both seasons, the RHEAS model was positively biased on the order of ~1.6 MT/ha. The overall median unbiased RMSE was 0.66 MT/ha. The RHEAS system shows skill at simulating extreme departures in anomalies, and a majority of the time (62.5%) the reported yields fall within the interquartile range of the simulations. SIGNIFICANCE: One of the most important areas of improvement for the next generation of agricultural data and models is to better understand and communicate the inherent uncertainties. This is especially critical in data-limited regions. Here we present a modeling system and its implementation that begins to address these concerns. We demonstrate the ability to simulate broad trends in yields at the county level for sub-annual yields with skills that commensurate previous national/annual level studies.

Crop model↗

One Stomatal Model to Rule Them All? Toward Improved Representation of Carbon and Water Exchange in Global Models

Abstract Stomatal conductance schemes that optimize with respect to photosynthetic and hydraulic functions have been proposed to address biases in land‐surface model (LSM) simulations during drought. However, systematic evaluations of both optimality‐based and alternative empirical formulations for coupling carbon and water fluxes are lacking. Here, we embed 12 empirical and optimization approaches within a LSM framework. We use theoretical model experiments to explore parameter identifiability and understand how model behaviors differ in response to abiotic changes. We also evaluate the models against leaf‐level observations of gas‐exchange and hydraulic variables, from xeric to wet forest/woody species spanning a mean annual precipitation range of 361–3,286 mm yr −1 . We find that models differ in how easily parameterized they are, due to: (a) poorly constrained optimality criteria (i.e., resulting in multiple solutions), (b) low influence parameters, (c) sensitivities to environmental drivers. In both the idealized experiments and compared to observations, sensitivities to variability in environmental drivers do not agree among models. Marked differences arise in sensitivities to soil moisture (soil water potential) and vapor pressure deficit. For example, stomatal closure rates at high vapor pressure deficit range between −45% and +70% of those observed. Although over half the new generation of stomatal schemes perform to a similar standard compared to observations of leaf‐gas exchange, two models do so through large biases in simulated leaf water potential (up to 11 MPa). Our results provide guidance for LSM development, by highlighting key areas in need for additional experimentation and theory, and by constraining currently viable stomatal hypotheses.

54 ENVIRONMENTAL SCIENCES↗

Differences Between OCO‐2 and GOME‐2 SIF Products From a Model‐Data Fusion Perspective

Space-borne retrievals of solar-induced chlorophyll fluorescence (SIF) over land surfaces have recently become a resource for studying and quantifying the broad scale dynamics of gross carbon uptake (gross primary productivity—GPP) across ecosystems. To prepare for the assimilation of SIF data in terrestrial biosphere models, we examine how differences between SIF products (due to differences in acquisition characteristics and processing chain) may affect the optimization of model parameters and the resultant GPP estimate. We compare recent daily mean SIF products (one from the Orbiting Carbon Observatory-2 [OCO-2] and two from the Global Ozone Monitoring Experiment–2 [GOME-2], GlobFluo [GF] and NASA-v28 [N28], missions), averaged at 0.5° × 0.5° spatial resolution and 16-day temporal resolution, at the biome level. Phase differences between these products are relatively small. A first-order correction of the difference in spectral sampling between the two instruments shows that OCO-2 and N28 are consistent in terms of magnitude and amplitude, while GF is twice as large as the others. Using a bias-blind toy data assimilation framework, we analyze how biases between SIF products, and between model and products, can be partially alleviated by optimizing the slope and intercept parameters of a linear GPP-SIF operator. As observation biases can transfer to biases in other optimized process-based parameters and to modeled carbon fluxes— thereby resulting in unidentified inaccurate parameter values—we argue that potential SIF biases should be treated cautiously in real-world experiments in order to achieve realistic and reliable future simulations.

Gross primary production↗

The Terrestrial Biosphere Model Farm

Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901–2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.

54 ENVIRONMENTAL SCIENCES↗

Influence of Vegetation on Simulation of the Water Balance and Hydrological Response to El Niño–Southern Oscillation in Western Tropical South America

Land surface models (LSMs) rely on vegetation parameters for use in hydrological and energy balance analysis, monitoring, and forecasting. This study examines the influence that vegetation representation in the Noah-Multiparameterization (Noah-MP) LSM has on hydrological simulations across the diverse climate zones of western tropical South America (WTSA), with specific consideration of hydrological variability associated with El Niño–Southern Oscillation (ENSO). The influence of model representation of vegetation on simulated hydrology is evaluated through three simulation experiments that use 1) satellite-derived constant MODIS; 2) satellite-derived time-varying MODIS; and 3) the Noah-MP dynamic leaf model. We find substantial differences in vegetation fields between these simulations, with the Noah-MP dynamic leaf model diverging significantly from satellite-derived vegetation fields in many ecoregions. Impacts on simulated hydrology were, however, found to be modest across climate zones, except for select extreme events. Also, although impacts on hydrology under ENSO-induced variability were small, we find that the Noah-MP dynamic leaf model simulates a positive relationship between rainfall and vegetation in humid ecoregions of WTSA, where satellite observations may indicate the opposite. The relatively small sensitivity of simulated hydrology to vegetation scheme suggests that the performance of hydrological monitoring and forecasting in WTSA that uses Noah-MP is largely unaffected by the choice of vegetation scheme, such that using a simple climatological default is generally no worse than adopting more complicated options. The presence of some differences between the time-varying and constant MODIS simulations for hydrologic extremes, however, indicates that time-varying MODIS configuration might be more suitable for hydrological hazards applications.

South America↗