Finite element models to represent seismic activity of the Indian plate

Quantification of seismic activity is one of the most challenging problems faced by earthquake engineers in probabilistic seismic hazard analysis. Currently, this problem has been attempted using empirical approaches which are based on the regional earthquake recurrence relations from the available earthquake catalogue. However, at a specified site of engineering interest, these empirical models are associated with large number of uncertainties due to lack of sufficient data. Due to these uncertainties, engineers need to develop mechanistic models to quantify seismic activity. A wide range of techniques for modeling continental plates provides useful insights on the mechanics of plates and their seismic activity. Among the different continental plates, the Indian plate experiences diffused seismicity. In India, although Himalaya is regarded as a plate boundary and active region, the seismicity database indicates that there are other regions in the Indian shield reporting sporadic seismic activity. It is expected that mechanistic models of Indian plate, based on finite element method, simulate stress fields that quantify the seismic potential of active regions in India. This article explores the development of a finite element model for Indian plate by observing the simulated stress field for various boundary conditions, geological and rheological conditions. The study observes that the magnitude and direction of stresses in the plate is sensitive to these conditions. The numerical analysis of the models shows that the simulated stress field represents the active seismic zones in India.     

空间大地测量测定板块运动新进展

主要阐述利用空间大地测量测定和研究板块运动的研究成果及最新进展，并与地质学方法进行了比较，结果表明空间大地测量方法测定板块运动由单一手段发展为多种手段的综合，最新的研究是利用ITRF₉₆、ITRF₉₇地球参考架的站坐标及速率建立新的板块运动模型，这些模型和相同时期的地质模型，整体上一致但仍存在着差异。对于太平洋板块所有大地测量测定的运动速率均大于地质学模型，极位置的经度也有系统偏差。     

Relating rapid plate-motion variations to plate-boundary forces in global coupled models of the mantle/lithosphere system: Effects of topography and friction

Recent high-resolution models of past plate motions and their comparison with plate motion models inferred from space geodetic techniques reveal a number of short-term variations in global plate velocities over the past 10 Myrs. Such variations serve as powerful probe into the nature and magnitude of plate boundary forces, because they are unlikely to originate from changes in mantle buoyancy forces, which evolve on longer time scales. Here we explore the constraints of the velocity record using a novel coupled modeling-approach of global neo-tectonic simulations combined with realistic plate driving forces obtained from mantle circulation models (MCMs) to arrive at simple global budgets of mantle, lithosphere and plate boundary forces. We focus on three plate boundary systems along the Nazca/South America plate margin, the Aleutian trench and the India/Australia plate boundary to show that gravitational spreading from high topography in the Andes and Tibet contributes substantially to the global plate tectonic force balance and that this contribution is sufficient to explain some 35% of recent velocity changes over the Earth's surface, including among others the observed 30% convergence reduction between the Nazca/South America plates. Our models make a number of specific predictions such as significant lateral variations in plate coupling forces along a given margin revealed by trench-parallel gravity and bathymetry anomalies and the occurrence of large earthquakes, as well as differences by as much as a factor of five from margin to margin. They also support the notion of a relatively young plate boundary separating the India and Australia plates, which has been previously suggested based on independent observations. Importantly, we find that the modeled Nazca/South America convergence reduction explains recent spreading-rate variations in the South Atlantic and South Pacific, which points to the importance of far field effects on the adjacent continents in explaining the spreading record of oceanic basins. Our numerical results demonstrate (a) that detailed budgets of forces acting upon plates can be obtained and (b) support the notion of strong forcing along weak plate boundaries.     

Kinematic model for out-of-sequence thrusting: Motion of two ramp-flat faults and the production of upper plate duplex systems

