全球五大板块的运动和形变──用卫星激光测距数据导出的站速度分析

根据德国卫星激光测距（ＳＬＲ）数据分析中心ＧＦＺ对１９８０年１月－１９９１年７月获取的ＳＬＲ观测数据处理后得到的４１个ＳＬＲ站的站速度，解算了北美、欧亚、太平洋、南美和澳大利亚板块之间的相对运动欧拉矢量，得到了第１个ＳＬＲ实测的板块运动模型ＳＰＭＭ１．与地学板块模型ＲＭ２和ＮＵＶＥＬ－１的比较指出，ＳＰＭＭ１大体上与地学模型一致，与ＮＵＶＥＬ－１更为接近；ＳＰＭＭ１的欧亚与北美板块相对运动欧拉极与ＮＵＶＥＬ－１的相应极很接近，但旋转速率明显偏小．还分析了各板块上ＳＬＲ站的局部形变特征．欧亚板块东部和西部存在１０－２０ｍｍ／ａ的相对运动；板块边界附近ＳＬＲ站的残差站速度基本上反映了该边界的构造形变特征．     

Determination of Euler parameters of Philippine Sea plate and the inferences

Euler vectors of 12 plates, including Philippine Sea plate (PH), relative to a randomly fixed Pacific plate(PA) were determined by inverting the 1122 data from NUVEL-1 global plate motion model, earthquake slip vectors along Philippine Sea plate boundary, and GPS observed velocities. Euler vectors of Philippine Sea plate relative to adjacent plates are also gained. Our results are well consistent with observed data and can satisfy the geological and geophysical constraints along the Caroline(CR)-PH and PA-CR boundaries. Deformation of Philippine Sea plate is also discussed by using the plate motion Euler parameters.     

东亚北部地区现代板块构造的运动学分析

从原欧亚板块划出阿穆利亚微板块,并建立一个包括太平洋、北美、欧亚、鄂霍次克及菲律宾海等板块在内的６板块系统,利用地震滑移矢量,转换断层走向,洋中脊扩展速率,以及最新的ＧＰＳ观测资料,采用Ｍｏｎｔｅ－Ｃａｒｌｏ盒子法,反演得到该６板块系统的欧拉运动矢量．本文得到的有关欧拉矢量为：ＡＭ-ＥＵ（６０．４２°Ｎ,１２３．２５°Ｅ,０．０２５°／Ｍａ）,ＡＭ-ＯＫ（５３．２０°Ｎ,１４１．９５°Ｅ,０．４７６°／Ｍａ）以及ＡＭ-ＰＨ（５０．８４°Ｎ,１５８.１３°Ｅ,１.２０４°／Ｍａ）,ＡＭ－ＥＵ的欧拉运动极点位于贝加尔隆起的东北及斯塔诺夫山脉的西北,并给出在贝加尔隆起有０．４－０．７ｍｍ／ａ的扩张,与从地质学证据得到的估计结果一致;由ＡＭ－ＯＫ及ＡＭ－ＰＨ欧拉运动矢量可以分别得到日本海东缘的６-１５ｍｍ／ａ及日本南海地槽带的５１－６６ｍｍ／ａ收敛速率,与最新的ＧＰＳ观测结果并不一致．因此,东亚北部地区存在独立的南中国微板块及其它亚尺度的微板块,可以更好地描述该地区现今板块构造的运动规律．     

Relative Motion Between the Rivera and North American Plates Determined from the Slip Directions of Earthquakes

