SEARCH · Engineering Papers

Results for “atmospheric 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Decadal Hybrid GCM Simulation Using Deep‐Learning‐Based Cloud and Convection Parameterization Generalized to a Warm Climate

A critical challenge for machine‐learning (ML) parameterization in global climate models (GCMs) is to achieve stable, accurate simulations under climates not seen during training. Previous studies have demonstrated promising offline performance and year‐long online stability in aquaplanet simulations but have encountered difficulties in real geography and under climate warming. Here we report that a GCM with real geography configuration using neural‐network‐based cloud and convection parameterization, trained exclusively with present‐day climate data, successfully performs a stable, decade‐long simulation of a warm climate with +4 K sea surface temperature (SST). The neural network (NN) is based on Han et al. (2023, https://doi.org/10.1029/2022ms003508 ) with additional inputs. The simulation captures the global precipitation distribution, surface temperatures, vertical atmospheric structures, and extreme precipitation very well, closely matching simulations from both the superparameterized CAM (SPCAM) and the conventional CAM5 in the warm climate without accuracy degradation compared to those in the baseline climate. Moreover, it produces a climate response to +4 K SST in atmospheric thermodynamic states and circulations similar to those from SPCAM and CAM5. Prognostic ablation tests on NN input variables show that the NN without convective memory as input suffers from numerical instability, and the NN without considering radiative variables and land fraction as input, or with reduced training samples produce less accurate results. To our knowledge, this is the first time an ML parameterization successfully achieves online extrapolation to a warm climate without using additional warm‐climate data for training. It demonstrates the potential of ML‐driven parameterizations for credible long‐term climate projections.

Atmosphere model↗

DOE OSTI2021

Characteristics of thermospheric infrasound predicted using ray tracing and weakly non-linear waveform analyses

Infrasonic signals refracted by thermal gradients in the rarefied upper atmosphere are modeled using a combination of ray tracing and weak shock theory to develop an understanding of thermospheric infrasound signals produced by energetic, transient sources. Canonical arrival structures in the form of u-wave signatures are identified for returns refracted at lower altitudes within the thermosphere, and possible multi-pathing produced by effective sound speed inflections are investigated to elucidate more complex arrival structures, which are found to be spatially localized. Furthermore, variability in the source characteristics is investigated and it is found that whereas some waveform phase information is lost due to finite amplitude effects, arrival characteristics are strongly dependent on the peak overpressure near the source. Variability in the propagation path is considered using archived atmospheric specifications and implies that despite uncertainties related to the dynamic and sparsely sampled nature of the atmosphere, thermospheric signatures might be useful in estimating the yield for explosive sources. Last, thermospheric arrivals from a failed rocket launch, as well as those from several large chemical explosions, are analyzed and it is found that qualitative trends match those predicted, and analyses here provide additional insight into such signatures.

42 ENGINEERING↗

DOE OSTI2021

Calculated Roughness Changes due to Clear-Cutting At Savannah River Site

At the Savannah River Site (SRS), multiple facilities are routinely producing and working with radiological constituents which require safety basis modeling by the National Nuclear Safety Administration. The models used to determine the safety basis employ roughness length as an input parameter which describes the roughness of the ground surface which influences momentum exchange near the surface and the subsequent generation of mechanical turbulence. Because roughness length is dependent on surface characteristics, including vegetation, any changes to the surface would have a subsequent impact on the roughness length and require a new safety basis analysis to be performed. An analysis of measurements prior to and after the dear-cuts and thinning was performed and roughness lengths were calculated for each of sixteen wind sectors. Not all sectors contained dear-cuts, so this provides comparison of natural variations in roughness. Based on the results of the estimated changes, consideration will be given to how these changes may affect atmospheric dispersion modeling. It was found that clear-cutting activities generally led to a reduction in roughness. Reductions in roughness length also occurred in some of the wind sectors containing only undisturbed forest, but those in sectors with dear-cuts were found to be slightly greater in magnitude than natural variation. Only in one instance was the change in roughness length found to exceed a decrease of 10%, indicating that there is no significant safety basis to the site as a whole based on previous clear-cutting activities.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗

DOE OSTI2021

Refined Telluric Absorption Correction for Low-resolution Ground-based Spectroscopy: Resolution and Radial Velocity Effects in the O{sub 2}A-band for Exoplanets and K i Emission Lines

Telluric correction of spectroscopic observations is either performed via standard stars that are observed close in time and airmass along with the science target, or recently growing in importance, by theoretical telluric absorption modeling. Both approaches work fine when the telluric lines are resolved, i.e., at a spectral resolving power larger than about 10,000, and it is sufficient to facilitate the detection of spectral features at lower resolution. However, a meaningful quantitative analysis also requires the reliable recovery of line strengths. Here, we show for the Fraunhofer A-band of molecular O{sub 2} that the standard telluric correction approach fails in this at lower spectral resolutions, as an example for the general problem. Doppler-shift-dependent errors of the restored flux may arise, which can amount to more than 50% in extreme cases, depending on the line shapes of the target spectral features. Two applications are discussed: the recovery of the O{sub 2} band in the reflected light of an Earth analog atmosphere, as facilitated potentially in the future using an orbiting starshade and a ground-based extremely large telescope; and the recovery of the intrinsic ratio of the K i lines in the post-nova V4332 Sgr tracing the optical depth of the emitting region, to exemplify the relevance using present-day instrumentation. We show how one should derive correction functions for the compensation of the error in dependence of radial velocity shift, spectral resolution, and target line-profile function by use of high-resolution atmospheric transmission modeling, which has to be solved for the individual case.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Meteorological characteristics of extreme ozone pollution events in China and their future predictions

Ozone (O 3 ) has become one of the most concerning air pollutants in China in recent decades. In this study, based on surface observations, reanalysis data, global atmospheric chemistry model simulations, and multi-model future predictions, meteorological characteristics conducive to extreme O 3 pollution in various regions of China are investigated, and their historical changes and future trends are analyzed. During the most severe O 3 polluted months, the chemical production of O 3 is enhanced under the hot and dry conditions over the North China Plain (NCP) in June 2018 and the Yangtze River Delta (YRD) in July 2017, while regional transport is the main reason for the severe O3 pollution over the Sichuan Basin (SCB) in July 2015 and the Pearl River Delta (PRD) in September 2019. Over the last 4 decades, the frequencies of high-temperature and low-relative-humidity conditions increased in 2000–2019 relative to 1980–1999, indicating that O 3 pollution in both the NCP and YRD has become more frequent under historical climate change. In the SCB and PRD, the occurrence of atmospheric circulation patterns similar to those during the most polluted months increased, together with the more frequent hot and dry conditions, contributing to the increases in severe O 3 pollution in the SCB and PRD during 1980–2019. In the future (by 2100), the frequencies of months with anomalous high temperature show stronger increasing trends in the high-forcing scenario (Shared Socioeconomic Pathway (SSP5-8.5)) compared to the sustainable scenario (SSP1-2.6) in China. It suggests that high anthropogenic forcing will not only lead to slow economic growth and climate warming but also likely result in environmental pollution issues.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Coupled Air Quality and Boundary-Layer Meteorology in Western U.S. Basins during Winter: Design and Rationale for a Comprehensive Study

The wintertime episodes of high aerosol concentrations occur frequently in urban and agricultural basins and valleys worldwide. These episodes often arise following development of persistent cold-air pools (PCAPs) that limit mixing and modify chemistry. While field campaigns targeting either basin meteorology or wintertime pollution chemistry have been conducted, coupling between interconnected chemical and meteorological processes remains an insufficiently studied research area. Gaps in understanding the coupled chemical–meteorological interactions that drive high-pollution events make identification of the most effective air-basin specific emission control strategies challenging. To address this, a September 2019 workshop occurred with the goal of planning a future research campaign to investigate air quality in western U.S. basins. Approximately 120 people participated, representing 50 institutions and five countries. Workshop participants outlined the rationale and design for a comprehensive wintertime study that would couple atmospheric chemistry and boundary layer and complex-terrain meteorology within western U.S. basins. Participants concluded the study should focus on two regions with contrasting aerosol chemistry: three populated valleys within Utah (Salt Lake, Utah, and Cache Valleys) and the San Joaquin Valley in California. This paper describes the scientific rationale for a campaign that will acquire chemical and meteorological datasets using airborne platforms with extensive range, coupled to surface-based measurements focusing on sampling within the near-surface boundary layer, and transport and mixing processes within this layer, with high vertical resolution at a number of representative sites. No prior wintertime basin-focused campaign has provided the breadth of observations necessary to characterize the meteorological–chemical linkages outlined here, nor to validate complex processes within coupled atmosphere–chemistry models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Meteoroid Stream of Comet C/1961 T1 (Seki) and Its Relation to the December ρ-Virginids and γ-Sagittariids

We answer the question of whether long-period comet C/1961 T1 (Seki) is the parent body of a meteor shower observed in Earth’s atmosphere. We model discrete parts of a meteoroid stream associated with this comet and follow the orbital evolution of its particles via a numerical integration of their orbits. The gravitational perturbations of major planets and the Poynting–Robertson effect, which influences the motion of the particles, are taken into account. We predict two showers which are related to C/1961 T1. One of these showers was identified with the December ρ-Virginids, #502. This identification confirms an earlier suggested relationship between the shower and the comet. The second predicted shower was identified with the γ-Sagittariids, #657. The radiant areas of both showers are symmetrical in respect to the apex of Earth’s heliocentric motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Constraints on Metastable Helium in the Atmospheres of WASP-69b and WASP-52b with Ultranarrowband Photometry

Infrared observations of metastable 2{sup 3}S helium absorption with ground- and space-based spectroscopy are rapidly maturing, as this species is a unique probe of exoplanet atmospheres. Specifically, the transit depth in the triplet feature (with vacuum wavelengths near 1083.3 nm) can be used to constrain the temperature and mass-loss rate of an exoplanet’s upper atmosphere. Here, we present a new photometric technique to measure metastable 2{sup 3}S helium absorption using an ultranarrowband filter (FWHM 0.635 nm) coupled to a beam-shaping diffuser installed in the Wide-field Infrared Camera on the 200 inch Hale Telescope at Palomar Observatory. We use telluric OH lines and a helium arc lamp to characterize refractive effects through the filter and to confirm our understanding of the filter transmission profile. We benchmark our new technique by observing a transit of WASP-69b and detect an excess absorption of 0.498% ± 0.045% (11.1σ), consistent with previous measurements after considering our bandpass. We then use this method to study the inflated gas giant WASP-52b and place a 95th percentile upper limit on excess absorption in our helium bandpass of 0.47%. Using an atmospheric escape model, we constrain the mass-loss rate for WASP-69b to be 5.25{sub −0.46}{sup +0.65}×10{sup −4} M{sub J} Gyr{sup −1} (3.32{sub −0.56}{sup +0.67}×10{sup −3} M{sub J} Gyr{sup −1}) at 7000 K (12,000 K). Additionally, we set an upper limit on the mass-loss rate of WASP-52b at these temperatures of 2.1×10{sup −4} M{sub J} Gyr{sup −1} (2.1×10{sup −3} M{sub J} Gyr{sup −1}). These results show that ultranarrowband photometry can reliably quantify absorption in the metastable helium feature.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Integrating Energy-Efficient Computing with Computational Research to Accelerate Energy Technology

NREL's computational sciences center hosts the largest high performance computing (HPC) capabilities dedicated to energy research while functioning as a living laboratory for energy-efficient computing. NREL's HPC capabilities support the research needs of the Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). In ten years of operation, HPC use in EERE-sponsored research has grown by a factor of 30, including work in electricity generation, energy efficiency, transportation, and energy system modeling. This paper analyzes this research portfolio, providing examples of individual use cases. The paper documents NREL's history of operating one of the world's most energy-efficient data centers while examining pathways to reduce economic and environmental impact beyond reduction of Power Usage Efficiency (PUE). This paper concludes by examining the unique opportunities created for accelerating improvements in data center efficiency created by combining an HPC system dedicated to energy research and a research program in energy-efficient computing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Benchmarking Cylindrical Blast Wave Theory Against the OSIRIS-REx Sample Return Capsule Reentry

