Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Multi-temporal shoreline change rates on Taku Atoll

Recent-past shoreline changes on reef islands are now subject to intensified monitoring via remote sensing data. Based on these data, rates of shoreline change calculated from long-term measurements (decadal) are often markedly lower than recent short-term rates (over a number of years). This observation has raised speculations about the growing influence of sea-level rise on reef island stability. This observation, however, can also be explained if we consider two basic principles of geomorphology and sedimentology. For Takú Atoll, Papua New Guinea, we show that natural shoreline fluctuations of dynamic reef islands have a crucial influence on the calculation of short-term rates of change. We analyze an extensive dataset of multitemporal shoreline change rates from 1943 to 2012 and find that differing rates between long- and short-term measurements consistently reflect the length of the observation interval. This relationship appears independent from the study era and indicates that reef islands were equally dynamic during the early periods of analysis, i.e. before the recent acceleration of sea-level rise. Consequently, we suggest that high rates of shoreline change calculated from recent short-term observations may simply result from a change in temporal scale and a shift from geomorphic equilibrium achieved over cyclic time towards an apparent disequilibrium during shorter periods of graded time. This new interpretation of short- and long-term shoreline change rates has important implications for the ongoing discussion about reef island vulnerability, showing that an observed jump from low to high rates of change may be independent from external influences, including but not limited to sea-level rise.

Automated Tools for Layer Recognition (BMPix) and Counting (PEAK)

We present tools for rapid and quantitative detection of sediment lamination. The BMPix tool extracts color and gray-scale curves from images at pixel resolution. The PEAK tool uses the gray-scale curve and performs, for the first time, fully automated counting of laminae based on three methods. The maximum count algorithm counts every bright peak of a couplet of two laminae (annual resolution) in a smoothed curve. The zero-crossing algorithm counts every positive and negative halfway-passage of the curve through a wide moving average, separating the record into bright and dark intervals (seasonal resolution). The same is true for the frequency truncation method, which uses Fourier transformation to decompose the curve into its frequency components before counting positive and negative passages. We applied the new methods successfully to tree rings, to well-dated and already manually counted marine varves from Saanich Inlet, and to marine laminae from the Antarctic continental margin. In combination with AMS14C dating, we found convincing evidence that laminations in Weddell Sea sites represent varves, deposited continuously over several millennia during the last glacial maximum. The new tools offer several advantages over previous methods. The counting procedures are based on a moving average generated from gray-scale curves instead of manual counting. Hence, results are highly objective and rely on reproducible mathematical criteria. Also, the PEAK tool measures the thickness of each year or season. Since all information required is displayed graphically, interactive optimization of the counting algorithms can be achieved quickly and conveniently.

Multi-temporal mapping of seagrass cover, species and biomass of the Eastern Banks, Moreton Bay, Australia, with links to shapefiles

The spatial and temporal dynamics of seagrasses have been studied from the leaf to patch (100 m**2) scales. However, landscape scale (> 100 km**2) seagrass population dynamics are unresolved in seagrass ecology. Previous remote sensing approaches have lacked the temporal or spatial resolution, or ecologically appropriate mapping, to fully address this issue. This paper presents a robust, semi-automated object-based image analysis approach for mapping dominant seagrass species, percentage cover and above ground biomass using a time series of field data and coincident high spatial resolution satellite imagery. The study area was a 142 km**2 shallow, clear water seagrass habitat (the Eastern Banks, Moreton Bay, Australia). Nine data sets acquired between 2004 and 2013 were used to create seagrass species and percentage cover maps through the integration of seagrass photo transect field data, and atmospherically and geometrically corrected high spatial resolution satellite image data (WorldView-2, IKONOS and Quickbird-2) using an object based image analysis approach. Biomass maps were derived using empirical models trained with in-situ above ground biomass data per seagrass species. Maps and summary plots identified inter- and intra-annual variation of seagrass species composition, percentage cover level and above ground biomass. The methods provide a rigorous approach for field and image data collection and pre-processing, a semi-automated approach to extract seagrass species and cover maps and assess accuracy, and the subsequent empirical modelling of seagrass biomass. The resultant maps provide a fundamental data set for understanding landscape scale seagrass dynamics in a shallow water environment. Our findings provide proof of concept for the use of time-series analysis of remotely sensed seagrass products for use in seagrass ecology and management.

Multi-proxy analyses of sediment core GeoB13801-2

Over the Uruguayan shelf and uppermost slope the coalescence of northward flowing Subantarctic Shelf Water and southward flowing Subtropical Shelf Water forms a distinct thermohaline front termed the Subtropical Shelf Front (STSF). Running in a SW direction diagonally across the shelf from the coastal waters at 32°S towards the shelf break at ca. 36°S, the STSF represents the shelf-ward extension of the Brazil-Malvinas Confluence zone. This study reconstructs latitudinal STSF shifts during the Holocene based on benthic foraminifera d18O and d13C, total organic carbon, carbonate contents, Ti/Ca, and grain-size distribution from a high-accumulation sedimentary record located at an uppermost continental-slope terrace. Our data provide direct evidence for: (1) a southern STSF position (to the South of the core site) at the beginning of the early Holocene (>9.4 cal ka BP) linked to a more southerly position of the Southern Westerly Winds in combination with restricted shelf circulation intensity due to lower sea level; (2) a gradual STSF northward migration (bypassing the core site towards the North) primarily forced by the northward migration of the Southern Westerly Winds from 9.4 cal ka BP onwards; (3) a relatively stable position of the front in the interval between 7.2 and 4.0 cal ka BP; (4) millennial-scale latitudinal oscillations close to 36°S of the STSF after 4.0 cal ka BP probably linked to the intensification in El Niño Southern Oscillation; and (5) a southward migration of the STSF during the last 200 years possibly linked to anthropogenic influences on the atmosphere.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2014S13, deployed during the CryoVEx2014 field campaign

