中国东南陆缘大地水准面的确定

利用我国东南地区１１８２２个重力数据、３０″×３０″数值地形模型和ＷＤＭ９４重力位模型,以Ｓｔｏｋｅｓ理论为基础,采用低通滤波方法处理远区影响,计算了中国东南陆缘２．５′×２．５′格网重力大地水准面,其结果与福建、山东的ＧＰＳ水准相比较,精度分别为±０．１２４ｍ和±０．０６２ｍ。     

A study on the combination of satellite, airborne, and terrestrial gravity data

Satellite gravity missions, such as CHAMP, GRACE and GOCE, and airborne gravity campaigns in areas without ground gravity will enhance the present knowledge of the Earths gravity field. Combining the new gravity information with the existing marine and ground gravity anomalies is a major task for which the mathematical tools have to be developed. In one way or another they will be based on the spectral information available for gravity data and noise. The integration of the additional gravity information from satellite and airborne campaigns with existing data has not been studied in sufficient detail and a number of open questions remain. A strategy for the combination of satellite, airborne and ground measurements is presented. It is based on ideas independently introduced by Sjöberg and Wenzel in the early 1980s and has been modified by using a quasi-deterministic approach for the determination of the weighting functions. In addition, the original approach of Sjöberg and Wenzel is extended to more than two measurement types, combining the Meissl scheme with the least-squares spectral combination. Satellite (or geopotential) harmonics, ground gravity anomalies and airborne gravity disturbances are used as measurement types, but other combinations are possible. Different error characteristics and measurement-type combinations and their impact on the final solution are studied. Using simulated data, the results show a geoid accuracy in the centimeter range for a local test area.     

Sensitivity of GPS and GLONASS orbits with respect to resonant geopotential parameters

The Center for Orbit Determination in Europe (CODE) has been involved in the processing of combined GPS/GLONASS data during the International GLONASS Experiment (IGEX). The resulting precise orbits were analyzed using the program SORBDT. Introducing one satellites positions as pseudo-observations, the program is capable of fitting orbital arcs through these positions using an orbit improvement procedure based on the numerical integration of the satellites orbit and its partial derivative with respect to the orbit parameters. For this study, the program was enhanced to estimate selected parameters of the Earths gravity field. The orbital periods of the GPS satellites are —in contrast to those of the GLONASS satellites – 2:1 commensurable (P _Sid:P _GPS) with the rotation period of the Earth. Therefore, resonance effects of the satellite motion with terms of the geopotential occur and they influence the estimation of these parameters. A sensitivity study of the GPS and GLONASS orbits with respect to the geopotential coefficients reveals that the correlations between different geopotential coefficients and the correlations of geopotential coefficients with other orbit parameters, in particular with solar radiation pressure parameters, are the crucial issues in this context. The estimation of the resonant geopotential terms is, in the case of GPS, hindered by correlations with the simultaneously estimated radiation pressure parameters. In the GLONASS case, arc lengths of several days allow the decorrelation of the two parameter types. The formal errors of the estimates based on the GLONASS orbits are a factor of 5 to 10 smaller for all resonant terms. AcknowledgmentsThe authors would like to thank all the organizations involved in the IGS and the IGEX campaign, in particular those operating an IGS or IGEX observation site and providing the indispensable data for precise orbit determination.     

Automatic Location of Swarm Earthquakes from Local Network Data

An advanced method of automated seismic phase picking and exact location and magnitude determination of swarm micro-earthquakes from local network data is presented. The phase picker is applied in two steps: first, S-wave groups are identified using a polarisation detector, and then corresponding P-wave groups are searched for. The times of maximum P- and S-amplitudes are then used as starting points for the determination of accurate P- and S-arrival times. The maximum S-wave amplitudes are utilised for determining local magnitudes. The whole procedure is checked by simultaneous preliminary hypocentre location providing estimates of local magnitudes and a compatibility check of the candidate P- and S-phases. The closest station to the earthquake cluster is used as a master, and the phase search at the remaining stations is governed by the P- and S-phases identified at the master station. Thanks to the use of apriori information on the approximate position of hypocentres, the procedure is also capable of picking the individual P- and S-phases of sequences of overlapping swarm events. The performance of the procedure was tested by comparison of the automatically and interactively created catalogues of the January 1997 NW-Bohemia micro-earthquake swarm. With stations located at epicentral distances between 0 and 20 km, the difference between hypocentre coordinates obtained by automatic and interactive processing did not exceed 80 m for 86% events. All events above magnitude 0.5 were identified, and the automatically determined polarity of first P-wave motion proved to be correct in 89% of them.     

