Efficient GOCE satellite gravity field recovery based on least-squares using QR decomposition

We develop and apply an efficient strategy for Earth gravity field recovery from satellite gravity gradiometry data. Our approach is based upon the Paige-Saunders iterative least-squares method using QR decomposition (LSQR). We modify the original algorithm for space-geodetic applications: firstly, we investigate how convergence can be accelerated by means of both subspace and block-diagonal preconditioning. The efficiency of the latter dominates if the design matrix exhibits block-dominant structure. Secondly, we address Tikhonov-Phillips regularization in general. Thirdly, we demonstrate an effective implementation of the algorithm in a high-performance computing environment. In this context, an important issue is to avoid the twofold computation of the design matrix in each iteration. The computational platform is a 64-processor shared-memory supercomputer. The runtime results prove the successful parallelization of the LSQR solver. The numerical examples are chosen in view of the forthcoming satellite mission GOCE (Gravity field and steady-state Ocean Circulation Explorer). The closed-loop scenario covers 1 month of simulated data with 5 s sampling. We focus exclusively on the analysis of radial components of satellite accelerations and gravity gradients. Our extensions to the basic algorithm enable the method to be competitive with well-established inversion strategies in satellite geodesy, such as conjugate gradient methods or the brute-force approach. In its current development stage, the LSQR method appears ready to deal with real-data applications.     

Ocean tide loading (OTL) displacements from global and local grids: comparisons to GPS estimates over the shelf of Brittany, France

In this paper we examine OTL displacements detected by GPS stations of a dedicated campaign and validate ocean tide models. Our area of study is the continental shelf of Brittany and Cotentin in France. Brittany is one of the few places in the world where tides provoke loading displacements of ∼10–12 cm vertically and a few cm horizontally. Ocean tide models suffer from important discrepancies in this region. Seven global and regional ocean tide models were tested: FES2004 corrected for K2, TPXO.7.0, TPXO.6.2, GOT00.2, CSR4.0, NAO.99b and the most recent regional grids of the North East Atlantic (NEA2004). These gridded amplitudes and phases of ocean tides were convolved in order to get the predicted OTL displacements using two different algorithms. Data over a period of 3.5 months of 8 GPS campaign stations located on the north coast of Brittany are used, in order to evaluate the geographical distribution of the OTL effect. We have modified and implemented new algorithms in our GPS software, GINS 7.1. GPS OTL constituents are estimated based on 1-day batch solutions. We compare the observed GPS OTL constituents of M₂, S₂, N₂ and K₁ waves with the selected ocean tide models on global and regional grids. Large phase-lag and amplitude discrepancies over 20° and 1.5 cm in the vertical direction in the semi-diurnal band of M₂ between predictions and GPS/models are detected in the Bay of Mont St-Michel. From a least squares spectral analysis of the GPS time-series, significant harmonic peaks in the integer multiples of the orbital periods of the GPS satellites are observed, indicating the existence of multipath effects in the GPS OTL constituents. The GPS OTL observations agree best with FES2004, NEA2004, GOT00.2 and CSR4.0 tide models.     

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.     

Accuracy assessment of lidar-derived digital elevation models

Despite the relatively high cost of airborne lidar-derived digital elevation models (DEMs), such products are usually presented without a satisfactory associated estimate of accuracy. For the most part, DEM accuracy estimates are typically provided by comparing lidar heights against a finite sample of check point coordinates from an independent source of higher accuracy, supposing a normal distribution of the derived height differences or errors. This paper proposes a new methodology to assess the vertical accuracy of lidar DEMs using confidence intervals constructed from a finite sample of errors computed at check points. A non-parametric approach has been tested where no particular error distribution is assumed, making the proposed methodology especially applicable to non-normal error distributions of the type usually found in DEMs derived from lidar. The performance of the proposed model was experimentally validated using Monte Carlo simulation on 18 vertical error data-sets. Fifteen of these data-sets were computed from original lidar data provided by the International Society for Photogrammetry and Remote Sensing Working Group III/3, using their respective filtered reference data as ground truth. The three remaining data-sets were provided by the Natural Environment Research Council's Airborne Research and Survey Facility lidar system, together with check points acquired using high precision kinematic GPS. The results proved promising, the proposed models reproducing the statistical behaviour of vertical errors of lidar using a favourable number of check points, even in the cases of data-sets with non-normally distributed residuals. This research can therefore be considered as a potentially important step towards improving the quality control of lidar-derived DEMs.     

An Efficient Reordering Prediction-Based Lossless Compression Algorithm for Hyperspectral Images

In this letter, we propose an efficient lossless compression algorithm for hyperspectral images; it is based on an adaptive spectral band reordering algorithm and an adaptive backward previous closest neighbor (PCN) prediction with error feedback. The adaptive spectral band reordering algorithm has some strong points. It can adaptively determine the range of spectral bands needed to be reordered, and it can efficiently find the optimum branches. Hyperspectral images have a large number of spectral bands, which express the same land cover structure and have high correlation. The adaptive backward PCN prediction with error feedback can sufficiently make use of this correlation. Experiments show that implementing both the reordering of the spectral bands before prediction and the prediction with error feedback improve compression performance     

