Altimeter Missions Research Articles

Potential of cold-atom airborne gravimetry to improve coastal gravity field and quasigeoid modelling

We investigate using the GIRAFE cold-atom gravimeter during an airborne gravity survey for improving gravity field and quasigeoid modelling. The study is conducted over the Bay of Biscay, France. Geoid/quasigeoid determination is usually a major challenge over such coastal areas due to scarce and inconsistent gravity data. In a first step, the GIRAFE dataset is analysed and compared with available surface gravity data as well as with global altimetry models from UCSD and DTU. The comparisons indicate that the DTU model is better than the UCSD model within around 10 km from the coastline. Furthermore, recent satellite altimeter missions significantly improve the altimetry models in coastal areas. A significant bias (− 4.00 mGal) in shipborne data is also found from this comparison. In a second step, eight quasigeoid solutions are calculated to evaluate the contribution of GIRAFE data. This contribution reaches 3 cm in terms of height anomaly for DTU21 while being much larger for UCSDv31 and shipborne data. Finally, the quasigeoid solutions are validated using GNSS-levelling data. The results indicate that using GIRAFE data improves by approximately 50% the quality of quasigeoid models over land near the coast. The highest accuracy, around 1 cm, is achieved when GIRAFE data are merged with refined gravity data. Importantly, the standard deviation is just 1.2 cm when compared with GNSS-levelling points if using only GIRAFE data over marine areas, which is very close to the 1 cm goal of geoid/quasigeoid model determination in modern geodesy. This study thus confirms the benefits of performing airborne gravity survey using quantum sensors.

CryoSat Long-Term Ocean Data Analysis and Validation: Final Words on GOP Baseline-C

ESA’s Earth explorer mission CryoSat-2 has an ice-monitoring objective, but it has proven to also be a valuable source of observations for measuring impacts of climate change over oceans. In this paper, we report on our long-term ocean data analysis and validation and give our final words on CryoSat-2’s Geophysical Ocean Products (GOP) Baseline-C. The validation is based on a cross comparison with concurrent altimetry and with in situ tide gauges. The highlights of our findings include GOP Baseline-C showing issues with the ionosphere and pole tide correction. The latter gives rise to an east–west pattern in range bias. Between Synthetic Aperture Radar (SAR) and Low-Resolution Mode (LRM), a 1.4 cm jump in range bias is explained by a 0.5 cm jump in sea state bias, which relates to a significant wave height SAR-LRM jump of 10.5 cm. The remaining 0.9 cm is due to a range bias between ascending and descending passes, exhibiting a clear north–south pattern and ascribed to a timing bias of +0.367 ms, affecting both time-tag and elevation. The overall range bias of GOP Baseline-C is established at −2.9 cm, referenced to all calibrated concurrent altimeter missions. The bias drift does not exceed 0.2 mm/yr, leading to the conclusion that GOP Baseline-C is substantially stable and measures up to the altimeter reference missions. This is confirmed by tide gauge comparison with a selected set of 309 PSMSL tide gauges over 2010–2022: we determined a correlation of R = 0.82, a mean standard deviation of σ=5.7 cm (common reference and GIA corrected), and a drift of 0.17 mm/yr. In conclusion, the quality, continuity, and reference of GOP Baseline-C is exceptionally good and stable over time, and no proof of any deterioration or platform aging has been found. Any improvements for the next CryoSat-2 Baselines could come from sea state bias optimization, ionosphere and pole tide correction improvement, and applying a calibrated value for any timing biases.

Steric heights of submesoscale processes and internal gravity waves in the subtropical northwestern Pacific and northern South China Sea as revealed by moored observations

