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.     

黑龙江气候变暖的时空变化特征

进入80年代以来，黑龙江气温明显变暖，近120年来年平均气温上升1.4℃，春季次之，气温突变在1990年前后，从年、各季平均气温变化趋势与全国北半球及全球气温变化进行比较，年及冬季我国增暖中心在黑龙江，春季增暖中心在内蒙古北部及黑龙江省；夏秋增暖不明显，但仍有上升趋势，近50年来，50～70年代黑龙江省及全国各大区域普遍增暖不明显，80年代开始黑龙江省增暖明显，三北地区次之，90年代黑龙江增暖更加明显，而全国及其他区域上升已达到80年代北方增温幅度，这里从黑龙江变暖的事实，突变时间，与全球、北半球、全国区域的气温及雪盖海温变化关系来揭示黑龙江气候变化的时空特点。     

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

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

Electrical aerosol spectrometer of Tartu University

The electrical aerosol spectrometer (EAS) of the parallel measuring principle at Tartu University is an efficient instrument for rapid measurement of the unstable size spectrum of aerosol particles. The measuring range from 10 nm to 10 μm is achieved by simultaneously using a pair of differential mobility analyzers with two different particle chargers. The particle spectrum is calculated and measurement errors are estimated in real time by using a least-squares method. Experimental calibration ensures reliability of measurement. The instrument is well suited for continuous monitoring of atmospheric aerosol.     

Multi‐objective optimal design of FLC driven hybrid mass damper for seismically excited structures

Structural vibration control using active or passive control strategy is a viable technology for enhancing structural functionality and safety against natural hazards such as strong earthquakes and high wind gusts. Both the active and passive control systems have their limitations. The passive control system has limited capability to control the structural response whereas the active control system depends on external power. The power requirement for active control of civil engineering structures is usually quite high. Thus, a hybrid control system is a viable solution to alleviate some of the limitations. In this paper a multi‐objective optimal design of a hybrid control system for seismically excited building structures has been proposed. A tuned mass damper (TMD) and an active mass driver (AMD) have been used as the passive and active control components of the hybrid control system, respectively. A fuzzy logic controller (FLC) has been used to drive the AMD as the FLC has inherent robustness and ability to handle the non‐linearities and uncertainties. The genetic algorithm has been used for the optimization of the control system. Peak acceleration and displacement responses non‐dimensionalized with respect to the uncontrolled peak acceleration and displacement responses, respectively, have been used as the two objectives of the multi‐objective optimization problem. The proposed design approach for an optimum hybrid mass damper (HMD) system, driven by FLC has been demonstrated with the help of a numerical example. It is shown that the optimum values of the design parameters of the hybrid control system can be determined without specifying the modes to be controlled. The proposed FLC driven HMD has been found to be very effective for vibration control of seismically excited buildings in comparison with the available results for the same example structure but with a different optimal absorber. Copyright © 2002 John Wiley & Sons, Ltd.     

Study on effects of damping in laminated rubber bearings on seismic responses for a ⅛ scale isolated test structure

The effects of damping in various laminated rubber bearings (LRB) on the seismic response of a ?‐scale isolated test structure are investigated by shaking table tests and seismic response analyses. A series of shaking table tests of the structure were performed for a fixed base design and for a base isolation design. Two different types of LRB were used: natural rubber bearings (NRB) and lead rubber bearings (LLRB). Three different designs for the LLRB were tested; each design had a different diameter of lead plug, and thus, different damping values. Artificial time histories of peak ground acceleration 0.4g were used in both the tests and the analyses. In both shaking table tests and analyses, as expected, the acceleration responses of the seismically isolated test structure were considerably reduced. However, the shear displacement at the isolators was increased. To reduce the shear displacement in the isolators, the diameter of the lead plug in the LLRB had to be enlarged to increase isolator damping by more than 24%. This caused the isolator stiffness to increase, and resulted in amplifying the floor acceleration response spectra of the isolated test structure in the higher frequency ranges with a monotonic reduction of isolator shear displacement. Copyright © 2002 John Wiley & Sons, Ltd.     

马岭油田北三区污水回注的实验研究

根据马岭油田北三区油田水性质，结合污水与地层水混合的静态实验结果分析，发现污水回注对地层可能引起不同程度的伤害。同时，模拟地层条件做岩心流动实验。将采出的污水回注到地层与地层水混合后，亦会使地层渗透率下降，且随着污水含量的增加，渗透率伤害程度增大。     

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.     

震源时间函数与震源破裂过程

简要介绍了现代数字地震学的两个重要概念。一个是震源时间函数，另一个是地震震源的时空破裂过程。首先，从地震断层位移表示定量出发，介绍震源时间函数和震源破裂过程的基本概念。然后，分别介绍两种从远场地震记录中提取震源时间函数和获取有限断层面上时空破裂过程图像的反演方法。     

Most likely locations of large earthquakes in the Kanto and Tokai areas, Japan, based on the local recurrence times

We examined the hypothesis that minima in local recurrence time, T_L, or equivalently maxima in local probability, P_L, may map asperities in the Kanto and Tokai areas of Japan, where the earthquake catalog of the National Research Institute for Earth Science and Disaster Prevention (NIED) is complete at the M=1.5 (M1.5) level. We mapped T_L (P_L) based on the a- and b-values of the nearest earthquakes within 20 km of every node of a grid spaced 0.01° for M7 target events. Only earthquakes within the top 33 km were used. The b-values increase strongly with depth, in several areas. Therefore, some of the T_L (P_L) anomalies are not revealed if data from the entire crustal seismogenic zone are mixed. Thus, we mapped T_L (P_L) separately for the top 15 km and the rest of the depth range, as well as for the entire seismogenic crust. The resulting T_L- and P_L-maps show that approximately 12% of the total area shows anomalously short recurrence times. Out of six shallow target events with M≥6.5 and which occurred since 1890, five are located within the anomalous areas with T_L <450 years. We interpret this to mean that areas with anomalously short T_L map asperities, which are more likely than other areas to generate future target events. The probability that this result is due to chance is vanishingly small. The great Kanto rupture of 1923 appears to have initiated in the most significant asperity we mapped in the study area. One anomaly is located in the northeastern part of the area of the proposed future rupture of the Tokai earthquake, and another one at its southwestern corner. The absolute values of T_L calculated are uncertain because they depend on the size of the volume used for the calculation.