Physical theory for near-bed turbulent particle suspension capacity.

The inability to capture the physics of solid-particle suspension in turbulent fluids in simple formulas is holding back the application of multiphase fluid dynamics techniques to many practical problems in nature and society involving particle suspension. We present a force balance approach to particle suspension in the region near no-slip frictional boundaries of turbulent flows. The force balance parameter Γ contains gravity and buoyancy acting on the sediment and vertical turbulent fluid forces; it includes universal turbulent flow scales and material properties of the fluid and particles only. Comparison to measurements shows that Γ = 1 gives the upper limit of observed suspended particle concentrations in a broad range of flume experiments and field settings. The condition of Γ > 1 coincides with the complete suppression of coherent turbulent structures near the boundary in direct numerical simulations of sediment-laden turbulent flow. Γ thus captures the maximum amount of sediment that can be contained in suspension at the base of turbulent flow, and it can be regarded as a suspension capacity parameter. It can be applied as a simple concentration boundary condition in modelling studies of the dispersion of particulates in environmental and man-made flows

Newly recognized turbidity current structure can explain prolonged flushing of submarine canyons

Seabed-hugging flows called turbidity currents are the volumetrically most important process transporting sediment across our planet and form its largest sediment accumulations. We seek to understand the internal structure and behavior of turbidity currents by reanalyzing the most detailed direct measurements yet of velocities and densities within oceanic turbidity currents, obtained from weeklong flows in the Congo Canyon. We provide a new model for turbidity current structure that can explain why these are far more prolonged than all previously monitored oceanic turbidity currents, which lasted for only hours or minutes at other locations. The observed Congo Canyon flows consist of a short-lived zone of fast and dense fluid at their front, which outruns the slower moving body of the flow. We propose that the sustained duration of these turbidity currents results from flow stretching and that this stretching is characteristic of mud-rich turbidity current systems. The lack of stretching in previously monitored flows is attributed to coarser sediment that settles out from the body more rapidly. These prolonged seafloor flows rival the discharge of the Congo River and carry ~2% of the terrestrial organic carbon buried globally in the oceans each year through a single submarine canyon. Thus, this new structure explains sustained flushing of globally important amounts of sediment, organic carbon, nutrients, and fresh water into the deep ocean

Electrical conductivity during incipient melting in the oceanic low-velocity zone

International audienceThe low-viscosity layer in the upper mantle, the asthenosphere, is a requirement for plate tectonics1. The seismic low velocities and the high electrical conductivities of the asthenosphere are attributed either to subsolidus, water-related defects in olivine minerals2, 3, 4 or to a few volume per cent of partial melt5, 6, 7, 8, but these two interpretations have two shortcomings. First, the amount of water stored in olivine is not expected to be higher than 50 parts per million owing to partitioning with other mantle phases9 (including pargasite amphibole at moderate temperatures10) and partial melting at high temperatures9. Second, elevated melt volume fractions are impeded by the temperatures prevailing in the asthenosphere, which are too low, and by the melt mobility, which is high and can lead to gravitational segregation11, 12. Here we determine the electrical conductivity of carbon-dioxide-rich and water-rich melts, typically produced at the onset of mantle melting. Electrical conductivity increases modestly with moderate amounts of water and carbon dioxide, but it increases drastically once the carbon dioxide content exceeds six weight per cent in the melt. Incipient melts, long-expected to prevail in the asthenosphere10, 13, 14, 15, can therefore produce high electrical conductivities there. Taking into account variable degrees of depletion of the mantle in water and carbon dioxide, and their effect on the petrology of incipient melting, we calculated conductivity profiles across the asthenosphere for various tectonic plate ages. Several electrical discontinuities are predicted and match geophysical observations in a consistent petrological and geochemical framework. In moderately aged plates (more than five million years old), incipient melts probably trigger both the seismic low velocities and the high electrical conductivities in the upper part of the asthenosphere, whereas in young plates4, where seamount volcanism occurs6, a higher degree of melting is expected