Snow height was measured by the Snow Depth Buoy 2014S13, an autonomous platform, drifting on Arctic sea ice, deployed during the CryoVEx2014 field campaign. The resulting time series describes the evolution of snow height as a function of place and time between 2014-03-30 and 2014-07-20 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on multi year ice. In addition to snow height, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow height occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow height may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow height, because no initial readings of absolute snow height are available.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2013S3, deployed at the Barneo ice camp 2013

Snow height was measured by the Snow Depth Buoy 2013S3, an autonomous platform, drifting on Arctic sea ice. This buoy was deployed at the Barneo ice camp 2013. The resulting time series describes the evolution of snow height as a function of place and time between 2013-04-09 and 2013-06-13 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on multi year ice. In addition to snow height, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow height occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow height may still be used for sea ice drift analyses.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2014S15, deployed during POLARSTERN cruise ARK-XXVIII/4 (PS87)

Snow height was measured by the Snow Depth Buoy 2014S15, an autonomous platform, drifting on Arctic sea ice, deployed during POLARSTERN cruise ARK-XXVIII/4 (PS87). The resulting time series describes the evolution of snow depth as a function of place and time between 2014-08-29 and 2014-12-31 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). The buoy was installed on multi year ice. In addition to snow depth, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Records without any snow depth may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow depth, because no initial readings of absolute snow depth are available.

Oceanographic investiations from the Sea of Ochotsk

The Sea of Okhotsk is a marginal sea of the Pacific Ocean, which is characterized by strong variations in the productivity and sediment supply due to sea ice transport and river input. Furthermore the variations in the hydrological cycle determine the formation of the SOIW (Sea of Okhotsk Intermediate Water) which plays an important role in the ventilation processes in the intermediate water of the N-Pacific. Isotope data measured on planktonic and benthic foraminifera, sedimentological and geochemical studies of sediment cores and surface samples from the Sea of Okhotsk are used to reconstruct the paleoceanography during the past 350.000 years. The dating and correlation of the sediments are based on oxygen isotope stratigraphy, absolute ages, magnetic susceptibility as well as a detailled tephrachronology of the entire basin. The sedimentation rates are characterized by temporal and spatial variations. The maximum sedimentation rate takes place at the continental slope off Sakhalin due to the input of the Amur River, the sea ice drift and the high productivity. The sedimentation rate in the eastern part of the Sea of Okhotsk is generelly high because of the influence of the nutrient-rich Kamchatka Current. In the central and northern parts of the Sea of Okhotsk, areas with low productivity and reduced terrestrial supply, the sedimentation rate is the lowest. The analyses of the surface sediment samples make it possible to characterize the (sub)- recent sediment supply and transportation processes. The bulk sediment measurements, isotope data and the accumulation rate of ice-rafted debris (IRD) show a dominant sea ice cover and a region with a high productivity as well as a high Amur River input in the western part of the sea. The eastern part of the Sea of Okhotsk, however, is marked by the predominance of warm and nutrient-rich water masses coming from the Kamchatka Current which restricts the sea ice cover. This is reflected in low content of ice-rafted debris and high productivity proxies as well as in isotope data. The deposits of the Sea of Okhotsk are characterized by terrestrial, biogenic and volcanogenic sediment input which varies temporally and spatially. Here, the sedimentation pattern is dominated by the terrestrial input. Bulk sediment measurements and sample analyses of the > 63 micron particle input make it possible to distinguish glacial and interglacial fluctuations. The sedimentation processes during glacial times are determined by a high content of ice-rafted debris, whereas the primary production is higher during interglacial periods. During the last glacial/interglacial cycle the IRD-distribution pattern indicates a strong sea ice transport in the western part and in large areas of the open sea in the eastern part of the Sea of Okhotsk with a relatively constant ice-drift system. The IRD flux in sediments of the oxygen isotope Stage 6 reflects a new sedimentation pattern in the eastern part of the sea. This high IRD accumulation rate indicates ice advances beyond the shelf margin and an iceberg transport from NE-E direction into the Sea of Okhotsk. The several large, brief, negative anomalies in d13C values of Neogloboquadrina pachyderma (s) show releases of methane from basin sediments which correspond to periods of relative sea level falls. The high sedimentation rates on the Sakhalin slope allow insights into the climatic history in Holocene and indicate shorter-scale variations oscillation in Stage 3, which correlate with the global climatic changes. These variations are described as Dansgaard-Oeschger cycles in Greenland ice cores and as Heinrich-Events in several marine sediment cores from the N-Atlantic.