The North American datum of 1983: Project methodology and execution

A new adjustment of the geodetic control networks in North America has been completed, resulting in a new continental datum—the North American Datum of 1983 (NAD 83). The establishment ofNAD 83 was the result of an international project involving the National Geodetic Survey of the United States, the Geodetic Survey of Canada, and the Danish Geodetic Institute (responsible for surveying in Greenland). The geodetic data in Mexico and Central America were collected by the Inter American Geodetic Survey and validated by the Defense Mapping Agency Hydrographic/Topographic Center. The fundamental task ofNAD 83 was a simultaneous least squares adjustment involving 266,436 stations in the United States, Canada, Mexico, and Central America. The networks in Greenland, Hawaii, and the Caribbean islands were connected to the datum through Doppler satellite and Very Long Baseline Interferometry (VLBI) observations. The computations were performed with respect to the ellipsoid of the Geodetic Reference System of 1980. The ellipsoid is positioned in such a way as to be geocentric, and its axes are oriented by the Bureau International de l'Heure Terrestrial System of 1984. The mathematical model for theNAD readjustment was the height-controlled three-dimensional system. The least squares adjustment involved 1,785,772 observations and 928,735 unknowns. The formation and solution of the normal equations were carried out according to the Helmert block method. [Authors' note:This article is a condensation of the final report of the NAD 83 project. The full report (Schwarz,1989) contains a more complete discussion of all the topics.]     

The geopotential value <Emphasis Type="Italic">W</Emphasis><Subscript>0</Subscript> for specifying the relativistic atomic time scale and a global vertical reference system

The TOPEX/Poseidon (T/P) satellite alti- meter mission marked a new era in determining the geopotential constant W ₀. On the basis of T/P data during 1993–2003 (cycles 11–414), long-term variations in W ₀ have been investigated. The rounded value W ₀ = 62636856.0 ± 0.5) m ² s ⁻² has already been adopted by the International Astronomical Union for the definition of the constant L _G = W ₀/c ² = 6.969290134 × 10⁻¹⁰ (where c is the speed of light), which is required for the realization of the relativistic atomic time scale. The constant L _G , based on the above value of W ₀, is also included in the 2003 International Earth Rotation and Reference Frames Service conventions. It has also been suggested that W ₀ is used to specify a global vertical reference system (GVRS). W ₀ ensures the consistency with the International Terrestrial Reference System, i.e. after adopting W ₀, along with the geocentric gravitational constant (GM), the Earth’s rotational velocity (ω) and the second zonal geopotential coefficient (J ₂) as primary constants (parameters), then the ellipsoidal parameters (a,α) can be computed and adopted as derived parameters. The scale of the International Terrestrial Reference Frame 2000 (ITRF2000) has also been specified with the use of W ₀ to be consistent with the geocentric coordinate time. As an example of using W ₀ for a GVRS realization, the geopotential difference between the adopted W ₀ and the geopotential at the Rimouski tide-gauge point, specifying the North American Vertical Datum 1988 (NAVD88), has been estimated.     

现代大地测量参考系统

概述现代大地测量参考系统的定义及不同参考系统之间的关系.主要讨论我国当代大地测量界常使用的3种用以表示几何位置的参考系统1980年国家大地坐标系、全球大地测量系统 1984 (WGS 84)、国际地球参考系统(ITRS),和一种用以表示物理位置,高程的参考系统1985年国家高程基准.并讨论大地测量中框架和基准的概念.     

A bias-free geodetic boundary value problem approach to height datum unification

A geodetic boundary value problem (GBVP) approach has been formulated which can be used for solving the problem of height datum unification. The developed technique is applied to a test area in Southwest Finland with approximate size of 1.5° × 3° and the bias of the corresponding local height datum (local geoid) with respect to the geoid is computed. For this purpose the bias-free potential difference and gravity difference observations of the test area are used and the offset (bias) of the height datum, i.e., Finnish Height Datum 2000 (N2000) fixed to Normaal Amsterdams Peil (NAP) as origin point, with respect to the geoid is computed. The results of this computation show that potential of the origin point of N2000, i.e., NAP, is (62636857.68 ± 0.5) (m²/s²) and as such is (0.191 ± 0.003) (m) under the geoid defined by W ₀ = 62636855.8 (m²/s²). As the validity test of our methodology, the test area is divided into two parts and the corresponding potential difference and gravity difference observations are introduced into our GBVP separately and the bias of height datums of the two parts are computed with respect to the geoid. Obtaining approximately the same bias values for the height datums of the two parts being part of one height datum with one origin point proves the validity of our approach. Besides, the latter test shows the capability of our methodology for patch-wise application.     

