Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.

Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Global and Regional Sea level budget components from GRACE and radar altimetry (2002-2014)

This dataset contains the result of a joint least squares inversion of GRACE and altimetry data. The results are evaluated in terms of sea level change for the global ocean as well as dedicated areas. In addition, some auxiliary data is provided to enable reproducibility of the results in Rietbroek et al. 2016, and a google Earth kmz file is provided which visualizes the trends derived from the inversion results.

Static GPS measurements and elevation model in the hinterland of Neumayer

As part of the CryoSat Cal/Val activities and the pre-site survey for an ice core drilling contributing to the International Partnerships in Ice Core Sciences (IPICS), ground based kinematic GPS measurements were conducted in early 2007 in the vicinity of the German overwintering station Neumayer (8.25° W and 70.65° S). The investigated area comprises the regions of the ice tongues Halvfarryggen and Søråsen, which rise from the Ekströmisen to a maximum of about 760 m surface elevation, and have an areal extent of about 100 km x 50 km each. Available digital elevation models (DEMs) from radar altimetry and the Antarctic Digital Database show elevation differences of up to hundreds of meters in this region, which necessitated an accurate survey of the conditions on-site. An improved DEM of the Ekströmisen surroundings is derived by a combination of highly accurate ground based GPS measurements, satellite derived laser altimetry data (ICESat), airborne radar altimetry (ARA), and radio echo sounding (RES). The DEM presented here achieves a vertical accuracy of about 1.3 m and can be used for improved ice dynamic modeling and mass balance studies.

Updated bathymetric chart of the East Aral Sea

The Aral Sea is located in an arid region with much sand and salt deposits in the surrounding barren open land. Hence, significant displacements of sediments into the lakebed by the action of wind, water, gravity, or snow are likely. The bathymetry of the lake was last observed in the 1960s, and within the last half century, the structure of the lakebed has changed. Based on satellite observations of the temporal changes of shoreline (Landsat optical remote sensing) and water level (multi-mission satellite altimetry) over more than one decade an updated bathymetric chart for the East Basin of the Aral Sea has been generated. During this time, the geometry of the shallow East Basin experienced strong fluctuations due to the occurrence of periods of drying and strong inflow. By the year 2014 the East Basin fell dry. The dynamic behavior of the basin allowed for estimating the lake's bathymetry from a series of satellite-based information. The river mouth made its impression in the present East Aral Sea, because its shrinking led to the inflow of much sediment into the lake's interior. In addition, salt deposits along the shorelines increased the corresponding elevation, a phenomenon that is more pronounced in the reduced lakebed because of increased salinity. It must be noted that height estimates from satellite altimetry were only possible down to a minimum elevation of 27 m above sea level due to a lack of reliable altimetry data over the largely reduced water surface. In order to construct a complete bathymetric chart of the lakebed of the East Aral Sea heights below 27 m were obtained solely from Landsat optical images following the SRB (Selected Region Boundary) approach as described by Singh et al. (2015).

The CMCC Eddy-permitting Global Ocean Physical Reanalysis (C-GLORS v5, 1980-2014), links to NetCDF files

The CMCC Global Ocean Physical Reanalysis System (C-GLORS) is used to simulate the state of the ocean in the last decades. It consists of a variational data assimilation system (OceanVar), capable of assimilating all in-situ observations along with altimetry data, and a forecast step performed by the ocean model NEMO coupled with the LIM2 sea-ice model. KEY STRENGTHS: - Data are available for a large number of ocean parameters - An extensive validation has been conducted and is freely available - The reanalysis is performed at high resolution (1/4 degree) and spans the last 30 years KEY LIMITATIONS: - Quality may be discontinuos and depend on observation coverage - Uncertainty estimates are simply derived through verification skill scores

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

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.

Improved digital elevation model (DEM) for Lake Eyre, Australia

Here we demonstrate the applicability of using altimetry data and Landsat imagery to provide the most accurate digital elevation model (DEM) of Australia's largest playa lake - Lake Eyre. We demonstrate through the use of geospatial techniques a robust assessment of lake area and volume of recent lake-filling episodes whilst also providing the most accurate estimates of area and volume for larger lake filling episodes that occurred throughout the last glacial cycle. We highlight that at a depth of 25 m Lake Mega-Eyre would merge with the adjacent Lake Mega-Frome to form an immense waterbody with a combined area of almost 35,000 km**2 and a combined volume of ~520 km**3. This would represent a vast water body in what is now the arid interior of the Australian continent. The improved DEM is more reliable from a geomorphological and hydrological perspective and allows a more accurate assessment of water balance under the modern hydrological regime. The results presented using GLAS/ICESat data suggest that earlier historical soundings were correct and the actual lowest topographic point in Australia is -15.6 m below sea level. The results also contrast nicely the different basin characteristics of two adjacent lake systems; Lake Eyre and Lake Frome.

