Impact of an 0.2 km3 Rock Avalanche on Lake Eibsee (Bavarian Alps,Germany) – Part II: Catchment Response to Consecutive Debris Avalanche and Debris Flow

The ~0.2 km³ Eibsee rock avalanche impacted Paleolake Eibsee and completely displaced its waters. This study analyses the lake impact and the consequences, and the catchment response to the landslide. A quasi-3D seismic reflection survey, four sediment cores from modern Lake Eibsee, reaching far down into the rock avalanche mass, nine radiocarbon ages, and geomorphic analysis allow us to distinguish the main rock avalanche event from a secondary debris avalanche and debris flow. The highly fluidized debris avalanche formed a megaturbidite and multiple swashes that are recorded in the lake sediments. The new calibrated age for the Eibsee rock avalanche of ~4080–3970 cal yr BP indicates a coincidence with rockslides in the Fernpass cluster and subaquatic landslides in Lake Piburg and Lake Plansee, and raises the possibility that a large regional earthquake triggered these events. We document a complex history of erosion and sedimentation in Lake Eibsee, and demonstrate how the catchment response and rebirth of the lake are revealed through the complementary application of geophysics, sedimentology, radiocarbon dating, and geomorphology. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

A new pseudocylindrical equal-area projection for adaptive composite map projections

The recently introduced adaptive composite map projection technique changes the projection to the geographic area shown on a map. It is meant as a replacement for the commonly used web Mercator projection, which grossly distorts areas when representing the entire world. The original equal-area version of the adaptive composite map projection technique uses the Lambert azimuthal projection for regional maps and three alternative projections for world maps. Adaptive composite map projections can include a variety of other equal-area projections when the transformation between the Lambert azimuthal and the world projections uses Wagner’s method. To select the most suitable pseudocylindrical projection, the distortion characteristics of a pseudocylindrical projection family are analyzed, and a user study among experts in the area of map projections is carried out. Based on the results of the distortion analysis and the user study, a new pseudocylindrical projection is recommended for extending adaptive composite map projections. The new projection is equal-area throughout the transformation to the Lambert azimuthal projection and has better distortion characteristics then small-scale projections currently included in the adaptive composite map projection technique.     

On the probability distribution of GNSS carrier phase observations

When processing observational data from global navigation satellite systems (GNSS), the carrier phase measurements are generally assumed to follow a normal distribution. Although full knowledge of the probability distribution of the observables is not required for parameter estimation, for example when using the least-squares method, the distributional properties of GNSS observations play a key role in quality control procedures, such as outlier and cycle-slip detection, in ambiguity resolution, as well as in the reliability assessment of estimation results. In addition, when applying GNSS positioning under critical observation conditions with respect to multipath and atmospheric effects, the validity of the normal distribution assumption of GNSS observables certainly comes into doubt. This paper illustrates the discrepancies between the normal distribution assumption and reality, based on a large and representative data set of GPS phase measurements covering a range of factors, including multipath impact, baseline length, and atmospheric conditions. The statistical inferences are made using the first through fourth sample moments, hypothesis tests, and graphical tools such as histograms and quantile–quantile plots. The results show clearly that multipath effects, in particular the near-field component, produce the dominant influence on the distributional characteristics of GNSS observables. Additionally, using surface meteorological data, considerable correlations between distributional deviations from normality on the one hand and atmospheric relative humidity on the other are detected.     

Mineral grain availability and pupal-case building by lotic caddisflies: Effects on case architecture, stability and building expenses

The retreat-making larvae of many lotic caddisflies build entirely new pupal cases with fine gravel and sand that they collect in the neighbourhood of the building place to fix it with silk to cobbles in swift flow (where finer sediments are generally rare). Previous field observations on Hydropsyche siltalai pupal cases illustrate that natural local resource limitations of the preferred grain fraction (2.5–3.15 mm) produced chained effects across other grain fractions, as the alternative use of more grains in the 1.6–2 mm fraction (an unlimited resource) induced an increased use of more grains in the 0.315–0.5 mm fraction (another unlimited resource). To examine the implications of these observations for H. siltalai, we used (1) mesocosms to created minor deviations in the availability of the natural grain size composition of the building material of pupal cases at otherwise carefully replicated natural stream habitat conditions and (2) recently developed technologies to assess many case characteristics so far ignored in studies of caddisfly cases. When the preferred coarser grains (2.5–3.15 mm) were unavailable, more grains with intermediate size (1.25–2.0 mm) were used (and not other, still available coarse grains) and fewer larvae built cases in groups, thereby not only loosing the benefits (lower costs for grain transport and silk) but also avoiding potential disadvantages associated with grouped cases (more aggressive encounters with conspecifics for rare building material, less flow exposure and thus reduced water renewal in the pupal chamber). Unavailability of 2.5–3.15-mm and 0.315–0.5-mm grains caused a reduction of larvae building in groups, more use of grains with intermediate size, changes of several other grain characteristics (e.g. number, circularity) and considerable investment into silk to maintain the case resistance. Finally, grain availability deviating most from that observed in nature (no grains of 2.5–3.15 mm and 1.6–2.0 mm) caused dramatic responses, as mortality increased so that fewer pupal cases were built, using typically more coarse grains so that many cases had an elevated resistance against crushing forces; in addition, many males had a retarded development, whereas female development was unaffected. Thus, the response of H. siltalai to any of the three types of grain limitations differed, illustrating an immense diversity to respond to grain-size shortage.     

