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

Table 5 Comparison of surface elevations by Seasat and optical survey along EGIG line

Surface elevation maps of the southern half of the Greenland subcontinent are produced from radar altimeter data acquired by the Seasat satellite. A summary of the processing procedure and examples of return waveform data are given. The elevation data are used to generate a regular grid which is then computer contoured to provide an elevation contour map. Ancillary maps show the statistical quality of the elevation data and various characteristics of the surface. The elevation map is used to define ice flow directions and delineate the major drainage basins. Regular maps of the Jakobshavns Glacier drainage basin and the ice divide in the vicinity of Crete Station are presented. Altimeter derived elevations are compared with elevations measured both by satellite geoceivers and optical surveying.

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.

(Table 1) Ice thickness at lake Nam Co, Tibet, during 2008-2009

Lake ice change is one of the sensitive indicators of regional and global climate change. Different sources of data are used in monitoring lake ice phenology nowadays. Visible and Near Infrared bands of imagery (VNIR) are well suited for the observation of freshwater ice change, for example data from AVHRR and MODIS. Active and passive microwave data are also used for the observation of lake ice, e.g., from satellite altimetry and radiometry, backscattering coefficient from QuickSCAT, brightness temperature (Tb) from SSM/I, SMMR, and AMSR-E. Most of the studies are about lake ice cover phenology, while few studies focus on lake ice thickness. For example, Hall et al. using 5 GHz (6 cm) radiometer data showed a good relationship between Tb and ice thickness. Kang et al. found the seasonal evolution of Tb at 10.65 GHz and 18.7 GHz from AMSR-E to be strongly influenced by ice thickness. Many studies on lake ice phenology have been carried out since the 1970s in cold regions, especially in Canada, the USA, Europe, the Arctic, and Antarctica. However, on the Tibetan Plateau, very little research has focused on lake ice-cover change; only a small number of published papers on Qinghai Lake ice observations. The main goal of this study is to investigate the change in lake ice phenology at Nam Co on the Tibetan Plateau using MODIS and AMSR-E data (monitoring the date of freeze onset, the formation of stable ice cover, first appearance of water, and the complete disappearance of ice) during the period 2000-2009.

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.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.

Improved geoid solution for the Weddell Sea region

We present a geoid solution for the Weddell Sea and adjacent continental Antarctic regions. There, a refined geoid is of interest, especially for oceanographic and glaciological applications. For example, to investigate the Weddell Gyre as a part of the Antarctic Circumpolar Current and, thus, of the global ocean circulation, the mean dynamic topography (MDT) is needed. These days, the marine gravity field can be inferred with high and homogeneous resolution from altimetric height profiles of the mean sea surface. However, in areas permanently covered by sea ice as well as in coastal regions, satellite altimetry features deficiencies. Focussing on the Weddell Sea, these aspects are investigated in detail. In these areas, ground-based data that have not been used for geoid computation so far provide additional information in comparison with the existing high-resolution global gravity field models such as EGM2008. The geoid computation is based on the remove-compute-restore approach making use of least-squares collocation. The residual geoid with respect to a release 4 GOCE model adds up to two meters and more in the near-coastal and continental areas of the Weddell Sea region, also in comparison with EGM2008. Consequently, the thus refined geoid serves to compute new estimates of the regional MDT and geostrophic currents.

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.

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