Flow and turbulence structure across the ripple–dune transition: an experiment under mobile bed conditions

Current knowledge of flow and turbulent processes acting across the sand bed continuum is still unable to unequivocally explain the mechanism(s) by which ripples become dunes. Understanding has been improved by comparative high-resolution studies undertaken over fixed bedforms at different stages in the continuum. However, these studies both ignore the role of mobile sediment and do not examine flow structure during the actual transition from ripples to dunes. The aims of the paper are: (i) to describe flow and turbulence characteristics acting above mobile bedforms at several stages across the transition; and (ii) to compare these data with those arising from experiments over fixed ripples and dunes. Laboratory experiments are presented that examine the turbulence structure across seven distinct stages of the transition from ripples to dunes. Single-point acoustic Doppler velocimeter sampling at three flow heights above a developing mobile boundary was undertaken. Time-averaged statistics and the instantaneous quadrant record reveal distinct changes in flow structure either side of the change from ripples to dunes. Initially, shear-related, high-frequency vortex shedding dominates turbulence production. This increases until two-dimensional (2D) dunes have formed. Thereafter, turbulence intensities and Reynolds stress decline and three-dimensional dunes exhibit values found over 2D ripples. This is the result of shear layer dampening which occurs when the topographically-accelerated downstream velocity increases at a faster rate than flow depth. Activity at reattachment increases due to high velocity fluid imparting high mass and momentum transfer at the bed and/or wake flapping. Suspended sediment may also play a role in turbulence dampening and bed erosion. Ejections dominate over sweeps in terms of event frequency but not magnitude. Strong relationships between inward interactions and sweeps, and ejections and outward interactions, suggest that mass and momentum exchanges are dependent upon activity in all four quadrants. The results contradict the notion present in most physical models that larger bedforms exhibit most shear layer activity. Consequently an improved model for the ripple–dune transition is proposed.     

Pleistocene speleothems of Mallorca: implications for palaeoclimate and carbonate diagenesis in mixing zones

The Pleistocene speleothems of Sa Bassa Blanca cave, Mallorca, are excellent indicators of palaeoclimate variations, and are samples that allow evaluation of the products and processes of mixing‐zone diagenesis in an open‐water cave system. Integrated stratigraphic, petrographic and geochemical data from a horizontal core of speleothem identified two main origins for speleothem precipitates: meteoric‐marine mixing zone and meteoric‐vadose zone. Mixing‐zone precipitates formed at and just below the water–air interface of cave pools during interglacial times, when the cave was flooded as a result of highstand sea‐level. Mixing‐zone precipitates include bladed and dendritic high‐Mg calcite, microporous‐bladed calcite with variable Mg content, and acicular aragonite; their presence suggests that calcium‐carbonate cementation is significant in the studied mixing‐zone system. Fluid inclusion salinities, δ¹³C and δ¹⁸O compositions of the mixing‐zone precipitates suggest that mixing ratio was not the primary control on whether precipitation or dissolution occurred, rather, the proximity to the water table and degassing of CO₂ at the interface, were the major controls on precipitation. Thus, simple two‐end‐member mixing models may apply only in mixing zones well below the water table. Meteoric‐vadose speleothems include calcite and high‐Mg calcite with columnar and bladed morphologies. Vadose speleothems precipitated during glacial stages when sea level was lower than present. Progressive increase in δ¹³C and δ¹⁸O of the vadose speleothems resulted from cooling temperatures and more positive seawater δ¹⁸O associated with glacial buildup. Such covariation could be considered as a valid alternative to models predicting invariant δ¹⁸O and highly variable δ¹³C in meteoric calcite. Glacio‐eustatic oscillations of sea‐level are recorded as alternating vadose and mixing‐zone speleothems. Short‐term climatic variations are recorded as alternating aragonite and calcite speleothems precipitated in the mixing zone. Fluid‐inclusion and stable‐isotope data suggest that aragonite, as opposed to calcite, precipitated during times of reduced meteoric recharge.     

Structure,composition, and location of organic matter in the enstatite chondrite Sahara 97096 (EH3)

