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.     

经验模式分解的尺度阈值滤波在变形分析中的应用

本文将信号处理领域的经验模式分解算法应用于变形信息的提取中,通过引入阈值函数,建立了基于经验模式分解的尺度阈值滤波模型,采用优化模型确定了经验模式分解的次数。分别通过模拟试验和实测数据与小波阈值法和多项式拟合法进行了比对,分析表明:在低噪声情况下,三种方法都有一定的滤波效果;在高噪声情况下,经验模式分解的尺度阈值滤波法具有与小波阈值去噪法相等的精度,而且瞬时强噪声识别能力更好,优于多项式拟合法。     

On the Importance of Geolithological Features for the Estimate of the Site Response: The Case of Catania Metropolitan Area (Italy)

The evaluation of seismic site response in the urban area of Catania was tackled by selecting test areas having peculiar lithological and structural features, potentially favourable to large local amplifications of ground motion. The two selected areas are located in the historical downtown and in the northern part of Catania where the presence of a fault is evident. Site response was evaluated using spectral ratio technique taking the horizontal- to-vertical component ratio of ambient noise. Inferences from microtremor measurements are compared with results from synthetic accelerograms and response spectra computed at all drillings available for this area. Such method is particularly suitable in urban areas where the nature of the outcropping geological units is masked by city growth and anthropic intervention on the surface geology. The microtremor H/V spectral ratios evaluated at soft sites located within the downtown profile tend to be smaller than that usually reported in the literature for such soils. A tendency for amplifications to peaks near 2 Hz is observed only in some sites located on recent alluvial deposits. Evidences for amplifications of site effects (frequency range 4–8 Hz) were observed in the sampling sites located on the fault, with a rapid decrease of spectral amplitude just a few tenths of metres away from the discontinuity. Numerical simulations evidenced the importance of geolithological features at depth levels even greater than 20–30 m. Besides this, the results strongly confirm the importance of the subsurface geological conditions, in the estimate of seismic hazard at urban scale.     

Assessing the Vibrational Frequencies of the Cathedral of Cologne (Germany) by Means of Ambient Seismic Noise Analysis

Ambient seismic noise measurements were conducted inside the Cathedral of Cologne (Germany) for assessing its frequencies of vibration and for checking whether these occur in the range where soil amplification is expected. If this is the case, damages may increase in case of an earthquake due to an increased structural response of the building. Analysis of the ratio between the horizontal and vertical components of the spectra recorded at stations located inside the building as well as the ratio between the corresponding components of the spectra recorded simultaneously inside the building and at a reference station placed in the basement of the cathedral indicated several modes of vibration. Facilitated by these results an assessment of the seismic vulnerability was attempted for a 2D ground motion scenario using the finite element method.     

Site Effects Study in the Serravalle Village (NE Italy) by Means of Spectral Ratios Analyses Relatively to Multiple Reference Sites

The area of Serravalle, sited in the northern part of the town of Vittorio Veneto (TV), NE Italy, has been the target of a seismic microzonation campaign. 10 seismic stations have been deployed for a 7 months period to record in continuous mode. Three stations were installed on bedrock outcrops and seven on sedimentary sites with variable cover thickness. Spectral analyses have been performed on the collected data-set using the Generalized Inversion Technique (GIT, e.g. Andrews, 1986). In particular, spectral ratios have been calculated for each station relatively to the average of the three reference, bedrock sites. The spectral ratios provide quantitative estimates of the seismic motion amplifications which occur in each of the monitored sites. Two sites show high values of amplification, 5 times larger than signal amplitude at the reference sites, in correspondence of well discernible peak frequencies of 5 Hz. Results for the other stations show smaller amounts of site amplification spreading over a broad range of frequencies. Sites where the highest amplifications were recorded all lie on the left bank of the Meschio River and in areas farther away from its outlet into the plain correlating with the presence of thick layers of Quaternary deposits.     

Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. We show that this potential-field estimation problem has deep-seated connections to Slepian's spatiospectral localization problem which seeks bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures such as bias, variance and mean squared error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of our method make the use of spherical Slepian functions in earth and planetary geodesy practical.     

Spectral Simulation and Conditioning to Local Continuity Trends in Geologic Modeling

Spectral simulation has gained application in building geologic models due to the advantage of better honoring the spatial continuity of petrophysical properties, such as reservoir porosity and shale volume. Distinct from sequential simulation methods, spectral simulation is a global algorithm in the sense that a global density spectrum is calculated once and the inverse Fourier transform is performed on the Fourier coefficient also only once to generate a simulation realization. The generated realizations honor the spatial continuity structure globally over the whole field instead of only within a search neighborhood, as with sequential simulation algorithms. However, the disadvantage of global spectral simulation is that it traditionally cannot account for the local information such as the local continuity trends, which are often observed in reservoirs and hence are important to be accounted for in geologic models. This disadvantage has limited wider application of spectral simulation in building geologic models. In this paper, we present ways of conditioning geologic models to the relevant local information. To account for the local continuity trends, we first scale different frequency components of the original model with local-amplitude spectrum ratios that are specific to the local trend. The sum of these scaled frequency components renders a new model that displays the desired local continuity trend. The implementation details of this new method are discussed and examples are provided to illustrate the algorithm.     

基于高光谱卫星遥感数据的UPDM分析方法

提出了一种用于处理多高光谱卫星数据的UPDM分析方法。研究结果证明．该方法应用干Landsat／TM(ETM^-)、Terra／MODIS和ADEOS—II／GLI等高光谱卫星传感器时．光谱重构均方根误差小于0．029适用于研究高光谱卫星遥感数据。     

单频GPS整周模糊度动态快速求解的研究

提出了一种动态快速求解整周模糊度的方法,即先对系数阵进行QR分解,然后通过矩阵变换使模糊度参数和位置参数分离,从而降低矩阵的维数,满足实时动态求解的要求,最后应用LAMBDA方法搜索模糊度。为验证该算法,用单频GPS接收机进行了静态定位和动态定位两种试验。结果表明,静态定位误差小于1cm,动态定位误差小于3cm。由此可见,该方法能够为动态用户快速而精确地实施GPS单频动态定位。