Static GPS measurements and elevation model in the vicinity of the EDML deep-drilling site

Interpretation of ice-core records requires accurate knowledge of the past and present surface topography and stress-strain fields. The European Project for Ice Coring in Antarctica (EPICA) drilling site (0.0684° E and 75.0025° S, 2891.7 m) in Dronning Maud Land, Antarctica, is located in the immediate vicinity of a transient and splitting ice divide. A digital elevation model is determined from the combination of kinematic GPS measurements with the GLAS12 data sets from the ICESat satellite. Based on a network of stakes, surveyed with static GPS, the velocity field around the EDML drilling site is calculated. The annual mean velocity magnitude of 12 survey points amounts to 0.74 m/a. Flow directions mainly vary according to their distance from the ice divide. Surface strain rates are determined from a pentagon-shaped stake network with one center point, close to the drilling site. The strain field is characterised by along flow compression, lateral dilatation, and vertical layer thinning.

Static GPS measurements in the vicinity of EDML and an elevation model for Dronning Maud Land

The knowledge of ice sheet surface topography and the location of the ice divides are essential for ice dynamic modeling. An improved digital elevation model (DEM) of Dronning Maud Land (DML), Antarctica, is presented in this paper. It is based on ground-based kinematic GPS profiles, airborne radar altimetry, and data of the airborne radio-echo sounding system, as well as spaceborne laser altimetry from NASA's Ice, Cloud and land Elevation Satellite (ICESat). The accuracy of ICESat ice sheet altimetry data in the area of investigation is discussed. The location of the ice divides is derived from aspect calculation of the topography and is verified with several velocity data derived from repeated static GPS measurements.

A 20 m spatial resolution seamless multisource Digital Elevation/Depth Model for Lizard Island, northern Great Barrier Reef

Hydrographers have traditionally referred to the nearshore area as the "white ribbon" area due to the challenges associated with the collection of elevation data in this highly dynamic transitional zone between terrestrial and marine environments. Accordingly, available information in this zone is typically characterised by a range of datasets from disparate sources. In this paper we propose a framework to 'fill' the white ribbon area of a coral reef system by integrating multiple elevation and bathymetric datasets acquired by a suite of remote-sensing technologies into a seamless digital elevation model (DEM). A range of datasets are integrated, including field-collected GPS elevation points, terrestrial and bathymetric LiDAR, single and multibeam bathymetry, nautical chart depths and empirically derived bathymetry estimations from optical remote sensing imagery. The proposed framework ranks data reliability internally, thereby avoiding the requirements to quantify absolute error and results in a high resolution, seamless product. Nested within this approach is an effective spatially explicit technique for improving the accuracy of bathymetry estimates derived empirically from optical satellite imagery through modelling the spatial structure of residuals. The approach was applied to data collected on and around Lizard Island in northern Australia. Collectively, the framework holds promise for filling the white ribbon zone in coastal areas characterised by similar data availability scenarios. The seamless DEM is referenced to the horizontal coordinate system MGA Zone 55 - GDA 1994, mean sea level (MSL) vertical datum and has a spatial resolution of 20 m.

Terrain model of the eastern slope of the Porcupine Seabight (50 m grid size)

The topography of the eastern margin of the Porcupine Seabight was surveyed in June 2000 utilizing swath bathymetry. The survey was carried out during RV Polarstern cruise ANT XVII/4 as part of the GEOMOUND project. The main objective was to map and investigate the seafloor topography of this region. The investigated area contains a variability of morphological features such as deep sea channels and giant mounds. The survey was planned and realized on the basis of existing data so as to guarantee the complete coverage of the margin. In order to achieve a resolution of the final digital terrain model (DTM) that meets the project demands, data processing was adjusted accordingly. The grid spacing of the DTM was set to 50 m and an accuracy better than 1% of the water depth was achieved for 96% of the soundings.

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.

Northeast Greenland - Digital Bathymetric Compilation, with links to DBM in ArcGIS format

A new digital bathymetric model (DBM) for the Northeast Greenland (NEG) continental shelf (74°N - 81°N) is presented. The DBM has a grid cell size of 250 m × 250 m and incorporates bathymetric data from 30 multibeam cruises, more than 20 single-beam cruises and first reflector depths from industrial seismic lines. The new DBM substantially improves the bathymetry compared to older models. The DBM not only allows a better delineation of previously known seafloor morphology but, in addition, reveals the presence of previously unmapped morphological features including glacially derived troughs, fjords, grounding-zone wedges, and lateral moraines. These submarine landforms are used to infer the past extent and ice-flow dynamics of the Greenland Ice Sheet during the last full-glacial period of the Quaternary and subsequent ice retreat across the continental shelf. The DBM reveals cross-shelf bathymetric troughs that may enable the inflow of warm Atlantic water masses across the shelf, driving enhanced basal melting of the marine-terminating outlet glaciers draining the ice sheet to the coast in Northeast Greenland. Knolls, sinks, and hummocky seafloor on the middle shelf are also suggested to be related to salt diapirism. North-south-orientated elongate depressions are identified that probably relate to ice-marginal processes in combination with erosion caused by the East Greenland Current. A single guyot-like peak has been discovered and is interpreted to have been produced during a volcanic event approximately 55 Ma ago.

