Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997–2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the “truth” signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.     

A new method for indentifying the model error of adjustment system

Some theory problems affecting parameter estimation are discussed in this paper. Influence and transformation between errors of stochastic and functional models is pointed out as well. For choosing the best adjustment model, a formula, which is different from the literatures existing methods, for estimating and identifying the model error, is proposed. On the basis of the proposed formula, an effective approach of selecting the best model of adjustment system is given. Project supported by the Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, Wuhan University (No. 905276031-04-10).     

Investigating the propagation mechanism of unmodelled systematic errors on coordinate time series estimated using least squares

The propagation of unmodelled systematic errors into coordinate time series computed using least squares is investigated, to improve the understanding of unexplained signals and apparent noise in geodetic (especially GPS) coordinate time series. Such coordinate time series are invariably based on a functional model linearised using only zero and first-order terms of a (Taylor) series expansion about the approximate coordinates of the unknown point. The effect of such truncation errors is investigated through the derivation of a generalised systematic error model for the simple case of range observations from a single known reference point to a point which is assumed to be at rest by the least squares model but is in fact in motion. The systematic error function for a one pseudo-satellite two-dimensional case, designed to be as simple but as analogous to GPS positioning as possible, is quantified. It is shown that the combination of a moving reference point and unmodelled periodic displacement at the unknown point of interest, due to ocean tide loading, for example, results in an output coordinate time series containing many periodic terms when only zero and first-order expansion terms are used in the linearisation of the functional model. The amplitude, phase and period of these terms is dependent on the input amplitude, the locations of the unknown point and reference point, and the period of the reference point's motion. The dominant output signals that arise due to truncation errors match those found in coordinate time series obtained from both simulated data and real three-dimensional GPS data.     

显式自由表面模型时间分裂格式的一个积分方案

讨论了有关显式自由表面模型的时间分裂格式中的分裂误差。在传统的积分方案中，由于采用蛙跳格式，斜压动量方程具有三层的时间离散结构，但是，正压动量方程在相应的时间层上只有二层的离散结构。因此，正压方程与斜压方程的时间离散结构是不相容的。分析结果表明，这是导致分裂误差的原因。本提出一种能保持时间离散结构相容的正压和斜压方程的时间积分方案。该积分方案的分裂误差很小因此特别适合于需要长期积分的气候研究。该方案在北太平洋的一个自由表面模型中得到了应用。     

EH-4系统观测资料的非远区场校正研究

EH－4系统实测的卡尼亚视电阻率在非远区场时应作必要的校正。基于导出的EH－4系统的卡尼亚电阻率ρxy和ρyx的表达式，采用多项式分段逼近的方法，建立EH－4系统测深数据的全频域的视电阻率算法。根据建立的算法对均匀大地、二层大地和三层大地模型的EH－4系统测深数据进行非远区场校正的结果表明：均匀大地时，全频域的视电阻率与真电阻率吻合很好，D型、K型和Q型断面的视电阻率曲线与MT法的接近，G型、H型和A型的视电阻率曲线得到很大改善。     

Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Static deformation of two welded monoclinic elastic half-spaces due to a long inclined strike-slip fault

Static deformation of two monoclinic elastic half-spaces in welded contact due to a long inclined strike-slip fault situated in one of the half-spaces is studied analytically and numerically. Closed-form algebraic expressions for the displacement at any point of the medium are obtained. The variation of the displacement at the interface with the horizontal distance from the fault is studied. The effect of anisotropy on the displacement field is examined. It is found that while the anisotropy of the source half-space has a significant effect on the displacement at the interface, the anisotropy of the other half-space has only a marginal effect.     

一种计算固体潮观测数据改正值的新方法

从标定前一天和标定当天的地形变前兆数据变化趋势入手，在分析当天数据结构的基础上采用最小二乘法原理，提出了一种计算数采改正值的新方法。经典型实例计算表明，该方法精度较高，效果明显。并对误差的来源作了进一步的分析。     

ARTIFICIAL NEURAL NETWORK CORRECTION OF ROUGH ERRORS OF OBSERVATIONS*

In the context of tower measured radiation datasets.following the correction principle meeting a diagnostic equation in data quality control and in terms of a technique for model construction on data and ANN (artificial neural network) retrieval for BP correction of radiation measurements with rough errors available,a BP model is presented.Evidence suggests that the developed model works well and is superior to a convenient multivariate linear regression model,indicating its wide applications.     

Precise revision of the garnet-muscovite geothermometer

The garnet-muscovite geothermometer was refined through empirical calibration by using natural rocks metamorphosed under the physical conditions of 238—1306 MPa and 490—700℃. Input temperatures and pressures were determined through simultaneously applying the garnet-biotite geothermometer and the garnet-biotite-plagioclase-quartz barometer, assuming that all FeO in muscovite and garnet be ferrous. Garnet was treated as the asymmetric quaternary solid solution, and muscovite as the symmetric binary solid solution. Input muscovite compositions include Fe atoms between 0.03—0.19 and Mg atoms between 0.04—0.16 on the basis of 11 oxygen atoms, and input garnet compositions include spessartine fractions between 0.01—0.289, grossular fractions between 0.028—0.273, and the Fe/Mg ratio between 3.387-18.986. The resulting garnet-muscovite geothermometer reproduces temperatures within (50℃ compared with the garnet-biotite thermometer. Total random error of ±37℃ of the new thermometer may stem from the pressure uncertainty of ±200 MPa, and uncertainties of ±5% of Fe and Mg components in muscovite, and ±5% of Fe, Mg, Mn and Ca components in garnet, altogether. When there exist 10%, 20%, 30%, 40% and 50% Fe3+ in muscovite, respectively, the computed garnet-muscovite temperatures will be 1—6℃, 2—12℃, 3—16℃, 5—24℃ and 7—29℃, respectively, lower than those obtained when assuming that all FeO be ferrous. The new garnet-muscovite geothermometer can efficiently reflect temperature change of typical prograde sequences and contact aureole rocks, and may be applied to low- to high-grade and low- to high-pressure metamorphic rocks.