Least squares spectral analysis and its application to superconducting gravimeter data analysis

Detection of a periodic signal hidden in noise is the goal of Superconducting Gravimeter (SG) data analysis. Due to spikes, gaps, datum shrifts (offsets) and other disturbances, the traditional FFT method shows inherent limitations. Instead, the least squares spectral analysis (LSSA) has showed itself more suitable than Fourier analysis of gappy, unequally spaced and unequally weighted data series in a variety of applications in geodesy and geophysics. This paper reviews the principle of LSSA and gives a possible strategy for the analysis of time series obtained from the Canadian Superconducting Gravimeter Installation (CGSI), with gaps, offsets, unequal sampling decimation of the data and unequally weighted data points.     

最小二乘谱及其在超导重力观测数据分析中的应用

运用投影理论和基于向量空间中最小二乘平差原理,介绍了最小二乘谱分析方法,给出了最小二乘谱的计算公式。利用加拿大超导重力观测数据实例,讨论了最小二乘谱分析的策略和步骤。     

Least-squares cross-wavelet analysis and its applications in geophysical time series

The least-squares wavelet analysis, an alternative to the classical wavelet analysis, was introduced in order to analyze unequally spaced and non-stationary time series exhibiting components with variable amplitude and frequency over time. There are a few methods such as cross-wavelet transform and wavelet coherence that can analyze two time series together. However, these methods cannot generally be used to analyze unequally spaced and non-stationary time series with associated covariance matrices that may have trends and/or datum shifts. A new method of analyzing two time series together, namely the least-squares cross-wavelet analysis, is developed and applied to study the disturbances in the gravitational gradients observed by GOCE satellite that arise from plasma flow in the ionosphere represented by Poynting flux. The proposed method also shows its outstanding performance on the Westford–Wettzell very long baseline interferometry baseline length and temperature series.     

On the accuracy assessment of celestial reference frame realizations

The scatter of the celestial pole offset (CPO) time-series obtained from very long baseline interferometry (VLBI) observations is used as a measure of the accuracy of celestial reference frame (CRF) realizations. Several scatter indices (SI), including some proposed for the first time, are investigated. The first SI is based on residual analysis of CPO series with respect to a free core nutation (FCN) model. The second group of SIs includes Allan deviation and its extensions, which allow the treatment of unequally weighted and multidimensional observations. Application of these criteria to several radio source catalogues (RSCs) showed their ability to perform a preliminary assessment of the quality of each RSC. The 2D Allan deviation estimate seems to be the most sensitive measure. The proposed extensions of Allan deviation, weighted and multidimensional, can also be used for the statistical analysis of other time-series.     

多星定轨条件下北斗卫星钟差的周期性变化

基于地面监测网的多星精密定轨可以同时解算出北斗卫星轨道和卫星钟差。由于轨道和钟差的耦合影响,卫星钟差时序难免会出现周期性波动。此外,受限于目前并不完善的北斗全球监测网络分布、系统导航文件缺失以及定轨后处理软件的设置问题,3类卫星的钟差均存在大量数据间断问题。本文利用适用于间断数据的谱分析方法,对多星定轨条件下的北斗卫星钟差数据进行了周期项提取,并利用周期项改进后的钟差预报模型评估了24h以内的预报精度。基于近一年的数据分析表明,北斗GEO卫星钟差3个主周期依次为12、24和8h,IGSO卫星钟差的3个主周期依次为24、12和8h,而MEO卫星钟差的3个主周期依次为12.91、6.44和24h。与改进前相比,周期项改进后的钟差预报模型将北斗卫星钟差在24h以内的预报精度提高了20%~40%。     

GPS变形监测系统中消除噪声的一种有效方法

GPS变形监测系统的观测资料可看成为与时间有关的数据序列，应用小波分析理论，研究了时间序列观测数据的误差消除问题。结果表明，借助于小波分解与重构，可有效地从强噪声干扰的观测数据序列中提取变形特征。该方法解决了传统处理技术对非平稳、非等时间间隔观测数据序列滤波的局限性问题。     

iGPS: IDL tool package for GPS position time series analysis

