SEARCH · Search NASA

Results for “surface heat flux”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Surface heat flux and its association with the MJO in the tropical western Pacific using ARM observations

The Madden-Julian oscillation (MJO) is a major atmospheric phenomenon in the tropics that moves eastward every 20 to 100 days. It brings heavy rain and strong winds, influencing extreme weather events far beyond the tropics – including flooding, hurricanes, tornadoes and heavy snow in the United States. Therefore, an improved understanding of the MJO is critical for enhancing weather forecasts and supporting better decision-making for communities, emergency managers, and the private sector. Past studies, based on short-term observations or long-term model data, have emphasized the dominant role of atmospheric humidity in driving the MJO. Our research, using long-term observations (2000-2014) from three U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) sites in the tropical western Pacific (Manus, Nauru, and Darwin) confirms those past findings. However, we also show that atmospheric temperature, especially in the mid and upper troposphere, also plays a key role in the MJO’s evolution from its quiet (suppressed) to convective (active) phase. These findings provide valuable insights into how the MJO evolves and may help evaluate and improve weather and climate models. Our study also examined how the exchange of heat between the Earth's surface and the atmosphere, called surface heat flux, interacts with the MJO. This flux includes components such as latent heat (related to evaporation), sensible heat (related to temperature difference between surface and atmosphere), and radiation. We found that modulation in MJO convection is well connected to surface heat flux over the tropical western Pacific. One surface heat flux component often overlooked is the sensible heat flux due to precipitation (Q P ). Since falling raindrops are typically cooler than the land or ocean surface, they can cool the surface and affect local weather. To test the impact of Q P on convection during MJO, we incorporated it into a weather model, and ran simulations for two periods: April 2009 (when the MJO was active) and June 2006 (when it was not). Including Q P reduced errors in simulating the daily peak of rainfall—by 83% during the strong MJO and 23% during the inactive phase. It also improved the timing of peak daily rainfall and reduced the overall precipitation error by about 10%. These improvements were especially noticeable during periods of heavy rainfall. Our results suggest that even small heat exchanges from rainfall can play a significant role in shaping local weather. Accounting for these effects, particularly in regions like the tropical islands, can lead to more accurate simulations of the MJO-associated precipitation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗

DOE OSTI2025

Cloud radiative effect dominates variabilities of surface energy budget in the dark Arctic

Climate models simulate a wide range of temperatures in the Arctic. Here we investigate one of the main drivers of changes in surface temperature: the net surface heat flux in the models. We show that in the winter months of the dark Arctic, there is a more than two-fold difference in the net surface heat fluxes among the models, and this difference is dominated by the downward infrared radiation from clouds. Owing to the small amount of water vapor in the winter Arctic, infrared radiation from clouds transmits more easily to the surface in the Arctic than at other latitudes, resulting in large cloud radiative effect at the surface. The dominant role of the cloud effect is also found in the transient variability of the net surface heat flux. Results demonstrate that accurate simulation of clouds is crucial for determining the net surface heat flux, which in turn affects surface temperature and sea ice properties in the Arctic.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Relative Contributions of Surface Parcel and Atmospheric Environment to Convective Potential during Interstorm Periods

Convective instability is affected by both surface heat fluxes into near-surface air parcels and the atmospheric profile encountered as these air parcels are lifted aloft. In this study, we investigate the evolution of buoyancy experienced by surface-based parcels, taking interstorm soil moisture drydowns as fundamental time periods. We use profile observations of atmospheric temperature and humidity from ground-based passive remote sensing and in situ observations available at the U.S. Department of Energy Atmospheric Radiation Measurement Southern Great Plains site in Oklahoma, United States, and examine drydown events during the warm seasons from 2017 to 2022. Results reveal new insights on the emergence and growth of convective available potential energy (CAPE) during drydowns. We observe that the peak convective potential appears in the afternoon and grows in the midtroposphere (68% of 66 drydown events), with the largest CAPE values toward the end of the drydown event. We also find that the near-surface parcel contributes more to the evolution of CAPE than the atmospheric environment (86% of drydown events). These findings highlight the importance of surface heat fluxes and moisture availability from surface-based parcels in modulating the evolution of atmospheric conditions favorable for convection, underscoring the interplay between surface processes and atmospheric dynamics in driving convective instability.

