An improved equation of latitude and a global system of graticule distance coordinates

Two innovations are presented for coordinate time-series computation. First, an improved solution is given to a century-old problem, that is the non-iterative computation of conventional geodetic (CG: latitude, longitude, height) coordinates from geocentric Cartesian (GC: x, y, z) coordinates. The accuracy is 1 nm for heights < 500 km and < 10^-15 rad for latitude at any point, terrestrial or outer space. This can be 3 orders of magnitude more accurate than published non-iterative methods. Secondly, CG time series are transformed into a practical system of "graticule distance" (GD: easting, northing, height) curvilinear coordinates that, unlike the commonly used system of topocentric Cartesian (TC: east, north, up) coordinates, is global in nature without arbitrary specification of GC reference coordinates for every geodetic station. Since 2011, Nevada Geodetic Laboratory has publicly produced time series for 20,000 GPS stations in GD form that have been cited by hundreds of studies. The GD system facilitates direct comparison of positions for co-located stations. Users of GD time series are able: (1) to resolve different historical station names that have been assigned to the same physical benchmark and (2) to resolve different physical benchmarks that have been assigned the same name. This benefits historical reconstruction of benchmark occupation and local site tie analysis for reference frame integrity. GD coordinates have archival value, in that inversion back to GC coordinates is practically exact. For geodetic stations, GD time series closely emulate TC time series with rates agreeing to 0.01 mm/yr, and so can be used interchangeably.