A new tool package written in Interactive Data Language (IDL) was developed for processing and analyzing daily continuous GPS position time series. This software package can read continuous GPS position time series in various formats, detect outliers, remove abnormal observation spans, locate epochs of offsets or post-seismic relaxation events, and estimate their amplitudes. It also provides functionalities for epoch statistics, site selection, periodic noise analysis, spatial filtering (by regional stacking), etc. This tool is referred to as iGPS.     

剔除VAV软件中扰动数据

新的潮汐分析软件VAV不但可以处理带有间断的数据,还可以采用迭代的方法不断地检测和去除扰动数据,提高潮汐参数的求解精度。介绍VAV软件的迭代原理,并利用VAV软件对武汉国际重力潮汐基准台和布鲁塞尔台站超导重力仪记录的重力潮汐数据进行分析处理,潮汐参数的标准偏差可减小30％以上。     

Troposphere gradients from the ECMWF in VLBI analysis

Modeling path delays in the neutral atmosphere for the analysis of Very Long Baseline Interferometry (VLBI) observations has been improved significantly in recent years by the use of elevation-dependent mapping functions based on data from numerical weather models. In this paper, we present a fast way of extracting both, hydrostatic and wet, linear horizontal gradients for the troposphere from data of the European Centre for Medium-range Weather Forecasts (ECMWF) model, as it is realized at the Vienna University of Technology on a routine basis for all stations of the International GNSS (Global Navigation Satellite Systems) Service (IGS) and International VLBI Service for Geodesy and Astrometry (IVS) stations. This approach only uses information about the refractivity gradients at the site vertical, but no information from the line-of-sight. VLBI analysis of the CONT02 and CONT05 campaigns, as well as all IVS-R1 and IVS-R4 sessions in the first half of 2006, shows that fixing these a priori gradients improves the repeatability for 74% (40 out of 54) of the VLBI baseline lengths compared to fixing zero or constant a priori gradients, and improves the repeatability for the majority of baselines compared to estimating 24-h offsets for the gradients. Only if 6-h offsets are estimated, the baseline length repeatabilities significantly improve, no matter which a priori gradients are used.     

An analysis of vertical deflections derived from high-degree spherical harmonic models

The theoretical differences between the Helmert deflection of the vertical and that computed from a truncated spherical harmonic series of the gravity field, aside from the limited spectral content in the latter, include the curvature of the normal plumb line, the permanent tidal effect, and datum origin and orientation offsets. A numerical comparison between deflections derived from spherical harmonic model EGM96 and astronomic deflections in the conterminous United States (CONUS) shows that correcting these systematic effects reduces the mean differences in some areas. Overall, the mean difference in CONUS is reduced from −0.219 arcsec to −0.058 arcsec for the south–north deflection, and from +0.016 arcsec to +0.004 arcsec for the west–east deflection. Further analysis of the root-mean-square differences indicates that the high-degree spectrum of the EGM96 model has significantly less power than implied by the deflection data. Received: 9 December 1997 / Accepted: 21 August 1998     

超导重力仪观测数据分析

超导重力仪作为目前研究时变重力场的主要仪器具有很高的灵敏度和稳定性、很低的噪音水平和漂移以及很宽的动态频率响应范围。探测时间序列中隐含的周期信号是超导重力仪观测数据分析的主要目标之一,通过快速傅里叶变换(FFT)和最小二乘谱分析方法分析超导重力数据中隐含信号的检测,并对两种方法结果进行了分析比较。并最终用小波分析的方法进行了检校。     

利用GNSS数据分析大港验潮站地壳沉降

青岛大港验潮站的地壳沉降关系到该站平均海平面的绝对变化,因而也就关系到我国高程基准面的变化。本文利用青岛GNSS基准站约10年的观测数据对该站的地壳沉降变化进行分析。首先将青岛GNSS基准站纳入由50个国际IGS站和43个国内陆态网络基准站组成的全球网中,进行单日松弛解和单日约束解解算,获得该站坐标时间序列。然后对该站垂向坐标时间序列进行分析,利用粗差探测、偏差探测、趋势项分析、频谱分析等方法对粗差、偏差、趋势项和周期项进行探测、分析,并通过时间序列模型估计获得时间序列中的周期项振幅和偏差估值。分析表明青岛GNSS基准站垂直方向近一段时间未发现存在显著性的地壳沉降变化,但受到比较明显的周年和半周年周期变化影响。结合青岛大港验潮站验潮数据分析结果得出结论:青岛大港验潮站平均海平面的绝对上升速率是1.62mm/a。     