The Surface Water and Ocean Topography (SWOT) altimeter mission provides a good opportunity to globally detect submesoscale processes (submesoscales) from high-resolution sea surface height (SSH) measurements. However, in addition to the balanced mesoscale eddies and submesoscales, the measured SSH also contains signals of unbalanced internal gravity waves (IGWs) including internal tides (ITs) and supertidal IGWs (IGWs with frequencies of 3–8 cycles per day). Therefore, a quantitative understanding of the contributions of the above multiscale dynamic processes to the SSH is a prerequisite for exploring submesoscales from the SWOT-derived SSH data. Here, this issue was investigated through analyzing one-year moored observation-derived near-surface steric heights (SHs) in the subtropical northwestern Pacific (NWP) and northern South China Sea (SCS) which have different dynamic backgrounds. In the sense of annual mean, the root-mean-squared (RMS) SHs of submesoscales (defined as processes with periods of 2–15 days here) and IGWs in the SCS are 1.30 and 5.54 cm, respectively, while they are only 0.70 and 3.19 cm in the NWP. Seasonally, the RMS SH of submesoscales is larger in winter than summer in both regions; with respect to the IGWs, their RMS SH is larger in summer than winter in the SCS but it is opposite in the NWP. Generally, the SHs of diurnal and semidiurnal ITs are dominated by the stationary component in both regions but the proportion of the nonstationary component is much larger in the NWP than the SCS. The sum of the SHs by nonstationary diurnal and semidiurnal ITs exceeds that of submesoscales except in winter in the SCS. Furthermore, we also found that the SHs of supertidal IGWs are always larger than those of submesoscales in both regions and both seasons. By decomposing the IGWs with frequencies of 1–8 cpd (i.e., cycles per day) into different baroclinic modes, we found that for all of the IGWs bands, the SHs are overwhelmingly attributed to the first mode. The results in this study will provide useful information for the application of SWOT data in regions with active IGWs.

DPOD2020: A DORIS extension of the ITRF2020 for precise orbit determination

As one of the tracking systems used to determine orbits of the altimeter mission satellites (such as TOPEX/Poseidon, Envisat, Jason-1/-2/-3, CryoSat-2, Saral/Altika, Sentinel-3A/-3B, HY-2A/C/D, Jason-CS/Sentinel-6A, SWOT), DORIS (Doppler Orbitography Radiopositionning Integrated by Satellite) allows to determine positions and velocities of tracking stations that define a stable reference for the estimation of the precise orbits and thus are fundamental for the quality of the altimeter data and derived mean sea level products. Due to the time evolution of the DORIS ground network, some stations included in the 2020 realization of the International Terrestrial Reference Frame (ITRF2020) have been decommissioned and since 2021.0 a few new stations were added to the tracking network. Therefore, to satisfy operational requirements for POD (Precise Orbit Determination) and routine delivery of geodetic products, the International DORIS Service (IDS) regularly updates the DPOD (DORIS terrestrial reference frame for Precise Orbit Determination). The DPOD solutions include mean positions and velocities of all the DORIS stations since 1993.0 derived from the stacking of the latest IDS weekly combined series aligned to the current ITRF. In this paper, we first present the stacking process of the DPOD2020 version 1.0. Then, we address the validation procedure of the DPO2020 including comparison with ITRF2020 and POD tests. For eighty percent of all the time segments of all the DORIS stations, the station position differences between DPOD2020 version 1.0 and ITRF2020 are smaller than ten millimeters. The major position differences between these two solutions are associated with the DORIS sites either localized in the South Atlantic Anomaly region or with time spans smaller than one year. Compared to DPOD2014, the DPOD2020 shows reduction of the main statistics of the DORIS-to-DORIS tie residuals (differences between the estimated and measured ties). The POD tests showed similar results for DPOD2020 and DPOD2014 for most of the altimetric satellites (TOPEX/Poseidon, Jason-1, CryoSat-2, Jason-3). In addition, we observed better POD results with DPOD2020 for the latest altimetric satellite Sentinel-3A as well as a slight degradation for Jason-2. That degradation was fully explained by slightly worst results for the stations localized in the South Atlantic Anomaly region.

Ocean 2D eddy energy fluxes from small mesoscale processes with SWOT

