Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean observation”

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 325 records · Page 18

Formation and fate of freshwater on an ice floe in the Central Arctic

The melt of snow and sea ice during the Arctic summer is a significant source of relatively fresh meltwater. The fate of this freshwater, whether in surface melt ponds or thin layers underneath the ice and in leads, impacts atmosphere–ice–ocean interactions and their subsequent coupled evolution. Here, we combine analyses of datasets from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (June–July 2020) for a process study on the formation and fate of sea ice freshwater on ice floes in the Central Arctic. Our freshwater budget analyses suggest that a relatively high fraction (58 %) is derived from surface melt. Additionally, the contribution from stored precipitation (snowmelt) outweighs by 5 times the input from in situ summer precipitation (rain). The magnitude and rate of local meltwater production are remarkably similar to those observed on the prior Surface Heat Budget of the Arctic Ocean (SHEBA) campaign, where the cumulative summer freshwater production totaled around 1 m during both. A relatively small fraction (10 %) of freshwater from melt remains in ponds, which is higher on more deformed second-year ice (SYI) compared to first-year ice (FYI) later in the summer. Most meltwater drains laterally and vertically, with vertical drainage enabling storage of freshwater internally in the ice by freshening brine channels. In the upper ocean, freshwater can accumulate in transient meltwater layers on the order of 0.1 to 1 m thick in leads and under the ice. The presence of such layers substantially impacts the coupled system by reducing bottom melt and allowing false bottom growth; reducing heat, nutrient, and gas exchange; and influencing ecosystem productivity. Regardless, the majority fraction of freshwater from melt is inferred to be ultimately incorporated into the upper ocean (75 %) or stored internally in the ice (14 %). Terms such as the annual sea ice freshwater production and meltwater storage in ponds could be used in future work as diagnostics for global climate and process models. For example, the range of values from the CESM2 climate model roughly encapsulate the observed total freshwater production, while storage in melt ponds is underestimated by about 50 %, suggesting pond drainage terms as a key process for investigation.

54 ENVIRONMENTAL SCIENCES↗

Amplified bottom water acidification rates on the Bering Sea shelf from 1970–2022

The Bering Sea shelf supports a highly productive marine ecosystem that is vulnerable to ocean acidification (OA) due to the cold, carbon-rich waters. Previous observational evidence suggests that bottom waters on the shelf are already seasonally undersaturated with respect to aragonite (i.e. Ω arag <1) and that OA will continue to increase the spatial extent, duration, and intensity of these conditions. Here, we use a regional ocean biogeochemical model to simulate changes in ocean carbon chemistry for the Bering Sea shelf from 1970–2022. Over this timeframe, model results suggest that surface Ω arag decreases by −0.043 per decade and surface pH by −0.014 per decade, comparable to observed global rates of OA. However, bottom water pH decreases at twice the rate of surface pH, while bottom [H + ] decreases at nearly 3 times the rate of surface [H + ]. This amplified bottom water acidification has emerged over the past 25 years and is likely driven by a combination of anthropogenic carbon accumulation and increasing primary productivity and subsurface respiration and remineralization. Due to this enhanced bottom water acidification, the spatial extent of bottom waters with Ω arag <1 has greatly expanded over the past 2 decades, along with pH conditions harmful to red king crab. Interannual variability in surface and bottom Ω arag , pH, and [H + ] has also increased over the past 2 decades, resulting in part from the increased physical climate variability. We also find that the Bering Sea shelf is a net annual carbon sink of 1.1–7.9 Tg C yr −1 , with the range resulting from the difference in the two different atmospheric forcing reanalysis products used. Seasonally, the shelf is a significant carbon sink from April–October but a somewhat weaker carbon source from November–March.

54 ENVIRONMENTAL SCIENCES↗

Observing Arctic Sea Ice

