Global trench migration velocities and slab migration induced upper mantle volume fluxes: Constraints to find an Earth reference frame based on minimizing viscous dissipation

Since the advent of plate tectonics different global reference frames have been used to describe the motion of plates and trenches. The difference in plate motion and trench migration between different reference frames can be substantial (up to 4 cm/yr). This study presents an overview of trench migration velocities for all the mature and incipient subduction zones on Earth as calculated in eight different global reference frames. Calculations show that, irrespective of the reference frame: (1) trench retreat always dominates over trench advance, with 62–78% of the 244 trench segments retreating; (2) the mean and median trench velocity are always positive (retreating) and within the range 1.3–1.5 cm/yr and 0.9–1.3 cm/yr, respectively; (3) rapid trench retreat is only observed close to lateral slab edges (< 1500 km); and (4) trench retreat is always slow far from slab edges (> 2000 km). These calculations are predicted by geodynamic models with a varying slab width, in which plate motion, trench motion and mantle flow result from subduction of dense slabs, suggesting that trench motion is indeed primarily driven by slab buoyancy forces and that proximity to a lateral slab edge exerts a dominant control on the trench migration velocity. Despite these four general conclusions, significant differences in velocities between such reference frames remain. It is therefore important to determine which reference frame most likely describes the true absolute velocities to get an understanding of the forces driving plate tectonics and mantle convection. It is here proposed that, based on fluid dynamic considerations and predictions from geodynamic modelling, the best candidate is the one, which optimises the number of trench segments that retreat, minimizes the trench–perpendicular trench migration velocity (v_T) in the centre of wide (> 4000 km) subduction zones, maximizes the number of retreating trench segments located within 2000 km of the closest lateral slab edge, minimizes the average of the absolute of the trench–perpendicular trench migration velocity (|v_T|) for all subduction zones on Earth, and minimizes the global upper mantle toroidal volume flux (_To) that results from trench migration and associated lateral slab migration (i.e. slab rollback or slab roll-forward). Calculations show that these conditions are best met in one particular Indo-Atlantic hotspot reference frame, where 75% of the subduction zones retreat, v_T in the centre of wide subduction zones ranges between − 3.5 and 1.8 cm/yr, 83% of the trench segments located within 2000 km of the closest lateral slab edge retreat, the average of |v_T| is 2.1 cm/yr, and _To = 456 km³/yr (lower limit) and 539 km³/yr (upper limit). Inclusion of all the incipient subduction zones on Earth results in slightly greater fluxes of 465 km³/yr (lower limit) and 569 km³/yr (upper limit). It is also found that this reference frame is close to minimizing the total sub-lithospheric upper mantle volume flux (_K) associated with motion of continental keels located below the major cratons. It is stressed, however, that _K is an order of magnitude smaller than _To, and thus of subordinate importance. In conclusion, the Indo-Atlantic hotspot reference frame appears preferable for calculating plate velocities and plate boundary velocities.     

Caribbean island-arc rifting and back-arc basin development in the Late Cretaceous: Geochemical, isotopic and geochronological evidence from Central Hispaniola

We present new regional petrologic, geochemical, Sr–Nd isotopic, and U–Pb geochronological data on the Turonian–Campanian mafic igneous rocks of Central Hispaniola that provide important clues on the development of the Caribbean island-arc. Central Hispaniola is made up of three main tectonic blocks—Jicomé, Jarabacoa and Bonao—that include four broad geochemical groups of Late Cretaceous mafic igneous rocks: group I, tholeiitic to calc-alkaline basalts and andesites; group II, low-Ti high-Mg andesites and basalts; group III, tholeiitic basalts and gabbros/dolerites; and group IV, tholeiitic to transitional and alkalic basalts. These igneous rocks show significant differences in time and space, from arc-like to non-arc-like characteristics, suggesting that they were derived from different mantle sources. We interpret these groups as the record of Caribbean arc-rifting and back-arc basin development in the Late Cretaceous. The> 90 Ma group I volcanic rocks and associated cumulate complexes preserved in the Jicomé and Jarabacoa blocks represent the Albian to Cenomanian Caribbean island-arc material. The arc rift stage magmatism in these blocks took place during the deposition of the Restauración Formation from the Turonian–Coniacian transition (~ 90 Ma) to Santonian/Lower Campanian, particularly in its lower part with extrusion at 90–88 Ma of group II low-Ti, high-Mg andesites/basalts. During this time or slightly afterwards adakitic rhyolites erupted in the Jarabacoa block. Group III tholeiitic lavas represent the initiation of Coniacian–Lower Campanian back-arc spreading. In the Bonao block, this stage is represented by back-arc basin-like basalts, gabbros and dolerite/diorite dykes intruded into the Loma Caribe peridotite, as well as the Peralvillo Sur Formation basalts, capped by tuffs, shales and Campanian cherts. This dismembered ophiolitic stratigraphy indicates that the Bonao block is a fragment of an ensimatic back-arc basin. In the Jicomé and Jarabacoa blocks, the mainly Campanian group IV basalts of the Peña Blanca, Siete Cabezas and Pelona–Pico Duarte Formation, represent the subsequent stage of back-arc spreading and off-axis non-arc-like magmatism, caused by migration of the arc toward the northeast. These basalts have geochemical affinities with the mantle domain influenced by the Caribbean plume, suggesting that mantle was flowing toward the NE, beneath the extended Caribbean island-arc, in response to rollback of the subducting proto-Caribbean slab.     

