Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Land Surface Modeling”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

BioRT‐HBV 1.0: A Biogeochemical Reactive Transport Model at the Watershed Scale

Abstract Reactive Transport Models (RTMs) are essential tools for understanding and predicting intertwined ecohydrological and biogeochemical processes on land and in rivers. While traditional RTMs have focused primarily on subsurface processes, recent watershed‐scale RTMs have integrated ecohydrological and biogeochemical interactions between surface and subsurface. These emergent, watershed‐scale RTMs are often spatially explicit and require extensive data, computational power, and computational expertise. There is however a pressing need to create parsimonious models that require minimal data and are accessible to scientists with limited computational background. To that end, we have developed BioRT‐HBV 1.0, a watershed‐scale, hydro‐biogeochemical RTM that builds upon the widely used, bucket‐type HBV model known for its simplicity and minimal data requirements. BioRT‐HBV uses the conceptual structure and hydrology output of HBV to simulate processes including advective solute transport and biogeochemical reactions that depend on reaction thermodynamics and kinetics. These reactions include, for example, chemical weathering, soil respiration, and nutrient transformation. The model uses time series of weather (air temperature, precipitation, and potential evapotranspiration) and initial biogeochemical conditions of subsurface water, soils, and rocks as input, and output times series of reaction rates and solute concentrations in subsurface waters and rivers. This paper presents the model structure and governing equations and demonstrates its utility with examples simulating carbon and nitrogen processes in a headwater catchment. As shown in the examples, BioRT‐HBV can be used to illuminate the dynamics of biogeochemical reactions in the invisible, arduous‐to‐measure subsurface, and their influence on the observed stream or river chemistry and solute export. With its parsimonious structure and easy‐to‐use graphical user interface, BioRT‐HBV can be a useful research tool for users without in‐depth computational training. It can additionally serve as an educational tool that promotes pollination of ideas across disciplines and foster a diverse, equal, and inclusive user community.

Sadayappan, Kayalvizhi↗
DOE OSTI2025

Data and code for: Combining eddy covariance towers, field measurements, and the MEMS 2 ecosystem model improves confidence in the climate impacts of bioenergy with carbon capture and storage

BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty.

BECCS↗
NASA NTRS2013

Visualization and Analysis of Multi-scale Land Surface Products via Giovanni Portals

Large volumes of MODIS land data products at multiple spatial resolutions have been integrated into the Giovanni online analysis system to support studies on land cover and land use changes,focused on the Northern Eurasia and Monsoon Asia regions through the LCLUC program. Giovanni (Goddard Interactive Online Visualization ANd aNalysis Infrastructure) is a Web-based application developed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC), providing a simple and intuitive way to visualize, analyze, and access Earth science remotely-sensed and modeled data.Customized Giovanni Web portals (Giovanni-NEESPI andGiovanni-MAIRS) have been created to integrate land, atmospheric,cryospheric, and societal products, enabling researchers to do quick exploration and basic analyses of land surface changes, and their relationships to climate, at global and regional scales. This presentation shows a sample Giovanni portal page, lists selected data products in the system, and illustrates potential analyses with imagesand time-series at global and regional scales, focusing on climatology and anomaly analysis. More information is available at the GES DISCMAIRS data support project portal: http:disc.sci.gsfc.nasa.govmairs.

Shen, Suhung↗
NASA NTRSDate not supplied

Improved Treatment of Snow Over Ice Sheets in the NASA GISS Climate Model: Towards Ice Sheet–Climate Coupling

Representing the interactions between ice sheets and climate is essential for more accurate prediction of climate change and sea level rise. Ice sheets interact with the overlying atmosphere via the accumulation of snow and its compaction into firn, then ice, as well as the melting of surface snow and ice and the creation of runoff water. Getting an adequate representation of heat transfer, compaction, and melting processes is essential for an accurate representation of snow on land ice in global climate models. We are implementing an improved snow model on top of land ice as part of an effort to couple the NASA GISS climate model with the PISM ice sheet model. The new snow model includes additional layers and processes that are not currently incorporated (e.g., liquid water retention, percolation and refreezing, and snow densification), and mass and energy transfer methods that are consistent with both static ice sheets (with implicit iceberg fluxes) and interactive ice sheets (with explicit dynamics). We are tuning the densification scheme of this snow model with temperature and density data from common FirnCover and SumUp observations at locations in the accumulation zone of Greenland, and we compare the resulting density profiles to other SumUp density profiles in Greenland and Antarctica. We will assess the impact of this new snow model in climate model simulations with a static ice sheet compared with the previous (simpler) 2-layer snow model. Finally, we aim to use the non-coupled simulations as a baseline to assess the impact of dynamic coupling with an interactive ice sheet model.