Weak shock theory based on cylindrical blast waves has been used to interpret meteor infrasound, but it has not been systematically benchmarked against a non-ablating hypersonic source with independently known parameters. The objective of this study is not to propose a new theoretical framework, but to evaluate the operational validity of the existing suite of blast radius formulations against a high-fidelity ground truth dataset. The OSIRIS-REx Sample Return Capsule reentry on 24 September 2023 provides such a benchmark because the capsule geometry, trajectory, and infrasound emission points are constrained from mission data and ray tracing, reducing source-side uncertainty associated with ablation. Using observations from 39 infrasound stations, this benchmarking study evaluates six published blast radius (${R}_{0}$ ) formulations and three weak-shock transition coefficients ( C ) within a stratified atmospheric propagation model to predict signal period and peak overpressure. The benchmarking identifies the Sakurai formulation as the best-performing formulation for non-ablating bodies, with the Jones/Plooster formulation performing comparably when a physically appropriate C is adopted. Sakurai and Jones/Plooster yield linear-period median absolute percentage residuals of 9% and 11%, respectively. The period predictions show only weak sensitivity to C at these propagation distances. The Mach-diameter approximation commonly used in meteor studies overestimates ${R}_{0}$ by more than a factor of 3 in the absence of ablation. Finally, these results establish a performance baseline for applying cylindrical blast wave theory to effectively non-ablating hypersonic bodies and demonstrate that the signal period is a robust observable for constraining ${R}_{0}$.

