High‐resolution bathymetric mapping reveals subaqueous glacial landforms in the Arctic alpine lake Tarfala,Sweden

In Arctic alpine regions, glacio‐lacustrine environments respond sensitively to variations in climate conditions, impacting, for example,glacier extent and rendering former ice‐contact lakes into ice distal lakes and vice versa. Lakefloors may hold morphological records of past glacier extent, but remoteness and long periods of ice cover on such lakes make acquisition of high‐resolution bathymetric datasets challenging. Lake Tarfala and Kebnepakte Glacier, located in the Kebnekaise mountains, northern Sweden, comprise a small, dynamic glacio‐lacustrine system holding a climate archive that is not well studied. Using an autonomous surface vessel, a high‐resolution bathymetric dataset for Lake Tarfala was acquired in 2016, from which previously undiscovered end moraines and a potential grounding line feature were identified. For Kebnepakte Glacier, structure‐from‐motion photogrammetry was used to reconstruct its shape from photographs taken in 1910 and 1945. Combining these methods connects the glacial landform record identified at the lakefloor with the centennial‐scale dynamic behaviour of Kebnepakte Glacier. During its maximum 20^th century extent, attained c. 1910, Kebnepakte Glacier reached far into Lake Tarfala, but had retreated onto land by 1945, at an average of 7.9 m year^–1. Copyright © 2019 John Wiley & Sons, Ltd.     

Bleaching of the post‐IR IRSL signal: new insights for feldspar luminescence dating

Post‐infrared (pIR) stimulated luminescence dating of sedimentary feldspar largely avoids the effects of anomalous fading that affect conventional infrared stimulated luminescence (IRSL) dating. However, optical resetting of pIR signals is more difficult than resetting the conventional IRSL signal, which may undermine the crucial assumption that pIR signals were effectively bleached upon deposition and burial of sediment grains. In this study, we quantify the bleaching properties of several pIR signals on various samples using laboratory‐simulated bleaching in full sunlight and water‐attenuated sunlight. Our data show that bleaching is most efficient under full spectrum conditions for all pIR signals and that pIR signals measured at elevated temperature are increasingly harder to bleach than IR and pIR signals measured at low temperature (e.g. IR at 50°C). All bleaching curves exhibit a very slow and steady decrease, indicating that a fixed un‐bleachable residual level cannot be reached within the 11 days of solar simulator exposure undertaken here. We show that the magnitude of a laboratory‐determined residual dose depends on the adopted bleaching protocol and cannot be used as a proxy for the dose that remains in the sample at the time of burial (remnant dose). Our data emphasize the importance of finding a balance between sufficient signal stability and a minimized contribution of a remnant dose when using pIR procedures for feldspar luminescence dating.     

A multi‐storage groundwater concept for the SWAT model to emphasize nonlinear groundwater dynamics in lowland catchments

Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatio‐temporal tracer variability in the glacier melt end‐member — How does it affect hydrograph separation results?

Geochemical and isotopic tracers were often used in mixing models to estimate glacier melt contributions to streamflow, whereas the spatio‐temporal variability in the glacier melt tracer signature and its influence on tracer‐based hydrograph separation results received less attention. We present novel tracer data from a high‐elevation catchment (17 km², glacierized area: 34%) in the Oetztal Alps (Austria) and investigated the spatial, as well as the subdaily to monthly tracer variability of supraglacial meltwater and the temporal tracer variability of winter baseflow to infer groundwater dynamics. The streamflow tracer variability during winter baseflow conditions was small, and the glacier melt tracer variation was higher, especially at the end of the ablation period. We applied a three‐component mixing model with electrical conductivity and oxygen‐18. Hydrograph separation (groundwater, glacier melt, and rain) was performed for 6 single glacier melt‐induced days (i.e., 6 events) during the ablation period 2016 (July to September). Median fractions (±uncertainty) of groundwater, glacier melt, and rain for the events were estimated at 49±2%, 35±11%, and 16±11%, respectively. Minimum and maximum glacier melt fractions at the subdaily scale ranged between 2±5% and 76±11%, respectively. A sensitivity analysis showed that the intraseasonal glacier melt tracer variability had a marked effect on the estimated glacier melt contribution during events with large glacier melt fractions of streamflow. Intra‐daily and spatial variation of the glacier melt tracer signature played a negligible role in applying the mixing model. The results of this study (a) show the necessity to apply a multiple sampling approach in order to characterize the glacier melt end‐member and (b) reveal the importance of groundwater and rainfall–runoff dynamics in catchments with a glacial flow regime.     

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.     

Investigating suspended particulate matter in coastal waters using the fractal theory

