Current horizontal strain field in Chinese mainland derived from GPS data

Introduction In the years when the reliable data could not be obtained and in the analysis of strain property and magnitude in history, the intensity, property and activity pattern of strain field were mainly inferred on the bases of geometric characters of surface traces and behaviors (especially the faults) as well as the characteristics of petrology (XIE, et al, 1993; Molnar, Tapponnier, 1975, 1977; Tapponnier, Molnar, 1977; FU, et al, 2000). However, they are the averaged results accumu…     

中国大陆及周边地区的水平应变场

推导并建立了块体的两种弹性运动方程：块体的整体旋转与均匀应变方程和块体的整体旋转与线性应变方程. 应用统计学原理,使用西域、青藏和华北块体上的GPS站速度数据,对这两种弹性运动方程与刚体运动方程模拟块体站速度的无偏性和有效性进行了统计检验. 检验结果表明,块体的整体旋转与线性应变方程是描述块体运动的最优模型. 将中国大陆划分为10个块体,应用块体的整体旋转与线性应变方程和块体上的GPS站速度估计了各个块体上的旋转与应变参数,按照1°×1°的间距计算了中国大陆及周边地区上1005个点的应变参数,分析了中国大陆及周边地区应变场的基本特征. 用本文方法得到的主压应变方向与地质学方法和测震学方法得到的主压应力轴方向具有很好的一致性(华南块体除外).     

Recent crustal deformation of İzmir, Western Anatolia and surrounding regions as deduced from repeated GPS measurements and strain field

To investigate contemporary neotectonic deformation in İzmir, Western Anatolia and in its neighborhood, a relatively dense Global Positioning System (GPS) monitoring network was established in 2001. Combination of three spatially dense GPS campaigns in 2001, 2003 and 2004 with temporally dense campaigns between 1992 and 2004 resulted in a combined velocity field representing active deformation rate in the region. We computed horizontal and vertical velocity fields with respect to Earth-centered, Earth-fixed ITRF2000, to Eurasia and to Anatolia as well.The rates of principal and shear strains along with rigid-body rotation rates were derived from velocity field. Results show east–west shortening between Karaburun Peninsula and northern part of İzmir Bay together with the extension of İzmir Bay in accordance with general extension regime of Western Anatolia and Eastern Agea. East–west shortening and north–south extension of Karaburun Peninsula are closely related to right-lateral faulting and a clockwise rotation. There exists a block in the middle of the peninsula with a differential motion at a rate of 3–5 ± 1 mm/year and 5–6 ± 1 mm/year to the east and south, respectively.As is in Western Anatolia, north–south extension is dominant in almost all parts of the region despite the fact that they exhibit significantly higher rates in the middle of the peninsula. Extensional rates along Tuzla Fault lying nearly perpendicular to İzmir Bay and in its west are maximum in the region with an extension rate of 300–500 ± 80–100 nanostrain/year and confirm its active state. Extensional rates in other parts of the region are at level of 50–150 nanostrain/year as expected in the other parts of Western Anatolia.     

鄂尔多斯块体周缘地区现今地壳水平运动与应变

位于青藏块体和华北块体之间的鄂尔多斯块体及其周缘地区是中国大陆构造活动最活跃的地区之一,从1300年至今,在块体周边断陷盆地和西南缘断裂带上发生了五次8级以上的地震.为了了解该地区现今地壳运动、应变状态以及断裂滑动分布,我们收集了中国大陆构造环境监测网络2009-2013年、国家GPS控制网、跨断陷盆地的8个GPS剖面等共527个流动站和32个连续站GPS观测数据,获得了高空间分辨率的地壳水平运动速度场,进一步用均匀弹性模型计算了应变率分布.结果表明,块体内部GPS站点向NEE方向运动,速度变化较小,应变率大多在（-1.0~1.0）×10-8/a之间;山西断陷带构造运动与变形最为强烈,盆地相对于鄂尔多斯块体为拉张变形,应变率为（1.0~3.0）×10-8/a,相对于东部山地则为挤压变形,应变率为（-2.0~-3.0）×10-8/a,盆地西侧断裂（如罗云山断裂、交城断裂）以拉张运动为主,拉张速率为2~3 mm·a-1,盆地东侧断裂主要以右旋缩短运动为主,速率为1~3 mm·a-1;河套断陷带西部的临河凹陷处于较强的张性应变状态,应变率为（2.0~3.0）×10-8/a;块体西南边缘处于压缩应变状态,应变率为（-1.0~-2.0）×10-8/a,六盘山断裂存在明显的地壳缩短运动,速率约为2.1 mm·a-1,速率在断裂附近逐渐减小,反映了断裂处于闭锁状态;相对于鄂尔多斯块体内部渭河断裂带为左旋运动,速率为1.0 mm·a-1,盆地处在弱拉张变形状态.     

