Evolution of a carbonate delta generated by gateway‐funnelling of episodic currents

Cool‐water carbonate sedimentation has dominated Mediterranean shelves since the Early Pliocene. Skeletal sand and gravel herein consist of remains of heterozoan organisms, which are susceptible to reworking due to weak early cementation in non‐tropical waters. This study documents the Lower Pleistocene carbonate wedge of Favignana Island (Italy), which prograded from a 5   km wide passage between two palaeo‐islands into a perpendicular, 10 to 15   km wide strait between the palaeo‐islands at one side and Sicily at the other during the Emilian highstand (1·6   Ma to 1·1   Ma). The clinoformed carbonate wedge, which is 50   m thick and 6   km long, formed by east/south‐east progradation of a platform on the submarine sill by currents that were funnelled between the two palaeo‐islands. Platform‐slope clinoforms evolved from initial aggradation (thin and low‐angle) into a progradation phase (thick and high‐angle). Both clinoform types are characterized by a bimodal facies stacking pattern defined by sedimentary structures created by: (i) subaqueous dunes associated with dilute subcritical currents; and (ii) upper‐flow‐regime bedforms associated with sediment‐laden supercritical turbidity currents. Focusing of episodic currents on the platform by funnelling between the islands controlled the downstream formation of a sediment body, here named carbonate delta. The carbonate delta interfingers with subaqueous dune deposits formed in the perpendicular strait. This study uses a reconstruction of bedform dynamics to unravel the evolution of this gateway‐related carbonate accumulation.     

An experimental investigation of sand–mud suspension settling behaviour: implications for bimodal mud contents of submarine flow deposits

The settling behaviour of particulate suspensions and their deposits has been documented using a series of settling tube experiments. Suspensions comprised saline solution and noncohesive glass‐ballotini sand of particle size 35·5 μm < d < 250 μm and volume fractions, φ_s, up to 0·6 and cohesive kaolinite clay of particle size d < 35·5 μm and volume fractions, φ_m, up to 0·15. Five texturally distinct deposits were found, associated with different settling regimes: (I) clean, graded sand beds produced by incremental deposition under unhindered or hindered settling conditions; (II) partially graded, clean sand beds with an ungraded base and a graded top, produced by incremental deposition under hindered settling conditions; (III) graded muddy sands produced by compaction with significant particle sorting by elutriation; (IV) ungraded clean sand produced by compaction and (V) ungraded muddy sand produced by compaction. A transition from particle size segregation (regime I) to suppressed size segregation (regime II or III) to virtually no size segregation (IV or V) occurred as sediment concentration was increased. In noncohesive particulate suspensions, segregation was initially suppressed at φ_s ～ 0·2 and entirely inhibited at φ_s ≥ 0·6. In noncohesive and cohesive mixtures with low sand concentrations (φ_s < 0·2), particle segregation was initially suppressed at φ_m ～ 0·07 and entirely suppressed at φ_m ≥ 0·13. The experimental results have a number of implications for the depositional dynamics of submarine sediment gravity flows and other particulate flows that carry sand and mud; because the influence of moving flow is ignored in these experiments, the results will only be applicable to flows in which settling processes, in the depositional boundary, dominate over shear‐flow processes, as might be the case for rapidly decelerating currents with high suspended load fallout rates. The ‘abrupt’ change in settling regimes between regime I and V, over a relatively small change in mud concentration (<5% by volume), favours the development of either mud‐poor, graded sandy deposits or mud‐rich, ungraded sandy deposits. This may explain the bimodality in sediment texture (clean ‘turbidite’ or muddy ‘debrite’ sand or sandstone) found in some turbidite systems. Furthermore, it supports the notion that distal ‘linked’ debrites could form because of a relatively small increase in the mud concentration of turbidity currents, perhaps associated with erosion of a muddy sea floor. Ungraded, clean sand deposits were formed by noncohesive suspensions with concentrations 0·2 ≤ φ_s ≤ 0·4. Hydrodynamic sorting is interpreted as being suppressed in this case by relatively high bed aggradation rates which could also occur in association with sustained, stratified turbidity currents or noncohesive debris flows with relatively high near‐bed sediment concentrations.     

Energy relations of density-current flows: an experimental investigation*

