3.4.3 Earth Rotation Model
Our Earth is the mother of all clocks. It has given us the time units of days, hours, minutes, and seconds that govern our lives. Not until the development of atomic clocks in the mid‐twentieth century were we able to observe the imperfections in our Earth‐clock due to shifting mass distributions caused by tectonic events, the transfers of angular momentum with the atmosphere and hydrosphere, and the gradual slowdown from the transfer of energy and angular momentum within the Earth–moon system due to tides.4 Despite these, we continue to use Earth rotation as our primary time reference, adding or subtracting leap seconds to atomic clocks to keep them synchronized to the rotation of the Earth. These time variations are significant for GNSS navigation, but not for inertial navigation.
The value of earthrate in the World Geodetic System 1984 (WGS 84) Earth model is
3.4.4 Gravity Models
Because an INS operates in a world with gravitational accelerations it is unable to sense and unable to ignore, it must use a reasonably faithful model of gravity.
Gravity models for the Earth include centrifugal acceleration due to rotation of the Earth as well as true gravitational accelerations due to the mass distribution of the Earth, but they do not generally include oscillatory effects such as tidal variations.
3.4.4.1 Gravitational Potential
Gravitational potential of a unit of mass is defined to be zero at a point infinitely distant from all massive bodies and to decrease toward massive bodies such as the Earth. That is, a point at infinity is the reference point for gravitational potential.
In effect, the gravitational potential at a point in or near the Earth is defined by the potential energy lost per unit of mass falling to that point from infinite altitude. In falling from infinity, potential energy is converted to kinetic energy,
3.4.4.2 Gravitational Acceleration
Gravitational acceleration is the negative gradient of gravitational potential. Potential is a scalar function, and its gradient is a vector. Because gravitational potential increases with altitude, its gradient points upward and the negative gradient points downward.
3.4.4.3 Equipotential Surfaces
An equipotential surface is a surface of constant gravitational potential. If the ocean and atmosphere were not moving, then the surface of the ocean at static equilibrium would be an equipotential surface. Mean sea level is a theoretical equipotential surface obtained by time‐averaging the dynamic effects. Orthometric altitude is measured along the (curved) plumbline.
WGS84 Ellipsoid
The WGS84 Earth model approximates mean sea level (an equipotential surface) by an ellipsoid of revolution with its rotation axis coincident with the rotation axis of the Earth, its center at the center of mass of the Earth, and its prime meridian through Greenwich. Its semimajor axis (equatorial radius) is defined to be 6 378 137 m, and its semiminor axis (polar radius) is defined to be 6 356 752.3142 m.
Geoid Models
Geoids are approximations of mean sea‐level orthometric height with respect to a reference ellipsoid. Geoids are defined by additional higher‐order shapes, commonly modeled by spherical harmonics of height deviations from an ellipsoid, as illustrated in Figure 3.8. There are many geoid models based on different data, but the more recent, most accurate models depend heavily on GPS data. Geoid heights deviate from reference ellipsoids by tens of meters, typically.
The WGS84 geoid heights vary about
3.4.4.4 Longitude and Latitude Rates
The second integral of acceleration in locally level coordinates should result in the estimated vehicle position. This integral is somewhat less than straightforward when longitude and latitude are the preferred horizontal location variables.
The rate of change of vehicle altitude equals its vertical velocity, which is the first integral of net (i.e. including gravity) vertical acceleration. The rates of change of vehicle longitude and latitude depend on the horizontal components of vehicle velocity, but in a less direct manner. The relationship between longitude and latitude rates and east and north velocities is further complicated by the oblate shape of the Earth.
Figure 3.8 Equipotential surface models for Earth.
The rates at which these angular coordinates change as the vehicle moves tangent to the surface will depend upon the radius of curvature of the reference surface model. Radius of curvature can depend on the direction of travel, and for an ellipsoidal model there is one radius of curvature for north–south motion and another radius of curvature for east–west motion.
Meridional Radius of Curvature
The radius of curvature for north–south motion is called the “meridional” radius of curvature, because north–south travel is along a meridian (i.e. line of constant longitude). For an ellipsoid of revolution, all meridians have the same shape, which is that of the ellipse that was rotated to produce the ellipsoidal surface model. The tangent circle with the same radius of curvature as the ellipse is called the “osculating circle” (osculating means “kissing”). As illustrated in Figure 3.10 for an oblate Earth model, the radius of the meridional osculating circle is smallest where the geocentric radius is largest (at the equator), and the radius of the osculating circle is largest where the geocentric radius is smallest (at the poles). The osculating circle lies inside or on the ellipsoid at the equator and outside or on the ellipsoid at the poles and passes through the ellipsoid surface for latitudes in between.