Kinematic models developed here suggest a bewildering array of structural styles can be generated during out-of-sequence thrusting. Many of these structures would be difficult to distinguish from a normally stacked thrust sequence and the process can produce younger-on-older faults that could easily be misinterpreted as normal faults. This paper considers a small subset of this problem within a large model space by considering structures that develop along a pair of ramp-flat faults that are moving simultaneously, or sequentially. Motion on the lower ramp warps the structurally higher fault due to fault-bend folding and when the fault ruptures through the warp it transfers a horse to the upper hanging wall. Continuity of the process generates what is referred to here as an “upper plate duplex” to distinguish the structure from a conventional duplex. Kinematic parameters are developed for two models within this general problem: 1) a system with a fixed ramp in the lower thrust, overridden by an upper thrust; and 2) a double-duplex system where a conventional duplex develops along the lower fault at the same time as an upper plate duplex is formed along the upper fault. The theory is tested with forward models using 2D Move software and these tests indicate different families of structural styles form in association with relative scaling of ramp systems, slip-ratio between faults, and aspect ratios of horse blocks formed in the upper-plate duplex. A first-order result of the analysis is that an upper plate duplex can be virtually indistinguishable from a conventional duplex unless the trailing branch lines of the horses are exposed or imaged; a condition seldom met in natural exposures. Restoration of an upper-plate duplex produces counterintuitive fault geometry in the restored state, and thus, restorations of upper plate duplexes that erroneously assume a conventional duplex model would produce restored states that are seriously in error. In addition, in most of the models some fault segments place younger rocks on older rocks which could be easily misinterpreted as normal fault systems. In some models younger-on-older juxtapositions are significant and if scaled to crustal scale would produce core-complex style structures that would be difficult to recognize as contractional features. Collectively, these observations imply that many areas where simultaneous contraction and extension are inferred may be entirely contractional with younger-on-older relationships generated by out-of-sequence thrust systems. Examples where this process may have occurred are in southwestern North America and the Moine thrust system and future studies should evaluate these systems in light of these models. Distinguishing upper plate duplex from conventional duplex is potentially important in economic evaluations of thrust systems because fluid migration paths would be very different in the two alternatives. The process may also be important in seismogenic mechanisms, particularly in subduction megathrusts, because faults warping faults could produce fault irregularities that would form transient asperities along the fault.     

A note on the relationship between dip-slip fault parameters and plate parameters in the two-dimensional modelling of preseismic deformation

Preseismic lithospheric deformation at a subduction zone can be modelled as dip-slip dislocation on an inclined fault or as flexure of a thin plate. Both these models predict a region of positive topography known as forebulge or outer rise. By matching the location and the magnitude of the forebulge, we derive useful relations between the dip-slip fault parameters and the plate parameters. In particular, we determine the width of a long dip-slip fault of given dip corresponding to a semi-infinite plate of given thickness. The displacement profiles of the two models are also compared.     

http://www.sciencedirect.com/science/article/pii/S167498711400108X

We present three 3D numerical models of deep subduction where buoyant material from an oceanic plateau and a plume interact with the overriding plate to assess the influence on subduction dynamics,trench geometry,and mechanisms for plateau accretion and continental growth.Transient instabilities of the convergent margin are produced,resulting in:contorted trench geometry;trench migration parallel with the plate margin;folding of the subducting slab and orocline development at the convergent margin;and transfer of the plateau to the overriding plate.The presence of plume material beneath the oceanic plateau causes flat subduction above the plume,resulting in a "bowed" shaped subducting slab.In plateau-only models,plateau accretion at the edge of the overriding plate results in trench migration around the edge of the plateau before subduction is re-established directly behind the trailing edge of the plateau.The plateau shortens and some plateau material subducts.The presence of buoyant plume material beneath the oceanic plateau has a profound influence on the behaviour of the convergent margin.In the plateau + plume model,plateau accretion causes rapid trench advance.Plate convergence is accommodated by shearing at the base of the plateau and shortening in the overriding plate.The trench migrates around the edge of the plateau and subduction is re-established well behind the trailing edge of the plateau,effectively embedding the plateau into the overriding plate.A slab window forms beneath the accreted plateau and plume material is transferred from the subducting plate to the overriding plate through the window.In all of the models,the subduction zone maintains a relatively stable configuration away from the buoyancy anomalies within the downgoing plate.The models provide a dynamic context for plateau and plume accretion in Phanerozoic accretionary orogenic systems such as the East China Orogen and the Central Asian Orogen(Altiads),which are characterised by accreted ophiolite complexes with diverse geochemical affinities,and a protracted evolution of accretion of exotic terranes including oceanic plateau and terranes with plume origins.     

