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.     

Distribution of vertical crustal movement rates in the Tohoku district, Japan, predicted by least-squares collocation

Least-squares prediction using an empirically deduced local covariance function was performed to investigate the temporal change in the rates of vertical crustal movements in the Tohoku district, Japan. Levelling data covering an area of approximately 450 × 275 km² observed between 1966 and 1995 were used and the results shown in the form of contour maps. Firstly we derived a covariance function of the rates of vertical crustal movement with a Gaussian form function. We used this function to estimate the spatial distribution of the rates of vertical crustal movements. By the present method, a steady tilt of the Tohoku district to the east, toward the Japan Trench and an areal uplift in the southwestern part were well reproduced. Moreover, a significant temporal change in vertical movement rates is clearly seen. Received: 15 July 1996      

Least-squares self-coherency analysis of superconducting gravimeter records in search for the Slichter triplet

We develop a new approach for the spectral analysis of the superconducting gravimeter data to search for the spheroidal oscillation ₁S₁ of the Earth solid inner core. The new method, which we call least-squares (LS) self-coherency analysis, is based on the product of the least-squares spectra of segments of the time series under consideration. The statistical foundation of this method is presented in the new least-squares product spectrum theorem that establishes rigorously confidence levels for detecting significant peaks. We apply this approach along with a number of other innovative ideas to a 6-year long gravity series collected at the Canadian Superconducting Gravimeter Installation (CSGI) in Cantley, Canada, by splitting it into 72 statistically independent monthly records. Each monthly record is analysed spectrally and all monthly LS spectra are multiplied to construct the self-coherency spectrum of the 6-year gravity series. The self-coherency spectrum is then used to detect significant peaks in the band 3-7 h at various significant levels with the aim to identify a triplet of periods associated with the rotational/ellipsoidal splitting of ₁S₁ (Slichter triplet). From all the Slichter periods predicted by various researchers so far, Smylie's triplet appears to be the most supported one, albeit very weakly, both, before and after the atmospheric pressure effect is removed from the series. Using the viscous splitting law [Smylie, D.E., 1992. The inner core translational triplet and the density near Earth's center. Science 255, 1678-1682] as guide, we can also see one interesting and statistically significant triplet with periods A = {4.261 h, 4.516 h, 4.872 h}, which changes slightly to A′ = {4.269 h, 4.516 h, 4.889 h} after the atmospheric pressure correction is applied to the gravity series.     

Mapping and initial dilution estimation of an ocean outfall plume using an autonomous underwater vehicle

A monitoring mission to study the shape and estimate initial dilution of the S. Jacinto outfall plume using an autonomous underwater vehicle (AUV) was performed on July 30, 2002. In order to reduce the uncertainty about plume location and to concentrate the vehicle mission only in the hydrodynamic mixing zone, outputs of a near-field prediction model, based on effective real-time in situ measurements of current speed and direction and density stratification, were opportunistically used to specify in real time the mission transects. The surface characteristics of the outfall plume were found to be influenced strongly by the relatively weak stratification and low current velocities. Dilution was estimated using a temperature–salinity (TS-) diagram with initial mixing lines between wastewater and ambient waters. Effluent dilutions were at least 30:1 in this study. In order to efficiently map the plume dispersion we applied the least-squares collocation method technique. Our results demonstrate that AUVs can provide high-quality measurements of physical properties of effluent plumes in a quite effective manner and valuable considerations about the initial mixing processes under real oceanic conditions can be further investigated.     

New analytical and numerical approaches for geopotential modeling

 The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential coefficients _n,m . In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal harmonic series is developed for the geopotential and then its coefficients _n,m ^e are transformed to the unknown _n,m . The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From two sets of _n,m (n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068) by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to apply the LS technique for evaluating geopotential coefficients. Received: 12 December 2000 / Accepted: 21 June 2001     

大连双井测汞震例分析

根据大连双井测汞震例，用最小二乘法原理，就汞测量预测地震进行了定理化研究，给出适宜本地区的经验预测公式，取得初步研究成果。     

A study on the combination of satellite, airborne, and terrestrial gravity data

Satellite gravity missions, such as CHAMP, GRACE and GOCE, and airborne gravity campaigns in areas without ground gravity will enhance the present knowledge of the Earths gravity field. Combining the new gravity information with the existing marine and ground gravity anomalies is a major task for which the mathematical tools have to be developed. In one way or another they will be based on the spectral information available for gravity data and noise. The integration of the additional gravity information from satellite and airborne campaigns with existing data has not been studied in sufficient detail and a number of open questions remain. A strategy for the combination of satellite, airborne and ground measurements is presented. It is based on ideas independently introduced by Sjöberg and Wenzel in the early 1980s and has been modified by using a quasi-deterministic approach for the determination of the weighting functions. In addition, the original approach of Sjöberg and Wenzel is extended to more than two measurement types, combining the Meissl scheme with the least-squares spectral combination. Satellite (or geopotential) harmonics, ground gravity anomalies and airborne gravity disturbances are used as measurement types, but other combinations are possible. Different error characteristics and measurement-type combinations and their impact on the final solution are studied. Using simulated data, the results show a geoid accuracy in the centimeter range for a local test area.     

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).     

基于平面波加速的VTI介质最小二乘逆时偏移

地震各向异性集中表现为速度各向异性,势必影响地震波运动学特征.传统声波逆时偏移(RTM)和最小二乘逆时偏移(LSRTM)没有考虑介质各向异性特征,导致反射波不能正确归位、同相轴出现扭曲及寻优速度慢或不收敛等,VTI介质逆时偏移(VTI-RTM)矫正了声波成像的不足,但仍存在低频干扰严重、中深部成像不佳、振幅保持差等缺陷.为此,本文首先实现了VTI介质最小二乘逆时偏移(VTI-LSRTM)方法,为了节省I/O及内存需求并提高效率,进一步引入平面波编码技术,提出了一种基于平面波加速的VTI介质最小二乘逆时偏移(VTI-PLSRTM)策略.在此基础上开展了简单模型及复杂Marmousi模型成像试验,并与标准逆时偏移剖面对比表明:本方法能够校正各向异性造成的相位畸变,且在迭代中自动压制串扰及低频噪声、补偿中深部能量,是一种兼具质量与效率的保幅成像策略;对速度误差的敏感性测试说明该方法需要相对正确的偏移速度及Thomsen参数模型.     

�����沨���������й˼������Ҫ�Ե�����

????????????????й??????????????????????LSQR??WLSQR??MLS?????????????????????????????????????????????沨?????????в???????????????????????????????????????????????