Structure of a dense release produced by varying initial conditions

Buoyancy driven flows such as gravity currents, present in nature and human made applications, are conveyors of particles or dissolved substances for long distances with clear implications for the environment. This transport depends on the triggering conditions of the current. Gravity currents are experimentally simulated under varying initial conditions by combining three different initial buoyancies and five volumes of dense fluid released. The horizontal and vertical structures of the gravity currents are analysed and it is shown that the variation on the initial configuration is conditioning for these. Vertical transport through the gravity current is influenced at the bottom by the solid wall over which the current flows, and at the upper interface by the contact with the ambient water. The relative contribution of shear stress at the bottom and at the upper interfaces are estimated and analysed in terms of the initial triggering conditions of the current; these two compete with the buoyancy, the driver of the current, determining mixing and entrainment. By using a proper parametrization, which accounts for both initial volume of release and location of the observation position relative to the lock, a relation between the resistance at the bottom and at the upper interfaces with the initial conditions of release (i.e. the lock-length) has been found; this is found to be independent of the initial density in the lock. The present study shows that the variation of the initial conditions have consequences on (1) the configuration of the currents and on (2) the hydrodynamics of the currents, including mass and momentum exchanges, which are in addition mutually dependent.     

Gravity currents in two-layer stratified media

An analytical, experimental and numerical study of boundary gravity currents propagating through a two-layer stratified ambient of finite vertical extent is presented. Gravity currents are supposed to originate from a lock-release apparatus; the (heavy) gravity current fluid is assumed to span the entire channel depth, H, at the initial instant. Our theoretical discussion considers slumping, supercritical gravity currents, i.e. those that generate an interfacial disturbance whose speed of propagation matches the front speed, and follows from the classical analysis of Benjamin (J Fluid Mech 31:209?C248, 1968). In contrast to previous investigations, we argue that the interfacial disturbance must be parameterized so that its amplitude can be straightforwardly determined from the ambient layer depths. Our parameterization is based on sensible physical arguments; its accuracy is confirmed by comparison against experimental and numerical data. More generally, measured front speeds show positive agreement with analogue model predictions, which remain strictly single-valued. From experimental and numerical observations of supercritical gravity currents, it is noted that this front speed is essentially independent of the interfacial thickness, ??, even in the limiting case where ?? = H so that the environment is comprised of a uniformly stratified ambient with no readily discernible upper or lower ambient layer. Conversely, when the gravity current is subcritical, there is a mild increase of front speed with ??. Our experiments also consider the horizontal distance, X, at which the front begins to decelerate. The variation of X with the interface thickness and the depths and densities of the ambient layers is discussed. For subcritical gravity currents, X may be as small as three lock lengths whereas with supercritical gravity currents, the gravity current may travel long distances at constant speed, particularly as the lower layer depth diminishes.     

A review of gravity currents formed by submerged single-port discharges in inland and coastal waters

In this paper, the authors review the current state of the science on the dynamics of gravity currents generated by positively and negatively buoyant jet discharges from submerged round outfalls (i.e., a point source) in inland and coastal waters. Specifically, this article focuses on describing gravity currents occurring at both the bottom boundary and the free surface of the receiving fluid. The manmade discharge operations generating both types of gravity currents and their significance to sustainability of the surrounding hydro-environment are first described. The authors then summarize the flow regimes characteristics of these discharges before becoming gravity currents and how those flow regimes influence the dynamics of the gravity currents. The gravity current dynamics in the calm receiving waters are then analyzed. This analysis is followed by an analysis of the influence of the hydrodynamic forces (e.g., currents, turbulence, waves) on the dynamics of gravity currents. Finally, the authors review quantitative modeling approaches for different forms of gravity current, and identify the current knowledge gaps and research needs.     

Turbulent entrainment into fluid mud gravity currents

