The rate of sedimentation from turbulent suspension: An experimental model with application to pyroclastic density currents and discussion on the grain‐size dependence of flow runout

Large‐scale experiments generating ground‐hugging multiphase flows were carried out with the aim of modelling the rate of sedimentation, of pyroclastic density currents. The current was initiated by the impact on the ground of a dense gas‐particle fountain issuing from a vertical conduit. On impact, a thick massive deposit was formed. The grain size of the massive deposit was almost identical to that of the mixture feeding the fountain, suggesting that similar layers formed at the impact of a natural volcanic fountain should be representative of the parent grain‐size distribution of the eruption. The flow evolved laterally into a turbulent suspension current that sedimented a thin, tractive layer. A good correlation was found between the ratio of transported/sedimented load and the normalized Rouse number of the turbulent current. A model of the sedimentation rate was developed, which shows a relationship between grain size and flow runout. A current fed with coarser particles has a higher sedimentation rate, a larger grain‐size selectivity and runs shorter than a current fed with finer particles. Application of the model to pyroclastic deposits of Vesuvius and Campi Flegrei of Southern Italy resulted in sedimentation rates falling inside the range of experiments and allowed definition of the duration of pyroclastic density currents which add important information on the hazard of such dangerous flows. The model could possibly be extended, in the future, to other geological density currents as, for example, turbidity currents.     

Contour current imprints and contourite drifts in the Bahamian archipelago

New data collected along the slopes of Little and Great Bahama Bank and the abyssal plain of the Bahama Escarpment provides new insights about contour current‐related erosive structures and associated deposits. The Bahamian slope shows abundant evidence of bottom current activity such as furrows, comet‐like structures, sediment waves and drifts. At a seismic scale, large erosion surfaces and main periods of drift growth resulted from current acceleration related to plate tectonic processes and progressive opening and closure of gateways and long‐term palaeoclimate evolution. At present‐day, erosion features and contourite drifts are either related to relatively shallow currents (<1000 m water depth) or to deep currents (>2500 m water depth). It appears that the carbonate nature of the drifts does not impact the drift morphology at the resolution addressed in the present study. Classical drift morphologies defined in siliciclastic environments are found, such as mounded, plastered and separated drifts. In core, contourite sequences show a bi‐gradational trend that resembles classical contourite sequences in siliciclastic deposits showing a direct relationship with a change in current velocity at the sea floor. However, in a carbonate system the peak in grain size is associated with increased winnowing rather than increased sediment supply as in siliciclastic environments. In addition, the carbonate contourite sequence is usually thinner than in siliciclastics because of lower sediment supply rates. Little Bahama Bank and Great Bahama Bank contourites contain open‐ocean input and slope‐derived debris from glacial episodes. Inner platform, platform edge and open ocean pelagic input characterize the classical periplatform ooze during interglacials. In all studied examples, the drift composition depends on the sea floor topography surrounding the drift location and the type of sediment supply. Carbonate particles are derived from either the slope or the platform in slope and toe of slope drifts, very deep contourites have distant siliciclastic sources of sediment supply. The recent discovery of the importance of a large downslope gravitary system along Bahamian slopes suggests frequent interactions between downslope and along‐slope (contour currents) processes. The interlayering of mass flow deposits and contourites at a seismic scale or the presence of surface structures associated with both contour currents and mass flow processes shows that both processes act at the same location. Finally, contour currents have an important impact on the repartition of deep‐water coral mounds. Currents can actively interact with mounds as a nutrient and oxygen supplier or have a passive interaction, with mounds solely being obstacles orienting erosion and deposition.     

Pyroclastic dune bedforms: macroscale structures and lateral variations. Examples from the 2006 pyroclastic currents at Tungurahua (Ecuador)

