EOT11a - Global Empirical Ocean Tide model from multi-mission satellite altimetry, with links to model results

EOT11a is a global (E)mpirical (O)cean (T)ide model derived in 2011 by residual analysis of multi-mission satellite (a)ltimeter data. EOT11a includes amplitudes and phases of the main astronomical tides M2, S2, N2, K2, 2N2, O1, K1, P2, and Q1, the non-linear constituent M4, the long period tides Mm and Mf, and the radiational tide S1. Ocean tides as well as loading tides are provided. EOT11a was computed by means of residual tidal analysis of multi-mission altimeter data from TOPEX/Poseidon, ERS-2, ENVISAT, and Jason-1/2, as far as acquired between September 1992 and April 2010. The resolution of 7.5'x7.5' is identical with FES2004 which was used as reference model for the residual tide analysis. The development of EOT11a was funded by the Deutsche Forschungsgemeinschaft (DFG) under grant BO1228/6-2.

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.

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.

Mesoscale ocean eddy dataset in the Lofoten Basin from Satellite Altimetry

The Lofoten Basin is the most eddy rich region in the Norwegian Sea. In this paper, the characteristics of these eddies are investigated from a comprehensive database of nearly two decades of satellite altimeter data (1995-2013) together with Argo profiling floats and surface drifter data. An automated method identified 1695/1666 individual anticyclonic/cyclonic eddies in the Lofoten Basin from more than 10,000 altimeter-based eddy observations. The eddies are found to be predominantly generated and residing locally. The spatial distributions of lifetime, occurrence, generation sites, size, intensity, and drift of the eddies are studied in detail. The anticyclonic eddies in the Lofoten Basin are the most long-lived eddies (>60 days), especially in the western part of the basin. We reveal two hotspots of eddy occurrence on either side of the Lofoten Basin. Furthermore, we infer a cyclonic drift of eddies in the western Lofoten Basin. Barotropic energy conversion rates reveals energy transfer from the slope current to the eddies during winter. An automated colocation of surface drifters trapped inside the altimeter-based eddies are used to corroborate the orbital speed of the anticyclonic and cyclonic eddies. Moreover, the vertical structure of the altimeter-based eddies is examined using colocated Argo profiling float profiles. Combination of altimetry, Argo floats, and surface drifter data is therefore considered to be a promising observation-based approach for further studies of the role of eddies in transport of heat and biomass from the slope current to the Lofoten Basin.

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.

Composition of sulfides from the Sierra Leone Fracture Zone (Central Atlantic Ocean)

An analysis of data on the location of hydrothermal fields, seismicity, and satellite altimetry evidences that in mid-ocean ridges with low spreading rate hydrothermal fields tend to be grouped in areas with generally low seismic activity and at intersections of discontinuities and rift zones. Based on this assumption, the Sierra Leone Fracture Zone was studied in 2000 during Cruise 22 of R/V Akademik Nikolaj Strakhov. A study of gabbrodolerite and dolerite showed that sulfide ore minerals in them were formed both by hydrothermal and magmatic processes. An analysis of melt inclusions demonstrated that magmatic complexes formed from a high-temperature (1210-1255°C) low-potassium melt of the N-MORB type. Investigations of fluid inclusions revealed that gabbro and dolerite formed under influence of an active hydrothermal system at temperature 205-226°C. Thus, the Sierra Leone Fracture Zone is considered to be perspective for a discovery of a new hydrothermal field.

Crustal uplift rates of 4 GPS stations near Jakobshavn Isbrae (2006-2009)

We analyze 2006-2009 data from four continuous Global Positioning System (GPS) receivers located between 5 and 150 km from the glacier Jakobshavn Isbrae, West Greenland. The GPS stations were established on bedrock to determine the vertical crustal motion due to the unloading of ice from Jakobshavn Isbrae. All stations experienced uplift, but the uplift rate at Kangia North, only 5 km from the glacier front, was about 10 mm/yr larger than the rate at Ilulissat, located only ~45 km further away. This suggests that most of the uplift is due to the unloading of the Earth's surface as Jakobshavn thins and loses mass. Our estimate of Jakobshavn's contribution to uplift rates at Kangia North and Ilulissat are 14.6 ± 1.7 mm/yr and 4.9 ± 1.1 mm/yr, respectively. The observed rates are consistent with a glacier thinning model based on repeat altimeter surveys from NASA's Airborne Topographic Mapper (ATM), which shows that Jakobshavn lost mass at an average rate of 22 ± 2 km**3/yr between 2006 and 2009. At Kangia North and Ilulissat, the predicted uplift rates computed using thinning estimates from the ATM laser altimetry are 12.1 ± 0.9 mm/yr and 3.2 ± 0.3 mm/yr, respectively. The observed rates are slightly larger than the predicted rates. The fact that the GPS uplift rates are much larger closer to Jakobshavn than further away, and are consistent with rates inferred using the ATM-based glacier thinning model, shows that GPS measurements of crustal motion are a potentially useful method for assessing ice-mass change models.

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).

(Table 3) Details of field-based ice and snow profile measurements from the SIMBA 2007 experiment

Although sea-ice extent in the Bellingshausen-Amundsen (BA) seas sector of the Antarctic has shown significant decline over several decades, there is not enough data to draw any conclusion on sea-ice thickness and its change for the BA sector, or for the entire Southern Ocean. This paper presents our results of snow and ice thickness distributions from the SIMBA 2007 experiment in the Bellingshausen Sea, using four different methods (ASPeCt ship observations, downward-looking camera imaging, ship-based electromagnetic induction (EM) sounding, and in situ measurements using ice drills). A snow freeboard and ice thickness model generated from in situ measurements was then applied to contemporaneous ICESat (satellite laser altimetry) measured freeboard to derive ice thickness at the ICESat footprint scale. Errors from in situ measurements and from ICESat freeboard estimations were incorporated into the model, so a thorough evaluation of the model and uncertainty of the ice thickness estimation from ICESat are possible. Our results indicate that ICESat derived snow freeboard and ice thickness distributions (asymmetrical unimodal tailing to right) for first-year ice (0.29 ± 0.14 m for mean snow freeboard and 1.06 ± 0.40 m for mean ice thickness), multi-year ice (0.48 ± 0.26 and 1.59 ± 0.75 m, respectively), and all ice together (0.42 ± 0.24 and 1.38 ± 0.70 m, respectively) for the study area seem reasonable compared with those values from the in situ measurements, ASPeCt observations, and EM measurements. The EM measurements can act as an appropriate supplement for ASPeCt observations taken hourly from the ship's bridge and provide reasonable ice and snow distributions under homogeneous ice conditions. Our proposed approaches: (1) of using empirical equations relating snow freeboard to ice thickness based on in situ measurements and (2) of using isostatic equations that replace snow depth with snow freeboard (or empirical equations that convert freeboard to snow depth), are efficient and important ways to derive ice thickness from ICESat altimetry at the footprint scale for Antarctic sea ice. Spatial and temporal snow and ice thickness from satellite altimetry for the BA sector and for the entire Southern Ocean is therefore possible.