Search NASA⌕ Search

SEARCH · Search NASA

Results for “ATMOSPHERIC TURBULENCE”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

CROCUS High-Frequency Measurements of CO₂, H₂O, Wind, and Temperature at University of Illinois Chicago

As of April 9, 2026: this dataset is currently being versioned to include data up to March 31, 2026. Once the versioning process is complete, new data files will be available for access. The dataset metadata will also be updated to reflect the data availability of the new data being versioned. Raw atmospheric measurements collected at the University of Illinois Chicago (UIC) as part of the Community Research on Climate and Urban Science (CROCUS) project. The data includes high-frequency measurements of carbon dioxide (CO₂) concentration, water vapor (H₂O) concentration, wind speed (U, V, W components), temperature, and atmospheric pressure. These raw data are recorded by the LI-7500DS Open Path CO₂/H₂O Analyzer and a sonic anemometer at a 10 Hz acquisition frequency, providing the necessary inputs for calculating fluxes of CO₂, H₂O, heat, and momentum. In addition to gas concentration measurements, the dataset includes diagnostic information from the instruments, including absorptance, sample and reference signals, and diagnostic values. The data were collected to study urban atmospheric conditions and contribute to flux calculations for urban climate research. These measurements form the basis for calculating 30-minute average fluxes of key atmospheric variables using the eddy covariance method. This dataset provides critical raw input data for researchers interested in atmospheric fluxes, urban air quality, and the interaction between the urban environment and atmospheric processes. The data is part of the U.S. Department of Energy’s Biological and Environmental Research (BER) program, under the CROCUS Urban Integrated Field Laboratory (UIFL) project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving the representation of shallow cumulus convection with the simplified-higher-order-closure–mass-flux (SHOC+MF v1.0) approach

Abstract. Parameterized boundary layer turbulence and moist convection remain some of the largest sources of uncertainty in general circulation models. High-resolution climate modeling aims to reduce that uncertainty by explicitly attempting to resolve deep moist convective motions. An example of such a model is the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) with a target global resolution of 3.25 km, allowing for a more accurate representation of complex mesoscale deep convective dynamics. Yet, small-scale planetary boundary layer turbulence and shallow convection still need to be parameterized, which in SCREAM is accomplished through the turbulent-kinetic-energy-based (TKE-based) simplified higher-order closure (SHOC) – a simplified version of the assumed-double-Gaussian-PDF (probability density function) higher-order-closure method. In this paper, we implement a stochastic-multiplume mass-flux (MF) parameterization of dry and shallow convection in SCREAM to go beyond the limitations of double-Gaussian-PDF closures and couple it to SHOC (SHOC+MF). The new parameterization implemented in a single-column model type version of SCREAM produces results for two shallow cumulus convection cases (marine and continental shallow convection) that agree well with the reference data from large-eddy simulations, thus improving the general representation of the thermodynamic quantities and their turbulent fluxes as well as cloud macrophysics in the model. Furthermore, SHOC+MF parameterization shows weak sensitivity to the vertical grid resolution and model time step.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗
DOE OSTI2020

Drizzle, Turbulence, and Density Currents Below Post Cold Frontal Open Cellular Marine Stratocumulus Clouds

Nine cases of postcold frontal marine stratocumulus clouds exhibiting open cellular mesoscale organization are analyzed to characterize the drizzle, turbulence, and density currents below them. Data collected by the vertically pointing Doppler cloud radar and multiple lidars part of the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site are used in these analyses. A total of 251 drizzle shafts passed over the site with 76 shafts sampled for more than 15 min by the vertically pointing instruments. On average the drizzle shafts were ~18 km wide with below-cloud drizzle water path of ~40 g m −2 , cloud base rain rate of 7.78 mm day −1 , and cloud base drizzle modal diameter of 393 μm. The widths of the drizzle shafts did not exhibit any relationship with the cloud base rain rate, below-cloud drizzle water path, and the total water removed from the cloud. The downdrafts in the lowest 500 m within the drizzle cells strengthened with increasing cloud base rain rates, while the updrafts within drizzle cells did not exhibit this behavior. On average, the surface air density, pressure, water vapor mixing ratio, and wind speed along the background wind during the drizzle shaft were higher by 1.88 g m −3 , 8.60 Pa, 0.24 g kg −1 , and 1.22 m s −1 , respectively, from their background values. On the other hand, the surface air temperature was lower by 0.43 K during the drizzle shaft than its background value. The coincident measurements highlight the drizzle-turbulence-surface coupling in these cloud systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Carbon Outgassing in the Antarctic Circumpolar Current Is Supported by Ekman Transport From the Sea Ice Zone in an Observation‐Based Seasonal Mixed‐Layer Budget