Abstract. We investigate ocean dynamics at different scales in the Agulhas Current system, a region of important interocean exchange of heat and energy. While ocean observations and some of the most advanced climate models capture the larger mesoscale dynamics (> 100 km), the smaller-scale fronts and eddies are underrepresented. The recently launched NASA–CNES Surface Water and Ocean Topography (SWOT) wide-swath altimeter mission observes the smaller ocean geostrophic scales down to 15 km in wavelength globally. Here we will analyse different eddy diagnostics in the Agulhas Current region and quantify the contributions from the larger mesoscales observable today and the smaller scales to be observed with SWOT. Surface geostrophic diagnostics of eddy kinetic energy, strain, and energy cascades are estimated from modelled sea surface height (SSH) fields of the Massachusetts Institute of Technology general circulation model (MITgcm) latitude–longitude polar cap (LLC4320) simulation subsampled at 1/10∘. In this region, the smaller scales (<150 km) have a strong signature on the horizontal geostrophic strain rate and for all eddy diagnostics in the Western Boundary Current and along the meandering Agulhas Extension. We investigate the horizontal cascade of energy using a coarse-graining technique, and we observe that the wavelength range where the inverse cascade occurs is biased towards larger mesoscale wavelengths with today’s altimetric sampling. We also calculate the projected sampling of the eddy diagnostics under the SWOT swaths built with the NASA–CNES simulator to include the satellite position and realistic noise. For the swaths, a neural network noise mitigation method is implemented to reduce the residual SWOT random error before calculating eddy diagnostics. In terms of SSH, observable wavelengths of 15 to 20 km are retrieved after neural network noise mitigation, as opposed to wavelengths larger than 40 km before the noise reduction.

In Situ Validation of Altimetry and CFOSAT SWIM Measurements in a High Wave Environment

Abstract While satellite altimeters have revolutionized ocean science, validation measurements in high wave environments are rare. Using geodetic Global Navigation Satellite System (GNSS) data collected from the Southern Ocean Flux Station (SOFS; −47°S, 142°E) since 2019, as part of the Southern Ocean Time Series (SOTS), we present a validation of satellite missions in this energetic region. Here we show that high rate GNSS observations at SOFS can successfully measure waves in the extreme conditions of the Southern Ocean and obtain robust measurements in all wave regimes [significant wave height (SWH) ranging from 1.5 to 12.6 m]. We find good agreement between the in situ and nadir altimetry SWH (RMSE = 0.16 m, mean bias = 0.04 m, and n = 60). Directional comparisons with the Chinese–French Ocean Satellite (CFOSAT) Surface Waves Investigation and Monitoring (SWIM) instrument also show good agreement, with dominant directions having an RMSE of 9.1° (n = 22), and correlation coefficients between the directional spectra ranging between 0.57 and 0.79. Initial sea level anomaly (SLA) estimates capture eddies propagating through the region. Comparisons show good agreement with daily gridded SLA products (RMSE = 0.03 m, and n = 205), with scope for future improvement. These results demonstrate the utility of high rate geodetic GNSS observations on moored surface platforms in highly energetic regions of the ocean. Such observations are important to maximize the geophysical interpretation from altimeter missions. In particular, the ability to provide collocated directional wave observations and SLA estimates will be useful for the validation of the recently launched Surface Water and Ocean Topography (SWOT) mission where understanding the interactions between sea state and sea surface height poses a major challenge.

Recovering crustal thickness of the Moon from an extended Vening Meinesz-Moritz inversion

Crustal thickness (Moho surface relief) is closely related to geotectonic, which is the predominant research topic in lunar geophysics. For lunar crustal depth inversion, to address the problem of downward continuation in traditional spectral domain methodologies, we adopt a solution from gravity disturbance (geophysical gravity anomaly). It is available to investigate crustal thickness at arbitrary surface outside the geoid. Using gravity data from the Gravity Recovery and Interior Laboratory (GRAIL) mission and terrain information recovered from the Lunar Orbiter Laser Altimeter (LOLA) mission, global crustal depth of the Moon is recovered by the extended Vening Meinesz-Moritz method. The first application result of this technique on the Moon is validated by four types of released models in the form of spherical harmonic expansion with 310°. Moreover, from the statistics, our lunar Moho solution demonstrates that the average crustal thickness is approximately 41.5 km. The maximum value of 83.1 km and the minimum value of 0.8 km are located in the area of high mountains and Moscoviense mascon, respectively.

On global wave height climatology and trends from multiplatform altimeter measurements and wave hindcast