海原构造区及其周缘上部地壳结构研究

基于甘肃省夏河县—陕西省靖边县剖面的8次人工地震初至波数据, 利用有限差分走时方法反演得到了沿该剖面长约650 km的上部地壳速度结构和结晶基底的深度分布． 反演结果显示: 海原构造区西侧的西秦岭—祁连山褶皱区上部地壳的横向非均匀性明显, 基底深度从1 km到5 km不等, 反映了褶皱区改造变形强烈的构造特征; 其东侧的鄂尔多斯盆地基底深度约为5—6 km, 其速度均匀、 稳定, 上地壳呈弱速度梯度特征; 海原构造区及海原弧形断裂带附近上部地壳的破坏变形最严重, 区内横向高低速相间分布． 综上可知, 海原构造区东西两侧上地壳结构的显著差异揭示了其结构复杂性的成因及其与地震活动性的关系．      

Solar terminator-related ionosphere derived from GPS TEC measurements: A case study

We present a case study of ionospheric wave activity during solar terminator crossing using GPS TEC measurements. As a basic tool, the spatial gradient of total electron content (TEC) has been used. We tested with positive result the hypothesis on anisotropic response to assumed one-dimensional solar terminator forcing. The approximate drift velocity of irregularities has been computed. The wavelet analysis gave an interesting insight into variable frequency content of TEC gradient time series. We also proposed a filtering with respect to spatial scale allowing for resolving the spatio-temporal ambiguity.     

Field of tectonic stresses from focal mechanisms of earthquakes and recent crustal movements from GPS measurements in China

Orientations of the principal axes of the tectonic stress field reconstructed from seismological data on focal mechanisms of earthquakes and strain fields determined from GPS measurements in China are compared. The data of GPS measurements used in the paper were obtained by the Crustal Movement Observation Network of China (about 1000 stations) in the period of 1998–2004. On the basis of information on the recent horizontal crustal motions, the strain field is calculated for the study territory by the finite element method. Calculations of the strain tensor using GPS data were carried out with a step of 1° in latitude and longitude. A catalog of earthquake focal mechanisms was used for the reconstruction of tectonic stress field components. Focal mechanisms of earthquakes were calculated with the use of seismological data on signs of first arrivals from the bulletin of the International Seismological Center. To estimate characteristics of the regional stress field, an approach based on the kinematic method proposed by O.I. Gushchenko was applied. The tectonic stress field was reconstructed in depth intervals of 0 < H < 35 km and 35 km < H < 70 km from data on focal mechanisms of earthquakes over the periods of 1998–2004 and 1985–2004. Comparison of directions of the principal strain axes at the surface (according to GPS measurements) and directions of the principal stress axes (reconstructed from focal mechanisms of earthquakes) showed their good convergence. Seismotectonic strains and GPS measurements coincide within a larger part of the territory. The coincidence is best in a depth interval of 0 < H < 35 km. Maximum misfit values are confined to areas of high 3-D gradients of strain axis directions and are possibly related to the structural heterogeneity of the region, zones with strains of the same type along both horizontal axes (compression or extension along all directions), or areas of small absolute values of recent horizontal movements. Areas with invariable directions of the stress axes are recognizable regardless of the depth of initial data. Good reproducibility of results obtained by two different methods made it possible to check the method of stress field reconstruction using data on focal mechanisms of earthquakes.     

A precise velocity field of tectonic deformation in China as inferred from intensive GPS observations

##正##Active tectonics, e.g., faulting, folding, and rifting are clearly manifested over the vast territory of China. The in-tensive crustal deformation with long-lived structures has given rise to numerous highly-elevated mountains in the western China with the Yangtze and Yellow rivers downstreaming to the east,     

Ambient tectonic shear stress field in Southern California and seismic hazard regions