GPS数据处理中的坐标系统一问题

针对CGCS2000坐标系的启用,探讨了GPS数据中存在的坐标系统一问题,这些问题主要有：GPS基线向量坐标框架的统一、基线向量历元的统一、基线向量影响的量级分析和起算坐标的兼容性分析,并给出了相应的数学模型及适用条件。     

Rapid vertical motions and formation of volcanic arc gaps: Plateau collision recorded in the forearc geological evolution (Costa Rica margin)

The collision of bathymetric features with modern convergent margins has been investigated with the full range of tools used in geosciences. Hence, a comprehensive picture exists of the characteristic effects of collision events on the evolution of convergent margins. In contrast, much less studies documented past collisions of bathymetric features with convergent margins, as colliding features were generally lost to subduction. The arc‐trench system of southern Central America provides modern and past textbook examples of active margin interaction with incoming bathymetric reliefs. Here, we propose a synthesis which combines basin and terrane analysis of the forearc of northern Costa Rica and takes up the challenge of documenting past episodes of plateau accretion to the active margin. As illustrated in modern examples, our study shows that kilometric uplift of the overriding plate and termination of the volcanic arc activity are the most profound effects of colliding/accreting oceanic plateaus. Kilometric uplift of the forearc is documented by short‐lived (ca. 3 m.y.) occurrences of shallow‐water deposits in an overall deep‐water forearc record. These shallow deposits contain material reworked from underlying sedimentary and basement lithologies. The development of spatial gaps in arc volcanism is deduced from the transition from arc‐derived turbidites to pelagic sediments. Eventually, end of the collision event is evidenced by the subsidence of the whole forearc to deep‐water environments. Basin subsidence is accompanied or followed by renewed volcanic arc activity and coeval arc‐derived sedimentation, which may occur 1–7 m.y. after plateau collision. These past episodes of plateau accretion are archetypal for the following reasons: (a) they may be studied in outcrop, whereas most of the modern collisions of plateaus largely occur underwater; (b) no tectonic or metamorphic imprint has significantly complicated the forearc geological record; (c) the colliding feature and the sedimentary response to its collision are both preserved in the forearc geology; (d) they may be used as analogues for any setting where a bathymetric feature is suspected to have caused rapid forearc uplift and cessation of the volcanic arc activity.     

台站插值技术对检测固体地球变化的影响分析

搜集了全球约4 000个GNSS测站实测数据(站坐标和速度场信息),对固体地球体积及其变化率进行计算分析。采用板块运动模型插值方法对台站进行插值,使全球测站分布更加均匀。利用Delaunay数学模型计算固体地球体积及其变化率的结果发现:基于板块运动模型进行台站插值可以有效地提高计算结果的误差精度,不会对地球最终形态变化结论产生影响,台站数量及其分布均匀程度对计算最终结果影响较大。     

Tracking the paleogene India-Arabia plate boundary

The location of the India-Arabia plate boundary prior to the formation of the Sheba ridge in the Gulf of Aden is a matter of debate. A seismic dataset crossing the Owen Fracture Zone, the Owen Basin, and the Oman Margin was acquired to track the past locations of the India-Arabia plate boundary. We highlight the composite age of the Owen Basin basement, made of Paleocene oceanic crust drilled on its eastern part, and composed of pre-Maastrichtian continental and oceanic crust overlaid by ophiolites emplaced in Early Paleocene on its western side. A major fossil transform fault system crossing the Owen Basin juxtaposed these two slivers of lithosphere of different ages, and controlled the uplift of marginal ridges along the Oman Margin. This transform system deactivated ∼40 Myrs ago, coeval with the onset of ultra-slow spreading at the Carlsberg Ridge. The transform boundary then jumped to the edge of the present-day Owen Ridge during the Late Eocene-Oligocene period, before seafloor spreading began at the Sheba Ridge. This migration of the plate boundary involved the transfer of a part of the Indian oceanic lithosphere formed at the Carlsberg Ridge to Arabia. This Late Eocene-Oligocene tectonic episode at the India-Arabia plate boundary is synchronous with a global plate reorganization event corresponding to geological events at the Zagros and Himalaya belts. The Owen Ridge uplifted later, in Late Miocene times, and is unrelated to any major migration of the India-Arabia boundary.     

Free vibrations of simply supported double plate on two models of elastic soils