Our understanding of Arctic sea ice and its wide-ranging influence is deeply rooted in observation. Advancing technologies have profoundly improved our ability to observe Arctic sea ice, document its processes and properties, and describe atmosphere-ice-ocean interactions with unprecedented detail. Yet, our progress toward better understanding the Arctic sea ice system is mired by the stark disparities between observations that tend to be siloed by method, scientific discipline, and application. This article presents a review and philosophical design for observing sea ice and accelerating our understanding of the Arctic sea ice system. We give a brief history of Arctic sea ice observations and showcase the 2018 melt season within the context of five observational themes: spatial heterogeneity, temporal variability, cross-disciplinary science, scalability, and retrieval uncertainty. We synthesize buoy data, optical imagery, satellite retrievals, and airborne measurements to demonstrate how disparate data sets can be woven together to transcend observational-scale issues. The results show that there are limitations to interpreting any single data set alone. However, many of these limitations can be surmounted by combining observations that cross spatial and temporal scales. We conclude the article with pathways toward coordinating across observational platforms in order to: (1) optimize the scientific, operational, and community return on observational investments, and (2) facilitate a richer understanding of Arctic sea ice and its role in the climate system.

54 ENVIRONMENTAL SCIENCES↗

Advancing Organized Convection Representation in the Unified Model: Implementing and Enhancing Multiscale Coherent Structure Parameterization

To address the effect of stratiform latent heating on meso- to large-scale circulations, an enhanced implementation of the Multiscale Coherent Structure Parameterization (MCSP) is developed for the Met Office Unified Model. MCSP represents the top-heavy stratiform latent heating from under-resolved organized convection in general circulation models. We couple the MCSP with a mass-flux convection scheme (CoMorph-A) to improve storm lifecycle continuity. The improved MCSP trigger is specifically designed for mixed-phase deep convective cloud, combined with a background vertical wind shear, both known to be crucial for stratiform development. We also test a cloud top temperature dependent convective-stratiform heating partitioning, in contrast to the earlier fixed partitioning. Assessments from ensemble weather forecasts and decadal simulations demonstrate that MCSP directly reduces cloud deepening and precipitation areas by moderating mesoscale circulations. Indirectly, it amends tropical precipitation biases, notably correcting dry and wet biases over India and the Indian Ocean, respectively. Remarkably, the scheme outperforms a climate model ensemble by improving seasonal precipitation cycle predictions in these regions. The scheme also improves Madden-Julian Oscillation (MJO) spectra, achieving better alignment with observational and reanalysis data by intensifying the simulated MJO over the Indian Ocean during phases 4 to 5. However, the scheme increases precipitation overestimation over the Western Pacific. Shifting from fixed to temperature-dependent convective-stratiform partitioning reduces the Pacific precipitation overestimation and further improves the seasonal cycle in India. Spatially correlated biases highlight the necessity for advances beyond deterministic approaches to align MCSP with environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

Stream and ocean hydrodynamics mediate partial migration strategies in an amphidromous Hawaiian goby

Abstract Partial migration strategies, in which some individuals migrate but others do not, are widely observed in populations of migratory animals. Such patterns could arise via variation in migratory behaviors made by individual animals, via genetic variation in migratory predisposition, or simply by variation in migration opportunities mediated by environmental conditions. Here we use spatiotemporal variation in partial migration across populations of an amphidromous Hawaiian goby to test whether stream or ocean conditions favor completing its life cycle entirely within freshwater streams rather than undergoing an oceanic larval migration. Across 35 watersheds, microchemical analysis of otoliths revealed that most adult Awaous stamineus were freshwater residents (62% of n = 316 in 2009, 83% of n = 274 in 2011), but we found considerable variation among watersheds. We then tested the hypothesis that the prevalence of freshwater residency increases with the stability of stream flows and decreases with the availability of dispersal pathways arising from ocean hydrodynamics. We found that streams with low variation of daily discharge were home to a higher incidence of freshwater residents in each survey year. The magnitude of the shift in freshwater residency between survey years was positively associated with predicted interannual variability in the success of larval settlement in streams on each island based on passive drift in ocean currents. We built on these findings by developing a theoretical model of goby life history to further evaluate whether mediation of migration outcomes by stream and ocean hydrodynamics could be sufficient to explain the range of partial migration frequency observed across populations. The model illustrates that the proportion of larvae entering the ocean and differential survival of freshwater‐resident versus ocean‐going larvae are plausible mechanisms for range‐wide shifts in migration strategies. Thus, we propose that hydrologic variation in both ocean and stream environments contributes to spatiotemporal variation in the prevalence of migration phenotypes in A. stamineus . Our empirical and theoretical results suggest that the capacity for partial migration could enhance the persistence of metapopulations of diadromous fish when confronted with variable ocean and stream conditions.