According to the fracture mechanics rupture model of earthquakes put forward by us, several equations to compute tectonic ambient shear stress value τ₀ have been derived [equations (1), (2), (3), (5)]. τ₀ values for intermediate and small earthquakes occurred in Chinese mainland and Southern California have been calculated by use of these equations. The results demonstrate that the level and distribution of τ₀ are closely related to the location where large earthquakes will occur, i.e. the region with higher level of τ₀ will be prone to occur large earthquakes and the region with lower level will usually occur small earthquakes. According to the spatial distribution of τ₀, the seismic hazard regions or the potential earthquake source regions can in some degree be determined. According to the variation of τ₀ with time, the large earthquake occurrence time can be roughly estimated. According to the distribution of τ₀ in Southern California and variation with time, three high stress level regions are determined, one (Goldfield area) of them is the present seismic hazard region. Contribution No. 98A02023, Institute of Geophysics, State Seismological Bureau, China     

Ambient tectonic shear stress field in Southern California and seismic hazard regions

ＡｍｂｉｅｎｔｔｅｃｔｏｎｉｃｓｈｅａｒｓｔｒｅｓｆｉｅｌｄｉｎＳｏｕｔｈｅｒｎＣａｌｉｆｏｒｎｉａａｎｄｓｅｉｓｍｉｃｈａｚａｒｄｒｅｇｉｏｎｓＰＥＩＳＨＡＮＣＨＥＮ（陈培善）ＬＥＩＸＩＡＯ（肖磊）ＴＯＮＧＸＩＡＢＡＩ（白彤霞）ＪＩＮＣＨ...     

基于Monte Carlo方法的由GPS观测计算地应变率的误差分析

在利用GPS观测计算地应变率中,应变率的误差分布在地球动力学研究中非常重要. 为解决这一问题,文中提出基于Monte Carlo技术,利用GPS观测的速度计算地应变率进行误差分析的方法;以青藏高原地区为例,利用带有误差的GPS观测数据进行了大量、独立的应变率重复计算,然后运用统计方法给出了应变率分量及其误差的分布. 结果显示青藏高原的喜马拉雅及中东部地区地应变率的东西向、南北向及剪应变率分量的误差相对较大,但整个青藏高原及附近地区应变率结果稳定可靠;研究结果进一步表明,由GPS观测计算应变率的误差主要决定于GPS观测速度的误差, 在变形小的区域（如中国东部地区）,目前GPS的观测精度还不能完全满足地应变率计算的要求.     

Preliminary results on kinematic model of tectonic blocks derived from high precision GPS observations in Southwest China

ＰｒｅｌｉｍｉｎａｒｙｒｅｓｕｌｔｓｏｎｋｉｎｅｍａｔｉｃｍｏｄｅｌｏｆｔｅｃｔｏｎｉｃｂｌｏｃｋｓｄｅｒｉｖｅｄｆｒｏｍｈｉｇｈｐｒｅｃｉｓｉｏｎＧＰＳｏｂｓｅｒｖａｔｉｏｎｓｉｎＳｏｕｔｈｗｅｓｔＣｈｉｎａＬＩＲＥＮＨＵＡＮＧ１）（黄立人）...     

中国大陆及邻区现代构造应力场分区

以“中国大陆地壳应力环境基础数据库”资料为基础,总结了中国大陆及邻区现代构造应力场的基本特征,提出了划分构造应力分区的原则,将中国大陆及邻区现代构造应力场划为四级应力区. 指出一级构造应力区主要受板块边界的几何特征和作用在边界上的力所控制,应力作用在大范围和长时间内保持均匀和稳定;二级构造应力区受控于区域块体间的相互作用,在较大区域范围内应力作用具有较强的相关性;三级构造应力区受控于区域内部块体间的相互作用,在一定区域范围内应力结构具有相似性;四级构造应力区受控于块体和断裂相互作用的影响,其应力性状(方向、强度、结构等)一致性较好. 论文还初步分析了控制中国大陆及邻区现代构造应力场分布的动力学环境.     

Shallow Seismic Velocity Structure of the Betic Cordillera (Southern Spain) from Modelling of Rayleigh Wave Dispersion

