华北地区地震环境噪声特征研究

利用华北流动地震台阵观测的垂直分向连续波形数据, 通过计算功率谱密度和相应的概率密度函数, 对华北地区地震环境噪声特征进行了分析研究. 结果表明, 东部平原和沉积盆地2 Hz以上的高频环境噪声水平与全球新高噪声模型（NHNM）相近, 周期3——18 s的平均噪声水平低于NHNM和新低噪声模型（NLNM）的平均值; 山区及西部高原的高频噪声水平明显低于NHNM, 周期1——18 s的噪声水平大多明显低于NHNM和NLNM的平均值; 不同区域18 s以上周期的噪声水平差异相对较小. 流动地震台阵部分台站的环境噪声存在明显的昼夜变化, 个别台站噪声水平明显高于周边台站, 表明这些台站受人类活动干扰较大. 不同台站的噪声水平分析表明, 将台站布设在摆坑内, 能在一定程度上降低高频和低频段的噪声水平. 台站环境噪声特征的研究结果可为流动地震台阵观测数据质量的定量评估, 观测期间的台站优化调整等提供重要依据.      

Two-dimensional nonlinear site response analysis of Adapazarı plain and predictions inferred from aftershocks of the Kocaeli earthquake of 17 August 1999

The objective of this study is to investigate the effects of local site conditions in the Adapazarı plain crossing the severely damaged central part of Sakarya during the 17th August earthquake. A two-dimensional model has been used to obtain the free-field motions in the valley and on surface formations. A preliminary geotechnical profile model developed from available but limited soil data was checked using recordings of some specific locations where a strong ground motion accelerometer array was in operation after the main earthquake. The range of validity of the model was assessed and modifications were made to compare with the actual recorded motions. The purpose of building such a model, which was used to study the aftershocks in Adapazarı, is to provide an improved database to be used for the design of structures in the city. This study also attempts to provide greater insight into the local site response phenomena through the use of a two-dimensional nonlinear analysis. Simulated site responses are in global agreement with the recorded data. In general agreement between the amplification functions of the computed and recorded data is satisfactory for the frequencies of engineering interest. Alluvial sites show amplification factors in the range of 4–6 in the frequencies between 0.6 and 2 Hz. A relatively shallow alluvial station, HASTAHANE, experienced considerable amplification for small aftershocks and deamplification for the strongest aftershock in frequencies between 0.5 and 2 Hz. Another station, TOYOTA, with significant deamplification characteristics and a clearly observed peak response frequency shift at the soil resonant frequency may have experienced soil shear failure during the strongest aftershock. Other alluvial stations, SEKER and GENC, show deamplification characteristics at 0.55 and at around 1 Hz. with the peak frequency shifts during the magnitude 5.8 aftershock. In general, site responses are larger in the alluvial basin for all aftershocks except the magnitude 5.8 event.     

大兴安岭造山带及两侧邻区莫霍面深度与泊松比分布特征

利用2009～2016年内蒙古自治区数字地震台网宽频带固定地震台站的远震波形数据,采用接收函数H-k算法获得23个基岩台站下方的莫霍面深度和泊松比,同时,收集并筛选出277个已有探测台阵和流动台站的接收函数研究结果,综合分析给出大兴安岭造山带及两侧邻区莫霍面深度、泊松比的分布特征。研究表明,研究区域的莫霍面在整体上呈现自东向西逐渐加深的特征,莫霍面深度为25.0～42.3km,平均约为33.5km。莫霍面最浅的区域为松辽盆地(深度为27.0～35.0km),最深的区域为大兴安岭重力梯级带以西地区(深度为41.0～42.3km)。研究区域泊松比为0.19～0.33,平均值为0.26,大于全球大陆地壳的平均值。泊松比高值异常区集中在火山岩区及具有较厚沉积层的盆地。台站所处位置的海拔与莫霍面深度之间具有较强的正相关性,艾里补偿模式在研究区成立,莫霍面起伏与区域地形地貌特征间具有显著的镜像关系。大兴安岭地区的莫霍面深度与泊松比间存在显著的反相关关系,而在松辽盆地及周缘地区未发现明显的规律性,这也意味着松辽盆地在构造演化过程中经历了更为复杂的地壳改造过程。     

