Statistical analysis of discrepancies in high precision levelling

An investigation was made of the behaviour of the variable (where ρ_ij are the discrepancies between the direct and reverse measurements of the height of consecutive bench marks and theR _ij are their distance apart) in a partial net of the Italian high precision levelling of a total length of about1.400 km. The methods of analysis employed were in general non-parametric individual and cumulative tests; in particular randomness, normality and asymmetry tests were carried out. The computers employed wereIBM/7094/7040. From the results evidence was obtained of the existence of an asymmetry in respect to zero of thex _ij confirming the well-known results given firstly by Lallemand. A new result was obtained from the tests of randomness which put in evidence trends of the mean values of thex _ij and explained some anomalous behaviours of the cumulative discrepancy curves. The extension of this investigation to a broader net possibly covering other national nets would be very useful to get a deeper insight into the behaviour of the errors in high precision levelling. Ad hoc programs for electronic computers are available to accomplish this job quickly. Presented at the 14^th International Assembly of Geodesy (Lucerne, 1967).     

The use of atmospheric dispersion in optical distance measurement

The development of lasers, new electro-optic light modulation methods, and improved electronic techniques have made possible significant improvements in the range and accuracy of optical distance measurements, thus providing not only improved geodetic tools but also useful techniques for the study of other geophysical, meteorological, and astronomical problems. One of the main limitations, at present, to the accuracy of geodetic measurements is the uncertainty in the average propagation velocity of the radiation due to inhomogeneity of the atmosphere. Accuracies of a few parts in ten million or even better now appear feasible, however, through the use of the dispersion method, in which simultaneous measurements of optical path length at two widely separated wavelengths are used to determine the average refractive index over the path and hence the true geodetic distance. The design of a new instrument based on this method, which utilizes wavelengths of6328 ? and3681 ? and3 GHz polarization modulation of the light, is summarized. Preliminary measurements over a5.3 km path with this instrument have demonstrated a sensitivity of3×10 ⁻⁹ in detecting changes in optical path length for either wavelength using1-second averaging, and a standard deviation of3×10 ⁻⁷ in corrected length. The principal remaining sources of error are summarized, as is progress in other laboratories using the dispersion method or other approaches to the problem of refractivity correction.     

Fehlertheoretische Maxwell-Boltzmann-Verteilung

Summary The probability to find an error vector in multiples of the Helmert-Maxwell-Boltzmann point error σ² δ_ij(δ_ij Kronecker symbol) is calculated. It is found that the probability is for σ39%, for2 σ86% and for3 σ99% in two dimensions, for σ20%, for2 σ74% and for3 σ97% in three dimensions. The fundamental Maxwell-Boltzmann-distribution is tabulated0,02 (0,02) 4,50.      

The geopotential to (14, 14) from a combination of satellite and gravimetric data

A set of2261 5°×5° mean anomalies were used alone and with satellite determined harmonic coefficients of the Smithsonian' Institution to determine the geopotential expansion to various degrees. The basic adjustment was carried out by comparing a terrestrial anomaly to an anomaly determined from an assumed set of coefficients. The (14, 14) solution was found to agree within ±3 m of a detailed geoid in the United States computed using1°×1° anomalies for an inner area and satellite determined anomalies in an outer area. Additional comparisons were made to the input anomaly field to consider the accuracy of various harmonic coefficient solutions. A by-product of this investigation was a new γ_E=978.0463 gals in the Potsdam system or978.0326 gals in an absolute system if −13.7 mgals is taken as the Potsdam correction. Combining this value of γ_E withf=1/298.25, KM=3.9860122·10 ²² cm ³ /sec ² , the consistent equatorial radius was found to be6378143 m.     

On the multivariate total least-squares approach to empirical coordinate transformations. Three algorithms