A detailed dispersion analysis of Rayleigh waves generated by local earthquakes and occasionally by blasts that occurred in southern Spain, was undertaken to obtain the shear-wave velocity structure of the region at shallow depth. Our database includes seismograms generated by 35 seismic events that were recorded by 15 single-component short-period stations from 1990 to 1995. All these events have focal depths less than 10 km and body-wave magnitudes between 3.0 and 4.0, and they were all recorded at distances between 40 and 300 km from the epicentre. We analysed a total of 90 source-station Rayleigh-wave paths. The collected data were processed by standard digital filtering techniques to obtain Rayleigh-wave group-velocity dispersion measurements. The path-averaged group velocities vary from 1.12 to 2.25 km/s within the 1.0-6.0 s period interval. Then, using a stochastic inversion approach we obtained 1-D shear-wave velocity–depth models across the study area, which were resolved to a depth of circa 5 km. The inverted shear-wave velocities range approximately between 1.0 and 3.8 km/s with a standard deviation range of 0.05–0.16 km/s, and show significant variations from region to region. These results were combined to produce 3-D images via volumetric modelling and data visualization. We present images that show different shear velocity patterns for the Betic Cordillera. Looking at the velocity distribution at various depths and at vertical sections, we discuss of the study area in terms of subsurface structure and S-wave velocity distribution (low velocity channels, basement depth, etc.) at very shallow depths (0–5 km). Our results characterize the region sufficiently and lead to a correlation of shear-wave velocity with the different geological units features.     

Circulation in a Southern Italy coastal basin: Modelling and field measurements

The purpose of the present work is to study the hydrodynamic aspects in the Mar Piccolo, a coastal basin located on the northern side of the Gulf of Taranto in the Ionian Sea (Italy), by means of mathematical modelling and field measurements. The latter were assessed during three surveys carried out in the spring–summer of 2002. Collected data have been utilized as input by the 3-D Princeton Ocean Model, which is a sigma coordinate, free surface ocean model which was developed in the late 1970s by Blumberg and Mellor. Simulations in baroclinic condition were forced by a homogeneous and stationary wind field, a simple tidal wave, a constant outflow and vertical stratification of temperature and salinity. A comparison was made between the mathematical modelling results and the field measurements collected during the surveys, in terms of velocity. It was observed that during small tides, when the wind effect prevails over the stratification effect, the best model results were obtained for the most superficial layer and that superficial patterns reproduced by the model are more sensitive to wind direction than to stratification. On the contrary, when the wind effect decreases or the thermohaline effect rises, best results occurred in deeper layers.     

基于实测与模拟的青海湖冰厚时空变化特征

湖冰厚度是湖泊在封冻期的重要物理参数,明晰其时空变化特征对于认识气候变暖背景下的湖冰响应规律具有重要的理论价值和现实意义.基于ERA5 Climate Reanalysis气温数据集、MODIS MOD09GQ数据产品和2019年湖冰钻孔测厚数据及雷达测厚数据,重建20002019年青海湖冰厚时间序列并分析其时空变化特...     

龙门山断裂带及其周边地区重力场 和岩石层力学特性研究

继2008年汶川M_S8.0地震之后, 2013年4月20日又发生了芦山M_S7.0地震, 两次地震的发震构造同属龙门山断裂带. 根据最新的重力和地形资料, 采用岩石层弹性板模型, 计算了龙门山断裂带及其周边地区的二维岩石层有效弹性厚度分布, 并从岩石层的力学特征分析了穿过两次地震震中位置的重力剖面特征; 结合以往在该地区的研究成果, 分析了岩石层的力学变形问题. 结果表明, 以龙门山为界, 四川盆地所在的扬子板块弹性厚度为33(±4) km, 龙门山西北的松潘—甘孜地块的弹性厚度为13(±4) km, 两侧岩石层存在明显的力学强度差异. 包括两次地震震中范围的龙门山断裂带南部区域的有效弹性厚度值小于北部地区, 说明该区域的岩石层更容易发生变形, 可以解释在构造上具备强震发生的岩石层动力学条件.     

河北及周边地区岩石圈磁场变化及震前异常特征

利用2019年度大华北秋季复测156个测点的流动地磁矢量资料,获得河北及周边地区岩石圈磁场半年变化图,分析归纳河北及周边地区岩石圈磁场的总体变化特征,并研究测区内2019年12月3日河北怀安M_S3.4地震和2019年12月5日河北丰南M_S4.5地震前岩石圈磁场局部变化和异常特征。结果表明:丰南M_S4.5地震发生在H矢量的转向区,总强度F负异常的高值区,D要素负异常的高梯度带和Z值变化的高值区;怀安M_S3.4地震则发生在H矢量幅值弱化且转向区,D要素零值线附近的正值区和Z值变化的低值区。