Reliable estimates of ocean surface waves and their long-term global trends are fundamentally important for numerous ocean engineering and geophysical applications. Nonetheless, it is scientifically and technically demanding to accurately identify such relatively small trends [O(0.1−1) cm yr−1] from wave data currently available. In 2022, Young and Ribal demonstrated that by changing the calibration method, multiplatform altimeter observations yield noticeable differences in wave height estimates and their trends. In this paper, we use our global wave hindcast as a baseline and confirm that observations calibrated relative to one another (i.e., “altimeter-altimeter” calibration) ultimately outperform traditionally “altimeter-buoy” calibrated data. The former leads to more consistent mean and extreme wave heights across simultaneously operating altimeters and consequently to a more homogeneous evolution of global sea states over the past several decades. The underlying reason is hypothesized to be that temporal inhomogeneities in long-term NDBC buoy records due to changes in buoy hardware and software over time introduce comparable inhomogeneities to long-term altimeter observations when the “altimeter-buoy” calibration is adopted. It is encouraging to observe that the global trends of mean and 90th-percentile wave heights over the period 1992-2019 from our hindcast and re-calibrated altimeter data are very close: approximately 0.3 cm yr−1. Nonetheless, the altimeter-based trend of 99th-percentile wave height, Hs99, is substantially-but-spuriously stronger mainly due to markedly growing availability of altimeter missions over time. Our work clearly demonstrates the power of third generation spectral wave models as diagnostics tools. Given the uncertainties embedded in buoy and satellite observations, wave modeling and global wave hindcasts can indeed guide these measurements. Further analysis of the underlying climate physics driving the interannual variability of ocean waves in wind-sea dominated regions is also presented.

Satellite Altimetry for Ocean and Coastal Applications: A Review

More than 30 years of observations from an international suite of satellite altimeter missions continue to provide key data enabling research discoveries and a broad spectrum of operational and user-driven applications. These missions were designed to advance technologies and to answer scientific questions about ocean circulation, ocean heat content, and the impact of climate change on these Earth systems. They are also a valuable resource for the operational needs of oceanographic and weather forecasting agencies that provide information to shipping and fishing vessels and offshore operations for route optimization and safety, as well as for other decision makers in coastal, water resources, and disaster management fields. This time series of precise measurements of ocean surface topography (OST)—the “hills and valleys” of the ocean surface—reveals changes in ocean dynamic topography, tracks sea level variations at global to regional scales, and provides key information about ocean trends reflecting climate change in our warming world. Advancing technologies in new satellite systems allows measurements at higher spatial resolution ever closer to coastlines, where the impacts of storms, waves, and sea level rise on coastal communities and infrastructure are manifest. We review some collaborative efforts of international space agencies, including NASA, CNES, NOAA, ESA, and EUMETSAT, which have contributed to a collection of use cases of satellite altimetry in operational and decision-support contexts. The extended time series of ocean surface topography measurements obtained from these satellite altimeter missions, along with advances in satellite technology that have allowed for higher resolution measurements nearer to coasts, has enabled a range of such applications. The resulting body of knowledge and data enables better assessments of storms, waves, and sea level rise impacts on coastal communities and infrastructure amongst other key contributions for societal benefit. Although not exhaustive, this review provides a broad overview with specific examples of the important role of satellite altimetry in ocean and coastal applications, thus justifying the significant resource contributions made by international space agencies in the development of these missions.

The CNES CLS 2022 Mean Sea Surface: Short Wavelength Improvements from CryoSat-2 and SARAL/AltiKa High-Sampled Altimeter Data

