Search NASA⌕ Search

SEARCH · Search NASA

Results for “evapotranspiration model”

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 487 records · Page 27

Application of a Machine Learning Algorithm in Generating an Evapotranspiration Data Product From Coupled Thermal Infrared and Microwave Satellite Observations

Land surface evapotranspiration (ET) is one of the main energy sources for atmospheric dynamics and a critical component of the local, regional, and global water cycles. Consequently, accurate measurement or estimation of ET is one of the most active topics in hydro-climatology research. With massive and spatially distributed observational data sets of land surface properties and environmental conditions being collected from the ground, airborne or space-borne platforms daily over the past few decades, many research teams have started to use big data science to advance the ET estimation methods. The Geostationary satellite Evapotranspiration and Drought (GET-D) product system was developed at the National Oceanic and Atmospheric Administration (NOAA) in 2016 to generate daily ET and drought maps operationally. The primary inputs of the current GET-D system are the thermal infrared (TIR) observations from NOAA GOES satellite series. Because of the cloud contamination to the TIR observations, the spatial coverage of the daily GET-D ET product has been severely impacted. Based on the most recent advances, we have tested a machine learning algorithm to estimate all-weather land surface temperature (LST) from TIR and microwave (MW) combined satellite observations. With the regression tree machine learning approach, we can combine the high accuracy and high spatial resolution of GOES TIR data with the better spatial coverage of passive microwave observations and LST simulations from a land surface model (LSM). The regression tree model combines the three LST data sources for both clear and cloudy days, which enables the GET-D system to derive an all-weather ET product. This paper reports how the all-weather LST and ET are generated in the upgraded GET-D system and provides an evaluation of these LST and ET estimates with ground measurements. The results demonstrate that the regression tree machine learning method is feasible and effective for generating daily ET under all weather conditions with satisfactory accuracy from the big volume of satellite observations.

54 ENVIRONMENTAL SCIENCES↗

Impact of the ongoing Amazonian deforestation on local precipitation: A GCM simulation study

Numerical simulation experiments were conducted to delineate the influence of in situ deforestation data on episodic rainfall by comparing two ensembles of five 5-day integrations performed with a recent version of the Goddard Laboratory for Atmospheres General Circulation Model (GCM) that has a simple biosphere model (SiB). The first set, called control cases, used the standard SiB vegetation cover (comprising 12 biomes) and assumed a fully forested Amazonia, while the second set, called deforestation cases, distinguished the partially deforested regions of Amazonia as savanna. Except for this difference, all other initial and prescribed boundary conditions were kept identical in both sets of integrations. The differential analyses of these five cases show the following local effects of deforestation. (1) A discernible decrease in evapotranspiration of about 0.80 mm/d (roughly 18%) that is quite robust in the averages for 1-, 2-, and 5-day forecasts. (2) A decrease in precipitation of about 1.18 mm/d (roughly 8%) that begins to emerge even in 1-2 day averages and exhibits complex evolution that extends downstream with the winds. (3) A significant decrease in the surface drag force (as a consequence of reduced surface roughness of deforested regions) that, in turn, affects the dynamical structure of moisture convergence and circulation. The surface winds increase significantly during the first day, and thereafter the increase is well maintained even in the 2- and 5-day averages.

Walker, G. K.↗

An Operational Product for Peatland Applications: The Version 7 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture estimates from April 2015 to present with a mean latency of 2.5 days from the time of observation. The L4_SM estimates are derived from the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment Land Surface Model (CLSM). The recently released Version 7 of the L4_SM product comprises a key advancement that is of interest for peatland applications. CLSM now includes the recently developed PEATCLSM hydrology module for peatlands and uses an updated global peatland distribution. In this presentation, we first give an overview of the operational product for the peatland community. We further show how the incorporation of PEATCLSM considerably improves the dynamics of water table depth, surface soil moisture and evapotranspiration in L4_SM Version 7 over Version 6 when evaluated against in situ measurements in peatlands. These improvements are also manifested in smaller Tb observation-minus-forecast residuals. Eventually, we provide three research examples in which the data assimilation product is used to gain insights into: peatland hydrological processes, the peat moisture dependency of the carbon cycle, and wildfire occurrence in peatlands.

Michel Bechtold↗

