Quantifying river incision into low‐relief surfaces using local and catchment‐wide 10Be denudation rates

Relief generation in non‐glaciated regions is largely controlled by river incision into bedrock but datable fluvial terraces that allow quantifying incision rates are not always present. Here we suggest a new method to determine river incision rates in regions where low‐relief surfaces are dissected by streams. The approach consists of three steps and requires the ¹⁰Be concentrations of a stream sediment sample and a regolith sample from the low‐relief surface. In the first step, the spatial distribution of ¹⁰Be surface concentrations in the given catchment is modelled by assuming that denudation rates are controlled by the local hillslope angles. The slope–denudation rate relation for this catchment is then quantified by adjusting the relation between slope angle and denudation rate until the average ¹⁰Be concentration in the model is equal to the one measured in the stream sediment sample. In the second step, curved swath profiles are used to measure hillslope angles adjacent to the main river channel. Third, the mean slope angle derived from these swath profiles and the slope–denudation relation are used to quantify the river incision rate (assuming that the incision rate equals the denudation rate on adjacent hillslopes). We apply our approach to two study areas in southern Tibet and central Europe (Black Forest). In both regions, local ¹⁰Be denudation rates on flat parts of the incised low‐relief surface are lower than catchment‐wide denudation rates. As the latter integrate across the entire landscape, river incision rates must exceed these spatially averaged denudation rates. Our approach yields river incision rates between ~15 and ~30 m/Ma for the Tibetan study area and incision rates of ~70 to ~100 m/Ma in the Black Forest. Taking the lowering of the low‐relief surfaces into account suggests that relief in the two study areas increases at rates of 10–20 and 40–70 m/Ma, respectively. Copyright © 2018 John Wiley & Sons, Ltd.     

Identification of extreme storm surges with high-impact potential along the German North Sea coastline

Planning and design of coastal protection rely on information about the probabilities of very severe storm tides and the possible changes that may occur in the course of climate change. So far, this information is mostly provided in the form of high percentiles obtained from frequency distributions or return values. More detailed information and assessments of events that may cause extreme damages or have extreme consequences at the coast are so far still unavailable. We describe and compare two different approaches that may be used to identify highly unlikely but still physically possible and plausible events from model simulations. Firstly, in the case when consistent wind and tide-surge data are available, different metrics such as the height of the storm surge can be derived directly from the simulated water levels. Secondly, in cases where only atmospheric data are available, the so called effective wind may be used. The latter is the projection of the horizontal wind vector on that direction which is most effective in producing surges at the coast. Comparison of events identified by both methods show that they can identify extreme events but that knowledge of the effective wind alone does not provide sufficient information to identify the highest storm surges. Tracks of the low-pressure systems over the North Sea need to be investigated to find those cases, where the duration of the high wind is too short to induce extreme storm tides. On the other hand, factors such as external surges or variability in mean sea level may enhance surge heights and are not accounted for in estimates based on effective winds only. Results from the analysis of an extended data set suggest that unprecedented storm surges at the German North Sea coast are possible even without taking effects from rising mean sea level into account. The work presented is part of the ongoing project “Extreme North Sea Storm Surges and Their Consequences” (EXTREMENESS) and represents the first step towards an impact assessment for very severe storm surges which will serve as a basis for development of adaptation options and evaluation criteria.     

Sparsity‐enhanced wavelet deconvolution

We propose a three‐step bandwidth enhancing wavelet deconvolution process, combining linear inverse filtering and non‐linear reflectivity construction based on a sparseness assumption. The first step is conventional Wiener deconvolution. The second step consists of further spectral whitening outside the spectral bandwidth of the residual wavelet after Wiener deconvolution, i.e., the wavelet resulting from application of the Wiener deconvolution filter to the original wavelet, which usually is not a perfect spike due to band limitations of the original wavelet. We specifically propose a zero‐phase filtered sparse‐spike deconvolution as the second step to recover the reflectivity dominantly outside of the bandwidth of the residual wavelet after Wiener deconvolution. The filter applied to the sparse‐spike deconvolution result is proportional to the deviation of the amplitude spectrum of the residual wavelet from unity, i.e., it is of higher amplitude; the closer the amplitude spectrum of the residual wavelet is to zero, but of very low amplitude, the closer it is to unity. The third step consists of summation of the data from the two first steps, basically adding gradually the contribution from the sparse‐spike deconvolution result at those frequencies at which the residual wavelet after Wiener deconvolution has small amplitudes. We propose to call this technique “sparsity‐enhanced wavelet deconvolution”. We demonstrate the technique on real data with the deconvolution of the (normal‐incidence) source side sea‐surface ghost of marine towed streamer data. We also present the extension of the proposed technique to time‐varying wavelet deconvolution.     

Localized shock-induced melting of sandstone at low shock pressures (<17.5 GPa): An experimental study