snow↗
NASA NTRS1977

Improved scene generator capability

Components of an integrated system having all functions for model preparation and generation, image computation and scene generation, and collision detection and reporting for the visual simulation of space shuttle maneuvers and its interaction with the manipular arm and other objects, as well as during earth landing exercises are described. High surface solution by range priority allows for correct views of dynamically moving objects from any viewpoint. Smooth shading, anti-aliasing and edge smoothing techniques provide a high quality picture. The collision detection system allows for simulation of impenetrable objects and docking maneuvers, and the modeling systems permit production of new data bases for desired simulations.

Source record↗
NASA NTRSDate not supplied

Updating Lunar Subsurface Geotechnical Properties Based Upon Nonlinear Winkler Modeling of Apollo and Lunokhod Measurements

Returning human and robotic explorers to the lunar surface has identified a need for a renewed understanding of lunar geotechnical properties related to landing on and exploring the lunar surface. We review fundamental lunar geotechnical properties such as bulk density, cohesion, friction angle and shear strength, and apply new data derived from Apollo and Lunokhod geotechnical measurements to update the understanding of lunar regolith behavior. The analysis derives a new, nonlinear relationship of lunar geotechnical properties as a function of depth.

Lunar Geotechnical engineering↗
NASA NTRSDate not supplied

Advances in Stereo Photogrammetry Capabilities for in-Situ Ejecta Measurements in Lunar Plume-Surface Interactions

Plume-surface interaction is a critical topic of interest as human spaceflight programs establish a renewed presence on the Moon and push onwards to Mars. The erosion of regolith and high-speed ejection of soil and rocks which occurs during a powered spacecraft landing on lunar or terrestrial bodies can pose serious risks to the vehicle during landing and to existing ground assets. Ground testing campaigns have provided data sets in relevant environments (e.g., low pressure) in tandem with predictive computational model validation. To-date, limited in-situ flight data exists to inform the study of the key phenomena of plume-surface interaction: plume-flow physics, soil erosion, and ejecta dynamics. In two upcoming lunar missions, erosion of the lunar regolith will be measured in-situ using the stereo photogrammetry technique as part of NASA’s Commercial Lunar Payload Services (CLPS) Program. The Stereo Cameras for Lunar Plume Surface Studies (SCALPSS) will use a multiple-camera system to measure the erosion of the lunar regolith due to the vehicles’ rocket engine plumes during descent and after landing. The application of stereo photogrammetry to accurately measure a realistic eroded shape has been validated in laboratory work ahead of the scheduled flights. In preparation for future CLPS, Artemis, or other lunar science missions, recent advancements have been made to measure the ejecta structure and particle dynamics in addition to the erosion of regolith using a SCALPSS-like photogrammetry system. The thickness and ejection angle of entrained soil and dust produced during plume-induced erosion can be used to determine how far debris will travel, and thus how it may impact nearby infrastructure or the landing vehicle itself. The present work highlights the development of this new flight instrumentation capability and its value to advancements in the plume-surface interaction measurement topic for lunar or planetary landings.

Olivia K. Tyrrell↗
NASA NTRSDate not supplied

Geologic Maps for Civil Engineering and Infrastructure Planning

The Artemis Program mandate for sustainable, long-term presence on the Moon requires learning to "Live Off the Land". Geologic maps enable a paradigm shift to model the development and scaling of infrastructure for robotic and human exploration inside a bounded system on the lunar surface. Relying on geologic maps produced from orbital data and robotic precursor missions, in-situ resource utilization (ISRU) will be required to produce infrastructure, building materials and consumables. The ISRU value chain must be mapped to identify linkages of resources, producers, processes and equipment capabilities for transforming excavated regolith into usable feedstocks. Coordination will be required to balance supply and demand across the whole system of lunar infrastructure and operations, at each stage of development.