Despite its importance for the global cycling of carbon, there are still large gaps in our understanding of the processes driving annual and seasonal carbon fluxes in the high-latitude Southern Ocean. This is due in part to a historical paucity of observations in this remote, turbulent, and seasonally ice-covered region. Here, we use autonomous biogeochemical float data spanning 6 full seasonal cycles and with circumpolar coverage of the Southern Ocean, complemented by atmospheric reanalysis, to construct a monthly climatology of the mixed layer budget of dissolved inorganic carbon (DIC). We investigate the processes that determine the annual mean and seasonal cycle of DIC fluxes in two different zones of the Southern Ocean—the Sea Ice Zone (SIZ) and Antarctic Southern Zone (ASZ). We find that, annually, mixing with carbon-rich waters at the base of the mixed layer supplies DIC which is, in the ASZ, either used for net biological production or outgassed to the atmosphere. In contrast, in the SIZ, where carbon outgassing and the biological pump are weaker, the surplus of DIC is instead advected northward to the ASZ. In other words, carbon outgassing in the southern Antarctic Circumpolar Current (ACC), which has been attributed to remineralized carbon from deep water upwelled in the ACC, is also due to the wind-driven transport of DIC from the SIZ. These results stem from the first observation-based carbon budget of the circumpolar Southern Ocean and thus provide a useful benchmark to evaluate climate models, which have significant biases in this region.

Sauvé, Jade↗
DOE OSTI2025

Using ground-based lidar data to investigate the water–vapor budget in the daytime atmospheric boundary layer

The moisture advection term in the water–vapor budget equation is investigated with a combination of a vertically-staring water–vapor lidar and Doppler lidar systems. These instruments make it possible to get the mean profile of moisture tendency and the latent heat flux (LHF) divergence. We use data of the Land–Atmosphere Feedback Experiment (LAFE) at the Atmospheric Radiation Measurement (ARM) Program’s Southern Great Plains (SGP) site, Oklahoma, USA, collected on 30 August 2017 between 15 and 24 UTC, which corresponds to 09 to 18 LT. The lidars provide turbulence resolving profiles of moisture and vertical wind fluctuations. The LHF profile is derived from the covariance of these moisture and vertical wind fluctuations. The mean boundary layer height z i is determined from the peak of the moisture variance. The results demonstrate that the combination of two remote sensing instruments can be applied for determining the dominant water–vapor budget terms, namely moisture tendency, latent heat flux divergence and moisture advection.

Advection↗
DOE OSTI2026

Eddy Correlation Flux Measurement System (ECOR) Instrument Handbook

The eddy correlation (ECOR) flux measurement system provides in situ, half-hour measurements of surface turbulent fluxes of momentum, sensible heat, latent heat, and carbon dioxide (CO 2 ) at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility sites. ARM’s ECORs also historically measured methane fluxes at Utqiaġvik and Oliktok Point, Alaska. The fluxes are obtained with the eddy covariance technique, which involves correlation of the vertical wind component with the horizontal wind component, the air temperature, the water vapor density, and the CO 2 concentration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Barriers and variable spacing enhance convective cooling and increase power output in solar PV plants