EPN coordinate time series monitoring for reference frame maintenance

The EUREF permanent network (EPN) is a network of continuously operating global navigation satellite system (GNSS) stations, primarily installed for reference frame maintenance. In order to ensure the long term reliability of the EPN products, a thorough station performance monitoring has been initiated and carried out in addition to the routine GNSS data management, processing and analysis. This paper addresses the main factors influencing the quality of the coordinate time series in a permanent GNSS network. Relevant examples, based on the EPN experience are given, the analysis strategy is introduced, the estimated coordinate offsets are published and the importance of this analysis for site velocity estimation is demonstrated. The results are derived from the analysis of the EPN weekly combined solutions covering the period from 1996 to 2003. Our target is the identification, interpretation and elimination of offsets and outliers present in the EPN coordinate time series in order to estimate reliable coordinates and velocities and consequently maintain a high quality kinematic reference network.

遥感定量反演地表参数对于气候变化的响应研究

遥感定量反演地表参数时间序列产品已被广泛应用于植被动态变化、全球气候变化、防灾减灾及环境保护等领域。由于卫星观测往往受到大气条件（如云、气溶胶、水汽等）以及传感器自身稳定性的影响等,许多由卫星观测反演得到的陆表产品,如归一化差值植被指数（NDVI）、叶面积指数（LAI）、地表温度（LST）、微波极化亮温（PDBT）等存在严重的时空不连续问题。为了获取时间序列上连续、空间上完整的地表参数遥感产品以满足长时序的陆面过程分析与建模的需求,目前已发展多种遥感时间序列重建模型。本文介绍了基于傅里叶变换的时间序列谐波分析（HANTS）方法,能够识别并去除受到云和大气影响的像元（噪声）,对原始时序数据进行时间插值来重建连续时间序列的数据,并针对其面向多种不同时空尺度的遥感反演地表参数以及在非洲、南美洲、欧洲、中国及印度等全球不同地区的应用研究进行了综述,包括植被动态变化对于气候变化及流域水循环过程的响应、干旱监测、基于土壤含水量饱和度时间序列分析的洪涝灾害易发区监测、遥感估算地表蒸散发时间尺度扩展等方面的研究,充分阐释了遥感时间序列产品在地气相互作用的各类研究领域的应用。     

顾及误差频谱特性的CHZ重力仪航空应用研究

给出了航空重力测量误差频域分析的方法,利用功率谱密度从频域分析了航空标量重力测量系统恢复重力场的能力及影响因素。介绍了CHZ重力仪的主要特点,并利用实测空中重力异常数据及机载GPS动态加速度数据,结合航空重力测量的频谱范围,分析了CHZ重力仪在不同阻尼系数下的动态性能。计算结果表明,采用合适的阻尼系数,CHZ重力仪能够被用于固定翼飞机的航空重力测量。     

基于 Markov 理论的加权非等距GM（1，1）预测优化模型

背景值的构造方法是影响加权非等距GM(1,1)预测模型的精度和适应性的关键因素。文中通过等分函数法构造新的背景值对传统的加权非等距GM(1,1)模型进行优化,优化后的模型使其同时适应于高增长指数序列和低增长指数序列,提高传统模型的预测精度和适应性能力。但是优化后的模型依然易受建模数据随机扰动影响。马尔科夫(Markov)模型具有削弱建模数据的随机扰动性的优势。基于此,将优化的加权非等距GM(1,1)模型和Markov理论有机结合,构建优化的加权非等距Markov-GM(1,1)预测模型。最后,结合秀山湖二期工程的变形实测数据,运用新陈代谢的计算模式进行预测验证。结果表明:优化的加权非等距Markov-GM(1,1)预测模型的拟合和预测精度都优于传统的加权非等距GM(1,1)预测模型,新的预测模型的适用性更强,具有实际的参考价值。     

GPS position time-series analysis based on asymptotic normality of M-estimation