The multivariate total least-squares (MTLS) approach aims at estimating a matrix of parameters, Ξ, from a linear model (Y−E _Y = (X−E _X ) · Ξ) that includes an observation matrix, Y, another observation matrix, X, and matrices of randomly distributed errors, E _Y and E _X . Two special cases of the MTLS approach include the standard multivariate least-squares approach where only the observation matrix, Y, is perturbed by random errors and, on the other hand, the data least-squares approach where only the coefficient matrix X is affected by random errors. In a previous contribution, the authors derived an iterative algorithm to solve the MTLS problem by using the nonlinear Euler–Lagrange conditions. In this contribution, new lemmas are developed to analyze the iterative algorithm, modify it, and compare it with a new ‘closed form’ solution that is based on the singular-value decomposition. For an application, the total least-squares approach is used to estimate the affine transformation parameters that convert cadastral data from the old to the new Israeli datum. Technical aspects of this approach, such as scaling the data and fixing the columns in the coefficient matrix are investigated. This case study illuminates the issue of “symmetry” in the treatment of two sets of coordinates for identical point fields, a topic that had already been emphasized by Teunissen (1989, Festschrift to Torben Krarup, Geodetic Institute Bull no. 58, Copenhagen, Denmark, pp 335–342). The differences between the standard least-squares and the TLS approach are analyzed in terms of the estimated variance component and a first-order approximation of the dispersion matrix of the estimated parameters.     

Preliminary results of the establishment of a marine geodetic point in the Pacific Ocean

In November 1968, a marine geodetic control point was established in the Pacific Ocean at a water depth of6,200 feet. The control point (reference point) consists of three underwater acoustic transponders, two of which are powered with lead-acid batteries and the third with an underwater radioisotope power source “URIPS” with a10- to20- year life expectancy. Four independent measuring techniques (LORAC airborne line-crossing, satellite, ship inertial, and acoustic techniques) were used to measure and determine the coordinates of the control point. Preliminary analysis of the acoustic and airborne data indicates that high accuracies can be achieved in the establishment of geodetic reference points at sea. Geodetic adjustment by the method of variation of coordinates yielded a standard point error of±50 to±66 feet in determining the unknown ship station. The original location of the ship station as determined by shipboard navigation equipment was off by about1,600 feet. Paper previously published in the Proceedings of the Second Marine Geodesy Symposium of the Marine Technology Society.     

The free nonlinear boundary value problem of physical geodesy

Summary Within potential theory of Poisson-Laplace equation the boundary value problem of physical geodesy is classified asfree andnonlinear. For solving this typical nonlinear boundary value problem four different types of nonlinear integral equations corresponding to singular density distributions within single and double layer are presented. The characteristic problem of free boundaries, theproblem of free surface integrals, is exactly solved bymetric continuation. Even in thelinear approximation of fundamental relations of physical geodesy the basic integral equations becomenonlinear because of the special features of free surface integrals.     

Shallow-water loading tides in Japan from superconducting gravimetry

Gravity observations from superconducting gravimeters are used to observe loading effects from shallow-water tides on the Japanese east and west coasts. Non-linear third-diurnal and higher-frequency shallow-water tides are identified in the tide-gauge observations from these coastal areas. The most energetic constituents in the tide gauge observations are also seen in the gravity observations due to their loading effects on the deformation of the Earth. Even though the shallow-water tides at the Japanese east coast have an amplitude of only a few millimetres, they are still able to generate a loading signal at gravity sites located several hundred kilometres inland. In particular, the S₃, S₄ and S₅ solar tides occur in both gravity and tide gauge observations. It is indicated that in other shelf regions with large shallow water tides, the shallow water loading signals account for a significant signal, which should be taken into account.Acknowledgement The authors would like to thank the Hydrographic and Oceanographic Department (Japan Coast Guard), Japan Meteorological Agency and Hokkaido Development Agency for access to the tide-gauge data. Also, the Global Geodynamic Project Information System and Data Center (GGP-ISDC) is acknowledged for providing the gravity data.     

A generalisation of the least-squares method for the adjustment of free networks

Summary The system of normal equations for the adjustment of a free network is a singular one. Therefore, a number of coordinates has to be fixed according to the matrix. The mean square errors and the error ellipses of such an adjustment are dependent on this choice. This paper gives a simple, direct method for the adjustment of free networks, where no coordinates need to be fixed. This is done by minimizing not only the sum of the squares of the weighted errorsV ^T PV=minimun but also the Euclidean norm of the vectorX and of the covariance matrixQ X ^T X=minimum trace (Q)=minimum This last condition is crucial for geodetic problems of this type.     

TRIGONOMETRICAL HEIGHTS AND THE COEFFICIENT OF TERRESTRIAL REFRACTION

Abstract