Abstract– The insoluble organic matter (IOM) of an unequilibrated enstatite chondrite Sahara (SAH) 97096 has been investigated using a battery of analytical techniques. As the enstatite chondrites are thought to have formed in a reduced environment at higher temperatures than carbonaceous chondrites, they constitute an interesting comparative material to test the heterogeneities of the IOM in the solar system and to constrain the processes that could affect IOM during solar system evolution. The SAH 97096 IOM is found in situ: as submicrometer grains in the network of fine‐grained matrix occurring mostly around chondrules and as inclusions in metallic nodules, where the carbonaceous matter appears to be more graphitized. IOM in these two settings has very similar δ¹⁵N and δ¹³C; this supports the idea that graphitized inclusions in metal could be formed by metal catalytic graphitization of matrix IOM. A detailed comparison between the IOM extracted from a fresh part and a terrestrially weathered part of SAH 97096 shows the similarity between both IOM samples in spite of the high degree of mineral alteration in the latter. The isolated IOM exhibits a heterogeneous polyaromatic macromolecular structure, sometimes highly graphitized, without any detectable free radicals and deuterium‐heterogeneity and having mean H‐ and N‐isotopic compositions in the range of values observed for carbonaceous chondrites. It contains some submicrometer‐sized areas highly enriched in ¹⁵N (δ¹⁵N up to 1600‰). These observations reinforce the idea that the IOM found in carbonaceous chondrites is a common component widespread in the solar system. Most of the features of SAH 97096 IOM could be explained by the thermal modification of this main component.     

Risk analyses for a water supply system: Occoquan Reservoir,Fairfax and Prince William Counties,Virginia, USA / Des analyses aléatoires d'un système d'approvisionnement en eau: Occoquan Reservoir,Fairfax et Prince William Counties,Virginia, USA

Abstract

This report demonstrates two techniques for evaluating risks in the operation of a water supply system. Both rely on reconstructed historical streamflow data to develop estimates of the probabilities of certain specific events occurring in the future. These techniques are applied to the Occoquan Reservoir which was experiencing an unprecedented low level of storage in the autumn of 1977. The two techniques are used respectively to evaluate the overall adequacy of the existing reservoir and to evaluate the risks in the 1977 crisis.

In the first technique, the general risk analysis model (GRAM), simulations of the reservoir's contents are carried out under a set of assumptions about withdrawal rates and emergency procedures. The results of the GRAM simulation for the Occoquan Reservoir are in the form of estimates of the probabilities that in any year contain emergency procedures will have to be invoked. These estimates are given for a range of rates of withdrawals and for four different stages of emergency actions. Also given are the estimated probabilities of entering emergency conditions in one year given that an emergency has occurred in the previous year. Due to the year-to-year persistence of low flows, these latter (conditional) probabilities are higher than the former (marginal) probabilities.

The second technique used is position analysis. In this procedure probability distributions of future storages are estimated under existing storage conditions and an assumed rate of withdrawal from the reservoir. The position analysis which was initialized at the 1 October 1977 conditions indicates that the probability of entering a Stage III emergency (the prohibition of all uses of water non-essential to life, health, and safety) in the autumn of 1977 or winter of 1978 was 10 per cent at that time. With a reduction in water use by 8 million gallons per day, however, this probability would have fallen to 4 per cent.

If a long reconstructed historical streamflow record is available to a water supply agency, then the agency will have the capability to undertake its own risk analyses. It can carry out comparisons of alternative operating policies by using techniques such as GRAM. It can also evaluate the short-term risks of different plans of operation during crisis situations by using position analysis.     

Capturing and modelling metre‐scale spatial facies heterogeneity in a Jurassic ramp setting (Central High Atlas,Morocco)

