GPS RTK测量技术的应用与体会

GPS RTK测量技术以快捷、精准的特点广泛地被测量人员所接受，但是在到量中仍存在一些问题需要研究。本文介绍了几种经常遇到的棘手问题及其解决方法。     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

北秦岭古聚会带壳幔再循环

以同构造期代表古洋壳残片的蛇绿岩及产于古岛弧的玄武岩为基础，通过Ｎｄ，Ｐｂ同位素与微量元素示踪及岩浆源区分析，揭示出北秦岭元古宙上地幔以强亏损（εＮｄ（ｔ）＋６．３～＋７．３）和高的Ｙｂ／Ｈｆ，Ｎｂ／Ｌａ和Ｔｈ／Ｌａ比值为特征，北秦岭地壳和上地幔明显具有Ｐｂ同位素比值高的特征．北秦岭丹凤群岛弧火山岩、二郎坪群弧后玄武岩以及松树沟蛇绿岩中变拉斑玄武岩εＮｄ（ｔ）、放射成因Ｐｂ同位素、Ｙ／Ｔｂ和Ｔｉ－ＭｇＯ研究表明，本区玄武岩存在两类性质不同的岩浆源．一类与亏损的北秦岭岩石圈上地幔源区有关；另一类与携带海洋沉积物的洋壳板块俯冲参与有关．由此，论证了北秦岭古聚会带壳幔之间物质再循环     

Testing an earthquake prediction algorithm

A test to evaluate earthquake prediction algorithms is being applied to a Russian algorithm known asM8 TheM8 algorithm makes intermediate term predictions for earthquakes to occur in a large circle, based on integral counts of transient seismicity in the circle. In a retroactive prediction for the period January 1, 1985 to July 1, 1991 the algorithm as configured for the forward test would have predicted eight of ten strong earthquakes in the test area. A null hypothesis, based on random assignment of predictions, predicts eight earthquakes in 2.87% of the trials. The forward test began July 1, 1991 and will run through December 31, 1997. As of July 1, 1995, the algorithm had forward predicted five out of nine earthquakes in the test area, which success ratio would have been achieved in 53% of random trials with the null hypothesis.     

Seafloor spreading, obduction and triple junction tectonics of the Mineoka Ophiolite, Central Japan

The Tertiary Mineoka ophiolite occurs in a fault zone at the intersection of the Honshu and Izu forearcs in central Japan and displays structural evidence for three major phases of deformation: normal and oblique-slip faults and hydrothermal veins formed during the seafloor spreading evolution of the ophiolite at a ridge-transform fault intersection. These structures may represent repeated changes in differential stress and pore-fluid pressures during their formation. The second series of deformation is characterized by oblique thrust faults with Riedel shears and no significant mineral veining, and is interpreted to have resulted from transpressional dextral faulting during the obduction of the ophiolite through oblique convergence and tectonic accretion. This deformation occurred at the NW corner of a TTT-type (trench–trench–trench) triple junction in the NW Pacific rim before the middle Miocene. The third series of deformation of the ophiolite is marked by contractional and oblique shear zones, Riedel shears, and thrust faults that crosscut and offset earlier structures, and that give the Mineoka fault zone its lenticular (phacoidal) fabric at all scales. This deformation phase was associated with the establishment and the southward migration of the TTT Boso triple junction and with the kinematics of oblique subduction and forearc sliver fault development. The composite Mineoka ophiolite hence displays rocks and structures that evolved during its complex geodynamic history involving seafloor spreading, tectonic accretion, and triple junction evolution in the NW Pacific Rim.     

Determining the maximum degree of harmonic coefficients in geopotential models by Monte Carlo methods

Random errors for the harmonic coefficients of a geopotential model are generated from the matrix of normal equations by a parallel computer applying the Gibbs sampler. This leads to random values for the harmonic coefficients. They are transformed by nonlinear, quadratic transformations to random values for the square roots of degree variances, of mean squares of geoid undulations and gravity anomalies. The expected values of these quantities are not equal to the values of these quantities computed by the estimated harmonic coefficients, due to correlations and errors in the estimation. By hypothesis tests estimated harmonic coefficients distorted by correlations and errors are detected. Applying the tests to the geopotential model ITG-CHAMP01 of the Institute of Theoretical Geodesy in Bonn it is concluded that above the degree 62 the harmonic coefficients cannot add any information to the geopotential model.     

Crustal and upper mantle seismic structure of the Australian Plate, South Island, New Zealand

Seismic reflection and refraction data were collected west of New Zealand's South Island parallel to the Pacific–Australian Plate boundary. The obliquely convergent plate boundary is marked at the surface by the Alpine Fault, which juxtaposes continental crust of each plate. The data are used to study the crustal and uppermost mantle structure and provide a link between other seismic transects which cross the plate boundary. Arrival times of wide-angle reflected and refracted events from 13 recording stations are used to construct a 380-km long crustal velocity model. The model shows that, beneath a 2–4-km thick sedimentary veneer, the crust consists of two layers. The upper layer velocities increase from 5.4–5.9 km/s at the top of the layer to 6.3 km/s at the base of the layer. The base of the layer is mainly about 20 km deep but deepens to 25 km at its southern end. The lower layer velocities range from 6.3 to 7.1 km/s, and are commonly around 6.5 km/s at the top of the layer and 6.7 km/s at the base. Beneath the lower layer, the model has velocities of 8.2–8.5 km/s, typical of mantle material. The Mohorovicic discontinuity (Moho) therefore lies at the base of the second layer. It is at a depth of around 30 km but shallows over the south–central third of the profile to about 26 km, possibly associated with a southwest dipping detachment fault. The high, variable sub-Moho velocities of 8.2 km/s to 8.5 km/s are inferred to result from strong upper mantle anisotropy. Multichannel seismic reflection data cover about 220 km of the southern part of the modelled section. Beneath the well-layered Oligocene to recent sedimentary section, the crustal section is broadly divided into two zones, which correspond to the two layers of the velocity model. The upper layer (down to about 7–9 s two-way travel time) has few reflections. The lower layer (down to about 11 s two-way time) contains many strong, subparallel reflections. The base of this reflective zone is the Moho. Bi-vergent dipping reflective zones within this lower crustal layer are interpreted as interwedging structures common in areas of crustal shortening. These structures and the strong northeast dipping reflections beneath the Moho towards the north end of the (MCS) line are interpreted to be caused by Paleozoic north-dipping subduction and terrane collision at the margin of Gondwana. Deeper mantle reflections with variable dip are observed on the wide-angle gathers. Travel-time modelling of these events by ray-tracing through the established velocity model indicates depths of 50–110 km for these events. They show little coherence in dip and may be caused side-swipe from the adjacent crustal root under the Southern Alps or from the upper mantle density anomalies inferred from teleseismic data under the crustal root.