Molodensky potential telluroid based on a minimum-distance map. Case study: the quasi-geoid of East Germany in the World Geodetic Datum 2000

 A potential-type Molodensky telluroid based upon a minimum-distance mapping is derived. With respect to a reference potential of Somigliana–Pizzetti type which relates to the World Geodetic Datum 2000, it is shown that a point-wise minimum-distance mapping of the topographical surface of the Earth onto the telluroid surface, constrained to the gauge W(P)=u(p), leads to a system of four nonlinear normal equations. These normal equations are solved by a fast Newton–Raphson iteration. Received: 7 February 2000 / Accepted: 23 October 2001     

W0: an estimate in the Finnish Height Datum N60, epoch 1993.4, from twenty-five GPS points of the Baltic Sea Level Project

Based upon a data set of 25 points of the Baltic Sea Level Project, second campaign 1993.4, which are close to mareographic stations, described by (1) GPS derived Cartesian coordinates in the World Geodetic Reference System 1984 and (2) orthometric heights in the Finnish Height Datum N60, epoch 1993.4, we have computed the primary geodetic parameter W ₀(1993.4) for the epoch 1993.4 according to the following model. The Cartesian coordinates of the GPS stations have been converted into spheroidal coordinates. The gravity potential as the additive decomposition of the gravitational potential and the centrifugal potential has been computed for any GPS station in spheroidal coordinates, namely for a global spheroidal model of the gravitational potential field. For a global set of spheroidal harmonic coefficients a transformation of spherical harmonic coefficients into spheroidal harmonic coefficients has been implemented and applied to the global spherical model OSU 91A up to degree/order 360/360. The gravity potential with respect to a global spheroidal model of degree/order 360/360 has been finally transformed by means of the orthometric heights of the GPS stations with respect to the Finnish Height Datum N60, epoch 1993.4, in terms of the spheroidal “free-air” potential reduction in order to produce the spheroidal W ₀(1993.4) value. As a mean of those 25 W ₀(1993.4) data as well as a root mean square error estimation we computed W ₀(1993.4)=(6 263 685.58 ± 0.36) kgal × m. Finally a comparison of different W ₀ data with respect to a spherical harmonic global model and spheroidal harmonic global model of Somigliana-Pizetti type (level ellipsoid as a reference, degree/order 2/0) according to The Geodesist's Handbook 1992 has been made. Received: 7 November 1996 / Accepted: 27 March 1997     

1980国家坐标系向某市地方独立坐标系的转换

改变国家大地基准参考椭球的参数,以新的椭球作为投影基准,选择合适的中央子午线,通过坐标相似变换的方法,实现1980国家坐标系(1980西安坐标系)向某市地方独立坐标系的转换,并对转换的内符合精度和外符合精度进行了对比分析。     

World Geodetic Datum 2000

 Based on the current best estimates of fundamental geodetic parameters {W ₀,GM,J ₂,Ω} the form parameters of a Somigliana-Pizzetti level ellipsoid, namely the semi-major axis a and semi-minor axis b (or equivalently the linear eccentricity ) are computed and proposed as a new World Geodetic Datum 2000. There are six parameters namely the four fundamental geodetic parameters {W ₀,GM,J ₂,Ω} and the two form parameters {a,b} or {a,ɛ}, which determine the ellipsoidal reference gravity field of Somigliana-Pizzetti type constraint to two nonlinear condition equations. Their iterative solution leads to best estimates a=(6 378 136.572±0.053)m, b=(6 356 751.920 ± 0.052)m, ɛ=(521 853.580±0.013)m for the tide-free geoide of reference and a=(6 378 136.602±0.053)m, b=(6 356 751.860±0.052)m, ɛ=(521 854.674 ± 0.015)m for the zero-frequency tide geoid of reference. The best estimates of the form parameters of a Somigliana-Pizzetti level ellipsoid, {a,b}, differ significantly by −0.39 m, −0.454 m, respectively, from the data of the Geodetic Reference System 1980. Received: 1 February 1999 / Accepted: 31 August 1999     

The GEOID96 high-resolution geoid height model for the United States