A new mean sea surface (MSS) was determined by focusing on the accuracy provided by exact-repeat altimetric missions (ERM) and the high spatial coverage of geodetic (or drifting) missions. The goal was to obtain a high-resolution MSS that would provide centimeter-level precision. Particular attention was paid to the homogeneity of the oceanic content of this MSS, and specific processing was also carried out, particularly on the data from the geodetic missions. For instance, CryoSat-2 and SARAL/AltiKa data sampled at high frequencies were enhanced using a dedicated filtering process and corrected from oceanic variability using the results of the objective analysis of sea-level anomalies provided by DUACS multi-missions gridded sea-level anomalies fields (MSLA). Particular attention was also paid to the Arctic area by combining traditional sea-surface height (SSH) with the sea levels estimated within fractures in the ice (leads). The MSS was determined using a local least-squares collocation technique, which provided an estimation of the calibrated error. Furthermore, our technique takes into account altimetric noises, ocean-variability-correlated noises, and along-track biases, which are determined independently for each observation. Moreover, variable cross-covariance models were fitted locally for a more precise determination of the shortest wavelengths, which were shorter than 30 km. The validations performed on this new MSS showed an improvement in the finest topographic structures, with amplitudes exceeding several cm, while also continuing to refine the correction of the oceanic variability. Overall, the analysis of the precision of this new CNES_CLS 2022 MSS revealed an improvement of 40% compared to the previous model, from 2015.

The Rise and Fall of Alaska and Yukon Glaciers Detected by TOPEX/Poseidon and Jason‐2 Altimeters Using a Novel Glacier‐Threshold Method

AbstractRadar altimetry has been used to monitor sea level changes and ice sheet elevation changes for decades. Over mountain glaciers, radar altimetry has limited applications due to contaminated waveforms caused by complex glacier surfaces and steep terrains. In this study, we develop a glacier‐threshold method (GTM) to determine glacier elevation changes in Alaska and Yukon. The GTM retracks waveforms, reduces terrain effect and detects invalid elevation observations from the TOPEX/Poseidon (T/P) and Jason‐2 (J2) altimeters, resulting in an average usable rate of 35% from original altimeter‐measured heights. The selected measurements are used to construct time series of glacier elevation changes over 1993–2002 (T/P) and 2008–2016 (J2) at 47 sites. A crossover analysis, validation by airborne laser altimetry observations, and comparisons with recent studies confirm the estimated glacier elevation changes. Our findings suggest that the thinning of glaciers in Alaska and Yukon has been ubiquitous in recent years. The site near Walsh Glacier has the highest thinning rate of −5.71 ± 0.09 m/yr, followed by Chitina Glacier at −4.51 ± 0.21 m/yr. Only 17% of sites show glacier thickening due to surges and mass accumulations. Using a sophisticated data processing algorithm like the GTM, we show that altimeter data from the repeat TOPEX‐Jason missions can be used to monitor long‐term glacier elevation changes at inaccessible spots, creating an additional value for altimeter missions originally purposed for monitoring long‐term sea level change.

Reconstructing Long-Term Arctic Sea Ice Freeboard, Thickness, and Volume Changes from Envisat, CryoSat-2, and ICESat-2

Satellite altimeters have been used to monitor Arctic sea ice (ASI) thickness for several decades, but whether the different altimeter missions (such as radar and laser altimeters) are in agreement with each other and suitable for long-term research needs to be investigated. To analyze the spatiotemporal characteristics of ASI, continuous long-term first-year ice, and multi-year ice of ASI freeboard, thickness, and volume from 2002 to 2021 using the gridded nadirization method from Envisat, CryoSat-2, and ICESat-2, altimeter data are comprehensively constructed and assessed. The influences of sea surface temperature (SST) and sea surface wind field (SSW) on ASI are also discussed. The freeboard/thickness and extent/area of ASI all varied seasonally and reached their maximum and minimum in April and October, March and September, respectively. From 2002 to 2021, the freeboard, thickness, extent, and area of ASI all consistently showed downward trends, and sea ice volume decreased by 5437 km3/month. SST in the Arctic rose by 0.003 degrees C/month, and the sea ice changes lagged behind this temperature variation by one month between 2002 and 2021. The meridional winds blowing from the central Arctic region along the eastern coast of Greenland to the North Atlantic each month are consistent with changes in the freeboard and thickness of ASI. SST and SSW are two of the most critical factors driving sea ice changes. This study provides new data and technical support for monitoring ASI and exploring its response mechanisms to climate change.

4DVarNet-SSH: end-to-end learning of variational interpolation schemes for nadir and wide-swath satellite altimetry

