Technical note: Two types of absolute dynamic ocean topography

Abstract. Two types of marine geoid exist with the first type being the average level
of sea surface height (SSH) if the water is at rest (classical
definition), and the second type being satellite-determined with the
condition that the water is usually not at rest. The differences between
the two are exclusion (inclusion) of the gravity anomaly and non-measurable
(measurable) in the first (second) type. The associated absolute dynamic
ocean topography (referred to as DOT), i.e., SSH minus marine geoid,
correspondingly also has two types. Horizontal gradients of the first type of
DOT represent the absolute surface geostrophic currents due to water being at
rest on the first type of marine geoid. Horizontal gradients of the second
type of DOT represent the surface geostrophic currents relative to flow on
the second type of marine geoid. Difference between the two is quantitatively
identified in this technical note through comparison between the first type
of DOT and the mean second type of DOT (MDOT). The first type of DOT is
determined by a physical principle that the geostrophic balance takes the
minimum energy state. Based on that, a new elliptic equation is derived for
the first type of DOT. The continuation of geoid from land to ocean leads to
an inhomogeneous Dirichlet boundary condition with the boundary values taking
the satellite-observed second type of MDOT. This well-posed elliptic equation
is integrated numerically on 1∘ grids for the world oceans with the
forcing function computed from the World Ocean Atlas (T, S) fields and
the sea-floor topography obtained from the ETOPO5 model of NOAA. Between the
first type of DOT and the second type of MDOT, the relative root-mean square
(RRMS) difference (versus RMS of the first type of DOT) is 38.6 % and the
RMS difference in the horizontal gradients (versus RMS of the horizontal
gradient of the first type of DOT) is near 100 %. The standard deviation
of horizontal gradients is nearly twice larger for the second type
(satellite-determined marine geoid with gravity anomaly) than for the first
type (geostrophic balance without gravity anomaly). Such a difference needs
further attention from oceanographic and geodetic communities, especially the
oceanographic representation of the horizontal gradients of the second type
of MDOT (not the absolute surface geostrophic currents).


Abstract. Two types of marine geoid exist with the first type being the average level of sea surface height (SSH) if the water is at rest (classical definition), and the second type being satellite-determined with the condition that the water is usually not at rest. The differences between the two are exclusion (inclusion) of the gravity anomaly and nonmeasurable (measurable) in the first (second) type. The associated absolute dynamic ocean topography (referred to as DOT), i.e., SSH minus marine geoid, correspondingly also has two types. Horizontal gradients of the first type of DOT represent the absolute surface geostrophic currents due to water being at rest on the first type of marine geoid. Horizontal gradients of the second type of DOT represent the surface geostrophic currents relative to flow on the second type of marine geoid. Difference between the two is quantitatively identified in this technical note through comparison between the first type of DOT and the mean second type of DOT (MDOT). The first type of DOT is determined by a physical principle that the geostrophic balance takes the minimum energy state. Based on that, a new elliptic equation is derived for the first type of DOT. The continuation of geoid from land to ocean leads to an inhomogeneous Dirichlet boundary condition with the boundary values taking the satellite-observed second type of MDOT. This well-posed elliptic equation is integrated numerically on 1 • grids for the world oceans with the forcing function computed from the World Ocean Atlas (T , S) fields and the sea-floor topography obtained from the ETOPO5 model of NOAA. Between the first type of DOT and the second type of MDOT, the relative root-mean square (RRMS) difference (versus RMS of the first type of DOT) is 38.6 % and the RMS difference in the horizontal gradients (versus RMS of the horizontal gradient of the first type of DOT) is near 100 %. The standard deviation of horizontal gradients is nearly twice larger for the second type (satellite-determined marine geoid with gravity anomaly) than for the first type (geostrophic balance without gravity anomaly). Such a difference needs further attention from oceanographic and geodetic communities, especially the oceanographic representation of the horizontal gradients of the second type of MDOT (not the absolute surface geostrophic currents).