Three series of density-current experiments were performed in a 5.76 m flume. In the first series, the flume was horizontal, and in the second and third, it was inclined with a positive slope and negative slope, respectively. Energy relations during successive stages of density-current movement were computed from observed data, which showed an appreciable frictional energy dissipation. The computed friction factors of our experimental density-flows were compared to the friction factors for pipe flows (Moody diagram), and while the calculated friction factor increases with increasing Reynold's number within the range of our experiments (Re 2 × 10³?2 × 10⁴), it is concluded that with increasing Reynold's number above about 5 × 10⁴ the friction factor decreases. For natural turbidity currents, the Moody diagram gives a reasonable estimate of the friction factor between the current and sediment bed. The value of the friction factor for the interface between the current and overlying water was found to be about 0.2 times the friction factor for the current and flume. However, due to errors inherent in measuring the depth of the current, a value of 0.4 would be more reasonable for density-currents in our range of Reynold's number. Friction tends to decrease the value of the dimensionless coefficient in Keulegan's law of saline front and to decrease the thickness of the flow. In contrast, the presence of a slope in the direction of flow tends to compensate the effect of friction. The angle θ^c that provides the potential energy to exactly offset the energy losses incurred during movement by the density-currents in our experiments has a calculated value of 31′. An empirical formula φ= 0.935θ—0·57 relating friction, in terms of the hydraulic gradient φ, to the slope angle θ was obtained. Since the thickness of the current can be computed from the relationship between φ and θ, we estimated the thickness of naturally occurring density-currents in Swiss lakes. The results suggest the applicability of our experimental results to small turbidity currents in nature. Our analysis further indicates that large turbidity currents have a small φ and can be expected to flow very long distances on a flat abyssal plain.     

Reverse flow in turbidity currents: the role of internal solitons

Pickering & Hiscott, (1985) have demonstrated amply the presence of reverse-flow units within the thick-bedded calcareous wacke (TCW) beds of the turbiditic Cloridorme Formation (Middle Ordovician, Gaspé Peninsula, Quebec, Canada). These reverse-flow units are underlain and overlain by units which reveal flow in the primary (obverse) direction. In this paper, a model is proposed for this reverse flow, based on the probable nature of the primary turbidity flow. It appears that the initial flow was highly elongated (thickness h? length L), with h～ 500 m, velocity U～ 2 m s^-1 and sediment concentration C～ 1·25%_o. The rate of momentum loss of the flow is estimated by means of a useful parameter which we call the ‘drag distance’, symbol d_D, defined by where h and L are the thickness and length of the flow, respectively; _cC_d is a combined drag coefficient representing friction on the bottom and at the upper interface; and _fC_d is a form-drag coefficient related to the shape and size of the head. d_D is the distance travelled by a current of constant h and L, flowing over a horizontal bottom and obeying a quadratic friction law, for an e-fold reduction in velocity. Simple considerations, confirmed by our own experiments (described in this paper), show that such an elongated turbidity current cannot be reflected as a whole from an adverse slope: when the nose of the current reaches the slope, it forms a hump, which surges backwards and sooner or later breaks up into a series of internal solitons. The latter, probably numbering 4–7, will cause reverse flow at a given point as they pass by, provided that the residual velocity in the tail is not too great. Flow in the original (obverse) direction will be re-established after the passage of the solitons. Quiescent periods in front of, between and behind the solitons, when soliton-associated currents cancelled out the residual obverse flow, would allow the deposition of thin mud-drapes. Additional flow reversals observed in a few of the TCW beds cannot be explained readily by the re-passage of solitons, since wave breaking at the ends of the basin would cause massive energy loss; internal seiches are the preferred explanation for these later reversals.     

The formation of large mudwaves by turbidity currents on the levees of the Toyama deep-sea channel, Japan Sea