The entrainment of ambient water into non-Newtonian fluid mud gravity currents was investigated in this study. Constant volume release gravity currents were generated in a lock-exchange tank for a wide range of experimental conditions. A technique similar to the so-called light attenuation technique was used to find the boundary of the current, allowing for the calculation of both temporal and bulk entrainment parameters (in terms of the temporal and bulk entrainment velocities, respectively). It was found that the temporal entrainment velocity is dependent on different parameters in the different propagation phases. The slumping phase begins with an adjustment zone (henceforth, non-established zone) in which the temporal entrainment velocity is not a function of the current front velocity, followed by the established zone in which the temporal entrainment velocity is a function of the current front velocity. This dependence of the temporal entrainment velocity on the current front velocity carries through to the inertia-buoyancy phase. As expected, temporal entrainment velocity in the viscous-buoyancy phase was negligible in comparison to average entrainment velocity in the other phases. It is observed that the temporal entrainment characteristics in the non-established zone is governed by the competition between the entrainment-inhibiting density stratification effects and the entrainment-favouring effects of the Kelvin–Helmholtz billows that are quantified by the Richardson number and the Reynolds number of the gravity current, respectively. In the established zone, Reynolds number effects were observed to dominate over Richardson number effects in dictating temporal entrainment characteristics. A parameterization for the temporal entrainment velocity for non-Newtonian fluid mud gravity currents is developed based upon the experimental observations. This study also found that the bulk entrainment characteristics for the non-Newtonian fluid mud gravity currents can be parameterized by the Newtonian bulk entrainment parameterizations that rely solely on a bulk Richardson number. Interestingly, it was found that the non-Newtonian characteristics of the gravity current have little to no effect on the entrainment of the Newtonian ambient fluid.     

Gravity currents and intrusions of stratified fluids into a stratified ambient

We consider high-Reynolds-number Boussinesq gravity currents and intrusions systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a rectangular lock; the height ratio of the fluids H, and the stratification parameter of the ambient S, are quite general. We develop a one-layer shallow-water (SW) model which is an extension of previously used and tested formulations for currents and intrusions of constant density. The internal stratification enters as a new dimensionless parameter, s ? [0,1]{\sigma \in [0,1]} . Analytical results are obtained for the initial “slumping” stage during which the speed of propagation is constant, and finite-difference solutions are presented for the more general time-dependent motion. Overall, this is a versatile and robust self-contained prediction tool, which reduces smoothly to the classical case when σ = 0. We show that, in general, the speed of propagation decreases when the internal stratification becomes more pronounced (σ increases). An interesting non-expected behavior was detected: when the stratification of the ambient is weak and moderate then the height of the current decreases with σ, but the opposite occurs when the stratification of the ambient is strong (S ≈ 1, including the case of an intrusion). Moreover, when the stratification of the ambient is strong a current with internal stratification may “run out” of driving power. We also consider the Benjamin-type steady state current with internal linear stratification in a non-stratified ambient, and show that an analytical solution exists, and that the maximal thickness decreases to below half-channel depth when σ increases.     

Wave mixing rise inferred from Lyapunov exponents

We study the horizontal surface mixing and the transport induced by waves in a coastal environment. A comparative study is addressed by computing the Lagrangian coherent structures, via Finite Size Lyapunov Exponents, that arise in two different numerical settings: with and without wave coupled to currents. In general, we observe that mixing is increased in the area due to waves. Besides, the methodology presented here is tested by deploying a set of eight Lagrangian drifters at different locations. This dynamical approach is shown as a valuable tool to extract information about transport, mixing and residence embedded in the Eulerian time dependent velocity fields obtained from numerical models.     

The flow of high-Reynolds axisymmetric gravity currents of a stratified fluid into a stratified ambient: shallow-water and box model solutions

We consider the axisymmetric flow (in a full cylinder or a wedge) of high-Reynolds-number Boussinesq gravity currents and intrusions systems in which both the ambient and the propagating “current” are linearly stratified. The main focus is on a current of fixed volume released from a cylinder lock; the height ratio of the fluids H, and the stratification parameter of the ambient S, are quite general. We develop a one-layer shallow-water model. The internal stratification enters as a new dimensionless parameter, ${\sigma \in [0, 1]}$ . In general, the time-dependent motion is obtained by standard finite-difference solutions; a self-similar analytical solution exists for S?= 0. We show that, in general, the speed of propagation decreases when the internal stratification becomes more pronounced (σ increases). We also developed a box-model approximation, and show that the resulting radius of propagation is in good agreement with the more rigorous shallow-water prediction.     

Mixing dynamics of turbidity currents interacting with complex seafloor topography

