Implementation of a wetting-and-drying model in simulating the Kennebec–Androscoggin plume and the circulation in Casco Bay

A high-resolution coastal ocean model was developed to simulate the temporal/spatial variability of the Kennebec–Androscoggin (K–A) river plume and the circulation in Casco Bay. The model results agree favorably with the moored and shipboard observations of velocity, temperature, and salinity. The surface salinity gradient was used to distinguish the plume from the ambient coastal water. The calculated plume thickness suggests that the K–A plume is surface trapped. Its horizontal scales correlate well with Q ^0.25, where Q is the volume discharge of the rivers. Directional spreading is affected by the wind with the upwelling favorable wind transporting the plume water offshore. Both the wind and the tide also enhance mixing in the plume. The inclusion of a wetting-and-drying (WAD) scheme appears to enhance the mixing and entrainment processes near the estuary. The plume becomes thicker near the mouth of the estuary, the outflow velocity of the plume is weaker, and the radius of the river plume shrinks. The flow field in the model run with the WAD is noisier, not only in shallow areas of Casco Bay but also in the plume and even on the shelf. We speculate that the WAD processes can affect much larger areas than the intertidal zones, especially via a river plume that feeds into a coastal current.     

The dual-beam mini-DOAS technique—measurements of volcanic gas emission,plume height and plume speed with a single instrument

The largest error in determining volcanic gas fluxes using ground based optical remote sensing instruments is typically the determination of the plume speed, and in the case of fixed scanning instruments also the plume height. We here present a newly developed technique capable of measuring plume height, plume speed and gas flux using one single instrument by simultaneously collecting scattered sunlight in two directions. The angle between the two measurement directions is fixed, removing the need for time consuming in-field calibrations. The plume height and gas flux is measured by traversing the plume and the plume speed is measured by performing a stationary measurement underneath the plume. The instrument was tested in a field campaign in May 2005 at Mt. Etna, Italy, where the measured results are compared to wind fields derived from a meso-scale meteorological model (MM5). The test and comparison show that the instrument is functioning and capable of estimating wind speed at the plume height.     

Particle recycling in volcanic plumes

We have developed a new theoretical model of an eruption column that accounts for the re-entrainment of particles as they fall out of the laterally spreading umbrella cloud. The model illustrates how the mass flux of particles in the plume may increase with height in the plume, by a factor as large as 2.5 because of this recycling. Three important consequences are that (1) the critical velocity required to generate a buoyant eruption column for a given mass flux increases, (2) the total height of rise of the column may decrease, and (3) we infer that in relatively wind-free environments, for eruption columns near the conditions for collapse, the recycling of particles may lead to an unsteady oscillating motion of the plume, which, in time, may lead to the formation of interleaved fall and flow deposits.     

Development and verification of deep-water blowout models

Modeling of deep-water releases of gas and oil involves conventional plume theory in combination with thermodynamics and mass transfer calculations. The discharges can be understood in terms of multiphase plumes, where gas bubbles and oil droplets may separate from the water phase of the plume and rise to the surface independently. The gas may dissolve in the ambient water and/or form gas hydrates––a solid state of water resembling ice. All these processes will tend to deprive the plume as such of buoyancy, and in stratified water the plume rise will soon terminate. Slick formation will be governed by the surfacing of individual oil droplets in a depth and time variable current. This situation differs from the conditions observed during oil-and-gas blowouts in shallow and moderate water depths. In such cases, the bubble plume has been observed to rise to the surface and form a strong radial flow that contributes to a rapid spreading of the surfacing oil. The theories and behaviors involved in deepwater blowout cases are reviewed and compared to those for the shallow water blowout cases.     

A trigger mechanism for the Lake Nyos disaster