atmospheric entry↗

DOE OSTI2023

Evaluating Shallow Convection Parameterization Assumptions With a q t – w Quadrant Analysis

Abstract Uncertainties associated with the response of shallow clouds to global warming remain a great challenge for climate projection. Due to the small size of these clouds, parameterizations are required to represent them in both current‐ and next‐generation climate models. We present a quantitative evaluation of several important assumptions used in both mass‐flux and Assumed‐PDF Higher‐Order Closure parameterizations of shallow convection. We use large‐eddy simulations (LESs) of four different shallow convection regimes as benchmarks and apply a q t (total moisture)– w (vertical velocity) quadrant analysis technique to identify the “coherent structures” of moist and dry up/down‐drafts. The statistics of these coherent structures and the environment are then used to evaluate commonly used assumptions in mass‐flux parameterizations related to (a) the convective downdrafts and (b) the subplume variability of convective updrafts. For quantitative evaluation of the PDF closure used in AP‐HOCs, we perform offline calculations with the PDF closure used in the Cloud Layers Unified By Bi‐normals and the Simplified Higher‐Order Closure schemes by supplying the closure with statistical quantities directly calculated from the LESs. Then, the impact of parameterization assumptions embedded in the PDF closure in the same two categories is analyzed and understood through comparisons of the output from the closure with the statistics of the moist and dry up/down‐drafts from the LESs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

