Antenna phase center calibration for precise positioning of LEO satellites

Phase center variations of the receiver and transmitter antenna constitute a remaining uncertainty in the high precision orbit determination (POD) of low Earth orbit (LEO) satellites using GPS measurements. Triggered by the adoption of absolute phase patterns in the IGS processing standards, a calibration of the Sensor Systems S67-1575-14 antenna with GFZ choke ring has been conducted that serves as POD antenna on various geodetic satellites such as CHAMP, GRACE and TerraSAR-X. Nominal phase patterns have been obtained with a robotic measurement system in a field campaign and the results were used to assess the impact of receiver antenna phase patterns on the achievable positioning accuracy. Along with this, phase center distortions in the actual spacecraft environment were characterized based on POD carrier phase residuals for the GRACE and TerraSAR-X missions. It is shown that the combined ground and in-flight calibration can improve the carrier phase modeling accuracy to a level of 4 mm which is close to the pure receiver noise. A 3.5 cm (3D rms) consistency of kinematic and reduced dynamic orbit determination solutions is achieved for TerraSAR-X, which presumably reflects the limitations of presently available GPS ephemeris products. The reduced dynamic solutions themselves match the observations of high grade satellite laser ranging stations to 1.5 cm but are potentially affected by cross-track biases at the cm-level. With respect to the GPS based relative navigation of TerraSAR-X/TanDEM-X formation, the in-flight calibration of the antenna phase patterns is considered essential for an accurate modeling of differential carrier phase measurements and a mm level baseline reconstruction.

Single-frequency precise point positioning with optimal filtering

The accuracy of standalone GPS positioning improved significantly when Selective Availability was turned off in May 2000. With the availability of various public GPS related products including precise satellite orbits and clocks, and ionosphere maps, a single-frequency standalone user can experience even a further improvement of the position accuracy. Next, using carrier phase measurements becomes crucial to smoothen the pseudorange noise. In this contribution, the most critical sources of error in single-frequency standalone positioning will be reviewed and different approaches to mitigate the errors will be considered. An optimal filter (using also carrier phase measurements) will be deployed. The final approach will then be evaluated in a decently long static test with receivers located in different regions of the world. Kinematic experiments have also been performed in various scenarios including a highly dynamic flight trial. The accuracy, in general, can be confirmed at 0.5 m horizontal and 1 m vertical, with static tests. Ultimate results demonstrate an accuracy close to 2 dm (95%) for the horizontal position components and 5 dm (95%) for the vertical in the flight experiment.

Tracking and orbit determination performance of the GRAS instrument on MetOp-A

The global navigation satellite system receiver for atmospheric sounding (GRAS) on MetOp-A is the first European GPS receiver providing dual-frequency navigation and occultation measurements from a spaceborne platform on a routine basis. The receiver is based on ESA’s AGGA-2 correlator chip, which implements a high-quality tracking scheme for semi-codeless P(Y) code tracking on the L1 and L2 frequency. Data collected with the zenith antenna on MetOp-A have been used to perform an in-flight characterization of the GRAS instrument with focus on the tracking and navigation performance. Besides an assessment of the receiver noise and systematic measurement errors, the study addresses the precise orbit determination accuracy achievable with the GRAS receiver. A consistency on the 5 cm level is demonstrated for reduced dynamics orbit solutions computed independently by four different agencies and software packages. With purely kinematic solutions, 10 cm accuracy is obtained. As a part of the analysis, an empirical antenna offset correction and preliminary phase center correction map are derived, which notably reduce the carrier phase residuals and improve the consistency of kinematic orbit determination results.

Analysis of high-frequency multipath in 1-Hz GPS kinematic solutions

High-frequency multipath would be problematic for studies at seismic or antenna dynamical frequencies as one could mistakenly interpret them as signals. A simple procedure to identify high-frequency multipath from global positioning system (GPS) time series records is presented. For this purpose, data from four GPS base stations are analyzed using spectral analyses techniques. Additional data, such as TEQC report files of L1 pseudorange multipath, are also used to analyze the high-frequency multipath and confirmation of the high-frequency multipath inferred from the phase records. Results show that this simple procedure is effective in identification of high-frequency multipath. The inferred information can aid interpretation of multipath at the GPS site, and is important for a number of reasons. For example, the information can be used to study GPS site selections and/or installations.

Computer algebra solution of the GPS N-points problem