Comparison of Deep Underwater Measurements and Radar Observations of Rainfall

Deep-water acoustical measurements of rainfall are compared to high-resolution ground radar observations for the first time. The measurements of underwater ambient sound were made from a subsurface mooring near Methoni, Greece, in 2004. The acoustical measurements were at 60-, 200-, 1000-, and 2000-m depths. Simultaneous ground-based polarimetric -band radar observations were made over the acoustic mooring. Comparisons show acoustic detection of rain events and storm structure that are in agreement with the radar observations. Results from a comparison between the underwater sound pressure level at different depths and the observed radar reflectivities are presented.     

High-Resolution Forecast Models of Water Vapor Over Mountains: Comparison With MERIS and Meteosat Data

Propagation delay due to variable tropospheric water vapor (WV) is one of the most intractable problems for radar interferometry, particularly over mountains. The WV field can be simulated by an atmospheric model, and the difference between the two fields is used to correct the radar interferogram. Here, we report our use of the U.K. Met Office Unified Model in a nested mode to produce high-resolution forecast fields for the 3-km-high Mount Etna volcano. The simulated precipitable-water field is validated against that retrieved from the Medium-Resolution Imaging Spectrometer (MERIS) radiometer on the Envisat satellite, which has a resolution of 300 m. Two case studies, one from winter (November 24, 2004) and one from summer (June 25, 2005), show that the mismatch between the model and the MERIS fields ( rms = 1.1 and 1.6 mm, respectively) is small. One of the main potential sources of error in the models is the timing of the WV field simulation. We show that long-wavelength upper tropospheric troughs of low WV could be identified in both the model output and Meteosat WV imagery for the November 24, 2004 case and used to choose the best time of model output.     

Best prediction in linear models with mixed integer/real unknowns: theory and application

In this contribution, we extend the existing theory of minimum mean squared error prediction (best prediction). This extention is motivated by the desire to be able to deal with models in which the parameter vectors have real-valued and/or integer-valued entries. New classes of predictors are introduced, based on the principle of equivariance. Equivariant prediction is developed for the real-parameter case, the integer-parameter case, and for the mixed integer/real case. The best predictors within these classes are identified, and they are shown to have a better performance than best linear (unbiased) prediction. This holds true for the mean squared error performance, as well as for the error variance performance. We show that, in the context of linear model prediction, best predictors and best estimators come in pairs. We take advantage of this property by also identifying the corresponding best estimators. All of the best equivariant estimators are shown to have a better precision than the best linear unbiased estimator. Although no restrictions are placed on the probability distributions of the random vectors, the Gaussian case is derived separately. The best predictors are also compared with least-squares predictors, in particular with the integer-based least-squares predictor introduced in Teunissen (J Geodesy, in press, 2006).     

Non-destructive estimation of potato leaf chlorophyll from canopy hyperspectral reflectance using the inverted PROSAIL model

Optimizing nitrogen (N) fertilization in crop production by in-season measurements of crop N status may improve fertilizer N use efficiency. Hyperspectral measurements may be used to assess crop N status by estimating leaf chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content. Trials were conducted with two potato cultivars under different N fertility rates (0–300 kg N ha⁻¹). Canopy reflectance, leaf area index (LAI) and leaf chlorophyll and N contents were measured. The PROSAIL model was able to predict leaf chlorophyll content with reasonable accuracy later in the growing season. The low estimation accuracy earlier in the growing season could be due to model sensitivity to non-homogenous canopy architecture and soil background interference before full canopy closure. Canopy chlorophyll content (leaf chlorophyll content × LAI) was predicted less accurately than leaf chlrophyll content due to the low estimation accuracy of LAI for values higher than 4.5.     

GPS sidereal filtering: coordinate- and carrier-phase-level strategies

Multipath error is considered one of the major errors affecting GPS observations. One can benefit from the repetition of satellite geometry approximately every sidereal day, and apply filtering to help minimize this error. For GPS data at 1 s interval processed using a double-difference strategy, using the day-to-day coordinate or carrier-phase residual autocorrelation determined with a 10-h window leads to the steadiest estimates of the error-repeat lag, although a window as short as 2 h can produce an acceptable value with > 97% of the optimal lag’s correlation. We conclude that although the lag may vary with time, such variation is marginal and there is little advantage in using a satellite-specific or other time-varying lag in double-difference processing. We filter the GPS data either by stacking a number of days of processed coordinate residuals using the optimum “sidereal” lag (23 h 55 m 54 s), and removing these stacked residuals from the day in question (coordinate space), or by a similar method using double-difference carrier-phase residuals (observational space). Either method results in more consistent and homogeneous set of coordinates throughout the dataset compared with unfiltered processing. Coordinate stacking reduces geometry-related repeating errors (mainly multipath) better than carrier-phase residual stacking, although the latter takes less processing time to achieve final filtered coordinates. Thus, the optimal stacking method will depend on whether coordinate precision or computational time is the over-riding criterion.