Search

This paper presents the RDSAR, SAR/SARin & FF-SAR altimetry processors available in the ESA Altimetry Virtual Lab (AVL) hosted on the EarthConsole® platform. An overview on processors and features as well as preliminary analyses using AVL output data are reported to demonstrate the quality of the ESA Altimetry Virtual Lab altimetry services. Soon additional processors from ESA research contracts will be added to the AVL portfolio to continue providing innovative solutions to the radar altimetry community.

Information on sea level and its temporal and spatial variability is of great importance for various scientific, societal, and economic issues. This article reports about a new sea level dataset for the North Sea (named North SEAL) of monthly sea level anomalies (SLAs), absolute sea level trends, and amplitudes of the mean annual sea level cycle over the period 1995–2019. Uncertainties and quality flags are provided together with the data. The dataset has been created from multi-mission cross-calibrated altimetry data preprocessed with coastal dedicated approaches and gridded with an innovative least-squares procedure including an advanced outlier detection to a 6–8 km wide triangular mesh. The comparison of SLAs and tide gauge time series shows good consistency, with average correlations of 0.85 and maximum correlations of 0.93. The improvement with respect to existing global gridded altimetry solutions amounts to 8 %–10 %, and it is most pronounced in complicated coastal environments such as river mouths or regions sheltered by islands. The differences in trends at tide gauge locations depend on the vertical land motion model used to correct relative sea level trends. The best consistency with a median difference of 0.04±1.15 mm yr−1 is reached by applying a recent glacial isostatic adjustment (GIA) model. With the presented sea level dataset, for the first time, a regionally optimized product for the entire North Sea is made available. It will enable further investigations of ocean processes, sea level projections, and studies on coastal adaptation measures. The North SEAL data are available at https://doi.org/10.17882/79673 (Müller et al., 2021).

EOT20 is the latest in a series of empirical ocean tide (EOT) models derived using residual tidal analysis of multi-mission satellite altimetry at DGFI-TUM. The amplitudes and phases of 17 tidal constituents are provided on a global 0.125∘ grid based on empirical analysis of seven satellite altimetry missions and four extended missions. The EOT20 model shows significant improvements compared to the previous iteration of the global model (EOT11a) throughout the ocean, particularly in the coastal and shelf regions, due to the inclusion of more recent satellite altimetry data as well as more missions, the use of the updated FES2014 tidal model as a reference to estimated residual signals, the inclusion of the ALES retracker and improved coastal representation. In the validation of EOT20 using tide gauges and ocean bottom pressure data, these improvements in the model compared to EOT11a are highlighted with the root sum square (RSS) of the eight major tidal constituents improving by ∼ 1.4 cm for the entire global ocean with the major improvement in RSS (∼ 2.2 cm) occurring in the coastal region. Concerning the other global ocean tidal models, EOT20 shows an improvement of ∼ 0.2 cm in RSS compared to the closest model (FES2014) in the global ocean. Variance reduction analysis was conducted comparing the results of EOT20 with FES2014 and EOT11a using the Jason-2, Jason-3 and SARAL satellite altimetry missions. From this analysis, EOT20 showed a variance reduction for all three satellite altimetry missions with the biggest improvement in variance occurring in the coastal region. These significant improvements, particularly in the coastal region, provide encouragement for the use of the EOT20 model as a tidal correction for satellite altimetry in sea-level research. All ocean and load tide data from the model can be freely accessed at https://doi.org/10.17882/79489 (Hart-Davis et al., 2021). The tide gauges from the TICON dataset used in the validation of the tide model, are available at https://doi.org/10.1594/PANGAEA.896587 (Piccioni et al., 2018a).