Present-day King George Island ice cap geometry measured by ground penetrating radar in 2008/09

We present a consistent data set for the ice thickness, the bedrock topography and the ice surface topography of the King George Island ice cap (Arctowski Icefield and the adjacent central part). The data set is composed of groundbased and airborne Ground Penetrating Radar (GPR) and differential GPS (DGPS) measurements, obtained during several field campaigns. The data set incorporates groundbased measurements in the safely accessible inner parts and airborne measurements in the heavily crevassed coastal areas of the ice cap. In particular, the inclusion of airborne GPR measurements with the 30MHz BGR-P30-System developed at the Institute of Geophysics (University of Münster) completes the picture of the ice geometry substantially. The compiled digital elevation model of the bedrock shows a rough, highly variable topography with pronounced valleys, ridges, and troughs. The mean ice thickness is approx. 238m, with a maximum value of approx. 400m in the surveyed area. Noticeable are bounded areas in the bedrock topography below sea level where marine based ice exists.

Downscaled land surface temperatures obtained by a particle smoothing approach, with links to NetCDF files

This dataset contains continuous time series of land surface temperature (LST) at spatial resolution of 300m around the 12 experimental sites of the PAGE21 project (grant agreement number 282700, funded by the EC seventh Framework Program theme FP7-ENV-2011). This dataset was produced from hourly LST time series at 25km scale, retrieved from SSM/I data (André et al., 2015, doi:10.1016/j.rse.2015.01.028) and downscaled to 300m using a dynamic model and a particle smoothing approach. This methodology is based on two main assumptions. First, LST spatial variability is mostly explained by land cover and soil hydric state. Second, LST is unique for a land cover class within the low resolution pixel. Given these hypotheses, this variable can be estimated using a land cover map and a physically based land surface model constrained with observations using a data assimilation process. This methodology described in Mechri et al. (2014, doi:10.1002/2013JD020354) was applied to the ORCHIDEE land surface model (Krinner et al., 2005, doi:10.1029/2003GB002199) to estimate prior values of each land cover class provided by the ESA CCI-Land Cover product (Bontemps et al., 2013) at 300m resolution . The assimilation process (particle smoother) consists in simulating ensemble of LST time series for each land cover class and for a large number of parameter sets. For each parameter set, the resulting temperatures are aggregated considering the grid fraction of each land cover and compared to the coarse observations. Miniminizing the distance between the aggregated model solutions and the observations allow us to select the simulated LST and the corresponding parameter sets which fit the observations most closely. The retained parameter sets are then duplicated and randomly perturbed before simulating the next time window. At the end, the most likely LST of each land cover class are estimated and used to reconstruct LST maps at 300m resolution using ESA CCI-Land Cover. The resulting temperature maps on which ice pixels were masked, are provided at daily time step during the nine-year analysis period (2000-2009).

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Holocene experiment with coupled atmosphere-ocean-model ECHAM5/MPI-OM

Changes in the Earth's orbit lead to changes in the seasonal and meridional distribution of insolation. We quantify the influence of orbitally induced changes on the seasonal temperature cycle in a transient simulation of the last 6000 years - from the mid-Holocene to today - using a coupled atmosphere-ocean general circulation model (ECHAM5/MPI-OM) including a land surface model (JSBACH). The seasonal temperature cycle responds directly to the insolation changes almost everywhere. In the Northern Hemisphere, its amplitude decreases according to an increase in winter insolation and a decrease in summer insolation. In the Southern Hemisphere, the opposite is true. Over the Arctic Ocean, decreasing summer insolation leads to an increase in sea-ice cover. The insulating effect of sea ice between the ocean and the atmosphere leads to decreasing heat flux and favors more "continental" conditions over the Arctic Ocean in winter, resulting in strongly decreasing temperatures. Consequently, there are two competing effects: the direct response to insolation changes and a sea-ice insulation effect. The sea-ice insulation effect is stronger, and thus an increase in the amplitude of the seasonal temperature cycle over the Arctic Ocean occurs. This increase is strongest over the Barents Shelf and influences the temperature response over northern Europe. We compare our modeled seasonal temperatures over Europe to paleo reconstructions. We find better agreements in winter temperatures than in summer temperatures and better agreements in northern Europe than in southern Europe, since the model does not reproduce the southern European Holocene summer cooling inferred from the paleo reconstructions. The temperature reconstructions for northern Europe support the notion of the influence of the sea-ice insulation effect on the evolution of the seasonal temperature cycle.

Grain size distribution and sea level changes in Gåsfjärden fjord

