Geodetic methods for detecting volcanic unrest: a theoretical approach

In this paper we study the application of different geodetic techniques to volcanic activity monitoring, using theoretical analysis. This methodology is very useful for obtaining an idea of the most appropriate (and efficient) monitoring method, mainly when there are no records of geodetic changes previous to volcanic activity. The analysis takes into account the crustal structure of the area, its geology, and its known volcanic activity to estimate the deformation and gravity changes that might precede eruptions. The deformation model used includes the existing gravity field and vertical changes in the crustal properties. Both factors can have a considerable effect on computed deformation and gravity changes. Topography should be considered when there is a steep slope (greater than 10°). The case study of Teide stratovolcano (Tenerife, Canary Islands), for which no deformation or gravity changes are available, is used as a test. To avoid considering topography, we worked at the lowest level of Las Cañadas and examined a smaller area than the whole caldera or island. Therefore, the results are only a first approach to the most adequate geodetic monitoring system. The methodology can also be applied to active areas where volcanic risk is not associated with a stratovolcano but instead with monogenetic scattered centers, especially when sites must be chosen in terms of detection efficiency or existing facilities. The Canary Islands provide a good example of this type of active volcanic areas, and we apply our model to the island of Lanzarote to evaluate the efficiency of the monitoring system installed at the existing geodynamic station. On this island topography is not important. The results of our study show clearly that the most appropriate geodetic volcano monitoring system is not the same for all different volcanic zones and types, and the particular properties of each volcano/zone must be taken into account when designing each system.     

A non-iterative and non-singular perturbation solution for transforming Cartesian to geodetic coordinates

The Cartesian-to-Geodetic coordinate transformation is re-cast as a minimization algorithm for the height of the Satellite above the reference Earth surface. Optimal necessary conditions are obtained that fix the satellite ground track vector components in the Earth surface. The introduction of an artificial perturbation variable yields a rapidly converging second order power series solution. The initial starting guess for the solution provides 3–4 digits of precision. Convergence of the perturbation series expansion is accelerated by replacing the series solution with a Padé approximation. For satellites with heights < 30,000 km the second-order expansions yields ~mm satellite geodetic height errors and ~10⁻¹² rad errors for the geodetic latitude. No quartic or cubic solutions are required: the algorithm is both non-iterative and non-singular. Only two square root and two arc-tan calculations are required for the entire transformation. The proposed algorithm has been measured to be ~41% faster than the celebrated Bowring method. Several numerical examples are provided to demonstrate the effectiveness of the new algorithm.     

浅谈厦门市大地基准的现代化

GPS定位技术的迅猛发展给测绘行业带来了革命性的变化,建于上世纪90年代初的厦门市首级控制网的升级换代摆上了议事日程。本文探讨厦门市三维、动态、高精度的现代大地基准建设及似大地水准面的精化。     

基于JSCORS的江苏省现代大地坐标框架的建设与维护

本文结合江苏省现代大地坐标框架的建设和JSCORS的运行维护,通过对江苏省原大地坐标框架的情况分析,介绍了江苏大地坐标框架的建设方案.结合区域稳定基点组的确定、框架维护和复测机制的建立、框架维护系统的建设等内容,提出江苏省大地坐标框架维护的解决方案.     

Alternatives to current GPS-RTK services and some implications for CORS infrastructure and operations

