Integer aperture bootstrapping: a new GNSS ambiguity estimator with controllable fail-rate

In this contribution, we introduce a new bootstrap-based method for Global Navigation Satellite System (GNSS) carrier-phase ambiguity resolution. Integer bootstrapping is known to be one of the simplest methods for integer ambiguity estimation with close-to-optimal performance. Its outcome is easy to compute due to the absence of an integer search, and its performance is close to optimal if the decorrelating Z-transformation of the LAMBDA method is used. Moreover, the bootstrapped estimator is presently the only integer estimator for which an exact and easy-to-compute expression of its fail-rate can be given. A possible disadvantage is, however, that the user has only a limited control over the fail-rate. Once the underlying mathematical model is given, the user has no freedom left in changing the value of the fail-rate. Here, we present an ambiguity estimator for which the user is given additional freedom. For this purpose, use is made of the class of integer aperture estimators as introduced in Teunissen (2003). This class is larger than the class of integer estimators. Integer aperture estimators are of a hybrid nature and can have integer outcomes as well as non-integer outcomes. The new estimator is referred to as integer aperture bootstrapping. This new estimator has all the advantages known from integer bootstrapping with the additional advantage that its fail-rate can be controlled by the user. This is made possible by giving the user the freedom over the aperture of the pull-in region. We also give an exact and easy-to-compute expression for its controllable fail-rate.     

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.     

K 近邻非参数回归概率预报技术及其应用

针对参数回归技术制作概率预报存在拟合好、但预报结果不稳定的现象, 提出了用K近邻非参数回归技术制作概率预报的新途径。K 近邻非参数回归技术包括历史样本数据库、近邻子集生成和优化以及预报量估计4 个主要部分。利用该技术进行了单要素概率预报(主要包括云量和降水)和多维联合概率预报(降水、总云量、风速和气温)试验, 并对试验结果进行了检验。实例研究结果表明:该文所给出的计算方案预报稳定性好, 准确率较高,具有良好的业务应用价值。     

改进的切线性模式对一个边界层模式变分资料同化的改善

大气、海洋中边界层物理过程的强非线性会对切线性模式近似（尤其是当积分时间长时）有较大的影响，从而给相应的4维变分资料同化问题的求解造成困难。本构造了Mellor-Yamada湍流封闭模式（level 2.5)的一个改进的切线性模式，相比通常的切线性模式和简化的切线性模式可以提高对非线性扰动的逼近。利用这个改进的切线性模式的伴随模式分别进行了1到7天的变分资料同化试验，得到了满意的结果。而用通常的伴随模式和简化的伴随模式都无法得到令人满意的结果。     

Seismic response of self‐centring hysteretic SDOF systems

The seismic response of single‐degree‐of‐freedom (SDOF) systems incorporating flag‐shaped hysteretic structural behaviour, with self‐centring capability, is investigated numerically. For a SDOF system with a given initial period and strength level, the flag‐shaped hysteretic behaviour is fully defined by a post‐yielding stiffness parameter and an energy‐dissipation parameter. A comprehensive parametric study was conducted to determine the influence of these parameters on SDOF structural response, in terms of displacement ductility, absolute acceleration and absorbed energy. This parametric study was conducted using an ensemble of 20 historical earthquake records corresponding to ordinary ground motions having a probability of exceedence of 10% in 50 years, in California. The responses of the flag‐shaped hysteretic SDOF systems are compared against the responses of similar bilinear elasto‐plastic hysteretic SDOF systems. In this study the elasto‐plastic hysteretic SDOF systems are assigned parameters representative of steel moment resisting frames (MRFs) with post‐Northridge welded beam‐to‐column connections. In turn, the flag‐shaped hysteretic SDOF systems are representative of steel MRFs with newly proposed post‐tensioned energy‐dissipating connections. Building structures with initial periods ranging from 0.1 to 2.0s and having various strength levels are considered. It is shown that a flag‐shaped hysteretic SDOF system of equal or lesser strength can always be found to match or better the response of an elasto‐plastic hysteretic SDOF system in terms of displacement ductility and without incurring any residual drift from the seismic event. Copyright © 2002 John Wiley & Sons, Ltd.     

