The third glacier inventory of North Tyrol, Austria, for 2006

A new, third, glacier inventory (GI3) is presented for North Tyrol, which is based on Airborne Laser Scanning for the year 2006. 447 glaciers covering 235 km**2 are included. 4 small, formerly neglected glaciers have been added. Basic quantities such as size, maximum, minimum and median elevation show large variances. Very recent glacier changes between the former inventory (GI2: 1998) and GI3 show a strong reduction in area (-8%) and mean thickness (-7 m). An asymmetry of mean maximum, minimum and median elevation is quantified with approximately 200 m higher values for south-exposed glaciers. Rates of changes are around 1% per year and 1 m per year between GI2 and GI3. The strongest volume losses occurred for glaciers between 5 and 10 km**2.

Digital elevation model (DEM) of the ice sheet in western Dronning Maud Land, Antarctica

In this paper, a new digital elevation model (DEM) is derived for the ice sheet in western Dronning Maud Land, Antarctica. It is based on differential interferometric synthetic aperture radar (SAR) from the European Remote Sensing 1/2 (ERS-1/2) satellites, in combination with ICESat's Geoscience Laser Altimeter System (GLAS). A DEM mosaic is compiled out of 116 scenes from the ERS-1 ice phase in 1994 and the ERS-1/2 tandem mission between 1996 and 1997 with the GLAS data acquired in 2003 that served as ground control. Using three different SAR processors, uncertainties in phase stability and baseline model, resulting in height errors of up to 20 m, are exemplified. Atmospheric influences at the same order of magnitude are demonstrated, and corresponding scenes are excluded. For validation of the DEM mosaic, covering an area of about 130,000 km**2 on a 50-m grid, independent ICESat heights (2004-2007), ground-based kinematic GPS (2005), and airborne laser scanner data (ALS, 2007) are used. Excluding small areas with low phase coherence, the DEM differs in mean and standard deviation by 0.5 +/- 10.1, 1.1 +/- 6.4, and 3.1 +/- 4.0 m from ICESat, GPS, and ALS, respectively. The excluded data points may deviate by more than 50 m. In order to suppress the spatially variable noise below a 5-m threshold, 18% of the DEM area is selectively averaged to a final product at varying horizontal spatial resolution. Apart from mountainous areas, the new DEM outperforms other currently available DEMs and may serve as a benchmark for future elevation models such as from the TanDEM-X mission to spatially monitor ice sheet elevation.

Crete fault scarp point cloud data

This paper assesses the along strike variation of active bedrock fault scarps using long range terrestrial laser scanning (t-LiDAR) data in order to determine the distribution behaviour of scarp height and the subsequently calculate long term throw-rates. Five faults on Cretewhich display spectacular limestone fault scarps have been studied using high resolution digital elevation model (HRDEM) data. We scanned several hundred square metres of the fault system including the footwall, fault scarp and hanging wall of the investigated fault segment. The vertical displacement and the dip of the scarp were extracted every metre along the strike of the detected fault segment based on the processed HRDEM. The scarp variability was analysed by using statistical and morphological methods. The analysis was done in a geographical information system (GIS) environment. Results show a normal distribution for the scanned fault scarp's vertical displacement. Based on these facts, the mean value of height was chosen to define the authentic vertical displacement. Consequently the scarp can be divided into above, below and within the range of mean (within one standard deviation) and quantify the modifications of vertical displacement. Therefore, the fault segment can be subdivided into areas which are influenced by external modification like erosion and sedimentation processes. Moreover, to describe and measure the variability of vertical displacement along strike the fault, the semi-variance was calculated with the variogram method. This method is used to determine how much influence the external processes have had on the vertical displacement. By combining of morphological and statistical results, the fault can be subdivided into areas with high external influences and areas with authentic fault scarps, which have little or no external influences. This subdivision is necessary for long term throw-rate calculations, because without this differentiation the calculated rates would be misleading and the activity of a fault would be incorrectly assessed with significant implications for seismic hazard assessment since fault slip rate data govern the earthquake recurrence. Furthermore, by using this workflow areas with minimal external influences can be determined, not only for throw-rate calculations, but also for determining samples sites for absolute dating techniques such as cosmogenic nuclide dating. The main outcomes of this study include: i) there is no direct correlation between the fault's mean vertical displacement and dip (R² less than 0.31); ii) without subdividing the scanned scarp into areas with differing amounts of external influences, the along strike variability of vertical displacement is ±35%; iii) when the scanned scarp is subdivided the variation of the vertical displacement of the authentic scarp (exposed by earthquakes only) is in a range of ±6% (the varies depending on the fault from 7 to 12%); iv) the calculation of the long term throw-rate (since 13 ka) for four scarps in Crete using the authentic vertical displacement is 0.35 ± 0.04 mm/yr at Kastelli 1, 0.31 ± 0.01 mm/yr at Kastelli 2, 0.85 ± 0.06 mm/yr at the Asomatos fault (Sellia) and 0.55 ± 0.05 mm/yr at the Lastros fault.