The 2 arc-minute × 2 arc-minute geoid model (GEOID96) for the United States supports the conversion between North American Datum 1983 (NAD 83) ellipsoid heights and North American Vertical Datum 1988 (NAVD 88) Helmert heights. GEOID96 includes information from global positioning system (GPS) height measurements at optically leveled benchmarks. A separate geocentric gravimetric geoid, G96SSS, was first calculated, then datum transformations and least-squares collocation were used to convert from G96SSS to GEOID96. Fits of 2951 GPS/level (ITRF94/NAVD 88) benchmarks to G96SSS show a 15.1-cm root mean square (RMS) around a tilted plane (0.06 ppm, 178^∘ azimuth), with a mean value of −31.4 cm (15.6-cm RMS without plane). This mean represents a bias in NAVD 88 from global mean sea level, remaining nearly constant when computed from subsets of benchmarks. Fits of 2951 GPS/level (NAD 83/NAVD 88) benchmarks to GEOID96 show a 5.5-cm RMS (no tilts, zero average), due primarily to GPS error. The correlated error was 2.5 cm, decorrelating at 40 km, and is due to gravity, geoid and GPS errors. Differences between GEOID96 and GEOID93 range from −122 to +374 cm due primarily to the non-geocentricity of NAD 83. Received: 28 July 1997 / Accepted: 2 September 1998     

Accuracy analysis for the NAD 83 geodetic reference system

The North American Datum of 1983 (NAD 83) provides horizontal coordinates for more than 250,000 geodetic stations. These coordinates were derived by a least squares adjustment of existing terrestrial and space-based geodetic data. For pairs of first order stations with interstation distances between 10km and 100km, therms discrepancy between distances derived fromNAD 83 coordinates and distances derived from independentGPS data may be suitably approximated by the empirical rulee=0.008 K^0.7 where e denotes therms discrepancy in meters and K denotes interstation distance in kilometers. For the same station pairs, therms discrepancy in azimuth may be approximated by the empirical rule e=0.020 K^0.5. Similar formulas characterize therms discrepancies for pairs involving second and third order stations. Distance and orientation accuracies, moreover, are well within adopted standards. While these expressions indicate that the magnitudes of relative positional accuracies depend on station order, absolute positional accuracies are similar in magnitude for first, second, and third order stations. Adjustment residuals reveal a few local problems with theNAD 83 coordinates and with the weights assigned to certain classes of observations.     

Ellipsoidal spectral properties of the Earth’s gravitational potential and its first and second derivatives

The spectral analysis of the Earth’s gravitational potential, its first and second derivatives is performed in spherical/ellipsoidal harmonics relative to the International Reference Sphere/International Reference Ellipsoid. The highlights of the diagrammatic approach are: (1) Up- and downward continuation of incremental gravity gradients, (2) Downward continuation of incremental gravity gradients (four tensor-valued harmonic functions) to the incremental gravity potential on the International Reference Figure, (3) Direct conversion of external incremental gravity gradients to geoidal undulation by means of the spherical/ellipsoidal Bruns Formula. The International Reference Ellipsoid was chosen as an equipotential surface in the Somigliana-Pizzetti reference potential field.     

新疆现代测绘基准体系建设的探讨

测绘是国民经济和社会发展的一项前期性、基础性和公益性工作,事关国家主权、国家安全和民族尊严,是推动信息化建设、促进经济社会可持续发展的基础保障[1]。全文提出新疆建设现代化的测绘基准体系的重要意义,阐述了新疆现代测绘基准体系建设结构和建设方案。     

Accuracy of OPUS solutions for 1- to 4-h observing sessions

We processed 30 consecutive days of Global Positioning System (GPS) data using the On-line Positioning Users Service (OPUS) provided by the National Geodetic Survey (NGS) to determine how the accuracy of derived three-dimensional positional coordinates depends on the length of the observing session T, for T ranging from 1 h to 4 h. We selected five Continuously Operating Reference Stations (CORS), distributed widely across the USA, and processed the GPS data for each with corresponding data from three of its nearby CORS. Our results support the current OPUS policy that recommends using a minimum of 2 h of static GPS data. In particular, 2 h of data yielded a root mean square error of 0.8, 2.1, and 3.4 cm in the north, east, and up components of the derived positional coordinates, respectively. Results drastically improve for solutions containing 3 h or more of GPS data.     

Geodesy and relativity

Relativity, or gravitational physics, has widely entered geodetic modelling and parameter determination. This concerns, first of all, the fundamental reference systems used. The Barycentric Celestial Reference System (BCRS) has to be distinguished carefully from the Geocentric Celestial Reference System (GCRS), which is the basic theoretical system for geodetic modelling with a direct link to the International Terrestrial Reference System (ITRS), simply given by a rotation matrix. The relation to the International Celestial Reference System (ICRS) is discussed, as well as various properties and relevance of these systems. Then the representation of the gravitational field is discussed when relativity comes into play. Presently, the so-called post-Newtonian approximation to GRT (general relativity theory) including relativistic effects to lowest order is sufficient for practically all geodetic applications. At the present level of accuracy, space-geodetic techniques like VLBI (Very Long Baseline Interferometry), GPS (Global Positioning System) and SLR/LLR (Satellite/Lunar Laser Ranging) have to be modelled and analysed in the context of a post-Newtonian formalism. In fact, all reference and time frames involved, satellite and planetary orbits, signal propagation and the various observables (frequencies, pulse travel times, phase and travel-time differences) are treated within relativity. This paper reviews to what extent the space-geodetic techniques are affected by such a relativistic treatment and where—vice versa—relativistic parameters can be determined by the analysis of geodetic measurements. At the end, we give a brief outlook on how new or improved measurement techniques (e.g., optical clocks, Galileo) may further push relativistic parameter determination and allow for refined geodetic measurements.     

