Effects of background rotation on the dynamics of multiphase plumes

Abstract We report laboratory results on bubble plumes released from a point source into a homogeneous liquid environment of depth $H$ and rotation rate $\varOmega$. The gas phase is characterised by the non-dimensional slip velocity $u_N=u_s/(BH^{-1})^{1/3}$, where $u_s$ is the slip velocity of the dominant bubble species and $B$ is the source buoyancy flux. The effects of the background rotation are characterised by the Rossby number $Ro=(B\varOmega ^{-3})^{1/4}/H$. We study the regimes $0.06\lessapprox u_N\lessapprox 0.36$ and $0.03\lessapprox Ro\lessapprox 0.3$. We establish that, from $\varOmega t\approx 2.3$, the growth of the maximum plume width $b_{max}$ is slowed down compared to the non-rotating $t^{3/4}$ power law, where $t$ is the time. At $\varOmega t\approx {\rm \pi}$, the plume axis starts to tilt laterally and causes a slowdown of the rise of the height $h_c$ of the plume silhouette centroid from the $t^{3/4}$ power law. These critical times do not depend on $u_N$. After $\varOmega t \approx {\rm \pi}$, the slip velocity counteracts the effects of rotation such that $h_c$ rises faster for larger $u_N$. The subsequent onset of the anticyclonic plume precession causes the disintegration of the rising plume front into vertical columnar structures for $Ro\lessapprox 0.15$. Once the plume reaches the free surface, the subsurface lateral dispersion of bubbles is increased compared to the non-rotating case. However, background rotation suppresses the lateral dispersion of bubbles at the free surface. We find that, asymptotically, the surface area $A$ affected by the bubbles scales as $(B\varOmega ^{-1})^{1/2}t$, and the proportionality factor reduces with an increasing $u_N$.


Introduction
The Deepwater Horizon (DwH) blowout in 2010 resulted in the formation of a long-lived turbulent multiphase plume from the damaged rig site at the ocean bottom. This incident has attracted considerable public and scientific attention. The DwH plume consisted of water containing oil droplets and gas bubbles, and both dispersed phases contributed comparable amounts to the DwH source buoyancy flux with the oil and gas volume flow rate ratio around 1: 0.8 (Socolofsky, Adams & Sherwood 2011). The internal dynamics of the DwH plume exhibited multifarious physical, chemical and biological processes, including oil droplet formation, break-up and coalescence, temperature-and pressure-dependent gas dissolution and hydrate formation as well as biodegradation (Socolofsky et al. 2016). Additionally, the DwH plume dynamics and the subsequent dispersal of oil both on and below the ocean surface was influenced by external factors such as the ocean density stratification, local currents, waves, winds and the Earth's rotation. This wide range of mechanisms on different length and time scales presents a challenge for modelling the hydrocarbon transport in the environment and for development of effective mitigation strategies.
Motivated by this prominent example of a multiphase plume, we seek to understand the interplay between the background rotation and the additional forcing introduced by the presence of a multiphase effluent in a plume. In this study, we consider a bubble plume released into a rotating unstratified environment and investigate in detail the behaviour of its gas phase. In particular, we focus on the changes of the plume structure caused by the background rotation compared to a multiphase plume in the non-rotating case: we study the rise speed and the lateral spreading of the plume, which are directly related to the subsurface mixing and dispersion of the multiphase effluent and the quantity of gas and oil reaching the surface. Characterising these effects is important for assessing the risks associated with deep-water blowouts and for oil spill mitigation purposes. For example, the correlations for the subsurface location, the amount and the composition of the effluent can be used as the initial conditions for far-field transport models (Socolofsky et al. 2016).
Disregarding the wide range of chemical reactions, phase transitions and biochemical processes, the multiphase component can be characterised by the so-called slip velocity u s , which is the velocity of an individual bubble (or an oil droplet) relative to the surrounding fluid and can lead to the separation between the gas bubbles (or oil droplets) and the entrained water plume. Asaeda & Imberger (1993) and Socolofsky & Adams (2005) investigated in laboratory experiments the effects of the slip velocity u s on the bubble plume structure and the formation of subsurface lateral intrusions in quiescent non-rotating stratified environments. Depending on the relative magnitude of the slip velocity u s to the characteristic plume velocity (BN) 1/4 in a stratified environment, where B is the source buoyancy flux and N is the buoyancy frequency of the stratification, three plume types can be identified that differ in the pattern of the lateral detrainment and the structure of secondary plumes. In the DwH case, it is estimated that bubbles rise faster than oil droplets so that u s can be based on the effective gas droplet diameter (Socolofsky et al. 2011). The non-dimensionalised slip velocity u Ns = u s /(BN) 1/4 can therefore be placed in the range between 0.6 and 1.5 (Fabregat Tomàs et al. 2017;Socolofsky et al. 2011) which, following the classification by Socolofsky & Adams (2005), renders the DwH plume a Type 1* plume with a single distinct subsurface intrusion layer. Understanding the effects of u N on the plume structure is therefore crucial for modelling the vertical and horizontal transport of oil and gas after deep-water blowouts.
The duration of the DwH spill of several months implies that the DwH plume is likely to have been affected by the Earth's rotation with angular velocity Ω E . The conventionally defined Rossby number Ro s = N/f for the DwH plume is moderate Ro s ≈ 10, 915 A2-2 where f = 2Ω E sin θ is the vertical component of the Coriolis frequency at latitude θ. However, previous numerical and experimental studies showed that even at such moderate Rossby numbers, we can expect the background rotation to alter significantly the plume dynamics if the plume is maintained for a sufficiently long time (Fabregat Tomàs et al. 2016;Frank et al. 2017). In particular, regardless of the Rossby number, the system rotation modifies the plume dynamics after one rotation period due to the conservation of the angular momentum (Frank et al. 2017). The DwH oil spill occurred at the latitude of θ ≈ 28.7 • at which the Earth's rotation period is approximately T = 50 h ≈ 2 d. During the spill duration of 87 days, the DwH oil plume is hence very likely to have been significantly impacted by the Earth's rotation with associated substantial changes in the plume dynamics.
Laboratory studies of single-phase plumes in a homogeneous rotating environment were conducted by Fernando, Chen & Ayotte (1998) and Goodman et al. (2004). Fernando et al. (1998) proposed expressions for the critical times when the plume rise height and width start to be affected by the background rotation: Ωt ≈ 2.4 and Ωt ≈ 5.5, respectively, where Ω is the rotation rate of the environment and t is the time. Sheremet (2004) studied the scenario when there is a finite angle between the direction of gravity and the vector of rotation and documented novel features in the plume development that are absent in the case of no misalignment, such as, for example, the plume flow direction being intermediate between the gravity acceleration and rotation vectors, as well as the formation of tilted plume structures. Helfrich & Battisti (1991) investigated the development of a single-phase plume from a point source in a stratified rotating environment and established that a cyclonic circulation forms in the vicinity of the source and an anticyclonic eddy at the height of the lateral intrusion (that exists provided that the ambient fluid is deep enough, i.e. H 5(BN −3 ) 1/4 ). Furthermore, they observed that the thickness of the intrusion layer grows with decreasing Rossby number and the plume acquires a cylindrical shape. Convection from extended sources in a homogeneous rotating environment was examined by Maxworthy & Narimousa (1994). They noted that beyond a critical vertical distance from the source, z ≈ 12.7(B a /f 3 ) 1/2 with B a being the buoyancy flux per unit area, the advancing three-dimensional turbulent front disintegrates by dividing into a field of columnar vortices with a characteristic diameter ∼ (B a /f 3 ) 1/2 and vertical propagation speed ∼ (B a /f ) 1/2 . The vortex generation from buoyant line sources was investigated in detail by Bush & Woods (1999) for both stratified and homogeneous rotating environments. In particular, in a homogeneous case if the line plume is affected by the background rotation before it impinges on the free water surface, Bush & Woods (1999) reported the formation of unsteady anticyclonic columns of characteristic radius ∼ B 1/3 l /f with B l the source buoyancy flux per unit length. Possibly most relevant to the current study, Frank et al. (2017) observed and characterised the anticyclonic precession of the plume axis for a single-phase plume from a point source in a homogeneous rotating environment: the plume exhibits a pronounced lateral tilt after approximately one rotation period and the precession frequency ω depends linearly on the background rotation speed as ω ≈ (0.4 ± 0.04)Ω. The phenomenon of the anticyclonic plume precession was also observed in experiments by Ma, Flynn & Sutherland (2019) who confirmed the linear relationship between the precession frequency and the background rotation, although they reported a slightly higher proportionality factor ω ≈ (0.52 ± 0.09)Ω.
These previous studies examine the effects of rotation on the dynamics of single-phase plumes and provide useful background for the present investigation. However, at the time of writing, there appears to be no experimental studies that additionally consider the presence of a multiphase plume in a rotating environment. A recent numerical simulation of bubble plumes in a rotating environment was conducted by Fabregat Tomàs et al. (2017).
Using a simple Eulerian-Eulerian Boussinesq model for the monodisperse gas phase, where the gas phase is treated as a Eulerian scalar field obeying a convection-diffusion equation and is characterised by a slip velocity u s , they studied the behaviour of a bubble plume at Rossby numbers and slip velocities comparable to the estimates for the DwH oil plume. Such an approach does not model directly the bubble-induced turbulence, but Fabregat Tomàs et al. (2017) stated that the simple slip velocity formulation sufficiently captures the increase in turbulence production observed in the presence of bubbles. The simulations were conducted for a stratified rotating environment and represented direct numerical simulations of the plume at laboratory scale at a Reynolds number of 7000 based on the characteristic plume velocity and length scales, (BN) 1/4 and (BN −3 ) 1/4 , respectively. The main finding of these simulations was that the presence of the gas phase amplifies the rotational effects and also changes the plume structure and dispersion compared to a single-phase plume in a rotating environment. In particular, compared to single-phase plumes, the presence of the gas phase leads to a greater decrease of the detrainment height, a thickening of the intrusion layer, and it elevates the near-source concentration of a passive tracer injected at the source of the buoyant plume.
The goal of the present study is to investigate experimentally the dynamics of a multiphase turbulent plume from a point source in a homogeneous environment under the combined effects of the background rotation and the presence of a buoyant gas phase. We generated several bubble size distributions and characterised them by the slip velocity u s of the bubble species that contributed most to the buoyancy flux. For the non-rotating homogeneous environment, a characteristic plume velocity scale is given by (B/H) 1/3 (Morton, Taylor & Turner 1956), where B is the plume buoyancy flux and H is the water depth, so that a non-dimensionalised slip velocity can be defined as u N = u s /(B/H) 1/3 . Note that u N is related to the dimensionless parameter M H = B/(4πα 2 Hu 3 s ) used by Asaeda & Imberger (1993) Bombardelli et al. (2007), u N can be expressed as u N = (4πα 2 D/H) −1/3 . In the experiments reported in this paper, u N ranges between 0.06 and 0.36. We also vary the Rossby number Ro = (BΩ −3 ) 1/4 /H in the range 0.03 Ro 0.3 by changing the buoyancy flux B and the system rotation rate Ω. We note that, in this study, the plume buoyancy flux B is solely due to the presence of bubbles released at the plume source, and we focus only on characterising the gas phase and do not analyse the liquid phase of the plume.
The present work aims to fill the gap left by the scarcity of experimental results on multiphase plumes in a rotating environment. Although stratification was a very important factor for the dynamics of the DwH plume, we consider here the simplified problem of a multiphase plume released into a homogeneous rotating environment. The present investigation is intended as a first step in the study of the largely unexplored research area on multiphase plumes in a rotating environment. We intend to include stratification in future experiments. We focus here on three particular aspects of the problem. First, we consider the initial rise characteristics of the bubble plume, such as the rise height and the lateral spreading of the gas phase, and how they depend on the interplay between the slip velocity and the rotation of the environment. Second, we briefly discuss the subsurface bubble distribution after the plume has risen through the entire water column and has reached the free surface. Last, as an externally accessible diagnostic, we consider the signature of the bubble plume at the free surface and how it evolves in time.
The paper is structured as follows. We describe our experimental set-up and the experimental procedure in § 2.1. The bubble source conditions (the buoyancy flux B and the bubble size distributions) are characterised in § 2.2. Section 2.3 explains our procedure for the data analysis. Section 3 presents our experimental results with § 3.1 briefly focusing on bubble plumes in a non-rotating environment before considering features specific to rotating plumes. Sections 3.2, 3.3, 3.4 and 3.5 explore the anticyclonic plume precession, the initial rise characteristics, the late-stage evolution and the surface footprint of the plume, respectively. The findings of this paper are summarised in § 4, where we also briefly discuss the relevance of our results for geophysical plumes and deep-water oil plumes such as the 2010 DwH oil spill plume.