地震速报参数不确定性的应急灾害损失快速评估模型

本文针对目前地震应急灾害损失快速评估中存在的问题,建立了考虑地震速报参数不确定性的灾害损失快速评估模型。并利用1990年来全国的81组速报震中与宏观震中数据。得到速报震中与宏观震中偏差的概率分布经验参数。     

高精度热电离质谱(TIMS)铀系法洞穴沉积物(石笋)年龄的研究

用高精度热电离质谱技术测定石笋的年龄,初步建立了热电离质谱铀系年龄的方法。 4个样品的测定结果复现性较好,铀含量、同位素比值及年龄值与标准值吻合,数据的精度大大提高,显示了质谱铀系年龄方法具有样品量少、时间短、精度高的优点。     

基于用户星精密轨道的TDRS卫星定轨

分析了TDRS卫星的轨道特性及传统的地基测距跟踪技术定轨精度不高的现状,研究了基于空基的用户星精密轨道的TDRS卫星定轨,解决了基于空基的一般GEO卫星定轨问题。     

Precise orbit determination for the GRACE mission using only GPS data

The GRACE (gravity recovery and climate experiment) satellites, launched in March 2002, are each equipped with a BlackJack GPS onboard receiver for precise orbit determination and gravity field recovery. Since launch, there have been significant improvements in the background force models used for satellite orbit determination, most notably the model for the geopotential. This has resulted in significant improvements to orbit accuracy for very low altitude satellites. The purpose of this paper is to investigate how well the orbits of the GRACE satellites (about 470 km in altitude) can currently be determined using only GPS data and based on the current models and methods. The orbit accuracy is assessed using a number of tests, which include analysis of orbit fits, orbit overlaps, orbit connecting points, satellite Laser ranging residuals and K-band ranging (KBR) residuals. We show that 1-cm radial orbit accuracy for the GRACE satellites has probably been achieved. These precise GRACE orbits can be used for such purposes as improving gravity recovery from the GRACE KBR data and for atmospheric profiling, and they demonstrate the quality of the background force models being used.     

Efficient satellite orbit modelling using pseudo-stochastic parameters

If the force field acting on an artificial Earth satellite is not known a priori with sufficient accuracy to represent its observations on their accuracy level, one may introduce so-called pseudo-stochastic parameters into an orbit determination process, e.g. instantaneous velocity changes at user-defined epochs or piecewise constant accelerations in user-defined adjacent time subintervals or piecewise linear and continuous accelerations in adjacent time subintervals. The procedures, based on standard least-squares, associated with such parameterizations are well established, but they become inefficient (slow) if the number of pseudo-stochastic parameters becomes large. We develop two efficient methods to solve the orbit determination problem in the presence of pseudo-stochastic parameters. The results of the methods are identical to those obtained with conventional least-squares algorithms. The first efficient algorithm also provides the full variance–covariance matrix; the second, even more efficient algorithm, only parts of it.     

Explicit formula for the geoid-quasigeoid separation

The explicit formula for the geoid-to-quasigeoid correction is derived in this paper. On comparing the geoidal height and height anomaly, this correction is found to be a function of the mean value of gravity disturbance along the plumbline within the topography. To evaluate the mean gravity disturbance, the gravity field of the Earth is decomposed into components generated by masses within the geoid, topography and atmosphere. Newton’s integration is then used for the computation of topography-and atmosphere-generated components of the mean gravity, while the combined solution for the downward continuation of gravity anomalies and Stokes’ boundary-value problem is utilized in computing the component of mean gravity disturbance generated by mass irregularities within the geoid. On application of this explicit formulism a theoretical accuracy of a few millimetres can be achieved in evaluation of the geoid-to-quasigeoid correction. However, the real accuracy could be lower due to deficiencies within the numerical methods and to errors within the input data (digital terrain and density models and gravity observations).