Atmosphere-land interaction↗

DOE OSTI2025

The Atlantic Meridional Overturning Circulation’s Response to CO2 Increase: Assessing the Roles of Surface Flux and Oceanic Advection Feedbacks

Abstract The Atlantic meridional overturning circulation (AMOC) is projected to slow down in climate models due to greenhouse gas emissions. However, the physical mechanisms determining the rate of the projected AMOC slowdown remain unclear. Accordingly, this study isolates the roles of oceanic advection and surface flux feedbacks that might accelerate or decelerate the AMOC’s weakening using carbon dioxide (CO 2 ) quadrupling simulations in the CESM1.2 model. Surface flux feedbacks are isolated in partially coupled experiments in which either all surface flux components or the momentum flux responses to AMOC’s weakening that might provide feedback are suppressed, while a tracer decomposition of ocean density anomalies isolates the advection feedbacks. Comparing the ocean density components in the experiments shows that the AMOC’s response is initially determined by CO 2 -induced anomalous surface heat fluxes, and afterward, feedbacks determine its response. In the fully coupled case, surface heat flux feedback strongly promotes AMOC slowdown and causes its near shutdown, while a weaker but active AMOC is maintained when the surface flux feedback is inhibited in the partially coupled case. The positive surface heat flux feedback works by canceling out the negative oceanic heat advection feedback on deep-water formation in the subpolar North Atlantic (SPNA). With the heat advection feedback thus reduced, the positive salinity advection feedback becomes the dominant contributor to SPNA density changes and deep-water formation. In the partially coupled case, negative ocean heat advection feedback and the CO 2 -induced subtropical Atlantic saline anomalies imported into the SPNA play a stabilizing role. The results highlight the importance of SPNA salinity gradients and gyre circulation strength in determining the AMOC’s response rate or recovery.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗

DOE OSTI2026

Concept of Dynamic Heat Insulation for Rotating Detonation Engines

This work introduces a new class of materials concept to dynamically reduce instantaneous heat fluxes in Rotating Detonation Engine (RDE) combustor chamber walls. The high-frequency and high-amplitude surface heat fluxes observed in RDEs arise from large instantaneous temperature differences between the detonation shockwave and chamber wall surface. These temperature gradients drive substantial energy loss and reduce the chamber gas pressure, ultimately limiting the cycle’s thermodynamic efficiency. This work introduces a concept for dynamically insulating the combustion chamber surfaces using surface layers or coatings with low thermal time scale. With such coatings, the surface temperature may follow the fluctuations of the cyclic detonation wave temperature, thus reducing the instantaneous heat flux therefore cycle-mean heat flux. To analyze these cyclic thermal phenomena, a one-dimensional analytical conduction solver was utilized with the capability to handle multilayered structures. Parametric modeling was performed using transient heat flux boundary conditions representative of a hydrogen–air RDE across a broad range of coating thermal properties and engine conditions. The coating effectiveness scaled with the product of thermal time constant and detonation wave frequency and the results were non-dimensionalized to guide future materials development. This strategy may offer benefits in increasing material survivability, reducing cooling requirements, and enhancing pressure gain.

heat transfer↗

DOE OSTI2026

Large Eddy Simulations of the Atmospheric Boundary Layer Over Satellite‐Sensed Sea Ice Maps

Surface heterogeneity in the marginal ice zone (MIZ) causes multiscale secondary atmospheric circulations that are challenging to model or observe. The absence or inadequate representation of these circulations in ocean‐atmosphere exchange schemes in climate models is partially responsible for the underestimation of Arctic sea ice loss. Observationally, such circulations obfuscate the interpretation of polar atmospheric chemistry measurements, among others. To address this open challenge, large‐eddy simulations are conducted over real‐world satellite‐sensed sea ice maps with an ice‐sea temperature contrast, as well as over idealizations of these maps that alter the ice pattern but conserve its fraction, showing that the ice fraction of a surface is not sufficient to predict the surface heat flux. In a second suite of simulations, three other heterogeneity metrics (representing the surface fragmentation, patch edge tortuosity, and patch size variability) are introduced to complement the ice fraction in describing the surface. Simulations varying these parameters suggest that they also significantly impact surface‐air interactions. A roughness contrast is then added to the surface temperature contrast, indicating that the contribution of roughness changes to the resulting atmospheric circulations is less pronounced than that of thermal heterogeneity. Based on these results, we illustrate, using a multi‐linear regression on these map features, that MIZ surface parameterizations in Earth Systems Models can be improved if they account for these various characteristics of the sea‐ice patterns.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Self-consistent modeling of tokamak edge plasma transport with lithium sources