Towards a direct collapse–load method of design for concrete frames subjected to severe ground motions

Most current methods of design for concrete structures under earthquake loads rely on highly idealized ‘equivalent’ static representations of the seismic loads and linear‐elastic methods of structural analysis. With the continuing development of non‐linear methods of dynamic analysis for the overload behaviour and collapse of complete concrete structures, a more direct and more accurate design procedure becomes possible which considers conditions at system collapse. This paper describes an evaluation procedure that uses non‐linear dynamic collapse–load analysis together with global safety coefficients. A back‐calibration procedure for evaluating the global safety coefficients is also described. The aim of this paper is to open up discussion of alternative methods of design with improved accuracy which are necessary to move towards a direct collapse–load method of design. Copyright © 2002 John Wiley & Sons, Ltd.     

海底天然气渗漏形成水合物量的线性动力学模型:以墨西哥湾GC185区块Bush高地为例

世界许多地区如墨西哥湾，存在着大量天然气渗漏并形成水合物的现象，出于科学，生态，气候和安全的需要。人类有必要弄清存在于海底以水合成形式存在的天然气的比例。因此，构建了水合物形成量的线性动力学模型，以分析墨西哥湾GC185块区Bush高地海底渗漏天然气的地质过程和行为。Bush高地渗漏的天然气来源于附近的Jolliet气藏，结果表明，渗漏天然气中约有9％的海底形成了水合物。     

Method of Data Reduction and Uncertainty Estimation for Platinum-Group Element Data Using Inductively Coupled Plasma-Mass Spectrometry

A data reduction method is described for determining platinum-group element (PGE) abundances by inductively coupled plasma-mass spectrometry (ICP-MS) using external calibration or the method of standard addition. Gravimetric measurement of volumes, the analysis of reference materials and the use of procedural blanks were all used to minimise systematic errors. Internal standards were used to correct for instrument drift. A linear least squares regression model was used to calculate concentrations from drift-corrected counts per second (cps). Furthermore, mathematical manipulations also contribute to the uncertainty estimates of a procedure. Typical uncertainty estimate calculations for ICP-MS data manipulations involve: (1) Carrying standard deviations from the raw cps through the data reduction or (2) calculating a standard deviation from multiple final concentration calculations. It is demonstrated that method 2 may underestimate the uncertainty estimate of the calculated data. Methods 1 and 2 do not typically include an uncertainty estimate component from a regression model. As such models contribute to the uncertainty estimates affecting the calculated data, an uncertainty estimate component from the regression must be included in any final error calculations. Confidence intervals are used to account for uncertainty estimates from the regression model. These confidence intervals are simpler to calculate than uncertainty estimates from method 1, for example. The data reduction and uncertainty estimation method described here addresses problems of reporting PGE data from an article in the literature and addresses both precision and accuracy. The method can be applied to any analytical technique where drift corrections or regression models are used.     

非线性地球物理反演方法:回顾与展望

自20世纪90年代以来，非线性地球物理反演方法已走向成熟，这些方法包括线性化迭代法，仿真淬火法，遗传算法及联合反演方法等，线性化迭代反演基于开放物理系统状态发生相变的原理，要进一步改善模型参数化的技术和迭代过程的自组织；仿真淬火法与遗传算法基于自然过程的指数率或生物演化的优生率，可以相互结合以提高解估计的分辨率与置信度；联合反演要结合岩石物理性质的统计规律，才能取得兼容地质与综合方法的应用效果，地球动力学中的反问题不仅涉及偏微分方程系数项的求解，而且涉及初始条件或初始边界的求解，对地球动力学作用过程研究特别重要。