So far, the direction and rate of relative motion between the Rivera and the North American plates (RIV-NAM) has been determined by the combination of two Euler poles: Rivera (RIV), with respect to Pacific (PAC), and PAC with respect to North America. Here, we estimate the relative motion of this plate pair (RIV-NAM) assuming that the horizontal projection of the direction of slip of the earthquakes occurring on the RIV-NAM boundaries reflect their relative plate motion. A catalog of earthquakes for which focal mechanisms are reported since 1976 is used in the analysis. Earthquakes were considered in the three segments of the RIV-NAM plate boundary: the subduction zone of the Rivera plate beneath the Jalisco block, the Tres Marias Escarpment and the events associated with the Tamayo Fracture Zone. The best fitting Euler pole is determined using a grid search of 64 potential poles. The slip direction predicted for each grid point is compared to the slip direction of the focal mechanisms of the earthquakes on the plate boundary. The best fitting Euler pole, determined in a root mean square sense (RMS), is located at 21.8°N, 107.6°W. A rate of rotation of 5.3°/year is estimated assuming the seismic earthquake cycle of the 1932 and 1995 great earthquakes represents a lower bound of the rate of plate motion in the subduction zone. The best fitting Euler pole shows that the subduction of the Rivera plate takes place in a direction perpendicular to the trench with a relative velocity of 4.3 cm/year, offshore Manzanillo. The rate of relative motion RIV-NAM decreases from SE to NW. North of approximately 21°N, the subduction of the Rivera plate becomes oblique to the trench and the relative velocity between the two plates decreases to an average of 1.9 cm/year. This slow rate of convergence may explain the rapid decrease of seismicity in the trench and the apparent absence of large earthquakes in this region. In the Tres Marias Escarpment, our best-fitting pole suggests that subduction stops, giving way to high-angle reverse faulting perpendicular to the Tres Marias Escarpment, in agreement with the reverse faulting earthquakes occurring here. To the north of 22.5°N, the slip predicted by the best-fitting pole suggests right-lateral faulting in a direction parallel to the Tamayo Fracture Zone, at a very low velocity (0.5–1.0 cm/year). The best fitting Euler pole determined here lies very close to the RIV-NAM plate boundary in the vicinity of the Tamayo Fracture Zone. This location of our best fitting Euler pole explains the low relative plate velocity, the relatively low level of seismic activity and the presence of a broad zone of deformation that accommodates the RIV-NAM motion.     

Early Tertiary rupture of the Pacific plate: 1700 km of dextral offset along the Emperor trough-Line Islands lineament

Between 67 and ～40 Ma ago a northwest-southeast-trending fracture system over 8000 km long split the Pacific plate and accumulated at least 1700 km of dextral offset between the east and west portions. This system, here named the Emperor fracture zone (EFZ) system, consisted of several segments, one along the present trace of the Emperor trough and another along the Line Islands, joined by short spreading ridges. The EFZ terminated at its northern end against the Kula-Pacific ridge, and at its southern end in a ridge-transform system, called the Emperor spreading system, which extended to the west, north of Australia.The finite angular velocity vector describing the relative motion between the East and West Pacific plates is ～0.6°/Ma about a pole at 36°N, 70°W. This vector, added to the known Early Tertiary motion of the Pacific plate with respect to the global hotspot reference frame, accounts in large part for the NNW trend of the Emperor seamount chain relative to the WNW Hawaiian trend, without violation of the integrity of the Antarctic plate. The Meiji-Emperor and Emperor-Hawaiian bends date, respectively, the initiation (～67 Ma ago) and cessation (～40 Ma ago) of seafloor spreading on the Emperor spreading system.The postulated Early Tertiary relative motion along the EFZ between the East and West Pacific plates explains (1) the present misalignment of the two sets of magnetic bights of the Pacific, (2) the abrupt truncation of eastern Pacific bathymetric lineaments against the Emperor trough and Line Islands, (3) the contrast in paleolatitude between the eastern and western Pacific as indicated by paleomagnetic and sedimentologic studies, and (4) the anomalous gravity signature of the central Hawaiian ridge that indicates that the ridge loaded thin hot lithosphere.     

Secular polar motion,crustal movements,and international latitude service observations

Summary The latitude variations observed by optical astrometry at five International Latitude Service (ILS) stations since the beginning of the century are analyzed to estimate the secular motion of the pole and the correction of the proper motions of the observed stars. The geophysical model of tectonic plate motions NUVEL-1 is used to correct the observed latitude drifts, thus referring the derived secular polar motion to the "no net rotation" terrestrial reference system. It is shown that one of the ILS stations, Ukiah (near the U.S. west coast), very probably exhibits an anomalous drift which differs significantly from the motion of both the Pacific and North American plates.     

利用海洋重力场变化分析洋底板块运动