A full-plate global reconstruction of the Neoproterozoic

Neoproterozoic tectonic geography was dominated by the formation of the supercontinent Rodinia, its break-up and the subsequent amalgamation of Gondwana. The Neoproterozoic was a tumultuous time of Earth history, with large climatic variations, the emergence of complex life and a series of continent-building orogenies of a scale not repeated until the Cenozoic. Here we synthesise available geological and palaeomagnetic data and build the first full-plate, topological model of the Neoproterozoic that maps the evolution of the tectonic plate configurations during this time. Topological models trace evolving plate boundaries and facilitate the evaluation of “plate tectonic rules” such as subduction zone migration through time when building plate models. There is a rich history of subduction zone proxies preserved in the Neoproterozoic geological record, providing good evidence for the existence of continent-margin and intra-oceanic subduction zones through time. These are preserved either as volcanic arc protoliths accreted in continent-continent, or continent-arc collisions, or as the detritus of these volcanic arcs preserved in successor basins. Despite this, we find that the model presented here still predicts less subduction (ca. 90%) than on the modern earth, suggesting that we have produced a conservative model and are likely underestimating the amount of subduction, either due to a simplification of tectonically complex areas, or because of the absence of preservation in the geological record (e.g. ocean-ocean convergence). Furthermore, the reconstruction of plate boundary geometries provides constraints for global-scale earth system parameters, such as the role of volcanism or ridge production on the planet's icehouse climatic excursion during the Cryogenian. Besides modelling plate boundaries, our model presents some notable departures from previous Rodinia models. We omit India and South China from Rodinia completely, due to long-lived subduction preserved on margins of India and conflicting palaeomagnetic data for the Cryogenian, such that these two cratons act as ‘lonely wanderers’ for much of the Neoproterozoic. We also introduce a Tonian-Cryogenian aged rotation of the Congo-São Francisco Craton relative to Rodinia to better fit palaeomagnetic data and account for thick passive margin sediments along its southern margin during the Tonian. The GPlates files of the model are released to the public and it is our expectation that this model can act as a foundation for future model refinements, the testing of alternative models, as well as providing constraints for both geodynamic and palaeoclimate models.     

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.     

Tectonic evolution of the triple junction off central Honshu for the past 1 million years

We discuss several models of the evolution of the trench-trench-trench triple junction off central Honshu during the past 1 m.y. on the basis of plate kinematics, morphology, gravity and seismic reflection profile data available for the area. The study area is characterized by large basins, 7–8 km deep on the inner lower trench slope on the Philippine Sea side and the deep (9 km) Izu-Bonin Trench to the east. Between the basins and the trench, there are 6–7 km-deep basement highs. The triple junction is unstable due to the movement of the Philippine Sea plate at a velocity of 3 cm/yr in WNW direction with respect to Eurasia (Northeast Japan), subparallel to the strike of the Sagami Trough. Generally we can expect the boundary area between the Philippine Sea and Pacific plates to be extended because the Pacific plate is unlikely to follow the retreating Philippine Sea plate due to the obstruction of the southeastern corner of Eurasia. The above peculiar morphology of the junction area could have resulted from this lack of stability. However, there are several possible ways to explain the above morphology.

Our gravity model across the trench-basement high-basin area shows that the basement highs are made of low-density materials (1.8–2 g/cm³). Thus we reject the mantle diapir model which proposes that the basement highs have been formed by diapiric injection of serpentinites between the retreating Philippine Sea plate and the Pacific plate.