Optimum strong-motion array geometry for source inversion

Optimum strong-motion array geometry for source inversion is quantitatively determined for each of three types of earthquake faults: strike—slip, dip-slip and offshore subduction thrust. The estimation is done by using the source inversion method^1,2 previously developed and based on Wolberg's prediction analysis.³ Only the far-field S waves are taken into account in source inversion; we test only the ability to estimate seismic moments of subfault elements. A comparison of the optimum array geometries obtained in this study with the ones proposed on the basis of an empirical judgment at the 1978 International Workshop on Strong-Motion Earthquake Instrument Arrays⁴ shows that the Workshop geometries are not the best ones as far as source studies are concerned. The most preferable array involves two different kinds of stations regardless of the target fault type: stations close to the fault and aligned parallel to its strike, and stations surrounding the fault area with good azimuthal coverage. The former stations resolve the later stage of the rupturing process while the latter stations resolve the earlier stage.     

The Aterno valley strong-motion array: seismic characterization and determination of subsoil model

The paper focuses on the strong motion array deployed in the upper Aterno River Valley, in the immediate outskirts north-west of the town of L’Aquila, which is part of the Italian Strong Motion Network operated by the Department of Civil Protection. The array is composed of six accelerometric stations located along a cross section of the valley. The importance of this array relies on the fact that a large amount of high-quality records were obtained during the 2009 L’Aquila seismic sequence, from both the mainshock and several aftershocks. These data are especially important to investigate site effects in sediment-filled valleys during moderate earthquakes in epicentral area because well-documented observational studies are very limited in the literature. However, the main drawback for the study of site effects in the Aterno valley is the lack of a detailed knowledge of the geometry of the valley, soil layering and dynamic properties of materials. The main motivation for this study stems from the need to provide a reliable subsoil model of the valley coupled with high-quality strong motion data. Based on the above, in the framework of S4 project, a major effort was undertaken to get a trustworthy cross section of the valley by an ad hoc investigation, comprising geological and geotechnical surveys as well as an extensive geophysical campaign, characterized by both active and passive measurements. These results were complemented by additional geological and geotechnical data available in the literature. By merging all the information acquired, a 2D subsoil model of the transversal section of the upper Aterno valley has been produced. The valley is characterised by an asymmetric shape with a shallower rock basement at the western edge of the valley that deepens at the valley centre. Moreover, based on the results of geophysical tests, representative V_s values were assigned to the different lithologic units forming the alluvial deposits filling the valley. Shear wave velocity is a fundamental parameter for ground response studies and it is also effective in classifying the accelerometric station from a seismic point of view. The 2D model may be therefore, considered a benchmark model for future studies of site effects. It will offer the possibility to examine site effects with a complex underlying geology and to validate the results of numerical simulations of site response analyses with the numerous observations from earthquake recordings, both for weak and strong ground motion conditions.     

Exploration of S-wave velocity profiles at strong motion stations in Eskisehir,Turkey, using microtremor phase velocity and S-wave amplification

We have explored 1D S-wave velocity profiles of shallow and deep soil layers over a basement at strong motion stations in Eskisehir Province, Turkey. Microtremor array explorations were conducted at eight strong motion stations in the area to know shallow 1D S-wave velocity models. Rayleigh wave phase velocity at a frequency range from 3 to 30 Hz was estimated with the spatial autocorrelation analysis of array records of vertical microtremors at each station. Individual phase velocity was inverted to a shallow S-wave velocity profile. Low-velocity layers were identified at the stations in the basin. Site amplification factors from S-wave parts of earthquake records that had been estimated at the strong motion stations by Yamanaka et al. (2017) were inverted to the S-wave velocities and Q-values of the sedimentary layers. The depths to the basement with an S-wave velocity of 2.2 km/s are about 1 km in the central part of the basin, while the basement becomes shallow as 0.3 km in the marginal part of the basin. We finally discussed the effects of the shallow and deep sedimentary layers on the 1D S-wave amplification characteristics using the revealed profiles. It is found that the shallow soil layers have no significant effects in the amplification at a frequency range lower than 3 Hz in the area.     

