Quantifying fine‐grained sediment sources in the River Axe catchment,southwest England: application of a Monte Carlo numerical modelling framework incorporating local and genetic algorithm optimisation

Increasing recognition of the deleterious environmental effects of excessive fine sediment delivery to watercourses means that reliable sediment source assessment represents a fundamental component of catchment planning targeting the protection of freshwater resources and their ecological integrity. Sediment tracing or fingerprinting approaches have been increasingly used to provide catchment scale sediment source information, but there is a need to continue refining existing procedures especially with respect to uncertainty analysis during mass balance modelling. Consequently, an updated Monte Carlo numerical modelling framework was devised and tested, incorporating both conventional and robust statistics coupled with random and Latin Hypercube Sampling (LHS) together with local and genetic algorithm (GA) optimisation. A sediment sourcing study undertaken in the River Axe catchment, southwest England, suggested that the use of robust statistics and LHS with GA optimisation generated the best performance with respect to predicting measured bed sediment geochemistry in six out of eight model applications. On this basis, the catchment‐wide average median sediment source contributions were predicted to be 38 ± 1% (pasture topsoils), 3 ± 1% (cultivated topsoils), 37 ± 1% (damaged road verges) and 22 ± 1% (channel banks/subsurface sources). Using modelling frameworks which provide users with flexibility to compare local and global optimisation during uncertainty analysis is recommended for future sediment tracing studies. Copyright © 2011 John Wiley & Sons, Ltd.     

The timing of migmatization in the northern Arabian–Nubian Shield: Evidence for a juvenile sedimentary component in collision‐related batholiths

Collision‐related granitoid batholiths, like those of the Hercynian and Himalayan orogens, are mostly fed by magma derived from metasedimentary sources. However, in the late Neoproterozoic calcalkaline (CA) batholiths of the Arabian–Nubian Shield (ANS), which constitutes the northern half of the East African orogen, any sedimentary contribution is obscured by the juvenile character of the crust and the scarcity of migmatites. Here, we use paired in situ LASS‐ICP‐MS measurements of U–Th–Pb isotope ratios and REE contents of monazite and xenotime and SHRIMP‐RG analyses of separated zircon to demonstrate direct linkage between migmatites and granites in the northernmost ANS. Our results indicate a single prolonged period of monazite growth at 640–600 Ma, in metapelites, migmatites and peraluminous granites of three metamorphic suites: Abu‐Barqa (SW Jordan), Roded (S Israel) and Taba–Nuweiba (Sinai, Egypt). The distribution of monazite dates and age zoning in single monazite grains in migmatites suggest that peak thermal conditions, involving partial melting, prevailed for c. 10 Ma, from 620 to 610 Ma. REE abundances in monazite are well correlated with age, recording garnet growth and garnet breakdown in association with the prograde and retrograde stages of the melting reactions, respectively. Xenotime dates cluster at 600–580 Ma, recording retrogression to greenschist facies conditions as garnet continued to destabilize. Phase equilibrium modelling and mineral thermobarometry yield P–T conditions of ~650–680°C and 5–7 kbar, consistent with either water‐fluxed or muscovite‐breakdown melting. The expected melt production is 8–10 vol.%, allowing a melt connectivity network to form leading to melt segregation and extraction. U–Pb ages of zircon rims from leucosomes indicate crystallization of melt at 610 ± 10 Ma, coinciding with the emplacement of a vast volume of CA granites throughout the northern ANS, which were previously considered post‐collisional. Similar monazite ages (c. 620 Ma) retrieved from the amphibolite facies Elat schist indicate that migmatites are the result of widespread regional rather than local contact metamorphism, representing the climax of the East African orogenesis.     

Holocene and Lateglacial relative sea‐level change in north‐west England: implications for glacial isostatic adjustment models