Ice thickness distribution and glacier bed of Hallstätter Gletscher

Hallstätter Glacier is the northernmost glacier of Austria. Appendant to the northern Limestone Alps, the glacier is located at 47°28'50'' N, 13°36'50'' E in the Dachstein-region. At the same time with its advance linked to the Little Ice Age (LIA), research on changes in size and mass of Hallstätter glacier was started in 1842 by Friedrich Simony. He observed and documented the glacier retreat related to its last maximum extension in 1856. In addition, Hallstätter Glacier is a subject to scientific research to date. In this thesis methods and results of ongoing mass balance measurements are presented and compared to long term volume changes and meteorological observations. The current mass balance monitoring programm using the direct glaciological method was started 2006. In this context, 2009 the ice thickness was measured with ground penetrating radar. The result are used with digital elevation models reconstucted from historical maps and recent digital elevation models to calculate changes in shape and volume of Hallstätter Glacier. Based on current meteorological measurements near the glacier and longtime homogenized climate data provided by HISTALP, time series of precipitation and temperature beginning at the LIA are produced. These monthly precipitation and monthly mean temperature data are used to compare results of a simple degree day model with the volume change calculated from the difference of the digital elevation models. The two years of direct mass balance measurements are used to calibrate the degree day model. A number of possible future scenarios are produced to indicate prospective changes. Within the 150-year-period between 1856 and 2007 the Hallstätter Glacier lost 1940 meters of its length and 2.23 km**2 in area. 37% of the initial volume of 1856 remained. This retreat came along with a change in climate. The application of a running avarage of 30 years shows an increase in precipitation of 18.5% and a warming of 1.3°C near the glacier between 1866 and 1993. The mass loss was continued in the hydrological years 2006/2007 and 2007/2008 showing mean specific mass balance of -376 mm and -700 mm, respectively. Applying a temperature correction for the different minimum elevations of the glacier, the degree day approach based on the two measured mass balances can reproduce sign and order of magnitude of the volume change of Hallstätter Glacier since 1856. Nevertheless, the relative deviation is significant. Future scenarios show, that 30% of the entire glacier volume remains after subtracting the elevation changes between the digital elevation models of 2002 and 2007 ten times from the surface of 2007. The past and present mass changes of Hallstätter Glacier are showing a retreating glacier as a consequence of rising temperatures. Due to high precepitation, increased with previous warming, the Hallstätter Glacier can and will exist in lower elevation compared to inner alpine glaciers.

(Table 1) Summary statistics of in situ sea-ice characteristics of Ice Station Belgica sites Fabra, Patria/Liège and Brussels in October 2007