54 ENVIRONMENTAL SCIENCES↗

Impacts of Synoptic-Scale Dynamics on Clouds and Radiation in High Southern Latitudes

High-latitudinal mixed-phase clouds significantly affect Earth's radiative balance. Observations of cloud and radiative properties from two field campaigns in the Southern Ocean and Antarctica were compared with two global climate model simulations. A cyclone compositing method was used to quantify “dynamics-cloud-radiation” relationships relative to the extratropical cyclone centers. Observations show larger asymmetry in cloud and radiative properties between western and eastern sectors at McMurdo compared with Macquarie Island. Most observed quantities at McMurdo are higher in the western (i.e., post-frontal) than the eastern (frontal) sector, including cloud fraction, liquid water path (LWP), net surface shortwave and longwave radiation (SW and LW), except for ice water path (IWP) being higher in the eastern sector. Here, the two models were found to overestimate cloud fraction and LWP at Macquarie Island but underestimate them at McMurdo Station. IWP is consistently underestimated at both locations, both sectors, and in all seasons. Biases of cloud fraction, LWP, and IWP are negatively correlated with SW biases and positively correlated with LW biases. The persistent negative IWP biases may have become one of the leading causes of radiative biases over the high southern latitudes, after correcting the underestimation of supercooled liquid water in the older model versions. By examining multi-scale factors from cloud microphysics to synoptic dynamics, this work will help increase the fidelity of climate simulations in this remote region.

58 GEOSCIENCES↗

The GFDL Variable-Resolution Global Chemistry-Climate Model for Research at the Nexus of US Climate and Air Quality Extremes

We present a variable-resolution global chemistry-climate model (AM4VR) developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) for research at the nexus of US climate and air quality extremes. AM4VR has a horizontal resolution of 13 km over the US, allowing it to resolve urban-to-rural chemical regimes, mesoscale convective systems, and land-surface heterogeneity. With the resolution gradually reducing to 100 km over the Indian Ocean, we achieve multi-decadal simulations driven by observed sea surface temperatures at 50% of the computational cost for a 25-km uniform-resolution grid. In contrast with GFDL's AM4.1 contributing to the sixth Coupled Model Intercomparison Project at 100 km resolution, AM4VR features much improved US climate mean patterns and variability. In particular, AM4VR shows improved representation of: precipitation seasonal-to-diurnal cycles and extremes, notably reducing the central US dry-and-warm bias; western US snowpack and summer drought, with implications for wildfires; and the North American monsoon, affecting dust storms. AM4VR exhibits excellent representation of winter precipitation, summer drought, and air pollution meteorology in California with complex terrain, enabling skillful prediction of both extreme summer ozone pollution and winter haze events in the Central Valley. AM4VR also provides vast improvements in the process-level representations of biogenic volatile organic compound emissions, interactive dust emissions from land, and removal of air pollutants by terrestrial ecosystems. We highlight the value of increased model resolution in representing climate–air quality interactions through land-biosphere feedbacks. AM4VR offers a novel opportunity to study global dimensions to US air quality, especially the role of Earth system feedbacks in a changing climate.

54 ENVIRONMENTAL SCIENCES↗