洋底板块运动是地球动力学和全球变化研究的重要内容.本文根据质量迁移与地球外部重力场变化的对应关系,利用不同时期测高资料推算的1995—2019全球海洋重力场变化结果,反演分析全球洋底板块运动特征.研究表明,板块汇聚边界、板块内无震海岭、海山群、断裂带等区域重力异常变化显著,而在板块离散边界无明显变化趋势;西南印度洋中脊、大西洋中脊、中印度洋中脊等地区重力异常垂直梯度变化显著,且在西太平洋俯冲带、部分海岭区域也存在明显变化,其空间分布与地形基本吻合.海洋重力场变化整体上准确反映了全球洋底板块构造运动.相较于重力异常变化反演结果,重力垂直梯度的变化能够更为准确地反映洋底板块运动特征,特别是在洋中脊区域,扩张速率越小,垂直重力梯度变化越显著.此外,详细讨论了测高海洋重力场不确定因素对洋底板块运动分析结果的影响,海面坡度改正是主要因素之一.     

南半球减速膨胀的定量分析

利用空间大地测量技术的长期观测资料，得出南半球纬线圈纬线长变化率和全球活动板块边缘扩张与汇聚运动速率，并与3Ma平均地质地磁模型NUVEL1A的估算结果进行比较：(1)空间大地测量技术测得南半球纬线圈纬线变化率均为正值；(2)南半球测站的垂向运动除了赤道附近几个测站下沉，其余91髎的台站全上升；（3）南半球相邻板块的现今汇聚和扩张运动速率均比3Ma平均地质模型NUVEL1A估值小，而北半球相邻板块的汇聚和扩张运动速率没有系统性的变化. 这些实测结果反映了南半球纬线圈方向在减速伸展，南北方向在减速拉伸，即南半球在减速膨胀.     

The present relative motion between Africa and Antarctica

Movement between the Africa and Antarctica plates is at present accomplished by sea-floor spreading on the Southwest Indian Ocean Ridge. This movement may be described in terms of an angular rotation vector. Bathymetric and magnetic observations from marine geophysical surveys near the Bouvet triple junction, at 52°S, 15°E and in the environs of the Prince Edward Islands are combined with spreading azimuths derived from earthquake fault plane solutions to define this vector. The rotation pole which describes the motion is located at 10.7°N, 40.9°W and the angular velocity is 1.44 × 10^?7 deg/yr. This pole is significantly different from some other poles obtained by global closure or vector addition. The possibility that the differences may be due to Africa being split into two plates is investigated but there would have to be convergence across the African Rift system for this possibility to be true. Closure of the vector velocity triangle around the Central Indian triple junction is checked by using the pole derived in this study and published poles and rates for the Africa/India and Antarctica/India motions to determine this triangle. The triangle is found to close when errors in the Africa/India and Antarctica/India motions are taken into account. This suggests that it is errors in the data that cause the differences between the observed and predicted poles.     

America-Eurasia plate boundary in eastern Asia and the opening of marginal basins

The opening of the Arctic Ocean during the past 55 Ma resulted in relative rotation of America with respect to Eurasia about a pole located in eastern Siberia, near the plate boundary. The extensional plate boundary enters deep inside the Eurasian continent up to the rotation pole. On the opposite side of the pole, on the Pacific side of the plate boundary, compressive tectonics are recorded along Sakhalin and Hokkaido. From the Oligocene to Middle Miocene, the relative movement was accommodated by strike-slip motion along Sakhalin and Hokkaido although the rotation pole was not located at a significatively different position from now. In this paper we explain this by independent motion of the southernmost tip of the American plate towards the Pacific margin which behaves as a free boundary. This oceanward motion resulted in an extension of the American plate giving rise to the wedge structure of the Okhotsk Sea. The Japan Sea opened as a pull-apart basin along the strike-slip boundary; finally the increasing extension in the Okhotsk Sea led to the opening of the oceanic Kuril Basin.     

Tectonics of the Easter plate

We present a new model for the Easter plate in which rift propagation has resulted in the formation of a rigid plate between the propagating and dying ridges. We use the distribution of earthquakes, eleven new focal mechanisms, and existing bathymetric and magnetic data to describe the tectonics of this area. Both the Easter-Nazca and Easter-Pacific Euler poles are sufficiently close to the Easter plate to cause rapid changes in rates and directions of motion along the boundaries. The east and west boundaries are propagating and dying ridges; the southwest boundary is a slow-spreading ridge and the northern boundary is a complex zone of convergent and transform motion.The Easter plate may reflect the tectonics of rift propagation on a large scale, where rigid plate tectonics requires boundary reorientation. We use simple schematic models to illustrate the general features and processes which occur at plates resulting from large-scale rift propagation.     