The development of mudwaves on the levees of the modern Toyama deep‐sea channel has been studied using gravity core samples combined with 3·5‐kHz echosounder data and airgun seismic reflection profiles. The mudwaves have developed on the overbank flanks of a clockwise bend of the channel in the Yamato Basin, Japan Sea, and the mudwave field covers an area of 4000 km². Mudwave lengths range from 0·2 to 3·6 km and heights vary from 2 to 44 m, and the pattern of mudwave aggradation indicates an upslope migration direction. Sediment cores show that the mudwaves consist of an alternation of fine‐grained turbidites and hemipelagites whereas contourites are absent. Core samples demonstrate that the sedimentation rate ranged from 10 to 14 cm ka^?1 on the lee sides to 17–40 cm ka^?1 on the stoss sides. A layer‐by‐layer correlation of the deposits across the mudwaves shows that the individual turbidite beds are up to 20 times thicker on the stoss side than on the lee side, whereas hemipelagite thicknesses are uniform. This differential accretion of turbidites is thought to have resulted in the pattern of upcurrent climbing mudwave crests, which supports the notion that the mudwaves have been formed by spillover turbidity currents. The mudwaves are interpreted to have been instigated by pre‐existing large sand dunes that are up to 30 m thick and were created by high‐velocity (10°ms^?1), thick (c. 500 m) turbidity currents spilling over the channel banks at the time of the maximum uplift of the Northern Japan Alps during the latest Pliocene to Early Pleistocene. Draping of the dunes by the subsequent, lower‐velocity (10^?1ms^?1), mud‐laden turbidity currents is thought to have resulted in the formation of the accretionary mudwaves and the pattern of upflow climbing. The dune stoss slopes are argued to have acted as obstacles to the flow, causing localized loss of flow strength and leading to differential draping by the muddy turbidites, with greater accretion occurring on the stoss side than on the lee slope. The two overbank flanks of the clockwise channel bend show some interesting differences in mudwave development. The mudwaves have a mean height of 9·8 m on the outer‐bank levee and 6·2 m on the inner bank. The turbidites accreted on the stoss sides of the mudwaves are 4–6 times thicker on the outer‐bank levee than their counterparts on the inner‐bank levee. These differences are attributed to the greater flow volume (thickness) and sediment flux of the outer‐bank spillover flow due to the more intense stripping of the turbidity currents at the outer bank of the channel bend. Differential development of mudwave fields may therefore be a useful indicator in the reconstruction of deep‐sea channels and their flow hydraulics.     

Subaqueous sediment density flows: Depositional processes and deposit types

Submarine sediment density flows are one of the most important processes for moving sediment across our planet, yet they are extremely difficult to monitor directly. The speed of long run‐out submarine density flows has been measured directly in just five locations worldwide and their sediment concentration has never been measured directly. The only record of most density flows is their sediment deposit. This article summarizes the processes by which density flows deposit sediment and proposes a new single classification for the resulting types of deposit. Colloidal properties of fine cohesive mud ensure that mud deposition is complex, and large volumes of mud can sometimes pond or drain‐back for long distances into basinal lows. Deposition of ungraded mud (T_E‐3) most probably finally results from en masse consolidation in relatively thin and dense flows, although initial size sorting of mud indicates earlier stages of dilute and expanded flow. Graded mud (T_E‐2) and finely laminated mud (T_E‐1) most probably result from floc settling at lower mud concentrations. Grain‐size breaks beneath mud intervals are commonplace, and record bypass of intermediate grain sizes due to colloidal mud behaviour. Planar‐laminated (T_D) and ripple cross‐laminated (T_C) non‐cohesive silt or fine sand is deposited by dilute flow, and the external deposit shape is consistent with previous models of spatial decelerating (dissipative) dilute flow. A grain‐size break beneath the ripple cross‐laminated (T_C) interval is common, and records a period of sediment reworking (sometimes into dunes) or bypass. Finely planar‐laminated sand can be deposited by low‐amplitude bed waves in dilute flow (T_B‐1), but it is most likely to be deposited mainly by high‐concentration near‐bed layers beneath high‐density flows (T_B‐2). More widely spaced planar lamination (T_B‐3) occurs beneath massive clean sand (T_A), and is also formed by high‐density turbidity currents. High‐density turbidite deposits (T_A, T_B‐2 and T_B‐3) have a tabular shape consistent with hindered settling, and are typically overlain by a more extensive drape of low‐density turbidite (T_D and T_C,). This core and drape shape suggests that events sometimes comprise two distinct flow components. Massive clean sand is less commonly deposited en masse by liquefied debris flow (D_CS), in which case the clean sand is ungraded or has a patchy grain‐size texture. Clean‐sand debrites can extend for several tens of kilometres before pinching out abruptly. Up‐current transitions suggest that clean‐sand debris flows sometimes form via transformation from high‐density turbidity currents. Cohesive debris flows can deposit three types of ungraded muddy sand that may contain clasts. Thick cohesive debrites tend to occur in more proximal settings and extend from an initial slope failure. Thinner and highly mobile low‐strength cohesive debris flows produce extensive deposits restricted to distal areas. These low‐strength debris flows may contain clasts and travel long distances (D_M‐2), or result from more local flow transformation due to turbulence damping by cohesive mud (D_M‐1). Mapping of individual flow deposits (beds) emphasizes how a single event can contain several flow types, with transformations between flow types. Flow transformation may be from dilute to dense flow, as well as from dense to dilute flow. Flow state, deposit type and flow transformation are strongly dependent on the volume fraction of cohesive fine mud within a flow. Recent field observations show significant deviations from previous widely cited models, and many hypotheses linking flow type to deposit type are poorly tested. There is much still to learn about these remarkable flows.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.     