Direct Numerical Simulations are employed to investigate the mixing dynamics of turbidity currents interacting with seamounts of various heights. The mixing properties are found to be governed by the competing effects of turbulence amplification and enhanced dissipation due to the three-dimensional topography. In addition, particle settling is seen to play an important role as well, as it affects the local density stratification, and hence the stability, of the current. The interplay of these different mechanisms results in the non-monotonic dependence of the mixing behavior on the height of the seamount. Regions of dilute lock fluid concentration generally mix more intensely as a result of the seafloor topography, while concentrated lock fluid remains relatively unaffected. For long times, the strongest mixing occurs for intermediate bump heights. Particle settling is seen to cause turbidity currents to mix more intensely with the ambient than gravity currents.     

Generation of internal waves by sheared turbulence: experiments

A series of experiments is presented that model the generation of internal gravity waves in the ocean by the forcing of turbulent eddies in the surface mixed layer.The experimental setup consists of a shallow mixed upper layer and a deep continuously stratified lower layer. A source of turbulence is dragged through the upper layer. Internal waves can freely propagate in the lower layer. The internal waves are measured using synthetic schlieren to determine the frequencies of the generated waves. Consistent with other studies, it is found that the characteristic frequencies of internal waves generated by turbulence are an approximate constant fraction of the buoyancy frequency.     

3D numerical modelling of turbidity currents

During floods, the density of river water usually increases due to a subsequent increase in the concentration of the suspended sediment that the river carries, causing the river to plunge underneath the free surface of a receiving water basin and form a turbidity current that continues to flow along the bottom. The study and understanding of such complex phenomena is of great importance, as they constitute one of the major mechanisms for suspended sediment transport from rivers into oceans, lakes or reservoirs. Unlike most of the previous numerical investigations on turbidity currents, in this paper, a 3D numerical model that simulates the dynamics and flow structure of turbidity currents, through a multiphase flow approach is proposed, using the commercial CFD code FLUENT. A series of numerical simulations that reproduce particular published laboratory flows are presented. The detailed qualitative and quantitative comparison of numerical with laboratory results indicates that apart from the global flow structure, the proposed numerical approach efficiently predicts various important aspects of turbidity current flows, such as the effect of suspended sediment mixture composition in the temporal and spatial evolution of the simulated currents, the interaction of turbidity currents with loose sediment bottom layers and the formation of internal hydraulic jumps. Furthermore, various extreme cases among the numerical runs considered are further analyzed, in order to identify the importance of various controlling flow parameters.     

Material dispersion by oceanic internal waves

Internal gravity waves that are generated in the open ocean have a universal frequency spectrum, called Garrett–Munk spectrum. By initializing internal waves that satisfy the Garrett–Munk spectrum in a non-hydrostatic numerical model, we investigate the material dispersion produced by these internal waves. Three numerical experiments are designed: Exp.-1 uses a linearly stratified fluid, Exp.-2 has an upper mixed layer, and Exp.-3 incorporates a circular front into the upper mixed layer. Resorting to neutrally buoyant particles, we investigate the dispersion in terms of metrics of the relative dispersion and finite-scale Lyapunov exponent (FSLE). Exp.-1 shows that the dispersion regime produced by these internal waves is between ballistic and diffusive based on relative dispersion, and is however ballistic according to FSLE. The maximum FSLE at scales of 100 m is about 5 day\(^{-1}\), which is comparable to that calculated using ocean drifters. Exp.-2 demonstrates that internal waves can generate flows and material dispersion in an upper mixed layer. However, when mixed layer eddies are present, as in Exp.-3, the dispersion in the mixed layer is controlled by the eddies. In addition, we show that inertial oscillations do not affect the relative dispersion, but impact FSLE at scales of inertial oscillations.     

Lock-exchange gravity currents over rough bottoms

In nature, density driven currents often flow over or within a bottom roughness: a sea breeze encountering tall buildings, a shallow flow encountering aquatic vegetation, or a dense oceanic current flowing over a rough bottom. Laboratory experiments investigating the mechanisms by which bottom roughness enhances or inhibits entrainment and dilution in a lock-exchange dense gravity current have been conducted. The bottom roughness has been idealized by an array of vertical, rigid cylinders. Both spacing (sparse vs. dense configuration) and height of the roughness elements compared with the height of the current have been varied. Two-dimensional density fields have been obtained. Experimental results suggest that enhancement of the entrainment/dilution of the current can occur due to two different mechanisms. For a sparse configuration, the dense current propagates between the cylinders and the entrainment is enhanced by the vortices generated in the wake of the cylindrical obstacles. For a dense configuration, the dense current rides on top of the cylinders and the dilution is enhanced by the onset of convective instability between the dense current above the cylinders and the ambient lighter water between the cylinders. For low values of the ratio of the cylinder to lock height \(\lambda \) the dense current behavior approaches that of a current over a smooth bottom, while the largest deviations from the smooth bottom case are observed for large values of \(\lambda \).     