Unique Temperature Trend Pattern Associated With Internally Driven Global Cooling and Arctic Warming During 1980–2022

Abstract Diagnosing the role of internal variability over recent decades is critically important for both model validation and projections of future warming. Recent research suggests that for 1980–2022 internal variability manifested as Global Cooling and Arctic Warming (i‐GCAW), leading to enhanced Arctic Amplification (AA), and suppressed global warming over this period. Here we show that such an i‐GCAW is rare in CMIP6 large ensembles, but simulations that do produce similar i‐GCAW exhibit a unique and robust internally driven global surface air temperature (SAT) trend pattern. This unique SAT trend pattern features enhanced warming in the Barents and Kara Sea and cooling in the Tropical Eastern Pacific and Southern Ocean. Given that these features are imprinted in the observed record over recent decades, this work suggests that internal variability makes a crucial contribution to the discrepancy between observations and model‐simulated forced SAT trend patterns.

54 ENVIRONMENTAL SCIENCES↗

Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.

54 ENVIRONMENTAL SCIENCES↗

A spatial evaluation of Arctic sea ice and regional limitations in CMIP6 historical simulations

The Arctic sea ice response to a warming climate is assessed in a subset of models participating in Phase 6 of the Coupled Model Intercomparison Project (CMIP6), using several metrics in comparison with satellite observations and results from the Pan-Arctic Ice Ocean Modeling and Assimilation System and the Regional Arctic System Model. Our study examines the historical representation of sea ice extent, volume, and thickness using spatial analysis metrics, such as the integrated ice-edge error, Brier score, and Spatial Probability Score. We find that the CMIP6 multi-model mean captures the mean annual cycle and 1979-2014 sea ice trends remarkably well. However, individual models experience a wide range of uncertainty in the spatial distribution of sea ice when compared against satellite measurements and reanalysis data. Our metrics expose common and individual regional model biases, which sea ice temporal analyses alone do not capture. We identify large ice edge and ice thickness errors in Arctic sub-regions, implying possible model specific limitations in or lack of representation of some key physical processes. We postulate that many of them could be related to the oceanic forcing, especially in the marginal and shelf seas, where seasonal sea ice changes are not adequately simulated. We therefore conclude that an individual model’s ability to represent the observed/reanalysis spatial distribution still remains a challenge. We propose the spatial analysis metrics as useful tools to diagnose model limitations, narrow down possible processes affecting them, and guide future model improvements critical to the representation and projections of Arctic climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The LAKE model input dataset for three Arctic lakes

This dataset contains meteorological data collected for three Arctic lakes and compiled to satisfy input requirements of the LAKE 2.0 model. The dataset was generated to act as a benchmarking dataset for future model-data inter-comparisons. The LAKE 2.0 model simulates temperatures within the water later and the sedimentary layer of a lake. The LAKE2.0. is an open-source code and available to download via this weblike http://tesla.parallel.ru/Viktor/LAKE/-/wikis/LAKE-model (last visit July 14, 2021). The meteorological data are required to simulate the surface energy balance at the surface of a lake. This dataset includes a compilation of the meteorological data pulled from multiple data streams, including National Oceanic and Atmospheric Administration (NOAA) climate data, Circumarctic Lakes Observation Network (CALON) data, and the United States Geological Survey (USGS) data. The data were compiled for three Arctic lakes: FoxDen (66.55877, -164.45670), Atqasuk (70.452497, -156.951984), and Toolik (68.63150, -149.60740). Each meteorological data is in comma-delimited format (file extension ‘.dat’) and includes eight columns: Temperature [K], Pressure [Pa], longwave downward radiation [W/m2], shortwave downward radiation [W/m2], “U” wind speed [m/s], ”V” wind speed [m/s], humidity [kg/kg], precipitation [m/s]. In addition to the meteorological data file, we included setup and driver files. The Toolik lake is the deepest out of three lakes and has inflowing and outflowing groundwater data. InflowOutflowREADME.txt has more information about inflow and outflow flies. The other two lakes are much shallower and modeled as a closed system (i.e. no water inflow or outflow).