Geologic Mapping↗
NASA NTRS2011

Algorithm for Autonomous Landing

Because of their small size, high maneuverability, and easy deployment, micro aerial vehicles (MAVs) are used for a wide variety of both civilian and military missions. One of their current drawbacks is the vast array of sensors (such as GPS, altimeter, radar, and the like) required to make a landing. Due to the MAV s small payload size, this is a major concern. Replacing the imaging sensors with a single monocular camera is sufficient to land a MAV. By applying optical flow algorithms to images obtained from the camera, time-to-collision can be measured. This is a measurement of position and velocity (but not of absolute distance), and can avoid obstacles as well as facilitate a landing on a flat surface given a set of initial conditions. The key to this approach is to calculate time-to-collision based on some image on the ground. By holding the angular velocity constant, horizontal speed decreases linearly with the height, resulting in a smooth landing. Mathematical proofs show that even with actuator saturation or modeling/ measurement uncertainties, MAVs can land safely. Landings of this nature may have a higher velocity than is desirable, but this can be compensated for by a cushioning or dampening system, or by using a system of legs to grab onto a surface. Such a monocular camera system can increase vehicle payload size (or correspondingly reduce vehicle size), increase speed of descent, and guarantee a safe landing by directly correlating speed to height from the ground.

Kuwata, Yoshiaki↗
DOE OSTI2020

Fire, dust, air and water: Improving aerosol biogeochemistry interactions in ACME. Final report

Aerosol-biogeochemistry interactions will have important consequences for human health, ecosystem services, and climate feedbacks over the next several decades. In this proposal we seek to link how human activities on land, including direct industrial emissions and landscape fires, impact the atmosphere and the biogeochemistry of the ocean and land. Our overall goal was to complement existing work at DOE labs to improve the capability and performance of the E3SM model in short-term, high-resolution projections of climate. Projections of ocean uptake of anthropogenic carbon are highly dependent on the details of the ocean physical and biogeochemistry parameterizations. Potentially equally important are the estimates of increases in iron and soluble nitrogen deposition in aerosols to the open ocean, but there are substantial uncertainties in these feedbacks. The E3SM land model simulates natural and anthropogenic fires, but improvement is needed to more accurately simulate fire responses to climate variability as well as fire contributions to terrestrial climate feedbacks. This proposal addressed four of the six priorities within the ACME/E3SM call by improving the representation of fire, industrial combustion sources, and desert dust and adding iron and soluble iron to the atmospheric chemistry (1). Our work focused on improving ocean biogeochemistry and the response of the oceans to anthropogenic atmospheric iron and nitrogen (3). We also improved the representation of land disturbance on fires and fire mediated carbon cycle feedbacks in the land model (4). Finally we had planned to couple the land, atmosphere and ocean biogeochemical systems (6) and explore how this coupling changes the feedbacks in the climate system. Our work under this proposal was separated into the following 5 tasks: i) improve representation of fire dynamics in both boreal and tropical ecosystems, ii) improve estimates of industrial emissions of iron and phosphorus, iii) improve representation of atmospheric iron, phosphorus and nitrogen transport, chemistry and deposition, iv) improve representation of iron, phosphorus and nitrogen deposition impacts on ocean biogeochemistry, and v) couple emissions into atmosphere, to deposition onto land and ocean biogeochemistry. For all 5 tasks, we conducted offline experiments of the historical period to enable comparison with available aerosol and land surface observations. We had planned to conduct fully coupled simulations, however, these will be conducted as part of the coupling in version 3 and 4 of the E3SM model, and thus have not been completed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Description of historical and future projection simulations by the global coupled E3SMv1.0 model as used in CMIP6

