Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Occurence and thickness of flood layers in varved sediments of Lake Ammersee

Microfacies analyses and X-ray fluorescence scanning (µ-XRF) at sub-mm resolution were conducted on the varved Mid- to Late Holocene interval of two sediment profiles from pre-alpine Lake Ammersee (southern Germany). The coring sites are located in a proximal (AS10prox) and distal (AS10dist) position towards the main tributary River Ammer, in 1.8 km distance from each other. To shed light on sediment distribution within the lake, particular emphasis was (1) the detection of intercalated detrital layers and their micro-sedimentological features, and (2) intra-basin correlation of these event deposits. Detrital layers were dated by microscopic varve counting, verified by accelerator mass spectrometry 14C dating of terrestrial plant macrofossils. Since ~5500 varve years (vyr) BP, in total 1573 detrital layers were detected in either one or both of the investigated sediment profiles. Based on their microfacies, geochemistry, and proximal-distal deposition pattern, detrital layers were interpreted as River Ammer flood deposits. Earlier studies on flood layer seasonality have proven that flood layer deposition occurs predominantly during spring and summer, the flood season at Lake Ammersee. Most prominent features of the record are the onset of regular flood layer deposition at ~5500 vyr BP in AS10prox and ~ 2800 vyr BP in AS10dist as well as three major increases in mean flood layer thickness at ~5500, 2800, and 400 vyr BP. Integrating information from both sediment profiles allowed to interpret these changes in terms of shifts towards higher mean flood intensity. Proposed triggering mechanisms are gradual reduction in Northern Hemisphere orbital summer forcing and superimposed centennial-scale solar activity minima. Likely responses to this forcing are enhanced equator-to-pole temperature gradients and changes in synoptic-scale atmospheric circulation. The consequences for the Ammersee region are more intense cyclones leading to extremer rainfall and flood events in spring and summer.