A computer algebra solution is applied here to develop and evaluate algorithms for solving the basic GPS navigation problem: finding a point position using four or more pseudoranges at one epoch (the GPS N-points problem). Using Mathematica 5.2 software, the GPS N-points problem is solved numerically, symbolically, semi-symbolically, and with Gauss–Jacobi, on a work station. For the case of N > 4, two minimization approaches based on residuals and distance norms are evaluated for the direct numerical solution and their computational duration is compared. For N = 4, it is demonstrated that the symbolic computation is twice as fast as the iterative direct numerical method. For N = 6, the direct numerical solution is twice as fast as the semi-symbolic, with the residual minimization requiring less computation time compared to the minimization of the distance norm. Gauss–Jacobi requires eight times more computation time than the direct numerical solution. It does, however, have the advantage of diagnosing poor satellite geometry and outliers. Besides offering a complete evaluation of these algorithms, we have developed Mathematica 5.2 code (a notebook file) for these algorithms (i.e., Sturmfel’s resultant, Dixon’s resultants, Groebner basis, reduced Groebner basis and Gauss–Jacobi). These are accessible to any geodesist, geophysicist, or geoinformation scientist via the GPS Toolbox () website or the Wolfram Information Center ().

Correction of apparent position shifts caused by GNSS antenna changes

Antenna changes at GNSS reference stations frequently produce discontinuities in the coordinate time series. These apparent position shifts are mainly caused by changes of carrier-phase multipath effects and different errors in the antenna phase center corrections. A monitoring method was developed and successfully tested, which requires additional GNSS observations from a local, temporary reference station. Changes of carrier-phase measurement errors due to the antenna change are determined and stored in L1 and L2 phase maps. These phase maps provide corrections to be applied either to the observation data obtained before the antenna change or to the observation data obtained after the antenna change. The observation corrections are able to remove coordinate discontinuities independent of the selected coordinate estimation algorithm.

Kalman-filter-based GPS clock estimation for near real-time positioning

In this article, an algorithm for clock offset estimation of the GPS satellites is presented. The algorithm is based on a Kalman-filter and processes undifferenced code and carrier-phase measurements of a global tracking network. The clock offset and drift of the satellite clocks are estimated along with tracking station clock offsets, tropospheric zenith path delay and carrier-phase ambiguities. The article provides a brief overview of already existing near-real-time and real-time clock products. The filter algorithm and data processing scheme is presented. Finally, the accuracy of the orbit and clock product is assessed with a precise orbit determination of the MetOp satellite and compared to results gained with other real-time products.

GPS scintillation over the European Arctic during the November 2004 storms

Small-scale irregularities in the background electron density of the ionosphere can cause rapid fluctuations in the amplitude and phase of radio signals passing through it. These rapid fluctuations are known as scintillation and can cause a Global Positioning System (GPS) receiver to lose lock on a signal. This could compromise the integrity of a safety of life system based on GPS, operating in auroral regions. In this paper, the relationship between the loss of lock on GPS signals and ionospheric scintillation in auroral regions is explored. The period from 8 to 14 November 2004 is selected for this study, as it includes both geomagnetically quiet and disturbed conditions. Phase and amplitude scintillation are measured by GPS receivers located at three sites in Northern Scandinavia, and correlated with losses of signal lock in receivers at varying distances from the scintillation receivers. Local multi-path effects are screened out by rejection of low-elevation data from the analysis. The results indicate that losses of lock are more closely related to rapid fluctuations in the phase rather than the amplitude of the received signal. This supports the idea, suggested by Humphreys et al. (2005) (performance of GPS carrier tracking loops during ionospheric scintillations. Proceedings Internationsl Ionospheric Effects Symposium 3–5 May 2005), that a wide loop bandwidth may be preferred for receivers operating at auroral latitudes. Evidence from the Imaging Riometer for Ionospheric Studies (IRIS) appears to suggest that, for this particular storm, precipitation of particles in the D/E regions may be the mechanism that drives the rapid phase fluctuations in the signal.

Analyzing the shear strength parameters of the Chiu-fen-erh-shan landslide: integrating strong-motion and GPS data to determine the best-fit accelerogram

This study adopts the Chiu-fen-erh-shan landslide as a case study for incorporating comprehensive accelerograph and global positioning system (GPS) data to determine the best-fit acceleration data for analyzing a rock avalanche. Previous investigations indicate that the distance from an accelerograph to a landslide site is crucial to determining the best-fit acceleration data to use in conducting a seismic analysis. Unfortunately, the Chiu-fen-erh-shan landslide and its nearest accelerograph station are located in different geological zones. Thus, GPS data were compared to the displacements derived from the accelerograms of nearby monitoring stations to help select the best accelerograph data. This emphasizes that a high density distribution of accelerographs and GPS installations are essential to acquire the best data for the seismic analysis, especially in complex geological zones. After applying the best-fit accelerogram to Newmark’s sliding model and an empirical displacement attenuation formula to back-calculate the shear strength parameters of the sliding surface, a cohesion of 0 kPa and friction angle of the sliding surface of 24.8° were found for this landslide.