Monte Carlo Simulations of Metasomatic Enrichment in the Lithosphere and Implications for the Source of Alkaline Basalts

One hypothesis for the origin of alkaline lavas erupted on oceanic islands and in intracontinental settings is that they represent the melts of amphibole-rich veins in the lithosphere (or melts of their dehydrated equivalents if metasomatized lithosphere is recycled into the convecting mantle). Amphibole-rich veins are interpreted as cumulates produced by crystallization of low-degree melts of the underlying asthenosphere as they ascend through the lithosphere. We present the results of trace-element modelling of the formation and melting of veins formed in this way with the goal of testing this hypothesis and for predicting how variability in the formation and subsequent melting of such cumulates (and adjacent cryptically and modally metasomatized lithospheric peridotite) would be manifested in magmas generated by such a process. Because the high-pressure phase equilibria of hydrous near-solidus melts of garnet lherzolite are poorly constrained and given the likely high variability of the hypothesized accumulation and remelting processes, we used Monte Carlo techniques to estimate how uncertainties in the model parameters (e.g. the compositions of the asthenospheric sources, their trace-element contents, and their degree of melting; the modal proportions of crystallizing phases, including accessory phases, as the asthenospheric partial melts ascend and crystallize in the lithosphere; the amount of metasomatism of the peridotitic country rock; the degree of melting of the cumulates and the amount of melt derived from the metasomatized country rock) propagate through the process and manifest themselves as variability in the trace-element contents and radiogenic isotopic ratios of model vein compositions and erupted alkaline magma compositions. We then compare the results of the models with amphibole observed in lithospheric veins and with oceanic and continental alkaline magmas. While the trace-element patterns of the near-solidus peridotite melts, the initial anhydrous cumulate assemblage (clinopyroxene ± garnet ± olivine ± orthopyroxene), and the modelled coexisting liquids do not match the patterns observed in alkaline lavas, our calculations show that with further crystallization and the appearance of amphibole (and accessory minerals such as rutile, ilmenite, apatite, etc.) the calculated cumulate assemblages have trace-element patterns that closely match those observed in the veins and lavas. These calculated hydrous cumulate assemblages are highly enriched in incompatible trace elements and share many similarities with the trace-element patterns of alkaline basalts observed in oceanic or continental setting such as positive Nb/La, negative Ce/Pb, and similiar slopes of the rare earth elements. By varying the proportions of trapped liquid and thus simulating the cryptic and modal metasomatism observed in peridotite that surrounds these veins, we can model the variations in Ba/Nb, Ce/Pb, and Nb/U ratios that are observed in alkaline basalts. If the isotopic compositions of the initial low-degree peridotite melts are similar to the range observed in mid-ocean ridge basalt, our model calculations produce cumulates that would have isotopic compositions similar to those observed in most alkaline ocean island basalt (OIB) and continental magmas after ~0·15 Gyr. However, to produce alkaline basalts with HIMU isotopic compositions requires much longer residence times (i.e. 1–2 Gyr), consistent with subduction and recycling of metasomatized lithosphere through the mantle. EM magmas cannot readily be explained without appealing to other factors such as a heterogeneous asthenosphere. These modelling results support the interpretation proposed by various researchers that amphibole-bearing veins represent cumulates formed during the differentiation of a volatile-bearing low-degree peridotite melt and that these cumulates are significant components of the sources of alkaline OIB and continental magmas. The results of the forward models provide the potential for detailed tests of this class of hypotheses for the origin of alkaline magmas worldwide and for interpreting major and minor aspects of the geochemical variability of these magmas

Regional patterns in the paragenesis and age of inclusions in diamond, diamond composition, and the lithospheric seismic structure of Southern Africa