The justification for the establishment of CORS networks was initially in support of geodesy and other geoscientific applications at the global and regional level. However, increasingly GPS CORS network operators have sought ways of making their network infrastructure the basis of a profitable business. This has arisen with the introduction of real-time centimeter-level accuracy services, carrier phase-based modes of operation generally referred to as GPS-RTK (real-time kinematic). One approach is to try to recruit a core group of users who are prepared to pay for the GPS-RTK services. But this is only feasible if the number of users, and the fees that are charged, are sufficient to generate a reasonable return-on-investment (ROI). This ROI (or at the very least “cost-recovery”) is important for many network operators in order that they may provide for the maintenance and upgrade of the CORS infrastructure. On the other hand, there are those who advocate that there is no need to recoup CORS investment, that the installed GPS receivers should be seen as public infrastructure in a similar manner to roads, bridges, etc. This paper discusses some new business and operational models for GPS-RTK services. These include models for the establishment and operation of CORS infrastructure, service provision, business cases, and options for value-added services beyond the standard GPS-RTK service. One concept is based on a “client–server” model. Currently GPS-RTK service providers have no control over the quality of the results computed by users. This makes it difficult for them to justify charging for their services. What if instead of broadcasting RTK corrections and placing the onus of obtaining a final solution on the user and his equipment, the user’s coordinates are determined by the service provider? Putting the computational effort on the server side will justify more easily the charging of users for a value-added product: an accurate and quality assured coordinate in the local reference frame. This paper describes the client–server concept as well as possible business models that may underpin such a service model. These models include some derived from mobile telephony and service/hospitality businesses. Furthermore, with the projected proliferation of independent, competitive GPS-RTK services, the concept of a GPS data or service “broker” is worth exploring.     

空间直角坐标计算大地坐标的抛物线逼近法

采用抛物线逼近法求解大地纬度和大地高,先计算空间点在椭球面上的子午面坐标,然后求解点的大地纬度和大地高。     

Altimetry–gravimetry problems with free vertical datum

 The definition and connection of vertical datums in geodetic height networks is a fundamental problem in geodesy. Today, the standard approach to solve it is based on the joint processing of terrestrial and satellite geodetic data. It is generalized to cases where the coverage with terrestrial data may change from region to region, typically across coastlines. The principal difficulty is that such problems, so-called altimetry–gravimetry boundary-value problems (AGPs), do not admit analytical solutions such as Stokes' integral. A numerical solution strategy for the free-datum problem is presented. Analysis of AGPs in spherical and constant radius approximation shows that two of them are mathematically well-posed problems, while the classical AGP-I may be ill posed in special situations. Received: 2 December 1998 / Accepted: 30 November 1999     

基于ArcGIS的80西安坐标系转换到2000国家坐标系的研究

2000中国大地坐标系(CGCS2000)是我国新启用的地心坐标系,但我国目前用以测图及工程规划、设计以及其他用途的地图坐标系一般又都是基于北京54坐标系或1980西安坐标系。如何将这些坐标系转换到2000国家坐标系是当前必须解决的问题。本文就西安80坐标系到2000国家坐标系之间的转换的基础理论和方法进行了研究,给出了基于ArcGIS环境下将1980西安坐标系转换到2000中国大地坐标系的有效解决办法。     

Transforming ellipsoidal heights and geoid undulations between different geodetic reference frames

Transforming height information that refers to an ellipsoidal Earth reference model, such as the geometric heights determined from GPS measurements or the geoid undulations obtained by a gravimetric geoid solution, from one geodetic reference frame (GRF) to another is an important task whose proper implementation is crucial for many geodetic, surveying and mapping applications. This paper presents the required methodology to deal with the above problem when we are given the Helmert transformation parameters that link the underlying Cartesian coordinate systems to which an Earth reference ellipsoid is attached. The main emphasis is on the effect of GRF spatial scale differences in coordinate transformations involving reference ellipsoids, for the particular case of heights. Since every three-dimensional Cartesian coordinate system ‘gauges’ an attached ellipsoid according to its own accessible scale, there will exist a supplementary contribution from the scale variation between the involved GRFs on the relative size of their attached reference ellipsoids. Neglecting such a scale-induced indirect effect corrupts the values for the curvilinear geodetic coordinates obtained from a similarity transformation model, and meter-level apparent offsets can be introduced in the transformed heights. The paper explains the above issues in detail and presents the necessary mathematical framework for their treatment. An erratum to this article can be found at