Investigations on the hierarchy of reference frames in geodesy and geodynamics

Problems related to reference directions were investigated. Space and time variant angular parameters are illustrated in hierarchic structures or towers. Using least squares techniques, model towers of triads are presented which allow the formation of linear observation equations. Translational and rotational degrees of freedom (origin and orientation) are discussed along with and the notion of length and scale degrees of freedom. According to the notion of scale parallelism, scale factors with respect to a unit length are given. Three-dimensional geodesy was constructed from the set of three base vectors (gravity, earth-rotation and the ecliptic normal vector). Space and time variations are given with respect to a polar and singular value decomposition or in terms of changes in translation, rotation, deformation (shear, dilatation or angular and scale distortions)

Positioning systems in Minkowski space-time: Bifurcation problem and observational data

In the framework of relativistic positioning systems in Minkowski space-time, the determination of the inertial coordinates of a user involves the {\em bifurcation problem} (which is the indeterminate location of a pair of different events receiving the same emission coordinates). To solve it, in addition to the user emission coordinates and the emitter positions in inertial coordinates, it may happen that the user needs to know {\em independently} the orientation of its emission coordinates. Assuming that the user may observe the relative positions of the four emitters on its celestial sphere, an observational rule to determine this orientation is presented. The bifurcation problem is thus solved by applying this observational rule, and consequently, {\em all} of the parameters in the general expression of the coordinate transformation from emission coordinates to inertial ones may be computed from the data received by the user of the relativistic positioning system.Comment: 10 pages, 7 figures. The version published in PRD contains a misprint in the caption of Figure 3, which is here amende

Is it possible to test directly General Relativity in the gravitational field of the Moon?

In this paper the possibility of measuring some general relativistic effects in the gravitational field of the Moon via selenodetic missions, with particular emphasis to the future Japanese SELENE mission, is investigated. For a typical selenodetic orbital configuration the post-Newtonian Lense-Thirring gravitomagnetic and the Einstein's gravitoelectric effects on the satellites orbits are calculated and compared to the present-day orbit accuracy of lunar missions. It turns out that for SELENE's Main Orbiter, at present, the gravitoelectric periselenium shift, which is the largest general relativistic effect, is 1 or 2 orders of magnitude smaller than the experimental sensitivity. The systematic error induced by the mismodelled classical periselenium precession due to the first even zonal harmonic J2 of the Moon's non-spherical gravitational potential is 3 orders of magnitude larger than the general relativistic gravitoelectric precession. The estimates of this work could be used for future lunar missions having as their goals relativistic measurements as well.Comment: Latex2e, 7 pages, no figures, ets2000.cls and art12.sty used. Major rewriting in introduction. References adde

Second Order Design of Geodetic Networks by the Simulated Annealing Method

The problem of determining the required precision in observations in order to obtain a desired precision in final parameters, classically known as the second order design problem, is revisited in this paper and proposed to be solved by the simulated annealing method. An example and a flexible implementation in MATLAB are given. © 2011 American Society of Civil Engineers.Baselga Moreno, S. (2011). Second Order Design of Geodetic Networks by the Simulated Annealing Method. Journal of Surveying Engineering. 137(4):167-173. doi:10.1061/(ASCE)SU.1943-5428.0000053167173137

Strain Rate Distribution in South‐Central Tibet From Two Decades of InSAR and GPS

The degree to which deformation and seismicity is focused on major mapped structures remains a key unknown in assessing seismic hazards and testing continental deformation models. Here we combine 208 Global Positioning System (GPS) velocities with 12‐track Interferometric Synthetic Aperture Radar (InSAR) rate maps to form high‐resolution velocity and strain rate fields for south‐central Tibet. Our results show that deformation is not evenly distributed across the region. We find a few zones with high strain rates, most notably the Yutian‐Zhongba strain rate zone. However, the average of the strain rates is similar within and outside the mapped fault zones. In addition, the slip rates are low on all the conjugate strike‐slip faults widespread in central Tibet. The observations are difficult to reconcile with time‐invariant block models or with continuum models that lack mechanisms for strain localization. Our results support arguments that the most robust estimates of seismic hazard should integrate seismicity catalogues, active fault maps, and geodetic strain rate models