Some new aspects of two-dimensional turbidity currents*

A theoretical consideration of two dimensional underflows and surge-type turbidity currents results in a general momentum equation. A number of formulae in current use are special cases of this equation, among which are the modified Chézy equation and Bagnold's criterion for autosuspension. Five dimensionless parameters are included: the Richardson number Ri (defined as the inverse square of the Froude number), the friction coefficient c_f, the slope β, the dimensionless settling velocity of the sediment V_s/u and the changes in flow height with distance dD/dx. The latter is mainly a measure of the dilution by entrainment of ambient water. For chalk powder experiments on surge type turbidity currents and on the initial front of continuous underflows the momentum equation is shown to be correct. Values for Ri range from about 1.5 at 0° slope to about 0.75 at 5° and are slightly to substantially lower than values from earlier authors. The two types of turbidity currents investigated show close similarity. A surprising attribute is their strong dilution even at very low-angle slopes. Pelitic sedimentation is possible from the upper, dilute part of the currents, graded intervals found at the base of turbidites can be explained as bedload deposits from the lowermost, concentrated layer of the current; hydraulic jumps are expected to be rare in surge-type turbidity currents and fronts of incipient underflows.     

Bedload transport and bed resistance associated with density and turbidity currents

Turbidity currents in the ocean are driven by suspended sediment. Yet results from surveys of the modern sea floor and turbidite outcrops indicate that they are capable of transporting as bedload and depositing particles as coarse as cobble sizes. While bedload cannot drive turbidity currents, it can strongly influence the nature of the deposits they emplace. This paper reports on the first set of experiments which focus on bedload transport of granular material by density underflows. These underflows include saline density flows, hybrid saline/turbidity currents and a pure turbidity current. The use of dissolved salt is a surrogate for suspended mud which is so fine that it does not settle out readily. Thus, all the currents can be considered to be model turbidity currents. The data cover four bed conditions: plane bed, dunes, upstream‐migrating antidunes and downstream‐migrating antidunes. The bedload transport relation obtained from the data is very similar to those obtained for open‐channel flows and, in fact, is fitted well by an existing relation determined for open‐channel flows. In the case of dunes and downstream‐migrating antidunes, for which flow separation on the lee sides was observed, form drag falls in a range that is similar to that due to dunes in sand‐bed rivers. This form drag can be removed from the total bed shear stress using an existing relation developed for rivers. Once this form drag is subtracted, the bedload data for these cases collapse to follow the same relation as for plane beds and upstream‐migrating antidunes, for which no flow separation was observed. A relation for flow resistance developed for open‐channel flows agrees well with the data when adapted to density underflows. Comparison of the data with a regime diagram for field‐scale sand‐bed rivers at bankfull flow and field‐scale measurements of turbidity currents at Monterey Submarine Canyon, together with Shields number and densimetric Froude number similarity analyses, provide strong evidence that the experimental relations apply at field scale as well.     

A comparison between Shields' threshold criterion and the movement of loosely packed gravel in a tidal channel

Observations of the threshold of movement of loosely packed gravel in a tidal current are described. For gravel with equivalent ‘spherical’ diameters D in the range 0.2 ?D? 5.0cm the critical friction velocity u*_c, corresponding to the initiation of sediment transport, is given by u*_c=7.0 D^0.2. At large values of D within the quoted range, the value u*_c is significantly lower than would be obtained by a Shields experiment (u*_c∞D^0.5). By comparing our values of u*_c with those obtained under well-controlled laboratory conditions, the discrepancy with Shields is shown to be due to the open spacing between, and exposure of, individual pebbles on the seabed. By comparing our results with those from upland gravel streams and flume experiments, it is suggested that Shields assumed an excessively large water depth to particle size ratio as a constraint within which the critical sediment entrainment number 0_c is valid.     

The stability of sublittoral, fine-grained sediments in a subarctic estuary