Estimation of Background Noise for International Monitoring System Seismic Stations

v--vThe prototype International Data Centre (IDC) in Arlington, Virginia has been acquiring data from seismic stations at locations designated in the Comprehensive Test-Ban Treaty for the International Monitoring System (IMS) since the start of 1995. A key characteristic of these stations is their background noise levels and their seasonal and diurnal variability. Since June 1997 an automated sample selection effort has collected over 700,000 individual noise sample spectra from 39 primary and 57 auxiliary stations. Monthly median and 5 and 95 percentile estimates have been calculated for each channel of every station. Compatibility of median spectra obtained for the same station and channel in the same month for two different years confirms the consistency of the noise-sampling algorithm used. A preliminary analysis of the results shows strong (more than a factor of two) seasonal variation at a quarter of all stations. Strong diurnal variations at half of the sites indicate that many of the selected sites are poorly located with respect to cultural noise sources. The results of this study are already being used to evaluate station quality, improve those processes that require background noise values, such as automatic association and requesting auxiliary station data, and to improve the estimation of station and network detection and location thresholds.     

江苏地电台站渗流方位角和极化方位角差异特征分析研究

利用江苏地电台站的地电场和地磁数据,计算各个台站的极化方位角和渗流方位角,其中,渗流方位角分别采用谐波振幅和峰谷值2种方法进行计算,并对计算结果进行分析。结果表明:(1)新沂台和南京台由于台站覆盖层较浅,极化方位角存在线性极化和较多的非线性极化现象;而高邮台和海安台由于覆盖层较厚,极化方位角多数存在非线性极化现象。(2)用峰谷值法与用谐波振幅法计算所得的渗流方位角,两者相差不大,与理想模型间的差距也较小;而极化方位角与模型差距较大。因此,用渗流方位角尝试进行地震映震关系研究是基本可靠的。(3)高邮台的渗流方位角与高邮-宝应M_S4.9地震间具有较好的对应关系。     

GPS Observation and Study on Recent Crustal Movement in the Northeastern Area of Pamir,Xinjiang

Calculation of repeated observation data at the densified GPS monitoring network in northeastern area of Pamir together with data from IGS stations in the peribphyery of the area yielded the movement rate of more than 40GPS station sites in the area,and,hence,the recent crustal deformation rate pattern and time series of fiductial GPS stations in the area were obtained.The result indicates that the principal movement direction of the GPS station sites is NNW,basically diagonal to the strike of Tianshan fold belt,i.e.a normal compression occurs in the Tianshan region.The movement pattern near Jiashi and its southwestern zone is some different from that of station sites in their surrounding areas,indicating a certain relation of tectonic deformation in Jiashi area to seismic activity during last years.The movement rate of station sites in the periphery of Taim basin less varies and its direction is basically consistent It indicates less or basically no deformation within Tarim basin.     

川西地区台阵环境噪声瑞利波相速度层析成像

2006年中国地震局地质研究所地震动力学国家重点实验室在川西地区(26°N~32°N,100°E~105°E)布设了由297台宽频带数字地震仪组成的流动观测台阵.利用该密集台阵29°N以北156个台站2007年1~12月份的地震环境噪声记录和互相关技术,我们得到了所有台站对的面波经验格林函数和瑞利波相速度频散曲线,并进一步反演得到了观测台阵下方2~35 s周期的瑞利波相速度分布图像．本文结果表明,观测台阵覆盖的川滇地块、松潘-甘孜地块和四川盆地的地壳速度结构存在显著差异,具体表现为：(1)短周期（2~8 s）相速度分布与地表构造特征相吻合,作为川滇地块、松潘-甘孜地块和四川盆地之间的边界断裂,龙门山断裂带和鲜水河断裂带对上述三个地块上地壳的速度结构具有明显的控制作用,四川盆地前陆低速特征表明相应区域存在较厚的（约10 km）沉积盖层;(2)中周期（12~18 s）相速度分布表明,川滇地块和松潘-甘孜地块中上地壳速度结构存在明显的不均匀横向变化,并形成了尺度不同且高、低速相间的分块结构,而四川盆地中地壳整体上已经表现出相对高速;(3)长周期（25~35 s）相速度分布表明,松潘-甘孜地块,特别是川滇地块中下地壳表现为广泛的明显低速异常,意味着它们的中下地壳相对软弱,而四川盆地的中下地壳呈现整体性的相对高速,意味着四川盆地具有相对坚硬的中下地壳,并且以汶川地震的震中为界,龙门山断裂带的地壳结构显示了北段为高速异常,南段为低速异常的分段特征．     