Abstract. This paper documents the experimental setup and general features of the coupled historical and future climate simulations with the first version of the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv1.0). The future projected climate characteristics of E3SMv1.0 at the highest emission scenario (SSP5-8.5) designed in the Scenario Model Intercomparison Project (ScenarioMIP) and the SSP5-8.5 greenhouse gas (GHG) only forcing experiment are analyzed with a focus on regional responses of atmosphere, ocean, sea ice, and land. Due to its high equilibrium climate sensitivity (ECS of 5.3 K), E3SMv1.0 is one of the Coupled Model Intercomparison Project phase 6 (CMIP6) models with the largest surface warming by the end of the 21st century under the high-emission SSP5-8.5 scenario. The global mean precipitation change is highly correlated with the global temperature change, while the spatial pattern of the change in runoff is consistent with the precipitation changes. The oceanic mixed layer generally shoals throughout the global ocean. The annual mean Atlantic meridional overturning circulation (AMOC) is overly weak with a slower change from ∼ 11 to ∼ 6 Sv (Sverdrup) relative to other CMIP6 models. The sea ice, especially in the Northern Hemisphere, decreases rapidly with large seasonal variability. We detect a significant polar amplification in E3SMv1.0 from the atmosphere, ocean, and sea ice. Comparing the SSP5-8.5 all-forcing experiment with the GHG-only experiment, we find that the unmasking of the aerosol effects due to the decline of the aerosol loading in the future projection period causes transient accelerated warming in the all-forcing experiment in the first half of the 21st century. While the oceanic climate response is mainly controlled by the GHG forcing, the land runoff response is impacted primarily by forcings other than GHG over certain regions, e.g., southern North America, southern Africa, central Africa, and eastern Asia. However, the importance of the GHG forcing on the land runoff changes grows in the future climate projection period compared to the historical period.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗
NASA NTRS2006

Aerosol-Water Cycle Interaction: A New Challenge in Monsoon Climate Research

Long recognized as a major environmental hazard, aerosol is now known to have strong impacts on both regional and global climate. It has been estimated that aerosol may reduce by up to 10% of the seasonal mean solar radiation reaching the earth surface, producing a global cooling effect that opposes global warming (Climate Change 2001). This means that the potential perils that humans have committed to global warming may be far greater than what we can detect at the present. As a key component of the Earth climate system, the water cycle is profoundly affected by the presence of aerosols in the atmosphere. Through the so-called direct effect , aerosol scatters and/or absorbs solar radiation, thus cooling the earth surface and changing the horizontal and vertical radiational heating contrast in the atmosphere. The heating contrast drives anomalous atmospheric circulation, resulting in changes in convection, clouds, and rainfall. Another way aerosol can affect the water cycle is through the so-called indirect effects, whereby aerosol increases the number of cloud condensation nuclei, prolongs life time of clouds, and inhibits the growth of cloud drops to raindrops. This leads to more clouds, and increased reflection of solar radiation, and further cooling at the earth surface. In monsoon regions, the response of the water cycle to aerosol forcing is especially complex, not only because of presence of diverse mix of aerosol species with vastly different radiative properties, but also because the monsoon is strongly influenced by ocean and land surface processes, land use, land change, as well as regional and global greenhouse warming effects. Thus, sorting out the impacts of aerosol forcing, and interaction with the monsoon water cycle is a very challenging problem. In this talk, I will offer some insights into how aerosols may impact the Asian monsoon based on preliminary results from satellite observations and climate model experiments. Specifically, I will discuss the elevated heat pump hypothesis, involving atmospheric heating by absorbing aerosols (dust and black carbon) over the southern slopes of the Himalayas, and feedback with the deep convection, in modifying monsoon water cycle over South .and East Asia. The role of aerosol forcing relative to those due to sea surface temperature and land surface processes, as well as observation requirements to verify such a hypothesis will also be discussed.

Lau, William K. M.↗
NASA NTRS2006

Aerosol-Water Cycle Interaction: A New Challenge in Monsoon Climate Research