Abstract The Archean lithospheric mantle beneath the Kaapvaal -Zimbabwe craton of Southern Africa shows F 1% variations in seismic P-wave velocity at depths within the diamond stability field (150 -250 km) that correlate regionally with differences in the composition of diamonds and their syngenetic inclusions. Seismically slower mantle trends from the mantle below Swaziland to that below southeastern Botswana, roughly following the surface outcrop pattern of the Bushveld-Molopo Farms Complex. Seismically slower mantle also is evident under the southwestern side of the Zimbabwe craton below crust metamorphosed around 2 Ga. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa, and Jwaneng have Re -Os isotopic ages that range from circa 2.9 Ga to the Proterozoic and show little correspondence with these lithospheric variations. However, silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane do show some regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity correlates with a greater proportion of eclogitic versus peridotitic silicate inclusions in diamond, a greater incidence of younger Sm -Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds whereas the converse is true for diamonds from higher velocity mantle. The oldest formation ages of diamonds indicate that the mantle keels which became continental nuclei were created by middle Archean (3.2 -3.3 Ga) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of sulfide inclusions that are eclogitic in the 2.9 Ga age population links late Archean (2.9 Ga) subduction-accretion events involving an oceanic lithosphere component to craton stabilization. These events resulted in a widely distributed younger Archean generation of eclogitic diamonds in the lithospheric mantle. Subsequent Proterozoic tectonic and magmatic event

Lessons learned from monitoring of turbidity currents and guidance for future platform designs

Turbidity currents transport globally significant volumes of sediment and organic carbon into the deep-sea and pose a hazard to critical infrastructure. Despite advances in technology, their powerful nature often damages expensive instruments placed in their path. These challenges mean that turbidity currents have only been measured in a few locations worldwide, in relatively shallow water depths (≪2 km). Here, we share lessons from recent field deployments about how to design the platforms on which instruments are deployed. First, we show how monitoring platforms have been affected by turbidity currents including instability, displacement, tumbling and damage. Second, we relate these issues to specifics of the platform design, such as exposure of large surface area instruments within a flow and inadequate anchoring or seafloor support. Third, we provide recommended improvements to improve design by simplifying mooring configurations, minimising surface area, and enhancing seafloor stability. Finally we highlight novel multi-point moorings that avoid interaction between the instruments and the flow, and flow-resilient seafloor platforms with innovative engineering design features, such as ejectable feet and ballast. Our experience will provide guidance for future deployments, so that more detailed insights can be provided into turbidity current behaviour, and in a wider range of settings

Land‐To‐Sea Sediment Fluxes From a Major Glacial Lake Outburst Flood Were Stepped Rather Than Instantaneous

Glacial lake outburst floods can transport large volumes of sediment. Where these floods reach the coastline, much of the particulate matter is delivered directly to the marine environment. It has been suggested that offshore deposits, specifically in fjord settings, may provide a faithful record of past outburst flood events. However, a lack of observations means that the mechanics and the timing of sediment transport offshore following a glacial lake outburst event remain poorly constrained. Here, we document the changes in sea surface sediment dynamics following the 28 November 2020 Elliot Lake outburst flood in British Columbia, which transported ∼4.3 × 106 m3 of sediment into an adjacent fjord (Bute Inlet) as a deep nepheloid layer directly following the event. However, analysis of sea surface turbidity using in situ measurements and satellite‐derived estimates reveals that changes in fjord‐head surface turbidity in the months following the major flood were surprisingly small. The highest measured sea surface turbidity instead occurred 5 months after the initial outburst flood. This delayed increase in seaward sediment flux coincided with the onset of the spring freshet, when the discharge of the rivers feeding Bute Inlet increases each year. We suggest that large quantities of sediment were temporarily stored within the river catchment and were only remobilized when river discharge exceeded a threshold level following seasonal snowmelt. Our results reveal a temporal disconnect, where onshore to offshore transfer of sediment is stepped following a glacial lake outburst flood, which could complicate the sedimentology of subsequent deposits