Pyroclastic currents are catastrophic flows of gas and particles triggered by explosive volcanic eruptions. For much of their dynamics, they behave as particulate density currents and share similarities with turbidity currents. Pyroclastic currents occasionally deposit dune bedforms with peculiar lamination patterns, from what is thought to represent the dilute low concentration and fluid‐turbulence supported end member of the pyroclastic currents. This article presents a high resolution dataset of sediment plates (lacquer peels) with several closely spaced lateral profiles representing sections through single pyroclastic bedforms from the August 2006 eruption of Tungurahua (Ecuador). Most of the sedimentary features contain backset bedding and preferential stoss‐face deposition. From the ripple scale (a few centimetres) to the largest dune bedform scale (several metres in length), similar patterns of erosive‐based backset beds are evidenced. Recurrent trains of sub‐vertical truncations on the stoss side of structures reshape and steepen the bedforms. In contrast, sporadic coarse‐grained lenses and lensoidal layers flatten bedforms by filling troughs. The coarsest (clasts up to 10 cm), least sorted and massive structures still exhibit lineation patterns that follow the general backset bedding trend. The stratal architecture exhibits strong lateral variations within tens of centimetres, with very local truncations both in flow‐perpendicular and flow‐parallel directions. This study infers that the sedimentary patterns of bedforms result from four formation mechanisms: (i) differential draping; (ii) slope‐influenced saltation; (iii) truncative bursts; and (iv) granular‐based events. Whereas most of the literature makes a straightforward link between backset bedding and Froude‐supercritical flows, this interpretation is reconsidered here. Indeed, features that would be diagnostic of subcritical dunes, antidunes and ‘chute and pools’ can be found on the same horizon and in a single bedform, only laterally separated by short distances (tens of centimetres). These data stress the influence of the pulsating and highly turbulent nature of the currents and the possible role of coherent flow structures such as Görtler vortices. Backset bedding is interpreted here as a consequence of a very high sedimentation environment of weak and waning currents that interact with the pre‐existing morphology. Quantification of near‐bed flow velocities is made via comparison with wind tunnel experiments. It is estimated that shear velocities of ca 0·30 m.s^?1 (equivalent to pure wind velocity of 6 to 8 m.s^?1 at 10 cm above the bed) could emplace the constructive bedsets, whereas the truncative phases would result from bursts with impacting wind velocities of at least 30 to 40 m.s^?1.     

Optimization of the Sefid-Roud Dam desiltation process using a sophisticated one-dimensional numerical model

Although water and soil conservation activities reduce reservoir sedimentation, it is inevitable that reservoirs fed by rivers transporting high amounts of sediment will experience sedimentation. The Ghezel-Ozan and Shah-Roud rivers, which flow to the Sefld-Roud reservoir dam, are both highly sediment-laden and transport significant amounts of sediment in both bed load and suspended load forms to the reservoir. Hence, it seems that the only practical way to remove the sediment from the reservoir is to flush it out using the Chasse method. In the present paper, field measurements of Chasse operation characteristics taken in previous years are presented, and a numerical model that simulates this process is introduced. After calibrating the model using field measured data, the calculated results （for reservoir pressure flushing and released sediment volume） of the numerical model were compared with other measured data for the same Chasse operation and the results agree well. Finally, using the numerical simulation results, the best approaches to ensure highly effective flushing while conserving reservoir water are presented （at least for the Sefid-Roud dam）. The operation of the bottom outlet gates, the shape of the output hydrograph, and the reservoir water level variation during flushing were optimized. In addition, the numerical model and related parameters, which need to be calibrated, are discussed.     

Numerical simulation of sedimentation processes in reservoirs as a function of outlet location

This paper builds on a recently published one-dimensional moving-boundary model of the coevolution of topset, foreset and bottomset in a reservoir that captures the dynamics of the internal muddy pond typical to reservoirs. This model was modified to account for different outlet locations at the reservoir＇s downstream end. This model considers a river carrying two sustained phases of sediments： coarse （sand） and fine （mud）. The coarse phase deposits in the topset and delta foreset, while the fine phase forms a dilute suspension of wash load in the river. As the river enters the reservoir, the muddy water plunges on the foreset to form a Froude-supercritical （purely depositional） turbidity current. This turbidity current emplaces the bottomset. The modified numerical model was tested against five laboratory experiments previously reported by the author. The model successfully locates the muddy-water/clear-water interface. In addition, modeled and measured bed deposits are in good agreement. Results clearly indicate that the location of the internal hydraulic jump plays a key role in the final bed deposit.     

Analytical study of the transverse distribution of along-channel and transverse residual flows in tidal estuaries

We present an analytical model to decompose complex along-channel and transverse residual flows into components induced by individual mechanisms. The model describes the transverse distribution of residual flows in tidally dominated estuaries. Scaling and perturbation techniques are used to obtain analytical solutions for residual flows over arbitrary across-channel bed profiles. The flows are induced by horizontal density gradients, tidal rectification processes, river discharge, wind, channel curvature and the earth's rotation. These rectification processes induce residual flows that are up-estuary to the right and down-estuary to the left of an estuarine channel (looking up-estuary in the northern hemisphere). The tidal rectification processes fundamentally change the transverse structure of along-channel residual flows in many tidal estuaries, as these processes cause the flows to be internally asymmetric about the mid-axis of the channel for relatively large tidal velocities, steep channels or narrow estuaries. In addition, velocity scales are derived from the analytical solutions to estimate the relative importance of the various residual flow mechanisms from estuarine parameters. A case study of a transect across the Upper Chesapeake Bay showed that important features of the residual flow observed in that transect are reproduced and explained by the analytical model. The velocity scales were able to identify the relevant residual flow mechanisms as well. The tidal rectification processes considered here result from advection of along-channel tidal momentum by Coriolis-induced transverse tidal currents.     