Vertical land motion (VLM) at the coast is a substantial contributor to relative sea level change. In this work, we present a refined method for its determination, which is based on the combination of absolute satellite altimetry (SAT) sea level measurements and relative sea level changes recorded by tide gauges (TGs). These measurements complement VLM estimates from the GNSS (Global Navigation Satellite System) by increasing their spatial coverage. Trend estimates from the SAT and TG combination are particularly sensitive to the quality and resolution of applied altimetry data as well as to the coupling procedure of altimetry and TGs. Hence, a multi-mission, dedicated coastal along-track altimetry dataset is coupled with high-frequency TG measurements at 58 stations. To improve the coupling procedure, a so-called “zone of influence” (ZOI) is defined, which confines coherent zones of sea level variability on the basis of relative levels of comparability between TG and altimetry observations. Selecting 20 % of the most representative absolute sea level observations in a 300 km radius around the TGs results in the best VLM estimates in terms of accuracy and uncertainty. At this threshold, VLMSAT-TG estimates have median formal uncertainties of 0.58 mm yr−1. Validation against GNSS VLM estimates yields a root mean square (rmsΔVLM) of VLMSAT-TG and VLMGNSS differences of 1.28 mm yr−1, demonstrating the level of accuracy of our approach. Compared to a reference 250 km radius selection, the 300 km zone of influence improves trend accuracies by 15 % and uncertainties by 35 %. With increasing record lengths, the spatial scales of the coherency in coastal sea level trends increase. Therefore, the relevance of the ZOI for improving VLMSAT-TG accuracy decreases. Further individual zone of influence adaptations offer the prospect of bringing the accuracy of the estimates below 1 mm yr−1.

In the context of the ESA Climate Change Initiative project, we are engaged in a regional reprocessing of high-resolution (20 Hz) altimetry data of the classical missions in a number of the world's coastal zones. It is done using the ALES (Adaptive Leading Edge Subwaveform) retracker combined with the X-TRACK system dedicated to improve geophysical corrections at the coast. Using the Jason-1 and Jason-2 satellite data, high-resolution, along-track sea level time series have been generated, and coastal sea level trends have been computed over a 14-year time span (from July 2002 to June 2016). In this paper, we focus on a particular coastal site where the Jason track crosses land, Senetosa, located south of Corsica in the Mediterranean Sea, for two reasons: (1) the rate of sea level rise estimated in this project increases significantly in the last 4–5 km to the coast compared to what is observed further offshore, and (2) Senetosa is the calibration site for the TOPEX/Poseidon and Jason altimetry missions, which are equipped for that purpose with in situ instrumentation, in particular tide gauges and a Global Navigation Satellite System (GNSS) antenna. A careful examination of all the potential errors that could explain the increased rate of sea level rise close to the coast (e.g., spurious trends in the geophysical corrections, imperfect inter-mission bias estimate, decrease of valid data close to the coast and errors in waveform retracking) has been carried out, but none of these effects appear able to explain the trend increase. We further explored the possibility that it results from real physical processes. Change in wave conditions was investigated, but wave setup was excluded as a potential contributor because the magnitude was too low and too localized in the immediate vicinity of the shoreline. A preliminary model-based investigation about the contribution of coastal currents indicates that it could be a plausible explanation of the observed change in sea level trend close to the coast.