The stretched basin model proposes that the basins have been formed by stretching of the Philippine Sea plate wedge. We estimated the extension to be about 10 km at the largest basin. We reconstructed the morphology at 1 Ma by moving the Philippine Sea plate 20 km farther to the east after closing the basins, and thus obtained 8 km depth of the 1 Ma trench, which is similar to that of the present Japan Trench to the north. Although this stretched basin model can explain the formation of the basins and the deep trench, other models are equally possible. For instance, the eduction model explains the origin of the basin by the eduction of the Philippine Sea basement from beneath the basement high, while the accretion model explains the basement highs by the accretion of the Izu-Bonin trench wedge sediments. In both of these models we can reconstruct the 1 Ma trench depth as about 8 km, similar to that of the stretched basin model.

The deformation of the basement of the basins constitutes the best criterion to differentiate between these models. The multi-channel seismic reflection profiles show that the basement of the largest basin is cut by normal faults, in particular at its eastern edge. This suggests that the stretched basin model is most likely. However, the upper part of the sediments shows that the basement high to the east has been recently uplifted. This uplift is probably due to the recent (0.5 Ma) start of accretion of the trench wedge sediments beneath this basement high.     

Multibeam bathymetry and sidescan imaging of the Rivera Transform–Moctezuma Spreading Segment junction, northern East Pacific Rise: New constraints on Rivera–Pacific relative plate motion

To better understand the recent motion of the Pacific plate relative to the Rivera plate and to better define the limitations of the existing Rivera–Pacific plate motion models for accurately predicting this motion, total-field magnetic data, multibeam bathymetric data and sidescan sonar images were collected during the BART and FAMEX campaigns of the N/O L'Atalante conducted in April and May 2002 in the area surrounding the Moctezuma Spreading Segment of the East Pacific Rise, located offshore of Manzanillo, Mexico, at 106°16′W, between 17.8°N and 18.5°N. Among the main results are: (1) the principle transform displacement zone of the Rivera Transform is narrow and well defined east of 107o15′W and these azimuths should be used preferentially when deriving new plate motion models, and (2) spreading rates along the Moctezuma Spreading Segment should not be used in plate motion studies as either seafloor spreading has been accommodated at more than one location since the initiation of seafloor spreading in the area of the Moctezuma Spreading Segment, or this spreading center is not a Rivera–Pacific plate boundary as has been previously assumed. Comparison of observed transform azimuths with those predicted by the best-fit poles of six previous models of Rivera–Pacific relative motion indicate that, in the study area, a significant systematic bias is present in the predictions of Rivera–Pacific motion. Although the exact source of this bias remains unclear, this bias indicates the need to derive a new Rivera–Pacific relative plate motion model.     

Eocene to Miocene geometry of the West Antarctic Rift System

Tectonic models for the Late Cretaceous/Tertiary evolution of the West Antarctic Rift System range from hundreds of kilometres of extension to negligible strike-slip displacement and are based on a variety of observations, as well as kinematic and geodynamic models. Most data constraining these models originate from the Ross Sea/Adare Trough area and the Transantarctic Mountains. We use a new Antarctic continental crustal-thinning grid, combined with a revised plate-kinematic model based on East Antarctic – Australia – Pacific – West Antarctic plate circuit closure, to trace the geometry and extensional style of the Eocene – Oligocene West Antarctic Rift from the Ross Sea to the South Shetland Trench. The combined data suggest that from chron 21 (48 Ma) to chron 8 (26 Ma), the West Antarctic Rift System was characterised by extension in the west to dextral strike-slip in the east, where it was connected to the Pacific – Phoenix – East Antarctic triple junction via the Byrd Subglacial Basin and the Bentley Subglacial Trench, interpreted as pullapart basins. Seismic-reflection profiles crossing the De Gerlache Gravity Anomaly, a tectonic scar from a former spreading ridge jump in the Bellingshausen Sea, suggest Late Tertiary reactivation in a dextral strike-slip mode. This is supported by seismic-reflection profiles crossing the De Gerlache Gravity Anomaly in the Bellingshausen Sea, which show incised narrow sediment troughs and vertical faults indicating strike-slip movement along a north – south direction. Using pre-48 Ma plate circuit closure, we test the hypothesis that the Lord Howe Rise was attached to the Pacific Plate during the opening of the Tasman Sea. We show that this plate geometry may be plausible at least between 74 and 48 Ma, but further work especially on Australian – Antarctic relative plate motions is required to test this hypothesis.     

