A New Model of Lithosphere Deformation Beneath the Okinawa Trough Based on Gravity Data

The Ryukyu trench-arc system can be divided into two types according to its subduction model. The normal subduction in the northern part of the Philippine Sea plate creates a hinge sedimentary wedge with large deformation at the collision front, while the oblique subduction in the southern part gives rise to a smaller accretion with small deformation than that in the northern part. The mechanisms that cause the distinction between these two types have been analysed and calculated by using gravity data based on the lithosphere rheology and the stress state of the lithosphere in the subduction boundary. The two types of subduction model are associated with the internal extension in the southern Okinawa Trough and the small extension in the northern part. The difference of the stress state between the two types of subduction model is also manifested in other tectonic features, such as topography, volcanic activity and crust movement. Modeling bathymetric and gravity data from this area suggests that the oblique subduction of low angle, together with smooth geometry of the overlying plate crust, results in small stress released on the south of the trench by the subduction plate. The intraplate faults in the southern Okinawa Trough behind the trench stand in surplus intensive stress. On the other hand, the normal subduction of high angle, together with strong undulation geometry of the overlying crust, results in more intensive stress released in the northern Ryukyu Trench than that in the south. The intraplate faults in the northern Okinawa Trough behind the northern Ryukyu Trench stand in small stress.     

Nodes and Field of Tourist Origins to Ancient Village——A Case Study of Huangcheng Village in Shanxi Province of China

Based on Hagget's theory of spatial structure,researches on the nodes and field of tourist origins to Huang-cheng Village in Shanxi Province of China have been explored.Nodal hierarchy structure of tourist origins is analyzed with cluster analysis and the gravity model.And field of tourist origins is analyzed with attraction radius index(R)and geography concentration index(G).In the field analysis,R and G of Huangcheng Village are compared with Xidi Village that is a world heritage located in Huangshan City of Anhui Province in China.According to comparison of loca-tions of two areas,influential factors for field area of Huangcheng Village are identified.It is concluded that:1)cluster analysis and gravity model can be complementary methods to each other for nodal hierarchy structure analysis of tour-ist origins;and 2)as far as location is concerned,the weak intensity effect of tourism resources in the tourist region is a major cause for explaining why tourist origins to Huangcheng Village are mainly its neighboring areas.Moreover,it is suggested that the regional effect of tourist resources should be regarded as a component of destination attractiveness when applying gravity model.     

Intercontinental response to variations in the Arabia-Eurasia Plate convergence,calcite e-twin evidence of the Kuhbanan Fault system,Central Iran

Variations in the plate convergence direction have generally reflected on the kinematics of the major fault zones developed in the intercontinental parts of the collision zones. The Kuhbanan Fault system is one of the most important intercontinental faults in the Arabia-Eurasia collision zone with a dextral strike slip mechanism. This fault system is composed of three fault strands including Kuhbanan, Bazargan, and Kerman Faults. Here we used calcite e-twin analysis of the vein samples developed in these fault zones to reconstruct deformation condition and the paleodifferential stress direction and magnitude at the Kuhbanan Fault system. Our results represent 190°C-200°C of the deformation temperature and related 5–6 ± 1 km of deformation depth at the Kuhbanan Fault system. Calculated paleodifferential stress magnitude in the Kuhbanan Fault system using e-twin parameters ranges from 169-196 MPa similar to the inner parts of the orogenic systems. Comparing the data set of the Kuhbanan Fault system with previous studies at the Zagros orogen represents an increase of deformation depth and paleodifferential stress magnitudes from the foreland of the Zagros orogen to its hinterland and finally northward around the Kuhbanan Fault system in the central Iranian blocks. We have also proposed a shift of the stress direction from paleo NE to recent N directed by comparing paleostress direction deduced from the calcite e-twin and recent stress direction calculated from GPS and earthquake data analysis around the Kuhbanan Fault system. This kinematics change due to plate reorganization is in agreement with the observed regional variation in convergence direction all in the Arabia-Eurasia collision zone.     

FEM simulation of formation of metamorphic core complex with ANSYS software