Sea state data are of major importance for climate studies, marine engineering, safety at sea and coastal management. However, long-term sea state datasets are sparse and not always consistent, and sea state data users still mostly rely on numerical wave models for research and engineering applications. Facing the urgent need for a sea state climate data record, the Global Climate Observing System has listed “Sea State” as an Essential Climate Variable (ECV), fostering the launch in 2018 of the Sea State Climate Change Initiative (CCI). The CCI is a programme of the European Space Agency, whose objective is to realise the full potential of global Earth observation archives established by ESA and its member states in order to contribute to the ECV database. This paper presents the implementation of the first release of the Sea State CCI dataset, the implementation and benefits of a high-level denoising method, its validation against in situ measurements and numerical model outputs, and the future developments considered within the Sea State CCI project. The Sea State CCI dataset v1 is freely available on the ESA CCI website (http://cci.esa.int/data, last access: 25 August 2020) at ftp://anon-ftp.ceda.ac.uk/neodc/esacci/sea_state/data/v1.1_release/ (last access: 25 August 2020). Three products are available: a multi-mission along-track L2P product (http://dx.doi.org/10.5285/f91cd3ee7b6243d5b7d41b9beaf397e1, Piollé et al., 2020a), a daily merged multi mission along-track L3 product (http://dx.doi.org/10.5285/3ef6a5a66e9947d39b356251909dc12b, Piollé et al., 2020b) and a multi-mission monthly gridded L4 product (http://dx.doi.org/10.5285/47140d618dcc40309e1edbca7e773478, Piollé et al., 2020c).

The ESA Earth Explorer CryoSat-2 was launched on 8 April 2010 to monitor the precise changes in the thickness of terrestrial ice sheets and marine floating ice. To do that, CryoSat orbits the planet at an altitude of around 720 km with a retrograde orbit inclination of 92∘ and a quasi repeat cycle of 369 d (30 d subcycle). To reach the mission goals, the CryoSat products have to meet the highest quality standards to date, achieved through continual improvements of the operational processing chains. The new CryoSat Ice Baseline-D, in operation since 27 May 2019, represents a major processor upgrade with respect to the previous Ice Baseline-C. Over land ice the new Baseline-D provides better results with respect to the previous baseline when comparing the data to a reference elevation model over the Austfonna ice cap region, improving the ascending and descending crossover statistics from 1.9 to 0.1 m. The improved processing of the star tracker measurements implemented in Baseline-D has led to a reduction in the standard deviation of the point-to-point comparison with the previous star tracker processing method implemented in Baseline-C from 3.8 to 3.7 m. Over sea ice, Baseline-D improves the quality of the retrieved heights inside and at the boundaries of the synthetic aperture radar interferometric (SARIn or SIN) acquisition mask, removing the negative freeboard pattern which is beneficial not only for freeboard retrieval but also for any application that exploits the phase information from SARIn Level 1B (L1B) products. In addition, scatter comparisons with the Beaufort Gyre Exploration Project (BGEP; https://www.whoi.edu/beaufortgyre, last access: October 2019) and Operation IceBridge (OIB; Kurtz et al., 2013) in situ measurements confirm the improvements in the Baseline-D freeboard product quality. Relative to OIB, the Baseline-D freeboard mean bias is reduced by about 8 cm, which roughly corresponds to a 60 % decrease with respect to Baseline-C. The BGEP data indicate a similar tendency with a mean draft bias lowered from 0.85 to −0.14 m. For the two in situ datasets, the root mean square deviation (RMSD) is also well reduced from 14 to 11 cm for OIB and by a factor of 2 for the BGEP. Observations over inland waters show a slight increase in the percentage of good observations in Baseline-D, generally around 5 %–10 % for most lakes. This paper provides an overview of the new Level 1 and Level 2 (L2) CryoSat Ice Baseline-D evolutions and related data quality assessment, based on results obtained from analyzing the 6-month Baseline-D test dataset released to CryoSat expert users prior to the final transfer to operations.

Coupled with the occurrence of regional/local sea level rise on urbanized coastal cities is the possibility of land subsidence that contaminates the measurement by the tide gauge (TG) sensors. Another technology that could possibly check the in-situ data from tide gauge is satellite altimetry. The sea surface height (SSH) measured from satellite altimeter is compared with the observed tide gauge sea level (TGSL) to detect vertical land motion (VLM). This study used satellite altimeter retracked products near the TG Stations in Tagbilaran, Bohol; Dumaguete, Negros Oriental; and Mambajao, Camiguin located in the vicinity of the South East Bohol Trench and Fault (SEBTF).Based on the results, the TG site in Tagbilaran is undergoing land subsidence. The rate of VLM is around 5 mm/year from 2009 to 2017. The same trend was manifested in the GNSS observed data in the PHIVOLCS monitoring station in Tagbilaran and the geodetic levelling done in the area. After the October 15, 2013 earthquake in Bohol, downward trends of around 27 mm/year and 17 mm/year were observed from GNSS measurements and SSH-TGSL difference respectively. These different rates may be due to the distance between the two sensors. The comparison between SSH and TGSL in Dumaguete showed small difference with a VLM rate of 1.8 mm/year. The difference in SSH-TGSL in Mambajao is quite large with a downward rate of 9.4 mm/year. This result needs to be further investigated for TG or TGBM instability or monitored for a possibility of land uplift.

A deeper knowledge about geostrophic ocean surface currents in the northern Nordic Seas supports the understanding of ocean dynamics in an area affected by sea ice and rapidly changing environmental conditions. Monitoring these areas by satellite altimetry results in a fragmented and irregularly distributed data sampling and prevents the computation of homogeneous and highly resolved spatio-temporal datasets. In order to overcome this problem, an ocean model is used to fill in data when altimetry observations are missing. The present study provides a novel dataset based on a combination of along-track satellite-altimetry-derived dynamic ocean topography (DOT) elevations and simulated differential water heights (DWHs) from the Finite Element Sea ice Ocean Model (FESOM) version 1.4. This innovative dataset differs from classical assimilation methods because it substitutes altimetry data with the model output when altimetry fails or is not available. The combination approach is mainly based on a principal component analysis (PCA) after reducing both quantities by their constant and seasonal signals. In the main step, the most-dominant spatial patterns of the modeled differential water heights as provided by the PCA are linked with the temporal variability in the estimated DOT from altimetry by performing a principal component synthesis (PCS). After the combination, the annual signal obtained by altimetry and a constant offset are re-added in order to reference the final data product to the altimetry height level. Surface currents are computed by applying the geostrophic flow equations to the combined topography. The resulting final product is characterized by the spatial resolution of the ocean model around 1 km and the temporal variability in the altimetry along-track derived DOT heights. The combined DOT is compared to an independent DOT product, resulting in a positive correlation of about 80 %, to provide more detailed information about short periodic and finer spatial structures. The derived geostrophic velocity components are evaluated by in situ surface drifter observations. Summarizing all drifter observations in equally sized bins and comparing the velocity components shows good agreement in spatial patterns, magnitude and flow direction. Mean differences of 0.004 m s−1 in the zonal and 0.02 m s−1 in the meridional component are observed. A direct pointwise comparison between the combined geostrophic velocity components interpolated onto the drifter locations indicates that about 94 % of all residuals are smaller than 0.15 m s−1. The dataset is able to provide surface circulation information within the sea ice area and can be used to support a deeper comprehension of ocean currents in the northern Nordic Seas affected by rapid environmental changes in the 1995–2012 time period. The data are available at https://doi.org/10.1594/PANGAEA.900691 (Müller et al., 2019).

The dynamic ocean topography (DOT) of the polar seas can be described by satellite altimetry sea surface height observations combined with geoid information as well as by ocean models. The altimetry observations are characterized by an irregular sampling and seasonal sea ice coverage complicating reliable DOT estimations. Models display various spatiotemporal resolutions but are limited to their computational and mathematical context and introduced forcing models. In the present paper, ALES+ retracked altimetry ranges and derived along-track DOT heights of ESA's Envisat and water heights of the Finite Element Sea Ice-Ocean Model (FESOM) are compared to investigate similarities and discrepancies. The goal of the present paper is to identify to what extent pattern and variability of the northern Nordic seas derived from measurements and model agree with each other, respectively. The study period covers the years 2003–2009. An assessment analysis regarding seasonal DOT variabilities shows good agreement and confirms the dominant impact of the annual signal in both datasets. A comparison based on estimated regional annual signal components shows 2–3 times stronger amplitudes of the observations but good agreement of the phase. Reducing both datasets by constant offsets and the annual signal reveals small regional residuals and highly correlated DOT time series (Pearson linear correlation coefficient of at least 0.67). The highest correlations can be found in areas that are ice-free and affected by ocean currents. However, differences are visible in sea-ice-covered shelf regions. Furthermore, remaining constant artificial elevations in the observational data can be attributed to an insufficient representation of the used geoid. In general, the comparison results in good agreement between simulated and altimetry-based descriptions of the DOT in the northern Nordic seas.

Sea level is a very sensitive index of climate change since it integrates the impacts of ocean warming and ice mass loss from glaciers and the ice sheets. Sea level has been listed as an essential climate variable (ECV) by the Global Climate Observing System (GCOS). During the past 25 years, the sea level ECV has been measured from space by different altimetry missions that have provided global and regional observations of sea level variations. As part of the Climate Change Initiative (CCI) program of the European Space Agency (ESA) (established in 2010), the Sea Level project (SL_cci) aimed to provide an accurate and homogeneous long-term satellite-based sea level record. At the end of the first phase of the project (2010–2013), an initial version (v1.1) of the sea level ECV was made available to users (Ablain et al., 2015). During the second phase of the project (2014–2017), improved altimeter standards were selected to produce new sea level products (called SL_cci v2.0) based on nine altimeter missions for the period 1993–2015 (https://doi.org/10.5270/esa-sea_level_cci-1993_2015-v_2.0-201612; Legeais and the ESA SL_cci team, 2016c). Corresponding orbit solutions, geophysical corrections and altimeter standards used in this v2.0 dataset are described in detail in Quartly et al. (2017). The present paper focuses on the description of the SL_cci v2.0 ECV and associated uncertainty and discusses how it has been validated. Various approaches have been used for the quality assessment such as internal validation, comparisons with sea level records from other groups and with in situ measurements, sea level budget closure analyses and comparisons with model outputs. Compared with the previous version of the sea level ECV, we show that use of improved geophysical corrections, careful bias reduction between missions and inclusion of new altimeter missions lead to improved sea level products with reduced uncertainties on different spatial and temporal scales. However, there is still room for improvement since the uncertainties remain larger than the GCOS requirements (GCOS, 2011). Perspectives on subsequent evolution are also discussed.