TChem-atm (v2.0.0): scalable performance-portable multiphase atmospheric chemistry

We present TChem-atm, a performance-portable approach that enables efficient simulation of chemically detailed and multiphase atmospheric chemistry on modern heterogeneous computing architectures. Unlike previous efforts that rely on architecture-specific code or focus exclusively on gas-phase chemistry, TChem-atm supports fully coupled gas–aerosol systems with execution across CPUs, NVIDIA GPUs, and AMD GPUs through the Kokkos programming model. It integrates the flexible multiphase capabilities of the Community Atmospheric Model Chemistry Package (CAMP) with the high-performance kinetic routines of TChem, and includes automatic Jacobian construction with support for a range of stiff ODE solvers. In a proof-of-concept integration with the particle-resolved model PartMC, TChem-atm reproduces the existing PartMC–CAMP implementation within solver tolerances and delivers substantial GPU speedups, especially for large particle populations. Performance benchmarks reveal substantial speedups on GPU platforms, particularly for large particle populations, with consistent results across hardware backends. TChem-atm enables performance-portable execution across CPUs and GPUs, though optimal efficiency may require modest architecture-specific tuning (e.g., team and vector sizes), with up to a twofold improvement on the NVIDIA H100. It directly supports sectional and particle-resolved host models, while modal aerosol schemes require minor adaptation to provide particle-scale quantities such as representative diameters. By enabling chemically detailed, multiphase simulations with performance portability and host-model flexibility, TChem-atm facilitates the incorporation of advanced chemistry into atmospheric models.