This study utilizes ANSYS to establish FEM＇s model of metamorphic core complex, and used thermal-structure analysis to simulate metamorphic core complex＇s temperature field and stress field. The meta- morphic core complex formation mechanism is discussed. The simulation results show that the temperature field change appearing as the earth surface＇ s temperature is the lowest, and the temperature of metamorphic core complex＇ s nucleus is the highest. The temperature field is higher along with depth increase, and the stress field change appearing as the biggest stress occurs in the nucleus. The next stress field occurs at the top of the cover.     

Cohesive strength and seismogenic stress pattern along the active basement faults of the Precordillera-Sierras Pampeanas ranges, western Argentina: An experimental analysis by means of numerical model

A two-dimensional finite element method （FEM） model that incorporates faults, elastic rock physical properties, topographical load due to gravity and far-field plate velocity boundary conditions was used to recognize the seismogenic stress state along the fold-and-thrust belt of the Precordillera-Sierras Pampeanas ranges of western Argentina. A plane strain model with nine experiments was presented here to examine the fault strength with two major rock phyical properties： cohesion and angle of internal friction. Mohr-Coulomb failure criterion with bulk rock properties were applied to analyse faults. The stress field at any point of the model was assumed to be comprised of gravitational and tectonic components. The analysis was focused to recognize the seismogenic shear strain concentrated in the internal-cristaline domain of the orogene shown by the modeling. Modeling results are presented in terms of four parameters, i. e., （i） distributions, orientations, and magnitudes of principal stresses （σ1 and σ3）, （ii） displacement vector1 （iii） strain distribution, and （iv） maximum shear stress （τmax） contour line within the model. The simulation results show that the compressive stress is distributed in and around the fault systems. The overall orientation of of σ1 is in horizontal directions, although reorientations do occur within some stress weaker parts, especially subsequent to the faults. A large-scale shear stress is accumulating along the active faults of Tapias-Villicum Fault （TVF）, Salinas-Berros Fault （SBF）, Ampacama-Niquizanga Fault （ANF） and Las Charas Fault （CF）, which could act as local stress and strain modulators to localize the earthquakesoccurrence.     

Development of a High-Resolution Coastal Circulation Model for the Ocean Observatory in Lunenburg Bay

An advanced ocean observatory has been established in Lunenburg Bay of Nova Scotia, Canada as part of an interdisciplinary research project of marine environmental prediction. The development of a high-resolution coastal circulation model is one of important components of the observatory. The model horizontal resolution is 60m and the vertical resolution is about 1 m. The coastal circulation model is used to simulate the semi-diurnal tidal circulation and associated nonlinear dynamics with the M2 forcing specified at the model open boundaries. The model is also used to simulate the storm-induced circulation in the bay during Hurricane Juan in September 2003, with the model forcing to be the combination of tides and remotely generated waves specified at the model open boundaries and wind stress applied at the sea surface. The model results demonstrate strong interactions between the local wind stress, tidal forcing, and remotely generated waves during this period. Comparison of model results with the surface elevation and current observations demonstrates that the coastal circulation model has reasonable skills in simulating the tidal and storm-induced circulation in the bay.     

THE DIVERGENCE-CONVERGENCE EFFECT OF SURFACE CURRENT AND ITS ACTION ON THE BOHAI SEA

The action of the wind field and the influence of topography can cause divergence or convergence of surface current. The existence of the divergence-convergence effect is proved and the dynamical significance of the divergent or convergent state and its link with many marine phenomena are pointed out. Divergence fields of surface current in the Bohai Sea in winter and summer are obtained by numerical modelling describing the divergence-convergence character of seasonally wind-driven current. The relation between the effect and seasonal marine phenomena is discussed. Study on the divergence-convergence effect of surface current (DCESC)can be an indirect method for testing the calculated results.     

Study on dynamic numerical simulation of mountain system in tibet plateau