Shock-induced recovery experiments were performed to investigate melt formation in porous sandstones in the low shock pressure regime between 2.5 and 17.5 GPa. The sandstone shocked at 2.5 and 5 GPa is characterized by pore closure, fracturing of quartz (Qtz), and compression and deformation of phyllosilicates; no melting was observed. At higher pressures, five different types of melts were generated around pores and alongside fractures in the sandstone. Melting of kaolinite (Kln), illite (Ill), and muscovite (Ms) starts at 7.5, 12, and 15 GPa, respectively. The larger the amount of water in these minerals (Kln ~14 wt%, Ill ~6–10 wt%, and Ms ~4 wt% H₂O), the higher the shock compressibility and the lower the shock pressure required to induce melting. Vesicles in the almost dry silicate glasses attest to the loss of structural water during the short shock duration of the experiment. The compositions of the phyllosilicate-based glasses are identical to the composition of the parental minerals or their mixtures. Thus, this study has demonstrated that phyllosilicates in shocked sandstone undergo congruent melting during shock loading. In experiments at 10 GPa and higher, iron melt from the driver plate was injected into the phyllosilicate melts. During this process, Fe is partitioned from the metal droplets into the surrounding silicate melts, which induced unmixing of silicate melts with different chemical properties (liquid immiscibility). At pressures between 7.5 and 15 GPa, a pure SiO₂ glass was formed, which is located as short and thin bands within Qtz grains. These bands were shown to contain tiny crystals of experimentally generated stishovite.     

Continuous monitoring of volcanic eruption dynamics: a review of various techniques and new results from a frequency-modulated radar Doppler system

 In situ measurement of volcanic eruption velocities is one of the great challenges left in geophysical volcanology. In this paper we report on a new radar Doppler technique for monitoring volcanic eruption velocities. In comparison with techniques employed previously (e.g., photographic methods or acoustic Doppler measurements), this method allows continuous recordings of volcanic eruptions even during poor visibility. Also, radar Doppler instruments are usually light weight and energy efficient, which makes them superior to other Doppler techniques based on laser light or sound. The proposed new technique was successfully tested at Stromboli Volcano in late 1996 during a period of low activity. The recorded data allow a clear distinction between particles rising from the vent and particles falling back towards the vent. The mean eruption velocity was approximately 10 m/s. Most of the eruptions recorded by radar were correlated to seismic recordings. The correlation between the magnitude of the volcanic shocks and the eruption force index defined in the paper may provide new insights into magma transport in the conduit. Received: 15 May 1998 / Accepted: 15 December 1998     

Microscopic evidence of stishovite generated in low‐pressure shock experiments on porous sandstone: Constraints on its genesis

It has been almost exactly half a century since the first synthesis of stishovite in shock experiments on quartz was reported, but its formation conditions during shock is still under debate. Here, we present direct transmission electron microscopic observation of stishovite within material recovered from high‐explosive shock experiments on porous sandstone shocked at 7.5 and 12.5 GPa. Our observations allow for new conclusions on the genesis of stishovite in a close‐to‐nature environment. The formation of stishovite in short‐time shock experiments proves that its crystallization is ultrafast (<1 μs). Crystals were found only embedded in amorphous veins indicating homogeneous nucleation. Crystallization from melt rather than from glass can be concluded from the observation of roundish, defect‐free crystals up to 150 nm in diameter embedded in nondensified glass. The formation of stishovite at 7.5 GPa is in accordance with the phase diagram of silica, if rapid undercooling is present that becomes only possible by the existence of small hot spots in an otherwise cold material, which is supported by transient heat calculation. The absence of coesite at 7.5 GPa suggests kinetic hindrance of its crystallization from melt and, thus, smaller critical cooling rates compared to stishovite where critical cooling rates are estimated to be as large as 10¹¹ K s^?1. While the amorphous veins containing stishovite represent unambiguously hot spots, no associated pressure amplification could be verified within these veins. The rapid liquidus crystallization of stishovite only in hot spots generated in porous material is an alternative formation mechanism to the widely accepted theory of solid–solid transition from quartz to stishovite and might represent the more general mechanism occurring in nature for low shock pressure events.     

Microfacies and depositional environment of Tertiary Tanjung Enim low rank coal, South Sumatra Basin, Indonesia

