skip to main content

Title: Turbulent Transport of Spray Droplets in the Vicinity of Moving Surface Waves

A common technique for estimating the sea surface generation functions of spray and aerosols is the so-called flux–profile method, where fixed-height concentration measurements are used to infer fluxes at the surface by assuming a form of the concentration profile. At its simplest, this method assumes a balance between spray emission and deposition, and under these conditions the concentration profile follows a power-law shape. It is the purpose of this work to evaluate the influence of waves on this power-law theory, as well as investigate its applicability over a range of droplet sizes. Large-eddy simulations combined with Lagrangian droplet tracking are used to resolve the turbulent transport of spray droplets over moving, monochromatic waves at the lower surface. The wave age and the droplet diameter are varied, and it is found that droplets are highly influenced both by their inertia (i.e., their inability to travel exactly with fluid streamlines) and the wave-induced turbulence. Deviations of the vertical concentration profiles from the power-law theory are found at all wave ages and for large droplets. The dynamics of droplets within the wave boundary layer alter their net vertical fluxes, and as a result, estimates of surface emission based on the flux–profile method can yield significant errors. In practice, the resulting implication is that the flux–profile method may unsuitable for large droplets, and the combined effect of inertia and wave-induced turbulence is responsible for the continued spread in their surface source estimates.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
American Meteorological Society
Date Published:
Journal Name:
Journal of Physical Oceanography
Page Range / eLocation ID:
p. 1789-1807
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We provide a first-principles analysis of the energy fluxes in the oceanic internal wave field. The resulting formula is remarkably similar to the renowned phenomenological formula for the turbulent dissipation rate in the ocean, which is known as the finescale parameterization. The prediction is based on the wave turbulence theory of internal gravity waves and on a new methodology devised for the computation of the associated energy fluxes. In the standard spectral representation of the wave energy density, in the two-dimensional vertical wavenumber–frequency (mω) domain, the energy fluxes associated with the steady state are found to be directed downscale in both coordinates, closely matching the finescale parameterization formula in functional form and in magnitude. These energy transfers are composed of a “local” and a “scale-separated” contributions; while the former is quantified numerically, the latter is dominated by the induced diffusion process and is amenable to analytical treatment. Contrary to previous results indicating an inverse energy cascade from high frequency to low, at odds with observations, our analysis of all nonzero coefficients of the diffusion tensor predicts a direct energy cascade. Moreover, by the same analysis fundamental spectra that had been deemed “no-flux” solutions are reinstated to the status of “constant-downscale-flux” solutions. This is consequential for an understanding of energy fluxes, sources, and sinks that fits in the observational paradigm of the finescale parameterization, solving at once two long-standing paradoxes that had earned the name of “oceanic ultraviolet catastrophe.”

    Significance Statement

    The global circulation models cannot resolve the scales of the oceanic internal waves. The finescale parameterization of turbulent dissipation, a formula grounded in observations, is the standard tool by which the energy transfers due to internal waves are incorporated in the global models. Here, we provide an interpretation of this parameterization formula building on the first-principles statistical theory describing energy transfers between waves at different scales. Our result is in agreement with the finescale parameterization and points out a large contribution to the energy fluxes due to a type of wave interactions (local) usually disregarded. Moreover, the theory on which the traditional understanding of the parameterization is mainly built, a “diffusion approximation,” is known to be partly in contradiction with observations. We put forward a solution to this problem, visualized by means of “streamlines” that improve the intuition of the direction of the energy cascade.

    more » « less
  2. Abstract

    Quantifying the influence of sea spray on air‐sea fluxes under high‐wind conditions is challenging due to a variety of factors. Among existing models, the so‐called bulk air‐sea flux model is commonly used in meteorological applications due to its simplicity, which often involves strong but untested assumptions on spray‐mediated heat fluxes and feedback effects. For example, a common assumption is to treat each droplet size as an independent contribution; that is, it does not interact with other sizes. Thus, the interactions between different size classes of spray are often neglected. In this study, we focus on the polydispersity of the spray size distribution and investigate the appropriateness of assuming an independent contribution from different spray size classes. We implement direct numerical simulations (DNS) with Lagrangian tracking of spray droplets. Based on DNS results, the bulk spray model fails to capture the interactions between different sizes that are observed directly from the droplet and feedback statistics in DNS. Thus, assuming independent contributions from spray droplets results in significant overestimates on the total heat fluxes. We further test different representative sizes of a spray size distribution. We find that the volume‐weighted representative size is capable of predicting the droplet‐modified temperature and humidity fields and generally captures the vertical profiles of spray‐mediated and interfacial heat fluxes. The results indicate that the computation of spray‐mediated fluxes can be simplified in large‐scale parameterizations.

    more » « less
  3. Abstract The Eulerian multifluid mathematical model is developed to describe the marine atmospheric boundary layer laden with sea spray under the high-wind condition of a hurricane. The model considers spray and air as separate continuous interacting turbulent media and employs the multifluid E –ϵ closure. Each phase is described by its own set of coupled conservation equations and characterized by its own velocity. Such an approach enables us to accurately quantify the interaction between spray and air and pinpoint the effect of spray on the vertical momentum transport much more precisely than could be done with traditional mixture-type approaches. The model consistently quantifies the effect of spray inertia and the suppression of air turbulence due to two different mechanisms: the turbulence attenuation, which results from the inability of spray droplets to fully follow turbulent fluctuations, and the vertical transport of spray against the gravity by turbulent eddies. The results of numerical and asymptotic analyses show that the turbulence suppression by spray overpowers its inertia several meters above wave crests, resulting in a noticeable wind acceleration and the corresponding reduction of the drag coefficient from the reference values for a spray-free atmosphere. This occurs at much lower than predicted previously spray volume fraction values of ∼10 −5 . The falloff of the drag coefficient from its reference values is more strongly pronounced at higher altitudes. The drag coefficient reaches its maximum at spray volume fraction values of ∼10 −4 , which is several times smaller than predicted by mixture-type models. 
    more » « less
  4. Abstract The quantification of pressure fields in the airflow over water waves is fundamental for understanding the coupling of the atmosphere and the ocean. The relationship between the pressure field, and the water surface slope and velocity, are crucial in setting the fluxes of momentum and energy. However, quantifying these fluxes is hampered by difficulties in measuring pressure fields at the wavy air-water interface. Here we utilise results from laboratory experiments of wind-driven surface waves. The data consist of particle image velocimetry of the airflow combined with laser-induced fluorescence of the water surface. These data were then used to develop a pressure field reconstruction technique based on solving a pressure Poisson equation in the airflow above water waves. The results allow for independent quantification of both the viscous stress and pressure-induced form drag components of the momentum flux. Comparison of these with an independent bulk estimate of the total momentum flux (based on law-of-the-wall theory) shows that the momentum budget is closed to within approximately 5%. In the partitioning of the momentum flux between viscous and pressure drag components, we find a greater influence of form drag at high wind speeds and wave slopes. An analysis of the various approximations and assumptions made in the pressure reconstruction, along with the corresponding sources of error, is also presented. 
    more » « less
  5. Abstract Upper-ocean turbulence is central to the exchanges of heat, momentum, and gasses across the air/sea interface, and therefore plays a large role in weather and climate. Current understanding of upper-ocean mixing is lacking, often leading models to misrepresent mixed-layer depths and sea surface temperature. In part, progress has been limited due to the difficulty of measuring turbulence from fixed moorings which can simultaneously measure surface fluxes and upper-ocean stratification over long time periods. Here we introduce a direct wavenumber method for measuring Turbulent Kinetic Energy (TKE) dissipation rates, ϵ , from long-enduring moorings using pulse-coherent ADCPs. We discuss optimal programming of the ADCPs, a robust mechanical design for use on a mooring to maximize data return, and data processing techniques including phase-ambiguity unwrapping, spectral analysis, and a correction for instrument response. The method was used in the Salinity Processes Upper-ocean Regional Study (SPURS) to collect two year-long data sets. We find the mooring-derived TKE dissipation rates compare favorably to estimates made nearby from a microstructure shear probe mounted to a glider during its two separate two-week missions for (10 −8 ) ≤ ϵ ≤ (10 −5 ) m 2 s −3 . Periods of disagreement between turbulence estimates from the two platforms coincide with differences in vertical temperature profiles, which may indicate that barrier layers can substantially modulate upper-ocean turbulence over horizontal scales of 1-10 km. We also find that dissipation estimates from two different moorings at 12.5 m, and at 7 m are in agreement with the surface buoyancy flux during periods of strong nighttime convection, consistent with classic boundary layer theory. 
    more » « less