Strong influence of convective heat transfer efficiency on the cooling benefits of green roof irrigation

Abstract By enhancing evapotranspiration (ET), green roofs provide cooling benefits for the urban environment and are recognized as a promising heat mitigation strategy. The evaporative cooling effects of green roofs strongly depend on the soil moisture conditions and thus irrigation may be needed to sustain the cooling benefits. It has been shown that the magnitude of cooling benefits offered by green roof irrigation varies spatially, but its controlling factors remain elusive. In this study, we combine a surface energy balance (SEB) model with global simulations generated by an improved Earth System Model to illustrate the key factors controlling the cooling benefits of green roof irrigation. We employ a simple irrigation scheme, which is only active when there is no ice in soil layers and when the soil moisture is below field capacity. As a result, most of the irrigation water leaves the green roof system via ET. We find that the magnitude and also the spatial variability of the cooling benefits of green roof irrigation are controlled by the irrigation amount, and a surface energy redistribution factor that encodes the efficiencies of different SEB components in transferring heat. Further analysis indicates that the enhancement of latent heat flux due to irrigation is largely balanced by the reduction of sensible heat flux on green roofs. Therefore, the amount of irrigation needed per unit decrease of green roof surface temperature is mainly controlled by the convective heat transfer efficiency. A lower convective heat transfer efficiency (e.g. under a lower wind speed) helps reduce the amount of irrigation needed per unit decrease of green roof surface temperature. This study highlights the importance of SEB in constraining the cooling benefits of green roof irrigation and provides valuable guidance for urban planners and policy makers in terms of heat mitigation and sustainable water management.

Wang, Linying (ORCID:0000000308769220)↗

Overland flow numerical model prediction, Lower Triangle Region in East River Watershed, Colorado, 3 days

This data package contains numerical simulation results of surface flow variables such as flow velocity and water depth in Lower Triangle Region in East River Watershed, Colorado. The surface flow is a consequence of a high intensity rainfall event with a total duration of 3 days, available at a resolution of 10 minutes. The results are computed on triangular multiresolution meshes with resolutions ranging from 10 meter to 80 meter. The data package also contains a simulation on a uniform triangular mesh with a resolution of 10 meter. The simulations consider surface flow only and neglect subsurface flow, infiltration, and evapotranspiration. The purpose of the data is to assess the quality of a mesh refinement strategy.

54 ENVIRONMENTAL SCIENCES↗

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Negative Feedback of Vegetation to Greenhouse Warming: A Modeling Approach

Several climate models indicate that in a 2 x CO2 environment, temperature and precipitation would increase and runoff would increase faster than precipitation. These models, however, did not allow the vegetation to increase its leaf density as a response to the physiological effects of increased CO2 and consequent changes in climate. Other assessments included these interactions but did not account for the vegetation down-regulation to reduce plant's photosynthetic activity and as such resulted in a weak vegetation negative response. When we combine these interactions in climate simulations with 2 x CO2, the associated increase in precipitation contributes primarily to increase evapotranspiration rather than surface runoff, consistent with observations, and results in an additional cooling effect not fully accounted for in previous simulations with elevated CO2. By accelerating the water cycle, this feedback slows but does not alleviate the projected warming, reducing the land surface warming by 0.6 C. Compared to previous studies, these results imply that long term negative feedback from CO2-induced increases in vegetation density could reduce temperature following a stabilization of CO2 concentration.

Bounous, L.↗

Examining the peatland shrubification-evapotranspiration feedback following multi-decadal water table manipulation