The erodibility of natural estuarine sediments was measured in sit along a longitudinal transect of Manitounuk Sound, Hudson Bay, using the benthic flume Sea Carousel. Sedimentation processes along the transect varied from continuous, rapid, post-glacial sedimentation in the inner Sound, to glacial outcrops and seabed reworking of the outer Sound. The grain size and physical bulk properties reflect changes in depositional environment and correlate with sediment erosion threshold stress (τ_c), erosion rate (E), erosion type and still-water mass settling rate. There was a steady increase in τ_c (0·8–2·0 Pa) with distance down the Sound in parallel with the decreasing sedimentation rate (0·003–0·001 m yr^?1) and increasing sediment bulk density (1650–2010 kg m^?3). The near-surface friction coefficient varied up to 68° in proportion to the clay content of post-glacial material. Glacial sediments were characterized by variable results and generally higher friction coefficients. Seabed erosion in Sea Carousel began with surface creep of loose aggregates, pellets and organic debris. This was followed by Type I bed erosion at rates that varied between 0·0002 and 0·0032 kg m^?2 s^?1 (mean 0·0015). Type I peak erosion rate was inversely related to applied bed shear stress (τ_o). Type II erosion succeeded Type I, often after a broad transitional period. Simulations of suspended sediment concentration in Sea Carousel were made using four commonly used erosion (E) algorithms. The best results were obtained using Krone's dimensionless ratio relationship: E=M(τ_o/τ_c-1). Simulations were highly sensitive to the definition of erosion threshold with sediment depth [τ_c(z)]. Small errors in definition of τ_c(z) caused large errors in the prediction of suspended sediment concentration which far exceeded differences between the methods tested.     

Three‐dimensional seismic analysis of sediment waves and related geomorphological features on a carbonate shelf exposed to large amplitude internal waves,Browse Basin region,Australia

This study analyses the three‐dimensional geometry of sedimentary features recorded on the modern sea floor and in the shallow subsurface of a shelf to upper slope region offshore Australia that is characterized by a pronounced internal wave regime. The data interpreted comprise an extensive, >12 500 km² industrial three‐dimensional seismic‐reflection survey that images the northern part of the Browse Basin, Australian North West Shelf. The most prominent seismic–morphological features on the modern sea floor are submarine terrace escarpments, fault‐scarps and incised channels, as well as restricted areas of seismic distortion interpreted as mass wasting deposits. Besides these kilometre‐scale sea floor irregularities, smaller bedforms were discovered also, including a multitude of sediment waves with a lateral extent of several kilometres and heights up to 10 m. These sedimentological features generally occur in extensive fields in water depths below 250 m mostly at the foot of submerged terraces, along the scarps of modern faults and along the shelf break between the outer shelf and the upper continental rise. Additional bedforms that characterize the more planar regions of the outer shelf are elongate, north‐west/south‐east oriented furrows and ridges. The formation of both sediment waves and furrow‐ridge systems requires flow velocities between 0·3 m sec^?1 and 1·5 m sec^?1, which could be generated by oceanic currents, gravity currents or internal waves. In the studied setting, these velocities can be best explained as being generated by bottom currents induced by internal waves, an interpretation that is discussed against oceanographic background data and modelling results. In addition to the documentation of three‐dimensional seismic–geomorphological features of the modern sea floor, it was also possible to map kilometre‐scale buried sediment wave fields in the seismic volume down to ca 500 ms two‐way‐time below the present sea floor, indicating the general potential for the preservation of such bedforms in the sedimentary record.     

An empirical model for the development and equilibrium morphology of current ripples in fine sand

A flume study on the development and equilibrium morphology of current ripples in fine sand (D₅₀ = 0·238 mm) was performed to extend an empirical model for current ripple stability in 0·095 mm sand to larger grain sizes. The results of the flume experiments agree with the very fine sand model that current ripple development from a flat bed is largely independent of flow velocity. At all flow velocities, ripples evolve from incipient, through straight, sinuous and non-equilibrium linguoid, to equilibrium linguoid plan morphology. The time needed to achieve an equilibrium linguoid plan form is related to an inverse power of flow velocity and ranges from several minutes to more than hundreds of hours. Average equilibrium height and length are 17·0 mm and 141·1 mm respectively. These values are about 20% larger than in very fine sand. Equilibrium ripple height and length are proportional to flow velocity near the stability field of dunes. In the same velocity range, a characteristic grouping of ripples with smaller ripples migrating on the upstream face of larger ripples was observed. Bed-form development shows a conspicuous two-phase behaviour at flow velocities < 0·49 m s^?1. In the first phase of development, ripple height and length increase along an exponential path, similar to that at higher flow velocities, thus reaching intermediate equilibrium values of 14·8 mm and 124·5 mm respectively. After some time, however, a second phase commences, that involves a rapid increase in bed-form size to the typical equilibrium values for 0·238 mm sand. A comparison with literature data shows that the results obtained for 0·238 mm sand agree reasonably well with other flume studies at similar grain size. Yet considerable variability in the relationships between ripple dimensions and flow strength ensues from, among others, underestimation of equilibrium time, shallow flow depths and differences in sediment texture.     