Impacts of Urbanization and Station-relocation on Surface Air Temperature Series in Anhui Province, China

Since the 1990s, many meteorological stations in China have passively “entered” cities, which has led to frequent relocation and discontinuity in observational records at many stations. To study the impacts of urbanization on surface air temperature series, 52 meteorological stations in Anhui Province were chosen based firstly on a homogeneity test of the time series, and then their surrounding underlying surfaces during different decades were identified utilizing Landsat Multispectral Scanner images from the 1970s, Landsat Thematic Mapper images from 1980s and 1990s, and Enhanced Thematic Mapper images after 2000, to determine whether or not the station “entered” city, and then these stations were categorized into three groups: urban, suburban, and rural using Landsat-measured land use/land cover (LULC) around the station. Finally, variations in annual mean air temperature (T _mean), maximum air temperature (T _max), and minimum air temperature (T _min) were analyzed in urban-type stations and compared to their surrounding rural-type stations. The results showed that, in Anhui Province over the past two decades, many rural stations experienced urbanization and changed into urban or suburban locations. This process is referred as the “city-entering” phenomena of stations. Consequently, many of the latest stations were relocated and moved to currently rural and suburban areas, which significantly influenced the continuity of observational records and the homogeneity of long-term trends. Based on homogeneous data series, the averaged annual T _mean, T _max, and T _min over Anhui Province increased at a rate of 0.407, 0.383 and 0.432 °C decade^?1 from 1970 to 2008. The strongest effect of urbanization on annual T _mean, T _max, and T _min trends occurred at urban stations, with corresponding contributions of 35.824, 14.286, and 45.161 % to total warming, respectively. This work provides convincing evidences that (1) urban expansion has important impacts on the evaluation of regional climate change, (2) high spatial resolution images of Landsat are very useful for selecting reference climate stations for evaluating the potential urban bias in the surface air temperature data in certain regions of the continents, and (3) meteorological observation adjustments of station-relocation-induced inhomogeneities are essential for the study of regional or global climate change.     

CORRELATION PATTERNS IN CHARACTERISTICS OF SPATIALLY VARIABLE SEISMIC GROUND MOTIONS

A methodology for the investigation of the spatial variation of seismic ground motions is presented; data recorded at the SMART-1 dense instrument array in Lotung, Taiwan, during Events 5 and 39 are used in the analysis. The seismic motions are modeled as superpositions of sinusoidal functions, described by their amplitude, frequency, wavenumber and phase. For each event and direction (horizontal or vertical) analysed, the approach identifies a coherent, common component in the seismic motions at all recording stations, and variabilities in amplitudes and phases around the common component sinusoidal characteristics, that are particular for each recording station. It is shown that the variations in both the amplitudes and the phases of the motions at the station locations around the common component characteristics contribute significantly to the spatially variable nature of the motions, and, furthermore, they are correlated: increase in the variability of the amplitudes of the motions recorded at individual stations around the common amplitude implies increase in the variability of the phases around the common phase. The dispersion range of the amplitude and phase variability around their corresponding common components appear also to be associated with physical parameters. The spatially variable arrival time delays of the waveforms at the stations due to their upward travelling through the site topography, in addition to the wave passage delays identified from signal processing techniques, constitute another important cause for the spatial variation of the motions; their consideration in the approach facilitates also the identification of the correlation patterns in the amplitudes and phases. © 1997 by John Wiley & Sons, Ltd.     

柴达木盆地东部地震地面运动放大效应