Each simulation algorithm, including Truncated Gaussian Simulation, Sequential Indicator Simulation and Indicator Kriging is characterized by different operating modes, which variably influence the facies proportion, distribution and association of digital outcrop models, as shown in clastic sediments. A detailed study of carbonate heterogeneity is then crucial to understanding these differences and providing rules for carbonate modelling. Through a continuous exposure of Bajocian carbonate strata, a study window (320 m long, 190 m wide and 30 m thick) was investigated and metre‐scale lithofacies heterogeneity was captured and modelled using closely‐spaced sections. Ten lithofacies, deposited in a shallow‐water carbonate‐dominated ramp, were recognized and their dimensions and associations were documented. Field data, including height sections, were georeferenced and input into the model. Four models were built in the present study. Model A used all sections and Truncated Gaussian Simulation during the stochastic simulation. For the three other models, Model B was generated using Truncated Gaussian Simulation as for Model A, Model C was generated using Sequential Indicator Simulation and Model D was generated using Indicator Kriging. These three additional models were built by removing two out of eight sections from data input. The removal of sections allows direct insights on geological uncertainties at inter‐well spacings by comparing modelled and described sections. Other quantitative and qualitative comparisons were carried out between models to understand the advantages/disadvantages of each algorithm. Model A is used as the base case. Indicator Kriging (Model D) simplifies the facies distribution by assigning continuous geological bodies of the most abundant lithofacies to each zone. Sequential Indicator Simulation (Model C) is confident to conserve facies proportion when geological heterogeneity is complex. The use of trend with Truncated Gaussian Simulation is a powerful tool for modelling well‐defined spatial facies relationships. However, in shallow‐water carbonate, facies can coexist and their association can change through time and space. The present study shows that the scale of modelling (depositional environment or lithofacies) involves specific simulation constraints on shallow‐water carbonate modelling methods.     

Magnetic mineral transport and sorting in the swash‐zone: northern Lake Erie,Canada

A combined field and laboratory study in northern Lake Erie has provided new insights into the origin and dynamics of heavy mineral placer deposits on beaches consisting primarily of non‐magnetic sediment. Work was conducted on the cross‐shore and longshore transport of heavy magnetic minerals using magnetic susceptibility and fluorescent paints to trace the movement, in the field, of samples of magnetic (magnetite) and non‐magnetic (quartz and calcite) grains, respectively. Laboratory experiments examined how the burial of small, dense magnetic minerals is affected by the grain size of the non‐magnetic host material, and how grain burial affects magnetic susceptibility measurements at the surface. The field experiments demonstrated that the magnetic mineral tracers were buried rapidly beneath coarser, non‐magnetic grains under low to moderate wave conditions, and subsequently were unable to move in the longshore or cross‐shore directions. The laboratory experiments showed that the magnetic susceptibility rapidly decreased with the rate and depth of burial of the magnetic minerals, and that magnetic grain burial was most effective beneath coarser rather than finer non‐magnetic sand and, for the latter sediments, under less rather than more energetic conditions. The results imply that magnetic mineral concentrations develop in this area through magnetic grain burial under fairly mild conditions, and subsequent settling, exposure and concentration in the upper swash zone during more energetic periods, when the non‐magnetic grains are eroded. It is probably during these erosional periods, when the magnetic minerals are exposed in fairly homogeneous deposits, that longshore and cross‐shore transport takes place.     

Climate‐controlled aggradation and cyclicity of continental loessic siliciclastic sediments in Asselian–Sakmarian cyclothems,Permian, Hugoton embayment,USA

Siliciclastic intervals in Lower Permian carbonate–siliciclastic cyclothems in western Kansas record climate control on facies progression, deposition and preservation. The 26 000 km² study area comprises seven marine‐continental (carbonate–siliciclastic) cyclothems caused by glacioeustasy. Core data and a three‐dimensional geological model provide a detailed view of the sub‐surface on a gently sloping ramp. Siliciclastic intervals in the cyclothems are fine‐grained red beds with extensive pedogenic features, indicating a continental origin. Bed geometry (sheet‐like deposits that thin to the east), lateral grading, grain size (very fine‐grained sand to silt) and grain angularity (sub‐angular to angular) suggest that the sediment is loess sourced from the west, probably the Ancestral Rocky Mountains. There is a repeated record of glacial‐cycle‐scale, climate‐controlled cyclicity within siliciclastic intervals that has not been recognized previously. Aeolian silt grain size coarsens upward towards the middle, then fines upward in each siliciclastic interval. When sea‐level was high (interglacial) and carbonate production flourished, aeolian sedimentation nearly ceased, suggesting increased vegetation and rainfall at the source. As sea‐level fell, fine‐grained siliciclastic sediments were deposited under relatively dry, but seasonally wet conditions on an exposed ramp. Laterally graded coarser grained siliciclastic sediments with diagnostic fabrics indicate drier conditions with seasonal rainfall during a continued relative fall in sea‐level. The coarsest siliciclastic sediments were deposited during the lowest sea‐level and driest conditions, but still with sufficient seasonal moisture to allow vegetative cover and bioturbation. Subsequent upward fining is correlated with sedimentological indications of wetter conditions during relative sea‐level rise. Unlike common sequence stratigraphic models that relate siliciclastic sediment accumulation to base‐level rise, continental deposits were preserved because plants and pedogenesis stabilized aeolian sediment. The aggradational landscape formed by this process had several metres of positive relief that reduced accommodation for overlying marine carbonate strata. Thus, this mechanism for continental siliciclastic aggradation has a significant effect on sequence stratigraphic architecture.     