Northern peatlands are globally important long-term sinks of carbon due to their predominantly saturated conditions. However, these ecosystems are expected to become drier with climate change, potentially leading to shrubification. As such, the response of the shrubification–evapotranspiration (ET) feedback may be of critical importance to future peatland energy, water and carbon dynamics. Here we examined the effect of multi-decadal peatland water table (WT) alteration at three adjacent sites with increasing depth to WT (WET, INTermediate, and DRY). In order to better understand the WT–shrubification–ET feedback, we measured peatland vegetation composition, microtopography and ET partitioning, where ET was measured at the ecosystem, microform, and leaf level using eddy covariance (EC), chambers and porometry, respectively. Averaged across microforms and WT treatments, there was a difference in the median measured leaf resistance (r leaf ) between plant functional types ranging from 213 s m -1 for erect dwarf shrubs, 325 s m -1 for graminoids/sedges, and 520 s m -1 for prostrate dwarf shrubs. Scaled based on LAI, the low r leaf of erect dwarf shrubs dominated hummocks, where sites with a higher proportion of hummocks had lower median canopy resistance (r v ) of 141, 133 and 130 s m -1 at the WET, INT and DRY sites respectively. Nevertheless, ET was highest at the WET site and similar between the INT and DRY sites, with greater evaporation from the moss surface at the WET site. Porometry and EC data along with a three-source model were used to independently assess the evaporative contribution from the moss surface, which ranged from 17% to 40%. For moderate and persistent changes in WT from land-use or climate change, our results suggest vegetation succession is minimal, but the microtopographic development and the concomitant differences in LAI for the various plant functional types is key to understanding changes in total ET and partitioning.

54 ENVIRONMENTAL SCIENCES↗

Energy‐Water Asynchrony Principally Determines Water Available for Runoff From Snowmelt in Continental Montane Forests

ABSTRACT Changes in the volume, rate, and timing of the snowmelt water pulse have profound implications for seasonal soil moisture, evapotranspiration (ET), groundwater recharge, and downstream water availability, especially in the context of climate change. Here, we present an empirical analysis of water available for runoff using five eddy covariance towers located in continental montane forests across a regional gradient of snow depth, precipitation seasonality, and aridity. We specifically investigated how energy‐water asynchrony (i.e., snowmelt timing relative to atmospheric demand), surface water input intensity (rain and snowmelt), and observed winter ET (winter AET) impact multiple water balance metrics that determine water available for runoff (WAfR). Overall, we found that WAfR had the strongest relationship with energy‐water asynchrony (adjusted r 2 = 0.52) and that winter AET was correlated to total water year evapotranspiration but not to other water balance metrics. Stepwise regression analysis demonstrated that none of the tested mechanisms were strongly related to the Budyko‐type runoff anomaly (highest adjusted r 2 = 0.21). We, therefore, conclude that WAfR from continental montane forests is most sensitive to the degree of energy‐water asynchrony that occurs. The results of this empirical study identify the physical mechanisms driving variability of WAfR in continental montane forests and are thus broadly relevant to the hydrologic management and modelling communities.

Webb, Ryan William↗

Daily Precipitation Frequency Distributions Impacts on Land-Surface Simulations of CONUS

Many precipitation-driven data products from land data assimilation systems support assessments of droughts, floods, and other societally-relevant land-surface processes.The accumulated precipitation used as input to these products has a significant impacton water budgets; however, the effects of daily distribution of precipitation on theseproducts are not well known. A comparison of the Integrated Multi-satellite Retrievalsfor GPM (IMERG) and Climate Hazards Group InfraRed Precipitation with Stationsversion 2 (CHIRPS2) rainfall products over the continentalUnited States (CONUS) wasperformed to quantify the impacts of the daily distributionof precipitation on biases anderrors in soil moisture, runoff, and evapotranspiration (ET). Since the total accumulatedprecipitation between the IMERG and CHIRPS product differed, a third precipitationproduct, CHIRPS-to-IMERG (CHtoIM), was produced that usedCHIRPS2 accumulatedprecipitation totals and the daily precipitation frequency distribution of IMERG. Thisnew product supported a controlled analysis of the impact ofprecipitation frequencydistribution on simulated hydrological fields. The CHtoIM had higher occurrences ofprecipitation in the 0–5 mm day−1range, with a lower occurrence of dry days, whichdecreased soil moisture and surface runoff in the land-surface model. The surface soillayer had a tendency to reach saturation more often in the CHIRPS2 simulations, wherethe number of moderate to heavy precipitation days (>5mm day−1) was increased. Usingthe blended CHtoIM product as input reduced errors in surface soil moisture by 5–15%when compared to Soil Moisture Active/Passive (SMAP) data.Similarly, ET errors werealso slightly decreased (∼2%) when compared to SSEBop data. Moderate changes indaily precipitation distributions had a quantifiable impact on soil moisture, runoff, andET. These changes usually improved the model when compared to other modeled andobservational datasets, but the magnitude of the improvements varied by region andtime of year.

Daniel P Sarmiento↗