Magnetic confinement fusion devices require effective heat and particle exhaust solutions on the divertor plates to operate sustainably, especially under reactor-relevant conditions. Liquid lithium divertors have been proposed to address two major challenges: control of excessive heat flux to plasma-facing components through vapor shielding and minimization of core plasma contamination from impurities. The National Spherical Torus Experiment-Upgrade (NSTX-U) will explore lithium as a divertor material due to its potential to meet both objectives. We present a self-consistent coupling framework between the plasma boundary transport code UEDGE and the lithium wall transport code Wall–Li to evaluate the feasibility and operational limits of lithium-based divertors. The model aims to optimize lithium sourcing levels to prevent core plasma contamination via fuel dilution while ensuring divertor protection through vapor shielding. This integrated framework, applicable to any tokamak with lithium sources, dynamically adjusts lithium sourcing based on local plasma conditions and surface temperature. The coupled model is tested using NSTX-like geometry and plasma conditions to assess its performance and reliability. Wall–Li calculates lithium fluxes from plasma-facing components, incorporating physical sputtering, thermally enhanced sputtering, and evaporation driven by surface temperature and ion flux. These fluxes are reintroduced into UEDGE as neutral lithium atoms, enabling simulation of their transport and distribution within the plasma. UEDGE computes plasma and neutral transport, surface heat flux, and iteratively feeds this information back to Wall–Li. A small time step is employed to ensure numerical stability and convergence, enabling accurate simulations over typical tokamak discharge durations. This integrated modeling approach provides a robust tool for identifying operational regimes that balance effective lithium sourcing with minimal core plasma contamination, offering critical insights for optimizing lithium-based divertor systems in current and future fusion devices.

Magnetic confinement fusion↗

DOE OSTI2025

CFD Simulations of Molten Salts Within Plasma-Facing Components with Twisted Tape Inserts Under One-Sided–Heated Conditions

For liquid immersion blankets (LIBs), the heat removal within different plasma-facing components (PFCs) requires substantial heat removal capabilities that could benefit from the usage of heat transfer enhancements. Of consideration are twisted tape insert heat transfer enhancements that are ribbon-like metal pieces inserted within cooling channels that increase fluid mixing, resulting in higher heat transfer and pressure drop. Unlike traditional applications, PFCs such as divertor coolant tubes experience a significantly higher heat flux on one side versus the other or entirely one-sided heating. These conditions are not well studied for molten salts where the higher Prandtl numbers characteristic of the fluids are significantly higher than water, gases, and liquid metals. As an effort to study this knowledge gap, large eddy simulations were done with Prandtl numbers representative of FLiBe at proposed operating conditions. In this study, the effects of applying a surface heat flux to one side of a tube versus the entire tube with a twisted tape insert were investigated for the heat transfer coefficient and effects on the temperature profile. The findings of the study suggest that for Prandtl number fluids above unity, the Nusselt number was largely unaffected by the change in boundary condition. Both boundary conditions showed good agreement with the Manglik and Bergles Nusselt number correlation for tight-fitting twisted tapes. Finally, the outcomes of this study will be useful for future design studies of LIBs to provide fundamental guidance implementing twisted tape inserts and their performance in one-sided–heated conditions.

molten salts↗

DOE OSTI2026

Urban Land Surface Effects on Summertime Clouds and Moist Convection in Houston Under Different Synoptic Conditions