Sediment waves in a modern high-energy glacilacustrine environment

A 4·7 km² field of sediment waves occurs in front of the Slims River delta in Kluane Lake, the largest lake in the Yukon Territory. Slims River heads in the Kaskawulsh Glacier, part of the St Elias Ice Field and discharges up to 400 m³ s^?1 of water with suspended sediment concentrations of up to 7 g l^?1. The 19 km long sandur of Slims River was created in the past 400 years since Kaskawulsh Glacier advanced and dammed the lake and the sandur has advanced into Kluane Lake at an average rate of 48 m a^?1. However, this rate is decreasing as flow is diverted from Slims River because of the retreat of the Kaskawulsh Glacier. The sandur and a road constructed on the delta remove coarse‐grained sediment, so the river delivers dominantly mud to the lake. Inflow during summer generates quasi‐continuous turbidity currents with velocities up to 0·6 m s^?1. The front of the delta consists of a plane surface sloping lakeward at 0·0188 (1·08°). A field of sediment waves averaging 130 m in length and 2·3 m in amplitude has developed on this surface. Slopes on the waves vary from ?0·067 (?3·83°, i.e. sloping in the opposite direction to the regional slope) to 0·135 (7·69°). The internal structure of the sediment waves, as documented by seismic profiling, shows that sedimentation on the stoss portion of the wave averages 2·7 times that on the lee portion. Rates of sediment accumulation in the wave field are about 0·3 m a^?1, so these lacustrine waves have formed in a much shorter period of time (less than 200 years) and are advancing upslope towards the delta much more quickly (1 to 2 m a^?1) than typical marine sediment waves. These waves formed on the flat surface of the lake floor, apparently in the absence of pre‐existing forms, and they are altered and destroyed as the wave field advances and the characteristics of the turbidity currents change.     

Time‐scale of deformation and intertectonic phases revealed by P–T–D–t relationships in the orogenic middle crust of the Orlica‐Śnieżnik Dome,Polish/Czech Central Sudetes

A section of the orogenic middle crust (Orlica‐?nie?nik Dome, Polish/Czech Central Sudetes) was examined to constrain the duration and significance of deformation (D) and intertectonic (I) phases. In the studied metasedimentary synform, three deformation events produced an initial subhorizontal foliation S1 (D₁), a subsequent subvertical foliation S2 (D₂) and a late subhorizontal axial planar cleavage S3 (D₃). The synform was intruded by pre‐, syn‐ and post‐D₂ granitoid sheets. Crystallization–deformation relationships in mica schist samples document I_1–2 garnet–staurolite growth, syn‐D₂ staurolite breakdown to garnet–biotite–sillimanite/andalusite, I_2–3 cordierite blastesis and late‐D₃ chlorite growth. Garnet porphyroblasts show a linear Mn–Ca decrease from the core to the inner rim, a zone of alternating Ca–Y‐ and P‐rich annuli in the inner rim, and a Ca‐poor outer rim. The Ca–Y‐rich annuli probably reflect the occurrence of the allanite‐to‐monazite transition at conditions of the staurolite isograd, whereas the Ca‐poor outer rim is ascribed to staurolite demise. The reconstructed P–T path, obtained by modelling the stability of parageneses and garnet zoning, documents near‐isobaric heating from ~4 kbar/485 °C to ~4.75 kbar/575 °C during I_1–2. This was followed by a progression to 4–5 kbar/580–625 °C and a subsequent pressure decrease to 3–4 kbar during D₂. Pressure decrease below 3 kbar is ascribed to I_2–3, whereas cooling below ~500 °C occurred during D₃. In the dated mica schist sample, garnet rims show strong Lu enrichment, oscillatory Lu zoning and a slight Ca increase. These features are also related to allanite breakdown coeval with staurolite appearance. As Lu‐rich garnet rims dominate the Lu–Hf budget, the 344 ± 3 Ma isochron age is ascribed to garnet crystallization at staurolite grade, near the end of I_1–2. For the dated sample of amphibole–biotite granitoid sheet, a Pb–Pb single zircon evaporation age of 353 ± 1 Ma is related to the onset of plutonic activity. The results suggest a possible Devonian age for D₁, and a Carboniferous burial‐exhumation cycle in mid‐crustal rocks that is broadly coeval with the exhumation of neighbouring HP rocks during D₂. In the light of published ages, a succession of telescoping stages with time spans decreasing from c. 10 to 2–3 Ma is proposed. The initially long period of tectonic quiescence (I_1–2 phase, c. 10 Ma) inferred in the middle crust contrasts with contemporaneous deformation at deeper levels and points to decoupled P–T–D histories within the orogenic wedge. An elevated gradient of ~30 °C km^?1 and assumed high heating rates of c. 20 °C Ma^?1 are explained by the protracted intrusion of granitoid sheets, with or without deformation, whereas fast vertical movements (2–3 Ma, D₂ phase) in the crust require the activity of deformation phases.     

