The concentration of microbes and other particulates is frequently enriched in the droplets produced by bursting bubbles. As a bubble rises to the ocean surface, particulates in the bulk liquid can be transported to the sea surface microlayer by attaching to the bubble’s interface. When the bubble eventually ruptures, a fraction of these particulates is often ejected into the surroundings in film droplets with a particulate concentration that is higher than in the liquid from which they formed. The precise mechanisms responsible for this enrichment are unclear, yet such enrichment at the ocean surface influences important exchange processes with the atmosphere. Here we provide evidence that drainage, coupled with scavenging, is responsible for the enrichment. By simultaneously recording the drainage and rupture effects with high-speed and standard photography, we directly measured the particulate concentrations in the thin film of a bubble cap at the moment before it ruptures. We observed that the enrichment factor strongly depends on the film thickness at rupture, and developed a physical model, based on scavenging and drainage, that is consistent with our observations. We have also demonstrated that this model is quantitatively consistent with prior observations of film drop enrichment, indicating its potential for a broader range of applications in the study of the sea surface microlayer and related phenomena.
Covering more than 70% of the Earth’s surface, the sea surface microlayer (SML) is critical to a range of transport processes involved in climate change and human health (Després et al., 2012; Cunliffe et al., 2013). For example, bacteria from the SML have been found in the upper atmosphere where they may serve as ice and cloud condensation nuclei, and toxins originating from red tide events are linked to respiratory irritation experienced by coastal residents (Woodcock 1948; Bauer et al., 2003; Kirkpatrick et al., 2004; Després et al., 2012; Cunliffe et al., 2013). In both of these cases, the microscopic particulates involved are transferred by the jet and film droplets formed from the rupture of bubbles entrained by breaking waves (Deane and Stokes, 2002; Lewis and Schwartz, 2004; Grythe et al., 2014; Veron, 2015). Furthermore, it has been observed that the particulate concentration in these droplets is often higher than the concentration in either the bulk water or SML, which amplifies the influence of these particulates (Blanchard and Syzdek, 1970; 1972; Angenent et al., 2005; Aller et al., 2005). Yet, for film drops, it is unclear how large an enrichment to expect; one previous study reports an enrichment factor between 10 and 20 (Blanchard and Syzdek, 1982), whereas another reports an enrichment factor of over 100 (Wangwongwatana et al., 1990).
Before reaching the ocean surface, entrained air bubbles scavenge particulates from the subsurface water during their rise, ultimately leading to the enrichment that defines the SML (Weber et al., 1983; Aller et al., 2005). As a bubble approaches the surface, it deforms it into a spherical cap (Toba, 1959). The shape of this cap remains constant until the bubble spontaneously ruptures, often fragmenting into film droplets. It is tacitly assumed that the concentration in the enriched film drops is the same as the concentration in the bubble film cap when it is freshly formed. However, because of natural surfactants, the bubble evolves after forming, draining and thinning until rupture (Lhuissier and Villermaux, 2012). Therefore, before rupture, many of the particles initially contained in the bubble cap likely drain back into the bulk liquid. This combination of draining and thinning makes it difficult to assess the impact of each factor on the final enrichment by solely collecting the film droplets produced, as has typically been done in the past (Blanchard and Syzdek, 1982; Wangwongwatana et al., 1990). Instead, our approach is to investigate the film contents immediately before rupture. Based on our findings, we have developed a physical model highlighting the roles of bubble scavenging and drainage in the film as it is rupturing and have demonstrated that our model is quantitatively consistent with the enrichment measurements from films drops collected in earlier studies.
To better understand the processes leading to an enriched bubble film, we performed a series of single bubble scavenging experiments. The experimental setup used to investigate this process is outlined in Figure 1. A cylindrical container was filled with a suspension of common yeast cells, Saccharomyces cerevisiae, prepared with tap water at various bulk concentrations Cb ranging from Cb ≈ 107 cells mL–1 to Cb ≈ 8 × 107 cells mL–1. We injected an air bubble of radius rb ≈ 2.5 mm (circled in Figure 1a) at the bottom of the yeast cell suspension. The newly formed bubble rose a distance H = 70 mm and established a spherical cap after reaching the free surface (identified with a dotted rectangle near the top of Figure 1a). The drainage and subsequent rupture processes were recorded simultaneously with high-speed and traditional cameras. To ensure that we were able to focus each camera on the bubble and its surface, we overfilled the plastic tube containing the suspension to form a convex meniscus at the top. This technique exploits the natural tendency of the bubble to “self center” in the tube under the influence of gravity.
Before the bubble ruptured, the film thinned and a fraction of the yeast cells initially contained in the bubble cap returned back to the bulk liquid by a combination of gravitational and surface tension forces. We captured this drainage, typically occurring over a time span of seconds, with a standard camera that was fitted with a microscope objective and recorded images at a rate of 30 frames per second. Figure 1b shows a segment of a bubble cap as it forms, rests on the interface for nearly two seconds, and then spontaneously pops (see Video S1). Note that the speckles in the images are the individual yeast cells. Additionally, these yeast cells serve as tracers for the flow, revealing that in this video the bubble is draining at approximately 0.5 mm s–1. These images were taken at a high enough resolution (1.3 microns per pixel) that the number of yeast cells per unit area, N, could be counted immediately before the film ruptured (Figure 1c). Because the film was curved, the bubble passed through the focal plane, resulting in a thin region in focus. Cell counts were performed only in this region.
The bubble rupture occurred when a small hole spontaneously initiated and rapidly expanded (see Figure 1d). The entire rupturing process, often occurring in less than t = 2 ms, was captured with a Photron FASTCAM SA5 high-speed camera at a rate of 60,000 frames per second.
Our choice in using yeast cells rather than other microbial particulates was deliberate. We needed a simple, safe system to make the next series of advances on this topic, though we recognize the potential future value of using marine microorganisms in seawater solution. Serratia marcescens, which has been used in the past (Blanchard and Syzdek, 1982), is classified as a human pathogen that requires specialized handling procedures, especially when aerosolized (Yu, 1979). Furthermore, the yeast cells had a measured particle size of rp = 3.2 µm (standard deviation of 0.4 µm), which is significantly larger than the previously used bacteria, allowing for easier visualization.
Results and Discussion
As might be expected, the number of yeast cells per area, N, counted in the region of interest at the time of rupture (Figure 1c) depended on the bulk concentration Cb, with higher bulk concentrations leading to higher number counts N (Figure 2a). In contrast, the film thickness at rupture h appeared to be independent of the bulk yeast concentration Cb over the range of concentrations tested (Figure 2b). Here we used the standard technique to measure the film thickness at rupture h by measuring the speed at which the nearly uniform film retracted (Pandit and Davidson, 1990; Lhuissier and Villermaux, 2012; Modini et al., 2013). The surface tension-driven retraction is regulated by inertia and is well approximated by the Taylor-Culick velocity, , where γ is the surface tension and ρ is the liquid density (Taylor, 1959; Culick, 1960). The thicknesses at rupture that we measured ranged from h = 5 µm to h = 60 µm. These values are larger than what might be expected in particulate-free bubble films (Spiel, 1998; Modini et al., 2013), and we suspect that the presence of particulates may act as nucleation sites that initiate rupture earlier in the drainage process. However, the physics that determine the thickness at which a water bubble spontaneous ruptures is currently unknown (Lhuissier and Villermaux, 2012).
By measuring both the particle surface concentration (N in cells cm–2) and film thickness h at rupture, we were able to calculate the volumetric cell concentration in the thin film Cf = N/h, which we report in cells per milliliter. Two noteworthy features become clear when the film concentration at rupture Cf is plotted as a function of the bulk concentration Cb (Figure 2c). First, there is no clear relationship between the two concentrations, and second, the film concentration Cf is categorically larger than the bulk concentration Cb. Specifically, the points in Figure 2c are all above the line denoting equivalence, indicating that the film became enriched.
The ratio of film concentration Cf to the bulk concentration Cb is defined as the enrichment factor EF at rupture. Further insight was gained by plotting this enrichment factor as a function of the rupture film thickness h (Figure 2d). For bubble films that were thick when they ruptured, the film concentration Cf is only slightly larger than the bulk concentrating Cb, resulting in an enrichment factor less than ≈ 5. However, if the bubble film was thinner when it ruptured, the enrichment factor tended to increase. For the thinnest films in our experiments, we measured enrichment factors near 20 (Figure 2d). When thin bubble caps of water drain, it is known that minuscule amounts of indigenous surfactants can rigidify the interface through Marangoni stresses due to surface tension gradients (Scriven and Sternling, 1960). This rigidity of the inner and outer surfaces may lead to the interior of the bubble draining at a faster rate. Indeed, we saw evidence of this phenomenon in our experiments (see Video S1) where, by taking advantage of our limited visual depth of field, we could see the interior of the film draining downwards while the inner and outer interfaces were, by comparison, stationary. In our experiments, these semi-rigid interfaces were likely caused by the natural surfactants introduced by the yeast.
We can estimate the thickness of the bubble film when it reaches the surface (Figure 1b) by recognizing that the drainage velocity and bubble thickness are coupled through conservation of mass. The film thickness over time decays exponentially if drainage velocity is constant. For the velocity, bubble size, and thickness at rupture, in Figure 1b, this model would predict a thickness of 75 µm at t = –1.9s before rupture. Note that there is evidence that the drainage may be better modeled with a power law (Lhuissier and Villermaux, 2012), although this model in comparison predicts the earlier thickness to be larger than a centimeter. Regardless of the precise drainage dynamics, the bubble thins as it drains until it ruptures at a thickness h that can be measured from the film retraction speed.
To develop a model for the enrichment factor in the spherical cap, we combined the measured values of film thickness at rupture h with the collision efficiency from scavenging theory Ecol, here defined as the fraction of particles contained in the volume swept out by the bubble that attach to its surface (See Figure 3 for model geometry). Here we use one of the simplest scavenging models, the Sutherland potential flow model Ecol = 3rp/rb, to estimate the number of particles brought to the surface by the rising bubble (Sutherland, 1948). For this model to be appropriate, the scavenged particles should be large enough to ignore diffusive effects (rp > 0.1 µm) and follow the streamlines around the bubble. The flow around the bubble should also be well-approximated as potential flow (Miettinen et al., 2010). Both of these conditions were met by our experiments. Additionally, we assume that any particle colliding with the rising bubble attaches (Figure 3a, b) and remains attached until bursting at the free surface. This collision mechanism is commonly referred to as interception, as the particle passes close enough to be “intercepted” by the rising bubble. Although there are more complex collision models that incorporate a large number of parameters, such as particle wettability, density mismatch, and the detaching of particles, our aim is to develop the simplest model that is able to capture the physics of the enrichment process. In addition, the flow conditions may also influence the collision model. For example, in highly turbulent environments, we anticipate the collision efficiency to require modification, typically for larger particles rp > 50 µm, to account for particle velocities that are uncorrelated to the flow conditions immediately surrounding the bubble (Abrahamson, 1975).
For simplicity, our model assumes that the scavenged particles are evenly spread across the bubble surface area. This condition leads to a uniform distribution of particles across the bubble cap area Af, an assumption that is consistent with our observations during data collection (see Figure 1b, c). The total number of particles in the bubble’s cap Nf is the sum of two contributions. The first contribution is from the particles scavenged during the bubble’s rise that attach to the interface Natt (Figure 3d) and the second contribution is from the particles contained in the film volume at rupture AfhCb. The number of particles scavenged by the bubble is determined by multiplying the chosen model for the collision efficiency by the number of particles encountered during the bubble’s rise . Dividing the total number of particles in the film (Nf = Natt + AfhCb) by the film volume at rupture Afh yields the film concentration Cf. Introducing a dimensionless grouping of parameters , that we have chosen to scale with the thickness, we finally normalize the film concentration Cf by the original bath concentration Cb, resulting in a prediction for the enrichment factor:
Re-plotting each of our individual yeast scavenging experiments from Figure 2 in terms of and the enrichment factor Cf/Cb results in a reasonable agreement with our model (Figure 4). However, like all models there are limitations, and we anticipate instances where our model may not yield good agreement. For example, the biotechnology industry routinely introduces additives, such as the non-ionic surfactant Pluronic® F-68, to prevent cells from attaching to bubbles, as rupturing bubbles are generally agreed to be the largest cause of damage to the cells being grown in production scale bioreactors (Hu et al., 2011). In this example, we would anticipate the film concentration to be approximately equal to the bulk concentration, as the primary method of initially trapping cells in the film has been inhibited, which effectively reduces the collision efficiency to zero. Alternatively, if the density of the scavenged particles were sufficiently large and the bubble long-lived, the particles may sediment, moving independently of the interface, and modify the predicted concentration.
Because we make the assumption that the yeast cells act as passive particles, we would anticipate similar results for inert, non-biological particles. To test this prediction, we ran a separate series of experiments using Polystyrene beads of rp = 3 µm. The same experimental procedure outlined in the Methods section was repeated. Additionally, the Polystyrene beads were suspended in the same water as the yeast after removing the yeast by centrifugation. This step ensured that any natural surfactants present in our original experiments would be maintained. The Polystyrene beads also followed the model reasonably well (Figure 4, blue squares). We note that our model tends to underestimate the observed enrichment factors measured in our setup. This underestimation may be due, in part, to the centering of the bubbles with a convex meniscus at the top of our cylindrical container. Specifically, the free surface was not continuously refreshed, leading to the concentration in the upper most layer of the suspension likely being higher than the bulk concentration (Blanchard and Syzdek, 1982).
We conclude our discussion by comparing our proposed model with two film droplet experiments reported in the literature. Because the thickness at rupture was not measured in either of these studies, we modeled the value of thickness to be equal to the particle size. The first study, by Blanchard and Syzdek (1982), injected air bubbles with radius rb = 0.85 mm that rose a height H = 19 mm through a suspension of Serratia marcescens bacteria. Because the bacteria are rod-like we use their equivalent radius rp = 0.65 µm (Weber et al., 1983) and assume the bubble ruptured at a thickness h = 1 µm, the bacterium’s larger dimension. Based on these estimations, we find reasonable agreement between the reported value of enrichment factor EF = 18 and our predicted value of EF = 12.
The second comparison study used latex spheres with particle radius rp = 0.28 µm (Wangwongwatana et al., 1990). Because of the particle size, this example approaches the applicability of our model, specifically that the particles were sufficiently large to ignore diffusive effects. Here the authors injected air bubbles of radius rb = 1.29 mm that rose a height of H = 600 mm to the free surface. Although the value of enrichment was not explicitly reported in the study, the authors state that the enrichment factor is well over 100×, compared to our expected value of 175.
Finally, we apply the same approximations to our data presented in Figure 4. Specifically, we compare the average enrichment factor over the course of our experiments, ignoring the variation between cases. Based on the crude assumption of the bubble rupturing at the thickness of the particle, h = 6.4 µm (the yeast diameter), and accounting for the rise distance in our experiments H = 70 mm, our model estimates an approximate overall enrichment factor of 11 versus our measured average EF value of 9. Taken together, these examples illustrate that our model provides a reasonable approximation for the expected enrichment for a wide range of scavenging conditions.
The findings presented in this paper suggest that a combination of both scavenging and drainage concentrate suspended particulates within a bubble cap prior to rupture of the bubble and thus provide a more detailed mechanism of film drop enrichment. Specifically, we provide a framework of the enrichment process as a combination of two steps. First, the bubble transports particles that attach to its interface during its rise to the surface (Weber et al., 1983). Second, after forming a spherical cap at the surface, the bubble drains and thins returning a portion of the scavenged particles back to the bulk. We find that the interior of the bubble cap film drains faster than the inner and outer surfaces. This difference in drainage increases the enrichment factor as the bubble thins. Beyond our bench-scale setup, we compared our model to previous studies and found quantitative agreement between the reported values and our values predicted by our model, supporting the notion that the processes contributing most to particle enrichment occur before the bubble ruptures. Previous studies have demonstrated that bubbles with thinner caps produce a larger number of film drops than thicker caps (Lhuissier and Villermaux, 2012). Additionally, the film drops that these thinner caps produce tend to be smaller and persist longer in the atmosphere (Lewis and Schwartz, 2004; Veron, 2015). Our results suggest that these large number of persistent droplets are also likely to be the most enriched. We anticipate our results may be important to processes involving enriched film drops such as may occur in the sea surface microlayer, with important implications to potential disease or toxin transfer nearshore and climate-related air-sea exchanges on a more global scale.
This work was supported by the National Science Foundation under Grant No. 1351466.
The authors have no competing interests to declare.
All authors contributed equally to this article.