Long recognized as a major environmental hazard, aerosol is now known to have strong impacts on both regional and global climate. It has been estimated that aerosol may reduce by up to 10% of the seasonal mean solar radiation reaching the earth surface, producing a global cooling effect that opposes global warming (Climate Change 2001). This means that the potential perils that humans have committed to global warming may be far greater than what we can detect at the present. As a key component of the Earth climate system, the water cycle is profoundly affected by the presence of aerosols in the atmosphere. Through the so-called "direct effect", aerosol scatters and/or absorbs solar radiation, thus cooling the earth surface and changing the horizontal and vertical radiational heating contrast in the atmosphere. The heating contrast drives anomalous atmospheric circulation, resulting in changes in convection, clouds, and rainfall. Another way aerosol can affect the water cycle is through the so-called "indirect effects", whereby aerosol increases the number of cloud condensation nuclei, prolongs life time of clouds, and inhibits the growth of cloud drops to raindrops. This leads to more clouds, and increased reflection of solar radiation, and further cooling at the earth surface. In monsoon regions, the response of the water cycle to aerosol forcing is especially complex, not only because of presence of diverse mix of aerosol species with vastly different radiative properties, but also because the monsoon is strongly influenced by ocean and land surface processes, land use, land change, as well as regional and global greenhouse warming effects. Thus, sorting out the impacts of aerosol forcing, and interaction with the monsoon water cycle is a very challenging problem. In this talk, I will offer some insights into how aerosols may impact the Asian monsoon based on preliminary results from satellite observations and climate model experiments. Specifically, I will discuss the "elevated heat pump" hypothesis, involving atmospheric heating by absorbing aerosols (dust and black carbon) over the southern slopes of the Himalayas, and feedback with the deep convection, in modifying monsoon water cycle over South and East Asia. The role of aerosol forcing relative to those due to sea surface temperature and land surface processes, as well as observation requirements to verify such a hypothesis will also be discussed.

Lau, William K. M.↗
NASA NTRS2002

Impact of Soil Moisture Initialization on Seasonal Weather Prediction

The potential role of soil moisture initialization in seasonal forecasting is illustrated through ensembles of simulations with the NASA Seasonal-to-Interannual Prediction Project (NSIPP) model. For each boreal summer during 1997-2001, we generated two 16-member ensembles of 3-month simulations. The first, "AMIP-style" ensemble establishes the degree to which a perfect prediction of SSTs would contribute to the seasonal prediction of precipitation and temperature over continents. The second ensemble is identical to the first, except that the land surface is also initialized with "realistic" soil moisture contents through the continuous prior application (within GCM simulations leading up to the start of the forecast period) of a daily observational precipitation data set and the associated avoidance of model drift through the scaling of all surface prognostic variables. A comparison of the two ensembles shows that soil moisture initialization has a statistically significant impact on summertime precipitation and temperature over only a handful of continental regions. These regions agree, to first order, with regions that satisfy three conditions: (1) a tendency toward large initial soil moisture anomalies, (2) a strong sensitivity of evaporation to soil moisture, and (3) a strong sensitivity of precipitation to evaporation. The degree to which the initialization improves forecasts relative to observations is mixed, reflecting a critical need for the continued development of model parameterizations and data analysis strategies.

Koster, Randal D.↗
DOE OSTI2025

Assessing methods in fusion and fitting for time series construction in remote sensing-based earth observations

This study evaluates the comparative performance of spatiotemporal fusion and time-series fitting methods for constructing high-spatiotemporal-resolution remote sensing time-series data. Due to in-class similarity of fusion methods and fitting methods, we employ the Fit-FC (Fitting, spatial Filtering, and residual Compensation) model as a representative fusion method and the linear harmonic fitting model as a representative fitting method. Both Fit-FC and the linear harmonic fitting are widely used for high-spatiotemporal-resolution time-series data construction, and we modify the original Fit-FC model to enable automatic time-series fusion. To ensure data representativeness, we use 3 years (2019–2021) of Harmonized Landsat and Sentinel-2 surface reflectance datasets and Terra MCD43A4 products. Eight experimental regions are selected worldwide to guarantee generalization of the comparative performance between fusion and fitting methods, covering diverse land-use types (cropland, developed land, forest, and grassland) and varying climatological conditions. Time-series of NDVI and surface reflectance are analyzed under both actual observations and simulated data-missing scenarios. The constructed time-series data reveals that (1) the modified Fit-FC and linear harmonic fitting model achieve excellent performance in constructing high-resolution time-series images; (2) the fusion method outperforms the fitting method in constructing time-series of NDVI and surface reflectance images in cropland-, forest-, and grassland-dominated regions; (3) both methods achieve comparable performance in developed-dominated regions; (4) the fusion method is more robust to missing data, and better captures abrupt phenological transitions under conditions of continuous missing data; (5) the fitting method is computationally more efficient, making it suitable for large-scale time-series image reconstruction. This study provides valuable insights for selecting optimal strategies to generate high-resolution time-series images across diverse application scenarios and lays a foundation for extensions to other vegetation indices or land surface variables.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2008

