SIMD correlator library for GNSS software receivers

Software receivers have had a discernable impact on the GNSS research community. Often such receivers are implemented in a compiled programming language, such as C or C++. A software receiver must emulate the digital signal processing (DSP) algorithms executed on dedicated hardware in a traditional receiver. The DSP algorithms, most notably correlation, have a high computational cost; this burden precludes many software receivers from running in real time. However, the computational cost can be lessened by utilizing single instruction multiple data (SIMD) operations found on modern ×86 processors. The following demonstrates how C/C++ compatible code can be written to directly utilize the SIMD instructions. First, an analysis is carried out to demonstrate why real time operation is not possible when using traditional C/C++ code is carried out. Secondly a tutorial outlines how to write and insert ×86 assembly, with SIMD operations, into C/C++ code. Performance gains achieved via SIMD operations are then demonstrated, and pseudo code outlines how SIMD operations can be employed to perform correlation. Finally, a C/C++ compatible SIMD enabled arithmetic library is added to the GPS Toolbox for use in software receivers.     

Modeling the length of day and extrapolating the rotation of the Earth

The stochastic behavior of the length of day (LOD) process is analyzed and is modeled within statistical accuracy on a time-scale ranging from weeks to millennia by a three-component model comprising a global Brownian motion process, decadal fluctuations, and a 50-day Madden–Julian oscillation. While the model is intended to be phenomenological, some possible physical models underlying the three components are speculated upon. The model is applied to estimate long-range extrapolation errors. For example, it predicts a standard error of 1 h in the clock-time correction ΔT for extrapolation by 1,500 years from 500 to 2000 BC.     

Reviewers of the Journal of Geodesy for Volume 79

Acknowledgement for Referees

Reviewers of the Journal of Geodesy for Volume 79     

Efficient satellite orbit modelling using pseudo-stochastic parameters

If the force field acting on an artificial Earth satellite is not known a priori with sufficient accuracy to represent its observations on their accuracy level, one may introduce so-called pseudo-stochastic parameters into an orbit determination process, e.g. instantaneous velocity changes at user-defined epochs or piecewise constant accelerations in user-defined adjacent time subintervals or piecewise linear and continuous accelerations in adjacent time subintervals. The procedures, based on standard least-squares, associated with such parameterizations are well established, but they become inefficient (slow) if the number of pseudo-stochastic parameters becomes large. We develop two efficient methods to solve the orbit determination problem in the presence of pseudo-stochastic parameters. The results of the methods are identical to those obtained with conventional least-squares algorithms. The first efficient algorithm also provides the full variance–covariance matrix; the second, even more efficient algorithm, only parts of it.     

Effect of signal contamination in matched-filter detection of the signal on a cluttered background

To derive a matched filter for detecting a weak target signal in a hyperspectral image, an estimate of the band-to-band covariance of the target-free background scene is required. We investigate the effects of including some of the target signal in the background scene. Although the covariance is contaminated by the presence of a target signal (there is increased variance in the direction of the target signature), we find that the matched filter is not necessarily affected. In fact, if the variation in plume strength is strictly uncorrelated with the variation in background spectra, the matched filter and its signal-to-clutter ratio (SCR) performance will not be impaired. While there is little a priori reason to expect significant correlation between the plume and the background, there usually is some residual correlation, and this correlation leads to a suppressing effect that limits the SCR obtainable even for strong plumes. These effects are described and quantified analytically, and the crucial role of this correlation is illustrated with some numerical examples using simulated plumes superimposed on real hyperspectral imagery. In one example, we observe an order-of-magnitude loss in SCR for a matched filter based on the contaminated covariance.     

ARGO upper salinity measurements: perspectives for L-band radiometers calibration and retrieved sea surface salinity validation

With the view of preparing the strategy for the calibration/validation of future L-band satellite radiometers, we examine the salinity variability recorded by Array for Real-Time Geostrophic Oceanography (ARGO) floats in the upper 10-m layer of the surface ocean. Using one year of ARGO measurements,we show that the surface salinity variability at ten days and 200-km scales is above /spl plusmn/0.1 psu for 30% of the drifters and that this variability is larger than 0.2 psu in tropical regions affected by strong river discharges and by precipitations, and in frontal areas characterized by strong mesoscale activity. Vertical gradient observed between 5-10-m depth is much lower than the horizontal variability but leads to systematic biases in the tropics. The South Pacific Ocean appears to be the less variable both vertically and horizontally.     

Plate kinematics of Nubia–Somalia using a combined DORIS and GPS solution

We have used up to 12 years of data to assess DORIS performance for geodynamics applications. We first examine the noise characteristics of the DORIS time-series of weekly station coordinates to derive realistic estimates of velocity uncertainties. We find that a combination of white and flicker noise best explains the DORIS time-series noise characteristics. Second, weekly solutions produced by the Institut Géographique National/Jet Propulsion Laboratory (IGN/JPL) DORIS Analysis Centre are combined to derive a global velocity field. This solution is combined with two independent GPS solutions, including 11 sites on Nubia and 5 on the Somalia plate. The combination indicates that DORIS horizontal velocities have an average accuracy of 3 mm/year, with best-determined sites having velocity accuracy better than 1 mm/year (one-sigma levels). Using our combined velocity field, we derive an updated plate kinematics model with a focus on the Nubia–Somalia area. Including DORIS data improves the precision of the angular velocity vector for Nubia by 15%. Our proposed model provides robust bounds on the maximum opening rates along the East African Rift (4.7–6.7 mm/year). It indicates opening rates 15 and 7% slower than values predicted by NUVEL-1A for the southern Atlantic Ocean and Indian Ocean, respectively. These differences are likely to arise from the fact that NUVEL-1A considered Africa as a single non-deforming plate, while here we use a more refined approach.