Coastal circulation and hydrography in the Gulf of Tehuantepec,Mexico, during winter

Winter observations of shelf and slope hydrography and currents in the inner Gulf of Tehuantepec are analysed from two field studies in 1989 and 1996 to specify the variability of near-shore conditions under varying wind stress. During the winter period frequent outbursts of ‘Norte’ winds over the central Gulf result in persistent alongshore inflows along both its eastern and western coasts. Wind-induced variability on time scales of several days strongly influences the shelf currents, but has greater effect on its western coast because of the generation and separation of anticyclonic eddies there. The steadier inflow (∼0.2 m s⁻¹) on the eastern shelf is evident in a strong down-bowing of shallow isosurfaces towards the coast within 100 km of shore, below a wedge of warmer, fresher and lighter water. This persistent entry of less saline (33.4–34.0), warmer water from the southeast clearly originates in buoyancy input by rivers along the Central American coast, but is augmented by a general shoreward tendency (0.2 m s⁻¹) in the southeastern Gulf. The resultant shallow tongue of anomalous water is generally swept offshore in the head of the Gulf and mixed away by the strong outflow and vertical overturning of the frequent ‘Norte’ events but during wind relaxations the warm, low-salinity coastal flow may briefly extend further west. In the head of the Gulf, flow is predominantly offshore (<0.2 m s⁻¹) as the alongshore component alternates eastward and westward in association with elevation or depression, respectively, of the pycnocline against the shore. More saline, open ocean water is introduced from the north-western side of the Gulf by the inflow along the west coast. During extended wind relaxations, the flow becomes predominantly eastward beyond the shelf while nearshore the coastally trapped buoyant inflow from the southeast penetrates across the entire head of the gulf at least as far as its western limit. On the basis of these and other recent observations, it seems that the accepted view of a broad, persistent Costa Rica Coastal Current (CRCC) is the result of averaging over many relatively sparse observations and that the instantaneous CRCC is a highly variable and convoluted flow around and between constantly changing eddies. The buoyancy-driven shelf current reported here forms a hitherto unrecognized, but major, component of this CRCC system.     

SST observations of upwelling induced by tidal straining in the Rhine ROFI

In this paper SST imagery and a three-dimensional numerical model of a river plume were employed to detect upwelling induced by tidal straining in the Rhine ROFI (region of fresh water influence). Previous studies have shown that the Rhine ROFI in the North Sea exhibits strong cross-shore density gradients that compete with tidal and wind mixing to establish stratification. During neap periods with low mixing energy an area measuring 30 km offshore by 100 km alongshore becomes stratified. When the ROFI is stratified strong cross-shore currents are observed, with surface currents rotating anti-cyclonically and bottom currents rotating cyclonically. The cross-shore currents interact with the cross-shore density gradients to produce a semi-diurnal cycle of stratification. Due to continuity requirements imposed by the proximity of the coast, the offshore-directed surface currents and onshore-directed bottom currents should lead to coastal upwelling.     

Features and formation processes of multiple deposition layers from the 2004 Indian Ocean Tsunami at Ban Nam Kem,southern Thailand

Multiple‐layered tsunami deposits have been frequently reported from coastal stratigraphic sequences, but the formation processes of these layers remain uncertain. A terrestrial sandy deposit formed by the 2004 Indian Ocean Tsunami was investigated at Ban Nam Kem, southern Thailand. Four internal layers induced by two tsunami waves were identified in the tsunami deposit. Sedimentary structures indicated that two units were formed by run‐up currents caused by the tsunami and the other two units were deposited by the backwash flows. Graded bedding was common in the layers, but inverse grading was observed at limited intervals on the surveyed transects. The characteristics of the multiple‐layered tsunami deposit vary remarkably over a very short distance (<1 m) in response to the local topography. Remarkable asymmetries in thickness and grain‐size distribution are recognized between the run‐up and backwash flow deposits. On the basis of the interpretation of sedimentary structures, the formation process of the multiple‐layered tsunami deposit observed in this study can be explained in a schematic model as the modification of the ideal tsunami sequence by local erosion and the asymmetric hydraulic properties of tsunami waves, such as the maximum shear velocity and the heterogeneity of the flow velocity field.