Lagrangian Trajectory Modeling of Lunar Dust Particles

Apollo landing videos shot from inside the right LEM window, provide a quantitative measure of the characteristics and dynamics of the ejecta spray of lunar regolith particles beneath the Lander during the final 10 [m] or so of descent. Photogrammetry analysis gives an estimate of the thickness of the dust layer and angle of trajectory. In addition, Apollo landing video analysis divulges valuable information on the regolith ejecta interactions with lunar surface topography. For example, dense dust streaks are seen to originate at the outer rims of craters within a critical radius of the Lander during descent. The primary intent of this work was to develop a mathematical model and software implementation for the trajectory simulation of lunar dust particles acted on by gas jets originating from the nozzle of a lunar Lander, where the particle sizes typically range from 10 micron to 500 micron. The high temperature, supersonic jet of gas that is exhausted from a rocket engine can propel dust, soil, gravel, as well as small rocks to high velocities. The lunar vacuum allows ejected particles to travel great distances unimpeded, and in the case of smaller particles, escape velocities may be reached. The particle size distributions and kinetic energies of ejected particles can lead to damage to the landing spacecraft or to other hardware that has previously been deployed in the vicinity. Thus the primary motivation behind this work is to seek a better understanding for the purpose of modeling and predicting the behavior of regolith dust particle trajectories during powered rocket descent and ascent.

Lane, John E.↗
DOE OSTI2022

Estimation of the Surface Fluxes for Heat and Momentum in Unstable Conditions with Machine Learning and Similarity Approaches for the LAFE Data Set

Abstract Measurements of three flux towers operated during the land atmosphere feedback experiment (LAFE) are used to investigate relationships between surface fluxes and variables of the land–atmosphere system. We study these relations by means of two machine learning (ML) techniques: multilayer perceptrons (MLP) and extreme gradient boosting (XGB). We compare their flux derivation performance with Monin–Obukhov similarity theory (MOST) and a similarity relationship using the bulk Richardson number (BRN). The ML approaches outperform MOST and BRN. Best agreement with the observations is achieved for the friction velocity. For the sensible heat flux and even more so for the latent heat flux, MOST and BRN deviate from the observations while MLP and XGB yield more accurate predictions. Using MOST and BRN for latent heat flux, the root mean square errors (RMSE) are 107 Wm $$^{-2}$$ - 2 and 121 Wm $$^{-2}$$ - 2 , respectively, as well as the intercepts of the regression lines are $$\approx 110$$ ≈ 110 Wm $$^{-2}$$ - 2 . For the ML methods, the RMSEs reduce to 31 Wm $$^{-2}$$ - 2 for MLP and 33 Wm $$^{-2}$$ - 2 for XGB as well as the intercepts to just 4 Wm $$^{-2}$$ - 2 for MLP and $$-1$$ - 1 Wm $$^{-2}$$ - 2 for XGB with slopes of the regression lines close to 1, respectively. These results indicate significant deficiencies of MOST and BRN, particularly for the derivation of the latent heat flux. In fact, in contrast to the established theories, feature importance weighting demonstrates that the ML methods base their improved derivations on net radiation, the incoming and outgoing shortwave radiations, the air temperature gradient, and the available water contents, but not on the water vapor gradient. The results imply that further studies of surface fluxes and other turbulent variables with ML techniques provide great promise for deriving advanced flux parameterizations and their implementation in land–atmosphere system models.

54 ENVIRONMENTAL SCIENCES↗