When the temperature of solar photovoltaic (PV) modules rises, efficiency drops and module degradation accelerates. Thus, it is beneficial to reduce module operating temperatures. Previous studies of solar power plants have illustrated that incoming flow characteristics, turbulent mixing, and array geometry can strongly impact convective cooling, as measured by the convective heat transfer coefficient h. In the fields of heat transfer and plant canopy flow, previous work has shown that system-scale arrangement modifications—e.g., variable spacing, barriers, or windbreaks—can passively alter the flow, enhance turbulent mixing, and influence convection. However, researchers have not yet explored how variable spacing or barriers might enhance convective cooling in solar power plants. Here, high-resolution large-eddy simulations model the air flow and heat transfer through solar power plant arrangements modified with missing modules and barrier walls. We then perform a control volume analysis to evaluate the net heat flux and compute h, which quantifies the influence of these spatial modifications on convective cooling and, thus, module temperature and power output. Installing barrier walls yields the greatest improvements, increasing h by 3.4%, reducing module temperature by an estimated 2.5 °C, and boosting power output by an estimated 1.4% on average. These findings indicate that incorporating variable spacing or barrier-type elements into PV plant designs can reduce module temperature and, thus, improve PV performance and service life.

14 SOLAR ENERGY↗
DOE OSTI2023

A Comparison of Adverse Meteorology Inputs for Emergency Management and Safety Basis Dispersion Modeling Applications at SRS

The MELCOR Accident Consequence Calculation System (MACCS), used for safety basis accident analyses at SRS, implements adjustment factors to consider the increased vertical turbulence associated with local surface roughness and the increased lateral turbulence associated with time-based plume meander. The two models used for SRS emergency planning radiological consequence assessment, i.e., Hotspot and Puff/Plume, currently do not have an option to input or apply diffusion adjustments for meander or surface roughness. Adverse meteorology conditions at SRS have been determined as E stability class and 1.3 meters per second wind speed for ground level releases. This is the input used for emergency planning models. Contrastingly safety basis analyses use a different input, representing adverse meteorology values as F stability class and wind speeds of 1.3 meters per second (for a ground level release). This results in an apparent discrepancy between the two approaches. In this report we document the differing modeling assumptions, discuss the logic behind each modeling approach used to represent local roughness and plume meander conditions, and numerically compare dispersion coefficients and concentration estimates obtained from each set of model inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Insights of Boundary Layer Turbulence Over the Complex Terrain of Central Himalaya from GVAX Field Campaign

Limited observations hinder understanding of turbulent characteristics in mountainous terrain resulting from heating or cooling of slopes, wind, vertical motions, and heat or moisture advection, which disperse aerosols and other pollutants over the region. In this study, the 1290 MHz radar wind profiler data are utilized to compute the boundary layer height (BLH), the refractive index structure constant (C n 2 ), and the energy dissipation rate (ε) over the central Himalayan site for the period of November 2011 to March 2012, from the intense Ganges Valley Aerosol Experiment (GVAX) field measurements. The radar wind profiler (RWP) based estimation of BLH and ε is validated against the radiosonde, representing the effectiveness of the datasets for further investigation. The strong seasonal variation of log C n 2 and log ε, with average values of ≈ -12 m -2/3 and -2 m 2 s -3 , respectively, is associated with the mountain-induced local circulations and stability in the atmospheric boundary layer. The weak stratification during weak flow is found to be responsible for deep mixing, particularly in the nocturnal boundary layer in spring. Furthermore, the level of cloud cover significantly impacts the strength of turbulence, with the highest cloud cover resulting in a substantial increase in log C n 2 (approximately -11 m -2/3 ) due to intense updraft and downdraft motions compared to clear skies. Additionally, the distribution of aerosol loading across the site, coupled with the behavior of BLH, atmospheric stability, and orographic-induced circulations, implies distinctive seasonal mechanisms for transporting aerosols toward the mountains. This study offers valuable insights into the diurnal and seasonal patterns of turbulent mixing and the mechanisms behind the transport of pollutants through boundary layer processes over the region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Large eddy simulations of offshore wind turbine wakes for two floating platform types