The effect of sampling scale on actualistic sandstone petrofacies

Empirical correlations between plate tectonic setting and sand/sandstone composition have been the basis for large scale petrological models. These models do not explicitly treat sampling scale. Four areas from the western USA with diverse tectonic settings and rock types provide a natural laboratory for sampling sand at three different scales: talus piles to small drainages (first order), streams and rivers draining mountain ranges (second order), and large rivers and marine environments (third order). Existing plate tectonic petrofacies models should only be applied to third order settings because the data were derived from studies of such settings. This is especially true in tectonic settings with diverse source rocks (e.g. continental rifts and transform settings). On the other hand, some settings, such as active magmatic arcs and foreland fold-thrust belts, provide uniform results at any sampling scale because of homogeneity of source rocks. The Rio Grande drainage area is especially complex, with diverse igneous, metamorphic and sedimentary source areas. Some components (e.g. basalt) are destroyed with minimal transport, whereas others (e.g. quartz) are relatively enriched with greater transport. In this complex continental rift setting, first and second order sand is diverse and heterogeneous due to input from tributaries. The Santa Clara River of southern California also has heterogeneous sand due to diverse source rocks in this transform setting. It is only after considerable homogenization and stabilization due to weathering and mixing with more stable components, and/or considerable transport, that homogeneous compositions are produced in these two settings. In contrast, the Cascade magmatic arc and the Canadian Rocky Mountain fold-thrust belt have uniform source rocks (dominantly volcanic in the former and dominantly sedimentary in the latter). Uniform sand composition that is unique to each of these tectonic settings results at any sampling scale in these two cases. Uniformity of data collection and analysis is essential for reproducible results. Use of the Gazzi-Dickinson point counting method allows direct comparison among source rocks (zero order samples), modern sand of any order and ancient sandstone of unknown provenance. Lack of recognition of the effect of sampling scale in the development of actualistic petrofacies models has led to incorrect rejection of many existing models. Third order sands are excellent predictors of plate tectonic setting, but first and second order sands can provide ambiguous plate tectonic interpretations in many settings. More complex actualistic petrofacies models based on diverse sampling scales are needed.     

大陆俯冲带热结构的数值模拟

俯冲带热结构是控制俯冲板块演化的最主要因素之一.前人通过建立解析模型和数值模型对大洋俯冲带热结构进行了一系列研究,发现俯冲板块年龄和俯冲速度是影响俯冲带热结构的关键因素.为了认识大陆俯冲带热结构,特别是理解数值模型结果与岩石学结果之间的差异,我们建立了二维大陆俯冲带运动学和动力学数值模型研究其热结构演化.模型结果显示,如果大陆俯冲板块的俯冲速度与角度和大洋板块一致的话,较低的大陆俯冲带初始温度导致其板块温度比大洋俯冲带低.但是,当大陆俯冲板块的初始温度较高,俯冲速度超慢并且考虑大陆地壳中的放射性元素生热时,模型得到的大陆俯冲带热结构能够解释通过高压和超高压变质岩得到的较热的俯冲板块温度.另一方面,如果俯冲板块与上覆板块存在动力学解耦作用,也能够得到较热的俯冲温压数据.      

http://www.sciencedirect.com/science/article/pii/S1674987114001625

Rock masses are commonly used as the underlying layer of important structures such as bridges, dams and transportation constructions. The success of a foundation design for such structures mainly depends on the accuracy of estimating the bearing capacity of rock beneath them. Several traditional numerical approaches are proposed for the estimation of the bearing capacity of foundations resting on rock masses to avoid performing elaborate and expensive experimental studies. Despite this fact, there still exists a serious need to develop more robust predictive models. This paper proposes new nonlinear prediction models for the ultimate bearing capacity of shallow foundations resting on non-fractured rock masses using a novel evolutionary computational approach, called linear genetic programming. A comprehensive set of rock socket, centrifuge rock socket, plate load and large-scaled footing load test results is used to develop the models. In order to verify the validity of the models, the sensitivity analysis is conducted and discussed. The results indicate that the proposed models accurately characterize the bearing capacity of shallow foundations. The correlation coefficients between the experimental and predicted bearing capacity values are equal to 0.95 and 0.96 for the best LGP models. Moreover, the derived models reach a notably better prediction performance than the traditional equations.     