Verification of ARMA identification for modelling temporal correlations of GNSS observations using the ARMASA toolbox

The classical least-squares (LS) algorithm is widely applied in practice of processing observations from Global Satellite Navigation Systems (GNSS). However, this approach provides reliable estimates of unknown parameters and realistic accuracy measures only if both the functional and stochastic models are appropriately specified. One essential deficiency of the stochastic model implemented in many available GNSS software products consists in neglecting temporal correlations of GNSS observations. Analysing time series of observation residuals resulting from the LS evaluation, the temporal correlation behaviour of GNSS measurements can be efficiently described by means of socalled autoregressive moving average (ARMA) processes. For a given noise realisation, a well-fitting ARMA model can be automatically estimated and identified using the ARMASA toolbox available free of charge in MATLAB® Central.In the preliminary stage of applying the ARMASA toolbox to residual-based modelling of temporal correlations of GNSS observations, this paper presents an empirical performance analysis of the automatic ARMA estimation tool using a large amount of simulated noise time series with representative temporal correlation properties comparable to the GNSS residuals. The results show that the rate of unbiased model estimates increases with data length and decreases with model complexity. For large samples, more than 80% of the identified ARMA models are unbiased. Additionally, the model error representing the deviation between the true data-generating process and the model estimate converges rapidly to the associated asymptotical value for a sufficiently large sample size with respect to the correlation length.     

Cold seeps in Monterey Bay, California: Geochemistry of pore waters and relationship to benthic foraminiferal calcite

An extensive geochemical and biogeochemical examination of CH₄ seeps in the Clam Flats area of Monterey Bay provides insight into the character of relationships between seep geochemistry and benthic foraminiferal geochemistry. The area is characterized by sulfide-rich fluids. Sulfide increases are associated with large increases in alkalinity, as well as small decreases in dissolved Ca and Mg. In addition, only small increases in NH₄ are observed, but values of δ¹³C of dissolved inorganic C are as low as −60‰ at shallow depths (<3 cm). These observations indicate that all these processes are related to the bacterial oxidation of CH₄, which is transported upward by slow seepage of pore fluids. The geochemistry of the pore fluids should be relevant to the geochemistry of the carbonate tests of living and dead foraminifera. However, a profound disequilibrium of approximately an order of magnitude occurs between the δ¹³C values of stained (cytoplasm-containing) foraminiferal carbonate and the C isotope values of ambient pore water dissolved inorganic C. Reasons are unclear for this isotopic disequilibrium, but have important implications for interpretations of foraminiferal carbonate as a paleoenvironmental proxy. Much fine scale work is needed to fully understand the relationships between the biogeochemistry of benthic foraminifera and the geochemistry of the pore waters where they live.     

Improving the representation of major landforms in analytical relief shading

Relief shading is the most common type of cartographic relief representation for print and digital maps. Manual relief shading results in informative and visually pleasing representations of terrain, but it is time consuming and expensive to produce. Current analytical relief shading can be created quickly, but the resulting maps are not as aesthetically appealing and do not show landscape features in an explicit manner. This article introduces an automated digital method that produces shaded relief with locally adjusted illumination directions to simulate the techniques and cartographic principles of manual relief shading. Ridgelines and valley lines are derived from a digital terrain model, vectorized, and used in a diffusion curve algorithm. A graph analysis generalizes the lines before using them for diffusion curve shading. The direction of illumination is adjusted based on the spatial orientation of ridgelines and valley lines. The diffusion curve shading is combined with standard analytical relief shading to create a final diffusion relief shading image. Similar to manual relief shading, major landforms and the structure of the terrain are more clearly shown in the diffusion relief shading. The presented method best highlights major landforms in terrain characterized by sharp, clearly defined ridges and valleys.     

A method to detect and characterize sub‐daily flow fluctuations

Hydro peaking causes an important environmental impact on running water ecosystems. Many affected rivers have a poor ecological status. In rivers affected by hydro peaking, the flow conditions are highly complex and difficult to grasp. To develop a general framework for detecting and characterizing sub‐daily flow fluctuations, we analysed more than 500 Austrian hydrographs, covering the whole range from unimpacted to heavily impacted rivers. Different fluctuation types could be identified according to the potential source: e.g. sub‐daily flow fluctuations caused by hydro peaking, rainfall or snow and glacier melt. Additionally, the term ‘hydro fibrillation’ was established, to indicate frequently occurring artificial fluctuations with comparably low intensities. An automatic procedure was used to detect frequency and intensity of each flow fluctuation. Using variables based on duration curves of flow fluctuation rates (ramping rates), amplitudes, flow ratios, durations and daily numbers of fluctuations, a predictive model (linear discriminant analysis) was fitted to classify hydrographs into predominant fluctuation regimes. This is the basis for a detailed investigation of present sub‐daily flow regimes and to analyse the differences between the regimes. Based on the results, we finally propose a general framework that enables a standardized assessment of flow fluctuations regarding event intensities and/or event timing. The proposed framework offers a standardized selection of particular flow fluctuations referring to increase and decrease events separately. The selection of specific flow fluctuations can be defined with respect to several research questions (e.g. ecologically relevant fluctuations), which offers a wide range of applications. Copyright © 2015 John Wiley & Sons, Ltd.