Isochron‐burial dating of glaciofluvial deposits: First results from the Swiss Alps

In this study, we use isochron‐burial dating to date the Swiss Deckenschotter, the oldest Quaternary deposits of the northern Alpine Foreland. Concentrations of cosmogenic ¹⁰Be and ²⁶Al in individual clasts from a single stratigraphic horizon can be used to calculate an isochron‐burial age based on an assumed initial ratio and the measured ²⁶Al/¹⁰Be ratio. We suggest that, owing to deep and repeated glacial erosion, the initial isochron ratio of glacial landscapes at the time of burial varies between 6.75 and 8.4. Analysis of 22 clasts of different lithology, shape, and size from one 0.5 m thick gravel bed at Siglistorf (Canton Aargau) indicates low nuclide concentrations: <20 000 ¹⁰Be atoms/g and <150 000 ²⁶Al atoms/g. Using an ²⁶Al/¹⁰Be ratio of 7.6 (arithmetical mean of 6.75 and 8.4), we calculate a mean isochron‐burial age of 1.5 ± 0.2 Ma. This age points to an average bedrock incision rate between 0.13 and 0.17 mm/a. Age data from the Irchel, Stadlerberg, and Siglistorf sites show that the Higher Swiss Deckenschotter was deposited between 2.5 and 1.3 Ma. Our results indicate that isochron‐burial dating can be successfully applied to glaciofluvial sediments despite very low cosmogenic nuclide concentrations. Copyright © 2017 John Wiley & Sons, Ltd.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

The probabilistic seismic hazard assessment of Germany—version 2016, considering the range of epistemic uncertainties and aleatory variability

The basic seismic load parameters for the upcoming national design regulation for DIN EN 1998-1/NA result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering (DIBt). This 2016 version of the national seismic hazard assessment for Germany is based on a comprehensive involvement of all accessible uncertainties in models and parameters and includes the provision of a rational framework for integrating ranges of epistemic uncertainties and aleatory variabilities in a comprehensive and transparent way. The developed seismic hazard model incorporates significant improvements over previous versions. It is based on updated and extended databases, it includes robust methods to evolve sets of models representing epistemic uncertainties, and a selection of the latest generation of ground motion prediction equations. The new earthquake model is presented here, which consists of a logic tree with 4040 end branches and essential innovations employed for a realistic approach. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. Seismic load parameters, for rock conditions of \(v_{S30}\) = 800 m/s, are calculated for three hazard levels (10, 5 and 2% probability of occurrence or exceedance within 50 years) and delivered in the form of uniform hazard spectra, within the spectral period range 0.02–3 s, and seismic hazard maps for peak ground acceleration, spectral response accelerations and for macroseismic intensities. Results are supplied as the mean, the median and the 84th percentile. A broad analysis of resulting uncertainties of calculated seismic load parameters is included. The stability of the hazard maps with respect to previous versions and the cross-border comparison is emphasized.     

Decadal evolution of tidal flats and channels in the Outer Weser estuary,Germany

This study focuses on the medium scale morphodynamics of the tidal flat and channel system Fedderwarder Priel, located in the Outer Weser estuary (Wadden Sea, Germany). Tidal channels and adjacent flats are highly dynamic systems whose morphologic evolution are driven by tidal, wind, and wave forcings. These coastal environments are an important ecosystem and react to changes in hydrodynamic conditions in various spatial and temporal scales. Based on annual medium-resolution digital elevation models from 1998 to 2016, we describe changes in the surface area over depth with hypsometries and use vertical dynamic trends in order to analyze and visualize the morphologic evolution of the Fedderwarder Priel and adjacent tidal channels. It is shown that several intertidal flats rise in the order of 1.3 to 5.6 cm/year. The findings indicate that the Outer Weser estuary was not in an equilibrium state for the investigated period, and tidal flats accreted with a rate exceeding mean sea level rise.     

The fusion crust of the Winchcombe meteorite: A preserved record of atmospheric entry processes