Díaz-Ibarra, Oscar Homero [Sandia National Laborat↗

DOE OSTI2024

Impacts of future nuclear power generation on the international monitoring system

Many countries are considering nuclear power as a means of reducing greenhouse gas emissions, and the IAEA (IAEA, 2022) has forecasted nuclear power growth rates up to 224% of the 2021 level by 2050. Nuclear power plants release trace quantities of radioxenon, an inert gas that is also monitored under international agreements as a signature of nuclear weapons tests. To better understand how nuclear energy growth (and resulting Xe emissions) could affect this global nonproliferation architecture, we modeled daily releases of radioxenon isotopes used for nuclear explosion detection in the International Monitoring System (IMS) that is part of the Comprehensive Nuclear Test-Ban Treaty: 131m Xe, 133 Xe, 133m Xe, and 135 Xe to examine the change in the number of radioxenon detections as compared to the 2021 detection levels. If a 40-station IMS network is used, the detections of 133 Xe in 2050 would range from 82% for the low-power scenario to 195% for the high-power scenario, compared to the detections in 2021. If an 80-station IMS network is used, the detections of 133 Xe in 2050 would range from 83% of the 2021 detection rate for the low-power scenario to 209% for the high-power scenario. Essentially no detections of 131m Xe and 133m Xe are expected. The high growth scenario could lead to a six-fold increase in 135 Xe detections, but the total number of detections is still small (on the order of 1 detection per day in the entire network).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗

DOE OSTI2025

Enhanced Simulation of Atmospheric Blocking in a High‐Resolution Earth System Model: Projected Changes and Implications for Extreme Weather Events

Atmospheric blocking is closely linked to the occurrence of extreme weather events. However, low-resolution Earth system models often underestimate the frequency of blocking, undermining confidence in future projections. Here, in this study, we use the high-resolution Community Earth System Model (CESM-HR; 25 km atm and 10 km ocean) to show that CESM-HR reduces biases in atmospheric blocking for both winter and summer, particularly for events lasting longer than 10 days. This improvement is partly due to reduced sea surface temperature biases at higher resolution. Additionally, applying a bias correction to the 500 hPa geopotential height further enhances blocking frequency simulations, highlighting the crucial role of the mean state. Under the Representative Concentration Pathway 8.5 scenario, CESM-HR projects a decrease in wintertime blocking over regions such as the Euro-Atlantic and Chukchi-Alaska, consistent with previous studies. In contrast, summer blocking is expected to become more frequent and persistent, driven by weakened zonal winds. The blocking center shifts from historical locations over Scandinavia and eastern Russia to central Eurasia, significantly increasing blocking over the Ural region. Summer blocking frequency over the Scandinavia-Ural region may eventually surpass historical winter blocking over the Euro-Atlantic. This increase in summer blocking could exacerbate summer heatwaves in a warming climate, making severe heatwaves, like those observed recently, more common in the future.

Atmospheric blocking↗

DOE OSTI2021

Benchmark estimate of the effect of anthropogenic emissions on the ocean surface

Abstract Investigations into the role of anthropogenic emissions in the occurrence of extreme weather often use a method that compares simulations of atmospheric climate models run under a factual scenario of historical boundary conditions observed during the period of the event against simulations run under a counterfactual scenario of what those boundary conditions might naturally have been over that same period in the absence of anthropogenic emissions. A particular requirement for this experiment design is an accurate estimation of ocean surface boundary conditions for use by the counterfactual natural simulations. Here we use output from the CMIP5 multi‐climate‐model archive to develop a robust estimate of sea surface temperatures and sea ice conditions for use in counterfactual natural simulations, intended as a benchmark estimate to facilitate comparison across climate models and across studies. This development includes tests to ensure that the final estimate is stable from year‐to‐year and stable against other perturbations to the methodology, as well as consideration of the strengths and weaknesses in comparison to other available attributable warming estimates. While this estimate is tailored specifically for the International CLIVAR C20C+ Detection and Attribution Project, it can be used by related projects as well.

Stone, Dáithí A.↗

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↗