Relative sea‐level (RSL) change is reconstructed for central Cumbria, UK, based on litho‐ and biostratigraphical analysis from the Lateglacial to the late Holocene. The RSL curve is constrained using ten new radiocarbon‐dated sea‐level index points in addition to published data. The sea‐level curve identifies a clear Lateglacial sea‐level highstand approximately 2.3 m OD at c. 15–17 k cal a BP followed by rapid RSL fall to below ?5 m OD. RSL then rose rapidly during the early Holocene culminating in a mid‐Holocene highstand of approximately 1 m OD at c. 6 k cal a BP followed by gradual fall to the present level. These new data provide an important test for the RSL predictions from glacial isostatic adjustment models, particularly for the Lateglacial where there are very little data from the UK. The new RSL curve shows similar broad‐scale trends in RSL movement predicted by the models. However, the more recent models fail to predict the Lateglacial sea level highstand above present reconstructed by the new data presented here. Future updates to the models are needed to reduce this mismatch. This study highlights the importance for further RSL data to constrain Lateglacial sea level from sites in northern Britain. Copyright © 2012 John Wiley & Sons, Ltd.     

Impacts of model initialization on an integrated surface water–groundwater model

Integrated hydrologic models characterize catchment responses by coupling the subsurface flow with land surface processes. One of the major areas of uncertainty in such models is the specification of the initial condition and its influence on subsequent simulations. A key challenge in model initialization is that it requires spatially distributed information on model states, groundwater levels and soil moisture, even when such data are not routinely available. Here, the impact of uncertainty in initial condition was explored across a 208 km² catchment in Denmark using the ParFlow.CLM model. The initialization impact was assessed under two meteorological conditions (wet vs dry) using five depth to water table and soil moisture distributions obtained from various equilibrium states (thermal, root zone, discharge, saturated and unsaturated zone equilibrium) during the model spin‐up. Each of these equilibrium states correspond to varying computation times to achieve stability in a particular aspect of the system state. Results identified particular sensitivity in modelled recharge and stream flow to the different initializations, but reduced sensitivity in modelled energy fluxes. Analysis also suggests that to simulate a year that is wetter than the spin‐up period, an initialization based on discharge equilibrium is adequate to capture the direction and magnitude of surface water–groundwater exchanges. For a drier or hydrologically similar year to the spin‐up period, an initialization based on groundwater equilibrium is required. Variability of monthly subsurface storage changes and discharge bias at the scale of a hydrological event show that the initialization impacts do not diminish as the simulations progress, highlighting the importance of robust and accurate initialization in capturing surface water–groundwater dynamics. Copyright © 2015 John Wiley & Sons, Ltd.     

Chronology,stratigraphy and palaeoenvironmental interpretation of a Late Pleistocene to mid‐Holocene cave accumulation on Kangaroo Island,South Australia

Chronological, sedimentological and geochemical analyses of a clastic infill from Kelly Hill Cave (5K1), Kangaroo Island, document a palaeoenvironmental record that spans from the Late Pleistocene to the middle Holocene. We AMS radiocarbon‐dated bone collagen and U–Th‐dated speleothem to determine that fossiliferous sediments were deposited between >20 ka and 7 ka ago. Most of the 15 sedimentary layers are dominated by sand‐ and silt‐sized quartz that is physically and geochemically comparable with surface soils in the Kelly Hill area. Late Pleistocene and Last Glacial Maximum strata are represented primarily by homogeneous, poorly sorted quartz‐rich sediments that contain little organic matter, but include a thin layer composed largely of silt‐sized clay pellets that resemble sediments deflated from playa lakes. Microstructures observed in petrographic slides indicate that, with the exception of one layer, all sediments experienced little reworking once deposited in the cave. Some layers display pedogenic microstructures such as redeposited clays and opaline silica infilling that indicate postdepositional modification; that is, cave‐floor soil development. Overlying Holocene‐aged sediments also consist mainly of quartz but have much greater organic matter content. Some of these sediments have been strongly influenced by re‐precipitated organic matter that appears to have been transported into the cave via vadose drip water. The presence of dissolved organic matter in soil/vadose waters suggests a high vegetation density and acidic soils, which are congruent with the more equitable climatic conditions characteristic of the Holocene. The sediments described here provide a valuable palaeoenvironmental record that will facilitate future interpretation of associated vertebrate fossils.     

An automated routing methodology to enable direct rainfall in high resolution shallow water models