Fusion crusts form during the atmospheric entry heating of meteorites and preserve a record of the conditions that occurred during deceleration in the atmosphere. The fusion crust of the Winchcombe meteorite closely resembles that of other stony meteorites, and in particular CM2 chondrites, since it is dominated by olivine phenocrysts set in a glassy mesostasis with magnetite, and is highly vesicular. Dehydration cracks are unusually abundant in Winchcombe. Failure of this weak layer is an additional ablation mechanism to produce large numbers of particles during deceleration, consistent with the observation of pulses of plasma in videos of the Winchcombe fireball. Calving events might provide an observable phenomenon related to meteorites that are particularly susceptible to dehydration. Oscillatory zoning is observed within olivine phenocrysts in the fusion crust, in contrast to other meteorites, perhaps owing to temperature fluctuations resulting from calving events. Magnetite monolayers are found in the crust, and have also not been previously reported, and form discontinuous strata. These features grade into magnetite rims formed on the external surface of the crust and suggest the trapping of surface magnetite by collapse of melt. Magnetite monolayers may be a feature of meteorites that undergo significant degassing. Silicate warts with dendritic textures were observed and are suggested to be droplets ablated from another stone in the shower. They, therefore, represent the first evidence for intershower transfer of ablation materials and are consistent with the other evidence in the Winchcombe meteorite for unusually intense gas loss and ablation, despite its low entry velocity.     

MODFLOW as a Configurable Multi-Model Hydrologic Simulator

MODFLOW 6 is the latest in a line of six “core” versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW-USG, the Newton-Raphson formulation in MODFLOW-NWT, and the support for partitioned stress boundaries in MODFLOW-CDSS. The flexible multi-model capability allows users to configure MODFLOW 6 simulations to represent the local-grid refinement (LGR) capabilities available in MODFLOW-LGR, the multi-species transport capabilities in MT3DMS, and the coupled variable-density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code.     

A case study of using Raman lidar measurements in high-accuracy GPS applications

This paper investigates the impact of rapid small-scale water vapor fluctuations on GPS height determination. Water vapor measurements from a Raman lidar are used for documenting the water vapor heterogeneities and correcting GPS signal propagation delays in clear sky conditions. We use data from four short observing sessions (6 h) during the VAPIC experiment (15 May–15 June 2004). The retrieval of wet delays from our Raman lidar is shown to agree well with radiosonde retrievals (bias and standard deviation (SD) were smaller than 1 and 2.8 mm, respectively) and microwave radiometers (from two different instruments, bias was 6.0/−6.6 mm and SD 1.3/3.8 mm). A standard GPS data analysis is shown to fail in accurately reproducing fast zenith wet delay (ZWD) variations. The ZWD estimates could be improved when mean post-fit phase residuals were removed. Several methodologies for integrating zenith lidar observations into the GPS data processing are also presented. The final method consists in using lidar wet delays for correcting a priori the GPS phase observations and estimating a scale factor for the lidar wet delays jointly with the GPS station position. The estimation of this scale factor allows correcting for a mis-calibration in the lidar data and provides in the same way an estimate of the Raman lidar instrument constant. The agreement of this constant with an independent determination using radiosonde data is at the level of 1–4%. The lidar wet delays were derived by ray-tracing from zenith pointing measurements: further improvement in GPS positioning is expected from slant path lidar measurements that would properly account for water vapor anisotropy.     

History of scoria-cone eruptions on the eastern shoulder of the Kenya–Tanzania Rift revealed in the 250-ka sediment record of Lake Chala near Mount Kilimanjaro

Reconstructions of the timing and frequency of past eruptions are important to assess the propensity for future volcanic activity, yet in volcanic areas such as the East African Rift only piecemeal eruption histories exist. Understanding the volcanic history of scoria-cone fields, where eruptions are often infrequent and deposits strongly weathered, is particularly challenging. Here we reconstruct a history of volcanism from scoria cones situated along the eastern shoulders of the Kenya–Tanzania Rift, using a sequence of tephra (volcanic ash) layers preserved in the ~250-ka sediment record of Lake Chala near Mount Kilimanjaro. Seven visible and two non-visible (crypto-) tephra layers in the Lake Chala sequence are attributed to activity from the Mt Kilimanjaro (northern Tanzania) and the Chyulu Hills (southern Kenya) volcanic fields, on the basis of their glass chemistry, textural characteristics and known eruption chronology. The Lake Chala record of eruptions from scoria cones in the Chyulu Hills volcanic field confirms geological and historical evidence of its recent activity, and provides first-order age estimates for seven previously unknown eruptions. Long and well-resolved sedimentary records such as that of Lake Chala have significant potential for resolving regional eruption chronologies spanning hundreds of thousands of years.     