Urban landscapes modify cloud formation and convection through complex thermodynamic and aerodynamic processes; however, their influence under different synoptic regimes remains poorly understood. This study investigates the impact of the Houston metropolitan area on summertime cloud cover and convective cell characteristics using a combination of satellite observations, radar data, and high-resolution process-based modeling. We isolate urban effects by comparing model simulations with realistic urban land cover against hypothetical scenarios where all urban areas are replaced by rural vegetation. Results reveal that Houston's land cover consistently enhances cloud fraction and convective activity relative to surrounding rural areas, altered by large-scale meteorological forcing. Under weakly forced conditions, enhanced surface heat flux primarily contributes to driving low-level convergence and vertical ascent, leading to over 8% cloud fraction increase between 2 and 6 km over the city. Under strongly forced conditions, urban influences manifest differently depending on cloud type. For non-convective clouds, the city acts as a barrier that decelerates and lifts moist southerly inflow, increasing low cloud cover over the urban core, while decreasing it downwind the city. For convective clouds, both synoptic ascent and urbanization modulate moisture redistribution and cloud structure, producing modest cloud enhancement over the city and slight suppression over the downwind area. Urbanization exerts small changes in the intensity of the convective cells; however, it significantly decreases their duration and traveling distance. This work highlights the importance of accounting for land surface heterogeneity in modeling clouds and precipitation and demonstrates that urban impacts on clouds are highly regime-dependent.