柴达木盆地是青藏高原东北部大型断陷山间盆地,该地区的流动观测记录了2008年11月10日发生于大柴旦附近的M_W6.3地震。和附近的基岩上的记录相比,盆地内部的记录显示出非常显著的地面运动放大效应,表现为峰值速度的增大、持续时间的延长,其呈现出长持续时间的后续震相。傅里叶频谱分析表明盆地内部显著的后续震相的频率和直达波相比较低,地面质点运动轨迹图显示后续震相为面波运动特征。为了解释地面运动的差异,构建二维模型,通过交错网格高阶有限差分方法计算了地震波在盆地内部的传播过程,结果显示盆地内部低速层的存在造成直达波的放大以及多次反射与转换,盆地边缘结构造成的波的相干叠加产生了强烈的次生面波,其低频、大振幅、长持续时间的特征是盆地内部地面运动放大的主要原因。     

On reproducibility of relative locations of recurrent earthquakes recovered from similarity of their waveforms at a single station

Relative locations of recurrent earthquakes can be recovered from values of cross-correlation coefficients between their waveforms at a single station, provided that an idea of Menke (1999) that the values depend inversely exponentially from separation between the sources is true. An algorithm developed in Gnyp (2013) based on this idea was applied here for recovering relative locations of the two groups of the 2005–2006 recurrent earthquakes in Mukacheve, Ukraine, from the sets of cross-correlation coefficients between their records at the four local network stations to find out whether the locations were in principle reproducible independently at each of the stations, at epicentral distances from about 6 to 64 km. The resulting sets of locations occurred in good agreement with each other, which may be considered as an experimental evidence for validity of the Menke’s dependence over the corresponding range of distances. The reproducibility of relative locations has been analyzed in terms of existence and uniqueness of the inverse problem’s solution, and in connection with data completeness and the influence of noise on it.     

Optimum strong-motion array geometry for source inversion—II

Optimum strong-motion array geometry for source inversions is again determined for each of three types of earthquake faults: strike-slip, dip-slip and offshore subduction thrust. The method is the same as employed in a previous study;¹ however, use of a complete Green's function in an elastic half-space provides better results for engineering practice. It is found that the complete Green's function is capable of stabilizing the accuracy of an inversion solution obtained using theoretical seismograms, regardless of the differences in array configuration. The optimum strong-motion array for a strike-slip fault is characterized by stations well distributed in azimuth, while the optimum array for a dip-slip event has stations arranged in a grid-shaped form. The array geometries obtained here are grossly similar to those in the previous study,¹ which were derived using only the far-field S waves, and are more consistent with those proposed at the 1978 International Workshop on Strong-Motion Earthquake Instrument Arrays.².     

用剪切波分裂研究台湾北部地壳各向异性

研究主要使用台湾北部13个地震台站记录到的1991年7月～2002年12月的波形数据,采用剪切波分裂SAM分析方法,对台湾北部地区的剪切波分裂特征进行了研究.发现位于宜兰盆地内的台站的快剪切波优势偏振方向为近E-W方向,而位于山脉(西部麓山带、雪山山脉和中央山脉)的台站的快剪切波优势偏振方向为NNE向或NE方向.位于海...     

三峡水库运行前后长江中下游干流及两湖的径流过程变化

为探索三峡水库运行前后长江中下游干流及两湖径流过程的变化及其驱动因素,利用宜昌、监利、大通、七里山、湖口共5个水文站的流量资料,分析了各站径流过程的变化特征及其成因。结论:(1)各站年径流量均减少,但除七里山站之外,其它各站减少比例均小于10%且变化不显著;(2)干流各站月径流量最大减幅发生在10月,而七里山站、湖口站分别发生在7月、4月;(3)干流各站月径流量最大增幅发生在3月,而七里山站、湖口站分别发生在1月、6月;(4)宜昌站,1—4月径流量增加是三峡水库入库径流增加和水库调度的共同作用结果,6—8月径流量减少的主因是三峡水库入库径流量减少,5、9、10月径流量变化的主因是三峡水库调度;(5)监利站,径流量的变化与宜昌站表现出高度的一致性,但冬季各月径流量的增幅均大于宜昌站;(6)大通站,4—6月径流量变化方向与湖口站一致,其它月份变化方向均与宜昌站一致。(7)七里山站,7月径流量减少的主因是洞庭湖流域来水减少,9、10月径流量减少的主要原因是荆江分流减少,但洞庭湖流域来水减少也是重要原因。(8)湖口站,4、5月径流量减少的主因是流域降水减少,9、10月径流量减少的主要原因是鄱阳...     