Abstract. The reconstruction of sea surface currents from satellite altimeter data is a key challenge in spatial oceanography, especially with the upcoming wide-swath SWOT (Surface Water and Ocean and Topography) altimeter mission. Operational systems, however, generally fail to retrieve mesoscale dynamics for horizontal scales below 100 km and timescales below 10 d. Here, we address this challenge through the 4DVarnet framework, an end-to-end neural scheme backed on a variational data assimilation formulation. We introduce a parameterization of the 4DVarNet scheme dedicated to the space–time interpolation of satellite altimeter data. Within an observing system simulation experiment (NATL60), we demonstrate the relevance of the proposed approach, both for nadir and nadir plus SWOT altimeter configurations for two contrasting case study regions in terms of upper ocean dynamics. We report a relative improvement with respect to the operational optimal interpolation between 30 % and 60 % in terms of the reconstruction error. Interestingly, for the nadir plus SWOT altimeter configuration, we reach resolved space–timescales below 70 km and 7 d. The code is open source to enable reproducibility and future collaborative developments. Beyond its applicability to large-scale domains, we also address the uncertainty quantification issues and generalization properties of the proposed learning setting. We discuss further future research avenues and extensions to other ocean data assimilation and space oceanography challenges.

An improved triple collocation-based integration of multiple gravity anomaly grids from satellite altimetry: Contribution of ICESat-2

Satellite altimetry has been widely used to determine the marine gravity fields. Different altimeter missions lead to different gravity grids with varying spatial resolutions and accuracies. Marine gravity accuracy in offshore waters is often low. Integrating marine gravity anomalies (MGA) derived from various altimeter missions are key to obtaining high-precision and high-resolution global marine gravity fields. Here an improved triple collocation (ITC) method is first proposed to establish a new 1´×1′ marine gravity anomaly in the South China Sea (SCS), named MGAITC, by merging the grids of MGA from the Haiyang-2A (HY-2A), CryoSat-2, and ICESat-2 altimeter measurements. A combined MGA (MGACOM) is also constructed by using the least squares collocation (LSC) method from three sets of altimeter-derived geoid gradients. The qualities of MGAITC, MGACOM and each altimeter-derived MGA grids are assessed using shipborne gravity and different global gravity models over waters of varying depths and distances from the coastline. The accuracy of the MGAITC gravity against the shipborne gravity is 5.52 mGal, which represents improvements of 7%, 9%, 13%, and 3% relative to the gravity grids derived only from ICESat-2, CryoSat-2, HY-2A and from MGACOM. In all assessments involving different sub-regions, seafloors, and distances from the coastline, MGAITC performs the best, and ICESat-2 yields the best single-mission MGA grid, especially in offshore waters. We find that ICESat-2 has a high potential to improve shallow-water gravity accuracy and the global marine gravity field can be improved by integrating gravity grids from different altimetry missions using ITC.

Global submesoscale diagnosis using along-track satellite altimetry

Abstract. The ocean's sea surface height (SSH) field is a complex mix of motions in geostrophic balance and unbalanced motions including high-frequency tides, internal tides, and internal gravity waves. Barotropic tides are well estimated for altimetric SSH in the open ocean, but the SSH signals of internal tides remain. The transition scale, Lt, at which these unbalanced ageostrophic motions dominate balanced geostrophic motions is estimated for the first time using satellite altimetry. Lt is critical to define the spatial scales above which surface geostrophic currents can be inferred from SSH gradients. We use a statistical approach based on the analysis of 1 Hz altimetric SSH wavenumber spectra to obtain four geophysical parameters that vary regionally and seasonally: the background error, the spectral slope in the mesoscale range, a second spectral slope at smaller scales, and Lt. The mesoscale slope and error levels are similar to previous studies based on satellite altimetry. The break in the wavenumber spectra to a flatter spectral slope can only be estimated in midlatitude regions where the signal exceeds the altimetric noise level. Small values of Lt are observed in regions of energetic mesoscale activity, while larger values are observed towards low latitudes and regions of lower mesoscale activity. These results are consistent with recent analyses of in situ observations and high-resolution models. Limitations of our results and implications for reprocessed nadir and future swath altimetric missions are discussed.