Angular velocity of Arabian plate from multi-year analysis of GNSS data

We compute the secular tectonic motion of the Arabian plate, based on two different sets of Global Navigation Satellite System (GNSS) observations. The first set is formed by continuous GNSS (cGNSS) observations from 10 stations that are being operated since 2000 and have a data time span of at least 3.5 years. The cGNSS is supplemented with a second set of GNSS observations from 11 episodic stations (eGNSS) acquired between 2003 and 2009 (having at least six reoccupations). In order to evaluate the robustness of the solutions, we have estimated three different solutions using different parameterizations concerning the error models as follows: (a) estimating the velocities based only on the cGNSS sites; (b) combining the cGNSS and eGNSS, where the errors for the eGNSS solutions were rescaled to be consistent with those estimated for the cGNSS solutions; and (c) giving the same weight to all cGNSS and eGNSS solutions. Our conclusions show that case (b) is the preferred solution when the uncertainties of the eGNSS solutions are multiplied by a scale factor (4.1× in this particular case), since the inclusion of the eGNSS solutions allows to have a larger number of points with a better spatial distribution. Finally, we compare our angular velocity model with other models, in particular the most recent global models: GEODVEL (geodetic) and MORVEL (geophysical). We find that GEODVEL provides inaccurate predictions for the Arabian tectonic plate, whereas the MORVEL model shows a slightly better agreement with our solution. Our proposed model shows a significant improvement to model the present-day kinematics of Arabian plate. Thus, it is suggested to be used instead of the global models. As corollary, we also advocate that global models for present-day tectonic motions should be carefully verified for each modeled plate, so that dedicated regional models can provide better results.     

多年冻土区路基路面变形及应力的数值分析

针对青藏公路路基下发育多年冻土融化盘的实际情况,选择两种模型,应用ABAQUS有限元分析软件,对冻土路基从修筑到开放交通过程中的路基路面位移及应力进行了分析.结果表明:冻土路基以融沉为主的变形,一般情况下以路中心下最大,变形呈凹形;当路基下融化盘偏移时,最大变形位置随之偏移;路面层底拉应力最大,对融沉变形反映敏感;路面顶部压应力最大值出现在轴载作用位置,面层应力对轴载反映敏感.计算模型断面尺寸、路基填料、路面结构等对青藏公路具有代表性,在3.6 m路基总高度条件下,无论路基下融化盘偏移与否,融化盘厚达0.5 m时路基顶部(路面层底)拉应力即达基层抗拉强度,显示路基融沉变形可能导致路基失稳及路面破坏,此时路基高度即达最大值.     

Refining accretionary orogen models for the Tasmanides of eastern Australia