Tibet Plateau in mountain system is becoming one of the focuses of global technique research, because its crust is marvelous thick, which is twice of the normal thickness of mountain system in lithosphere, and its rapid raise from the Quaternary.By using a finite element analysis software ABAQUS,the numeric analysis has been carried out and presented in this paper for the lithosphere stress field.It is the first time to use the displacement loading in the simulation of Tibet Plateau. During the analysis,the deformed elements are used to simulate the structure band, and friction mechanism is used to model the fracture band. The boundary conditions are given according to the boundary displacements around the Plateau. The stress and displacement distributions are obtained for the geological evolution of the plateau, which are consistent with P axial orientations of the seismic origin mechanism and the measures principle stress orientations. The analysis is also given for the dynamic lithosphere evolution of the Mountain System in the Tibet Plateau.     

Joint inversion of magnetotelluric and fullwaveform seismic data based on alternating cross-gradient structural constraints

Magnetotelluric(MT)inversion and seismic inversion are important methods for the interpretation of subsurface exploration data,but separate inversion of MT and seismic produces ambiguous and non-unique results due to various factors.In order to achieve accurate results,the authors propose a joint inversion method of two-dimensional MT and seismic data in the frequency domain.The finite element method is used for numerical simulation of electromagnetic data in the forward modelling,and the Gauss-Newton method is used for the inversion.The 9-point-finite-difference method is used to solve the seismic wave field in the acoustic wave equation,and the inverse problem of seismic data is solved by full waveform inversion with a conjugate gradient,a simple and fast method.Cross gradient functions are used to provide constraint structure between resistivity and velocity parameters to carry out the joint inversion.The joint inversion algorithm is tested by double-rectangular model synthesis data,and the accuracy of the algorithm is verified.The results show that the joint inversion results are better than those from separate inversion.The algorithm is applied to a geophysical model of a metalliferous deposit in Jinchuan and is compared with the separate inversion results.It shows that the results obtained with joint inversion are much closer to the real model.     

Improvement of the one-dimensional Vertical Advection-diffusion Model in Seawater

The classical 1-D vertical advection-diffusion model was improved in this work.The main advantages of the improved model over the previous one are:1) The applicable condition of the 1-D model is made clear in the improved model,in that it is substantively applicable only to a vertical domain on which two end-member water masses are mixing.2) The substitution of parameter f(z) in the equation of the classical 1-D model with end-member fraction f1^z makes the model more precisely and easily solved.3) All the terms in the improved model equation have specific physical meanings,which makes the model easily understood.Practical application of the improved model to predict the vertical profiles of dissovled oxygen and micronutrients in abyssal ocean water of the North Pacific proved that the improvement of the 1-D advection-diffusion model is successful and practicable.     

A new model of lithosphere deformation beneath the Okinawa Trough based on gravity data

The Ryukyu trench-arc system can be divided into two types according to its subduction model. The normal subduction in the northern part of the Philippine Sea plate creates a hinge sedimentary wedge with large deformation at the collision front, while the oblique subduction in the southern part gives rise to a smaller accretion with small deformation than that in the northern part. The mechanisms that cause distinction between these two types have been analysed and calculated by using gravity data based on the lithosphere rheology and the stress state of the lithosphere in the subduction boundary. The two types of subduction model are associated with the internal extension in the southern Okinawa Trough and the small extension in the northern part. The difference of the stress state between the two types of subduction model is also manifested in other tectonic features, such as topography, volcanic activity and crust movement. Modeling bathymetric and gravity data from this area suggests that the oblique subduction of low angle, together with smooth geometry of the overlying plate crust, results in small stress released on the south of the trench by the subduction plate. The intraplate faults in the southern Okinawa Trough behind the trench stand in surplus intensive stress. On the other hand, the normal subduction of high angle, together with strong undulation geometry of the overlying crust, results in more intensive stress released in the northern Ryukyu Trench than that in the south. The intraplate faults in the northern Okinawa Trough behind the northern Ryukyu Trench stand in small stress.     

����GPS���ݷ������ձ���������ɱ���齻�紦��Ӧ����Ƚ�

??????????????????????????GPS???????????(1997-01~1999-12)??????Hori???????????????????齻????????????????????????????:????????????С??????????????????????????????7 cm/a??50%?????????????10~30 km??????????????????齻??????????????Yabuki-Matu??ura????????????????????????????Hori???????????б?????????????????????????????????????????????????в?????????Ч????????????????????????????????Hori?????????????????????????????????????????????????????????????????????????     

