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.     

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.     

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.     

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.     

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.     

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.     

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.     

Paired 10Be sampling of polished bedrock and erratic boulders to improve dating of glacial landforms: an example from the Western Alps

Cosmogenic nuclide dating of glacial landforms may lead to ambiguous results for ice retreat histories. The persistence of significant cosmogenic concentrations inherited from previous exposure may increase the apparent exposure ages for polished bedrocks affected by limited erosion under ice and for erratic boulders transported by glaciers and previously exposed in high-altitude rock walls. In contrast, transient burying by moraines, sediments and snow decreases the apparent exposure age. We propose a new sampling strategy, applied to four sites distributed in the Arc and Arve valleys in the Western Alps, to better constrain the factors that can bias exposure ages associated with glacial processes. We used the terrestrial cosmogenic nuclide ¹⁰Be (TCN) to estimate the exposure time from paired sampling of depth profiles in polished bedrock and on overlying erratic boulders. For a given sampling site, the exposure ages for both the polished bedrock and boulder are expected to be the same. However, in six cases out of seven, boulders had significantly higher ¹⁰Be surface concentrations than those of the associated polished surfaces. In present and past glacial processes, the ¹⁰Be distribution with depth for boulders and bedrocks implies the presence of an inheritance concentration of ¹⁰Be. Our study suggests that ¹⁰Be concentrations in erratic boulders and in polished bedrocks provide maximum and minimum exposure ages of the glacial retreat, respectively. © 2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Regression-based distribution mapping for bias correction of climate model outputs using linear quantile regression

Stochastic Environmental Research and Risk Assessment - Impact models are a major source of information for quantifying the consequences of future climate change for humans and the environment. To...     

Hydraulic-Head Formulation for Density-Dependent Flow and Transport

Density-dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic-head formulation can facilitate adaptation of existing constant-density groundwater flow codes to include density-driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density-dependent terms as a compartmentalized “correction” to constant-density calculations already performed by the code. The hydraulic-head formulation also accommodates complexities such as unconfined groundwater flow and Newton-Raphson solution schemes more readily than the freshwater-head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic-head formulation in MODFLOW.