Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

Sedimentary evolution of a Late Pleistocene temperate red algal reef (Coralligène) on Rhodes,Greece: correlation with global sea‐level fluctuations

Autochthonous red algal structures known as coralligène de plateau occur in the modern warm‐temperate Mediterranean Sea at water depths from 20 to 120 m, but fossil counterparts are not so well‐known. This study describes, from an uplifted coastal section at Plimiri on the island of Rhodes, a 450 m long by 10 m thick Late Pleistocene red algal reef (Coralligène Facies), interpreted as being a coralligène de plateau, and its associated deposits. The Coralligène Facies, constructed mainly by Lithophyllum and Titanoderma, sits unconformably upon the Plio‐Pleistocene Rhodes Formation and is overlain by a Maerl Facies (2 m), a Mixed Siliciclastic‐Carbonate Facies (0·2 m) and an Aeolian Sand Facies (2·5 m). The three calcareous facies, of Heterozoan character, are correlated with established members in the Lindos Acropolis Formation in the north of the island, while the aeolian facies is assigned to the new Plimiri Aeolianite Formation. The palaeoenvironmental and genetic‐stratigraphic interpretations of these mixed siliciclastic‐carbonate temperate water deposits involved consideration of certain characteristics associated with siliciclastic shelf and tropical carbonate shelf models, such as vertical grain‐size trends and the stratigraphic position of zooxanthellate coral growths. Integration of these results with electron spin resonance dates of bivalve shells indicates that the Coralligène Facies was deposited during Marine Isotope Stage 6 to 5e transgressive event (ca 135 to 120 ka), in water depths of 20 to 50 m, and the overlying Maerl Facies was deposited during regression from Marine Isotope Stage 5e to 5d (ca 120 to 110 ka), at water depths of 25 to 40 m. The capping Aeolian Sand Facies, involving dual terrestrial subunits, is interpreted as having formed during each of the glacial intervals Marine Isotope Stages 4 (71 to 59 ka) and 2 (24 to 12 ka), with soil formation during the subsequent interglacial periods of Marine Isotope Stages 3 and 1, respectively. Accumulation rates of about 0·7 mm year^?1 are estimated for the Coralligène Facies and minimum accumulation rates of 0·2 mm year^?1 are estimated for the Maerl Facies. The existence of older red algal reefs in the Plimiri region during at least Marine Isotope Stages 7 (245 to 186 ka) and 9 (339 to 303 ka) is inferred from the occurrence of reworked coralligène‐type lithoclasts in the basal part of the section and from the electron spin resonance ages of transported bivalve shells.     

Palaeoenvironmental significance of Late Permian palaeosols in the South‐Eastern Iberian Ranges,Spain

The Late Permian (Wuchiapingian) Alcotas Formation in the SE Iberian Ranges consists of one red alluvial succession where abundant soil profiles developed. Detailed petrographical and sedimentological studies in seven sections of the Alcotas Formation allow six different types of palaeosols, with distinctive characteristics and different palaeogeographical distribution, to be distinguished throughout the South‐eastern Iberian Basin. These characteristics are, in turn, related to topographic, climatic and tectonic controls. The vertical distribution of the palaeosols is used to differentiate the formation in three parts from bottom to top showing both drastic and gradual vertical upwards palaeoenvironmental changes in the sections. Reconstruction of palaeoenvironmental conditions based on palaeosols provides evidence for understanding the events that occurred during the Late Permian, some few millions of years before the well‐known Permian‐Triassic global crisis.     

Carbonate aeolianites of western Saurashtra,India: experimental decipherment of the depositional mechanisms

Granular carbonate deposits of Late Pleistocene to Early Holocene age, commonly referred to as ‘miliolite limestone’, occur in a linear belt, parallel to the southern coast of Saurashtra, India. In the present study area these carbonate deposits are found in select valleys between ridges and mounds of pyroclastic material present in the Deccan trap plateau. Two different depositional histories have been proposed for these sediments. The presence of marine bioclasts led to the postulation of a marine origin for these deposits. The second school of thought propounded redeposition of the coastal sediments by aeolian processes. Although a few features could not be explained by the proposed aeolian model, critical comparison of these two views favoured the aeolian origin. The mode of occurrence, lithological and structural attributes, and microscopic evidence presented here, also support a possible aeolian origin for these deposits. Experimental observation indicates that these carbonate aeolianites represent backflow deposits, which accumulated because of the flow separation caused by the presence of topographic highs. The conspicuous concave‐up geometry of the deposit conformed to the shape of the separation bulb. In view of the inferred depositional mechanism, the disposition of the deposits and the signature of the palaeoflow direction suggest that the carbonate particles were derived from the north‐western coast of Saurashtra by strong south‐easterly winds. Massive granular carbonates with outsized basement clasts appear to be the product of avalanching of granular material from the higher contours because of oversteepening of the primary deposit.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