The South Sumatra basin is among the most important coal producing basins in Indonesia. Results of an organic petrography study on coals from Tanjung Enim, South Sumatra Basin are reported. The studied low rank coals have a mean random huminite reflectance between 0.35% and 0.46% and are dominated by huminite (34.6–94.6 vol.%). Less abundant are liptinite (4.0–61.4 vol.%) and inertinite (0.2–43.9 vol.%). Minerals are found only in small amounts (0–2 vol.%); mostly as iron sulfide.Based on maceral assemblages, the coals can be grouped into five classes: (1) humotelinite-rich group, (2) humodetrinite-rich group, (3) humocollinite-rich group, (4) inertinite-rich group and (5) humodetrinite–liptinite-rich group. Comparing the distribution of maceral assemblages to the maceral or pre-maceral assemblages in modern tropical domed peat in Indonesia reveals many similarities. The basal section of the studied coal seams is represented typically by the humodetrinite–liptinite-rich group. This section might be derived from sapric or fine hemic peat often occurring at the base of modern peats. The middle section of the seams is characterized by humotelinite-rich and humocollinite-rich groups. The precursors of these groups were hemic and fine hemic peats. The top section of the coal seams is typically represented by the humodetrinite-rich or inertinite-rich group. These groups are the counterparts of fibric peat at the top of the modern peats. The sequence of maceral assemblages thus represents the change of topogenous to ombrogenous peat and the development of a raised peat bog.A comparison between the result of detailed maceral assemblage analysis and the paleodepositional environment as established from coal maceral ratio calculation indicates that the use of coal maceral ratio diagrams developed for other coal deposits fails to deduce paleo-peat development for these young tropical coals. In particular, mineral distribution and composition should not be neglected in coal facies interpretations.     

The behavior of Mg, Fe, and Ni during the replacement of olivine by orthopyroxene: experiments relevant to mantle metasomatism

Replacement of olivine by orthopyroxene is a frequently observed phenomenon in mantle metasomatism. In order to study element redistribution in SiO₂ metasomatism we synthesised orthopyroxene reaction rims at the contacts between forsterite-rich olivine and quartz. The orthopyroxene rims grew from the original quartz-olivine interface into both directions implying counterdiffusion of iron/magnesium and silicon. Following local equilibrium partitioning the X_Fe is lower in the orthopyroxene than in the reactant olivine at the olivine-orthopyroxene replacement front. The resulting local iron excess is compensated by formation of orthopyroxene with a higher X_Fe at the quartz-orthopyroxene interface, which is out of equilibrium with the reactant olivine. This is facilitated through short circuit diffusion along grain boundaries within the orthopyroxene rim. Due to the low capacity of orthopyroxene to accommodate Ni, this component is forced to diffuse back into the olivine producing a Ni enriched zone ahead of the replacement front. This leads to Ni contents in the orthopyroxene rim, which are higher than what is expected in equilibrium with the unaltered olivine. Taking quartz as a proxy for a silica rich fluid or liquid metasomatising agent, we conclude that the overall element fractionation between olivine and the silica rich phase may deviate from equilibrium partitioning so that the Fe and Ni concentrations in the orthopyroxene which is in contact with quartz are higher than in equilibrium with the reactant olivine. This indicates that kinetic fractionation is important for the chemical evolution of both the mantle rocks and the metasomatising agents.     

Potential decadal predictability and its sensitivity to sea ice albedo parameterization in a global coupled model

Decadal prediction is one focus of the upcoming 5th IPCC Assessment report. To be able to interpret the results and to further improve the decadal predictions it is important to investigate the potential predictability in the participating climate models. This study analyzes the upper limit of climate predictability on decadal time scales and its dependency on sea ice albedo parameterization by performing two perfect ensemble experiments with the global coupled climate model EC-Earth. In the first experiment, the standard albedo formulation of EC-Earth is used, in the second experiment sea ice albedo is reduced. The potential prognostic predictability is analyzed for a set of oceanic and atmospheric parameters. The decadal predictability of the atmospheric circulation is small. The highest potential predictability was found in air temperature at 2?m height over the northern North Atlantic and the southern South Atlantic. Over land, only a few areas are significantly predictable. The predictability for continental size averages of air temperature is relatively good in all northern hemisphere regions. Sea ice thickness is highly predictable along the ice edges in the North Atlantic Arctic Sector. The meridional overturning circulation is highly predictable in both experiments and governs most of the decadal climate predictability in the northern hemisphere. The experiments using reduced sea ice albedo show some important differences like a generally higher predictability of atmospheric variables in the Arctic or higher predictability of air temperature in Europe. Furthermore, decadal variations are substantially smaller in the simulations with reduced ice albedo, which can be explained by reduced sea ice thickness in these simulations.     

Comparison of precipitation in the regional climate model BALTIMOS to radar observations

Observational data and simulations of the regional climate system Baltic integrated model system (BALTIMOS) were used to study precipitation in the Baltic Sea and its drainage basin with a special focus on the diurnal cycle. The study includes a general evaluation of BALTIMOS precipitation, showing that BALTIMOS has too many light rain events causing an overestimation of the total annual precipitation amount. The diurnal cycle as well as its spatial distribution was analysed. BALTIMOS captures the broad characteristics: a significant diurnal variability with an afternoon peak above land and weak variability with a nocturnal peak above sea. An algorithm to distinguish between frontal and convective precipitation was applied to examine the diurnal cycle more thoroughly. The local solar time of maximum rain in summer is about 1 to 2 h earlier in BALTIMOS than in radar observations of precipitation.