Extreme hydroclimates amplify the biophysical effects of advanced green-up in temperate China

Vegetation phenology modulates climate by altering energy and water exchange between the land and atmosphere. However, how extreme hydroclimatic conditions modify these phenology-climate feedbacks is still poorly understood. Here, in this study, we used a land–atmosphere-coupled Weather Research and Forecasting model to explore the impacts of advanced green-up on air temperature under different hydroclimate conditions across temperate China and to Mechanistic analysis elucidate the underlying biophysical mechanisms. By imposing a 14-day earlier green-up in line with recent satellite observations, we found that under mean climate conditions, an earlier leaf-out induces immediate surface cooling of 0.14 °C during green-up and a lagging 0.02 °C warming during senescence averaged for temperate China. Extremely humid conditions amplify the cooling effects to 0.18 °C during green-up, extending this cooling into the senescence period. Conversely, under extremely arid conditions, earlier green-up cools air temperature by 0.09 °C, and amplified senescence warming to 0.16 °C. Mechanism analysis revealed that evapotranspiration-dominated non-radiative processes drive immediate cooling during green-up while radiation and circulation process dominates the delayed but opposite warming/cooling effects during senescence in extremely arid/humid hydroclimates. Given the projected continuation of warming trends and increased frequency of extreme climatic events, it is imperative to incorporate the biophysical effects of vegetation phenology into local climate adaptation strategies.

54 ENVIRONMENTAL SCIENCES↗

Flow Directions and Ages of Subsurface Water in a Salt Marsh System Constrained by Isotope Tracing

Salt marshes are dynamic hydrologic systems where terrestrial groundwater, terrestrial surface water, and seawater mix due to bi-directional flows and pressure gradients. Due to the counteracting terrestrial and marine forcings that control these environments, we do not comprehensively understand water fluxes in these complex coastal systems. To understand the water sources, flow directions, and velocities in salt marsh porewater, we employed a combination of geochemical tracers and analytical models across a hillslope-to-salt marsh continuum in a salt marsh experiencing daily inundation of estuarine surface water (SW) from tides and mixing of fresh seasonal groundwater. We used tritium ( 3 H) as a hydrologic tracer to assess porewater ages and stable water isotope ($\delta$ 2 H and $\delta$ 18 O) analyses to separate isotopically distinct estuarine and terrestrial groundwater across different depths and landscape positions in the study transect. We employed electrical conductivity to constrain the role of source mixing and evapotranspiration in salt marsh hydrology. Salinity and stable isotopes revealed that transpiration, rather than evaporation, increased subsurface water salinity to concentrations above estuarine SW during summer. Elevated salinity at depth indicated that salt marsh subsurface water is recharged during the dry growing season. Seasonal recharge patterns drive long-term deep subsurface water dynamics across the salt marsh, with 3 H ages of 3–7 years, and daily tidal cycles drive short-term shallow porewater dynamics with 3 H ages of 0 ± 3.6 years. Furthermore, our conceptual understanding of the spatiotemporal changes in SW-subsurface water interactions at the terrestrial-marine interface quantifies the hydrological constraints we are missing to improve our understanding of biogeochemical cycles within the salt marsh.

54 ENVIRONMENTAL SCIENCES↗

Improvements in lake volume predictions using Landsat data

A cumulative error in the water balance budget for Lake Okeechobee produces a one million acre-foot discrepancy in the predicted water volume over a 4-year period. The major source of error appears to be complex shoreline marshes that comprise 20 percent of the lake surface. The water balance budget model presently treats these marshes as open water. Using Landsat data, the vegetation in the lake's littoral zone was classified multispectrally to provide a data base for determining water budget information. First, the acreage of a given plant species in the littoral zone was obtained with satellite data. Second, the surface area occupied by plants (which therefore could not be considered open water) was used to adjust the vegetation acreage giving an effective water surface. Based on this information, more detailed representations of evapotranspiration and total water surface (and hence total lake volume) could be provided to the water balance budget computation.

Gervin, J. C.↗

Geocoding and stereo display of tropical forest multisensor datasets