利用加密GNSS数据和震源机制解分析川滇块体南部现今地壳活动特性

利用多源项目获得的补充加密GNSS观测资料计算研究区高空间分辨率的GNSS水平运动速度场和应变率场,采用地震波形资料,通过CAP方法求解震源机制解,在此基础上使用阻尼区域应力反演方法分析川滇地块区域构造应力场的空间分布特征。结合主要断裂的活动特性,综合分析川滇地区的地壳形变特征。结果表明：1）块体浅部的最大主压应力与地表的最大主压应变率由SSE向转变为近SN向,呈现出较好的一致性,GNSS观测至少可反映20 km以内的地壳形变;2）川滇块体南部的东向滑移和顺时针旋转,可能受青藏高原推挤、华南块体阻挡及印度板块与欧亚板块之间北向运动速率自西向东递减而形成的右旋剪切拖拽作用的顺时针力偶的综合影响;3）川滇块体在综合力偶作用下SE向挤出和顺时针旋转的同时,受到走滑逆冲断裂带的吸收转换,使得青藏高原物质SE向挤出有限。     

Plate convergence in the Indo-Pacific region

The Indo-Pacific convergence region is the best target to solve the teo remaining challenge s of the plate tectonics theory,i.e.,subduction initiation and the driving force of plate tectonics.Recent studies proposed that the Izu-Bonin subduction initiation belongs to spontaneous initiation,which implies that it started from extension,followed by low angle subduction.Numerical geodynamic modeling suggests that the initiation of plate subduction likely occurred along a transform fault,which put the young spreading ridge in direct contact with old oceanic crust.This,however,does not explain the simultaneous subduction initiation in the west Pacific region in the Cenozoic.Namely,the subduction initiations in the Izu-BoninMariana,the Aleutian,and the Tonga-Kermadec trenches are associated with oceanic crusts of different ages,yet they occurred at roughly the same time,suggesting that they were all triggered by a maj or change in the Pacific plate.Moreover,low angle subduction induces compression rather than extension,which requires external compression forces.Given that the famous Hawaiian-Emperor bending occurred roughly at the same time with the onset of westward subductions in the west Pacific,we propose that these Cenozoic subductions were initiated by the steering of the Pacific plate,which are classified as induced initiation.Induced subduction initiation usually occurs in young ocean basins,forming single-track subduction.The closure s of Neo-Tethys Oceans were likely triggered by plume s in the south,forming northward subductions.Interestingly,the Indian plate kept on moving northward more than 50 Ma after the collision between the Indian and Eurasian continents and the break-off of the subducted oceanic slab attached to it.This strongly suggests that slab pull is not the main driving force of plate tectonics,whereas slab sliding is.     

日本至中国东北区域构造应力场分布与Benioff带的形态关系研究

收集2 291个来源于NIED F-net宽频带地震台网中心M_W≥3.0地震的震源机制解,采用MSATSI构造应力反演方法,得到反映应力相对大小的R值的量θ和3个主应力轴的方位角与倾伏角,同时反演出每个网格的极射赤平投影P、T轴。以3D应力反演为导向,采用矩阵网格分区的方法,以纬度、经度、深度为坐标,将发生在每一区域的震源机制解投射到相应的矩阵区中,并反演出每个矩阵网格点的应力场变化及分布情况。反演结果表明,与板块交接的浅部区域不仅受到太平洋板块的挤压俯冲,同时也受北美板块的斜插作用;东北区域的深震为Benioff带俯冲作用的结果。     

Modern movement and deformation in the South China Sea shown by GPS measurements and numerical simulation

To better understand the crustal deformation of the South China Sea Basin, we produce a mechanically consistent 2-dimensional model for observing regional velocity field in the South China Sea (SCS). We studied the dominating regional tectonic stress field by geodetic measurements and finite element analysis, the spatial variations of velocity field and strain field, and relative movements among different blocks, using a 2-dimensional model describing crustal deformation of the South China Sea Basin. Strain results show that the SCS is extending at present. The western part of SCS is opening gradually in NWSE direction from its northern margin to the south, but the eastern part of SCS is opening gradually from its central part to the north and south. In addition, we analyzed the plate kinematics to the deformation of the SCS, using a two-dimensional finite element model. Our simulations results are well explained by available geodetic data. The movement of SCS is resulted from interactions among Indian Plate, Pacific Plate, Philippine Sea Plate, and Eurasian Plate.     