Visualization of the internal flow properties and the material exchange interface in an entraining viscous Newtonian gravity current

In order to simulate a simple entraining geophysical flow, a viscous Newtonian gravity current is released from a reservoir by a dam-break and flows along a rigid horizontal bed until it meets a layer of entrainable material of finite depth, identical to the current. The goal is to examine the entrainment mechanisms by observing the interaction between the incoming flow and the loose bed. The sole parameter varied is the initial volume of the gravity current, thus altering its height and velocity. The gravity current plunges or spills into the entrainable bed and the velocity of the flow front becomes linear with time. The bed material is directly affected: motion is generated in the fluid far downstream of, and in that lying beneath the encroaching front. Shear bands are identified, separating horizontal flow downstream from flow with a strong vertical component close to the step. Downstream of the step the flow is horizontal and stratified, with no slip on the bottom boundary and very low shear near the surface. Between these two regions may lie transitional zones with linear velocity profiles, separated by horizontal bands of high shear; the number of transitional zones in the cross-section varies with the initial volume of the dam-break.     

Experimental investigation of cross-shore sandbar volumes

Sandbars are critical to the cross-shore movement of sediment. Prediction of cross-shore sandbar volumes requires knowledge about the functional relationship of sediment transport rate conditions with waves, currents, base slope, sediment property and water depth. In this study, experiments on cross- shore sediment transport were carried out in a laboratory wave channel for initial base slopes of 1/8, 1/10 and 1/15. Using regular waves with different deep-water wave steepness generated by a pedal-type wave generator, bar volumes caused by cross-shore sediment transport are investigated for beach materials with the medium diameter of d₅₀?=?0.25, 0.32, 0.45, 0.62 and 0.80 mm. A non-dimensional equation for sandbar volume was obtained by using linear and non-linear regression methods through the experimental data and was compared with previously developed equations in the literature. The results have shown that the experimental data fitted well to the proposed equation with respect to the previously developed equations.     

Numerical simulation of particle-driven gravity currents

Particle-driven gravity currents frequently occur in nature, for instance as turbidity currents in reservoirs. They are produced by the buoyant forces between fluids of different density and can introduce sediments and pollutants into water bodies. In this study, the propagation dynamics of gravity currents is investigated using the FLOW-3D computational fluid dynamics code. The performance of the numerical model using two different turbulence closure schemes namely the renormalization group (RNG) ${k-\epsilon}$ scheme in a Reynold-averaged Navier-Stokes framework (RANS) and the large-eddy simulation (LES) technique using the Smagorinsky scheme, were compared with laboratory experiments. The numerical simulations focus on two different types of density flows from laboratory experiments namely: Intrusive Gravity Currents (IGC) and Particle-Driven Gravity Currents (PDGC). The simulated evolution profiles and propagation speeds are compared with laboratory experiments and analytical solutions. The numerical model shows good quantitative agreement for predicting the temporal and spatial evolution of intrusive gravity currents. In particular, the simulated propagation speeds are in excellent agreement with experimental results. The simulation results do not show any considerable discrepancies between RNG ${k-\epsilon}$ and LES closure schemes. The FLOW-3D model coupled with a particle dynamics algorithm successfully captured the decreasing propagation speeds of PDGC due to settling of sediment particles. The simulation results show that the ratio of transported to initial concentration C _o /C _i by the gravity current varies as a function of the particle diameter d _s . We classify the transport pattern by PDGC into three regimes: (1) a suspended regime (d _s is less than about 16 μm) where the effect of particle deposition rate on the propagation dynamics of gravity currents is negligible i.e. such flows behave like homogeneous fluids (IGC); (2) a mixed regime (16 μm < d _s <40 μm) where deposition rates significantly change the flow dynamics; and (3) a deposition regime (d _s ?> 40 μm) where the PDGC rapidly loses its forward momentum due to fast deposition. The present work highlights the potential of the RANS simulation technique using the RNG ${k-\epsilon}$ turbulence closure scheme for field scale investigation of particle-driven gravity currents.     