Ionospheric climatology derived from gps occultation observations made by the ionospheric occultation experiment

The ionospheric occultation experiment (IOX) is a GPS occultation sensor with an ionospheric mission focus. IOX measurements of GPS L1 and L2 carrier phase during Earth limb views of setting GPS satellites are used together with the Abel transform to determine vertical profiles of electron density from which F-region peak parameters are determined. Data from a four and a half month period beginning in November 2001 are statistically binned and compared with a climatological model. To account for potential errors in interpretation that could arise from violation of the Abel transform assertion of spherical symmetry, the data are compared to both the climatology and to statistics of simulated ionospheric inversions using the climatological model. General characteristics of the climatology are reproduced by the occultation data. However, several significant discrepancies between the model and the data are observed during this near-solar maximum time period. In particular, average mid-latitude daytime densities are shown to be higher than the climatological prediction and the height of F2 layer in the post-sunset equatorial region is underestimated by up to 150 km.

---

Use of C-Band frequencies for satellite navigation: benefits and drawbacks

Although not considered for the first generation of European Galileo satellites, the use of C-Band frequencies for navigation purposes may be taken into account for a future generation of Galileo. For this reason, a frequency band of 20 MHz bandwidth (5,010–5,030 MHz) has been allocated in the course of the World Radio Communications Conference 2000 held in Istanbul, Turkey. The use of C-Band navigation signals offers both advantages and drawbacks. One example is the ionospheric path delay which is inversely proportional to the (squared) carrier frequency and is therefore significantly smaller at C-Band. On the other hand, the use of C-Band frequencies results in increased attenuation effects such as free space loss or rainfall attenuation. It is therefore necessary to provide a detailed analysis of the effects of C-Band frequencies on the navigation process. In order to assess the feasibility of using C-Band frequencies, various aspects of signal propagation and satellite signal tracking at C-Band are examined in the context of this article. In particular, aspects like free space loss, atmospheric effects, foliage attenuation, code and carrier tracking performance, code noise, phase noise and multipath are discussed with respect to their performance at C-Band. In order to allow comparison with the current GPS system, the performance at C-Band is compared to the L-Band performance under similar or identical conditions. The results of this analysis will finally be discussed with respect to their impact on satellite payload and receiver design.

---

Fast direct GPS P-Code acquisition

GPS P-Code has a higher chipping rate, better accuracy, and anti-jamming property than C/A code. Traditionally, GPS P-Code acquisition depends on handover from C/A code. This potentially needs long acquisition time. Moreover, when C/A code is not available, it is no longer possible to acquire GPS P-Code through handover from C/A code. The purpose of this paper is to describe a new overlap average method to facilitate hardware design of fast direct P-Code acquisition. It allows the rapid code phase search to acquire GPS P-Code signals, and also decreases the hardware resource requirement. The small size FFT in the proposed methods is very promising for fast FPGA hardware system design using FFT cores. The simulation results and theoretical analysis are included demonstrating the overall performance of the proposed method.

Critical network infrastructure analysis: interdiction and system flow

Effective management of critical network infrastructure requires the assessment of potential interdiction scenarios. Optimization approaches have been essential for identifying and evaluating such scenarios in networked systems. Although a primary function of any network is the distribution of flow between origins and destinations, the complexity and difficulty of mathematically abstracting interdiction impacts on connectivity or flow has been a challenge for researchers. This paper presents an optimization approach for identifying interdiction bounds with respect to connectivity and/or flow associated with a system of origins and destinations. Application results for telecommunications flow are presented, illustrating the capabilities of this approach.

---

GPS sensitivity analysis applied to non-permanent deformation control networks

This paper illustrates the surveys and the results obtained in an experiment whose goal is to evaluate the Global Positioning System (GPS) sensitivity and accuracy for deformation control on non-permanent network of different extensions. To this aim a high-precision device was properly built to set up known displacements along three orthogonal axes of a GPS antenna. One of the antennas in the considered GPS networks was moved according to centimeter and sub-centimeter displacements; after careful GPS data processing, it was evaluated whether these simulated deformations were correctly a posteriori detected and at which probability level. This experiment was carried out both on a local (baselines ranging between 3 and 30 km) and on a regional (baselines ranging between 300 and 600 km) GPS network. The results show that in the local network it is possible to identify the displacements at a level of 10 mm in height and at a level of 5 mm in horizontal position. The analysis of the regional network showed that it is fundamental to investigate new strategies to model the troposphere; in fact, it is necessary to improve the precision of the height in order to correctly identify displacements lower than 60–80 mm; on the contrary, horizontal displacements can be evidenced at the level of 20 mm. Received: 27 April 1998 / Accepted: 14 December 1998     

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.