The catastrophic release of carbon dioxide gas from Lake Nyos on 21 August 1986 is discussed in the context of the buoyancy reversal instability. Originally proposed by Randall (1980) [Randall, D.S., 1980. Conditional instability of the first kind upside-down. Journal of Atmospheric Sciences 37: 125–130.] and Deardoff (1980) [Deardoff, J.W., 1980. Cloud-top entrainment instability. Journal of Atmospheric Sciences 37: 131–147.] for the `cloud-top entrainment instability' of stratocumulus clouds, the buoyancy reversal instability has been studied experimentally in water tank experiments by Shy and Breidenthal (1990) [Shy, S.S., Breidenthal, R.E., 1990. Laboratory experiments on the cloud-top entrainment instability. Journal of Fluid Mechanics 214: 1–15.], who identified three criteria for instability. The initial disturbance must be sufficiently large, so that its Reynolds number is above the mixing transition, its Richardson number must be less than one to achieve overturning and mixing, and the buoyancy reversal parameter must be greater than a critical value, of order one. The implications and applicability of these criteria to Lake Nyos are discussed. The criterion for the Reynolds number is easily satisfied for typical wind velocities in the Lake Nyos region. Similarly, the Richardson number based on incident turbulence is easily less than unity, and therefore satisfy the second criterion for instability. In the case of Lake Nyos, the continuous release of carbon dioxide at the bottom of the lake increases the value of the buoyancy reversal parameter until it reaches its critical value, at which point an explosion occurs. This instability provides a plausible trigger for the 1986 explosion. After each explosion, the buoyancy reversal parameter returns to below its critical value, only to slowly rise again over time, as CO₂ continues to enter the lake, setting the stage for the next explosion. Future explosions may be avoided if the value of the buoyancy reversal parameter is prevented from approaching its critical value by artificial mixing at the thermocline, such as with a bubble plume.     

A model study of tide- and wind-induced mixing in the Columbia River Estuary and plume

A numerical simulation of circulation in the Columbia River estuary and plume during the summer of 2004 is used to explore the mixing involved as river water is transformed into shelf water. The model is forced with realistic river flow, tides, wind stress, surface heat flux, and ocean boundary conditions. Simulated currents and water properties on the shelf near the mouth are compared with records from three moorings (all in 72 m of water) and five CTD sections. The model is found to have reasonable skill; statistically significant correlations between observed and modeled surface currents, temperature, and salinity are all 0.42–0.72 for the mooring records. Equations for the tidally averaged, volume-integrated mechanical energy budget (kinetic and potential) are derived, with attention to the effects of: (i) Reynolds averaging, (ii) a time varying volume due to the free surface, and (iii) dissipation very close to the bottom. It is found that convergence of tidal pressure work is the most important forcing term in the estuary. In the far field plume (which has a volume 15 times greater than that of the estuary), the net forcing is weaker than that in the estuary, and may be due to either tidal currents or wind stress depending on the time period considered. These forcings lead to irreversible mixing of the stratification (buoyancy flux) that turns river water into shelf water. This occurs in both the plume and estuary, but appears to be more efficient (17% vs. 5%), and somewhat greater (4.2 MW vs. 3.3 MW), in plume vs. estuary. This demonstrates the importance of both wind and tidal forcing to watermass transformation, and the need to consider the estuary and plume as part of a single system.     

The spread of subducted lithospheric material along the mid-mantle boundary

Studies of phase transitions in silicate minerals at high temperatures and pressures suggest that the bulk density of subducted lithosphere at the mid-mantle boundary is intermediate between the densities of the upper and lower mantle. We argue that, if this is the case, then the lithospheric material will intrude along the mid-mantle boundary driven by buoyancy forces resulting from the compositional density differences between the intrusion and its surroundings. The rate of spread of the intrusion is given by a balance between these buoyancy forces and the viscous resistance of the mantle to motion. Using results from our recent studies of the fluid mechanics of such viscous gravity currents, we find that lithospheric material can propagate between one thousand and two thoudand kilometres in a hundred million years and can cover the entire boundary in one to six billion years. This spreading may be reflected in the global distribution of the isotopic characteristics of oceanic basalts.     

Gravity currents down canyons: effects of rotation

The flow of dense water in a V-shaped laboratory-scale canyon is investigated by using a non-hydrostatic numerical ocean model with focus on the effects of rotation. By using a high-resolution model, a more detailed analysis of plumes investigated in the laboratory (Deep-Sea Res I 55:1021–1034 2008) for laminar flow is facilitated. The inflow rates are also increased to investigate plume structure for higher Reynolds numbers. With rotation, the plumes will lean to the side of the canyon, and there will be cross-canyon geostrophic currents and Ekman transports. In the present study, it is found that the cross-canyon velocities are approximately 5 % of the down-canyon velocities over the main body of the plume for the rotational case. With rotation, the flow of dense water through the body of the plume and into the plume head is reduced. The plume head becomes less developed, and the speed of advance of the head is reduced. Fluid parcels near the top of the plume will to a larger extent be left behind the faster flowing dense core of the plume in a rotating system. Near the top of the plume, the cross-canyon velocities change direction. Inside the plume, the cross-flow is up the side of the canyon, and above the interface to the ambient there is a compensating cross-flow down the side of the canyon. This means that parcels of fluid around the interface become separated. Parcels of fluid around the interface with small down-canyon velocity components and relative large cross-canyon components will follow a long helix-like path down the canyon. It is found that the entrainment coefficients often are larger in the rotational experiments than in corresponding experiments without rotation. The effects of rotation and higher inflow rates on the areal patterns of entrainment velocities are demonstrated. In particular, there are bands of higher entrainment velocities along the lateral edges of the plumes in the rotational cases.     

Temperature, density and buoyancy fluxes in “black smoker” plumes, and the criterion for buoyancy reversal

Measurements of the temperature and composition of effluent from vents on the sea floor can be used to deduce the in-situ density of this fluid, which is required for calculations of flow in the chimneys and through their porous walls. This density is, however, not directly relevant when calculating the buoyancy flux in the plume above a smoker. It is the asymptotic buoyancy flux, following extensive dilution with seawater, which is required when estimating the height of rise of plumes in a stably stratified ocean, and when calculating the criterion for reversal of buoyancy due to non-linear mixing effects. The results of mixing calculations show that the effluent from hydrothermal vents on the sea floor will exhibit reversing buoyancy if the ejected fluid has a temperature of 300°C and a salinity greater than 8 wt.% NaCl. If the temperature of the effluent is 200°C the salinity required for reversing buoyancy falls to 5.5 wt.% NaCl. Measurements of temperature and salinities of sea-floor hydrothermal fluid suggest that fluids with the characteristics required to form reversing plumes are ejected at the sea floor. The possibility that reversing plumes may be found has important implications for the formation of massive sulfide deposits.     

Circulation induced by subglacial discharge in glacial fjords: Results from idealized numerical simulations

The flow caused by the discharge of freshwater underneath a glacier into an idealized fjord is simulated with a 2D non-hydrostatic model. As the freshwater leaves horizontally the subglacial opening into a fjord of uniformly denser water it spreads along the bottom as a jet, until buoyancy forces it to rise. During the initial rising phase, the plume meanders into complex flow patterns while mixing with the surrounding fluid until it reaches the surface and then spreads horizontally as a surface seaward flowing plume of brackish water. The process induces an estuarine-like circulation. Once steady-state is reached, the flow consists of an almost undiluted buoyant plume rising straight along the face of the glacier that turns into a horizontal surface layer thickening as it flows seaward. Over the range of parameters examined, the estuarine circulation is dynamically unstable with gradient Richardson number at the sheared interface having values of <1/4. The surface velocity and dilution factors are strongly and non-linearly related to the Froude number. It is the buoyancy flux that primarily controls the resulting circulation with the momentum flux playing a secondary role.     

The movement of oil under non-breaking waves

The combined effects of wave kinematics, turbulent diffusion, and buoyancy on the transport of oil droplets at sea were investigated in this work using random walk techniques in a Monte Carlo framework. Six hundred oil particles were placed at the water surface and tracked for 500 wave periods. A dimensionless formulation was presented that allowed reporting distances in terms of the wave length and times in terms of the wave period. Stokes' drift was, expectedly, the major mechanism for horizontal transport. It was also found that plumes that have large terminal rise velocities move faster forward but spread less than those that have small terminal rise velocities. The increase in wave slope (or wave steepness) caused an increase in transport and spreading of the plume. Our results supported treating the oil as completely mixed vertically in a layer near the surface. In the horizontal direction, buoyant plumes had spreading coefficients that are essentially constant after about 200 wave periods. But neutrally buoyant plumes had horizontal spreading coefficients that increased with time (for the simulation time of 500 wave periods). Techniques for generalizing the results for a wide range of wave parameters were presented.     

Box and Gaussian plume models of the exhaust composition evolution of subsonic transport aircraft in- and out of the flight corridor

A box and a Gaussian plume model including gas-phase photochemistry and with plume dispersion parameters estimated from the few available plume observations are proposed and used for evaluation of photochemical transformations of exhausts from a single subsonic transport aircraft. The effects of concentration inhomogeneities in the plume cross section on the photochemical sources and sinks in the plume are analyzed for various groups of compounds. The influence of these inhomogeneities on the rate and on the mass of ambient air entrainment into the plume are studied also by comparing the box and the Gaussian plume model simulations during the first hours of their life. Due to the enterance of HO_X and NO_X from ambient air into the plume with rates varying from the wind shear and turbulence conditions, the rate of emitted NOX oxidation in the plume is dependent on these and also on the background concentration levels of HO_X and NO_X.     

The effects of river discharge and seasonal winds on the shelf off southeastern South America

The Río de la Plata waters form a low salinity tongue that affects the circulation, stratification and the distributions of nutrients and biological species over a wide extent of the adjacent continental shelf. The plume of coastal waters presents a seasonal meridional displacement reaching lower latitudes (28°S) during austral winter and 32°S during summer. Historical data suggests that the wind causes the alongshore shift, with southwesterly (SW) winds forcing the plume to lower latitudes in winter while summer dominant northeasterly (NE) winds force its southward retreat. To establish the connection between wind and outflow variations on the distribution of the coastal waters, we conducted two quasi-synoptic surveys in the region of Plata influence on the continental shelf and slope of southeastern South America, between Mar del Plata, Argentina and the northern coast of Santa Catarina, Brazil. We observed that: (A) SW winds dominating in winter force the northward spreading of the plume to low latitudes even during low river discharge periods; (B) NE winds displace the plume southward and spread the low salinity waters offshore over the entire width of the continental shelf east of the Plata estuary. The southward retreat of the plume in summer leads to a volume decrease of low salinity waters over the shelf. This volume is compensated by an increase of Tropical waters, which dominate the northern shelf. The subsurface transition between Subantarctic and Subtropical Shelf Waters, the Subtropical Shelf Front, and the subsurface water mass distribution, however, present minor seasonal variations. Along shore winds also influence the dynamics and water mass variations along the continental shelf area. In areas under the influence of river discharge, Subtropical Shelf Waters are kept away from the coastal region. When low salinity waters retreat southward, NE winds induce a coastal upwelling system near Santa Marta Cape. In summer, solar radiation promotes the establishment of a strong thermocline that increases buoyancy and further enhances the offshore displacement of low salinity waters under the action of NE winds.     

Pumped-storage hydropower schemes: integral plume models for pumpback simulations

Integral plume rise formulations for submerged jets in reservoirs are analysed for implementation in models of pumped-storage hydropower schemes. Calculations of entrainment require adequate representation of length scales, cross-sectional profiles, the entrainment coefficient and the zone of flow establishment (ZFE). A new method of determining the values of the flux variables at the end of the ZFE is given and recommendations are made for improved representations of the merging of a row of point diffusers.     

A model study of the Copper River plume and its effects on the northern Gulf of Alaska

A three-level nested Regional Ocean Modeling System was used to examine the seasonal evolution of the Copper River (CR) plume and how it influences the along- and across-shore transport in the northern Gulf of Alaska (NGoA). A passive tracer was introduced in the model to delineate the growth and decay of the plume and to diagnose the spread of the CR discharge in the shelf, into Prince William Sound (PWS) and offshore. Furthermore, a model experiment with doubled discharge was conducted to investigate potential impacts of accelerated glacier melt in future climate scenarios. The 2010 and 2011 simulation revealed that the upstream (eastward) transport in the NGoA is negligible. About 60 % of the passive tracer released in the CR discharge is transported southwestward on the shelf, while another one third goes into PWS with close to 60 % of which exiting PWS to the shelf from Montague Strait. The rest few percent is transported across the shelf break and exported to the GoA basin. The downstream transport and the transport into PWS are strongly regulated by the downwelling-favorable wind, while the offshore transport is related to the accumulation of plume water in the shelf, frontal instability, and the Alaskan Stream. It takes weeks in spring for the buoyancy to accumulate so that a bulge forms outside of the CR estuary. The absence of strong storms as in the summer of 2010 allows the bulge continue growing to trigger frontal instability. These frontal features can interact with the Alaskan Stream to induce transport pulses across the shelf break. Alternatively as in 2011, a downwelling-favorable wind event in early August (near the peak discharge) accelerates the southwestward coastal current and produces an intense downstream transport event. Both processes result in fast drains of the buoyancy and the plume content, thereby rapid disintegration of the plume in the shelf. The plume in the doubled discharge case can be two to three times in size, which affects not only the magnitude but also the timing of certain transport events. In particular, the offshore transport increases by several folds because the plume appears to be more easily entrained by the seaward flow along the side of Hinchinbrook Canyon.     

Richardson number profiles in laboratory experiments applied to shallow seas

Abstract

A unified analysis is given of the critical conditions for the onset of stratification due to either a vertical or a horizontal buoyancy flux, with tidal or wind stirring.

The critical conditions for the onset of stratification with a horizontal buoyancy flux are found to be of the form of ratios of the tidal slope, or wind setup, to the equivalent surface slope due to the lateral density gradient. These ratios, which are easily determined from sea data, indicate that the profiles of critical flux Richardson Number, averaged over the stirring cycle, are similar to those inferred from the laboratory experiments of Hopfinger and Linden (1982) in which there is zero mean shear turbulence with a stabilising buoyancy flux, and also that the efficiency for the conversion of kinetic energy to potential energy for tidal stirring is similar to that for wind stirring.

The observed much greater efficiency for wind stirring, compared with tidal stirring with a vertical buoyancy flux, is also consistent with the existence of flux Richardson Number profiles in the sea similar to those occurring in the corresponding laboratory experiments. Using the solution of the turbulent kinetic energy equation for the water column, the relative importance of the production of turbulent kinetic energy, and its diffusion by turbulence are assessed, and the critical conditions for the onset of stratification with a vertical buoyancy flux are shown to reduce the classical Simpson—Hunter form.     

Upscaling of CO2 vertical migration through a periodic layered porous medium: The capillary-free and capillary-dominant cases

We present an upscaled model for the vertical migration of a CO₂ plume through a vertical column filled with a periodic layered porous medium. This model may describe the vertical migration of a CO₂ plume in a perfectly layered horizontal aquifer. Capillarity and buoyancy are taken into account and semi-explicit upscaled flux functions are proposed in the two following cases: (i) capillarity is the main driving force and (ii) buoyancy is the only driving force. In both cases, we show that the upscaled buoyant flux is a bell-shaped function of the saturation, as in the case of a homogeneous porous medium. In the capillary-dominant case, we show that the upscaled buoyant flux is the harmonic mean of the buoyant fluxes in each layer. The upscaled saturation is governed by the continuity of the capillary pressure at the interface between layers. In the capillary-free case, the upscaled buoyant flux and upscaled saturation are determined by the flux continuity condition at the interface. As the flux is not continuous over the entire range of saturation, the upscaled saturation is only defined where continuity is verified, i.e. in two saturation domains. As a consequence, the upscaled buoyant flux is described by a piecewise continuous function. Two analytical approximations of this flux are proposed and this capillary-free upscaled model is validated for two cases of heterogeneity. Upscaled and cell averaged saturations are in good agreement. Furthermore, the proposed analytical upscaled fluxes provide satisfactory approximations as long as the saturation set at the inlet of the column is in a range where analytical and numerical upscaled fluxes are close.     

Numerical modeling of the Amazon River plume

Marine circulation above the northern Brazilian continental shelf is subject to energetic forcing factors of various origins: high water buoyancy fluxes induced by the Amazon River freshwater discharge, a strong coastal current associated with a mesoscale current (North Brazil Current (NBC)), a forcing by semidiurnal tide and by Northeast or Southeast trade winds according to the season. Using a three-dimensional (3-D) hydrodynamic numerical model (MOBEEHDYCS), and realistic bathymetry and coastline of the northern Brazilian shelf, this paper aims at studying the influence of some specific physical processes on the morphology of the Amazon plume. The very large volume discharge (180 000 m³/s on average) and the weak effect of Coriolis force are additional characteristics of the studied system, which induce a particular dynamics. The various forcing factors are successively introduced into the model in order to simulate and to determine their respective influences upon the plume extent and the hydrodynamics at the shelf scale. Simulation reveal that the coastal current is at the origin of the permanent northwestward Amazon plume extension while wind effect can either reinforce or moderate this situation. The tide intervenes also to modify the position of the salinity front: a horizontal migration of salinity front is observed under its action.     

Dynamics of medium-intensity dense water plumes in the Arkona Basin, Western Baltic Sea

In this study, the dynamics of medium-intensity inflow events over Drogden Sill into the Arkona Sea are investigated. Idealised model simulations carried out with the General Estuarine Transport Model suggest that most of the salt transport during such inflow events occur north of Kriegers Flak, a shoal with less than 20 m water depth surrounded by water depths of more than 40 m. This assumption about the pathway is supported by recent ship-based observations in the Arkona Sea during a medium-intensity inflow event. The propagation of a saline bottom plume could be observed during several days after having passed Drogden Sill. In the area north of Kriegers Flak the plume was about 10 m thick, and propagated with more than 0.5 m s⁻¹ and a salinity of up to 20 psu (with ambient water salinity being 8 psu) eastwards. Although the model simulations were idealised, the structural agreement between the observation and model result was good. The structure and pathways of these medium-intensity inflow events are of specific interest due to the plans for erecting extensive offshore wind farms in the Arkona Sea which may under certain circumstances lead to increased entrainment of ambient water into the bottom plumes.