Introduction
Let the coordinates (x, y, z) be in zonal, latitudinal, and vertical directions. The absolute dynamic ocean topography (hereafter referred to as DOT)D is the sea surface height (SSH, waves and tides filtered out) relative to the marine geoid (i.e., the equipotential surface), where S is the SSH;N is the marine geoid height above the reference ellipsoid ( Fig. 1).D is an important signal in oceanography andN is of prime interest in geodesy. Equation (1) is also applicable if defined relative to the center of the Earth. The geoid heightN (x, y) and other associated measurable quantities such as gravity anomaly g(x, y) are related to the anomaly of the gravitational potential V (x, y, z) to a first approximation by the well-known Brun's formula (e.g., Hofmann-Wellenhof and Moritz, 2005), where g = 9.81 m s −2 , is the globally mean normal gravity, which is usually represented by g 0 in geodesy. The gravity anomaly is the vertical derivative of the potential where the anomaly of the gravity potential V satisfies the Laplace equation The vertical deflection is the slope of the geoid which connects to the gravity anomaly by Equation (6) links the vertical gravity gradient to the horizontal Laplacian of the marine geoid heightN and serves as the basic principle in the satellite marine geodesy. Sincê D is the difference of the two large fields S andN (2 orders of magnitude larger thanD), it is extremely sensitive to any error in either S orN -even a 1 % error in either field can lead to error inD that is of the same order of magnitude asD itself (Wunsch and Gaposchkin, 1980;Bingham et al., 2008).
Before satellites came into operation, S was measured from sparse surveying ships and tide gauge stations located along irregular local coastlines. However,N was not easy to observe. Without satellite measurements, the marine geoid is defined as the average level of SSH if the water is at rest and denoted here by N, which is called the classical marine geoid (or first type of marine geoid, Fig. 1a). The first type of marine geoid can be taken as a stand-alone concept in oceanography since it is on the basis of the hypothesis (mean SSH when the water at rest) without using the gravity anomaly. In this framework, the geostrophic balance and hydrostatic balance are used for large-scale (i.e., scale > 100 km) processes. Here (u g , v g ) are geostrophic current components; f is the Coriolis parameter; (p,ρ) are in situ pressure and density, respectively, which can be decomposed intô Here, ρ 0 = 1025 kg m −3 is the characteristic density; (ρ, p) are horizontally uniform with ρ vertically increasing with depth (stable stratification) where n(z) is the buoyancy frequency (or called the Brunt-Vaisala frequency); (p, ρ) are anomalies of pressure and density. Near the ocean surface, it is common to use the characteristic density and corresponding pressure (p 0 , ρ 0 ) to represent (p,ρ). Vertical integration of Eq. (8) from N to S after replacing (p,ρ) by (p 0 , ρ 0 ) in Eqs. (7) and (8) leads to is the first type of DOT. Since the first type of marine geoid (N) is defined as the average level of SSH if the water is at rest, the horizontal gradient of D represents the absolute surface geostrophic currents. After satellites came into operation, SSH has been observed with uniquely sampled temporal and spatial resolutions by high-precision altimetry above a reference ellipsoid (not geoid; Fu and Haines, 2013). Two Gravity Recovery and Climate Experiment (GRACE) satellites, launched in 2002, provide data to compute the marine geoid (called the GRACE Gravity Model, GGM; see website: http://www.csr.utexas. edu/grace/ (last access: 31 July 2018); Tapley et al., 2003;Shum et al., 2011). In addition, the European Space Agency's GOCE mission data, along with the GRACE data, have produced the best mean gravity field or the geoid model at a spatial scale longer than 67 km half-wavelength (or spherical harmonics completed to degree 300). This marine geoid is the solution of Eq. (6), where N * is the satellite-determined marine geoid from the measurable gravity anomaly g, and called the second type of marine geoid (Fig. 1b), which is different from N, defined by Eq. (13). Correspondingly, the second type of DOT is defined by where N * (t) changes with time due to temporally varying gravity anomaly g. Thus, comparison between the firsttype and second-type geoids should be conducted between N and N * . Here, N * is the temporally mean of N * (t). As for DOT, the first type of DOT (D) should be compared to the second type mean DOT (MDOT), The oceanic conditions at N and N * are different: water is at rest on N (see Eq. 13), but in motion on N * . The oceanographic community ignores such a difference, also treating horizontal gradients of the second type of DOT as the absolute surface geostrophic currents. For example, the second type of MDOT (D * ) data are posted at the NASA/JPL website: https://grace.jpl.nasa.gov/data/get-data/ dynamic-ocean-typography/ (last access: 31 July 2018); its horizontal gradients are also taken as the absolute surface geostrophic currents.
A question arises: do the horizontal gradients of the second type of MDOT (D * ) represent the absolute surface geostrophic currents? This paper will answer the question using the temporally averaged SSH and marine geoid satellite altimetric and gravimetric measurements from NASA (i.e., the second type of MDOT (D * )), and solving a new elliptic equation of D numerically. Given (S, N * , D) leads to the answer of the question.
The rest of the paper is outlined as follows. Section 2 describes the change in DOT due to the change in marine geoid from the first to the second type. Section 3 describes geostrophic currents and energy related to the first type of DOT. Section 4 presents the governing equation of the first type of DOT with the boundary condition at the coasts. Section 5 shows the numerical solution for the world oceans. Section 6 evaluates the change in global DOT from first to second type with oceanographic implications. Section 7 concludes the study.

Change in DOT from first to second type
The second type of MDOT (D * ) data are downloaded from the NASA/JPL website: https://grace.jpl.nasa.gov/data/ get-data/dynamic-ocean-typography/. This dataset is a subtraction of a second type of marine geoid of GRACE (Bingham et al., 2011) from a mean (1993 to 2006) altimetric sea surface. The change in marine geoid from first (N ) to second (N * ) type is represented by Correspondingly, the change in DOT is given by where Eqs. (12) and (15) are used. D is of interest in oceanography. N is of interest in geodesy. Equation (17) shows that the key issue to evaluate D is to determine D (i.e., first type of DOT). Conservation of potential vorticity for a dissipation-free fluid does not apply precisely to sea water where the density is a function not only of temperature and pressure but also of the dissolved salts. The effect of salinity on density is very important in the distribution of water properties. However, for most dynamic studies the effect of the extra state variable is not significant and the conservation of potential vorticity is valid (Veronis, 1980). Based on the conservation of the potential vorticity, the geostrophic current reaches the minimum energy state (Appendix A). Due to the minimum energy state, an elliptic partial differential equation for D is derived with coefficients containing sea-floor topography H , and forcing function containing temperature and salinity fields.
If D is negligible in comparison to D, the change in marine geoid from N to N * does not change the oceanographic interpretation of absolute DOT, i.e., the horizontal gradients of D * also represent the absolute surface geostrophic currents. If D is not negligible, the horizontal gradient of D * does not represent the absolute surface geostrophic currents.

Geostrophic currents and energy
Equation (9) implies, Using the first type of marine geoid N, the horizontal gradients of D lead to the absolute surface geostrophic currents (see Eqs. 11 and 13). Integration of the thermal wind relation, from the ocean surface to depth z leads to depth-dependent geostrophic currents, 4 P. C. Chu: Technical Note: Two types of absolute dynamic ocean topography are the baroclinic geostrophic currents. Here, f = 2 sin(ϕ) is the Coriolis parameter; = 2π/(86 400 s) is the mean Earth rotation rate; and ϕ is the latitude. The volume-integrated total energy, i.e., sum of kinetic energy of the geostrophic currents and the available potential energy (Oort et al., 1989), for an ocean basin (W ) is given by Substitution of Eqs. (21) and (22) into Eq. (23) leads to 4 Governing equation of D For a given density field, the second integration in the right side of Eq. (24) is known. The geostrophic currents taking the minimum energy state provides a constraint for D, The three-dimensional integration Eq. (25) over the ocean basin is conducted by where R is the horizontal area of the water volume, H is the water depth. Thus, Eq. (25) becomes where the parameters (X, Y ) are given by which represent vertically integrated baroclinic geostrophic currents scaled by the factor f/g (unit: m). Here, Eq. (18) is used (i.e., horizontal gradient of in situ density is the same as that of density anomaly). The Euler-Lagragian equation of the functional Eq. (27) is given by Substitution of Eq. (28) into Eq. (31) gives an elliptic partial differential equation (i.e., the governing equation) for the first type of DOT (i.e., D), where a = 6370 km, is the mean Earth radius. The geostrophic balance does not exist at the Equator. The Coriolis parameter f needs some special treatment for low latitudes. In this study, f is taken as 2 sin(5π/180) if latitude is between 10 • N and 0 • ; and as −2 sin(5π/180) if latitude is between 0 • and 10 • S. Let be the coastline of the ocean basin. Continuation of the geoid from land to oceans gives which leads to  Here, N l is the geoid over land. The boundary condition Eq. (36) can be rewritten as which is the boundary condition of D.

Difference between the two DOTs
The first-type global DOT (D i,j ; Fig. 2a) is the numerical solution of the elliptic Eq. (32) with the boundary condition Eq. (37). The second type global MDOT (D * i,j ; Fig. 2b) is downloaded from the NASA/JPL website: https://grace.jpl. nasa.gov/data/get-data/dynamic-ocean-typography/ (last access: 31 July 2018). The difference between the two DOTs, is evident in the world oceans (Fig. 2c). Here, (i, j ) denote the horizontal grid point. The relative root-mean square (RRMS) of D is given by where M = 38 877 is the number of total grid points. Both D and D * have positive and negative values. The arithmetic mean values (0.524 and −3.84 cm) are much smaller than the RMS mean values. They are 1 order of magnitude smaller than the corresponding standard deviations (54.9 and 71.2 cm; see Fig. 2d, e). The magnitude of D and D * are represented by their root-mean squares, which are close to their standard deviations. Histograms for D i,j (Fig. 2d) and D * i,j (Fig. 2e) are both non-Gaussian and negatively skewed. The major difference between the two is the single modal form of D i,j with a peak at around 20 cm and the bimodal form of D * i,j with a high peak at around 30 cm and a low peak at −140 cm. The statistical parameters are different, e.g., mean value and standard deviation are 0.524 and 54.9 cm for D i,j , and −3.84 and 71.2 cm for D * i,j . Skewness and kurtosis are −0.83 and 3.01 for D i,j , and −0.87 and 2.80 for D * i,j .
Horizontal gradients of the DOT, (∂D i,j /∂x, ∂D i,j /∂y) and (∂D * i,j /∂x, ∂D * i,j /∂y), have oceanographic significance (related to the geostrophic currents). They are calculated using the central difference scheme at inside-domain grid points and the first order forward/backward difference scheme at grid points next to the boundary. The difference in global ∂D i,j /∂x (Fig. 3a) and ∂D * i,j /∂x (Fig. 3b) is evident with much smaller-scale structures in ∂D * i,j /∂x. The difference between the two gradients ( Fig. 3c), has the same order of magnitude as the gradients themselves with the relative RRMS of (∂D/∂x), which implies that the non-surface latitudinal geostrophic current component of the second type of MDOT has the same order of magnitude as the surface latitudinal geostrophic current component of the first type of DOT. Histograms for ∂D i,j /∂x (Fig. 3d) and ∂D * i,j /∂x (Fig. 3e) are near symmetric with mean values around (−1.29, −0.78)×10 −8 and standard deviations (2.69, 4.95)×10 −7 . The standard deviation of ∂D * i,j /∂x is almost twice that of ∂D i,j /∂x. Similarly, the difference in global ∂D i,j /∂y (Fig. 4a) and ∂D * i,j /∂y (Fig. 4b) is evident with much smaller-scale structures in ∂D * i,j /∂y. The difference between the two gradients ( Fig. 4c), has the same order of magnitude as the gradients themselves with the relative root-mean square (RRMS) of (∂D/∂y), which implies that the non-surface zonal geostrophic current component of the second type of MDOT has the same order of magnitude as the surface zonal geostrophic current component of the first type of DOT. Histograms for ∂D i,j /∂y (Fig. 4d) and ∂D * i,j /∂y (Fig. 4e) are also near symmetric with the mean values around (2.32, 1.18)×10 −7 and standard deviations (1.20, 2.44)×10 −6 . The standard deviation of ∂D * i,j /∂y is twice that of ∂D i,j /∂y. The denominators of Eqs. (41) and (43) represent the magnitude of the horizontal gradients of the first type of DOT.

Conclusions
The change in marine geoid from classically defined (first type, stand-alone concept in oceanography) to satellitedetermined (second type, stand-alone concept in marine geodesy) largely affects oceanography. With the classically defined marine geoid (average level of SSH if the water is at rest), the horizontal gradients of the first type of DOT represent the absolute surface geostrophic currents. With the satellite-determined (second type) marine geoid by Eq. (6), the horizontal gradients of the second type of MDOT do not represent the absolute surface geostrophic currents. The difference between the two types of DOT represents an additional component to the absolute surface geostrophic currents.
With conservation of potential vorticity, geostrophic balance represents the minimum energy state in an ocean basin where the mechanical energy is conserved. A new governing elliptic equation of the first type of DOT is derived with water depth (H ) in the coefficients and the three-dimensional temperature and salinity in the forcing function. This governing elliptic equation is well posed. Continuation of the geoid from land to ocean leads to an inhomogeneous Dirichlet boundary condition.
Difference between the two types of DOT is evident with a RRMS difference of 38.6 %. Horizontal gradients (representing geostrophic currents) of the two types of DOTs are different with much smaller-scale structures in the second type absolute DOT. The RRMS difference is near 1.0 in both (x, y) components of the DOT gradient, which implies that the non-absolute surface geostrophic currents identified from the second type has the same order of magnitude as the absolute surface geostrophic currents identified by the first type of DOT. The notable difference between the two types of DOT raises more questions in oceanography and marine geodesy: is there any theoretical foundation to connect the classical marine geoid (stand-alone concept in oceanography using the principle of surface geostrophic currents without g) to the satellite-determined marine geoid (stand-alone concept in marine geodesy using g without the principle of surface geostrophic currents)? How can the satellite-determined marine geoid using the gravity anomaly ( g) be conformed to the basic physical oceanography principle of surface geostrophic currents? What is the interpretation of the horizontal gradients of the second type of MDOT (D * )? Is there any evidence or theory to show (u g (N * ) = 0, v g (N * ) = 0) similar to Eq. (13)? More observational and theoretical studies are needed in order to solve those problems. The main challenge for oceanographers is how to use the satellite altimetry observed SSH such as the Surface Water and Ocean Topography (SWOT, https://swot.jpl.nasa.gov/, last access: 31 July 2018) to infer the ocean general circulations at the surface. A new theoretical framework rather than the geostrophic constraint needs to be established.
The GOCE satellite-determined data-only geoid model is more accurate and with higher resolution than GRACE. The change of GRACE to the GOCE geoid model may increase the accuracy of the calculation of the second type of DOT. However, such a replacement does not solve the fundamental problem presented here, i.e., incompatibility between satellite-determined marine geoid using the gravity anomaly ( g) and the classical marine geoid (mean SSH when the water at rest) on the basis of the basic physical oceanography principle of surface geostrophic currents.
Finally, the mathematical framework described here (i.e., the elliptic Eq. 32 with boundary condition Eq. 37) may lead to a new inverse method for calculating three-dimensional absolute geostrophic velocity from temperature and salinity fields since the surface absolute geostrophic velocity is the solution of Eq. (32). This will be a useful addition to the existing β-spiral method (Stommel and Schott, 1977), box model (Wunsch, 1978), and P -vector method (Chu, 1995;Chu et al., 1998Chu et al., , 2000. P. C. Chu: Technical Note: Two types of absolute dynamic ocean topography 9 Appendix A: Geostrophic balance as a minimum energy state in an energy conserved basin In large-scale motion (small Rossby number) with the Boussinesq approximation, the linearized potential vorticity ( ) is given by where, ρ 0 = 1025 kg m −3 is the characteristic density. Without the frictional force and 0 horizontally integrated buoyancy flux at the surface and bottom, the energy (including kinetic and available potential energies) is conserved in a threedimensional ocean basin (V ) The two terms of J are kinetic energy and available potential energy.
To show the geostrophic balance taking the minimum energy state for a given linear PV (see Eq. A1), the constraint is incorporated by extremizing the integral (see also in Vallis, 1992;Chu, 2018) where µ(x, y, z) is the Lagrange multiplier, which is a function of space. If it were a constant, the integral would merely extremize energy subject to a given integral of PV, and rearrangement of PV would leave the integral unaltered. Extremization of the integral Eq. (A4) gives the three Euler-Lagrange equations, where K is in the integrand appearing in Eq. (A4). Substitution of K into Eqs. (A5), (A6), (A7) leads to which shows that (u, v) = (u g , v g ) have the minimum energy state.
With the given boundary condition (Eq. 37, see Sect. 4) and forcing function (Eq. B4), the first type of DOT at the grid points can be computed from steps n to n+1, 2 1 + cos 2 φ j D (n+1) i,j Such iteration continues until the RRMS difference reaches the criterion, where M = 38 877, is the total number of the grid points on the ocean surface.