Impact of systematic errors on precise long-baseline kinematic GPS positioning

Many kinematic GPS applications rely on high accuracy, which usually requires the ambiguities to be fixed. Normally, a reference station in the rover’s vicinity is needed for successful ambiguity resolution. Alternatively, a network surrounding the rover and allowing one to derive area correction parameters is needed. Unfortunately, both approaches are not feasible in certain situations. This paper is a contribution to precise kinematic positioning over long baselines. Atmospheric refraction becomes critical in the error budget, but progress has been made to use numerical weather models to derive tropospheric corrections, for instance. The spatial correlation of both ionospheric and tropospheric propagation delays is investigated in this paper and special attention is paid on the systematic error behavior of tropospheric refraction. The principles developed are applied to an extended reliability test of the ambiguities. Finally, it is demonstrated in positioning experiments that kinematic positioning retrieval with fixed ambiguities is actually possible for baselines between 150 and 300 km with an accuracy of approximately 2 cm in post-mission processing.

---

A swash–backwash model of the single epidemic wave

While there is a large literature on the form of epidemic waves in the time domain, models of their structure and shape in the spatial domain remain poorly developed. This paper concentrates on the changing spatial distribution of an epidemic wave over time and presents a simple method for identifying the leading and trailing edges of the spatial advance and retreat of such waves. Analysis of edge characteristics is used to (a) disaggregate waves into ‘swash’ and ‘backwash’ stages, (b) measure the phase transitions of areas from susceptible, S, through infective, I, to recovered, R, status (S → I → R) as dimensionless integrals and (c) estimate a spatial version of the basic reproduction number, R ₀. The methods used are illustrated by application to measles waves in Iceland over a 60-year period from 1915 to 1974. Extensions of the methods for use with more complex waves are possible through modifying the threshold values used to define the start and end points of an event.

Hierarchical spatial interaction among the Italian regions: a nonlinear relative dynamics approach

Regional interactions and spillover effects should be considered as important factors in growth analysis of regional economies. Using modified versions of the Dendrinos–Sonis model, this paper analyses the spatial hierarchical system of Italy. The interaction among Italian regions is considered at three different levels of spatial aggregation, the NUTS-1, NUTS-2 and NUTS-3 levels. Compared to more popular spatial econometric approaches, the Dendrinos–Sonis model and its extensions provide greater flexibility in the way interaction between regions is handled but the results strongly depend on the choice of a reference region.

Non-linear motions of Australian geodetic stations induced by non-tidal ocean loading and the passage of tropical cyclones

We investigate daily and sub-daily non-tidal oceanic and atmospheric loading (NTOAL) in the Australian region and put an upper bound on potential site motion examining the effects of tropical cyclone Yasi that crossed the Australian coast in January/February 2011. The dynamic nature of the ocean is important, particularly for northern Australia where the long-term scatter due to daily and sub-daily oceanic changes increases by 20–55 % compared to that estimated using the inverted barometer (IB) assumption. Correcting the daily Global Positioning System (GPS) time series for NTOAL employing either a dynamic ocean model or the IB assumption leads to a reduction of up to 52 % in the weighted scatter of daily coordinate estimates. Differences between the approaches are obscured by seasonal variations in the GPS precision along the northern coast. Two compensating signals during the cyclone require modelling at high spatial and temporal resolution: uplift induced by the atmospheric depression, and subsidence induced by storm surge. The latter dominates ( \(>\) 135 %) the combined net effect that reaches a maximum of 14 mm, and 10 mm near the closest GPS site TOW2. Here, 96 % of the displacement is reached within 15 h due to the rapid transit of cyclones and the quasi-linear nature of the coastline. Consequently, estimating sub-daily NTOAL is necessary to properly account for such a signal that can be 3.5 times larger than its daily-averaged value. We were unable to detect the deformation signal in 2-hourly GPS processing and show that seasonal noise in the Austral summer dominates and precludes GPS detection of the cyclone-related subsidence.     

Nonlinear regional economic dynamics: continuous-time specification,estimation and stability analysis

This paper presents an innovative approach to the study of regional economic dynamics within a nonlinear continuous-time econometric framework—a generalized specification of the Lotka–Volterra system of equations. This specification, which accounts for interdependent behavior of three industrial sectors and spillover effects of activities in neighboring regions, is employed in an analysis of five Italian regions between 1980 and 2003. For these regions, we report estimation results, characterize the varying systems dynamics, analyze the models’ local and global stability properties, and determine via sensitivity analyses which structural features appear to exert the greatest influence on these properties.