The growing prospect for large farms of floating offshore wind turbines requires a better understanding of wake effects for floating turbines. In this work, large eddy simulations with an actuator line model are used to study the wake of the NREL 5 MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCwind semi-submersible platforms. The simulations are carried out in the Simulator fOr Wind Farm Applications (SOWFA) coupled with OpenFAST for the platform and turbine motion. The wake location, deficit, and turbulence levels are compared for the two floating platforms and equivalent fixed-turbine cases. The effects of neutral versus stable atmospheric conditions are also compared. Most notably, floating-turbine wakes are deflected upwards compared to fixed-turbine wakes, because of mean platform pitch. The spar wake deflects upwards more than the semi-submersible, while the stable atmosphere increases this vertical deflection compared to the neutral. The time-varying rotor motions do not significantly affect the mid-to-far wake, though the stable atmosphere shows larger fixed-floating differences in horizontal wake fluctuations.

17 WIND ENERGY↗
DOE OSTI2023

Surface impacts and associated mechanisms of a moisture intrusion into the Arctic observed in mid-April 2020 during MOSAiC

Distinct events of warm and moist air intrusions (WAIs) from mid-latitudes have pronounced impacts on the Arctic climate system. We present a detailed analysis of a record-breaking WAI observed during the MOSAiC expedition in mid-April 2020. By combining Eulerian with Lagrangian frameworks and using simulations across different scales, we investigate aspects of air mass transformations via cloud processes and quantify related surface impacts. The WAI is characterized by two distinct pathways, Siberian and Atlantic. A moist static energy transport across the Arctic Circle above the climatological 90th percentile is found. Observations at research vessel Polarstern show a transition from radiatively clear to cloudy state with significant precipitation and a positive surface energy balance (SEB), i.e., surface warming. WAI air parcels reach Polarstern first near the tropopause, and only 1–2 days later at lower altitudes. In the 5 days prior to the event, latent heat release during cloud formation triggers maximum diabatic heating rates in excess of 20 K d -1 . For some poleward drifting air parcels, this facilitates strong ascent by up to 9 km. Based on model experiments, we explore the role of two key cloud-determining factors. First, we test the role moisture availability by reducing lateral moisture inflow during the WAI by 30%. This does not significantly affect the liquid water path, and therefore the SEB, in the central Arctic. The cause are counteracting mechanisms of cloud formation and precipitation along the trajectory. Second, we test the impact of increasing Cloud Condensation Nuclei concentrations from 10 to 1,000 cm -3 (pristine Arctic to highly polluted), which enhances cloud water content. Resulting stronger longwave cooling at cloud top makes entrainment more efficient and deepens the atmospheric boundary layer. Finally, we show the strongly positive effect of the WAI on the SEB. This is mainly driven by turbulent heat fluxes over the ocean, but radiation over sea ice. The WAI also contributes a large fraction to precipitation in the Arctic, reaching 30% of total precipitation in a 9-day period at the MOSAiC site. However, measured precipitation varies substantially between different platforms. Therefore, estimates of total precipitation are subject to considerable observational uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Role of Ocean and Atmosphere Variability in Scale-Dependent Thermodynamic Air-Sea Interactions

This study investigates the influence of oceanic and atmospheric processes in extratropical thermodynamic air-sea interactions resolved by satellite observations (OBS) and by two climate model simulations run with eddy-resolving high-resolution (HR) and eddy-parameterized low-resolution (LR) ocean components. Here, in this work, spectral methods are used to characterize the sea surface temperature (SST) and turbulent heat flux (THF) variability and co-variability over scales between 50 and 10,000 km and 60 days to 80 years in the Pacific Ocean. The relative roles of the ocean and atmosphere are interpreted using a stochastic upper-ocean temperature evolution model forced by noise terms representing intrinsic variability in each medium, defined using climate model data to produce realistic rather than white spectral power density distributions. The analysis of all datasets shows that the atmosphere dominates the SST and THF variability over zonal wavelengths larger than ~2,000–2,500 km. In HR and OBS, ocean processes dominate the variability of both quantities at scales smaller than the atmospheric first internal Rossby radius of deformation (R 1 , ~600–2,000 km) due to a substantial ocean forcing coinciding with a weaker atmospheric modulation of THF (and consequently of SST) than at larger scales. The ocean forcing also induces oscillations in SST and THF with periods ranging from intraseasonal to multidecadal, reflecting a red spectrum response to ocean forcing similar to that driven by atmospheric forcing. Such features are virtually absent in LR due to a weaker ocean forcing relative to HR.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