54 ENVIRONMENTAL SCIENCES↗

The AtmoSOFAR Channel: First Direct Observations of an Elevated Acoustic Duct

The Sound Fixing and Ranging (SOFAR) channel in the ocean allows for low frequency sound to travel thousands of kilometers, making it particularly useful for detecting underwater nuclear explosions. Suggestions that an elevated SOFAR-like channel should exist in the stratosphere date back over half a century and imply that sources within this region can be reliably sensed at vast distances. However, this theory has not been supported with evidence of direct observations from sound within this channel. Here we show that an infrasound sensor on a solar hot air balloon recorded the first infrasound detection of a ground truth airborne source while within this acoustic channel, which we refer to as the AtmoSOFAR channel. Our results support the existence of the AtmoSOFAR channel, demonstrate that acoustic signals can be recorded within it, and provide insight into the characteristics of recorded signals. Results also show a lack of detections on ground-based stations, highlighting the advantages of using balloon-borne infrasound sensors to detect impulsive sources at altitude.

54 ENVIRONMENTAL SCIENCES↗

Advancing the Understanding of Cloud Microphysical Processes and Aerosol Indirect Effects in High-Latitude Mixed-Phase Clouds by Linking ARM Measurements with Climate Model Simulations (Final Report)

The key objectives of this project were to advance our understanding of cloud microphysical characteristics and aerosol indirect effects on mixed-phase clouds in high latitudes. To improve the representation of ice and mixed-phase clouds in Earth System Models (ESMs), we propose an integrated observation and modeling study of cloud macro- and microphysical properties, including spatial heterogeneities, mass partitioning between ice crystals and supercooled liquid water, effects of ice nucleating particles (INPs), and efficiency of secondary ice production (SIP), etc. Specifically, we took four main approaches in this project: (1) examining macro- and microphysical properties of ice and mixed-phase clouds based on in-situ and ground-based observations from multiple field campaigns funded by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program, including the Mixed-Phase Arctic Cloud Experiment (M-PACE), Indirect and Semi-Direct Aerosol Campaign (ISDAC), Ice Nucleating Particle Sources at Oliktok Point (INPOP), ARM West Antarctic Radiation Experiment (AWARE), Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS), and Macquarie Island Cloud and Radiation Experiment (MICRE); (2) evaluating the DOE Energy Exascale Earth System Model (E3SM) simulations based on observations, particularly for ice and mixed-phase cloud microphysical properties; (3) examining the impacts of INPs on ice and mixed-phase clouds. Specifically, a series of comparisons were conducted using observations over the Arctic, Southern Ocean, and Antarctica, including comparisons between the lower and higher southern latitudes as well as comparisons between the northern and southern hemispheres. In addition, aerosol indirect effects from distinct sources of dust particles were examined; and (4) investigating the impacts of SIP. Ultimately, these results helped to improve cloud microphysics and aerosol-cloud interaction parameterizations in the E3SM model. Overall, the project provided improved understanding regarding various factors, including thermodynamic, dynamic, and aerosol conditions, on the micro- and macrophysical properties of ice and mixed-phase clouds in the high latitudes. Resulting analysis helped to provide an improved physical basis for refining the current cloud microphysics parameterizations related to ice and mixed-phase clouds in E3SM.

54 ENVIRONMENTAL SCIENCES↗

Relative Contributions of Atmospheric, Oceanic, and Coupled Processes to North Pacific and North Atlantic Variability

