A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Temperature measured in permafrost borehole at ICDP site 5011-3 in El'gygytgyn Crater, 2008-2011

This study focuses on the temperature field observed in boreholes drilled as part of interdisciplinary scientific campaign targeting the El'gygytgyn Crater Lake in NE Russia. Temperature data are available from two sites: the lake borehole 5011-1 located near the center of the lake reaching 400 m depth, and the land borehole 5011-3 at the rim of the lake, with a depth of 140 m. Constraints on permafrost depth and past climate changes are derived from numerical simulation of the thermal regime associated with the lake-related talik structure. The thermal properties of the subsurface needed for these simulations are based on laboratory measurements of representative cores from the quaternary sediments and the underlying impact-affected rock, complemented by further information from geophysical logs and data from published literature. The temperature observations in the lake borehole 5011-1 are dominated by thermal perturbations related to the drilling process, and thus only give reliable values for the lowermost value in the borehole. Undisturbed temperature data recorded over more than two years are available in the 140 m deep land-based borehole 5011-3. The analysis of these observations allows determination of not only the recent mean annual ground surface temperature, but also the ground surface temperature history, though with large uncertainties. Although the depth of this borehole is by far too insufficient for a complete reconstruction of past temperatures back to the Last Glacial Maximum, it still affects the thermal regime, and thus permafrost depth. This effect is constrained by numerical modeling: assuming that the lake borehole observations are hardly influenced by the past changes in surface air temperature, an estimate of steady-state conditions is possible, leading to a meaningful value of 14 ± 5 K for the post-glacial warming. The strong curvature of the temperature data in shallower depths around 60 m can be explained by a comparatively large amplitude of the Little Ice Age (up to 4 K), with low temperatures prevailing far into the 20th century. Other mechanisms, like varying porosity, may also have an influence on the temperature profile, however, our modeling studies imply a major contribution from recent climate changes.

Particle positions of numerical experiments 1 to 4

In mixed sediment beds, erosion resistance can change relative to that of beds composed of a uniform sediment because of varying textural and/or other grain-size parameters, with effects on pore water flow that are difficult to quantify by means of analogue techniques. To overcome this difficulty, a three-dimensional numerical model was developed using a finite difference method (FDM) flow model coupled with a distinct element method (DEM) particle model. The main aim was to investigate, at a high spatial resolution, the physical processes occurring during the initiation of motion of single grains at the sediment-water interface and in the shallow subsurface of simplified sediment beds under different flow velocities. Increasing proportions of very fine sand (D50=0.08 mm) were mixed into a coarse sand matrix (D50=0.6 mm) to simulate mixed sediment beds, starting with a pure coarse sand bed in experiment 1 (0 wt% fines), and proceeding through experiment 2 (6.5 wt% fines), experiment 3 (10.5 wt% fines), and experiment 4 (28.7 wt% fines). All mixed beds were tested for their erosion behavior at predefined flow velocities varying in the range of U 1-5=10-30 cm/s. The experiments show that, with increasing fine content, the smaller particles increasingly fill the spaces between the larger particles. As a consequence, pore water inflow into the sediment is increasingly blocked, i.e., there is a decrease in pore water flow velocity and, hence, in the flow momentum available to entrain particles. These findings are portrayed in a new conceptual model of enhanced sediment bed stabilization.

Benthic oxygen flixes on the Washington shelf and slope

Benthic oxygen fluxes calculated from in situ microelectrode profiles arc compared with benthic flux chamber O2 uptake measurements on a transect of eight stations across the continental shelf and three stations on the slope of Washington State. Station depths ranged from 40 to 630 m and bottom-water oxygen concentrations were 127-38 µM. The fluxes measured by the two methods were similar on the slope, but on the shelf, the chamber flux exceeded the microelectrode flux by as much as a factor of 3-4. We attribute this difference to pore-water irrigation, a process which apparently accounts for the oxidation of a significant amount of organic C in the continental shelf sediments. Combining our diffusive flux data with other data demonstrates clearly that the bottomwater oxygen concentration must play some significant role in determining the sedimentary oxygen consumption rate. Numerical simulation of the microelectrode 0, profiles suggests that roughly half the diffusive 0, flux must be consumed within - 1 mm of the sediment surface. If this conclusion is correct, then the magnitude of the diffusive flux depends both on the bottom-water oxygen concentration and on the supply rate of labile C to the sediment surf'ace.

Simulated transition to agropastoralism in the Indus valley 7500-3000 BC

The Indus Valley Civilization (IVC) was one of the first great civilizations in prehistory. This bronze age civilization flourished from the end of the fourth millennium BC. It disintegrated during the second millennium BC; despite much research effort, this decline is not well understood. Less research has been devoted to the emergence of the IVC, which shows continuous cultural precursors since at least the seventh millennium BC. To understand the decline, we believe it is necessary to investigate the rise of the IVC, i.e., the establishment of agriculture and livestock, dense populations and technological developments 7000-3000 BC. Although much archaeologically typed information is available, our capability to investigate the system is hindered by poorly resolved chronology, and by a lack of field work in the intermediate areas between the Indus valley and Mesopotamia. We thus employ a complementary numerical simulation to develop a consistent picture of technology, agropastoralism and population developments in the IVC domain. Results from this Global Land Use and technological Evolution Simulator show that there is (1) fair agreement between the simulated timing of the agricultural transition and radiocarbon dates from early agricultural sites, but the transition is simulated first in India then Pakistan; (2) an independent agropas- toralism developing on the Indian subcontinent; and (3) a positive relationship between archeological artifact richness and simulated population density which remains to be quantified.