Mr. A. J. Morley has contributed a series of articles in the Review (E.S.R., iv, 23, 16; iv, 25, 136 and vi, 40, 76) on the adjustment of trigonometrical levels and the evaluation of the coefficient of terrestrial refraction with a view to ascertaining how other Colonies and Dominions deal with these problems. This object is very commendable as several problems concerning both the observational and theoretical sides arise in height determinations, regarding which there is not much guidance in the usual treatises on the subject.     

On the night errors in astronomical determinations

Summary The discrepancy between precision and accuracy in astronomical determinations is usually explained in two ways: on the one hand by ostensible large refraction anomalies and on the other hand by variable instrumental errors which are systematic over a certain interval of time and which are mainly influenced by temperature.In view of the research of several other persons and the author’s own investigations, the authors are of the opinion that the large night-errors of astronomical determinations are caused by variable, systematic instrumental errors dependent on temperature. The influence of refraction anomalies is estimated to be smaller than 0″.1 for most of the field stations. The possibility of determining the anomalous refraction from the observations by the programme given by Prof. Pavlov and Anderson has also been investigated. The precision of the determination of the anomalous refraction is good as long as no other systematic error working in a similar way is present.The results, which are interpreted as an effect of the anomalous refraction by Pavlov and Sergijenko, could also be interpreted as a systematic instrumental error. It is furthermore maintained thatthe latitude and longitude of a field station can be determined in a few hours of one night if the premisses given in [3, p.68]are kept. It has been deplored that the determination of the azimuth has not been given the necessary attention. It is therefore proposed to intensify the research on this problem. The profession has been called upon to acquaint itself better with the valuable possibilities of astronomical determinations and to apply them in a useful and appropriate manner. At the same time, attention has been called to the possibility of improving astronomical determinations with regard to accuracy as well as effectiveness.     

Cartography in Advertising

Abstract

This short paper attempts to pursue an analytical account of the cartography in present day advertising. Particular reference is made to the printed advertisements in newspapers, brochures and periodicals available in Britain.     

A statistical filter for geodetic observations

A method for filtering of geodetic observationwhich leaves the final result normally distributed, is presented. Furthermore, it is shown that if you sacrifice100.a% of all the observations you may be (1−β).100% sure that a gross error of the size Δ is rejected. Another and, may be intuitively, more appealing method is presented; the two methods are compared and it is shown why Method 1 should be preferred to Method 2 for geodetic purposes. Finally the two methods are demonstrated in some numerical examples.     

Physical Considerations of the Storage of Printed Cartographic Materials in Libraries,or, My Life in a Map Case

Abstract

The author offers advice on the storage needs of sheet maps, folded maps, oversized volumes, wall maps, globes, relief models in library and archive settings.     

ATTRACTION OF THE “COMPENSATION”: A NOTE

Abstract

The proof of the attraction of a uniformly thin vertical block in the last issue (No. 24, p. 87) does not appear as satisfactory as the textbook would suggest. It is proposed to attempt here a more direct explanation, independent of substitutional expedients.     

THE ADJUSTMENT OF A FIRST-ORDER LEVELLING NETWORK

Abstract

In the E.S.R. January and April numbers of 1955, Vol. xiii, Nos. 95 and 96, Mr. Hsuan-Loh Su described the “Adjustment of a Level Net by Successive Approximations and by Electrical Analogy”. It does not seem to be as generally known as it should be that the rigid least square solution can be greatly simplified by utilizing the electrical analogy and solving by Kirchhoff's method. The method as detailed below has been in use for over 40 years.     

LEAST SQUARES ADJUSTMENTS OF TRIANGULATION

Abstract

Mr. Rainsford's paper in E.S.R. No. 78 establishes as thoroughly as one can wish by computational example that there is little to choose in relative accuracy between direction adjustment and angle adjustment, and that the deciding factor is really ease and speed of computation.     

CONDITION FOR THE CONSTRUCTION OF A CONFORMAL PROJECTION

Abstract

In the Empire Survey Review for October 1938 (iv, 30, 480) a simple demonstration of the condition to be satisfied for conformal representation was given. This condition may be expressed by the equation w = f(z), where w and z are complex variables representing corresponding points in the w-plane and z-plane respectively, and f(z) is an analytic function of z.     

EXPANSION OF PAPER AS AFFECTING BEARING

Abstract

The problem of the influence of the paper itself on bearings taken off maps or plans has been recently discussed by G. T. M. in this Review (iii, 22, 479). The following alternative method of solving the problem may be of interest; the same figure and symbols are used as in the article mentioned.