Patterns of sea surface temperature (SST) variability over the northern oceans arise from a combination of atmospheric, oceanic, and coupled processes. In this work, we use a novel methodology and a suite of observations to quantify the processes contributing to the dominant patterns of interannual SST variability over these regions. We decompose the upper ocean heat content tendency associated with such dominant patterns into contributions from different heat fluxes: (a) atmospherically driven, (b) surface feedbacks, and (c) oceanic. We find that in the subtropics, cloud radiative flux, turbulent heat flux, and residual oceanic processes each contributes substantially to North Pacific SST variability, whereas turbulent heat flux primarily induces North Atlantic SST variability. Cloud radiative fluxes therefore provide a major source of interannual SST variability in the North Pacific but not in the North Atlantic. In midlatitudes, SST fluctuations over the northern oceans are driven by the combination of turbulent and oceanic heat fluxes.

54 ENVIRONMENTAL SCIENCES↗

A Global High‐Resolution Mesoscale Convective System Database Using Satellite‐Derived Cloud Tops, Surface Precipitation, and Tracking

Abstract A new methodology is developed to construct a global (60°S–60°N) long‐term (2000–2019) high‐resolution (∼10‐km h) mesoscale convective system (MCS) database by tracking MCS jointly using geostationary satellite infrared brightness temperature ( T b ) and precipitation feature (PF) characteristics from the Integrated Multi‐satellitE Retrievals for GPM precipitation data sets. Independent validation shows that the satellite‐based MCS data set is able to reproduce important MCS statistics derived from ground‐based radar network observations in the United States and China. We show that by carefully considering key PF characteristics in addition to T b signatures, the new method significantly improves upon previous T b ‐only methods in detecting MCSs in the midlatitudes for all seasons. Results show that MCSs account for over 50% of annual total rainfall across most of the tropical belt and in selected regions of the midlatitudes, with a strong seasonality over many regions of the globe. The tracking database allows Lagrangian aspects such as MCS lifetime and translational speed and direction to be analyzed. The longest‐lived MCSs preferentially occur over the subtropical oceans. The land MCSs have higher cloud‐tops associated with more intense convection, and oceanic MCSs have much higher rainfall production. While MCSs are observed in many regions of the globe, there are fundamental differences in their dynamic and thermodynamic structures that warrant a better understanding of processes that control their evolution. This global database provides significant opportunities for observational and modeling studies of MCSs, their characteristics, and roles in regional and global water and energy cycles, as well as their hydrologic and other impacts.

54 ENVIRONMENTAL SCIENCES↗

Predictive links between microbial communities and biological oxygen utilization in the Arctic Ocean

Microbial metabolism influences rates of net community production (NCP), exerting a direct biological control on marine oxygen and carbon fluxes. In the Arctic, it is increasingly important to understand and quantify this process, as ecological and oceanographic conditions shift due to changing climate. Here, we describe potential ecological links between pelagic microbial diversity and an NCP precursor, biological oxygen utilization, using machine learning and paired observations of community structure and metabolic activity from a seasonally and spatially variable transect of the Arctic Ocean (2019–2020 MOSAiC Expedition). Community structure was determined using 16S (prokaryotic) and 18S (eukaryotic) rRNA gene amplicon sequencing, and metabolic activity was derived from ΔO 2 /Ar. Using self-organizing maps, we identified clear successional patterns in observed microbial community structure that were seasonally driven in the upper ocean and vertically stratified with depth. Metabolic activity was also stratified, with a primarily net heterotrophic water column (median −1.5% biological oxygen saturation), excepting periodic oxygen supersaturation (maximum: 13.6%) within the mixed layer. Using DNA sequences as predictor variables, we then constructed a random forest regression model that reliably reconstructed biological oxygen concentrations (root mean squared error = 4.14 μmol kg −1 ). Top predictors from this model were from heterotrophic (bacteria) or potentially mixotrophic (dinoflagellate) taxa. These analyses highlight biologically driven diagnostic tools that can be used to expand biogeochemical datasets and improve the microbial perspectives and metabolisms represented in ecological models of net productivity and carbon flux in a changing Arctic Ocean.

Chamberlain, Emelia J. [Univ. of San Diego, San Di↗

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