In this paper the free vibrations of simply supported rectangular plates, resting on two different models of soils, are considered. The first model called Hetényi, by the name of the deviser, assumes a continuous plate, embedded in a Winkler‐type soil to realize the foundation partial continuity. The upper plate rests on a multiple layer characterized by a Winkler‐type soil, a continuous plate and another Winkler‐type soil in sequence. The two Winkler‐type soils have different modules. In the second model, the plate will be embedded in two layers of soils, whose behavior is similar to that of the Pasternak–Kerr‐type soil. The two models have been already used for the study of the double beam system. The free motions, in both cases, are described by a homogeneous set of partial differential equations, based on Kirchhoff–Love theory. Next, the homogeneous equations of motion are solved by using the classical Navier method. The free frequencies and associated vibration mode shapes of double plate system are found and numerical examples are illustrated to compare the two models. Copyright © 2008 John Wiley & Sons, Ltd.     

Large‐scale structures on Gazelle Peninsula,New Britain: Implications for the evolution of the New Britain arc

Geological mapping of fault systems on the Gazelle Peninsula, eastern New Britain arc, combined with a reinterpretation of existing sea floor data indicate that faults previously thought to be a possible location of the boundary between the North and South Bismarck Plates, do not appear to be directly related to the plate boundary spreading centres and transform faults in the 3.5 Ma Manus Basin. Structure on the Gazelle Peninsula is dominated by the Mediva Fault (new name) and the Wide Bay Fault System, both north‐northwest trending, deep‐seated features. The Mediva Fault, an element of the Baining Mountain Horst and Graben Zone, is an extensional structure which has focused Middle Miocene intrusive activity, controlled Mio‐Pliocene sedimentation in the central Gazelle Peninsula, and displaced Quaternary volcanic deposits. The Wide Bay Fault System has been active since at least the Late Oligocene. One hundred kilometres of sinistral strike‐slip motion is likely on this fault since at least the late Middle Miocene, moving the Gazelle Peninsula in a north‐northwest direction with respect to the remainder of New Britain. The nature and timing of movements along these two major structures indicate that some other major tectonic process has operated (and presently continues) in this region of the New Britain arc to create these structures.     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

Evidence for multiple burial–partial exhumation cycles from the Onodani eclogites in the Sambagawa metamorphic belt,central Shikoku,Japan

Eclogites from the Onodani area in the Sambagawa metamorphic belt of central Shikoku occur as layers or lenticular bodies within basic schists. These eclogites experienced three different metamorphic episodes during multiple burial and exhumation cycles. The early prograde stage of the first metamorphic event is recorded by relict eclogite facies inclusions within garnet cores (X_Sps 0.80–0.24, X_Alm 0–0.47). These inclusions consist of relatively almandine‐rich garnet (X_Sps 0.13–0.24, X_Alm 0.36–0.45), aegirine‐augite/omphacite (X_Jd 0.08–0.28), epidote, amphiboles (e.g. actinolite, winchite, barroisite and taramite), albite, phengite, chlorite, calcite, titanite, hematite and quartz. The garnet cores also contain polyphase inclusions consisting of almandine‐rich garnet, omphacite (X_Jd 0.27–0.28), amphiboles (e.g. actinolite, winchite, barroisite, taramite and katophorite) and phengite. The peak P–T conditions of the first eclogite facies metamorphism are estimated to be 530–590 °C and 19–21 kbar succeeded by retrogression into greenschist facies. The second prograde metamorphism began at greenschist facies conditions. The peak metamorphic conditions are defined by schistosity‐forming omphacites (X_Jd ≤ 49) and garnet rims containing inclusions of barroisitic amphibole, phengite, rutile and quartz. The estimated peak metamorphic conditions are 630–680 °C and 20–22 kbar followed by a clockwise retrograde P–T path with nearly isothermal decompression to 8–12 kbar. In veins cross‐cutting the eclogite schistosity, resorbed barroisite/Mg‐katophorite occurs as inclusions in glaucophane which is zoned to barroisite, suggesting a prograde metamorphism of the third metamorphic event. The peak P–T conditions of this metamorphic event are estimated to be 540–600 °C and 6.5–8 kbar. These metamorphic conditions are correlated with those of the surrounding non‐eclogitic Sambagawa schists. The Onodani eclogites were formed by subduction of an oceanic plate, and metamorphism occurred beneath an accretionary prism. These high‐P/T type metamorphic events took place in a very short time span between 100 and 90 Ma. Plate reconstructions indicate highly oblique subduction of the Izanagi plate beneath the Eurasian continent at a high spreading rate. This probably resulted in multiple burial and exhumation movements of eclogite bodies, causing plural metamorphic events. The eclogite body was juxtaposed with non‐eclogitic Sambagawa schists at glaucophane stability field conditions. The amalgamated metamorphic sequence including the Onodani eclogites were exhumed to shallow crustal/surface levels in early Eocene times (c. 50 Ma).