Compression of the North Hemisphere derived from space geodesy

Introduction With the development of global tectonics and overall detections for global tectonics with multi-geophysical methods, ones can roundly study on the geological tectonics of sampling and magnetic stripe image, so as to summarize and interpret the geometrical and kinematical charac-teristics for the distribution of the ocean and the land, and spreading state of the global tectonics in a global scale. From a comprehensive view, the South and North hemispheres are clearly unsym-metrical…     

Statistical tests of additional plate boundaries from plate motion inversions

We have investigated the application of the F-ratio test, a standard statistical technique, to the results of relative plate motion inversions. The method tests whether the improvement in fit of the model to the data resulting from the addition of another plate to the model is greater than that expected purely by chance. This approach appears to be useful in determining whether additional plate boundaries are justified. We confirm previous results favoring separate North American and South American plates with a boundary located between 30°N and the equator. Using Chase's global relative motion data, we show that in addition to separate West African and Somalian plates, separate West Indian and Australian plates, with a best-fitting boundary between 70°E and 90°E, can be resolved. These results are generally consistent with the observation that the Indian plate's internal deformation extends somewhat westward of the Ninetyeast Ridge. The relative motion pole is similar to Minster and Jordan's and predicts the NW-SE compression observed in earthquake mechanisms near the Ninetyeast Ridge.     

南大西洋洋中脊板块构造运动过程:古地磁的约束

海底磁异常的形态与洋中脊两侧板块的微运动或变形密切相关.因此,这方面的研究可为确定板块运动的演化历史、小尺度的动力学过程以及洋中脊分段的机制等提供重要约束.本文对南大西洋一段洋中脊(31°S—34.5°S)两侧的磁异常的偏度进行了系统研究.结果表明研究区域内扩张方向并不总是垂直于洋中脊走向,并且研究区域不同剖面的扩张方向也不一致,具体表现为从北向南,平均扩张方向逐渐增加,依次为33.6°±5.3°、62.8°±13.0°以及94.3°±8.0°.这表明洋中脊的倾斜扩张机制具有复杂性,初步解释应该与转换断层的剪切应力增加有关.深部辉长岩层倾斜和扩张速率不对称性对海底磁异常偏度的影响值得深入研究.另外,由北向南确定的欧拉极向东移动,表明洋中脊两侧的板块在6.5 Ma期间存在剧烈形变.     

Global tectonics and the plate motion obtained from the ITRF97 station velocity vectors

~~Global tectonics and the plate motion obtained from the ITRF97 station velocity vectors@马宗晋 @任金卫 @张进~~     

Hotspots,mantle convection and plate tectonics: A synthetic calculation

A model is developed that unifies vigorous hotspots with global-scale mantle convection and plate tectonics. The convection dynamics are assumed to generate flow patterns that emerge as closely packed polygonal cells in approaching the asthenosphere, and whose geometry is completely determined by a defining set of vigorous hotspots. Overlying viscously coupled rigid plates are driven with unique velocities (Euler vectors) at which the area integral of the shear forces is zero; these velocities are dynamically stable. The computed plate velocities, resulting from convection based on 15 hotspots, are compared with the velocities of plate motion models AM1-2 (Minster andJordan, 1978) and HS-NUVEL1 (Gripp andGordon, 1990), which combine transform fault geometries, magnetic anomalies and seismic data. The comparison shows a striking agreement for a majority of the plates. Geophysical implications of this numerical exercise are discussed.     

Pacific Plate motion and undulations in geoid and bathymetry