(Table 1) Lake elevation change derived from 4 or more years of measurements

In this study, ICESat altimetry data are used to provide precise lake elevations of the Tibetan Plateau (TP) during the period of 2003-2009. Among the 261 lakes examined ICESat data are available on 111 lakes: 74 lakes with ICESat footprints for 4-7 years and 37 lakes with footprints for 1 -3 years. This is the first time that precise lake elevation data are provided for the 111 lakes. Those ICESat elevation data can be used as baselines for future changes in lake levels as well as for changes during the 2003-2009 period. It is found that in the 74 lakes (56 salt lakes) examined, 62 (i.e. 84%) of all lakes and 50 (i.e. 89%) of the salt lakes show tendency of lake level increase. The mean lake water level increase rate is 0.23 m/year for the 56 salt lakes and 0.27 m/year for the 50 salt lakes of water level increase. The largest lake level increase rate (0.80 m/year) found in this study is the lake Cedo Caka. The 74 lakes are grouped into four subareas based on geographical locations and change tendencies in lake levels. Three of the four subareas show increased lake levels. The mean lake level change rates for subareas I, II, III, IV, and the entire TP are 0.12, 0.26, 0.19, -0.11, and 0.2 m/year, respectively. These recent increases in lake level, particularly for a high percentage of salt lakes, supports accelerated glacier melting due to global warming as the most likely cause.

Geochemical and isotopic composition of basalts from the axial Mid-Atlantic Ridge, South Atlantic

One of the essential problems of oceanic tectonics is estimation of the influence of plumes of the deep hot mantle on processes in the axial spreading zone. Areas of two giant (St. Helena and Tristan da Cunha) plumes in the Mid-Atlantic Ridge (MAR) rift zone (South Atlantic) are characterized by the effusion of basalts that differ from typical depleted riftogenic tholeiites by anomalously high contents of lithophile components and specific isotopic compositions. Moreover, the rift valley floor with basalt effusion is significantly uplifted above the adjacent sectors of the rift. The formation of the St. Helena Seamount located 400 km east of the MAR axis is related to magmatism that is active to this day. St. Helena Island is a member of the structural ensemble of large volcanic seamounts (Bonaparte, Bagration, and Kutuzov). Like St. Helena Island, each seamount incorporates a series of smaller rises of different morphologies and dimensions. Thus, a system of subparallel series of NE-trending (~45°) rises extend from the seamount ensemble to the African continent. According to the plate tectonics concept, the seamount series represent hotspots related to a deep mantle plume that can be projected onto the present-day St. Helena Island area (St. Helena plume). At the same time, the inferred topographic map based on satellite altimetry data shows that the seamount series also extend along the opposite southwestern direction (~225°) toward the axial MAR and even intersect the latter structure. This fact cannot be explained by the hotspot hypothesis, which suggests stationary positions of plumes relative to the mobile oceanic plate. In the course of Cruise 10 of the R/V Akademik Ioffe (2002), detailed geological and geophysical investigations were carried out at the junction of one structural series with the MAR rift zone located near the Martin Vaz Fracture Zone (Martin Vaz test area, 19°-20° S). The present communication is devoted to the study of lithology, geochemistry, and isotopy of basalts dredged at the test area.

Hydrography and biogeochemistry during S.A. Agulhas II cruise Agu011, in the SW Indian Ocean, April 2014

Several open-ocean mesoscale features, a "young" warm-core (or anti-cyclonic) eddy at 52°S, an "older" warm-core eddy at 57.5°S, as well as an adjacent cold-core (or cyclonic) eddy at 56°S, were surveyed during a R/V S.A. Agulhas II cruise in April 2014. The main aim of the survey was to obtain hydrographical and biogeochemical profile data for contrasting open-ocean eddies in the Southern Ocean, that will be suitable for comparison and modelling of their heat, salt and nutrient characteristics, and the changes that occur in these properties as warm-core eddies migrate from the polar front southwards into the Southern Ocean. A total of 18 CTD stations were occupied in a sector south of the South-West Indian Ridge, along three transects crossing several mesoscale features identified from satellite altimetry data prior to the cruise.