Response of unconfined turbidity current to deep‐water fold and thrust belt topography: Orthogonal incidence on solitary and segmented folds

Sea‐floor topography of deep‐water folds is widely considered to have a major impact on turbidity currents and their depositional systems, but understanding the flow response to such features was limited mainly to conceptual notions inspired by small‐scale laboratory experiments. High‐resolution three‐dimensional numerical experiments can compensate for the lack of natural‐scale flow observations. The present study combines numerical modelling of thrusts with fault‐propagation folds by Trishear3D software with computational fluid dynamics simulations of a natural‐scale unconfined turbidity current by MassFlow‐3D? software. The study reveals the hydraulic and depositional responses of a turbidity current (ca 50 m thick) to typical topographic features that it might encounter in an orthogonal incidence on a sea‐floor deep‐water fold and thrust belt. The supercritical current (ca 10 m sec^?1) decelerated and thickened due to the hydraulic jump on the fold backlimb counter‐slope, where a reverse overflow formed through current self‐reflection and a reverse underflow was issued by backward squeezing of a dense near‐bed sediment load. The reverse flows were re‐feeding sediment to the parental current, reducing its waning rate and extending its runout. The low‐efficiency current, carrying sand and silt, outran a downslope distance of >17 km with only modest deposition (<0·2 m) beyond the fold. Most of the flow volume diverted sideways along the backlimb to surround the fold and spread further downslope, with some overspill across the fold and another hydraulic jump at the forelimb toe. In the case of a segmented fold, a large part of the flow went downslope through the segment boundary. Preferential deposition (0·2 to 1·8 m) occurred on the fold backlimb and directly upslope, and on the forelimb slope in the case of a smaller fold. The spatial patterns of sand entrapment revealed by the study may serve as guidelines for assessing the influence of substrate folds on turbiditic sedimentation in a basin.     

The stability of a remediated bed in Hamilton Harbour,Lake Ontario,Canada