Preliminary results from studies of nocturnal bioluminescence with subsurface moored photometers

Spatial and temporal patterns of bioluminescent flashes were recorded from fall 1982 through Spring 1983 by photometers moored offshore in Scripps Canyon, La Jolla, California, USA. From depths between 8 and 90 m, realtime data were transmitted by cable to a laboratory on land approximately one mile (1.7 km) away. In addition, temperature, depth and current velocity and direction were monitored either in real time by direct coupling to a laboratory-based system, or by internal data storage-systems that were retrieved at regular intervals and subsequently analyzed. Our studies showed that our field station is largely uncoupled from wave action effects usually associated with luminescence measurements made from ships. Bioluminescent activity varied greatly both during a single night and between different nights. Vertical profiling of the water column between 8 and 90 m showed evidence of vertical migration, patchiness of distribution and large-scale spatial differences in total bioluminescent activity. Currents had a major impact on patterns of bioluminescent activity; however, sometimes high levels of luminescence were recorded in complete absence of currents. Diel cycles, organism patchiness, the level of downwelling ambient light, and currents appeared to interact in controlling the levels and patterns of bioluminescence.     

Estimating Climate Resilience for Conservation across Geophysical Settings

Conservationists need methods to conserve biological diversity while allowing species and communities to rearrange in response to a changing climate. We developed and tested such a method for northeastern North America that we based on physical features associated with ecological diversity and site resilience to climate change. We comprehensively mapped 30 distinct geophysical settings based on geology and elevation. Within each geophysical setting, we identified sites that were both connected by natural cover and that had relatively more microclimates indicated by diverse topography and elevation gradients. We did this by scoring every 405 ha hexagon in the region for these two characteristics and selecting those that scored >SD 0.5 above the mean combined score for each setting. We hypothesized that these high‐scoring sites had the greatest resilience to climate change, and we compared them with sites selected by The Nature Conservancy for their high‐quality rare species populations and natural community occurrences. High‐scoring sites captured significantly more of the biodiversity sites than expected by chance (p < 0.0001): 75% of the 414 target species, 49% of the 4592 target species locations, and 53% of the 2170 target community locations. Calcareous bedrock, coarse sand, and fine silt settings scored markedly lower for estimated resilience and had low levels of permanent land protection (average 7%). Because our method identifies—for every geophysical setting—sites that are the most likely to retain species and functions longer under a changing climate, it reveals natural strongholds for future conservation that would also capture substantial existing biodiversity and correct the bias in current secured lands.     

Mechanisms of cross-shelf transport of crab megalopae inferred from a time series of daily abundance

From 20 July 1982 to 19 July 1984, crab megalopae were trapped daily from the Scripps Institution of Oceanography pier. Pachygrapsus crassipes, Portunus xantusii, Cancer spp., Hemigrapsus spp., and Majid crab megalopae dominated the catch. Yearly, seasonal, and daily variations in the magnitude of the catch were observed. Yearly and seasonal variations were probably due to a strong El Niño event that occurred during the study and to the timing of spawning and duration of the larval phase, respectively. Daily variation was correlated with oceanographic processes that can transport larvae to shore. Catch of some taxon during some seasons correlated with wind stress suggesting that transport was wind driven. The correlations were, however, weak and the sign of the correlation varied between years. The maximum daily tidal range was significantly correlated (cross-correlations and cross Fourier analysis) to both daily seawater temperature anomalies (surface and bottom) and daily catch of crab megalopae in all taxa enumerated. Significant correlations between tidal range and temperature anomalies suggest that temperature anomalies were primarily due to the shoreward transport of warm and cold water by the internal tides. The consistent and relatively strong relationship between tidal range and catch of megalopae (the cross-Fourier analysis suggests that from 20 to >90% of the variation in catch can be attributed to variation in the tidal range) suggest that much of the shoreward transport of megalopae was via the internal tides. Shoreward transport of larvae by internal tides may be due to internal cold bores or convergences over large tidally generated internal waves (solitons). Peak catches of megalopae, however, were often not associated with cold anomalies suggesting that transport was due to moving convergences over internal waves.