Analysis of mapping coverage obtained from spacecraft

An exact analysis of the coverage obtained by spacecraft using cross-track scanning and nadir-centered conical imaging, under imposed viewing obliqueness and resolution requirements, is presented. In addition to exact expressions for the area acquired and the area acquisition rate, envelope theory is introduced to obtain the boundary of the imaged area. These expressions are relatively compact, allowing rapid machine computation. The effects of the sun phase angle, and of imaging system limitations are also examined. The Galileo mission encounter with Callisto is used as a numerical example, from which certain general conclusions are drawn regarding optimal imaging trajectories.     

Observations with the Visual and Infrared Mapping Spectrometer (VIMS) during Cassini's flyby of Jupiter

The Cassini Visual and Infrared Mapping Spectrometer (VIMS) is an imaging spectrometer covering the wavelength range 0.3-5.2 μm in 352 spectral channels, with a nominal instantaneous field of view of 0.5 mrad. The Cassini flyby of Jupiter represented a unique opportunity to accomplish two important goals: scientific observations of the jovian system and functional tests of the VIMS instrument under conditions similar to those expected to obtain during Cassini's 4-year tour of the saturnian system. Results acquired over a complete range of visual to near-infrared wavelengths from 0.3 to 5.2 μm are presented. First detections include methane fluorescence on Jupiter, a surprisingly high opposition surge on Europa, the first visual-near-IR spectra of Himalia and Jupiter's optically-thin ring system, and the first near-infrared observations of the rings over an extensive range of phase angles (0-120°). Similarities in the center-to-limb profiles of H⁺₃ and CH₄ emissions indicate that the H⁺₃ ionospheric density is solar-controlled outside of the auroral regions. The existence of jovian NH₃ absorption at 0.93 μm is confirmed. Himalia has a slightly reddish spectrum, an apparent absorption near 3 μm, and a geometric albedo of 0.06±0.01 at 2.2 μm (assuming an 85-km radius). If the 3-μm feature in Himalia's spectrum is eventually confirmed, it would be suggestive of the presence of water in some form, either free, bound, or incorporated in layer-lattice silicates. Finally, a mean ring-particle radius of 10 μm is found to be consistent with Mie-scattering models fit to VIMS near-infrared observations acquired over 0-120° phase angle.     

新世纪大洋钻探"三角鼎立"中的欧洲联合体

欧洲为了与美国、日本在IODP中处于同等重要的地位，成立了由ODP成员参加的欧洲大洋钻探科学指导委员会（ESCOD）和JEODI组织，并开展了IODP之前的一系列准备活动。在其中的APLACON会议上，论证了IODP中第三类特定钻探平台的必要性，讨论了该类平台在某些研究范围内的重要作用，包括研究地球历史中的极端气候（北极海区）、重建白垩纪－第三纪的水文地理学、气候的快速变化、沉积盆地形成与过程以及固体地球过程等。     

Methodology and use of tensor invariants for satellite gravity gradiometry

Although its use is widespread in several other scientific disciplines, the theory of tensor invariants is only marginally adopted in gravity field modeling. We aim to close this gap by developing and applying the invariants approach for geopotential recovery. Gravitational tensor invariants are deduced from products of second-order derivatives of the gravitational potential. The benefit of the method presented arises from its independence of the gradiometer instrument’s orientation in space. Thus, we refrain from the classical methods for satellite gravity gradiometry analysis, i.e., in terms of individual gravity gradients, in favor of the alternative invariants approach. The invariants approach requires a tailored processing strategy. Firstly, the non-linear functionals with regard to the potential series expansion in spherical harmonics necessitates the linearization and iterative solution of the resulting least-squares problem. From the computational point of view, efficient linearization by means of perturbation theory has been adopted. It only requires the computation of reference gravity gradients. Secondly, the deduced pseudo-observations are composed of all the gravitational tensor elements, all of which require a comparable level of accuracy. Additionally, implementation of the invariants method for large data sets is a challenging task. We show the fundamentals of tensor invariants theory adapted to satellite gradiometry. With regard to the GOCE (Gravity field and steady-state Ocean Circulation Explorer) satellite gradiometry mission, we demonstrate that the iterative parameter estimation process converges within only two iterations. Additionally, for the GOCE configuration, we show the invariants approach to be insensitive to the synthesis of unobserved gravity gradients.     

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.     