Are carbonate barrier islands mobile? Insights from a mid to late-Holocene system,Al Ruwais,northern Qatar

Barrier islands are important landforms in many coastal systems around the globe. Studies of modern barrier island systems are mostly limited to those of siliciclastic realms, where the islands are recognized as mobile features that form on transgressive coastlines and migrate landward as sea-level rises. Barrier islands of the ‘Great Pearl Bank’ along the United Arab Emirates coast are the best-known carbonate examples. These Holocene islands, however, are interpreted to be anchored by older deposits and immobile. The mid-Holocene to late-Holocene depositional system at Al Ruwais, northern Qatar, provides an example of a mobile carbonate barrier island system, perhaps more similar to siliciclastic equivalents. Sedimentological and petrographic analyses, as well as ¹⁴C-dating of shells and biogenic remains from vibracored sediments and surface deposits, show that after 7000 years ago a barrier system with a narrow back-barrier lagoon formed along what is now an exposed coastal zone, while, contemporaneously, a laterally-extensive coral reef was forming immediately offshore. After 1400 years ago the barrier system was forced to step ca 3 km seaward in response to a sea-level fall of less than 2 m, where it re-established itself directly on the mid-Holocene reef. Since that time, the barrier has retreated landward as much as 1000 m to its current position, exposing previously-deposited back-barrier lagoonal sediment at the open-coast shoreline. In modern neritic warm-water carbonate settings mobile barrier island systems are rare. Their construction and migration may be inhibited by reef formation, early cementation, and the relative inefficiency of sourcing beach sediments from open carbonate shelves. Carbonate barrier island systems likely formed more commonly during geological periods when ramps and unrimmed shelves predominated and in calcite seas, when meteoric cementation was minimized as a result of initial calcitic allochem mineralogy. As with their siliciclastic analogues, however, recognition of the influence of these transient landforms in the rock record is challenging.     

Comparison of environmental tracers including organic micropollutants as groundwater exfiltration indicators into a small river of a karstic catchment

Understanding groundwater–surface water (GW–SW) interactions is vital for water management in karstic catchments due to its impact on water quality. The objective of this study was to evaluate and compare the applicability of seven environmental tracers to quantify and localize groundwater exfiltration into a small, human-impacted karstic river system. Tracers were selected based on their emission source to the surface water either as (a) dissolved, predominantly geogenic compounds (radon-222, sulphate and electrical conductivity) or (b) anthropogenic compounds (predominantly) originating from wastewater treatment plant (WWTP) effluents (carbamazepine, tramadol, sodium, chloride). Two contrasting sampling approaches were compared (a) assuming steady-state flow conditions and (b) considering the travel time of the water parcels (Lagrangian sampling) through the catchment to account for diurnal changes in inflow from the WWTP. Spatial variability of the concentrations of all tracers indicated sections of preferential groundwater inflow. Lagrangian sampling techniques seem highly relevant for capturing dynamic concentration patterns of WWTP-derived compounds. Quantification of GW inflow with the finite element model FINIFLUX, based on observed in-stream Rn activities led to plausible fluxes along the investigated river reaches (0.265 m³ s⁻¹), while observations of other natural or anthropogenic environmental tracers produced less plausible water fluxes. Important point sources of groundwater exfiltration can be ascribed to locations where the river crosses geological fault lines. This indicates that commonly applied concepts describing groundwater–surface water interactions assuming diffuse flow in porous media are difficult to transfer to karstic river systems whereas concepts from fractured aquifers may be more applicable. In general, this study helps selecting the best suited hydrological tracer for GW exfiltration and leads to a better understanding of processes controlling groundwater inflow into karstic river systems.