Experimental set-up
We conducted 144 experiments on bubble plumes discharged at the base of a homogeneous rotating saltwater environment. A schematic of the experimental set-up is shown in figure 1. The experimental tank, located on a rotating table and described in detail in Frank et al. (2017), had an octagonal base with a circumscribed circle diameter of 1 m and a maximum depth of 1.2 m. The rotation rate Ω used was typically varied in the range 0-1 rad s −1 , with an incremental step of 0.1 rad s −1 . Typically, the tank was filled to a depth of H = 0.6 m.
The bubbles were produced by means of electrolysis. The bubble source generator, manufactured in house, consisted of a tube, with a 15 mm inner diameter and a 50 mm length, which was connected to a circular plastic base (figure 1c). Six layers of electrodes, each comprising multiple 127 μm diameter platinum wires arranged in a semi-random manner, spanned the generator tube at different heights. These electrodes were connected to the power supply such that there were three layers of cathodes alternating with three layers of anodes so that both hydrogen and chlorine gas bubbles were produced. Small holes were cut through the tube wall next to where it connected to the plastic base to form a chimney that emits a two-phase flow of saltwater and bubbles. The bubble generator was hidden underneath a false bottom in the tank except for the top approximately 30 mm, which protruded through a 40 mm diameter hole in the centre of the false bottom.
The bubble generator was connected to a computer-controlled power supply. The number and the size of the bubbles produced depended on the current I flowing through the bubble source and the salt concentration (and hence the density ρ a ) of the ambient saltwater. For our experiments, we used two different saltwater densities ρ a ≈ 1008 kg m −3 and 1179 kg m −3 and four different currents I = 0.3, 0.5, 0.7, 0.9 A.  Table 1 summarises our experimental parameters. Note that the gas production rate Q used to calculate the source buoyancy flux B was measured by collecting the gas bubbles in a measuring cylinder that was initially filled with water. This is explained in detail in § 2.2 and appendix A. The flow was recorded using two JAI Spark SP-5000C-CXP2 5-megapixel cameras that were fixed in the rotating frame of reference: one camera was attached to the sidewall of the tank and the other recorded the top view (figure 1a,b). We used the side illumination by means of two red LED panels for the side-view recordings of the plume. For the top-view recordings, a 2 mm thick horizontal laser sheet created by a 32 mW 520 nm Flexpoint laser, fitted with a 45 • Powell lens, was positioned approximately 15 mm below the water surface. The centrifugal force, caused by the rotation of the tank, produced a curvature of the water surface, especially for high values of Ω, so that the water level was depressed at the tank centre. As a consequence, the laser sheet was closer to the water surface in the centre of the tank than at the tank edges. Unfortunately, no simultaneous side-view and top-view recordings could be made because the light from the side LED panels would have interfered with the laser sheet and polluted the top-view videos. Instead, each experiment was repeated for the same experimental conditions twice: once recording with the side-view camera, and once with the top-view camera. The side-view camera was equipped with a wide angle 6 mm c-mount lens that allowed visualisation of the entire vertical extent of the plume from the source to the free surface.
The experimental procedure was as follows. The tank was spun up from rest for at least eight spin-up times T s = H(νΩ) −1/2 /4 (Greenspan & Howard 1963), where ν is the kinematic viscosity of water, until it reached solid body rotation. The video recording was started and, immediately, the bubble source was switched on with the selected current. The plume was captured for 300 s at 10 frames per second for the side-view and at 5 frames per second for the top-view videos. The exposure time was set to 100 and 200 ms, respectively. We note that for the source buoyancy fluxes B used, the rising times of the plume through the water column were of the order of magnitude of 20 s. Furthermore, the fixed recording time of 300 s corresponded to approximately 5 rotation periods for Ω = 0.1 rad s −1 and 50 rotation periods for Ω = 1 rad s −1 . In practice, we used only the first 120 s of our recordings for the data analysis. At the end of an experiment, the bubble source was switched off and a new rotation rate for the tank was selected. After an appropriate amount of time that was sufficient for the remaining bubbles to disappear and for the water in the tank to equilibrate in a new solid body rotation state, a minimum of four spin-up times T s , a new experiment could be started. We changed the water in the tank approximately every 10 experiments to avoid any yellow colouring of the water in the tank caused by potential secondary chemical reactions of other impurities. No evidence of a significant gas dissolution in the plume was observed during an experimental run.