ABSTRACT In situ measurements of lakebed sediment erodibility were made on three sites in Hamilton Harbour, Lake Ontario, using the benthic flume Sea Carousel. Three methods of estimating the surface erosion threshold (τ_c(0)) from a Carousel time series were evaluated: the first method fits measures of bed strength to eroded depth (the failure envelope) and evaluates threshold as the surface intercept; the second method regresses mean erosion rate (E_m) with bed shear stress and solves for the floc erosion rate (E_f) to derive the threshold for E_m = E_f = 1 × 10^?5 kg m^?2 s^?1; the third method extrapolates a regression of suspended sediment concentration (S) and fluid transmitted bed shear stress (τ₀) to ambient concentrations. The first field site was undisturbed (C) and acted as a control; the second (W) was disturbed through ploughing and water injection as part of lakebed treatment, whereas the third site (OIP) was disturbed and injected with an oxidant used for remediation of contaminated sediment. The main objectives of this study were: (1) to evaluate the three different methods of deriving erosion threshold; (2) to compare the physical behaviour of lacustrine sediments with their marine estuarine counterparts; and (3) to examine the effects of ploughing and chemical treatment of contaminated sediment on bed stability. Five deployments of Sea Carousel were carried out at the control site. Mean erosion thresholds for the three methods were: τ_c(0) = 0·5 (±0·06), 0·27 (±0·01) and 0·34 (±0·03) Pa respectively. Method 1 overpredicted bed strength as it was insensitive to effects in the surface 1–2 mm, and the fit of the failure envelope was also highly subjective. Method 2 exhibited a wide scatter in the data (low correlation coefficients), and definition of the baseline erosion rate (E_f) is largely arbitrary in the literature. Method 3 yielded stable (high correlation coefficients), reproducible and objective results and is thus recommended for evaluation of the erosion threshold. The results of this method correlated well with sediment bulk density and followed the same trend as marine counterparts from widely varying sites. Mass settling rates, expressed as a decay constant, k, of S(t), were strongly related to the maximum turbidity at the onset of settling (S_max) and were also in continuity with marine counterparts. Thus, it appears that differences in salinity had little effect on mass settling rates in the examples presented, and that biological activity dominated any effects normally attributable to changes in salinity. Bedload transport of eroded aggregates (2–4 mm in diameter) took place by rolling below a mean tangential flow velocity (U_y) of 0·32 ms^?1 and by saltation at higher velocities. Mass transport as bedload was a maximum at U_y = 0·4 ms^?1, although bedload never exceeded 1% of the suspended load. The proportion of material moving as bedload was greatest at the onset of erosion but decreased as flow competence increased. Given the low bulk density and strength of the lakebed sediment, the presence of a bedload component is notable. Bedload transport over eroding cohesive substrates should be greater in estuaries, where both sediment density and strength are usually higher. Significant differences between the ploughed and control sites were apparent in both the erosion rate and the friction coefficient (φ), and suggest that bed recovery after disruption is rapid (< 24 h). τ_c(0) increased linearly with time after ploughing and recovered to the control mean value within 3 days. The friction coefficient was reduced to zero by ploughing (diagnostic of fluidization), but increased linearly with time, regaining control values within 6 days. No long‐term reduction in bed strength due to remediation was apparent.     

The impact of mobile disarticulated shells ofCerastoderma edulis on the abrasion of a cohesive substrate

An annular laboratory flume was used to investigate the effect of mobile cockle shells on the erosion of a cohesive sediment bed. A standard clay bed was created and shells of differing sizes placed upon it. Flow in the flume was increased in increments and the onset of motion and the transport patterns of the cockles were monitored. The release of bed material to the water column was monitored and compared to controls made in the absence of shells (due only to the flow). The shells moved as bedload; first as surface creep (sliding) and then by rolling. The onset velocity of motion (U_c) of the shells was found to be directly related to the settling rate (W_s) in still water. The fluid-induced stresses did not cause any detectable erosion of the bed. The addition of even a single shell induced significant erosion rates (E). The erosion was found to be the result of abrasion rather than corrasion, as the shells never entered into saltation. There was a linear increase in erosion rate with increasing shell size, and an exponential increase in the suspended sediment concentration with time. The drag coefficients (C_d) for settling in traction were calculated. The ratio of the drag forces acting on the shells when settling and moving as traction was found to equal to 1/tan(ф) where ф is the friction angle.     

Estimation of drag coefficient over the western desert sector of the Indian summer monsoon trough

In the estimation of momentum fluxes over land surfaces by the bulk aerodynamic method, no unique value of the drag coefficient (C _D) is found in the literature. The drag coefficient is generally estimated from special observations at different parts of the world. In this study an attempt is made to estimate drag coefficient over the western desert sector of India using data sets of Monsoon Trough Boundary Layer Experiment (MONTBLEX) during the summer monsoon season of 1990. For this purpose, the fast and slow response data sets obtained simultaneously from a 30 m high micro-meteorological tower at Jodhpur are used. All the observations used in this study are confined to a wind speed regime of 2.5–9.0 ms⁻¹. A comparison of momentum fluxes computed by eddy correlation (direct estimation) with profile and bulk aerodynamic (C _D = 3.9 × 10⁻³, Garratt, 1977) methods revealed that though the nature of variation of the fluxes by all these methods is almost similar, both the indirect methods give an under-estimated value of the fluxes. The drag coefficient is estimated as a function of wind speed and surface stability by a multiple regression approach. An average value of the estimated drag coefficient is found to be of the order of 5.43 × 10⁻³. The estimated value ofC _D is validated with a set of independent observations and found to be quite satisfactory. The recomputed momentum fluxes by bulk aerodynamic method using the estimated drag coefficient are in close agreement with the directly estimated fluxes.