Mean Dynamic Topography Research Articles

The static gravity field model DGM-1S from GRACE and GOCE data: computation, validation and an analysis of GOCE mission’s added value

We present a global static model of the Earth’s gravity field entitled DGM-1S based on GRACE and GOCE data. The collection of used data sets includes nearly 7 years of GRACE KBR data and 10 months of GOCE gravity gradient data. The KBR data are transformed with a 3-point differentiation into quantities that are approximately inter-satellite accelerations. Gravity gradients are processed in the instrumental frame. Noise is handled with a frequency-dependent data weighting. DGM-1S is complete to spherical harmonic degree 250 with a Kaula regularization being applied above degree 179. Its performance is compared with a number of other satellite-only GRACE/GOCE models by confronting them with (i) an independent model of the oceanic mean dynamic topography, and (ii) independent KBR and gravity gradient data. The tests reveal a competitive quality for DGM-1S. Importantly, we study added value of GOCE data by comparing the performance of satellite-only GRACE/GOCE models with models produced without GOCE data: either ITG-Grace2010s or EGM2008 depending on which of the two performs better in a given region. The test executed based on independent gravity gradients quantifies this added value as 25–38 % in the continental areas poorly covered with terrestrial gravimetry data (Equatorial Africa, Himalayas, and South America), 7–17 % in those with a good coverage with these data (Australia, North America, and North Eurasia), and 14 % in the oceans. This added value is shown to be almost entirely related to coefficients below degree 200. It is shown that this gain must be entirely attributed to gravity gradients acquired by the mission. The test executed based on an independent model of the mean dynamic topography suggests that problems still seem to exist in satellite-only GRACE/GOCE models over the Pacific ocean, where noticeable deviations between these models and EGM2008 are detected, too.

GNSS-R Derived Centimetric Sea Topography: An Airborne Experiment Demonstration

The results of two airborne experiments performed to test the precision and the relative accuracy of the conventional Global Navigation Satellite Systems Reflectometry (GNSS-R) technique employing only the C/A code are presented. The first and the second experiments demonstrate, respectively, a 17 cm precision for a 500 m flight altitude with a 8 km along-track spatial resolution, and a 6 cm precision for a 3000 m flight altitude with a 6.6 km along-track spatial resolution. In both, the Relative Mean Dynamic Topography (RMDT) is compared with results derived from traditional radar altimetry provided by Jason-2. The Root Mean Square (RMS) of the RMDT difference between both measurement systems is 48 cm for the first flight, and 198 cm for the second flight. During the second flight, the feasibility of the proposed technique to measure the sea slopes is demonstrated by superposing over the aircraft ground track the measured sea surface height with the geoid undulations, which are about 1 meter.

The apparent British sea slope is caused by systematic errors in the levelling-based vertical datum