We analysed long-term variations in grain-size distribution in sediments from Gåsfjärden, a fjord-like inlet on the south-west Baltic Sea, and explored potential drivers of the recorded changes in sediment grain-size data. Over the last 5.4 thousand years (ka), the relative sea level decreased 17 m in the study region, caused by isostatic land uplift. As a consequence, Gåsfjärden has been transformed from an open coastal setting into a semi-closed inlet surrounded on the east by numerous small islands. To quantitatively estimate the morphological changes in Gåsfjärden over the last 5.4 ka and to further link the changes to our grain-size data, a digital elevation model (DEM)-based openness index was calculated. In the period between 5.4 and 4.4 ka BP, the inlet was characterised by the largest openness index. During this interval, the highest sand contents (~0.4 %) and silt/clay ratios (~0. 3) in the sediment sequence were recorded, indicating relatively high bottom water energy. After 4.4 ka BP, the average sand content was halved to ~0.2 % and the silt/clay ratios showed a significant decreasing trend over the last 4 ka. These changes are found to be associated with the gradual embayment of Gåsfjärden as represented in the openness index. The silt/clay ratios exhibited a delayed and slower change compared with the sand contents, which further suggest that finer particles are less sensitive to changes in hydrodynamic energy. Our DEM-based coastal openness index has proved to be a useful tool for interpreting the sedimentary grain-size record.

The glacier inventory of Salzburg, Austria, 2007/2009

A third glacier inventory (GI3) is presented for the province of Salzburg where 173 glaciers are located in the seven mountain ranges: Ankogel (47°4'N, 13°14'E), Glockner, Granatspitz, Sonnblick (Goldberg), Hochkönig, Venediger and Zillertal (47°8'N, 12°7'E). The basis for the new GI3 are orthophotos of 2007 and 2009 and the digital elevation model (DEM) of the southern part of Salzburg. On the basis of former inventories, area- and volume changes have been calculated. The biggest relative loss of glacier area per mountain range was found in the Ankogel range and on Hochkönig as a result of the disrupted structure of their small and thin glaciers. In terms of absolute values, the largest changes took place in the Glockner- and Venediger range with an area loss of -10.1 km**2 and -9.7 km**2 during the period between GI1 (1969) and GI3 (2007/2009), respectively. Volume changes have been calculated for nearly half of the glacier area in Salzburg, where DEMs were available. The Glockner, Granatspitz and Sonnblick mountain ranges showed a volume loss of -0.481 km**3 which corresponds to a mean thickness change of -10.5 m. An extrapolation of these changes to all of the 173 glaciers in Salzburg results in a loss of about 1.04 km**3 between GI1 and GI3 and 0.44 km**3 between GI2 and GI3. Overall annual changes in the province of Salzburg between GI2 and GI3 were higher than between GI1 and GI2 and show likewise changes such as those of Tyrol.

Sediment volume variations inferred from DGPS measurements within Bossons glacier proglacial area (Mont-Blanc massif, France)

This dataset presents Differential Global Positioning System data (DGPS) acquired within the Bossons glacier proglacial area. Bossons glacier is a rapidly retreating glacier and its proglacial area is deglaciated for ~30 years. Bossons stream is one of the outlets of the subglacial drainage system. It starts as a 800 m steep cascade reach, then flows through an area with gentler slope : the Plan des Eaux (PdE). PdE is a 300 m long, 50 m wide proglacial alluvial plain with an increasing channel mobility in the downstream direction but decreasing slope gradient and incision. As it may act a sediment trap, studying periglacial and proglacial erosion processes in the Bossons catchment requires to quantify PdE sediment volume evolution. A several meter-sized block located within Bossons proglacial area was set up as GPS base : its location was measured by one antenna (Topcon Hyper Pro) by performing 600 consecutive measurements throughout one day. A second antenna (Topcon Hyper Pro) was then used to measure XYZ location of points in the proglacial area with a ~2 m grid. Radio communication between the two antennas allowed differential calculations to be automatically carried out on field using the Topcon FC-250 hand controller. This methodology yields 3 cm XY and 1.5 cm Z uncertainties. DGPS data have been acquired through 10 campaigns from 2004 to 2014; campaigns from 2004 to 2008 cover a smaller area than those from 2010 to 2014. Digital Elevation Model (DEM) have been interpolated from DGPS data and difference between two DEMs yields deposited and eroded volume within PdE. Maps of PdE volume variation between two campaigns show that incision mainly occurs in the upper and lower sections where as deposition dominates in the middle section. Deposition, denudation and net rate (deposition rate - denudation rate) are calculated by normalizing volumes by DEM areas. Deposition dominates results with a mean net rate of 29 mm/yr. However, strong inter-annual variability exists and some years are dominated by denudation : -36 mm/yr and -100 mm/yr for 2006 and 2011, respectively. Nonetheless, oldest campaigns (2004 to 2008) were carried out on the lower part part of the alluvial plain and ruling them out to keep only complete DEM (2010 to 2014) yields a mean net rate of ~15 mm/yr. This results is coherent with field observations of both strong deposition (e.g. flood deposits) and strong erosion (e.g. 30 cm incision) evidences. Bossons glacier proglacial area is thus dynamic with year-to-year geormorphological changes but may leans toward increasing its mean elevation through a deposition dominated system.