Consistent realization of Celestial and Terrestrial Reference Frames

The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005–2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of \(\varDelta \hbox {UT1}\) results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several \(\upmu \)as. Moreover, the CRF is influenced by up to \(50~\upmu \)as if the station coordinates and EOP are dominated by the satellite techniques.     

National height datum,the Gauss–Listing geoid level value w0 and its time variation 0 (Baltic Sea Level Project: epochs 1990.8, 1993.8, 1997.4)

 A methodology for precise determination of the fundamental geodetic parameter w ₀, the potential value of the Gauss–Listing geoid, as well as its time derivative ₀, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w ₀ and ₀ values (w ₀=62 636 855.75 ± 0.21 m²/s², ₀=−0.0099±0.00079 m²/s² per year, w ₀/&γmacr;=6 379 781.502 m,₀/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m²/s²) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1) the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums of countries around the Baltic Sea. Received: 14 August 2000 / Accepted: 19 June 2001     

顾及测量方式的地球自转参数差异性分析

随着甚长基线干涉测量(VLBI)、卫星激光测距(SLR)、激光测月(LLR)、全球卫星导航系统(GNSS)、多里斯系统(DORIS)等多种空间大地测量手段的使用,地球自转参数(ERP)的测量精度不断提高,为航天器导航、深空探测等诸多领域提供了高精度的国际天球参考系(ICRS)和国际地表参考系统(ITRS)之间的转换参数. 以国际地球自转与参考系服务发布的C04序列为基础序列,选取500天ERP序列,分析不同测量手段得到的ERP数据的误差分布情况,为研究利用不同数据之间的一致性进行精度检核的可行性及精度水平提供数据基础,同时也为ERP预报提供更多的数据选择.      

FJCORS的构建及其在控制测量中的应用

具有实时定位服务功能的连续运行卫星定位服务系统（CORS）是当代GPS发展的热点之一。从系统组成、技术指标等方面详细介绍了FJCORS,并给出了一种应用FJCORS和区域大地水准面新型控制测量方法。     

Modelling position dependent errors in gravimetric deflections of the vertical

Comparisons of gravimetric and astrogeodetic deflections of the vertical in the Australian region indicate that the former are affected by position dependent systematic errors, even after orientation onto the Australian Geodetic Datum. These are probably due to errors in the predicted mean anomalies for gravimetrically unsurveyed oceanic regions to the east, south and west of the continent. Deflection component residuals (astrogeodetic minus oriented gravimetric) at 83 control stations are made the observables in a set of observation equations, based on the Vening Meinesz equations, from which pseudocorrections to the mean anomalies for a set of arbitrarily selected surface elements are computed. These pseudocorrections compensate for prediction errors in much larger unsurveyed regions. Their effects on individual deflection components are calculated using the Vening Meinesz equations. Statistical tests indicate that pseudocorrections computed for four large offshore elements and six smaller elements in unsurveyed areas produce corrections to the gravimetric deflections which make the ξ and η components in seconds of arc consistent with normally distributed populations N (0.00, 0.70²).     

格尔木市高精度数字高程基准模型研究

大地水准面（数字高程基准）为国家高程基准的建立与维持提供了全新的思路。然而,受限于地形、重力数据等原因,高原地区高精度数字高程基准模型的建立一直是大地测量领域的难题。本文以格尔木地区为例,探讨了高原地区高精度数字高程基准模型的建立方法。首先,基于重力和地形数据,由第二类Helmert凝集法计算了格尔木重力似大地水准面。在计算中,考虑到高原地形对大地水准面模型的影响,采用了7.5″×7.5″分辨率和高精度的地形数据来恢复大地水准面短波部分的方法,以提高似大地水准面的精度。然后,利用球冠谐调和分析方法将GNSS水准与重力似大地水准面联合,建立了格尔木高精度数字高程基准模型。与实测的67个高精度GNSS水准资料比较,重力似大地水准面的外符合精度为3.0 cm,数字高程基准模型的内符合精度为2.0 cm。