Zonal-Harmonics Perturbations

A recently developed iterative technique for solving the perturbed Kepler problem is explained, in a step-by-step detail, using the simple example of oblateness perturbations. The results are then extended to deal with higher-degree zonal harmonics.     

各向异性组合滤波法反演陆地水储量变化

地球时变重力场模型反演陆地水储量变化已为全球气候变化研究作出巨大贡献,考虑到时变重力场模型球谐系数中存在相关性,其高阶次项具有较大的误差,需采用最优的滤波方法进行空间平滑。本文提出一种新的各向异性组合滤波方法,其基本思想是将改进的高斯滤波法与均方根(root mean square,RMS)滤波法组合,即对球谐系数的低阶次采用改进的高斯滤波法,而高阶次采用RMS滤波法。首先分析了最新的GRACE RL05系列时变重力场模型系数误差特性,基于全球水储量变化反演结果,分析比较了高斯滤波、改进的高斯滤波、RMS滤波和DDK滤波与本文提出的组合滤波法的有效性及精度,并利用模型结果进行了验证,计算结果表明,组合滤波法的中误差最小。研究结果表明,本文提出的组合法相比于先前的滤波方法,可有效地过滤高阶次的噪声,消除南北条带误差,同时减少信号泄漏,提高信噪比,保留更多有效的地球物理信号,进而提高反演精度。     

京津冀地区地磁场球冠谐分析

2002年在京津冀地区进行了45个测点的地磁三分量测量,对测量资料进行通化处理,通化时间为2002年5月5日16～18时（世界时）. 通化后的观测均方差分别优于1.5nT（地磁场总强度F）,0.5′(磁偏角D和磁倾角I). 将国际参考地磁场（IGRF2000）作为地磁正常场,建立了京津冀地区地磁异常场的球冠谐模型（BTHASCH）. 球冠极的空间位置坐标为39.5°N和117.0°E,球冠半角为4°. 在模型计算过程中,球冠谐函数的截断阶数分别取为1～10. 经综合比较, 最终采用的截断阶数为5. 建立了京津冀地区参考地磁场的球冠谐模型（BTHGRF）.根据模型,绘制了京津冀地区地磁异常场图(ΔX、ΔY、ΔZ、ΔF、ΔD、ΔI)和京津冀地区地磁图(X、Y、Z、F、D、I).     

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.     

Modelling of mantle postglacial relaxation in axisymmetric geometry with a composite rheology and a glacial load interpolated by adjusted spherical harmonics analysis

Although studies on glacial isostatic adjustment usually assume a purely linear rheology, we have previously shown that mantle relaxation after the melting of Laurentide ice sheet is better described by a composite rheology including a non-linear term. This modelling is, however, based on axially symmetric geometry and glacial forcing derived from ICE-3G and suffers from a certain amount of arbitrariness in the definition of the ice load. In this work we apply adjusted spherical harmonics analysis to interpolate the ice thicknesses of ICE-3G and ICE-1 glaciological models. This filters out the non-axisymmetric components of the ice load by considering only the zonal terms in the spherical harmonics expansion. The resulting load function is used in finite-element simulation of postglacial rebound to compare composite versus purely linear rheology. Our results confirm that composite rheology can explain relative sea level (RSL) data in North America significantly better than a purely linear rheology. The performance of composite rheology suggests that in future investigations, it may be better to use this more physically realistic creep law for modelling mantle deformation induced by glacial forcing.     

Nonlinear Shallow Water Theories for Coastal Waves

Ocean waves entering the near-shore zone undergo nonlinear and dispersive processes. This paper reviews nonlinear models, focusing on the so-called Serre equations. Techniques to overcome their limitations with respect to the phase speed are presented. Nonlinear behaviours are compared with theoretical results concerning the properties of Stokes waves. In addition, the models are tested against experiments concerning periodic wave transformation over a bar topography and of the shoaling of solitary waves on a beach.