The efficacy of robust M-estimators is a well-known issue when dealing with observational blunders. When the number of observations is considerably large—long time series for instance—one can take advantage of the asymptotic normality of the M-estimation and compute reasonable estimates for the unknown parameters of interest. A few leading M-estimators have been employed to identify the most likely functional model for GPS coordinate time series. This includes the simultaneous detection of periodic patterns and offsets in the GPS time series. Estimates of white noise, flicker noise, and random walk noise components are also achieved using the robust M-estimators of (co)variance components, developed in the framework of the least-squares variance component estimation (LS-VCE) theory. The method allows one to compute confidence interval for the (co)variance components in asymptotic sense. Simulated time series using white noise plus flicker noise show that the estimates of random walk noise fluctuate more than those of flicker noise for different M-estimators. This is because random walk noise is not an appropriate noise structure for the series. The same phenomenon is observed using the results of real GPS time series, which implies that the combination of white plus flicker noise is well described for GPS time series. Some of the estimated noise components of LS-VCE differ significantly from those of other M- estimators. This reveals that there are a large number of outliers in the series. This conclusion is also affirmed by performing the statistical tests, which detect (large) parts of the outliers but can also leave parts to be undetected.     

Detecting discontinuities in GNSS coordinate time series with STARS: case study,the Bologna and Medicina GPS sites

Global navigation satellite systems (GNSS) data are a fundamental source of information for achieving a better understanding of geophysical and climate-related phenomena. However, discontinuities in the coordinate time series might be a severe limiting factor for the reliable estimate of long-term trends. A methodological approach has been adapted from Rodionov (Geophys Res Lett 31:L09204, 2004; Geophys Res Lett 31:L12707, 2006) and from Rodionov and Overland (J Marine Sci 62:328–332, 2005) to identify both the epoch of occurrence and the magnitude of jumps corrupting GNSS data sets without any a priori information on these quantities. The procedure is based on the Sequential t test Analysis of Regime Shifts (STARS) (Rodionov in Geophys Res Lett 31:L09204, 2004). The method has been tested against a synthetic data set characterized by typical features exhibited by real GNSS time series, such as linear trend, seasonal cycle, jumps, missing epochs and a combination of white and flicker noise. The results show that the offsets identified by the algorithm are split into 48 % of true-positive, 28 % of false-positive and 24 % of false-negative events. The procedure has then been applied to GPS coordinate time series of stations located in the southeastern Po Plain, in Italy. The series span more than 15 years and are affected by offsets of different nature. The methodology proves to be effective, as confirmed by the comparison between the corrected GPS time series and those obtained by other observation techniques.     

Polar motion modeling,analysis, and prediction with time dependent harmonic coefficients

A time dependent amplitude model was proposed for the analysis and prediction of polar motion time series. The formulation was implemented to analyze part of the new combined solution, EOP (IERS) C 04, daily polar motion time series of 14 years length using a statistical model with first order autoregressive disturbances. A new solution approach, where the serial correlations of the disturbances are eliminated by sequentially differencing the measurements, was used to estimate the model parameters using weighted least squares. The new model parsimoniously represents the 14-year time series with 0.5 mas rms fit, close to the reported 0.1 mas observed pole position precisions for the x and y components. The model can also predict 6 months into the future with less than 4 mas rms prediction error for both polar motion components, and down to sub mas for one-step ahead prediction as validated using a set of daily time series data that are not used in the estimation. This study is dedicated to the memory of Prof. Urho Uotila (1923–2006) whose teaching of “Adjustment Computations” over the years influenced so much, so many of us who had the privilege of being his students.     

顾及钟差变化率的GPS卫星钟差预报法

提出了一种基于钟差变化率拟合建模的卫星钟差预报方法。以附加周期项的线性或二次多项式作为基础模型对钟差变化率序列进行拟合,最优估计卫星钟差的趋势项系数,然后直接使用精密定轨得到的相应时刻的卫星钟差计算预报初始时刻的基准项系数,来建立卫星钟差的预报模型。以IGS发布的快速星历(IGR)的卫星钟差为试验数据,对GPS星座中各种型号的所有卫星钟差进行预报。结果表明:本文方法3、6、12与24h的预报精度分别可达0.43、0.58、0.90与1.47ns,相比于传统的基于钟差拟合的预报方法,精度分别提高69.3%、61.8%、50.5%与37.2%;与IGS发布的超快速星历(IGU)的预报钟差相比,钟差精度分别提高15.7%、23.7%、27.4%与34.4%。