the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The impact of entrained air on ocean waves
Juan M. Restrepo
Alex Ayet
Luigi Cavaleri
We make a physical–mathematical analysis of the implications that the presence of a large number of tiny bubbles may have, when present, on the thin upper layer of the sea. In our oceanographic example, the bubbles are due to intense rain. It was found that the bubbles increase momentum dissipation in the near surface and affect the surface tension force. For short waves, the implications of increased vorticity are momentum exchanges between wave and mean flow and modifications to the wave dispersion relation. For the direct effect we have analyzed, the implications are estimated to be nonsignificant when compared to other processes of the ocean. However, we hint at the possibility that our analysis may be useful in other areas of research or practical application.
Our interest is in the consequences of intense rain on the smallscale roughness of the ocean surface. More specifically, we analyze how the presence of tiny, raininduced air bubbles in the thin upper layer of the ocean affects the highfrequency tail of the gravity wave spectrum, including capillary waves. Although tiny, both in length and height, with respect to the majestic motion of the surface, these waves are important for both physical and practical reasons. Physically, they control a large part of the exchanges between ocean and atmosphere, and they act as roughness elements for atmospheric turbulence and modulate the interaction between longer waves, currents, and wind (e.g., Ayet et al., 2020). At the same time, their presence is what allows the scatterometers to measure, globally, the winds on the ocean (Kudryavtsev et al., 1999).
The disappearance of capillary sea waves with rain is not new. Known to mariners since their early attempts (reduced smallscale roughness leads also to reduced breaking), the subject has also been studied in recent times by, e.g., Le Méhauté and Khangaonkar (1990) and Cavaleri et al. (2018). The former authors focused on the mechanical action of rain on the dynamics of the wave orbital motion. The latter authors' attention, taking the wavelets attenuation as a de facto evidence, focused more on the implications for air–sea interactions. In this work, we try to clarify one of the reasons why rain attenuates surface wavelets. Granted, if strong enough, the raindrops' mechanical action and induced turbulence, we focus our attention on the raininduced presence of a large number of tiny air bubbles just below the surface (a few centimeters) and on how they affect the wave motion. To avoid the dominance of other processes, we exclude stormy conditions, considering either only swell or calm sea cases, but with enough wind to lead to waves on the surface. As for rain, just to frame the order of magnitude, under an intense rain with, e.g., 100 mm h^{−1}, assuming 4 mm diameter raindrops (volume $\mathrm{3}\times {\mathrm{10}}^{\mathrm{6}}$ m^{3}), this translates into 100 drops per square meter per second. Wolf (2001) suggests that, under these conditions, air injected into the sea can be a few cubic centimeters per square meter per second. Compared to other mechanisms that introduce air into the sea, one can conclude that rain contribution is very small. In general, this is true of quantitative terms, and we consider the process interesting from the physical point of view and with potential application in other practical fields.
On this basis, we make a physical–mathematical analysis of the ensuing processes. In Sect. 2, we propose a model that suggests how the entrainment of air leads to an increase in the effective viscosity of the upper ocean. We then use results from ocean acoustics (see Oguz and Prosperetti, 1990 and Prosperetti and Oguz, 1993) to compute the distribution of air bubbles for a given rain rate (Sect. 3). How changes of effective viscosity affect gravity waves is taken up in Sect. 4. It is found that, under vigorous rains, the air concentration due to air bubbles near the surface is capable of producing a damping effect akin to the (idealized) contamination of a free surface. A discussion of the results and conclusions and possible extensions of the present work appear in Sects. 5 and 6.
To determine the impact of air bubbles on the thin upper level of the ocean, we perform homogenization on the Navier–Stokes equations. Homogenization is a wellestablished technique in modeling transport in complex media, endowed with statistical homogeneity in the media, as viewed at large scales (see Babuška, 1976; Bensousan et al., 1978; Cioranescu and Donato, 1999).
We assume that air bubbles are distributed uniformly in the transverse direction within the vicinity of the ocean surface. For our specific example, we explore how the presence of bubbles affects the effective viscosity of the fluid averaged over a cell Ω, of size ℓ^{3} (subwave scale), over which the distributions of the density and viscosity of the combined water and air bubble mixture are statistically stationary. The typical subwave speed is u_{Ω}. The ratio of the bubble radii to the averaging ℓ defines, for us, a small parameter ϵ≪1.
Velocity and position are denoted by $\mathit{u}=(u,v,w)$ and $\mathit{r}=(x,y,z)$, respectively. The free surface is denoted by η, and z=0 corresponds to the quiescent sea level. Gravity is $g\widehat{\mathit{z}}$, and the vector $\widehat{\mathit{z}}$ points upwards, along increasing z. Velocity is scaled by u_{Ω}, length by ℓ, time by $\mathrm{\ell}/{u}_{\mathrm{\Omega}}$, density by ρ_{w}, and the density of water pressure by ${\mathit{\rho}}_{\mathrm{w}}{u}_{\mathrm{\Omega}}^{\mathrm{2}}$. We then define a Reynolds number $\mathit{\alpha}={u}_{\mathrm{\Omega}}\mathrm{\ell}/{\mathit{\nu}}_{\mathrm{w}}$, where ν_{w} is the kinematic viscosity of water. We also define a Froude number $\mathrm{1}/\sqrt{\mathit{\gamma}}={u}_{\mathrm{\Omega}}/\sqrt{g\mathrm{\ell}}$. The scaling leads to the following:
where Π is pressure and Ξ the stress tensor $\mathrm{\Xi}=\mathbf{\nabla}\mathit{u}+[\mathbf{\nabla}\mathit{u}{]}^{\top}$. D is the proportionality tensor.
Let $\mathbf{R}=(X,Y,Z)$ be the largescale space variable, such that $\mathbf{\nabla}=\mathbf{\nabla}+\mathit{\u03f5}{\mathbf{\nabla}}_{\mathbf{R}}$, and assume slow time ${\partial}_{t}=\mathit{\delta}{\partial}_{T}$. Also, assume that $\mathrm{\Pi}(\mathit{r},\mathbf{R},T)=pz+{\mathit{\u03f5}}^{\mathrm{2}}{p}_{\mathrm{0}}+{\mathit{\u03f5}}^{\mathrm{4}}{p}_{\mathrm{1}}+\mathrm{\dots}$, $\mathit{u}=\mathit{\u03f5}({\mathit{u}}_{\mathrm{0}}+\mathit{\u03f5}{\mathit{u}}_{\mathrm{1}}+{\mathit{\u03f5}}^{\mathrm{2}}{\mathit{u}}_{\mathrm{2}}\mathrm{\dots})$, and $\mathit{\eta}=\mathit{\alpha}({\mathit{\eta}}_{\mathrm{0}}+\mathit{\u03f5}{\mathit{\eta}}_{\mathrm{1}}+{\mathit{\u03f5}}^{\mathrm{2}}{\mathit{\eta}}_{\mathrm{2}}\mathrm{\dots})$. The orders are α=𝒪(ϵ^{2}), δ=𝒪(ϵ), and $\mathit{\gamma}=\mathcal{O}\left({\mathit{\u03f5}}^{\mathrm{1}}\right)$.
In the following, we derive the fluid mechanics equations for averaged dynamic quantities (see the related work of Caflisch et al., 1985), appropriate at wave spatiotemporal scales, for example. Our goal is to find effective equations valid in a composite media, under the assumptions made above. The equations and the dynamic quantities will agree with the nonaveraged ones, when the media has a single species density ρ_{w} and the tensor D=ν_{w}δ_{ij}. Here we introduce the inhomogeneity as tiny air bubbles, uniformly distributed in the transverse direction in a matrix of water. The source of the bubbles is rain, and the connection between the bubble density in the near surface and the rain rate will be suggested in a later section. Anticipating some of the results, the momentum and continuity equations for the averaged dynamic quantities are similar to the point versions. However, material properties, such as the density and the tensor D are different for the averaged cells, quantitatively affecting the momentum predictions.
Collecting by orders in ϵ, the momentum equation reads, as follows:

𝒪(ϵ).
$$\mathbf{\nabla}\cdot \left[D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{0}}\right]=\mathrm{0}.$$At this microscale, pressure gradients in the cell are negligible, as are variations in velocity within the cell. We are also assuming stationary conditions. Upon integrating in r and invoking periodicity, it is clear that ${\mathrm{\Xi}}_{\mathrm{0}}={\mathrm{\Xi}}_{\mathrm{0}}(\mathbf{R},t,T)$ and, thus, ${\mathit{u}}_{\mathrm{0}}={\mathit{u}}_{\mathrm{0}}(\mathbf{R},t,T)$.

𝒪(ϵ^{2}).
$$\begin{array}{}\text{(2)}& \mathbf{\nabla}\cdot \left[D\right(\mathit{r}\left)\phantom{\rule{0.125em}{0ex}}{\mathrm{\Xi}}_{\mathrm{1}}\right]+\mathbf{\nabla}\cdot \left[D\right(\mathit{r}\left)\phantom{\rule{0.125em}{0ex}}{\mathrm{\Xi}}_{\mathrm{0}}\right]+{\mathbf{\nabla}}_{\mathbf{R}}\cdot \left[D\right(\mathit{r}\left)\phantom{\rule{0.125em}{0ex}}{\mathrm{\Xi}}_{\mathrm{0}}\right]=\mathrm{0}.\end{array}$$The last term above is zero (we made use of $\mathbf{\nabla}\cdot (\mathrm{\nabla}\mathrm{\Xi}{)}^{T}=\mathbf{\nabla}(\mathbf{\nabla}\cdot \mathrm{\Xi})$). Integration by parts of Eq. (2) and using periodicity, ${\mathrm{\Xi}}_{\mathrm{1}}(\mathit{r},\mathbf{R},t,T)={\mathrm{\Xi}}_{\mathrm{0}}+{D}^{\mathrm{1}}\left(\mathit{r}\right)C(\mathbf{R},t,T)$, where C is a tensor that is independent of r. Periodicity and integration in r over Ω of every differential term in Ξ_{1} implies that, in the following:
$$\begin{array}{}\text{(3)}& \mathrm{0}=\mathrm{\Omega}{\mathrm{\Xi}}_{\mathrm{0}}+\underset{\mathrm{\Omega}}{\int}{D}^{\mathrm{1}}\left(\mathit{r}\right)\mathrm{d}\mathit{r}\phantom{\rule{0.125em}{0ex}}C(\mathbf{R},t,T).\end{array}$$Solving for the tensor from this equation and integrating with respect to r, one obtains the following:
$$C(\mathbf{R},t,T)={\left[{\displaystyle \frac{\mathrm{1}}{\mathrm{\Omega}}}\underset{\mathrm{\Omega}}{\int}{D}^{\mathrm{1}}\left(\mathit{r}\right)\mathrm{d}\mathit{r}\right]}^{\mathrm{1}}{\mathrm{\Xi}}_{\mathrm{0}}.$$We define the tensor as follows:
$$\begin{array}{}\text{(4)}& \mathit{\chi}={\left[{\displaystyle \frac{\mathrm{1}}{\mathrm{\Omega}}}\underset{\mathrm{\Omega}}{\int}{D}^{\mathrm{1}}\left(\mathit{r}\right)\mathrm{d}\mathit{r}\right]}^{\mathrm{1}}.\end{array}$$Then, we returning to Eq. (2), as follows:
$$\begin{array}{}\text{(5)}& {\mathrm{\Xi}}_{\mathrm{1}}={\mathrm{\Xi}}_{\mathrm{0}}+{D}^{\mathrm{1}}\left(\mathit{r}\right)\mathit{\chi}{\mathrm{\Xi}}_{\mathrm{0}},\end{array}$$which prescribes the firstorder correction to the stress tensor, which could be obtained explicitly if the zeroth order is known (D^{−1}(r) is an observed quantity). Taking the divergence of Eq. (5), as follows:
$$\begin{array}{}\text{(6)}& {\mathbf{\nabla}}_{\mathbf{R}}\cdot \left(D{\mathrm{\Xi}}_{\mathrm{1}}\right)={\mathbf{\nabla}}_{\mathbf{R}}\cdot \left(D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{0}}\right)+{\mathbf{\nabla}}_{\mathbf{R}}\cdot \left(\mathit{\chi}{\mathrm{\Xi}}_{\mathrm{0}}\right).\end{array}$$ 
𝒪(ϵ^{3}). The momentum balance reads as follows:
$$\begin{array}{}\text{(7)}& \begin{array}{rl}{\displaystyle \frac{\partial {\mathit{u}}_{\mathrm{0}}}{\partial T}}+{\mathit{u}}_{\mathrm{0}}& \cdot {\mathbf{\nabla}}_{\mathbf{R}}{\mathit{u}}_{\mathrm{0}}+{\mathit{u}}_{\mathrm{0}}\cdot \mathbf{\nabla}{\mathit{u}}_{\mathrm{1}}+{\mathbf{\nabla}}_{\mathbf{R}}{p}_{\mathrm{0}}\\ & =+\mathbf{\nabla}\cdot \left[D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{1}}\right]+{\mathbf{\nabla}}_{\mathbf{R}}\cdot \left[D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{0}}\right]\\ & \phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}+\mathbf{\nabla}\cdot \left[D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{2}}\right]+{\mathbf{\nabla}}_{\mathbf{R}}\cdot \left[D\right(\mathit{r}\left){\mathrm{\Xi}}_{\mathrm{1}}\right].\end{array}\end{array}$$
Using Eq. (6) in the last term in Eq. (7) and averaging all quantities over Ω, (i.e., integrating over r using periodicity in the cells), leads to the following:
The homogenized incompressibility condition is as follows:
At this point, we have a set of (closed) mean field equations for momentum and mass conservation, valid at the wave scales (R,T). The mean field velocity and pressure are well defined, assuming homogeneity conditions at the small scale are valid. The effective tensor of proportionality in the stress tensor term χ is formally obtained by knowing the microscale composition, density distribution, and dynamic viscosities. In practice, it is obtained by examining the stress/strain relation for the composite fluid.
Since the above derivation is not explicitly connected to rain, one can imagine that the same approach applies to similar situations where, for whichever reason, a large number of bubbles is distributed in the enclosing medium. We will touch on this point further in the final discussion.
The tensor χ takes the value of the ocean, changing it when bubbles appear for whatever reason. In what follows, we limit ourselves to the simplest possible case, i.e., χ=ν_{w}δ_{ij}, when bubbles are not present, where δ_{ij} is the 3dimensional space Kronecker delta. For our specific example, we use Eq. (4), assuming a matrix consisting of a homogeneous concentration of air bubbles, with diffusion constant ν_{a}, in a background fluid. In this very simple case, the enhanced value of χ becomes the following:
where ${N}_{\mathit{\nu}}=\frac{{\mathit{\nu}}_{\mathrm{w}}}{{\mathit{\nu}}_{\mathrm{a}}}$, and Θ connotes the volumetric ratio of air to water. For a given Θ>0, N_{ν}<1 leads to an increased effective momentum dissipation, and when N_{ν}>1, the effective dissipation is lower. For air or oil bubbles, N_{ν}<1.
With a fixed ratio ${N}_{\mathit{\nu}}=\mathrm{1}/\mathrm{15}$, which is roughly the ratio of viscosities of water and air, Fig. 1 shows how the effective dissipation coefficient K changes as a function of the volumetric ratio Θ. The relationship is nearly linear for the small volume fractions of the oceanic case.
If we adopt the specific form of χ as in Eq. (10), the dissipation term in Eq. (8) becomes $K{\mathbf{\nabla}}_{\mathbf{R}}^{\mathrm{2}}{\mathit{u}}_{\mathrm{0}}$ when bubbles are present and ${\mathit{\nu}}_{\mathrm{w}}{\mathbf{\nabla}}_{\mathbf{R}}^{\mathrm{2}}{\mathit{u}}_{\mathrm{0}}$, when there are no air bubbles. Equation (10) leads us to conclude that the effective momentum dissipation of water with air bubbles will be different from (actually larger than) ν_{w}. Physically, there is an increase in the rotational component of the flow and an associated increase in energy dissipation.
In this section, we assume a statistical distribution for raindrops and derive the consequent density of air bubbles in the ocean surface layer.
3.1 Rain distribution
The density of raindrops is assumed to follow the Marshall–Palmer distribution (see Marshall and Palmer, 1948). The density N_{r} of raindrops of radius r (meters) per unit volume (“number” m^{−3} m^{−1}) is as follows:
where ${N}_{\mathrm{0}}=\mathrm{8}\times {\mathrm{10}}^{\mathrm{6}}$ m^{−3}m^{−1}, $\mathrm{\Lambda}\left(R\right)=\mathrm{8.2}{R}^{\mathrm{0.21}}\times {\mathrm{10}}^{\mathrm{3}}$ (meter), and R (millimeters per hour) is the rain rate. Using N_{r}, we can then define the drop rate density (DRD), which describes the rate of falling raindrops of radius r per surface area (“number” m^{−2} s^{−1} m^{−1}). Namely, in the following:
where w_{r} is the terminal velocity of drops in the air (see Medwin et al., 1992). The terminal velocity is computed following the thirdorder polynomial estimates of Dingle and Lee (1972).
3.2 Bubble production
To complete our quantitative example, we need to relate the number of bubbles in the upper layer of the ocean to rain. Falling raindrops generate subsurface air bubbles in the neighborhood of the sea surface (see Prosperetti and Oguz, 1993 for a review). Small raindrops produce small air bubbles with a very narrow distribution of radii. Intermediate raindrops, with a radius between 0.55–1.1 mm, were shown not to produce air bubbles (see Medwin et al., 1992). For these, the impinging raindrops do not have the kinetic energy necessary to produce the requisite conical crater and jetting of the sea surface that engulfs air. Large raindrops, with radii larger than 1.1 mm, create a crater and a canopy on the sea surface, which, by collapsing, produces a downward liquid jet at the bottom of the crater, followed by the generation of an air bubble. At high rain rates, the larger raindrops are responsible for the bulk of the gas injection, creating bubbles with a varied distribution of radii.
The production of air bubbles by rain can be measured by acoustic means (see Prosperetti and Oguz, 1993). Oguz and Prosperetti (1990) classify air bubble production by falling raindrops in two regimes. Small raindrops, of radii between 0.41 and 0.55 mm, create type I air bubbles. These have a radius of 0.22 mm and are created at the bottom of a conical crater created on the ocean surface by the impinging raindrop. These air bubbles have a narrow acoustic spectrum with a spectral peak at 14 Hz. The acoustic spectrum for these type I air bubbles was found to be insensitive to the rain rate. Medwin et al. (1990) observed that, when rain hits the ocean surface at an angle, owing to strong winds or very steep wave conditions, the acoustic spectrum peak shifts downward, and there is a broadening of the spectrum at higher frequencies. The type I air bubbles do not contribute significantly to total submerged gas volume, and their contribution will be ignored in what follows.
Raindrops with a radius greater than 1.1 mm produce type II air resonant bubbles of varying radius. The relationship between the raindrop radius r (meters) and the peak acoustic emission frequency f_{0} (kilohertz) due to trapped air bubbles was empirically determined (see Medwin et al., 1992) as follows:
The relation between the nearsurface air bubble radius a (meters) and the peak acoustic emission follows from the Rayleigh–Plesset equation (Leighton, 1994, p. 306). For the large bubbles considered here, it can be simplified to Minnaert's formula (see Minnaert, 1933; Plesset and Prosperetti, 1977) as follows:
where γ=1.4 is the ratio of the specific heat of the bubble gas, P is the ambient pressure (surrounding the bubble), and ρ_{0}=1030 kg m^{−3} is the density of sea water. Close to the surface, the ambient pressure is approximately the surface pressure $P=\mathrm{1.01}\times {\mathrm{10}}^{\mathrm{5}}$ Pa and Eq. (14), further simplified as ${f}_{\mathrm{0}}=\mathrm{3.25}\times {\mathrm{10}}^{\mathrm{3}}/a$ (see Medwin and Clay, 1997).
The relation between the entrained air bubble radius a(r) (meters) and the incident raindrop radius r, for large raindrops, thus, reads, for $r>\mathrm{1.1}\times {\mathrm{10}}^{\mathrm{3}}$ m, as follows:
As an example, raindrops of radius 1.1–2.3 mm produce type II air bubbles of radius 0.2–1.3 mm, respectively.
Small raindrops (r<0.55 mm) almost always produce type I air bubbles. However, this is not the case for larger raindrops. The distribution of the number of air bubbles produced, as a function of the raindrop radius r, was found by Medwin et al. (1992) and is depicted in Fig. 2 as a probability distribution.
3.3 Bubble distribution in the wave boundary layer
Tying back to Sect. 2, by finding the bubble distribution in the wave boundary layer, we can estimate the dimensionless volume fraction Θ. More precisely, the volume fraction reads $\mathrm{\Theta}=\int V\left(a\right)N\left(a\right)\mathrm{d}a$ for bubbles of radius a, with density N(a) (“number” m^{−3} m^{−1}) that occupy a volume $V\left(a\right)=(\mathrm{4}/\mathrm{3})\mathit{\pi}{a}^{\mathrm{3}}$. The aim of the final part of this section is to link Θ to the rain precipitation rate.
In doing so, we will make some approximations. First, bubbles are injected at a very limited depth by raindrops, and this depth may vary with bubble size (see Ho et al., 2000). Second, a homogeneous air bubble distribution in the layer between the depth of injection and the surface is an approximation, as the accumulation of bubbles near the surface can be highly complex due to buoyant forces, damping, and the background turbulence (see, e.g., Merlivat and Memerly, 1983). We will specify how this is taken into account.
The bubble distribution N(a) can be described by an advection–diffusion equation (see Woolf and Thorpe, 1991 and references therein and, for a similar case for oil drops, Moghimi et al., 2018). This equation describes a balance between the advection of bubbles by the vertical fluid velocity, the attenuation of bubble density (i.e., diffusion, which includes the shrinking of bubble radius during their lifetime; see Merlivat and Memerly, 1983), the production of bubbles by rain, and a sink term that captures the loss of air bubbles bursting at the surface or being dispersed by oceanic turbulence.
Solutions of the advection–diffusion equation necessitate the specification of the velocity and the dispersion, generating nonstationary descriptions that might also include mixing due to wave turbulence (see Restrepo et al., 2015; Moghimi et al., 2018).
In the following, the advection–diffusion equation is simplified to estimate the air bubble concentration in a thin layer close to the surface, similar to those in Keeling (1993). First, we assume that the air bubble concentration can be linearly related to a bulk bubble distribution (i.e., averaged over the thin layer), which we, thus, define as c_{e}N(a). Second, a steady state is considered in which the bulk bubble distribution follows from a balance between the incoming bubble flux due to rain and the outgoing bubble flux due to upward bubble advection and bursting at the surface, as follows:
where P_{r}(r) is the probability of the production of an air bubble by a raindrop of radius r (see Fig. 2), and W(a) is the upward vertical ascent speed of bubbles of radius a. From this equation, the bubble distribution reads $N\left(a\right)={w}_{r}{N}_{r}(r,R){P}_{r}/\left[{c}_{e}W\right(a\left)\right]$. We assume, as an estimate for the bubble upward speed W(a), the following form:
where $\mathit{\tau}=\mathrm{72.8}\times {\mathrm{10}}^{\mathrm{3}}$ N m^{−1} is the air–sea surface tension (see Clift et al., 1978, p. 172). This form, predicted in Mendelson (1967), is valid for ellipsoidal bubbles with a radius greater than 0.65 mm and whose upward path is not rectilinear.
Referring to the above approximations, the factor c_{e} can first be interpreted as a depth scaling parameter accounting for differences between the bulk estimate and the nonhomogeneous bubble distribution N. A further approximation is that the effective vertical velocity of bubbles can be significantly lower than Eq. (17) by up to 60 % due to the presence of contaminants at the bubble surface (see Clift et al., 1978; Fig. 7.3). Hence, the factor c_{e} is also to be interpreted as a decrease in the bubble vertical velocity, which then reads c_{e}W(a).
Making use of Eq. (12), $N\left(a\right)=(\text{DRD}\times {P}_{r})\left(r\right)/{c}_{e}W\left(a\right)$. The volumetric ratio $\mathrm{\Theta}=\int V\left(a\right)N\left(a\right)\mathrm{d}a$, for bubbles of radius a with density N(a) that occupy volume $V\left(a\right)=(\mathrm{4}/\mathrm{3})\mathit{\pi}{a}^{\mathrm{3}}$ then reads as follows:
where the limits of integration encompass the volumetric contributions of type II bubbles. The obtained values of volume fraction are consistent with data from the laboratory experiments of Ho et al. (2000) (see their Fig. 7), which suggests that, for rain rates of 114 mm per hour, the volume fraction near the air–water interface is between 10^{−6} and 10^{−5}, corresponding to c_{e} of the order of 0.1.
The subject of mathematical models for wave dissipation that reproduce the observed dissipation rates observed in laboratory experiments is considered in the study by Henderson et al. (2015). The focus is on the dissipation effects of small waves due to surface contamination with air above it. Several wave damping models are compared to data. In this work, a model by Jenkins and Jacobs (1997) is compared to experimental data of the damping of small waves. The model posits a small (boundary) layer sitting over a deep water layer. The twolayer model yields expressions valid under several limits. Our work can contribute to this work by defining how the equations in the upper layer are modified by the presence of air bubbles and how their density relates to rain, if this is the mechanism for their generation.
In what follows, we will invoke a far less general, but more consistent, way of addressing the question on how the presence of rain bubbles affects water waves. The derivation of infinitesimal amplitude gravity waves dynamics, starting from Navier–Stokes, appears in Lamb (1926), article 349. We revert to dimensional quantities in what follows. The solution u_{0} satisfies the linear version of Eqs. (8) and (9) and the two stress conditions at the surface. Namely, in the following:
where ϕ_{0} is the mean field velocity potential, and ρ_{*} is the density of water (no rain) or the reduced density (i.e., raining, with air bubbles in the water), and time (at wave scales) is T. At some distance ℓ below the surface, the conditions are as follows:
By combining the boundary and nonboundary solutions, the approximate solution of the system, assuming a vanishingly small upper layer thickness and vanishing velocity at depth, is as follows:
where ψ is the streamfunction, and σ is the angular frequency. The solution (23) represents infinitesimal progressive waves traveling in the X direction.
The dispersion relation for the waves is as follows:
showing that changes in the surface tension due to the presence of air bubbles injected by the rain impacts the dispersion relation. The effective density is as follows:
which, in our simplified accounting, would be ${\mathit{\rho}}^{*}=\mathrm{\Theta}{\mathit{\rho}}_{\mathrm{a}}+(\mathrm{1}\mathrm{\Theta}){\mathit{\rho}}_{\mathrm{w}}$.
The effective wave dissipation 2K appears in the exponential factor exp [−2Kk^{2}T]. When the rainstorm is sufficiently intense (but the raindrops are not too large), the wave dissipation 2K will change from 2ν_{w} to a higher value, depending on how much air is injected by the impinging raindrops.
The dependence of the effective wave dissipation on the rain rate, using Eqs. (10) and (18), is depicted in Fig. 3. The surface tension effect increases as well, affecting the dispersion relation. These effects are manifest primarily in capillary, highfrequency wave components. The bubble model derived in Sect. 3 has one free parameter, c_{e}, which controls the upward flux of bubbles. The solid line shows the wave dissipation for c_{e}=1, i.e., the situation for which the bulk estimation of the upward flux of bubbles is assumed to be correct. The dashed line shows the wave dissipation corresponding to the more realistic value, c_{e}=0.1, which matches the measured volume fraction of Ho et al. (2000). With respect to c_{e}=1, it corresponds to a situation where the residence time of the bubbles in the upper ocean layer is increased, due to, e.g., nearsurface turbulence.
In both cases, the increase in effective wave damping due to the injection of air, at least as estimated by simple considerations, is small compared to the values reported in the experiments in Peirson et al. (2013) (and in Tsimplis, 1992), further proof of what was stated in the introduction – i.e., that, for large rain rates and, as in the wave tank experiments, large raindrops, the mechanical effect is dominant for the immediate, albeit limited, attenuation of water waves. In these instances, using very high rain rates, the effective dissipation was found to increase by 3–10 times when rain is present, compared to without rain.
Figure 3 also shows the sensitivity of the model to changes in ${\mathit{\nu}}_{\mathrm{w}}/{\mathit{\nu}}_{\mathrm{a}}$, which is much lower than to changes in the bubbles volume fraction, i.e., the tuning parameter c_{e}. The gray shading shows the intervals obtained by varying the air viscosity inside air bubbles (ν_{a}) from 10^{−5} to 10^{−4} (the lower and upper parts of the shading, respectively) for a fixed c_{e} of 0.1. This range of variation in ν_{a} corresponds to variations in the effective kinematic viscosity of air due to the mechanical effect of raindrops, as observed in Harrison and Veron (2017).
Stimulated by an oceanographic problem, we have developed a physical–mathematical approach to analyze how the presence of a large number of small bubbles affects the characteristics of the containing liquid. Given the general method and having specified the necessary conditions for its application, we focused on a specific example, namely the presence of air bubbles in the thin upper layer of the ocean, and their impact on ocean surface waves. We first derived the (upscaled) equations for mass and momentum conservation. The upscaled velocity and pressure are averages over subwave scales over which the material with heterogeneities assumes a stationary distribution. These equations revealed that the effective (upscaled) momentum dissipation is increased when bubbles are present with respect to the nobubble case.
Under homogenization conditions, the model yields a dynamic equation that has a higher effective diffusion constant with bubbles than without bubbles. This translates into an increased kinematic viscosity of the air–water mixture, which is directly related to an increase in the damping of surface waves. The source for the tiny uniformly distributed bubbles is, in our case, rain, whose small drops generate the bubble distribution (large drops lead to splashing and turbulence in the upper layer, which is a different process). For typical rain rates, we found that the resulting wave damping due to bubbles is not significant with respect to other damping sources, such as turbulenceinduced dissipation due to the impact of large raindrops. We also found that the presence of bubbles affects the capillary waves' dispersion relation through a change in the effective surface tension felt by the waves.
The increase in wave damping (related to kinematic viscosity) when bubbles are present may appear counterintuitive at first sight, but it results from kinematic viscosity being the ratio of the dynamic viscosity and the density of the fluid of interest. As bubbles are injected into water, the homogenized kinematic viscosity decreases, but so does the homogenized fluid density, with the overall effect being an increase in the homogenized kinematic viscosity. Physically, the presence of small bubbles at the typically small void ratios of nonbreaking waves enhances the momentum transfer between the waves and the mean flow. The enhanced wave rotational component results in higher viscous effects. Moreover, since the density enters explicitly in the definition of the surface tension force, when the volume fraction of air to water is not zero, then the average density of the medium decreases, and, thus, with this decrease comes an increase in surface tension forces. Finally, the viscous effects and changes in the surface tension also lead to changes in the dispersion relation and, as a consequence, in the wave group velocity as well.
Many simplifications were made in the formulation of the rain/wave model. In the momentum equation, we used the most simplistic model for the tensor D(r); the consequence is meek momentum exchanges between the rotational and irrotational components of the velocity, and this argues for a more vigorous research effort in determining more realistic models for the tensor. In particular, the tensor should account for all sources of smallscale heterogeneity besides air bubbles. The model for the mechanism that ties the rain, the bubble presence in the layer, and the momentum and mass conservation equations below the free surface has three critical parameters, namely the rain rate R, the airtowater volume fraction Θ (through the flux strength parameter, c_{e}, associated with the upward flux of air bubbles from their injection depth to the surface), and the ratio of diffusivities of water and air, N_{ν}. Given those approximations, the basic result is that the volume fraction of air to water would have to be exceptionally large for the effect to be significant when compared to other damping mechanisms. However, granted the limited effect in the analyzed case, larger amounts of bubbles or, mutatis mutandis, of buoyant organic or hydrocarbons are certainly possible. In these cases that we have not explored, the damping effect can indeed be relevant and, within the specified approximations, quantifiable following the procedure we have outlined.
A last simplification is related to the fact that waves occur in a salty ocean. The densities of the fresh water due to rain and the salt water of the ocean are different. Both fresh water and air bubbles are, thus, affecting the composition of the ocean in the near surface and, thus, the local density and the tensor D(r). The change in the water composition will also affect some aspects of the bubble distribution model presented earlier. The more complex composition will quantitatively affect the results, but it is presumed not to affect the qualitative conclusions presented. Another effect that has been ignored here is a full consideration of the water–air interface stress conditions when rain is present. In Veron and Meussiens (2016), a model is presented for how the impact of rain affects the boundary conditions (see their Eqs. 2.11, 2.14, and, for waves, 3.4), i.e., the pressure at the free surface is composed of the atmospheric pressure plus a raininduced pressure. The raininduced pressure is trivially space averaged in order to enter the homogenized set of momentum/condition equations presented here. The consequences on the waves enter through Bernoulli's equation and affect the dispersion of the waves. These issues are presented in Veron and Meussiens (2016) and not repeated here.
There are several processes in which the mixture of a liquid containing a large number of bubbles or droplets has different characteristics. Examples are a frothy ocean, whatever the content, or cavitating flows, as it happens in ship propellers or in pressurized flows. An example we came across while preparing this article describes how foams are used to ameliorate unwanted ship motion due to sloshing of their holding tank contents (see Denkov et al., 2005; Kim et al., 2007) and the stabilizing effect of free surface sloshing of bubbly drinks (see Sauret et al., 2015). These foamy cases are nontrivial extensions of the homogenization procedure we present in this paper. Consideration of the chemistry, compressible effects, and topology is required, and the intricate formulation of the stress conditions at the interfaces would need to be derived. Nevertheless, the procedure may play a constructive role in formulating a mean field description of the dissipation and the consequent attenuation of the foam on sloshing motions.
We itemize our main conclusions here, as follows:

A general methodology has been developed to quantify the effect of smallscale bubbles on the properties of the containing liquid. Averaged dynamic quantities can be defined, and for these, averaged momentum and mass conservation equations of general applicability can be derived.

A model was derived that permits an estimate of the amount of air injected by rain.

The presence of tiny air bubbles, which we assume are uniformly distributed, changes the physical characteristics of the containing liquid, affecting the momentum balances – the momentum dissipation of irrotational motions is enhanced due to an increase in the rotational component of the flow.

We have also specifically addressed the impact of the presence of air bubbles on gravity waves. We found that air bubbles increase the effective wave damping and, hence, that rain has a damping effect on waves. Changes in the density near the free surface will also affect the surface tension. Since the generation of bubbles by rain is small in the natural setting, the damping effect is small, akin to the damping of waves due to surface contamination. An increase in wave damping has implications for momentum exchanges between the waves and the mean flow. Wave damping changes and surface tension changes impact the wave dispersion relation and the group velocity. These effects are primarily relevant to the dynamics of small and capillary waves.

We have cited a number of situations where, granted certain conditions, our approach can be applied for a general estimate of the overall increased dissipation.
All data are synthetic and generated by the expressions in the paper.
LC brought the unresolved problem of the mechanism, whereby rain calms rough seas, to the attention of the other authors. All authors developed the homogenization model for momentum exchanges near the ocean surface when impurities (in this case air) were present. All authors participated in the model for the connection between rain rates and the presence of resulting air bubbles in the upper ocean. All authors participated in equal measure in the preparation of the paper.
The authors declare that they have no conflict of interest.
The authors wish to thank the Kavli Institute for Theoretical Physics (KITP), at the University of California, Santa Barbara, for their hospitality and for supporting this research project. The KITP is supported, in part, by the National Science Foundation (NSF). Alex Ayet acknowledges the support from the DGA and the French Brittany Regional Council. Juan M. Restrepo acknowledges the support from the NSF/OCE grant. Alex Ayet would like to thank Bertrand Chapron for the insightful comments. We also acknowledge the comments and suggestions from the reviewers, which helped improve the clarity of the paper.
This research has been supported by the National Science Foundation Directorate for Geosciences (NSF/OCE; grant no. 1434198), the NSF (grant no. PHY1748958) and the DGA (grant no. D0456JE075).
This paper was edited by Joachim Peinke and reviewed by two anonymous referees.
Ayet, A., Chapron, B., Redelsperger, J.L., Lapeyre, G., and Marié, L.: On the Impact of Long WindWaves on NearSurface Turbulence and Momentum Fluxes, Bound.Lay. Meteorol., 174, 465–491, 2020. a
Babuška, I.: Homogenization Approach In Engineering, in: Computing Methods in Applied Sciences and Engineering, edited by: Glowinski, R. and Lions, J. L., Springer, Berlin, Heidelberg, Lecture Notes in Economics and Mathematical Systems, vol. 134, https://doi.org/10.1007/9783642859724_8, 1976. a
Bensousan, A., Lions, J. L., and Papanicolaou, G.: Asymptotic Analysis for Periodic Structures, NorthHolland Pub. Co., Amsterdam; New York, ISBN 9780821853245, 1978. a
Caflisch, R. E., Miksis, M. J., Papanicolaou, G., and Ting, L.: Wave propagation in bubbly liquids at finite volume fraction, J. Fluid Mech., 160, 1–14, 1985. a
Cavaleri, L., Baldock, T., Bertotti, L., Langodan, S., Olfateh, M., and Pezzutto, P.: What a Sudden Downpour Reveals About Wind Wave Generation, Procedia IUTAM, 26, 70–80, 2018. a
Cioranescu, D. and Donato, P.: An Introduction to Homogenization, vol. 17, Oxford Lecture Series in Mathematics and Its Applications, Oxford University Press, Oxford, 272 pp., 1999. a
Clift, R., Grace, J. R., and Weber, M. E.: Bubbles, drops, and particles, Courier Corporation, Dover Publications, Mineola, NY, 381 pp., ISBN13: 9780486445809, 1978. a, b
Denkov, N. D., Subramanian, V., Gurovich, D., and Lips, A.: Wall slip and viscous dissipation in sheared foams: Effect of surface mobility, Colloid. Surface. A, 263, 129–145, 2005. a
Dingle, N. and Lee, Y.: Terminal fallspeeds of raindrops, J. Appl. Meteorol., 11, 877–879, 1972. a
Harrison, E. L. and Veron, F.: Nearsurface turbulence and buoyancy induced by heavy rainfall, J. Fluid Mech., 830, 602–630, 2017. a
Henderson, D., Rajan, G., and Segur, H.: Dissipation of NarrowBanded Surface Water Waves, in: Hamiltonian Partial Differential Equations and Applications, edited by: Guyenne, P., Nicholls, D., and Sulem, C., Springer, Fields Institute Communications, 75, 163–183, 2015. a
Ho, D., Ascher, W., Bliven, L. F., Schlosser, P., and Gordan, E.: On mechanisms of raininduced airwater gas exchange, J. Geophys. Res., 105, 24045–24057, 2000. a, b, c
Jenkins, A. and Jacobs, S.: Wave damping by a thin layer of viscous fluid, Phys. Fluids, 15, 1256–1264, 1997. a
Keeling, R. F.: On the role of large bubbles in airsea gas exchange and supersaturation in the ocean, J. Mar. Res., 51, 237–271, 1993. a
Kim, Y., Nam, B. W., Kim, D. W., and Kim, Y. S.: Study on coupling effects of ship motion and sloshing, Ocean Eng., 34, 2176–2187, 2007. a
Kudryavtsev, V., Makin, V., and Chapron, B.: Coupled sea surfaceatmosphere model: 2. Spectrum of short wind waves, J. Geophys. Res., 104, 7625–7639, 1999. a
Lamb, H.: Hydrodynamics, 1st edn., Cambridge University Press, Cambridge, 1926. a
Leighton, T.: The acoustic bubble, Academic press, London, https://doi.org/10.1016/B9780124419209.500018, 1994. a
Le Méhauté, B. and Khangaonkar, T.: Dynamic interaction of intense rain with water waves, J. Phys. Oceanogr., 20, 1805–1812, 1990. a
Marshall, J. S. and Palmer, W. M.: : The distribution of raindrops with size, J. Meteorol., 5, 165–166, 1948. a
Medwin, H. and Clay, C. S.: Fundamentals of acoustical oceanography, Academic press, Boston, ISBN13: 9780124875708, 1997. a
Medwin, H., Kurgan, A., and Nystuen, J. A.: Impact and bubble sound from raindrops at normal and oblique incidence, J. Acoust. Soc. Am., 88, 413–418, 1990. a
Medwin, H., Nystuen, J. A., Jacobus, P. W., Ostwald, L. H., and Snyder, D. E.: The anatomy of underwater rain noise, J. Acoust. Soc. Am., 92, 1613–1623, 1992. a, b, c, d, e
Mendelson, H. D.: The prediction of bubble terminal velocities from wave theory, AIChE Journal, 13, 250–253, 1967. a
Merlivat, L. and Memerly, L.: Gas Exchange Across an AirWater Interface: Experimental Results and Modeling of Bubble Contribution to Transfer, J. Geophys. Res., 88, 707–724, 1983. a, b
Minnaert, M.: XVI. On musical airbubbles and the sounds of running water, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 16, 235–248, https://doi.org/10.1080/14786443309462277, 1933. a
Moghimi, S., Ramírez, J., Restrepo, J. M., and Venkataramani, S.: Mass Exchange Dynamics of Surface and Subsurface Oil in ShallowWater Transport, Ocean Model., 128, 1–12, 2018. a, b
Oguz, H. N. and Prosperetti, A.: Bubble entrainment by the impact of drops on liquid surfaces, J. Fluid Mech., 219, 143–179, 1990. a, b
Peirson, W. L., Beyá, J. F., Banner, M. L., Peral, J. S., and Azarmsa, S. A.: Raininduced attenuation of deepwater waves, J. Fluid Mech., 724, 5–35, 2013. a
Plesset, M. S. and Prosperetti, A.: Bubble dynamics and cavitation, Annu. Rev. Fluid Mech., 9, 145–185, https://doi.org/10.1146/annurev.fl.09.010177.001045, 1977. a
Prosperetti, A. and Oguz, H. N.: The impact of drops on liquid surfaces and the underwater noise of rain, Annu. Rev. Fluid Mech., 25, 577–602, 1993. a, b, c
Restrepo, J. M., Ramírez, J. M., and Venkataramani, S. C.: An Oil Fate Model for Shallow Waters, J. Mar. Sci. Eng., 3, 1504–1543, https://doi.org/10.3390/jmse3041504, 2015. a
Sauret, A., Boulogne, F., Cappello, J., Dressaire, E., and Stone, H. A.: Damping of liquid sloshing by foams: from everyday observations to liquid transport, Phys. Fluids, 27, 022103, https://doi.org/10.1063/1.4907048, 2015. a
Tsimplis, M. N.: The Effect of Rain in Calming the Sea, J. Phys. Oceanogr., 22, 404–412, 1992. a
Veron, F. and Meussiens, L.: A kinetic model for particle surface interaction applied to rain falling on water waves, J. Fluid Mech., 796, 767–787, 2016. a, b
Woolf, D. K. and Thorpe, S. A.: Bubbles and the airsea exchange of gases in nearsaturation conditions, J. Mar. Res., 49, 435–466, 1991. a
 Abstract
 Introduction
 The dynamical approach and the presence of bubbles
 Air entrainment due to rain
 The effect of rain on gravity and capillary waves
 Discussion and conclusions
 Summary
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 The dynamical approach and the presence of bubbles
 Air entrainment due to rain
 The effect of rain on gravity and capillary waves
 Discussion and conclusions
 Summary
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References