The well-known southwest-to-northeast younging of stratigraphy over a present-day cross strike distance of >1500 km in the southern Tasmanides of eastern Australia has been used to argue for models of accretionary orogenesis behind a continually eastwards-rolling paleo-Pacific plate. However, these accretionary models need modification, since the oldest (ca 530 Ma) outcrops of Cambrian supra-subduction zone rocks occur in the outboard New England Orogen, now ～900 km east of the next oldest (520–510 Ma) supra-subduction zone rocks. This is not consistent with simple, continuous easterly rollback. Instead, the southern Tasmanides contain an early history characterised by a westwards-migrating margin between ca 530 and ca 520 Ma, followed by rapid eastwards rollback of the paleo-Pacific plate from 520 to 502 Ma that opened a vast backarc basin ～2000 km across that has never been closed. From the Ordovician through to the end of the Carboniferous, the almost vertical stacking of continental margin arcs (within a hundred kilometres of each other) in the New England Orogen indicates a constant west-dipping plate boundary in a Gondwana reference frame. Although the actual position of the boundary is inferred to have undergone contraction-related advances and extension-related retreats, these movements are estimated to be ～250 km or less. Rollback in the early Permian was never completely reversed, so that late Permian–Triassic to Cretaceous arcs lie farther east, in the very eastern part of eastern Australia, with rifted fragments occurring in the Lord Howe Rise and in New Zealand. The northern Tasmanides are even more anomalous, since they missed out on the middle Cambrian plate boundary retreat seen in the south. As a result, their Cambrian-to-Devonian history is concentrated in a ～300 km wide strip immediately west of Precambrian cratonic Australia and above Precambrian basement. The presence in this narrow region of Ordovician to Carboniferous continental margin arcs and backarc basins also implies a virtually stationary plate boundary in a Gondwana frame of reference. This bipolar character of the Tasmanides suggests the presence of a segmented paleo-Pacific Plate, with major transform faults propagating into the Tasmanides as tear faults that were favourably oriented for the formation of local supra-subduction zone systems and for subsequent intraplate north–south shortening. In this interpretation of the Tasmanides, Lower–Middle Ordovician quartz-rich turbidites accumulated as submarine fan sequences, and do not represent multiple subduction complexes developed above subduction zones lying behind the plate boundary. Indeed, the Tasmanides are characterised by the general absence of material accreted from the paleo-Pacific plate and by the dominance of craton-derived, recycled sedimentary rocks.     

非对称多层薄板折算力学模量的讨论

为使有限单元法中的板壳单元更适合于分析地学问题,本文放弃前人对多层薄板的几何和弹性对称性限制,根据等效原则,重新推导了非对称多层薄板的中性面位置,进而给出了用各分层弹性模量和厚度表示的折合抗弯刚度及弹性模量。并结合两个三层地壳模型例题讨论了地壳各分层弯曲和薄膜应力与等效应力的差异.      

Exhumation mechanisms of melt‐bearing ultrahigh pressure crustal rocks during collision of spontaneously moving plates

A series of 2D petrological–thermomechanical numerical experiments was conducted to: (i) characterize the variability of exhumation mechanisms of ultrahigh pressure metamorphic (UHPM) rocks during collision of spontaneously moving plates and (ii) study the possible geodynamic effects of melting at ultrahigh pressure conditions for the exhumation of high‐temperature–ultrahigh pressure metamorphic (HT–UHPM) rocks. To this end, the models include fluid‐ and melt‐induced weakening of rocks. Five distinct modes of exhumation of (U)HPM rocks associated with changes in several parameters in the models of plate collision and continent subduction are identified as follows: vertical crustal extrusion, large‐scale crustal stacking, shallow crustal delamination, trans‐lithospheric diapirism, and channel flow. The variation in exhumation mechanisms for (U)HPM rocks in numerical models of collision driven by spontaneously moving plates contrasts with the domination of the channel flow mode of exhumation in a majority of the published results from numerical models of collision that used a prescribed plate convergence velocity and/or did not include fluid‐ and melt‐induced weakening of rocks. This difference in the range of exhumation mechanisms suggests that the prescribed convergence velocity condition and the neglect of fluid‐ and melt‐related weakening effects in the earlier models may inhibit development of several important collisional processes found in our experiments, such as slab breakoff, vertical crustal extrusion, large‐scale stacking, shallow crustal delamination and relamination, and eduction of the continental plate. Consequently, the significance of channel flow for the exhumation of UHPM rocks may have been overstated based on the results of the earlier numerical experiments. In addition, the results from this study extend over a larger proportion of the high‐temperature range of P–T conditions documented from UHPM rocks, including those retrieved from HT–UHPM rocks, than the results of experiments from previous numerical models. In particular, the highest peak metamorphic temperatures (up to 1000 °C) are recorded in the case of the vertical crustal extrusion model in which subducted continental crust is subjected to a period of prolonged heating by asthenospheric mantle abutting the continental side of the vertically hanging slab. Nonetheless, some extreme temperature conditions which have been suggested for the Kokchetav and Bohemian massifs, perhaps up to 1100–1200 °C, are still to be achieved in experiments using numerical models.