Three ice type regimes at Ice Station Belgica (ISB), during the 2007 International Polar Year SIMBA (Sea Ice Mass Balance in Antarctica) expedition, were characterized and assessed for elevation, snow depth, ice freeboard and thickness. Analyses of the probability distribution functions showed great potential for satellite-based altimetry for estimating ice thickness. In question is the required altimeter sampling density for reasonably accurate estimation of snow surface elevation given inherent spatial averaging. This study assesses an effort to determine the number of laser altimeter 'hits' of the ISB floe, as a representative Antarctic floe of mixed first- and multi-year ice types, for the purpose of statistically recreating the in situ-determined ice-thickness and snow depth distribution based on the fractional coverage of each ice type. Estimates of the fractional coverage and spatial distribution of the ice types, referred to as ice 'towns', for the 5 km**2 floe were assessed by in situ mapping and photo-visual documentation. Simulated ICESat altimeter tracks, with spot size ~70 m and spacing ~170 m, sampled the floe's towns, generating a buoyancy-derived ice thickness distribution. 115 altimeter hits were required to statistically recreate the regional thickness mean and distribution for a three-town assemblage of mixed first- and multi-year ice, and 85 hits for a two-town assemblage of first-year ice only: equivalent to 19.5 and 14.5 km respectively of continuous altimeter track over a floe region of similar structure. Results have significant implications toward model development of sea-ice sampling performance of the ICESat laser altimeter record as well as maximizing sampling characteristics of satellite/airborne laser and radar altimetry missions for sea-ice thickness.

Topographic model of Greenland (elevation lines in PanMap and ASCII format)

In this paper, a new high-resolution elevation model of Greenland, including the ice sheet as well as the ice free regions, is presented. It is the first published full coverage model, computed with an average resolution of 2 km and providing an unprecedented degree of detail. The topography is modeled from a wide selection of data sources, including satellite radar altimetry from Geosat and ERS 1, airborne radar altimetry and airborne laser altimetry over the ice sheet, and photogrammetric and manual map scannings in the ice free region. The ice sheet model accuracy is evaluated by omitting airborne laser data from the analysis and treating them as ground truth observations. The mean accuracy of the ice sheet elevations is estimated to be 12-13 m, and it is found that on surfaces of a slope between 0.2° and 0.8°, corresponding to approximately 50% of the ice sheet, the model presents a 40% improvement over models based on satellite altimetry alone. On coastal bedrock, the model is compared with stereo triangulated reference points, and it is found that the model accuracy is of the order of 25-35 m in areas covered by stereo photogrammetry scannings and between 200 and 250 m elsewhere.

Digital surface model, hillshade and orthophoto of the world's largest fossil oyster reef, links to GeoTIFFs

The world's largest fossil oyster reef, formed by the giant oyster Crassostrea gryphoides and located in Stetten (north of Vienna, Austria) is studied by Harzhauser et al., 2015, 2016; Djuricic et al., 2016. Digital documentation of the unique geological site is provided by terrestrial laser scanning (TLS) at the millimeter scale. Obtaining meaningful results is not merely a matter of data acquisition with a suitable device; it requires proper planning, data management, and postprocessing. Terrestrial laser scanning technology has a high potential for providing precise 3D mapping that serves as the basis for automatic object detection in different scenarios; however, it faces challenges in the presence of large amounts of data and the irregular geometry of an oyster reef. We provide a detailed description of the techniques and strategy used for data collection and processing in Djuricic et al., 2016. The use of laser scanning provided the ability to measure surface points of 46,840 (estimated) shells. They are up to 60-cm-long oyster specimens, and their surfaces are modeled with a high accuracy of 1 mm. In addition to laser scanning measurements, more than 300 photographs were captured, and an orthophoto mosaic was generated with a ground sampling distance (GSD) of 0.5 mm. This high-resolution 3D information and the photographic texture serve as the basis for ongoing and future geological and paleontological analyses. Moreover, they provide unprecedented documentation for conservation issues at a unique natural heritage site.

Raw and normalized X-ray fluorescence (XRF) scannings of IODP Expedition 321, Site U1338

We used X-ray fluorescence (XRF) scanning on Site U1338 sediments from Integrated Ocean Drilling Program Expedition 321 to measure sediment geochemical compositions at 2.5 cm resolution for the 450 m of the Site U1338 spliced sediment column. This spatial resolution is equivalent to ~2 k.y. age sampling in the 0-5 Ma section and ~1 k.y. resolution from 5 to 17 Ma. Here we report the data and describe data acquisition conditions to measure Al, Si, K, Ca, Ti, Fe, Mn, and Ba in the solid phase. We also describe a method to convert the data from volume-based raw XRF scan data to a normalized mass measurement ready for calibration by other geochemical methods. Both the raw and normalized data are reported along the Site U1338 splice.