Impacts of massive pumping on flood processes in polders—with a case study on Lannihu basin in the Dongting Lake District,China

Flooding risk in polders is dictated by not only rainfall, topography, and land use, but also massive pumping. Unfortunately, existing models are inadequate for resolving floods as water transfer due to pumping is insufficiently accounted for. Here an improved hydrological model (MGB-MP) is proposed under the framework of the large-scale hydrological model (MGB) based on the principle of water balance, explicitly incorporating massive pumping within a polder and also out to external rivers. The proposed model is calibrated and validated for the Lannihu basin, a typical polder with an area of 1353 km² and 126 pumping stations in the Dongting Lake District, China and surrounded by Xiangjiang River and Zishui River. The model performs fairly well, with Nash-Sutcliffe efficiencies concerning water levels over 0.76 for the calibration and over 0.73 for the validation. The model is applied to the Lannihu basin under different pumping station settings and rainfall scenarios to unravel how and to what extent massive pumping affects the flood processes as characterized by water levels and discharge hydrographs. It is shown that massive pumping considerably alters the discharge hydrographs and accordingly leads to substantial decrease in the water levels of rivers, which are independent unit-polders, due to water transfer between unit-polders within the basin and out of the basin. The closer the unit-polders are to pumping stations, the more the water levels in unit-polders decrease. The water levels in unit-polders away from a pumping station is affected by the pumping station capacity to a greater extent than the pumping station's threshold water level for initiating pumping.     

Potential Occurrence of MTBE and BTEX in Groundwater Resources of Amman–Zarqa Basin,Jordan

This study investigates potential occurrence, distribution, and sources of the newly added gasoline oxygenate, methyl‐tert‐butyl ether (MTBE) and the petroleum derivatives benzene, toluene, ethylbenzene, and xylenes called collectively, BTEX, in Jordan's heavily populated Amman–Zarqa Basin (AZB). It presents the first data on the levels of MTBE and BTEX in the aquifers of this basin. One hundred and seventy‐nine (179) groundwater wells were sampled near petrol service stations, oil refinery storage tanks, car wrecks, bus stations, and chemical industries at different locations in the basin. Headspace GC and purge and trap GC–MS were utilized to determine the target substances in the samples. Concentrations of BTEX varied between no‐detection (minimum) for all of them to 6.6 µg/L (maximum) for ethylbenzene. MTBE was found in few samples but none has exceeded the regulated levels; its concentrations ranged between no‐detection to 4.1 µg/L. However, though the contamination levels are very low they should be considered alarming.     

A dense array experiment for the observation of waveform perturbations

We analyze observed three-component seismograms recorded by a dense array of 21 short-period receivers (f>2 Hz) installed on a homogeneous Quaternary plain with 1.5×3.5 km² horizontal dimensions. The experiment was conducted in the district of Caille, in SE France for 3 months. The observed waveform data are related to diffraction and refraction of the direct wavefield by different geological structures and topographic morphologies surrounding the array. The origin of the coda energy following the energetic Pg or Sg arrivals is either influenced by heterogeneities close to and inside of the valley, for instance thrust of Jurassic limestones over Cretaceous marls in the northern and southern limits of the valley, or by heterogeneities located outside of the valley such as the interface unconformity between crystalline basement and sedimentary cover in the northern border of the Maures–Tanneron ranges. The main energetic arrivals of the seismograms arrive at the array with a definite range of azimuths that is not correlated with the event back-azimuth. Moreover, these arrivals do not appear with azimuths related to the characteristic NS and EW valley orientations, and thus may be related to heterogeneities located outside of the valley.