Ocean Dynamics - Suspended particulate matters (SPM) in coastal waters were investigated with an approach combining suspended particulate matter concentrations (SPMCs) measured by an optical...     

Do landslides follow landslides? Insights in path dependency from a multi-temporal landslide inventory

Landslides are a major category of natural disasters, causing loss of lives, livelihoods and property. The critical roles played by triggering (such as extreme rainfall and earthquakes), and intrinsic factors (such as slope steepness, soil properties and lithology) have previously successfully been recognized and quantified using a variety of qualitative, quantitative and hybrid methods in a wide range of study sites. However, available data typically do not allow to investigate the effect that earlier landslides have on intrinsic factors and hence on follow-up landslides. Therefore, existing methods cannot account for the potentially complex susceptibility changes caused by landslide events. In this study, we used a substantially different alternative approach to shed light on the potential effect of earlier landslides using a multi-temporal dataset of landslide occurrence containing 17 time slices. Spatial overlap and the time interval between landslides play key roles in our work. We quantified the degree to which landslides preferentially occur in locations where landslides occurred previously, how long such an effect is noticeable, and how landslides are spatially associated over time. We also investigated whether overlap with previous landslides causes differences in landslide geometric properties. We found that overlap among landslides demonstrates a clear legacy effect (path dependency) that has influence on the landslide affected area. Landslides appear to cause greater susceptibility for follow-up landslides over a period of about 10  years. Follow-up landslides are on average larger and rounder than landslides that do not follow earlier slides. The effect of earlier landslides on follow-up landslides has implications for understanding of the landslides evolution and the assessment of landslide susceptibility.     

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.     

Prograded foredunes of Western Australia's macro‐tidal coast – implications for Holocene sea‐level change and high‐energy wave impacts

The Holocene evolution of the Canning Coast of Western Australia has largely been overlooked so far mainly due to its remoteness and low population density. We report on new data from a sequence of foredunes inside the macro‐tidal Admiral Bay, 110 km southwest of Broome. Based on sediment cores, differential global positioning system (dGPS)‐based elevation transects, and stratigraphical analyses on outcrops of the relict foredunes, we aim at reconstructing Holocene coastal changes and relative sea levels (RSLs), as well as identifying and dating imprints of extreme‐wave events. Sedimentary analyses comprise the documentation of bedding structures, foraminiferal content and macrofaunal remains, grain size distribution, and organic matter. The chronological framework is based on 26 carbon‐14 accelerator mass spectrometry (¹⁴C‐AMS) datings. Marine flooding of the pre‐Holocene surface landward of the 2.5 km‐wide foredune barriers occurred 7400–7200 cal bp , when mangroves colonized the area. After only 200–400 years, a high‐energy inter‐tidal environment established and prevailed until c. 4000 cal bp , before turning into the present supralittoral mudflat. During that time, coastal regression led to beach progradation and the formation of aligned foredunes. Drivers of progradation were a stable RSL or gradual RSL fall after the mid‐Holocene and a positive sand budget. The foredunes overlie upper beach deposits located up to >2 m above the present upper beach level and provide evidence for a higher mid‐Holocene RSL. Discontinuous layers of coarse shells and sand are intercalated in the foredunes, indicating massive coastal flooding events. One such layer was traced over three dune ridges and dated to c. 1700–1550 cal bp . However, it seems that most tropical cyclones induce net erosion rather than deposition at aligned foredunes and thus, they are only suitable for reconstructing temporal variability if erosional features or sedimentation reliably tied to these events can be identified and dated accurately. Copyright © 2014 John Wiley & Sons, Ltd.     

A Machine Learning Approach for Improving the Detection Capabilities at 3C Seismic Stations

We apply and evaluate a recent machine learning method for the automatic classification of seismic waveforms. The method relies on Dynamic Bayesian Networks (DBN) and supervised learning to improve the detection capabilities at 3C seismic stations. A time-frequency decomposition provides the basis for the required signal characteristics we need in order to derive the features defining typical “signal” and “noise” patterns. Each pattern class is modeled by a DBN, specifying the interrelationships of the derived features in the time-frequency plane. Subsequently, the models are trained using previously labeled segments of seismic data. The DBN models can now be compared against in order to determine the likelihood of new incoming seismic waveform segments to be either signal or noise. As the noise characteristics of seismic stations varies smoothly in time (seasonal variation as well as anthropogenic influence), we accommodate in our approach for a continuous adaptation of the DBN model that is associated with the noise class. Given the difficulty for obtaining a golden standard for real data (ground truth) the proof of concept and evaluation is shown by conducting experiments based on 3C seismic data from the International Monitoring Stations, BOSA and LPAZ.