Topographic signatures of progressive glacial landscape transformation

The Pleistocene glaciations left a distinct topographic footprint in mountain ranges worldwide. The geometric signature of glacial topography has been quantified in various ways, but the temporal development of landscape metrics has not been traced in a landscape evolution model so far. However, such information is needed to interpret the degree of glacial imprint in terms of the integrated signal of temporal and spatial variations in erosion as a function of glacial occupation time. We apply a surface process model for cold-climate conditions to an initially fluvial mountain range. By exploring evolving topographic patterns in model time series, we determine locations where topographic changes reach a maximum and where the initial landscape persists. The signal of glacial erosion, expressed by the overdeepening of valleys and the steepening of valley flanks, develops first at the glacier front and migrates upstream with ongoing glacial erosion. This leads to an increase of mean channel slope and its variance. Above steep flanks and head-walls, however, the observed mean channel slope remains similar to the mean channel slope of the initial fluvial topography. This leads to a characteristic turning point in the channel slope–elevation distribution above the equilibrium line altitude, where a transition from increasing to decreasing channel slope with elevation occurs. We identify this turning point and a high channel slope variance as diagnostic features to quantify glacial imprint. Such features are abundant in glacially imprinted mid-latitude mountain ranges such as the Eastern Alps. By analysing differently glaciated parts of the mountain range, we observe a decreasing clarity of this diagnostic morphometric property with decreasing glacial occupation. However, catchments of the unglaciated eastern fringe of the Alps also feature turning points in their channel slope–elevation distributions, but in contrast to the glaciated domain, the variance of channel slope is small at all elevation levels.     

Factors Controlling Total Concentration and Aqueous Speciation of Aluminium in an Acidic Headwater Stream of the Bavarian Forest National Park: a Modelling Approach

The purpose of the present paper is to analyse factors controlling total concentration and aqueous speciation of aluminium in the Große Ohe River, using a thermodynamic equilibrium model and a mixing approach. A model compound for humic substances is derived on the basis of the relation between anion deficit and the organic carbon content in the river as well as literature data. An equilibrium speciation model for aluminium is set up, considering this model compound and relevant inorganic solutes. Although the model cannot be verified directly, its results may be viewed as qualitatively correct. Applying the model to measured stream water samples highlights that aqueous speciation of aluminium is mainly controlled by the pH value and discharge and that free aluminium concentrations reach clearly toxic levels during acidic episodes. Comparing measured concentrations of sulfate and H⁺ and calculated concentrations of Al³⁺ with solubility curves of gibbsite like minerals and jurbanite clearly shows that total aluminium concentrations are not controlled by equilibria with these mineral phases alone. The observed relationship can be better explained from a mixture of two distinct waters, representing lowflow and highflow chemistry, and the resulting equilibrium concentrations. This indicates that total aluminium concentration, in particular during high discharge events, is mainly controlled by the mixture of waters with differing chemistry and flowpaths.     

Measurements of translation, rotation and strain: new approaches to seismic processing and inversion

We propose a novel approach to seismic tomography based on the joint processing of translation, strain and rotation measurements. Our concept is based on the apparent S and P velocities, defined as the ratios of displacement velocity and rotation amplitude, and displacement velocity and divergence amplitude, respectively. To assess the capability of these new observables to constrain various aspects of 3D Earth structure, we study their corresponding finite-frequency kernels, computed with a combination of spectral-element simulations and adjoint techniques. The principal conclusion is that both the apparent S and P velocities are generally sensitive only to small-scale near-receiver structure, irrespective of the type of seismic wave considered. It follows that knowledge of deeper Earth structure would not be required in tomographic inversions for local structure based on the new observables. In a synthetic finite-perturbation test, we confirm the ability of the apparent S and P velocities to directly detect both the location and the sign of shallow lateral velocity variations.     

Layer-Based Modelling of the Earth’s Gravitational Potential up to 10-km Scale in Spherical Harmonics in Spherical and Ellipsoidal Approximation

Global forward modelling of the Earth’s gravitational potential, a classical problem in geophysics and geodesy, is relevant for a range of applications such as gravity interpretation, isostatic hypothesis testing or combined gravity field modelling with high and ultra-high resolution. This study presents spectral forward modelling with volumetric mass layers to degree 2190 for the first time based on two different levels of approximation. In spherical approximation, the mass layers are referred to a sphere, yielding the spherical topographic potential. In ellipsoidal approximation where an ellipsoid of revolution provides the reference, the ellipsoidal topographic potential (ETP) is obtained. For both types of approximation, we derive a mass layer concept and study it with layered data from the Earth2014 topography model at 5-arc-min resolution. We show that the layer concept can be applied with either actual layer density or density contrasts w.r.t. a reference density, without discernible differences in the computed gravity functionals. To avoid aliasing and truncation errors, we carefully account for increased sampling requirements due to the exponentiation of the boundary functions and consider all numerically relevant terms of the involved binominal series expansions. The main outcome of our work is a set of new spectral models of the Earth’s topographic potential relying on mass layer modelling in spherical and in ellipsoidal approximation. We compare both levels of approximations geometrically, spectrally and numerically and quantify the benefits over the frequently used rock-equivalent topography (RET) method. We show that by using the ETP it is possible to avoid any displacement of masses and quantify also the benefit of mapping-free modelling. The layer-based forward modelling is corroborated by GOCE satellite gradiometry, by in-situ gravity observations from recently released Antarctic gravity anomaly grids and degree correlations with spectral models of the Earth’s observed geopotential. As the main conclusion of this work, the mass layer approach allows more accurate modelling of the topographic potential because it avoids 10–20-mGal approximation errors associated with RET techniques. The spherical approximation is suited for a range of geophysical applications, while the ellipsoidal approximation is preferable for applications requiring high accuracy or high resolution.