A data-driven approach to local gravity field modelling using spherical radial basis functions

We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.     

Geodetic intersection on the ellipsoid

Through each of two known points on the ellipsoid a geodesic is passing in a known azimuth. We solve the problem of intersection of the two geodesics. The solution for the latitude is obtained as a closed formula for the sphere plus a small correction, of the order of the eccentricity of the ellipsoid, which is determined by numerical integration. The solution is iterative. Once the latitude is obtained, the longitude is determined without iteration.     

Spin rate estimation of sounding rockets using GPS wind-up

Carrier phase wind-up is a well-known effect that arises from the relative rotation between a transmitting and receiving antenna. In GPS measurements at L1 frequency, this effect translates into an error of 19.029 cm per full relative rotation of antennas. Since this effect is independent of the satellite elevation for pure rotation about the antenna boresight axis, it is usually absorbed by the clock estimation in navigation algorithms. Therefore, the impact of wind-up is usually neglected for applications that do not require accuracies to the cm level like RTK. However, in receiving platforms with high rotation rate, the accumulated wind-up value can be important and actually be larger than receiver noise or even ionospheric variations. Therefore, in such scenarios, the wind-up contribution can be isolated and used as a source of information to compute the spin rate of such platforms using an appropriate combination of GPS observables. This work shows some results of a coarse, yet simple, approach to monitor the rotation angle and spin-rate of spin stabilized sounding rockets flown by DLR.     

An approach to aid INS using time-differenced GPS carrier phase (TDCP) measurements

This paper presents a simple and effective approach that incorporates single-frequency, L1 time-differenced GPS carrier phase (TDCP) measurements without the need of ambiguity resolution techniques and the complexity to accommodate the delayed-state terms. Static trial results are included to illustrate the stochastic characteristics and effectiveness of the TDCP measurements in controlling position error growth. The formulation of the TDCP observation model is also described in a 17-state tightly-coupled GPS/INS iterative, extended Kalman filter (IEKF) approach. Preliminary land vehicle trial results are also presented to illustrate the effectiveness of the TDCP which provides sub-meter positional accuracies when operating for more than 10 min.     

On the Use of CloudSat and MODIS Data for Estimating Hurricane Intensity

This letter presents preliminary results concerning the use of new observations from the A-Train Constellation for testing a new technique of remotely sensing hurricane intensity from space based on modeling a hurricane as a balanced, convectively neutral vortex. The key observational requirements are simultaneous, accurate measurements of cloud-top height, cloudtop temperature, and cloud profiling information across the center of the storm, although there are ways to bypass the need for cloud-top temperature. In this letter, the Moderate Resolution Imaging Spectroradiometer onboard Aqua provides an estimation of the cloud-top temperature, and the near-simultaneous CloudSat observations provide the essential cloud-top height and cloud profiling information. Initial results indicate that the new technique is a promising method for estimating storm intensity when compared post facto to the best track database. Potential uncertainties and room for further refinement of the technique are discussed.     

Using a Grid-Based Filter to Solve InSAR Phase Unwrapping

This letter presents a phase-unwrapping (PU) algorithm for synthetic aperture radar interferometry based on a grid-based filter. The proposed PU algorithm, which is based on state-space techniques, simultaneously performs noise filtering and PU. The formulation of this technique provides independence from noise statistics and is not constrained by the nonlinearity of the problem. Results obtained with synthetic data show a significant improvement with respect to other conventional PU algorithms in some situations.     

Income distribution dynamics and cross-region convergence in Europe

This paper presents a continuous version of the model of distribution dynamics to analyse the transition dynamics and implied long-run behaviour of the EU-27 NUTS-2 regions over the period 1995–2003. It departs from previous research in two respects: first, by introducing kernel estimation and three-dimensional stacked conditional density plots as well as highest density regions plots for the visualisation of the transition function, based on Hyndman et al. (J Comput Graph Stat 5(4):315–336, 1996), and second, by combining Getis’ spatial filtering view with kernel estimation to explicitly account for the spatial dimension of the growth process. The results of the analysis indicate a very slow catching-up of the poorest regions with the richer ones, a process of shifting away of a small group of very rich regions, and highlight the importance of geography in understanding regional income distribution dynamics.