The Meissl scheme for the geodetic ellipsoid

We present a variant of the Meissl scheme to relate surface spherical harmonic coefficients of the disturbing potential of the Earth's gravity field on the surface of the geodetic ellipsoid to surface spherical harmonic coefficients of its first- and second-order normal derivatives on the same or any other ellipsoid. It extends the original (spherical) Meissl scheme, which only holds for harmonic coefficients computed from geodetic data on a sphere. In our scheme, a vector of solid spherical harmonic coefficients of one quantity is transformed into spherical harmonic coefficients of another quantity by pre-multiplication with a transformation matrix. This matrix is diagonal for transformations between spheres, but block-diagonal for transformations involving the ellipsoid. The computation of the transformation matrix involves an inversion if the original coefficients are defined on the ellipsoid. This inversion can be performed accurately and efficiently (i.e., without regularisation) for transformation among different gravity field quantities on the same ellipsoid, due to diagonal dominance of the matrices. However, transformations from the ellipsoid to another surface can only be performed accurately and efficiently for coefficients up to degree and order 520 due to numerical instabilities in the inversion

Closed-Form transformation between geodetic and ellipsoidal coordinates

We present formulas for direct closed-form transformation between geodetic coordinates(Φ, λ, h) and ellipsoidal coordinates (β, λ, u) for any oblate ellipsoid of revolution.These will be useful for those dealing with ellipsoidal representations of the Earth's gravityfield or other oblate ellipsoidal figures. The numerical stability of the transformations for nearpolarand near-equatorial regions is also considered

Scalar and vector Slepian functions, spherical signal estimation and spectral analysis

It is a well-known fact that mathematical functions that are timelimited (or spacelimited) cannot be simultaneously bandlimited (in frequency). Yet the finite precision of measurement and computation unavoidably bandlimits our observation and modeling scientific data, and we often only have access to, or are only interested in, a study area that is temporally or spatially bounded. In the geosciences we may be interested in spectrally modeling a time series defined only on a certain interval, or we may want to characterize a specific geographical area observed using an effectively bandlimited measurement device. It is clear that analyzing and representing scientific data of this kind will be facilitated if a basis of functions can be found that are "spatiospectrally" concentrated, i.e. "localized" in both domains at the same time. Here, we give a theoretical overview of one particular approach to this "concentration" problem, as originally proposed for time series by Slepian and coworkers, in the 1960s. We show how this framework leads to practical algorithms and statistically performant methods for the analysis of signals and their power spectra in one and two dimensions, and, particularly for applications in the geosciences, for scalar and vectorial signals defined on the surface of a unit sphere.Comment: Submitted to the 2nd Edition of the Handbook of Geomathematics, edited by Willi Freeden, Zuhair M. Nashed and Thomas Sonar, and to be published by Springer Verlag. This is a slightly modified but expanded version of the paper arxiv:0909.5368 that appeared in the 1st Edition of the Handbook, when it was called: Slepian functions and their use in signal estimation and spectral analysi

Planejamento de redes geodésicas resistentes a múltiplos outliers

Ao se planejar o levantamento de uma rede geodésica, deseja-se que as observações a serem realizadas e as coordenadas dos pontos a serem estimadas atendam critérios de precisão e confiabilidade pré-estabelecidos de acordo com os objetivos do projeto. Na etapa de pré-análise, antes mesmo da coleta das observações, é possível estimar a precisão e confiabilidade da rede, estipulando uma geometria/configuração para a mesma e a precisão esperada para as observações. O objetivo deste artigo é apresentar o planejamento de uma rede geodésica que atenda critérios de precisão e confiabilidade, considerando a possível existência de dois ou mais erros não detectados nas observações, bem como a influência (simultânea) destes erros sobre os parâmetros (coordenadas ajustadas dos vértices). Além da revisão teórica, experimentos foram realizados em uma rede GNSS, onde foram estipulados critérios de precisão e confiabilidade considerando a existência de até duas observações contaminadas por erros (outliers), de maneira simultânea. O planejamento da rede foi feito por meio do método da tentativa e erro. Depois do processamento dos dados e do ajustamento da rede, se verificou que os critérios de precisão e confiabilidade que foram estipulados na etapa de pré-análise foram devidamente obtidos