Study on Recovering the Earth＇s Potential Field Based on GOCE Gradiometry

Given the second radial derivative Vrr（P） ｜δs of the Earth＇s gravitational potential V（P） on the surface δS corresponding to the satellite altitude, by using the fictitious compress recovery method, a fictitious regular harmonic field rrVrr（P）^＊ and a fictitious second radial gradient field V：（P） in the domain outside an inner sphere Ki can be determined, which coincides with the real field V（P） in the domain outside the Earth. Vrr^＊（P）could be further expressed as a uniformly convergent expansion series in the domain outside the inner sphere, because rrV（P）^＊ could be expressed as a uniformly convergent spherical harmonic expansion series due to its regularity and harmony in that domain. In another aspect, the fictitious field V^＊（P） defined in the domain outside the inner sphere, which coincides with the real field V（P） in the domain outside the Earth, could be also expressed as a spherical harmonic expansion series. Then, the harmonic coefficients contained in the series expressing V^＊（P） can be determined, and consequently the real field V（P） is recovered. Preliminary simulation calculations show that the second radial gradient field Vrr（P） could be recovered based only on the second radial derivative V（P）｜δs given on the satellite boundary. Concerning the final recovery of the potential field V（P） based only on the boundary value Vrr （P）｜δs, the simulation tests are still in process.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.     

A parametric study on the impact of satellite attitude errors on GOCE gravity field recovery

In the recent design of the Gravity field and steady-state Ocean Circulation Explorer (GOCE) satellite mission, the gravity gradients are defined in the gradiometer reference frame (GRF), which deviates from the actual flight direction (local orbit reference frame, LORF) by up to 3–4°. The main objective of this paper is to investigate the effect of uncertainties in the knowledge of the gradiometer orientation due to attitude reconstitution errors on the gravity field solution. In the framework of several numerical simulations, which are based on a realistic mission configuration, different scenarios are investigated, to provide the accuracy requirements of the orientation information. It turns out that orientation errors have to be seriously considered, because they may represent a significant error component of the gravity field solution. While in a realistic mission scenario (colored gradiometer noise) the gravity field solutions are quite insensitive to small orientation biases, random noise applied to the attitude information can have a considerable impact on the accuracy of the resolved gravity field models.     

TRMM卫星微波成像仪分级产品及其反演降水算法

文章叙述了获取定量降水信息的意义，简要介绍了对热带测雨卫星TRMM(Tropical Rainfall Measurement Mission)的仪器、美国国家宇航局提供的微波成像仪TMI(TRMM Microwave Imager)分级产品。对比了物理方法和经验方法反演降水的特点，并对一些经验方法以及倾斜对流系统对反演降水的影响、动态聚类分析、神经网络反演方法的研究成果进行了介绍。     

On the effect of a low viscosity asthenosphere on the temporal change of the geoid—A challenge for future gravity missions

New satellite technology to measure changes in the Earth’s gravity field gives new possibilities to detect layers of low viscosity inside the Earth. We used density models for the Earth mantle based on slab history as well as on tomography and fitted the viscosity by comparison of predicted gravity to the new CHAMP gravity model. We first confirm that the fit to the observed geoid is insensitive to the presence of a low viscosity anomaly in the upper mantle as long as the layer is thin ( 200 km) and the viscosity reduction is less than two orders of magnitude. Then we investigated the temporal change in geoid by comparing two stages of slablet sinking based on subduction history or by advection of tomography derived densities and compared the spectra of the geoid change for cases with and without a low viscosity layer, but about equal fit to the observed geoid. The presence of a low viscosity layer causes relaxation at smaller wavelength and thus leads to a spectrum with relatively stronger power in higher modes and a peak around degrees 5 and 6. Comparing the spectra to the expected degree resolution for GRACE data for a 5 years mission duration shows a weak possibility to detect changes in the Earth’s gravity field due to large scale mantle circulation, provided that other causes of geoid changes can be taken into account with sufficient accuracy. A discrimination between the two viscosity cases, however, demands a new generation of gravity field observing satellites.