Migration of the Ganga river and its implication on hydro-geological potential of Varanasi area, U.P., India

Borehole data reveals that during Late Quaternary, the Ganga river was non-existent in its present location near Varanasi. Instead, it was flowing further south towards peripheral craton. Himalayan derived grey micaceous sands were being carried by southward flowing rivers beyond the present day water divide of Ganga and mixed with pink arkosic sand brought by northward flowing peninsular rivers. Subsequently, the Ganga shifted to its present position and got incised. Near Varanasi, the Ganga river is flowing along a NW-SE tectonic lineament. The migration of Ganga river is believed to have been in response to basin expansion caused due to Himalayan tectonics during Middle Pleistocene times. Multi-storied sand bodies generated as a result of channel migration provide excellent aquifers confined by a thick zone of muddy sediments near the surface. Good quality potable water is available at various levels below about 70 m depth in sandy aquifers. Craton derived gravelly coarse-to-medium grained sand forms the main aquifer zones of tens of meter thickness with enormous yield. In contrast, the shallow aquifers made up of recycled interfluve silt and sandy silt occur under unconfined conditions and show water-level fluctuation of a few meters during pre-and post-monsoon periods.     

A modified DRASTIC model for Siting Confined Animal Feeding Operations in Williams County,Ohio, USA

Three of DRASTIC’s parameters (Depth to Water, Soil Media, and Topography) were modified and another parameter was added (land use/land cover) to the model to determine the potential impact on groundwater from Confined Animal Feeding Operations (CAFO) manure lagoon settings and manure application as fertilizer. Williams County is a mostly agricultural county located in northwest Ohio, USA. It currently has three CAFOs, all dairy, with the possibility of the construction of a multi-million chicken egg CAFO in the near future. A Geographic Information System (GIS) was utilized to modify the Ohio Department of Natural Resources (ODNR) DRASTIC map for the county to fully assess the county-wide pollution potential of CAFOs. The CAFO DRASTIC map indicates that almost half of Williams County has elevated groundwater pollution potential. The rest of the county, primarily the southeast corner, has lower CAFO groundwater pollution potential. Future CAFO development within the county should focus on the southeastern portion of the county where the groundwater table is deeper, and the aquifer is composed of shale substrate with low hydraulic conductivity. The CAFO DRASTIC results are intended to be used as a screening tool and are not to replace site-specific hydrogeologic investigations.     

Potential for localized groundwater contamination in a porous pavement parking lot setting in Rhode Island

The control of polluted surface runoff and the assessment of possible impacts on groundwater is a concern at the local and regional scale. On this background, a study investigates possible impacts of organic and inorganic pollutants (including bacteria) originating from a permeable asphalt parking lot on the water quality immediately beneath it. The functioning of the permeable pavement, including clogging and restricted vertical percolation, was also evaluated. Four nested sample ports (shallow and deep) were installed below low- and high-traffic areas, including one port outside the parking lot. At least initially there was a good hydraulic connection between the parking surface and the shallow sample ports. The presence of a geotextile layer at the base of the parking lot structure, however, was identified in lab tests as one factor restricting vertical percolation to the deeper ports. Clogging of the permeable surface was most pronounced in heavy traffic areas and below snow pile storage areas. Corroborated by high electric conductivity and chloride measurements, sand brought in by cars during winter was the principal cause for clogging. No bacteria or BOD were found in percolating water. Polycyclic aromatic hydrocarbons (PAH) were present at concentrations near minimum detection limit. Nutrients (nitrate and phosphate) were being leached into the ground via the permeable parking lot surface at annual flux rates of 0.45–0.84 g/m²/year. A multi-species tracer test demonstrated a retention capacity of the permeable parking lot structure of >90% for metals and 27% for nutrients, respectively.     

Low-pressure,water-assisted anatexis of basic dykes in a contact metamorphic aureole,Fuerteventura (Canary Islands): oxygen isotope evidence for a meteoric fluid origin

Migmatites produced by low-pressure anatexis of basic dykes are found in a contact metamorphic aureole around a pyroxenite–gabbro intrusion (PX2), on Fuerteventura. Dykes outside and inside the aureole record interaction with meteoric water, with low or negative δ¹⁸O whole-rock values (+0.2 to −3.4‰), decreasing towards the contact. Recrystallised plagioclase, diopside, biotite and oxides, from within the aureole, show a similar evolution with lowest δ¹⁸O values (−2.8, −4.2, −4.4 and −7.6‰, respectively) in the migmatite zone, close to the intrusion. Relict clinopyroxene phenocrysts preserved in all dykes, retain typically magmatic δ¹⁸O values up to the anatectic zone, where the values are lower and more heterogeneous. Low δ¹⁸O values, decreasing towards the intrusion, can be ascribed to the advection of meteoric water during magma emplacement, with increasing fluid/rock ratios (higher dyke intensities towards the intrusion acting as fluid-pathways) and higher temperatures promoting increasing exchange during recrystallisation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.