Absolute Calibration of the Chinese HY-2B Altimetric Mission with Fiducial Reference Measurement Standards

This research and collaboration work aims at the calibration and validation (Cal/Val) of the Chinese HY-2B satellite altimeter based upon two permanent Cal/Val facilities: (1) the China Altimetry Calibration Cooperation Plan in Qingdao, Bohai Sea and the Wanshan islands, China and (2) the permanent facility for altimetry calibration established by the European Space Agency in Crete, Greece. The HY-2B satellite altimeter and its radiometer have been calibrated and monitored using uniform, standardized procedures, as well as protocols and best practices, and they also built upon trusted and indisputable reference standards at both Cal/Val infrastructures in Europe and China. The HY-2B altimeter is thus monitored in a coordinated, absolute, homogeneous, long-term and worldwide manner. Calibration of altimeters is accomplished by examining satellite observations in open seas against reference measurements. Comparisons are established through precise satellite positioning, water level observations, GPS buoys and reference models (geoid, mean dynamic topography, earth tides, troposphere and ionosphere), all defined at the Cal/Val sites. In this work, the final uncertainty for the altimeter bias will be attributed to several individual sources of uncertainty, coming from observations in water level, atmosphere, absolute positioning, reference surface models, transfer of heights from Cal/Val sites to satellite observations, etc. Through this project, the procedures, protocols and best practices, originally developed in the course of the ESA FRM4ALT project, are updated, upgraded and followed at both Cal/Val facilities in Europe and China. All in all, the HY-2B satellite altimeter observes the sea level quite well and within its specifications.

Random, Environmental, and Representativeness Errors in Ocean Remote Sensing Versus In Situ Data: An Example of Wave Heights From Altimeters

Understanding the errors and uncertainties is important for the applications of oceanic remote sensing data products. In this study, an error modelling method is presented to decompose the error of ocean remote sensing data into three parts: random errors, environmental errors, and representativeness errors. Random errors refer to the zero-mean error uncorrelated to any external variable. Environmental errors are errors systematically dependent on some external or environmental variables. Representativeness errors are errors coming from the difference with the defined truth with respect to measured spatial-temporal ranges. These error components can be estimated with the collocations of only two systems if continuous observations are available from both systems. The methods of estimating and suppressing the magnitudes of these error components are also presented and are demonstrated using altimeter significant wave heights from the CCI-sea state dataset including eleven altimeter missions. For most altimeters, the three error components are within the same order of magnitude. After suppressing these error terms using the presented methods, the root-mean-square errors between altimeter and in-situ measurements of SWHs can be decreased for 25%~45% for these altimeters. The understanding of the contributions of different error components in remote sensing data products can guide the direction to improve data accuracy since different error components need to be suppressed in different ways, and it also provides the reference for applications such as data assimilation of models.

Baroclinic Tidal Cusps from Satellite Altimetry

Abstract The spatially averaged frequency spectrum of sea level has been computed at 4 cycle-per-year resolution and a Nyquist frequency of 0.5 cycles per hour using dual-satellite crossover data from the Jason and CryoSat-2 satellite altimeter missions. The novelty of the analysis is that it reveals unambiguous peaks due to high-frequency tidal signals, even after removing the predicted barotropic tide, without the usual aliasing caused by altimeter sampling. The tidal continuum, that is, a tidal cusp, is present in the spectrum in the diurnal and semidiurnal tidal bands, and a Lorentzian model spectrum has been fit within each band to identify the properties of the non-phase-locked tidal variability. An interesting feature of the semidiurnal tidal continuum is the unambiguous presence of an inner and an outer band, characterized by different Lorentzian bandwidths of roughly (180 day)−1 and (30 day)−1. Considering different latitude ranges, it is clear that the tidal continuum is most prominent in the range from −30° to 30° latitude. Within this range, it is found that 1.05-cm2 variance is associated with the semidiurnal continuum, and slightly less than half of this variance, 0.41 cm2, is associated with the slower, (180 day)−1 bandwidth, variability. The ratio of non-phase-locked to total baroclinic variability is about 62% in this latitude band, a value that is consistent with previous model-based estimates for this quantity. Quantification of the properties of the tidal continuum poleward of 30° latitude is not possible with the present data, due to the small size of the tidal signal compared to the mesoscale variability and other sources of noise.