Mapping Submerged Aquatic Vegetation with GIS in the Caloosahatchee Estuary: Evaluation of Different Interpolation Methods

This article evaluates different spatial interpolation methods for mapping submerged aquatic vegetation (SAV) in the Caloosahatchee Estuary, Florida. Data used for interpolation were collected by the Submersed Aquatic Vegetation Early Warning System (SAVEWS). The system consists of hydro-acoustic equipment, which operates from a slow-moving boat and records bottom depth, seagrass height, and seagrass density. This information is coupled with geographic location coordinates from a Global Positioning System (GPS) and stored together in digital files, representing SAV status at points along transect lines. Adequate spatial interpolation is needed to present the SAV information, including density, height, and water depth, as spatially continuous data for mapping and for comparison between seasons and years. Interpolation methods examined in this study include ordinary kriging with five different semivariance models combined with a variable number of neighboring points, the inverse distance weighted (IDW) method with different parameters, and the triangulated irregular network (TIN) method with linear and quintic options. Interpolation results were compared with survey data at selected calibration transects to examine the suitability of different interpolation methods. Suitability was quantified by the determination coefficient (R2) and the root-mean-square error (RMSE) between interpolated and observed values. The most suitable interpolation method was identified as the one yielding the highest R2 value and/or the lowest RMSE value. For different geographic conditions, seasons, and SAV parameters, different interpolation methods were recommended. This study identified that kriging was more suitable than the IDW or TIN method for spatial interpolation of all SAV parameters measured. It also suggested that transect data with irregular spatial distribution patterns such as SAV parameters are sensitive to interpolation methods. An inappropriate interpolation method such as TIN can lead to erroneous spatial representation of the SAV status. With a functional geographic system and adequate computing power, the evaluation and selection of interpolation methods can be automated and quantitative, leading to a more efficient and accurate decision.     

Impact histories of angrites,eucrites, and their parent bodies

Abstract– Eucrites, which are probably from 4 Vesta, and angrites are the two largest groups of basaltic meteorites from the asteroid belt. The parent body of the angrites is not known but it may have been comparable in size to Vesta as it retained basalts and had a core dynamo. Both bodies were melted early by ²⁶Al and formed basalts a few Myr after they accreted. Despite these similarities, the impact histories of the angrites and eucrites are very different: angrites are very largely unshocked and none are breccias, whereas most eucrites are breccias and many are shocked. We attribute the lack of shocked and unbrecciated angrites to an impact, possibly at 4558 Myr ago—the radiometric age of the younger angrites—that extracted the angrites from their original parent body into smaller bodies. These bodies, which may have had a diameter of approximately 10 km, suffered much less impact damage than Vesta during the late heavy bombardment because small bodies retain shocked rocks less efficiently than large ones and because large bodies suffer near‐catastrophic impacts that deposit vastly more impact energy per kg of target. Our proposed history for the angrites is comparable to that proposed by Bogard and Garrison (2003) for the unbrecciated eucrites with Ar‐Ar ages of 4.48 Gyr and that for unbrecciated eucrites with anomalous oxygen isotopic compositions that did not come from Vesta. We infer that the original parent bodies of the angrites and the anomalous eucrites were lost from the belt when the giant planets migrated and the total mass of asteroids was severely depleted. Alternatively, their parent bodies may have formed in the terrestrial planet region and fragments of these bodies were scattered out to the primordial Main Belt as a consequence of terrestrial planet formation.