走滑断裂活动导致地应力解耦的机理研究——以龙门山断裂带东北段为例

2008年汶川MS8.0地震后,沿龙门山断裂带系列水压致裂地应力测量结果表明,断裂带东北段(茂县-绵竹连线东北侧)北川、平武及广元等地浅表层最大主应力方向与区域应力场相比发生了25°~69°不等的偏转,亦称地应力解耦。初步分析认为,龙门山断裂带东北段的右旋走滑活动直接导致了区域内最大主应力方向的解耦,且最大主应力方位向偏转角度与断裂右旋走滑活动量正相关。为合理解释这一现象,探讨走滑断裂活动导致地应力解耦的机理,本文采用二维有限元数值模拟方法,选取汶川MS8.0地震断裂(龙门山中央断裂)上的同震右旋走滑位移作为边界条件,首先模拟了龙门山断裂带东北段右旋走滑活动产生的扰动应力特征,之后结合断裂带附近典型钻孔地应力测量结果,叠加该扰动应力作用,从力学机理上深入剖析断裂右旋走滑活动产生的扰动应力对区域应力场(或初始应力场)的影响,并对其附近不同构造部位最大主应力方向偏转现象进行定量解释。初步结果表明：1) 走滑断裂活动导致地应力解耦主要是通过断裂滑动产生的扰动应力对区域应力场的附加影响实现的;2) 在固定的断裂活动长度和摩擦强度情况下,走滑断裂活动导致区域最大主应力方向偏转角度基本与活动断裂走滑位移量呈正比;3) 受汶川MS8.0地震影响,龙门山断裂带大致以茂县-绵竹为界,西南侧的宝兴-映秀-汶川一带最大主应力方向偏转角度很小,基本与区域应力场一致,而东北侧的北川-平武-广元一带则发生角度不等的偏转,3个地区最大主应力方向平均偏转角度分别为65°、56°和29°,该角度与龙门山断裂带东北段同震右旋走滑位移量呈正相关。     

论双侧造山带——以桐柏—大别造山带为例

从双侧造山带的几何学和力学模式出发,详细地讨论了桐柏—大别造山带的构造特征。在综合地质和地球物理资料基础上,提出一个解释桐柏—大别造山带漂浮体及花状构造的运动学模式。强调指出,具有两个不同变形楔状体的桐柏—大别造山带,是由华北及扬子克拉通间在印支—燕山期斜向碰撞造成的。     

۰®���Ѵ�����Ӧ��״̬�о�����

???????????????????????????????????о???????????????????????????????????????????????????????????????NEE??????????ж???EW??????????????NEE?????????????????????????????????????????????????????????????????о?????????????????????????????????????о????????     

青藏高原东南缘构造应力场重构

采用有限元方法,针对青藏高原东南缘建立更细致、更精确的三维有限元弹性模型。选取9种不同的应力边界条件,分别进行优化分析后处理,将对应台站形变模拟值与GPS实测值进行误差分析,最终选取最佳方案作为古构造应力场。结果表明,青藏高原东南缘4 Ma BP的古应力场主要起源于中国大陆周围板块的相互作用,特别是印度板块NNE向强烈碰撞作用,成为中国大陆尤其是西南部青藏高原地区构造应力场最主要的动力来源,控制各个块体相互作用的方式和运动格局。青藏高原东南缘古应力场主要包括几个力源：西北部青藏高原侧向挤压造成的WE向应力约105 MPa;西南部直接来自于印度板块的NE-WS向应力约70 MPa;南部NS向作用力33 MPa;东南部扬子块体侧向NW-SE阻挡力56 MPa;北东部受扬子块体强烈EW向阻挡力90 MPa。这些力源共同作用于青藏高原东南缘,形成现今复杂应力场。