Internal tides off the Amazon shelf during two contrasted seasons: interactions with background circulation and SSH imprints

Abstract. The Amazon shelf break is a key region for internal tide (IT) generation. It also shows a large seasonal variation in circulation and associated stratification. This study, based on a high-resolution model (1/36∘) explicitly forced by tide, aims to better characterize how the ITs vary between two contrasted seasons. During the season from March to July (MAMJJ) the currents and mesoscale eddies are weak while the pycnocline is shallower and stronger. From August to December (ASOND) mean currents and mesoscale eddies are strong, and the pycnocline is deeper and weaker than in MAMJJ. For both seasons, semi-diurnal M2 ITs are generated on the shelf break mainly between the 100 and 1000 m isobath in the model. South of 2∘ N, the conversion from barotropic to baroclinic tide is more efficient in MAMJJ than in ASOND. Local dissipation of the coherent M2 at the generation sites is higher in MAMJJ (30 %) than in ASOND (22 %), because higher modes are favourably generated (modes 2 and 3), making the internal wave packet more dispersive. The remaining fraction (70 %–80 %) propagates away from the generation sites and mainly dissipates locally every ∼ 100 km, which corresponds to the mode 1 reflection beams. About 13 %, 30 %, and 40 % of the M2 coherent IT dissipates at the first, second, and third beams. M2 coherent baroclinic flux propagates more northward during MAMJJ while it seems to be blocked at 6∘ N during ASOND. There is no intensified dissipation of the coherent M2 that could explain the disappearance of the coherent flux. In fact, the flux at this location becomes more incoherent because of strong interaction with the currents. This has been shown in the paper using 25 h mean snapshots of the baroclinic flux that shows branching and stronger eastward deviation of the IT when interacting with mesoscale eddies and stratification during ASOND. Finally, we evaluated sea surface height (SSH) frequency and wavenumber spectra for subtidal (f<1/28 h−1), tidal (1/28 h−1 < f<1/11 h−1), and supertidal (f>1/11 h−1) frequencies. Tidal frequencies explain most of the SSH variability for wavelengths between 250 and 70 km. Below 70 km, the SSH is mainly incoherent and supertidal. The length scale at which the SSH becomes dominated by unbalanced (non-geostrophic) IT was estimated to be around 250 km. Our results highlight the complexity of correctly predicting IT SSH in order to better observe mesoscale and sub-mesoscale eddies from existing and upcoming altimetric missions, notably the Surface Water Ocean Topography (SWOT) mission.

Cloud-based framework for inter-comparing submesoscale-permitting realistic ocean models

Abstract. With the increase in computational power, ocean models with kilometer-scale resolution have emerged over the last decade. These models have been used for quantifying the energetic exchanges between spatial scales, informing the design of eddy parametrizations, and preparing observing networks. The increase in resolution, however, has drastically increased the size of model outputs, making it difficult to transfer and analyze the data. It remains, nonetheless, of primary importance to assess more systematically the realism of these models. Here, we showcase a cloud-based analysis framework proposed by the Pangeo project that aims to tackle such distribution and analysis challenges. We analyze the output of eight submesoscale-permitting simulations, all on the cloud, for a crossover region of the upcoming Surface Water and Ocean Topography (SWOT) altimeter mission near the Gulf Stream separation. The cloud-based analysis framework (i) minimizes the cost of duplicating and storing ghost copies of data and (ii) allows for seamless sharing of analysis results amongst collaborators. We describe the framework and provide example analyses (e.g., sea-surface height variability, submesoscale vertical buoyancy fluxes, and comparison to predictions from the mixed-layer instability parametrization). Basin- to global-scale, submesoscale-permitting models are still at their early stage of development; their cost and carbon footprints are also rather large. It would, therefore, benefit the community to document the different model configurations for future best practices. We also argue that an emphasis on data analysis strategies would be crucial for improving the models themselves.