Concern about the future of tropical forests has led to a demand for geocoded multisensor databases that can be used to assess forest structure, deforestation, thermal response, evapotranspiration, and other parameters linked to climate change. In response to studies being conducted at the Braulino Carrillo National Park, Costa Rica, digital satellite and aircraft images recorded by Landsat TM, SPOT HRV, Thermal Infrared Multispectral Scanner, and Calibrated Airborne Multispectral Scanner sensors were placed in register using the Landsat TM image as the reference map. Despite problems caused by relief, multitemporal datasets, and geometric distortions in the aircraft images, registration was accomplished to within + or - 20 m (+ or - 1 data pixel). A digital elevation model constructed from a multisensor Landsat TM/SPOT stereopair proved useful for generating perspective views of the rugged, forested terrain.

Welch, R.↗

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗

Quantifying the long-term changes of terrestrial water storage and their driving factors

Global warming is expected to cause changes in terrestrial water storage (TWS) across the land surface, with widespread impacts on ecosystems and society. Although extensive research has been performed to analyze TWS changes and possible drivers during the post-2000 period, longer-term evolution of TWS and associated environmental forcings remain relatively unexplored. In this study, we evaluated the performance of the Energy Exascale Earth System model (E3SM) land model ELM version 1 (ELM v1) in simulating global TWS, and used factorial simulations of ELMv1 to quantify global TWS changes and their drivers during 1948–2012. We found that ELM’s agreed best with existing satellites and reconstruction datasets in temperate regions unaffected by irrigation. Biome- and climate zone-averaged TWS mainly increased at rates between 0 and 10 mm/year over 1948–2012, but the second half of that period saw smaller positive trends than the first half or even negative trends. Climate change explained >80 % of the TWS trends across most biomes and climate zones, followed by land use and land cover change. The physiological and phenological effects of CO2 primarily induced noticeable TWS trends in the more humid biomes and climate zones across different latitudes. In contrast, nitrogen deposition and aerosol deposition generally had smaller and negative impacts across the biomes and climate regions. Among the meteorological drivers analyzed, the long-term average imbalance between precipitation (P), evapotranspiration (E), and runoff (Q) contributed >50 % of the TWS trends in most biomes and climate zones, with nonlinearity being induced by spatially heterogenous changes in E/P and Q/P ratios. The accumulated detrended anomalies in P, E, and Q also often contributed substantially, while the trends difference between P, E, and Q contributed little. Together, these findings unveiled an intensification of the global TWS and its diverse patterns of climate change and different non-withdrawal human-induced alterations, contributing to a more comprehensive understanding and projection of the global water cycle.

54 ENVIRONMENTAL SCIENCES↗

A Hybrid Biophysical‐Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy ( LE ) and sensible heat ( H ) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R 2 = 0.81–0.94) and H (R 2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

evapotranspiration↗

NGEE Arctic CO2, CH4 and Energy Eddy-Covariance (EC) Flux Tower Auxiliary Measurements, Utqiagvik (Barrow), Alaska, 2012 - 2019

Measurements began in September 2012 and are generally made throughout the snow-free growing season starting from April (pre-thawing) until November (freeze-up) when the instrumentation is taken down for maintenance and instrument calibration. CO2, CH4 and energy fluxes using the Eddy Covariance (EC) technique (Baldocchi 2003) and meteorological measurements are reported as 30-minute averages. The tower site is registered with AmeriFlux as US-NGB. See http://ameriflux.lbl.gov/sites/siteinfo/US-NGB for more information. Reported data include: (1) additional/extra meteorological and eddy covariance CO2, CH4 and energy flux data, (2) binned spectra for the three wind components, the sonic temperature and gas densities together with the binned cospectra for covariances of w (vertical wind component) and gas densities, (3) binned ogives (cumulative (co)spectra) for the three wind components, the sonic temperature and gas densities, and (4) tower site footprint matrices. Meteorological and eddy covariance CO2, CH4 and energy flux data (AmeriFlux BASE data product) are available from AmeriFlux https://doi.org/10.17190/AMF/1436326. The tower site is adjacent to the NGEE Arctic Tram system with several complementary time-series measurements https://doi.org/10.5440/1558782.This dataset supersedes the previously released dataset "Eddy-Covariance and auxiliary measurements, NGEE-Barrow, 2012-2013" DOI:10.5440/1124200 as this new product includes data for the longer time period processed using consistent methods for the entire record. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) near Utqiagvik (Barrow), Alaska and 2) multiple areas in the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