Previous studies have shown that the Pacific geoid and gravity fields exhibit lineated anomalies, trending approximately in the direction of absolute plate motion over the underlying mantle. Because the undulations obliquely cross fracture zones they have often been attributed a convective origin. Recently, lithospheric boudinage caused by diffuse extension has been proposed as a possible mechanism. We have examined the undulations in the free-air anomalies, geoid and bathymetry over a portion of the Pacific Plate to determine quantitatively how the undulations are related to plate motion. We compare the observed data to an axisymmetric, sinusoidal undulation defined in an arbitrary frame of reference; in particular, we seek the north pole of this reference frame that maximizes the correlation between data and model. Poles that are close to the Pacific hotspot pole represent copolar undulations possibly related to plate motion. The distance between the best-fitting poles and the hotspot pole is determined as a function of undulation wavelength and reveals several minima (with distance < 10°) for discrete geoid wavebands centered on wavelengths of 160 km, 225 km, 287 km, 400 km, 660 km, 850 km, 1000 km and 1400 km. Bathymetry data have copolar bathymetric expressions as well, giving an implied admittance of 2–3 m/km. The most co-polar geoid/bathymetry undulations (with poles within 2–3° of the average Pacific Euler pole) have wavelengths of 280 km and 1050 km, respectively. The latter could have a convective origin or be related to the spacing of hotspot swells. The former may reflect lithospheric boudinage formed in response to diffuse extension, but could also have a dynamic origin since flexural dampening may only have attenuated the bathymetric amplitude by 50% or less. Radiometric dating of volcanic ridges found in the troughs of prominent gravity lineations gives ages that correlate well with documented changes in Pacific and Indo/Australian Plate motion, suggesting the ridges formed in response to intermittent plate boundary stresses and not as a direct consequence of small-scale convection or diffuse extension.     

Petrologic consequences of rift propagation on oceanic spreading ridges

The production of anomalously differentiated lava compositions at several mid-ocean spreading centers can be attributed to magmatic processes associated with propagating rifts. The degree of differentiation attained by magmas beneath oceanic spreading ridges depends mainly on the balance between cooling rate and the supply rate of new magma to shallow chambers. Low supply rates and moderate cooling rates allow advanced degrees of closed-system fractionation to occur. High supply rates result in open systems in which magma compositions are buffered by frequent replenishment with new hot magma. Propagating rift tips are a special class of ridge-transform intersection in which the balance between cooling and supply rates is conducive to the development of advanced degrees of differentiation over an expanded length of ridge. This balance is affected by the spreading rate, the propagation rate of the rift, the length of the bounding transform and proximity to hotspots. Maximum compositional variability and maximum degree of differentiation occur within 50 km of propagating rift tips and subsequently diminish with increasing distance. Rifts that propagate through plates in directions approximating their absolute motion relative to the lower mantle are characterized by the presence of anomalously differentiated lavas over longer ridge segments than are rifts that propagate against their absolute motion. Geochemical anomalies may persist, though changing in degree and extent, for several million years on ridge segments that stop propagating. The concept of “magnetic telechemistry” is generally supported by our study, but in the vicinity of hotspots, magnetic anomaly amplitude may be controlled more by bathymetric and/or thickened magnetic layer effects than by geochemistry.     

Inversion of magnetic anomalies and sea-floor spreading in the Cayman Trough

For some time, sea-floor spreading has been hypothesized for the Mid-Cayman Rise based on inferences from seismicity, heat flow, topography and plate geometry. Here we present magnetic anomaly inversions from which a reasonable record of sea-floor spreading emerges. We obtain total opening rates of 20 ± 2 mm/yr for 0–2.4 m.y. B.P. and 40 ± 2 mm/yr for 2.4–6.0 m.y. B.P. Data on the west flank extend the half-opening rate of 20 mm/yr back to 8.3 m.y. B.P. Spreading has been very nearly symmetric. These new observations place important constraints on plate tectonic reconstructions by defining the relative motion between the North American and Caribbean plates. They also shed some light on sea-floor spreading processes in which the spreading center is a secondary feature in the sense that it is over an order of magnitude shorter than the adjoining transform faults.     

Classical focal mechanism techniques for body waves

The initial motion of primary body waves and polarization directions of secondary body waves have been applied successfully to the study of the mechanism operating at the earthquake focus. Equal area plots of these body wave characteristics resulted in radiation patterns that were compared to patterns due to theoretical focal mechanism sources. Such an approach indicated that a double couple force is the source mechanism operating at earthquake foci. This can be physically represented by faulting at the focus. This seems likely because of the relation between earthquakes and fault motion in shallow earthquakes. The possibility of other mechanisms operating in deep focus earthquakes has not been ruled out.The technique of solving for fault plane orientation/motion in equal area plots ofP-wave first motion andS-wave polarization, has been a powerful tool in areas where fault motions cannot be directly observed in studying tectonics. Such an approach has been used to test the theory of global tectonics. This approach has resulted in the confirmation of the suggested mechanics of plate motion, and the results of plate motion, such as spreading of the sea floor from mid-oceanic ridges and underthrusting of lithospheric plates at the sites of oceanic trenches.