Liu, Ye [Pacific Northwest National Laboratory (PN↗

DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Unsteady Land-Sea Breeze Circulations in the Presence of a Synoptic Pressure Forcing

Unsteady land-sea breezes (LSBs) that result from time-varying surface temperature contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, M g , oriented from sea to land (α = 0°) or land to sea (α = 180°). Large eddy simulations reveal the development of four distinctive regimes, depending on the joint interaction between M g , α, and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the velocity between the unsteady and the corresponding steady scenarios at Δθ = Δθ max , are found to be significant and to extend 2 hr longer for α = 0° compared to α = 180°. These diurnal dynamics result in nonequilibrium conditions that are significantly affected by the flow history, and that behave differently over the two patches for the different α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land to sea (α = 180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the turbulence-mean flow equilibrium response. Finally, the vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution and can thus be used to identify these regimes using k-means clustering.

58 GEOSCIENCES↗

DOE OSTI2025

Reduction of Pacific Double‐ITCZ Bias by Convection Parameterization in NCAR CESM2.2

Abstract The impact of convective closure on the double‐ITCZ bias in the NCAR CESM2.2 is investigated in this study. The standard CESM2.2 simulates a remarkable double‐ITCZ bias in the central and eastern Pacific, especially in boreal winter and spring. Modifications to the closure in convection parameterization scheme greatly reduce the double‐ITCZ bias in all seasons, demonstrating that convection parameterization can substantially influence the double‐ITCZ bias in CESM2.2. Further analyses suggest that convection parameterization can modulate the tropical atmosphere‐ocean feedback processes, through which it influences the SST in the southern ITCZ region and hence the double‐ITCZ bias. The changes in the upper ocean temperature advection induced by modified convective closure plays important roles in reducing the warm SST bias and double‐ITCZ precipitation bias in the southern ITCZ region. The modified convective closure improves the low‐level cloud and shortwave cloud radiative forcing in the southeastern Pacific. However, surface heat flux plays only a limited role in reducing warm SST bias and double ITCZ bias because the impacts of shortwave radiation changes are largely canceled by changes in longwave radiation and latent heat flux.

Song, Xiaoliang [Scripps Institution of Oceanograp↗

DOE OSTI2026

Thermodynamics-guided machine learning model for predicting convective boundary layer height and its multi-site applicability

Accurate estimation of convective boundary layer height (CBLH) is vital for weather, climate, and air quality modeling. Machine learning (ML) shows promise in CBLH prediction, but input parameter selection often lacks physical grounding, limiting generalizability. This study introduces a novel ML framework for CBLH prediction, integrating thermodynamic constraints and the diurnal CBLH cycle as an implicit physical guide. Boundary layer growth is modeled as driven by surface heat fluxes and atmospheric heat absorption represented with the low tropospheric stability, using the diurnal cycle as input and output. TPOT and AutoKeras are employed to select optimal models, validated against Doppler lidar-derived CBLH data, achieving an R 2 of 0.84 across untrained years. Comparisons of eddy covariance (ECOR) and energy balance Bowen ratio (EBBR) flux measurements show the same prediction capability. Models trained on the ARM SGP C1 site with ECOR data and tested at E37 and E39 yield R 2 values of 0.79 and 0.81, respectively, demonstrating their adaptability. The ML model trained with all sites' data slightly enhances the performance compared with ML models trained over single-site data. The interquartile range for predicted CBLH is consistently narrower than that for DL-derived CBLH, reflecting lower variability in predicted CBLH compared to DL-derived CBLH, which is influenced by additional factors, which are not well represented with the model inputs. The model's generalizability across multiple sites at the ARM SGP site demonstrates its potential for transfer to greater distances, offering a scalable approach for enhancing boundary layer parameterization in atmospheric models.

Chu, Yufei [Stony Brook Univ., NY (United States)]↗

DOE OSTI2025

A 1 km soil moisture dataset over eastern CONUS generated by assimilating SMAP data into the Noah-MP land surface model

An improved fine-scale soil moisture (SM) dataset at 1 km grid spacing, covering much of the eastern continental US, was generated by assimilating 9 km Soil Moisture Active Passive (SMAP) SM data into the v4.0.1 Noah-MP land surface model. With 12 ensemble members, the assimilation was carried out using the ensemble Kalman filter algorithm within NASA's Land Information System. The SM analysis for 2016 was fully validated against in situ observations from four different networks and compared with four other existing datasets. Results indicate that this SM analysis surpasses other datasets in top-layer SM distribution, including a machine-learning-based product, despite all SM estimates being less heterogeneous than observed. The analysis of anomalous errors suggests that large similarity in intrinsic errors is likely due to overlapping data sources among the selected SM datasets. More detailed evaluations were performed over two geographic areas. The observations collected by the Atmospheric Radiation Measurement facility in Oklahoma suggest that soil temperature and surface heat fluxes are concurrently simulated with good accuracy. Investigation into the 2016 southeastern US drought response further indicates drier conditions and higher evapotranspiration estimates compared to GLEAMv4.1. Notably, large errors are associated with grids having clay soil textures, underscoring the need for refined model treatments for specific soil types to further improve SM estimates. The dataset is publicly available on Zenodo at https://doi.org/10.5281/zenodo.14370563 (Tai et al., 2024).

Tai, Sheng-Lun [Pacific Northwest National Laborat↗

DOE OSTI2024

Atmospheric Variability Drives Anomalies in the Bering Sea Air–Sea Heat Exchange

Abstract High latitudes, including the Bering Sea, are experiencing unprecedented rates of change. Long-term Bering Sea warming trends have been identified, and marine heatwaves (MHWs), event-scale elevated sea surface temperature (SST) extremes, have also increased in frequency and longevity in recent years. Recent work has shown that variability in air–sea coupling plays a dominant role in driving Bering Sea upper-ocean thermal variability and that surface forcing has driven an increase in the occurrence of positive ocean temperature anomalies since 2010. In this work, we characterize the drivers of the anomalous surface air–sea heat fluxes in the Bering Sea over the period 2010–22 using ERA5 fields. We show that the surface turbulent heat flux dominates the net surface heat flux variability from September to April and is primarily a result of near-surface air temperature and specific humidity anomalies. The airmass anomalies that account for the majority of the turbulent heat flux variability are a function of wind direction, with southerly (northerly) wind advecting anomalously warm (cool), moist (dry) air over the Bering Sea, resulting in positive (negative) surface turbulent flux anomalies. During the remaining months of the year, anomalies in the surface radiative fluxes account for the majority of the net surface heat flux variability and are a result of anomalous cloud coverage, anomalous lower-tropospheric virtual temperature, and sea ice coverage variability. Our results indicate that atmospheric variability drives much of the Bering Sea upper-ocean temperature variability through the mediation of the surface heat fluxes during the analysis period. Significance Statement A long-term ocean warming trend and a recent increase in marine heatwaves in the Bering Sea have been identified. Previous work showed that anomalies in the exchange of heat between the ocean and the atmosphere were the primary driver of Bering Sea temperature variability, but the processes responsible for the heat exchange anomalies were unknown. In this work, we show that the atmosphere is the primary driver of anomalies in the Bering Sea air–sea heat exchange and therefore plays an important role in altering the thermal state of the Bering Sea. Our results highlight the importance of understanding more about how the ocean and the atmosphere interact at high latitudes and how this relationship will be affected by future climate change.

Hayden, Emily E. [a College of Earth, Ocean, and A↗