Recent high profile flood events have highlighted the need for hydraulic models capable of simulating pluvial flooding in urban areas. This paper presents a constant velocity rainfall routing scheme that provides this ability within the LISFLOOD‐FP hydraulic modelling code. The scheme operates in place of the shallow water equations within cells where the water depth is below a user‐defined threshold, enabling rainfall‐derived water to be moved from elevated features such as buildings or curbstones without causing instabilities in the solution whilst also yielding a reduction in the overall computational cost of the simulation. Benchmarking against commercial modelling packages using a pluvial and point‐source test case demonstrates that the scheme does not impede the ability of LISFLOOD‐FP to match both predicted depths and velocities of full shallow water models. The stability of the scheme in conditions unsuitable for traditional two‐dimensional hydraulic models is then demonstrated using a pluvial test case over a complex urban digital elevation model containing buildings. Deterministic single‐parameter sensitivity analyses undertaken using this test case show limited sensitivity of predicted water depths to both the chosen routing speed within a physically plausible range and values of the depth threshold parameter below 10 mm. Local instabilities can occur in the solution if the depth threshold is >10 mm, but such values are not required even when simulating extreme rainfall rates. The scheme yields a reduction in model runtime of ~25% due to the reduced number of cells for which the hydrodynamic equations have to be solved. Copyright © 2012 John Wiley & Sons, Ltd.     

Bioavailability of jarosite for stimulating acid mine drainage attenuation

Biological reduction of iron-sulfate minerals, such as jarosite, has the potential to contribute to the natural attenuation of acid mine drainage (AMD) sites. Previous studies of AMD attenuation at Davis Mine, an abandoned pyrite mine in Rowe Massachusetts, provided evidence of iron and sulfate reduction by indigenous bacteria. Jarosite is a large component of the sediment at Davis Mine and may play a role in AMD attenuation. In this study, microcosms were constructed with groundwater and sediment from Davis Mine and amended with glycerol, nitrogen and phosphorus (GNP) and naturally formed natrojarosite. Over time, higher total iron, sulfate, pH and sodium concentrations and lower oxidation–reduction potentials were observed in microcosms amended with GNP and jarosite, compared with unamended microcosms and killed controls. Geochemical modeling predicted jarosite precipitation under microcosm conditions, suggesting that abiotic processes were unlikely contributors to jarosite dissolution. SEM imaging at the jarosite surface showed microbial attachment. Microbial community composition analysis revealed a shift to higher populations of Clostridia, which are known to reduce both iron and sulfate. The results show that jarosite may be utilized as an electron acceptor by iron and/or sulfate reducing bacteria at Davis Mine and its presence may aid in the attenuation of AMD.     

Fulgurite morphology: a classification scheme and clues to formation

Fulgurites are natural glasses formed by cloud-to-ground lightning. Several different morphologies of fulgurites have been reported in previous studies, including sand fulgurites, rock fulgurites, and clay fulgurites. Herein, we examine sand, clay, and caliche fulgurites and demonstrate that these differ systematically in their morphology. We further use morphological features to constrain properties of fulgurite-forming lightning strikes. We classify fulgurites into four types of morphologies with an additional minor type. Type I fulgurites are sand fulgurites consisting of thin, glass walls; type II fulgurites are clay fulgurites, consisting of thick, melt rich walls; type III fulgurites are caliche fulgurites, consisting of thick, glass poor walls; and type IV fulgurites are rock fulgurites, consisting of glasses with walls consisting of surrounding, unmelted rock. Fulgurite morphology shows that the energy of fulgurite-forming strikes is between 1 and 30 MJ/m of fulgurite formed, suggests heating rates in the order of 1,000 K/s, and lightning channel thicknesses of about 1 mm diameter. Lightning generates mixtures of at least two components in most fulgurites: an SiO₂ glass identified as lechatelierite and a groundmass of more varied composition. In addition to these four primary types, a fifth type—droplet fulgurites—is morphologically dissimilar from the other types, but is compositionally related to the type II or IV fulgurites. Additionally, two fulgurites, both from York County, Pennsylvania, USA, showed the reduction of iron to iron metal with an assortment of Fe–Ti and Si–P compounds with stoichiometry that ranges from nearly pure Fe metal to FeSi. These metal silicides include stoichiometric Fe₃Si, Fe₂Si, and Fe₅Si₃, and possibly Fe₈Si₃ and Fe₇Si₃, and provide a terrestrial source for these phases, which are typically associated with extraterrestrial material.