The spirit-levelling–based British vertical datum (Ordnance Datum Newlyn) implies a south–north apparent slope in mean sea level of up to 53 mm deg –1 latitude, due to the datum falling on heading northwards. Although this apparent slope has been investigated since the 1960s, explanations of its origin have remained inconclusive. It has also been suggested that, rather than a slope, the British vertical datum includes a step of about 240 mm affecting all sites north of about 53°N. In either case, the British vertical datum may be of limited use for any study requiring accurate heights or changes in heights, such as testing geoid models, groundwater and hydrocarbon extraction, the calibration and validation of satellite-based digital terrain models, and the unification of vertical datums internationally. Within the last decade, however, based on an apparent reduction in the slope to only −12 mm deg –1 latitude with respect to recent geoid models, it has been claimed that the British vertical datum does provide a physically meaningful surface for use in scientific applications. In this paper, we reinvestigate the presence of apparent south–north sea slopes around Britain and reported slopes in the vertical datum, using the EGM2008 global gravitational model, together with mean sea level and GPS data from British tide gauges, GPS ellipsoidal heights of 178 fundamental benchmarks across mainland Britain, and vertical deflection observations at 192 stations. We demonstrate a south–north slope in the British vertical datum of −(20–25) mm deg –1 latitude with respect to both mean sea level (corrected for the ocean's mean dynamic topography and the inverse barometer response to atmospheric pressure loading) and the EGM2008 quasigeoid model, while EGM2008 is shown to exhibit a negligible slope of (2 ± 4) mm deg –1 with respect to mean sea level. It is clear, therefore, that the slope can only arise from systematic errors in the levelling, although we are unable to isolate their exact origin. Using an offset detection method based on a penalized likelihood maximization using the Schwarz Information Criterion, we do not detect a step in the vertical datum affecting all sites north of 53°N, but do identify regional distortions that we attribute to the inhomogeneity in both the levelling data used and the least squares adjustment procedures used to realize the datum. We conclude that the British vertical datum remains unsuitable for scientific purposes.

Multiple jets in the Malvinas Current

The velocity structure of the Malvinas Current is described based on the analysis of high‐resolution hydrographic data and direct current observations. The data show that though the current width exceeds 150 km, the flow is concentrated in two relatively narrow (~10–20 km) jets. Within these cores, the direct observations indicate surface velocities exceeding 0.5 m.s−1. Surface drifter, satellite‐derived mean dynamic topography, and sea surface temperature data suggest that the high‐velocity jets are also ubiquitous features of the time mean circulation. Both jets appear to be continuous features extending more than 900 km along the western slope of the Argentine Basin. These jets closely follow the 200 and 1400 m isobaths. Additional high‐velocity cores are apparent in direct current measurements and hydrographic observations, but these features are weaker and not continuous along the slope. Though the Malvinas Current transport is mostly barotropic, baroclinic jets are also identified in relative geostrophic velocity sections. The baroclinic jets are colocated with the barotropic jets. Our results suggest that the main Malvinas Current core is located over a relatively flat portion of the bottom, referred to as the Perito Moreno terrace. This observation is in agreement with recent seismic and geological evidence suggesting that in geological time scales the Malvinas Current played a key role in the configuration of the bottom sediments over the western slope of the Argentine Basin.

Analysis of a relative offset between vertical datums at the North and South Islands of New Zealand

The leveling networks realized by 13 different local vertical datums were jointly readjusted at the South and North Islands of New Zealand. The relation between these two leveling networks and the Word Height System was then defined using GPS-leveling data and the EGM08 global geopotential model. In this study, we investigate the relative offset between these two vertical datum realizations. This is done based on comparison of the geometric geoid/quasigeoid heights (obtained from GPS and newly adjusted leveling data) with the regional gravimetric geoid/quasigeoid solutions. Moreover, oceanographic and geodetic models of mean dynamic topography (MDT) are used to assess the relative offset between these two vertical datum realizations through the analysis of regional spatial variations of mean sea level (MSL). The comparison of GPS-leveling data with regional gravimetric solutions re- veals large systematic distortions (exceeding several deci- meters across New Zealand) between the geometric and gravimetric geoid/quasigeoid heights attributed mainly to systematic errors within regional gravimetric solutions. The presence of a significant offset between the vertical datum realizations at the North and South Islands is not confirmed. The MSL difference between tide gauges in Wellington and Dunedin of ∼24 cm is estimated based on the analysis of MDT models.

Gravity of the Arctic Ocean from satellite data with validations using airborne gravimetry: Oceanographic implications

Precise mappings of sea surface topography, slope, and gravity of the Arctic Ocean are derived from altimeter data collected by Envisat and ICESat. Both altimeters measured instantaneous sea surface height at leads in the sea ice. To reduce contamination by ice‐freeboard signal and tracker noise in Envisat height data, a retracking of the waveform data was performed. Analogous reprocessing of ICESat data was also done. Arctic mean sea surfaces (MSSs) were computed from Envisat data spanning 2002–2008 and ICESat data spanning 2003–2009. Farrell et al. (2012) used these “ICEn” MSSs to estimate mean dynamic topography (MDT). These same Envisat and ICESat data are used, in sea‐surface‐slope form, to compute the ARCtic Satellite‐only (ARCS‐2) altimetric marine gravity field. ARCS‐2 extends north to 86°N and uses GRACE/GOCE gravity data (GOCO02S) for its long‐wavelength (>260 km) components. Use of Envisat data improves the spatial resolution over that of existing Arctic marine gravity fields in many areas. ARCS‐2's spatial resolution aids in tracing tectonic fabric—e.g., extinct plate boundaries—over broad areas of the Arctic basin whose tectonic origin remains a mystery. ARCS‐2's precision is validated using NASA 2010/2011 Operation IceBridge (OIB) airborne gravimetry. ARCS‐2 and OIB gravity along with ICEn‐MSS results are employed to locate short‐wavelength errors approaching 1 m in current Arctic marine geoids (EGM2008). Precise OIB airborne gravity corroborates that such errors in current geoid/gravity models are widespread in Arctic areas lacking accurate surface gravity data. These geoid errors limit the spatial resolution at which MDT can be mapped.

Evaluation of global monitoring and forecasting systems at Mercator Océan

Abstract. Since December 2010, the MyOcean global analysis and forecasting system has consisted of the Mercator Océan NEMO global 1/4° configuration with a 1/12° nested model over the Atlantic and the Mediterranean. The open boundary data for the nested configuration come from the global 1/4° configuration at 20° S and 80° N. The data are assimilated by means of a reduced-order Kalman filter with a 3-D multivariate modal decomposition of the forecast error. It includes an adaptive-error estimate and a localization algorithm. A 3-D-Var scheme provides a correction for the slowly evolving large-scale biases in temperature and salinity. Altimeter data, satellite sea surface temperature and in situ temperature and salinity vertical profiles are jointly assimilated to estimate the initial conditions for numerical ocean forecasting. In addition to the quality control performed by data producers, the system carries out a proper quality control on temperature and salinity vertical profiles in order to minimise the risk of erroneous observed profiles being assimilated in the model. This paper describes the recent systems used by Mercator Océan and the validation procedure applied to current MyOcean systems as well as systems under development. The paper shows how refinements or adjustments to the system during the validation procedure affect its quality. Additionally, we show that quality checks (in situ, drifters) and data sources (satellite sea surface temperature) have as great an impact as the system design (model physics and assimilation parameters). The results of the scientific assessment are illustrated with diagnostics over the year 2010 mainly, assorted with time series over the 2007–2011 period. The validation procedure demonstrates the accuracy of MyOcean global products, whose quality is stable over time. All monitoring systems are close to altimetric observations with a forecast RMS difference of 7 cm. The update of the mean dynamic topography corrects local biases in the Indonesian Throughflow and in the western tropical Pacific. This improves also the subsurface currents at the Equator. The global systems give an accurate description of water masses almost everywhere. Between 0 and 500 m, departures from in situ observations rarely exceed 1 °C and 0.2 psu. The assimilation of an improved sea surface temperature product aims to better represent the sea ice concentration and the sea ice edge. The systems under development are still suffering from a drift which can only be detected by means of a 5-yr hindcast, preventing us from upgrading them in real time. This emphasizes the need to pursue research while building future systems for MyOcean2 forecasting.

Comparison of Geostrophic Current Estimation around the Korean Peninsula from Remote Sensing-based MSS Models

ABSTRACT Cho, J. M., Lee, D. H., Lee, M. R., Cho, J. H., Yun, H. S., Jeong, T. J., Jang, A. and Kim, T. Y., 2013. Comparison of Geostrophic Current Estimation around the Korean Peninsula from Remote Sensing-based MSS models Altimetry satellites like Geosat, ERS-1/2, Envisat and Topex-Poseidon, which have been developed and deployed since the 1970s, have collected immense arrays of data on the sea, and the data thus gathered have been used in oceanography as well as in a variety of scientific disciplines, such as geophysics and geodesy. The mean dynamic ocean topography or Sea Surface Topography (SST), obtained by calculating the difference between the Mean Sea Surface (MSS) and the geoid surface, plays a critical role in determining the earth's shape and geoid as well as in the precision modeling of the earth's gravity field, ocean currents, and tidal waves. In the study, the MSS_CNES_CLS10 (Centre National d'Etudes Spatiales Collectre Localisation Satellites 2010) and DTU10MSS (Danmarks Tekniske Universi...

Assessment of optimally filtered recent geodetic mean dynamic topographies

AbstractRecent geoids from the Gravity Recovery and Climate Experiment (GRACE) and the Gravity field and steady state Ocean Circulation Explorer satellite mission (GOCE) contain useful short‐scale information for the construction of a geodetic ocean mean dynamic topography (MDT). The geodetic MDT is obtained from subtracting the geoid from a mean sea surface (MSS) as measured by satellite altimetry. A gainful use of the MDT and an adequate assessment needs an optimal filtering. This is accomplished here by defining a cutoff length scale dmax for the geoid and applying a Gaussian filter with half‐width radius r on the MDT. A series of MDTs (GRACE, GOCE, and combined satellite‐only (GOCO) solutions) is tested, using different sets of filter parameters dmax and r. Optimal global and regional dependent filter parameters are estimated. To find optimal parameters and to assess the resulting MDTs, the geostrophic surface currents induced by the filtered geodetic MDT are compared to corrected near‐surface currents obtained from the Global Drifter Program (GDP). The global optimal cutoff degree and order (d/o) dmax (half‐width radius r of the spatial Gaussian filter) is 160 (1.1°) for GRACE; 180 (1.1–1.2°) for 1st releases of GOCE (time‐ and space‐wise methods) and GOCO models; and 210 (1.0 degree) for 2nd and 3rd releases of GOCE and GOCO models. The cutoff d/o is generally larger (smaller) and the filter length smaller (larger) for regions with strong, small‐scale (slow, broad scale) currents. The smallest deviations from the drifter data are obtained with the GOCO03s geoid model, although deviations of other models are only slightly higher.

The contribution of eddies to striations in absolute dynamic topography

AbstractDistinct 4 year averages of absolute dynamic topography reveal striations in all ocean basins during 1993–2008. Striations are alternating mesoscale jet‐like structures observed in time‐averaged zonal geostrophic velocity, . They are characterized by speeds O(1 cm s–1) and are nominally separated by 200 km in the meridional direction. Similar patterns have been observed in sea level anomaly, mean dynamic topography, and Argo float measurements. Use of a tracked‐eddy database in concert with a contour identification and eddy removal algorithm demonstrates that eddies are a dominant source of striations in in the South Pacific (20°S–50°S, 200°E–280°E). Eddies with lifetimes ≥ 4 weeks account for 46–57% of the variance in and correlation coefficients between total and eddy‐only are 0.90–0.93. Attention is given to the ability of the algorithm to correctly identify eddies and suggests that a more appropriate bound on the variance due to eddies is ∼ 30–70%. This permits the existence of latent zonal jets and/or β‐plumes. Additional findings of the study include (1) a large number of eddies having a broad range of amplitudes and scales contribute most to the eddy‐induced patterns and (2) the standard deviation of does not decay inversely with averaging period as proposed by a model of random eddies.

Accuracy of the Preliminary GOCE Geoid Models from an Oceanographic Perspective

In the framework of the ESA High Processing Facility (HPF), eight gravity models have been computed from the GOCE data since the beginning of the mission in 2009. The accuracy of these models for oceanographic application is assessed and compared with other geoid models (i.e., ITG-GRACE 2010S, EGM2008, EIGEN6C). For that purpose, ocean Mean Dynamic Topography (MDT) estimates are computed by subtracting the different geoids from an altimetric Mean Sea Surface. The MDT assessments were carried out by analyzing their associated geostrophic surface currents at different maximum harmonic degrees and by comparing them with independent in situ oceanographic data. Compared with GRACE-only geoid models, GOCE-only geoid models are in better agreement with independent data at scales smaller than 150 km while they show similar performances at scales larger than 200 km. Regional comparisons show that at 125 km the HPF releases compare well to independent observations with significant improvement with the last releases that use more GOCE data. At 100 km the western boundary currents are well resolved by the last releases, but the satellite-only geoid models need to be further improved to resolve smaller scales and reach the targeted 100 km resolution with higher accuracy.

Reconstruction of mean dynamic topography of the Black Sea for altimetry measurements

In this study a new approach for reconstructing the Mean Dynamic Topography of the Black Sea is applied. Constant observations (SVP measurements), drifters, and data of vertical sounding of the temperature and salt content together with measurements of sea level anomalies received from Topex/Poseidon mission satellite data were used. The absolute sea level received by altimetry data using the mean dynamic topography received during work was compared to the dynamic level received according to independent marine surveys. The comparison showed that the method represented in the study permits one to define more exactly the dynamic topography of the Black Sea when compared with the studies of previous authors. The results of this study will be useful to reconstruct the areas of the geostrophic currents according to satellite altimetry.

Global Surface Geostrophic Ocean Currents Derived from Satellite Altimetry and GOCE Geoid

The surface geostrophic currents (SGC) can be derived via the principle of geostrophy from the dynamic height of the ocean, or the deviation of the true, variable sea surface height with respect to the Earth's static geoid, both of which can be measured by geodetic means. Here we calculate the Mean Dynamic Topography (MDT) by subtracting the geoid height determined by the GOCE satellite mission from the Mean Sea Surface Topography (MSST) derived from multi-satellite ocean altimetry (T/P, Jason 1/2, ERS-1/2, GEOSAT). Results for SGC are compared with those obtained from a GRACE-based mean geoid, as well as with the mean circulation patterns from measurements done by in situ drifter buoys and from simulations of the ECCO Ocean General Circulation Model. We found GOCE-based geoid solution clearly leads to significant improvements in the spatial resolution of SGC globally except in the Equatorial band where special filtering may be needed, with current velocities and spatial patterns closest to the in situ measurements of currents, compared with the GRACE-based results or ECCO model simulations that give significantly weaker Values with lower spatial resolution.

Towards worldwide height system unification using ocean information

AbstractWe describe the application of ocean levelling to worldwide height system unification. The study involves a comparison of ‘geodetic’ and ‘ocean’ approaches to determination of the mean dynamic topography (MDT) at the coast, from which confidence in the accuracy of stateof- the-art ocean and geoid models can be obtained. We conclude that models are consistent at the sub-decimetre level for the regions that we have studied (North Atlantic coastlines and islands, North American Pacific coast and Mediterranean). That level of consistency provides an estimate of the accuracy of using the ocean models to provide an MDT correction to the national datums of countries with coastlines, and thereby of achieving unification. It also provides a validation of geoid model accuracy for application to height system unification in general. We show how our methods can be applied worldwide, as long as the necessary data sets are available, and explain why such an extension of the present study is necessary if worldwide height system unification is to be realised.

Nonlinear diffusion filtering of data on the Earth’s surface

The paper deals with data filtering on closed surfaces using linear and nonlinear diffusion equations. We define a surface finite-volume method to approximate numerically parabolic partial differential equations on closed surfaces, namely on a sphere, ellipsoid or the Earth’s surface. The closed surface as a computational domain is approximated by a polyhedral surface created by planar triangles and we construct a dual co-volume grid. On the co-volumes we define a weak formulation of the problem by applying Green’s theorem to the Laplace–Beltrami operator. Then the finite-volume method is applied to discretize the weak formulation. Weak forms of elliptic operators are expressed through surface gradients. In our numerical scheme we use a piece-wise linear approximation of a solution in space and the backward Euler time discretization. Furthermore, we extend a linear diffusion on surface to the regularized surface Perona–Malik model. It represents a nonlinear diffusion equation, which at the same time reduces noise and preserves main edges and other details important for a correct interpretation of the real data. We present four numerical experiments. The first one has an illustrative character showing how an additive noise is filtered out from an artificial function defined on a sphere. Other three examples deal with the real geodetic data on the Earth’s surface, namely (i) we reduce a stripping noise from the GOCE satellite only geopotential model up to degree 240, (ii) we filter noise from the real GOCE measurements (the component $$T_{zz})$$ , and (iii) we reduce a stripping noise from the satellite only mean dynamic topography at oceans. In all experiments we focus on a comparison of the results obtained by both the linear and nonlinear models presenting advantages of the nonlinear diffusion.

The impact of mean dynamic topography on a sea-level anomaly assimilation in the South China Sea based on an eddy-resolving model

The sea-level anomaly (SLA) from a satellite altimeter has a high accuracy and can be used to improve ocean state estimation by assimilation techniques. However, the lack of an accurate mean dynamic topography (MDT) is still a bothersome issue in an ocean data assimilation. The previous studies showed that the errors in MDT have significant impacts on assimilation results, especially on the time-mean components of ocean states and on the time variant parts of states via nonlinear ocean dynamics. The temporal-spatial differences of three MDTs and their impacts on the SLA analysis are focused on in the South China Sea (SCS). The theoretical analysis shows that even for linear models, the errors in MDT have impacts on the SLA analysis using a sequential data assimilation scheme. Assimilation experiments, based on EnOI scheme and HYCOM, with three MDTs from July 2003 to June 2004 also show that the SLA assimilation is very sensitive to the choice of different MDTs in the SCS with obvious differences between the experimental results and observations in the centre of the SCS and in the vicinity of the Philippine Islands. A new MDT for assimilation of SLA data in the SCS was proposed. The results from the assimilation experiment with this new MDT show a marked reduction (increase) in the RMSEs (correlation coefficient) between the experimental and observed SLA. Furthermore, the subsurface temperature field is also improved with this new MDT in the SCS.

The north‐south tilt in the Australian Height Datum is explained by the ocean's mean dynamic topography

Using geodetic and oceanographic data, we show that the apparent north‐south slope between the Australian Height Datum (AHD) and the geoid is caused almost completely by the ocean's time‐mean dynamic topography (MDT). This is because the AHD was constrained to zero height at local mean sea level at multiple tide gauges around the Australian continent. Using MDT models and corrected leveling data, almost all of the apparent north‐south slope can be removed from the AHD. An auxiliary observation is that a satellite‐only MDT model based on only around one year of GOCE data generates results commensurate with geodetic, oceanographic and combined MDT models.

The contribution of local gravimetric geoid models to the calibration of satellite altimetry data and an outlook of the latest GOCE GGM performance in Gavdos

The use of geoid heights has been one of the available methodologies utilized for the independent calibration/validation of altimeters on-board satellites. This methodology has been employed for long in the Gavdos dedicated cal/val facility (Crete, Greece), where calibration results for the Jason satellites have been estimated, both for ascending and descending passes. The present work gives a detailed overview of the methodology followed in order to estimate a high-resolution and accuracy gravimetric geoid model for the wider Gavdos area, in support of the on-going calibration work. To estimate the geoid model, the well-known remove-compute-restore method is used while residual geoid heights are estimated through least-squares collocation so that associated errors are determined as well. It is found that the estimated formal geoid errors from LSC along passes 018 and 109 of Jason satellites, used for the bias estimation, range between ±0.8–1.6cm. The so-derived geoid heights are employed in the determination of the Jason-2 altimeter bias for all available cycles (cycles 1-114, spanning the period from July 2008 to August 2011) together with the RioMed DOT model. From the results acquired the Jason-2 bias has been estimated to be +196.1±3.2mm for pass 109 and +161.9±5.1mm for pass 018. Within the same frame, the GOCE/GRACE-based geopotential model GOCO02s has been used to estimate the mean dynamic ocean topography and the steady-state circulation in the area around Gavdos. The so-derived DOT model was used to estimate the Jason-2 bias in an effort to evaluate the performance of satellite-only geoid models and investigate whether their spatial resolution and accuracy provides some improvement w.r.t. traditional local gravimetric geoids. From the results acquired with geoid heights from GOCO02s, the estimated Jason-2 bias deviates significantly from that of the local gravimetric model, which can be attributed to a possible mean offset and the low resolution of GOCE-based GGMs. On the other hand, when the newly estimated GOCE-based DOT was employed with geoid heights from the local gravimetric geoid model, the Jason-2 bias has been estimated to be +185.1±3.2mm for pass 109 and +130.2±5.1mm for pass 018.

High resolution dynamic ocean topography in the Southern Ocean from GOCE

A mean dynamic ocean topography (MDT) has been computed using a high resolution GOCE (Gravity field and steady-state Ocean Circulation Explorer) gravity model and a new mean sea surface obtained from a combination of satellite altimetry covering the period 1992 October till 2010 April. The considered gravity model is GO-CONS-GCF-2-TIM-R3, which computes geoid using 12 months of GOCE gravity field data. The GOCE gravity data allow for more detailed and accurate estimates of MDT. This is illustrated in the Southern Ocean where the commission error is reduced from 20 to 5cm compared to the MDT computed using the GRACE gravity field model ITG-Grace2010. As a result of the more detailed and accurate MDT, the calculation of geostrophic velocities from the MDT is now possible with higher accuracy and spatial resolution, and the error estimate is about 7 cms−1 for the Southern Ocean.

Performance of GOCE and GRACE-derived mean dynamic topographies in resolving Antarctic Circumpolar Current fronts

Presently, two satellite missions, Gravity Recovery and Climate Experiment (GRACE) and Gravity field and steady-state Ocean Circulation Explorer (GOCE), are making detailed measurements of the Earth’s gravity field, from which the geoid can be obtained. The mean dynamic topography (MDT) is the difference between the time-averaged sea surface height and the geoid. The GOCE mission is aimed at determining the geoid with superior accuracy and spatial resolution, so that a more accurate MDT can be estimated. In this study, we determine the mean positions of the Antarctic Circumpolar Current fronts using the purely geodetic estimates of the MDT constructed from an altimetric mean sea surface and GOCE and GRACE geoids. Overall, the frontal positions obtained from the GOCE and GRACE MDTs are close to each other. This means that these independent estimates are robust and can potentially be used to validate frontal positions obtained from sparse and irregular in situ measurements. The geodetic frontal positions are compared to earlier estimates as well as to those derived from MDTs based on satellite and in situ measurements and those obtained from an ocean data synthesis product. The position of the Sub-Antarctic Front identified in the GOCE MDT is found to be in better agreement with the previous estimates than that identified in the GRACE MDT. The geostrophic velocities derived from the GOCE MDT are also closer to observations than those derived from the GRACE MDT. Our results thus show that the GOCE mission represents an improvement upon GRACE in terms of the time-averaged geoid.