2.2.
Measuring the source buoyancy flux and the bubble size distribution In order to characterise the bubble plumes, we needed to determine both the source buoyancy flux B (which is governed by the current I through the bubble generator) and the size distribution of the bubbles produced in order to estimate the bubble slip velocity u s . The methods used to measure B and u s are described in detail in appendices A and B, respectively.
The gas production rate Q was found not to depend significantly on the salinity of the ambient water and to scale approximately linearly with the source current I. The source buoyancy flux B is inferred from the gas production rate Q as where g = 9.81 m s −2 is the acceleration due to gravity. Table 1 includes the approximate source buoyancy fluxes B that were achieved for each experiment. The bubble size distribution of the bubble generator was measured using a procedure described in detail in appendix B. The probability distribution of the bubble radii r is reasonably well approximated by a log-normal distribution f r (r), i.e. ln f r (r) ∼ N (μ, σ ), where N(μ, σ ) is a normal distribution with the mean value μ and the standard deviation σ . Table 2 lists the values of the associated parameters μ and σ for the fitted log-normal distributions. The radius corresponding to the peak of each distribution, denoted by r m , is calculated as exp(μ − σ 2 ).
The buoyancy distribution function f b (r) is calculated from the bubble radius distribution function f r (r) by multiplying the latter by r 3 and normalising it (2. 2) The bubble radii corresponding to the peaks of the buoyancy distribution functions are labelled as r mb and their values are recorded in table 2. We characterise each bubble size distribution by the slip velocity u s based on the bubble radius r mb containing the largest contribution to the buoyancy flux. The slip velocity u s is calculated from the balance of the buoyancy force and the friction force acting on a bubble (Stokes' law) 4 3 πr 3 mb gρ a = 6πμ dyn r mb u s , where μ dyn is the dynamic viscosity of saltwater. Table 2 includes the values calculated for u s and also for the non-dimensionalised bubble slip velocity u N . Note that the Reynolds number Re b based on the bubble radius of 50 mm, the slip velocity of 0.5 cm s −1 and the kinematic viscosity of water of 10 −6 m 2 s −1 can be estimated as Re b ≈ 0.25 which means that the flow around the bubble is laminar and (2.3) applies.

Image analysis
We analyse the experimental data by means of a formalised automated procedure for the post-processing of image sequences. In the following, we present this algorithm and, in particular, discuss our definition of the plume extent.
The recorded side-view raw image sequence shows the plume as a bright structure against a black background. The bubbles reflect the light from the side LED panels into the camera: the more bubbles that are present at a certain location (integrated in the direction normal to the camera field of view), the greater the light intensity and the brighter this location appears in an image. Appendix C shows an example of the experimental images of the plume at each stage of the post-processing. The camera recorded the raw images in the 8-bit system. We note that in the 8-bit system where 'black' is 0 and 'white' is 255, a typical mean light intensity of the raw background image in our experiments was around 9.5. A typical temporal variation in the light intensity of a background pixel (that was not affected by the plume at any time) was at most 1 or 2 in the 8-bit system. The first image of a recorded sequence, when the plume is not yet started, is used to remove the background from the subsequent frames through a simple transformation: the photographic negative (i.e. the image, normalised by the saturated intensity of 255, then subtracted from one) of every frame in the sequence is divided by the photographic negative of the first framê Here, I is the raw intensity of the image, I 0 is the intensity of the first frame andÎ is the intensity of the transformed image. This yields a sequence of images where the plume now appears as a dark structure against a white background. Upon a close inspection of these new transformed images, we found that the white background possesses a very uniform intensityÎ. The light intensityÎ of most background pixels assumes values in the range of 0.996−1 where 'white' is considered as 1 and 'black' as 0. Thus, we convert these transformed images into binary images of the plume where every pixel with the light intensityÎ < 0.995 is considered to contain bubbles. We now formalise our definition of the plume extent: the plume is defined to be the largest connected object in a binary image -effectively, the two-dimensional projection or silhouette of the plume. At each time step, the search for the largest object with pixels connected by any one of their eight nearest neighbours is performed using MATLAB ® 2016a.
Our choice of the plume as the largest connected object ensures that individual large bubbles that escape due to their slip velocity are not considered to be a part of the plume. In other words, a region is within the plume if the number of bubbles in this region is high enough and there is a gapless path between the imaged pixel of this region and the plume source. We note that the light intensity in our images is integrated in the direction normal to the camera field of view, so 'a region' is a two-dimensional projection of a three-dimensional object. In practice, this may lead to some uncertainty around the boundary where an edge region may or may not be considered to be a part of the plume depending on the orientation of the projection. However, this higher-order distinction would require either information of both gas and liquid-phase velocity field or a complete three-dimensional mapping of the plume structure.
Our definition of the plume extent is based on the distribution of the bubbles alone, since the light intensity distribution is linked to the bubble distribution, rather than on the liquid-phase velocity field. It accounts only for the presence of the bubble phase and not for the liquid phase since no information on the latter can be extracted from our experimental data. Under field conditions in real oceans, our definition of the plume extent would correspond to using a probe that measures the chemical composition of water and determines the pollutant concentration (for example, oil or gas). Then we define a location x to be within the plume if the pollutant concentration (accumulated along a horizontal line with a chosen direction through x) surpasses a certain threshold and there is a gapless path within the plume from x to the plume source. If the concentration accumulated along a horizontal line through x is too low -too little light is scattered into the camera in our experiments -then x is not regarded as being part of the plume. This definition is sensible because first, we have a direct access only to the light intensity field and not to the velocity field. Second, the radial profile of the vertical velocity for a plume in a rotating environment is no longer expected to follow a well-defined Gaussian distribution at each height, but exhibits a rather complicated structure with the flow reversal around the plume edges (see further discussion in § 3.3.3 and Fernando et al. 1998). Furthermore, we do not expect any azimuthal bias based on the view point of the camera, owing to the axisymmetric geometry of the set-up and the turbulent nature of the plume flow. Thus, our plume definition shows representative behaviours, especially when experiments are regarded in their ensemble.
The binary images are subsequently used for determining the edges of the plume in MATLAB ® 2016a. The detected edges are stored in a three-dimensional array with two spatial and one time dimension. In each pixel row of the image, the leftmost pixel where an edge has been detected is considered to be the left bounding edge of the plume and the rightmost edge pixel is defined as the right bounding edge of the plume. The region between the left and the right bounding edges is now defined as the extent of the plume and at each height, the width of the plume b is the distance between the left and the right bounding edge.
A similar procedure is used in the post-processing of the top-view images.

Bubble plume in a non-rotating environment
Before investigating the combined effects of the bubble slip velocity and the ambient rotation rate on the plume dynamics, we briefly discuss the role of the bubble slip velocity for the plume in a non-rotating environment. McDougall (1978) modelled a bubble plume as consisting of an inner core in which the majority of the bubbles are concentrated and which is surrounded by an annular flow of the entrained fluid. Milgram (1983) described the bubble plume as a single plume of bubbles and water with different Gaussian profiles for the velocity and density defect which were based on the plume radius and the smaller gas-containing radius, respectively. We reiterate that our imaging technique only captures the gas phase. Thus, hereafter 'the plume' only refers to the gas phase. Figure 2 shows average images of the plume for B = 50 cm 4 s −3 for two values of the non-dimensionalised slip velocity u N = 0.07 and 0.27. The averaging was performed over 60 s (600 frames) after the plume has reached the free water surface. The conical shape of the plume (bubble phase) appears visually thinner for a larger slip velocity u N (which corresponds to a smaller ambient density ρ a , see table 2). We formalise this by considering the vertical region between 10d 0 and 25d 0 above the nozzle (indicated by vertical double arrows in figure 2), where d 0 denotes the inner diameter of the bubble source. In each pixel row, a Gaussian is fitted to the light intensity distribution using a least-squares procedure. The standard deviations σ b of the fitted Gaussians as a function of height z above the nozzle are displayed in figure 3 for different u N . We observe that σ b ∼ z so that a linear fit σ b = α b z can be made: the coefficient α b as a function of u N is shown as the inset in figure 3 with the red dashed line showing the value for the entrainment constant α = 0.12 of a single-phase plume (see e.g. Morton et al. 1956;Papanicolaou & List 1988). This confirms our qualitative observation that the plume width reduces with increasing u N . We note that the coefficient α b cannot be regarded as the entrainment constant of the plume since the entrainment is defined for the liquid phase which we do not visualise.
The next feature we investigate is the effect of the bubble slip velocity on the initial rise height of the plume. For a single-phase plume from a point source with a source buoyancy flux B in a homogeneous non-rotating environment, the rise height h of the plume front is expected to scale as where t is the time after the start of the plume (Turner 1962).  The rise height h of the plume at each time point is calculated as the distance from the source to the lowest height near the top of the plume where the width of the plume first vanishes (based on the criterion described in § 2.3). Figure 4 shows a log-log plot of h/H as a function of the non-dimensionalised time (B/H 4 ) 1/3 t. We set t = 0 as the time when the plume starts. With the gas phase present, the height of the plume grows slightly faster than the scaling t 3/4 for homogeneous plumes. We fit the rise height h with a power law of the form at β up to the height of 25d 0 and find that the power coefficient β varies between approximately ∼t 0.75 and ∼t 0.87 . Based on the entire set of data presented in this paper and by considering the repeatability of experiments (see also § 3.3.1 and, in particular, figure 6), the error associated with the power coefficients is estimated to be ±0.03. Thus, an increase in the power-law coefficient with an increasing u N is statistically significant. Here, we stress that we define the plume rise height as the vertical extent of the largest connected structure in the bubble plume, as determined by our criterion described in § 2.3. For the series E, F, G and H, we observed that some individual bubbles separated from the plume front and rose faster than the plume. Had the height of the plume been based on the rise height of these separated bubbles, then the scaling would follow approximately h/H ∼ t, a speed comparable with that of the individual bubbles in isolation.
We can attribute the value of β being larger than 3/4 to the bubble slip velocity and our definition of the plume. Specifically, we consider only the bubble phase of the plume and define a region to be within the plume if the concentration of bubbles in that region is high enough. It would be interesting to see what happens to the liquid phase upon starting the bubble plume and what power law the rising plume front follows if it is defined by means of a passive (dye) tracer injected at the source. Predicting the value of β and how it depends on u N would be valuable in order to understand better the impact of bubbles on the plume dynamics. Likewise, investigating further the dependence of α b on u N would be useful for understanding the dynamics of these multiphase plumes. A careful study of these problems, however, would require a larger data set for a wider range of slip velocities than we currently have. This would also require measurements of the velocity field of the 0.10 0.05 0.20 0.60 1.00 liquid phase of these plumes. The data presented in this section are primarily included as a base non-rotating case (Ω = 0 rad s −1 ) to allow for comparison between the rotating and the non-rotating cases and for determining the effects of background rotation discussed in the following sections.

Precession frequency of the bubble plumes in a rotating environment
The bubble plumes exhibited essentially the same lateral deflection and precession phenomenon that was reported in Frank et al. (2017) for single-phase plumes: after approximately one rotation period, the plume was deflected laterally and then started to precess in the anticyclonic direction. When viewed from above, the tank was rotating in the counter-clockwise direction, so the direction of the observed precession was clockwise relative to the tank. We use a similar procedure to that of Frank et al. (2017), based on the Hilbert-Huang transform, to extract the precession rate of the plume axis. Figure 5 shows the extracted mean precession rateω of the bubble plume as a function of Ω. We observe thatω is a linear function of Ω with a proportionality factor 0.41 ± 0.05, which is a similar value to single-phase plumes (ω ≈ (0.4 ± 0.04)Ω, Frank et al. (2017) and ω ≈ (0.52 ± 0.09)Ω, Ma et al. 2019). The error estimate for the proportionality factor is based on the 95 % confidence interval for the coefficient estimate of the linear regression model. The independence ofω from B was already shown experimentally in the previous study (Frank et al. 2017). Here, we confirm that the precession rate is also independent of the bubble slip velocity u N . This is to be expected from the theoretical spinning-top toy model developed by Frank et al. (2017). The theoretical estimate (see (5) in that study) for the precession frequencyω relies on the value for the entrainment coefficient as one of the main factors. Bubble plumes considered in this paper may be less turbulent at the source and lazier, and possess an entrainment coefficient smaller than 0.1. However, changing the value of the entrainment coefficient to 0.08 instead of 0.1 used in Frank et al. (2017) results only in a marginally different proportionality factor of 0.43 instead of 0.35, that would be expected from the theoretical estimate. Such a difference is unlikely to be observed Rotation rate Ω in rad s -1 Figure 5. Extracted precession rateω of bubble plumes as a function of the rotation rate Ω for different source buoyancy fluxes B and slip velocities u N . The precession rateω scales linearly with Ω, with a proportionality factor of 0.41 ± 0.05, similar to single-phase plumes, and appears to be independent of u N . experimentally in the proportionality factor for the precession frequency, especially when taking into consideration the scatter of the data presented in figure 5.

Initial rise of the bubble plume in a rotating environment
In this section we investigate the combined effects of the slip velocity u N and the rotation rate Ω on the initial rise characteristics of the bubble plume, before the plume reaches the free water surface. For deep-water blowouts, or naturally occurring oceanic plumes in general, the only information readily available are the ocean surface observations. Therefore, here we first investigate the time delay between the initiating of the plume at the bottom and the appearance of its signature on the water surface. We also discuss the distribution of the plume fluid inside the water column by considering the rise height of the plume silhouette centroid and the disintegration of the first plume front.

Rise height of the plume front
The rise height h of the plume, determined as described in § 3.1, is plotted in figure 6 for six experimental series (see table 1) with t = 0 and h = 0 defined as before. We non-dimensionalise h as h/(BΩ −3 ) 1/4 and the time t as Ωt, similarly to the scalings introduced by Fernando et al. (1998) which coincides with (3.1) for the rise height of a single-phase plume in a non-rotating environment. Fernando et al. (1998) found that the rise height of the plume slows down from the 3/4 power law for Ωt 2.4. Close inspection of the rise height of bubble plumes presented in figure 6 does not reveal any observable change in the ascent of the plume front at any time within an experiment. Moreover, for a given source buoyancy flux B and a non-dimensional slip velocity u N , the curves for all Ω collapse within the margin of experimental error and appear to follow a The discrepancy between the slowing down of the single-phase plume initial rise beyond a critical time reported by Fernando et al. (1998) and the apparent absence of such a deviation in our data for bubble plumes -for any slip velocity u N -might be explained by the following remarks. First, Fernando et al. (1998) used the planar laser-induced fluorescence (PLIF) technique to visualise the plume. A sheet of laser light creates a two-dimensional plane through the plume and may obscure the three-dimensional nature of a plume in a rotating environment: the maximal rise height may be offset from the

Plume fingers
Plume fingers Figure 7. Snapshots of the plume at t = 20 s after its start for series A. The top row, showing the unaltered plume images, may misleadingly give the impression that the rise height of the plume is significantly lowered at high Ω because gas bubbles are increasingly concentrated around the source. However, the corresponding binary images, displayed in the bottom row, confirm that when the plume vertical extent is defined by a more sensitive concentration threshold, then its rise height does not change with Ω beyond a measuring error: (a) vertical line above the source due to the plume precession. Also, the disintegration of the first plume front, which will be discussed in § 3.3.4, might additionally contribute to this offset. Second, Fernando et al. (1998) define the extent of the plume to be the contour line of the 5 % of the maximum fluorescent dye concentration whereas we use a finer threshold of 1 % change in the (transformed) light intensityÎ. The rotation of the environment changes the distribution of the plume fluid within the water column, increasing the bubble concentration around the source. Hence, it is possible that the contour lines for different concentrations change in various ways: the 5 % contour line may be slowing down while the 1 % contour line does not. Last, the Rossby number based on the bubble radius and the bubble slip velocity is large (u s r mb −1 /Ω = O(10 2 )), meaning that the rise of an individual bubble (relative to the fluid containing it) is unaffected by the rotation of the environment. Since the rise of the plume front is governed by the rise of the largest (relatively isolated) bubbles, this implies that it should not be affected by the rotation, as shown by the results presented in figure 6. Figure 7 presents different plume snapshots taken at a fixed time after their starts for different rotation rates Ω. We use here dimensional times in order to compare the plume rise height to the non-rotating case Ω = 0 rad s −1 for which the non-dimensional time Ωt cannot be meaningfully defined. If in the rotating case the plume rises as (Ωt) 3/4 (see (3.2)), then Ω cancels and we are left with h ∼ B 1/4 t 3/4 , which is the same as in the non-rotating case. So, in that case we would expect to see the same vertical rise height of the plume for any Ω and the same t. The unaltered plume images in the top row show clearly the increased accumulation of the gas bubbles around the source for higher rotation rates. This may create a visual impression that the plume rises slower for higher rotation rates. However, the binary images in the bottom row demonstrate that the plume possesses the same vertical extent within the margin of accuracy, although the concentration of the gas bubbles at the plume front is greatly reduced with increasing Ω.
This result may have the following implications. If a sensor detecting the effluent concentration (e.g. gas bubbles or oil drops) was placed at the water surface, then the time delay between the start of the plume at the (ocean) bottom and the detection time of the plume by the probe would not differ between rotating and non-rotating environments. However, the initial effluent concentrations reaching the water surface in the rotating case are expected to be lower than in the non-rotating case, see figure 7.

Rise height of the plume silhouette centroid
The plume images presented in figure 7 demonstrate that, although the rise velocity of the plume front does not appear to be affected by the rotation rate Ω, the overall shape of the plume changes considerably. Here, we investigate the behaviour of other characteristic properties of the plume shape during its initial rise: the centroid of the plume silhouette and the plume width.
As a characteristic property in the vertical direction, we consider the temporal evolution of the silhouette centroid height h c of the plume. For each frame, the centroid height is calculated as follows. First, the plume shape is identified by extracting the largest connected structure from the binary images (see § 3.1). The centroid height is then calculated as is a two-dimensional vector and the index p runs over all pixels that are located inside the plume envelope P identified from the binary images. Here,Î is the transformed intensity of the pixel and z( p) is the vertical distance of the pixel from the bottom of the image. By definition, the centroid height h c of the plume silhouette calculated by (3.3) is a quantity which is purely based on the light intensity of the recorded images. It allows the tracking of the visually perceived changes to the plume shape. As this definition of h c does not depend on the image intensity being proportional to the bubble concentration, the temporal evolution of h c can be meaningfully discussed for any functional relation between the bubble concentration and the recorded light intensity. However, the centroid height h c can be linked to the physical height of the centre of mass (deficit) of the bubble plume. We have not measured the exact functional relation between the scattered (and recorded) light intensity and the bubble concentration although it can be surmised to be linear in the first approximation. Instead, we can calculate theoretically the height of the centre of mass (deficit) of the self-similar plume and compare it to h c . For the self-similar plume for which the buoyancy flux is constant with height, the mass deficit is the same at every height. Thus, for a self-similar plume, the height of the centre of mass (deficit) is located at h/2 where h is the rise height of the plume. We plotted h c as a function of time for our non-rotating experiments (not shown here) and confirmed that h c follows closely h/2. Figure 8 shows the mean non-dimensionalised centroid height h c /(BΩ −3 ) 1/4 as a function of the non-dimensionalised time Ωt for the six experimental series in a log-log plot. The averaging is performed over ten experiments in each series. Figure 25 in appendix D displays the data for each experiment separately. Initially, for each series, the data collapse onto a power-law curve (see also figure 25). However, h c subsequently deviates from this power-law growth and rises slower for Ωt 3. A least-squares power-law fit of the form for Ωt > 5. The exponent β 1 of (3.4) is β 1 ≈ 0.75 as is expected for plumes unaffected by the background rotation, see (3.2). The non-dimensionalised deviation time Φ c from the initial power law is determined by calculating the intersection point between the two power-law fits (3.4) and (3.5). The critical deviation time is plotted in figure 9(a) as a function of u N , while Φ c appears to be independent of the slip velocity u N . Furthermore, the deviation from the initial power law occurs for Ωt π and we estimate Φ c = 3.5 ± 0.5. Frank et al. (2017) reported that the lateral deviation of the plume can be clearly observed for π < Ωt < 2π. This suggests that the main reason for the slowdown of the rate of increase in time of h c is the onset of the anticyclonic plume precession. Section 3.2 showed that u N does not affect the plume precession, so our measurement of a constant Φ c is in agreement with that fact.
To characterise the growth of h c for Ωt > Φ c we fit the data with (3.5). The power exponent β 2 of these curves is shown in figure 9(b) as a function of u N . The plot suggests that β 2 rises with increasing u N , which indicates that the bubble slip velocity counteracts the slowdown of the h c growth induced by the lateral tilting of the plume and the onsetting precession. When the plume is initially tilted, the bubbles can rise vertically through the slanted upper side boundary of the plume envelope due to their slip velocity, thus reducing the inclination angle ψ from the vertical of the buoyancy defining the plume and, as a result, opposing the slowdown of the h c growth. We conjecture that ψ → 0 as u N → ∞ so that for very large u N the change in the power exponent from β 1 to β 2 may be undetectable. Additionally, bubbles with a large u N are less likely to be carried downwards by a downflow that develops around the plume edges (see § 3.3.3), which again would oppose the reduction in the rate of increase in time of h c .
As was already discussed in § 3.3.1, Fernando et al. (1998) observed the deviation of the plume rise height for single-phase plumes from the initial power law for Ωt ≈ 2.4, which is less than the measured Φ c in our experiments. We link Φ c to the onset of the plume precession whereas Fernando et al. (1998) did not report any such phenomena.

Growth of the maximum plume width
To investigate the effects of the background rotation and the slip velocity on the plume in the horizontal direction, we consider the maximum plume width b max . We extract the plume width for each time frame and each height from the binary images (see § 3.1). The data are relatively noisy because the data for the maximum plume width are instantaneous measurements. Averaging the data for each series over ten rotation rates (shown in figure 10) shows that initially the data collapse onto a power-law curve with an exponent very close to 3/4. This finding agrees with our expectation that the plumes are not affected by rotation at initial times. Subsequently, we observe a deviation from the 3/4 power law with the maximum plume width growing more slowly.  Figure 12. Images of the plume for series C and Ro ≈ 0.12 (Ω = 0.3 rad s −1 ) for four consecutive times. The plume still possesses a conical shape for Ωt = 1.8 with the maximum plume width located close to the plume front. The plume gradually starts to change its shape around Ωt ≈ 2.4. For Ωt ≈ 3.6, the plume shape becomes more cylindrical with a uniform width: the dashed lines around the plume shape are drawn to guide the eye. At the same time, we observe the onset of the lateral deflection of the plume, which is indicated with the red arrow: (a) Ωt ≈ 1.8; (b) Ωt ≈ 2.4; (c) Ωt ≈ 3; and (d) Ωt ≈ 3.6.
is determined in a similar way to h c by making a power-law fit to the average data for Ωt < 2 and a second fit, to the data in the range Ωt 4. Figure 11 presents the critical deviation time Φ b for different u N . We note that Φ b is not affected by u N and find Φ b ≈ 2.3 ± 0.8. Note that the critical time Φ b < π, which means that the plume width growth changes before the onset of the lateral plume deflection and the subsequent plume precession. Figure 12 illustrates the images of the plume (series C, Ro ≈ 0.12, Ω = 0.3 rad s −1 ) for four non-dimensionalised times. It confirms that the plume shape starts to undergo visible changes, from the conical shape at Ωt = 1.8 to the more cylindrical shape at Ωt = 3.6, at around Ωt ≈ 2.4 as suggested by Fernando et al. (1998). A close examination of the plume images for several series and rotation rates reveals that around the non-dimensionalised time Φ b , a downflow at the plume edges starts to develop (see also movie 1 in the supplementary material available at https://doi.org/10.1017/jfm. 2020.1181). This observation is in agreement with the study by Fernando et al. (1998) who explored saltwater plumes discharged vertically downwards into a homogeneous rotating environment and reported the beginning upflow around the boundaries of the plume around Ωt ≈ 2.4. They explained this phenomenon by suggesting that the rotation inhibited the horizontal entrainment flow. The plume then entrains fluid vertically which leads to a flow reversal around the plume edges. The same physical reason applies for the observed downflow around the edges of the bubble plumes: here, however, the bubbles with a large u N can escape the vertical downdrafts. Fernando et al. (1998) observed a slowdown in the plume front propagation for Ωt > 2.4. In our case, as was discussed in § 3.3.1, the plume ascent of the topmost connected edge continues unabated. It is currently not clear whether this disparity is a result of differences in visualisation approaches or whether this is the consequence of the bubble slip velocity. Further, Fernando et al. (1998) did not notice any significant change in the plume width growth until Ωt ≈ 5.5. After that time, their experimental results suggest that the maximum plume width remains constant. In our experiments, as is illustrated in figure 10, the growth of b max continues at a slower rate for Ωt 2.3. However, b max does not attain a constant maximum value but persistently increases with time. In fact, the starting downflow around the plume edge at Ωt 2.3 as well as the constant supply of the source buoyancy flux, provides a first hint that the vertical buoyancy flux may no longer be the same at every plume height: we conjecture that the buoyancy is accumulated around the source for Ωt 2.3, an idea that will be discussed further in the next subsection. This accumulation, a process which cannot persist for a long time for a continuously injected plume, may also be one of the contributing factors to the onset of the lateral plume deflection which may be observed for Ωt π and the ensuing plume precession. We conclude this section by examining the maximum plume width b max at one particular time instant t H when the plume reaches the free water surface, see figure 13. This time instant t H corresponds to different non-dimensional times Ωt for different experimental series because of different B and u N . We note that for a single-phase plume in a non-rotating environment, the maximum plume width at t H will always be constant around 2αH (Morton et al. 1956;Papanicolaou & List 1988), where α = 0.12 is the entrainment coefficient. Owing to self-similarity, it will not depend on how long it takes the plume to rise through the entire water column. We also note that for a fixed series, t H (as a dimensional time) will not vary significantly for different Ω because the plume front rise does not appear to be affected by the rotation (see § 3.3.1). The data in figure 13 are noisy because they depict instantaneous values of b max . Nevertheless, they allow us to recognise a general trend that the maximum plume width reduces for smaller Ro = (BΩ −3 ) 1/4 /H and larger u N . We note that the data are plotted with the log-log axes to improve the visual separation of individual data points for smaller Ro. The reduction with a decreasing Ro can be linked to the rotation inhibiting the horizontal entrainment discussed above. The reduction with u N is consistent with the data in figure 3 that show that the bubble phase of the plume becomes narrower for an increasing u N . We note that in the rotating environment, the maximum plume width is not necessarily found close to the free water surface, as will become obvious in the next subsection. Understanding the maximum plume width when the plume reaches the free water surface is, therefore, important for assessing the maximum subsurface horizontal extent of the polluted region. Maxworthy & Narimousa (1994) studied the turbulent convection of saltwater into a homogeneous rotating fluid of lesser density from horizontally extended sources (i.e. possessing a radius R so that R/H ≈ 1 where H is the water depth). Their experimental results suggest that the advancing turbulent three-dimensional front becomes rotationally dominated after a critical depth and disintegrates by dividing into a field of columnar vortices. These structures possess a characteristic radius and a constant vertical propagating velocity which scale with the surface buoyancy flux (defined as the total buoyancy flux per unit area) and the rotation rate of the system. Similarly, Bush & Woods (1999) reported that the front of a rising line plume in a homogeneous rotating environment undergoes an instability generating Taylor column structures whose horizontal dimensions again uniquely depend on the source buoyancy flux (per unit length) and the rotation rate of the environment. We observe a similar behaviour for point-source plumes. Under the conditions outlined below, the rising plume front loses its bell-shaped profile and disintegrates into a collection of columnar plume fingers, which can be seen in figures 7(c) and 7(d). The formation of the plume fingers appears to be at least partially related to the deflection of the plume and its subsequent precession. The mechanism for the onset of the formation of plume fingers is illustrated in figure 14. Initially, the plume is unaffected by the background rotation and discharges vertically upwards, figure 14(a). After Ωt ≈ π, the plume axis starts to deflect by tilting laterally, figure 14(b). However, the tilting does not occur uniformly along the entire plume height, the plume rather changes its discharge direction by flowing through the sides of its former shape envelope. This process leaves behind a patch of plume fluid, designated as parcel 1, which continues to rise vertically along the plume axis. The plume continues to change its direction and thereby sheds parcels of fluid in different directions, as can be observed in figures 14(c) and 14(d). Since the reason for the change in the discharge direction is the plume precession rather than the plume moving with the ambient fluid, parcel 2 and all subsequent parcels are not significantly swirled around the axis with the precessing plume. The parcels of plume fluid that have been shed rise vertically and become stretched into columnar vortices by the rotation of the environment: parcel 1 is transformed into a plume finger in figure 14(d). We note that the fluid parcels mentioned here do not appear to be surrounded by bubble-free fluid but remain topologically connected to the plume, and the plume source. However, they are also no longer directly supplied by the source fluid because of the changed plume discharge direction. In the simplest situation, the plume would be expected to form a helical structure due to the anticyclonic precession. In reality, the process of the plume precession is not completely steady (possibly due to interaction with the fingers), meaning that the plume can discharge for a longer period in some directions, which leads to the formation of the fluid parcels visible in figure 14. Furthermore, the occurrence of a downflow around the plume edges (see § 3.3.3 and Fernando et al. 1998) may additionally disrupt the formation of a spiralling plume shape. Note that the plume fingers do not appear to have a constant width as they rise: instead they are stretched and become thinner ( figure 15). This observed decrease in finger width appears to be different from the vortex columns forming in experiments with extended sources. The latter possess a characteristic radius as was reported by Maxworthy & Narimousa (1994). Also, it is not immediately clear how fast the plume fingers rise. We have not measured the upwards propagating velocity explicitly but our measurements of the first plume front rise height, which can be regarded as the rise height of the longest plume finger for small Ro, suggest that this velocity decreases in time. This again would be different from columnar structures observed by Maxworthy & Narimousa (1994) which have a constant propagating velocity. The noted discrepancies may originate in the mechanism of the plume finger formation that can be attributed to the plume precession rather than to the wave-like instability of an advancing plume front: the individual plume fingers are not continuously fed by buoyancy flux whereas columnar structures reported by Maxworthy & Narimousa (1994) are. We also note that the shedding of fluid parcels that are eventually stretched to plume fingers occurs within one or two cycles of the plume precession after its onset.

Disintegration of the first plume front
We can obtain an estimate for the range of Ro for which we can observe the formation of the plume fingers by noting that a pronounced lateral tilting of the plume is present at Ωt ≈ 2π (Frank et al. 2017). Thus, the first plume finger will start to form if the rise height of the plume is less than the water depth after that time leading to c(Bt 3 ) 1/4 = c(2 3 π 3 BΩ −3 ) 1/4 H. (3.8) The proportionality factor of approximately c ≈ 1.8 was reported in several previous studies (see Turner 1962;Fernando et al. 1998;Sheremet 2004). This is equivalent to Ro 0.14. (3.9) Another way to provide an estimate can be based upon arguing that plume fingers need to form if the expected Rossby radius of deformation is less than the width of the plume when it reaches the free water surface. In the threshold situation when the plume is still unaffected by the background rotation during its rise, its velocity scales as (B/H) 1/3 and the expected Rossby radius of deformation is therefore which can be transformed to Ro 0.1 3/4 ≈ 0.18.
(3.11) Sheremet (2004) performed experiments on single-phase plumes in a rotating environment when the gravity vector and the rotation axis are misaligned. In their experiment C2 (Ro ≈ 0.042 based on the parameter data provided), they noted that the incoming plume disintegrates into turbulent billows and diffuse tilted (because of the misalignment between the gravity and the rotation axis) columns, which eventually lead to the formation of tilted 'ink walls'. Sheremet (2004) does not discuss the mechanism for the formation of columnar structures nor their propagation characteristics. However, they provide an estimate for the distance beyond which the columnar structures appear to originate, by considering when the background rotation suppresses the three-dimensional turbulence and using experimental data by Fernando et al. (1998). For the case z c = H this translates to Ro 0.15. (3.13) All three threshold estimates (3.9), (3.11) and (3.13) are very close and, indeed, in our experiments with small u N we can observe the onset of the plume finger formation around Ro ≈ 0.15. This is shown in figure 16. For smaller Ro, the plume front clearly disintegrates into columnar structures (series A and series C in figure 16). An increasing slip velocity appears to counteract the formation of distinctive plume fingers and the plume acquires a more pronounced helical shape (series H in figure 16). We also note that snapshots for series H in figure 16 bear a resemblance to the images showing volume renderings of the gas fraction presented by Fabregat Tomàs et al. (2017). However, a direct comparison is difficult since that numerical study was conducted for a stratified environment. We also note here that for small slip velocities, the effects of the bubble separation from the entrained water plume are negligible and thus, the gas phase of the multiphase plume may be regarded as an approximation for a single-phase plume. Consequently, also for single-phase plumes we expect to see a similar disintegration of the rising plume front into plume fingers due to the mechanism of the plume precession.

Late-stage evolution of the plume
After the plume has reached the free surface, gas bubbles (or oil droplets) continue to accumulate throughout the height of the water column and spread laterally. This late-stage evolution of the plume is largely unexplored so that the pollutant distributions inside the water column, their concentration and the horizontal extent of the affected region, are mostly unknown.
In this section, we present the results using dimensional units rather than in a systematic re-scaled form. At this point we have been unable to find any simple underlying scaling laws. Consequently, the results discussed here are mostly of qualitative nature. Figure 17 presents four snapshots of the plume at t = 70 s for different rotation rates and slip velocities. We observe that the plume structure is highly complex and asymmetrical. For low Rossby numbers (figures 17b and 17d), sharply bounded vertically elongated structures are formed. These structures are similar to the 'ink walls' reported by Sheremet (2004). It is currently unknown how the temporal evolution of these structures will depend on the source parameters, the rotation rate and the bubble slip velocity. We expect here a complex flow pattern with a competition between different mechanisms such as an upflow due to buoyancy, downflow due to the downward entrainment, bubble slip velocity changing the buoyancy and leading to an accumulation of smaller bubble sizes inside the water column as well as formation of Stewartson and Ekman layers. The discussion of this late-stage evolution of the plume presented here is by no means exhaustive and warrants further in-depth studies. However, it relates to an interesting question for deep ocean blowouts: how does the horizontal extent of the contaminated underwater region evolve in time? Figure 18 presents the plots of the plume width b at each height z through the water column as a function of time (see colour bars) for different series and rotation rates. The plume width b and the vertical distance z have been normalised by the water depth H. The time t is experimental time (in seconds) with t = 0 defined here as the instance when the plume reaches the free water surface. Note that we use a very simple non-dimensionalisation of the spatial variable by H; this choice is arbitrary and not motivated by any scaling law. At every height, the plume width b has been determined by subtracting the location of the left bounding edge from the right bounding edge of the plume (see § 2.3). As is expected, in the non-rotating cases (left column of figure 18), the plume width b retains a temporally constant conical shape throughout the water column. There is a visible outflow region next to the free surface where the plume width b grows in time for z/H ≥ 0.8. In the rotating cases, we can distinguish two regions depending on the distance from the source. In the vicinity of the source, the width of the contaminated region b increases with the distance from the nozzle. The growth of the plume width b with height z/H in this near-source field is larger than the non-rotating case. In this region the plume precession occurs, which enhances the lateral dispersion of the bubbles. Further downstream, for a fixed instant in time t, the horizontal spread of the plume appears to remain approximately constant with height z/H. This is the region in which sharply bounded 'ink walls' are formed. Frank et al. (2017) reported that the vertical height of the precession region scales as (3.14) In figure 18, black dashed lines indicate the expected transition height using the scaling (3.14) between the near-field spreading region and the region of constant width. For small slip velocities, series A, C and D, (3.14) provides a reasonably good estimate for the transition height. For larger slip velocities, series E, G and H, the vertical extent of the near-field spreading region appears to be larger than given by (3.14), which is in line with the fact that the slip velocity enhances the vertical transport of the bubbles. In particular, for Ω = 0.1 rad s −1 , the scaling (3.14) yields values of 0.57H, 0.  The plume width plots presented in figure 18 suggest that in the constant-width region, for a fixed Ω, the plume width slightly increases with the buoyancy flux; compare plots for the series A, C and D (B = 30 cm 4 s −3 , B = 70 cm 4 s −3 , B = 90 cm 4 s −3 , respectively) as well as series E, G and H (B = 30 cm 4 s −3 , B = 70 cm 4 s −3 , B = 90 cm 4 s −3 , respectively). Similarly, for a fixed source buoyancy flux, the plume width b appears to decrease with an increasing Ω. However, it was not possible to determine any conclusive scaling laws based on our experimental data. Moreover, we find the bubble slip velocity does not have any significant effect on the temporal evolution of the plume shape at this late stage after the plume has reached the water surface.
In summary, this behaviour leads to the following pattern in the underwater bubble distribution in the rotating case. In the near-field region, up to a distance of approximately 2.6(BΩ −3 ) 1/4 above the source, the lateral dispersion of the bubbles is increased by the plume precession compared to the non-rotating case. However, the lateral extent of the bubbles in this near-field region remains approximately constant with time. Note, however, that our experiments were conducted for a finite time, so we cannot make any definitive statements about the long-time asymptotics. Further downstream, in the constant-width region, the initial lateral extent of the bubbles may be lower than in the non-rotating case. However, as time progresses, this region thickens and, eventually, the lateral extent of the bubbles will again be large compared to the non-rotating case. Finally, close to the water surface, the absence of a distinct outflow region for Ro 0.22 suppresses the lateral dispersion of bubbles for all times. Indeed, we expect that the affected area on the water surface will be smaller than in the non-rotating environment, which will be discussed in the next section.

Affected area on the water surface
For oceanic plumes, the only immediate and easily accessible information are the observations of the effluent spreading on the ocean surface. In this section, we discuss the temporal evolution of the surface area A affected by the multiphase effluent. This area is calculated from the processed binary images of the top-view recordings (see § 2.3) as where x = (x, y) is a two-dimensional position vector in the plane of the surface and S is the entire area of a top-view image. The index p runs over all pixels in the image.
Here, I b is the transformed intensity in a binary image, taking the value I b = 0 for black pixels (affected by bubbles) and I b = 1 for white pixels (not affected by bubbles). The intensity I b of binary images is calculated from the intensityÎ of transformed images in an analogous manner to the processing introduced in § 2.3. The diagnostic A as defined in (3.15) is a direct measure of the surface area affected by the presence of bubbles. Fernando et al. (1998) studied the lateral spreading of a single-phase plume from a point source in a rotating environment after its impingement on a solid boundary. They considered the case when the plume starts to feel the presence of the background rotation only after it reaches the boundary which, based on their experimental results, implies Ro = (BΩ −3 ) 1/4 /H > 0.3. In that regime, lateral spreading continues until it reaches a critical radius r c ≈ 1.75(BH −1 ) 1/3 /Ω that is proportional to the Rossby deformation radius. Subsequently, the laterally spreading front becomes baroclinically unstable and breaks down into a number of anticyclonic eddies. Figure 19. Snapshots from the top-view recordings showing the bubbles spreading on the free surface for series DT for Ro = 0.29, 0.086, 0.051 (left to right) and t = 1, 5, 9 s (top to bottom). The time is counted from the moment that we detect the impingement of the plume on the free water surface. We observe that for low Ro, there is no clear initial lateral spreading phase of the plume and that the affected surface area is very patchy even immediately after the plume reaches the free surface: Most of the experiments in the present study are in the range Ro < 0.3. Figure 19 shows snapshots from the top-view recordings of the plume for series DT (table 1) and for three Rossby numbers Ro = 0.29, 0.086, 0.051 (Ω = 0.1, 0.5, and 1 rad s −1 , respectively). The elapsed time t is measured from the instant when the plume first reaches the free water surface. In practice, we find the frame corresponding to t = 0 by manually inspecting our top-view recordings for each experiment and selecting the first frame in which we can observe the presence of the plume at the water surface.
Qualitatively, for experiments at very low Ro, we could no longer observe a clear lateral spreading phase after the plume impinges on the free water surface, as is shown in the second and third columns of figure 19. This is in line with the observations from § 3.4 shown in figure 18 that an outflow region can be detected for Ro ≈ 0.29 but not for Ro 0.22. In figure 19, for Ro 0.1, the affected surface area A appears to be incoherent and patchy, which can be attributed to the disintegration of the plume rising front into columnar plume fingers as was discussed in § 3.3.4.
For the case when the plume is affected by the background rotation only after it impinges on the free water surface, Ro 0.3 (Fernando et al. 1998), the velocity scale of the plume at the surface is u ∼ (B/H) 1/3 , which also corresponds to the frontal propagation velocity of the outflow (Fernando et al. 1998). Thus, immediately after the impingement, the radius of the deflected flow (i.e. the outflow) is expected to grow as ∼ ut and, hence, For the range of Ro when the plume is significantly affected by rotation during its initial rise, the region of constant plume width extends through the water column for z 2.6(BΩ −3 ) 1/4 (see § 3.4). Thus, by scaling arguments, we non-dimensionalise the surface area as A/(BΩ −3 ) 1/2 (note that for bubble plumes B and Q are directly related, see (2.1)) and the time as Ωt. Once the region of constant plume width is formed, the buoyancy influx into this plume column is the source buoyancy flux B and is constant in time. If we assume that the reduced gravity in that column is uniform, i.e. the bubble concentration is uniform, then its volume needs to scale linearly in t. Since the vertical extent of the constant-width region does not appear to vary in time (see § 3.4), the cross-sectional area A (which is the same as the surface area, but we note that it does not need to be circular, see figure 19) of the column needs to grow linearly in time. Thus, we conjecture that (3.17) Figure 20 shows the measured affected surface area A for series CT, DT, FT, GT and HT for different rotation rates Ω. The colour legend for the plotted lines is consistent throughout the figure and is shown in the right panels. The left panels of the figure display the plots of A as a function of time. We observe that, generally, rotation suppresses the affected surface area A. For the non-rotating case (black line), the initial affected surface growth appears to be consistent with the expectation A ∼ t 2 , which is fitted as a red dotted line for 0 s < t < 10 s. In experiments at Ω = 0.1 rad s −1 the lines initially follow the affected area growth of the non-rotating case (Ro = 0.27 for series CT and GT, Ro = 0.29 for series DT and HT, Ro = 0.25 for series FT). This is consistent with our observation that there is an initial outflow region close to the water surface for this range of Ro, see figure 18. As time progresses, the curves for Ω = 0.1 rad s −1 deviate from the non-rotating case. For higher rotation rates Ro 0.17, the measured affected area curves almost immediately deviate from the non-rotating case.
The data presented in figure 20 reveal that the growth of A in the rotating case is subject to characteristic oscillations around a rising trend, especially for lower rotating rates. When the time axis is non-dimensionalised as Ωt (middle panels), we observe that the first slowdown in the temporal evolution of A occurs for 2 Ωt 4. The onset of the lateral plume tilting (and the subsequent plume precession) occurs for Ωt ≈ π, so the first oscillation may be the (time-delayed) surface signature of the plume precession onset.
As time progresses, A exhibits a linear growth in the rotating case. This is demonstrated in the right panels of figure 20 in which we plot the ratio of A/(BΩ −3 ) 1/2 and Ωt as a function of the progressing time t. Asymptotically, we recognise that A(BΩ −1 ) −1/2 t −1 tends to a constant value (which appears to vary with u N ), which confirms our prediction (3.17).
We also observe the trend that an increasing slip velocity u N decreases A: this can be, for example, noted in the right panels of figure 20 where the asymptotic value is lowered for larger u N (series GT and HT). Also in the non-rotating case (left panel), the slip velocities presented here appear to stop the growth of the surface area because bubbles escape out of the water rather than spread laterally.

Conclusions and discussion
The aim of this study was to investigate the combined effects of the background rotation and the bubble slip velocity u N on the dynamics of a multiphase plume. We conducted experiments on bubble plumes, produced by electrolysis, that were discharged into a homogeneous rotating saltwater environment. A consistent image analysis of the recorded image sequences by means of formalised algorithms was performed in order to study the initial rise characteristics of the plume, its late-stage behaviour and its behaviour at the free surface.
Summarising our experimental results, the following timeline for the initial bubble plume development can be established. During the initial rise of the plume, we observe a downflow around the plume edges starting at around Ωt ≈ 2.3. At this point, a change in the temporal evolution of the maximum plume width occurs. The growth of the maximum plume width is slowed down in comparison to the non-rotating t 3/4 power law, albeit the width continues to increase without an apparent upper limit. The background rotation and the downflow initiated around the plume edges presumably cause a change in the plume entrainment, inducing a lateral deflection of the plume, tilting it away from the vertical at Ωt ≈ π. This deflection subsequently develops into an organised anticyclonic plume precession. The lateral plume tilt causes a slowdown of the rise of the plume silhouette centroid which deviates from the t 3/4 power law for Ωt π. The bubble slip velocity counteracts the effects of rotation such that, after the tilting occurs, the centroid rises faster for larger u N . The rise height of the plume front, however, does not appear to be affected by the background rotation at all and coincides with the rise height in a non-rotating environment. Its temporal evolution seems to be governed by bubbles which rise faster for larger u N .
This timeline for a multiphase plume is different than for a single-phase plume, as reported by Fernando et al. (1998). They also suggested the occurrence of the downflow after Ωt ≈ 2.4. However, their observations suggested that beyond that time the evolution of the plume height decelerated below the t 3/4 power law found for a non-rotating environment. The difference in the timelines between our experiments and the experiments of Fernando et al. (1998) may arise from the fact that individual bubbles around the plume front possess a very large Rossby number based on their size and are thus not directly affected by the background rotation. Additionally, the flow visualisation by means of PLIF may have obscured the three-dimensional precessing plume structure with the maximum rise height being offset from the vertical line above the source. Subsequently, Fernando et al. (1998) observed a change in the plume width growth only around Ωt ≈ 5.5 and the maximum plume width was constant after that time. However, they did not account for the anticyclonic plume precession (Frank et al. 2017) in their data analysis and did not consider the changes in the centroid rise.
During its initial rise through the water column, the front of the bubble plume may lose its coherent structure and disintegrate by dividing into a collection of columnar fingers. This observation is similar to the reported disintegration of advancing turbulent fronts from extended sources (Maxworthy & Narimousa 1994) and line plumes (Bush & Woods 1999). However, for the point sources used in this study, we notice that the mechanism for the formation of the plume fingers is related to the onset of the anticyclonic plume precession. This allows us to provide an estimate Ro 0.15 which defines the regime where we can observe the formation of the plume fingers inside a water column. The formation of plume fingers may imply that the buoyancy flux is no longer conserved with height, leading to the accumulation of buoyancy around the source. This accumulation is also related to the fact that the early-stage evolution of the plume is unsteady and not self-similar. Since the anticyclonic plume precession is not modified by the presence of bubbles, we expect the same mechanism to also work for single-phase plumes, especially since our bubble plumes for very small slip velocities can be regarded as an approximation for the single-phase plumes with the same buoyancy flux. We notice that an increasing u N may counteract the formation of plume fingers by organising the bubble phase into a spiralling structure.
Once the bubble plume has risen through the entire water column and reached the free surface, we investigated the temporal evolution of the effluent spread in the interior of the water column. The majority of the experiments reported in this study are in the range Ro < 0.3 with only a few at Ro ≈ 0.3. Generally, at these rotation rates, two regions in the plume structure can be recognised. The near-source or the precession region, where the plume width expands with the vertical distance from the source z and the region where the plume acquires a cylindrical shape and its width remains constant with z (the Taylor columns or 'ink walls' reported by Sheremet (2004)). The estimate for the vertical extent of the precession region 2.6(BΩ −3 ) 1/4 (Frank et al. 2017) provides a reasonable approximation for the transition height between these two regions for small u N but underestimates the transition height for larger u N (based on our experiments for u N 0.2). Additionally, for bubble plumes at Ro ≈ 0.3 there is a third region present close to the water surface where the plume starts to spread laterally upon reaching the surface. This outflow may inhibit the clear formation of the constant-width region. As time progresses, the plume width in the near-source region does not vary significantly. In contrast, further downstream, the width of the affected region grows smoothly in time uniformly across the whole vertical extent of the 'ink wall' region. This growth appears to be faster for an increasing B and slower for an increasing Ω, but no simple underlying scaling laws could be established based on the available set of the experimental data.
Summarising, after the plume reaches the water surface, the following changes in the subsurface pollutants spreading compared to the non-rotating environment can be observed. In the near-source region, the lateral dispersion of pollutants is spatially increased due to the plume precession (see also Frank et al. (2017)). For the heights z such that 2.6(BΩ −3 ) 1/4 z 0.8H, the initial dispersion of pollutants may be first lower than in the non-rotating environment but the lateral width of the contaminated region grows in time and will eventually surpass the non-rotating case, see figure 18 for z/H 0.8. Finally, close to the water surface, for z/H 0.8, the absence of an outflow region implies that the lateral dispersion of pollutants for Ro 0.3 will be reduced compared to the non-rotating environment, see figure 18 for z/H 0.8 and figure 20.
On the water surface, we can observe the signature of the changed subsurface effluent dispersion. For Ro 0.3, there is an initial spreading phase of the plume which follows closely the affected area growth in the non-rotating environment before deviating and slowing down. For Ro 0.3 (we note here that there is no sharp transition), the affected area growth is immediately smaller than in the non-rotating environment. There are characteristic oscillations around a rising trend in the affected area growth: the first dip occurs between 2 < Ωt < 4 which leads us to conclude that this may be the surface signature of the onsetting plume precession (the lateral tilting of the plume occurs for Ωt ≈ π). Asymptotically, the affected area A scales as (BΩ −1 ) 1/2 t. An increasing slip velocity u N additionally suppresses the growth of the affected area A due to bubbles escaping out of water. Fernando et al. (1998) considered only the case Ro 0.3. Therefore, we can conclude that an increasing slip velocity u N of the multiphase effluent counteracts the effects of the background rotation in the vertical direction but amplifies them in the horizontal direction.
We now return to one of the main motivations for this work and briefly discuss the relevance of our results to the DwH spill. In our laboratory study, we neglected several factors present in real oceans and could explore only a certain parameter range. First, we considered only a homogeneous rotating environment with no stratification. So, the results presented here might be relevant for plumes released into homogenised oceanic regions that for example are associated with the regions of deep oceanic convection (Steffen & D'Asaro 2002). The present study is intended as a first step in the largely unexplored research area on multiphase plumes discharged into a rotating environment. We intend to include the effect of stratification in a future study which will allow us to make better comparisons to the DwH oil plume.
The parameter range we explored, Ro < 0.3, is slightly outside of the Rossby number for the DwH plume. For the DwH plume, it can be estimated Ro s ≈ 10 or Ro = (BΩ −3 ) 1/4 /H ≈ 1.2, where we use the parameters B = 0.48 m 4 s −3 , Ω = 2π/86 400 s × sin(28.7) • ≈ 4 × 10 −5 rad s −1 , H = 1500 m (Socolofsky et al. 2011). However, several of the results we report here -such as the onset of the plume precession, the increased subsurface lateral dispersion in the rotating environment and the asymptotic behaviour of affected surface area -are independent of the Rossby number Ro = (BΩ −3 ) 1/4 /H defined in terms of the plume initial rise. Thus, we expect these findings to also apply for plumes at higher Rossby numbers such as the DwH plume. The Earth's rotation period at the latitude of the DwH oil spill is approximately 2 days, so the plume precession period (based on our results for the homogeneous rotating environment, see also Frank et al. 2017) can be expected to be approximately 4 days. However, we recognise that ambient conditions such as inertial oscillations and background currents are likely to have an impact on the precession signal. The explored parameter range Ro < 0.3 should also be applicable to natural gas or oil seeps with B ≈ 10 −3 − 10 −2 m 4 s −3 . Using these lower values for B (as well as Ω = 2π/86 400 s ≈ 10 −4 rad s −1 , H = 1500 m), we obtain Ro = (BΩ −3 ) 1/4 /H ≈ 0.1 − 0.26. However, we realise that a coherent bubble plume may not develop from a natural seep and the bubbles may continue to rise individually through the water column with their slip velocity (Roemer et al. 2012). Experimentally accessing a regime of a higher Rossby number is challenging. A significant reduction in the rotation rate Ω (even if achievable in a controlled way) would require a larger tank so that the vertical boundaries do not disturb the outflow on the water surface. A reduction in the water depth H would at some point interfere with the finite dimensions of the plume source and increase the influence of viscosity.
We note that in our homogeneous experiments the non-dimensional slip velocity u N = u s /(BH −1 ) 1/3 ranges from 0.06 to 0.36. For the DwH plume, released into a stratified environment, the parameter u Ns = u s /(BN) 1/4 is estimated to be between 0.6 and 1.5. Using the DwH values for the slip velocity of the gas bubbles u s ≈ 0.21 m s −1 (based of the effective gas bubble diameter of 2 cm), the source buoyancy flux B = 0.48 m 4 s −3 and the trap height h t = 300 m (Socolofsky et al. 2011), we can estimate u N = u s /(Bh −1 t ) 1/3 ≈ 1.8. However, during the DwH oil spill large amounts of dispersants were applied that altered the oil and potentially the bubble size distributions and led to smaller slip velocities. A reduction of the bubble size just by a factor of 3, would lead to a reduction of the non-dimensional slip velocity by a factor of approximately 10, which is inside the range that we examined. This highlights the need to study the behaviour of bubble plumes over a wide range of non-dimensional slip velocities.
It was not possible to achieve larger values of u N in our experiments. Producing larger bubbles through electrolysis required further reducing the ambient water density ρ a that in turn implied using higher voltages to maintain the same currents (and the source buoyancy fluxes). However, it was not possible to reliably start a bubble plume in our experiments if a voltage of more than 15 V was needed. Likewise, significantly reducing B below 30 cm 4 s −3 did not yield a steady bubble plume since the power supply that was available could not maintain a steady value for small currents I. Finally, the water depth of H = 60 cm was the maximum that could be used while allowing the entire plume to be within the field of view of the camera.
We also note that our electrolysis technique for generating bubbles produced a log-normal distribution rather than the bubbles being monodisperse. This in fact mirrors the real conditions of the DwH plume, although it hampers the modelling efforts. The actual droplet size distribution generated by DwH is a major focus of ongoing research, but in the absence of a definitive model is often predicted to be log-normally distributed (Socolofsky et al. 2016). One solution to this problem would be to swap from bubbles to particles with a very narrow size distribution (e.g. Mingotti & Woods 2019). However, such experiments will introduce slightly different physics at the boundary where the particles 'settle' (in principle, particles that were less dense could be used) and would require a larger volume flux through the nozzle to maintain a well-mixed suspension prior to release.
Lastly, we note that in this study we chose to focus only on the behaviour of the gas phase of the bubble plume and did not consider the liquid phase. For oceanic multiphase plumes, the primary question is to understand the dynamics of the multiphase component since it is the oil and gas that cause the most damaging effects for the environment and not the entrained seawater plume. We plan to investigate the liquid phase of a multiphase plume in a rotating environment in a future study.

Pump
Power supply Figure 21. Sketch of the set-up for measurements of the bubble size distribution.
The tank was filled with the homogeneous saltwater of different densities ρ a that were used in the rotating bubble plume experiments up to the height of about 55 cm. The bubble generator was placed at the bottom of the tank. Bubble plumes for four different currents I = 0.3, 0.5, 0.7, 0.9 A (approximately corresponding to source buoyancy fluxes B ≈ 30, 50, 70, 90 cm 4 s −3 , respectively) were discharged into the tank. A specifically constructed acrylic (Perspex) tube with a rectangular cross-section of 4 cm × 2 cm and a length of 50 cm was fixed to a metallic frame around the tank and inserted into the water so that its lower end was positioned approximately 20 cm above the outlet of the bubble source. In such an arrangement the upper half of the tube was protruding above the water surface. It was connected via a pipe to a gear pump and water with bubbles was pumped through the tube at a constant flow rate of 200 or 300 ml min −1 . The same water was re-added into the tank through another pipe connected to the gear pump and fixed sufficiently far away from the bubble source so as not to disturb the bubble plume. The flow rates were deemed small enough compared to the total volume of the tank (approximately 100 l) to not create any significant background circulation.
The JAI Spark 5 MP camera was attached to the metallic frame around the tank and with a field of view through the broader side of the acrylic tube. It was equipped with two camera lenses: a Cannon 5X macro lens and microscope 10X lens. The focal distances of the camera lenses were 4 cm and approximately 2 cm, respectively, so that the camera had to be placed very close to the Perspex tube wall. The flow through the tube was back illuminated with a projector. A piece of polyester tracing film was attached at the rear wall of the Perspex tube to create a more uniform light distribution. The videos of the bubble flow through the tube were recorded for five minutes at a frame rate of 10 frames/second (3000 frames in total). The exposure time for each frame was 1 ms.
The JAI Spark SP-5000C-CXP2 5-megapixel camera possesses a Lince5M sensor with 2560 × 2048 pixels (width × height) and a cell size of 5.0 μm × 5.0 μm. This corresponds to the physical size of the sensor of 1.28 cm × 1.24 cm. The 5X (a) ( b) ( c) Figure 22. Panel (a) shows the microscope slide with ruled etchings recorded with the 10X lens. The distance between two adjacent shortest ticks is 10 μm and corresponds to 20 pixels in the horizontal direction, which confirms the resolution of 0.5 μm per pixel. Panels (b) and (c) are examples of recorded bubble images for the 5X and 10X lens, respectively. The corresponding parameters are B ≈ 30 cm 4 s −3 in both images as well as ρ a ≈ 1.008 g cm −3 and ρ a ≈ 1.18 g cm −3 , respectively. Bubbles that are found to be in focus are circled in blue.
magnification lens imaged a region 2.56 mm × 2.48 mm (with the resolution of 1.0 μm × 1.0 μm per pixel) at the focal plane and the 10X lens magnifies a region 1.28 mm × 1.24 mm (with the resolution of 0.5 μm × 0.5 μm per pixel). Figure 22(a) shows a microscope slide with ruled etchings that was imaged with the 10X lens. Due to the large aperture of these lenses, the depth of field was limited to 0.1 cm for the 5X lens, resulting in an error of less than 1 % for bubbles located slightly off the focal plane. The flow rate of 200 ml min −1 was used and the camera was equipped with the 10X lens whenever the bubble source was placed in the ambient water of density ρ a = 1.179 g cm −3 . To measure the bubble size distribution in the ambient solution of density ρ a = 1.008 g cm −3 , the gear pump was set to the flow rate of 300 ml min −1 and the 5X lens was used.
A bubble which is in focus of the camera appears on a recorded image as a dark circle with a sharp boundary and a white dot in the centre against a lighter background. The frames of the captured videos were processed as individual images. The pumped flow of 200 or 300 ml min −1 created a vertical velocity in the tube of approximately 4 or 6 mm s −1 , respectively. We processed every third frame of our recordings to avoid any potential double counting of bubbles. For each processed image, we used the imfindcircles function of MATLAB ® 2016a to determine the position and the radius of recorded bubbles. The accuracy of the imfindcircles function is limited for circles of radius less than 5 pixels. To reduce the likelihood of errors due to noise, we set the minimum radius of bubbles to be found to 15 pixels (which corresponds to 15 μm for the 5X lens and 7.5 μm for the 10X lens). The imfindcircles function is based on a circular Hough transform algorithm that requires the user to specify a lower threshold on the gradient of the light intensity at the edge of a potential circle. We found empirically that setting this threshold to 0.35 for the recordings with the 5X lens and to 0.1 for the recordings with the 10X lens was an acceptable choice for the algorithm to reliably find the bubbles that were visually perceived as being in focus of the camera. Figures 22(b) and 22(c) show examples of the images with the bubbles for the 5X and 10X lens, respectively.
In each of the recorded videos, a sample of at least 2000 bubbles was found (some of the samples counted as many as 6000-8000 bubbles). The probability distribution of the bubble radii r is reasonably well approximated by a log-normal distribution f r (r), i.e. ln f r (r) ∼ N (μ, σ ), where N(μ, σ ) is a normal distribution with the mean value μ and the standard deviation σ . Figure 23(a,b) shows the fitted log-normal probability distribution functions f r (r) of the bubble radii as solid lines for different source buoyancy fluxes B (colours) and for ambient water densities ρ a ≈ 1.179 g cm −3 (23a) and  Table 2 also includes the values for the bubble radii r mb that contain most of the buoyancy flux and the corresponding slip velocities u s .
1.008 g cm −3 (23b). The measured distributions of the bubble radii are plotted as symbols in the same panels. Thereby, a measured distribution of the bubble radii is obtained by counting the number of bubbles in a bin divided by the total number of bubbles in a sample and divided by the width of the bin. We use bins of approximately 1 μm width for the 5X lens and of 0.5 μm width for the 10X lens. Table 2 lists the values of μ and σ for the fitted log-normal distributions. The bubble radius corresponding to the peak for each distribution is denoted by r m and is determined as exp(μ − σ 2 ). Figure 23(c,d) displays the buoyancy distribution functions f b (r) as solid lines for different B and ρ a values. The buoyancy distribution function f b (r) is calculated from the bubble radius distribution function f r (r) by multiplying the latter by r 3 and normalising it (solid lines) The measured buoyancy distributions are also shown as symbols in the same panels. The maximum locations of the buoyancy distribution functions are labelled as r mb and their values are recorded in table 2. Appendix C. Image analysis Figure 24 shows an example of the experimental images of the plume at each stage of the post-processing for B = 50 cm 4 s −3 , ρ a ≈ 1.008 g cm −3 (series F) and Ω = 0 rad s −1 . A raw image of the plume is displayed in figure 24(a). The transformed plume image with the removed background (see (2.4)) is presented in figure 24(b). Figure 24(c) shows the corresponding binary image if the threshold is set toÎ < 0.995. In particular, the binary image enhances the visibility of bubbles around the plume tip. Figure 24(d) shows the detected edges for the plume image. The detected edges are displayed in red. In each pixel row of the image, the leftmost pixel where an edge has been detected is considered to be the left bounding edge of the plume and the rightmost edge pixel is defined as the right bounding edge of the plume. The region between the left and the right bounding edges is now defined as the extent of the plume and at each height, the width of the plume is the distance between the left and the right bounding edge. The left and the right bounding edges are also shown in figure 24(d) in green and blue, respectively.
Appendix D. Rise height of the plume silhouette centroid experiments, in a log-log plot. Initially, for each series the data collapse onto a power-law curve (see also figure 8). However, h c subsequently deviates from this power-law growth and rises slower for Ωt 3.