Three dimensional velocity structure of the Koyna-Warna region using local earthquake tomography

A dense seismic network (～100 stations) was operated in the Koyna-Warna region from January 2010 to May 2010, that allow us to collect 400 high-quality local earthquake data of magnitude less than 4. In this region, the fault structure and tectonic setting that accommodate the induced seismicity is not well understood. To investigate the seismotectonics of the region, we have inverted 7826 P- and 7047 S-P arrival times for 3-D Vp and Vp/Vs tomographic models along with hypocenters parameters in the region. Although, Dixit et al. (2014) have performed 3-D local earthquake tomography with double-difference tomography code using catalog differential time data. In this paper, Simulps14 code on the same data set is applied. For better approach P arrival time and S-P travel times are inverted directly for Vp, Vp/Vs variations and earthquake locations. High Vp ～5.9 to 6.5 and low to high Vp/Vs ～1.69-1.74 imaged in the hypocenter region. These features interpreted as a fluid bearing rock mass under high pore pressure. It is also observed that below the trap basement form a local topography depression between the Koyna and Warna Reservoirs. To the South of the Warna reservoir, intense seismic activity defines a major cluster of ～ 5 km width at 3 to7 km deep, located under the trap, where the basement is deepening. Such regions are inferred to be associated with the seismically active faults zones. The obtain velocity anomalies are reliable down to a depth of 10 km. This is also confirmed by the analysis of three resolution parameters viz. Hit count, Derivative Weight sum (DWS) and Resolution Diagonal Elements (RDE).     

Method for determination of impending earthquake focal parameters based on daylight surface displacement data

Doklady Earth Sciences -     

岩土工程中的设计地震动参数确定

与岩土工程相关的设计地震动参数的规定分散在各行业的工程设计规范之中,比较了主要技术规范中设计地震动参数规定之间的差别,讨论了在岩土工程的勘察与设计中确定地震动参数的方法和内容,分析和研究了确定设计地震动参数时不同规范之间需要协调的主要问题,指出了在目前岩土工程中的设计地震动参数确定方面的需要重点关注的参数选择的合理性,提出了具有地震工程专业基础（注册地震安评工程师）的注册土木工程师（岩土）确定具体设计地震动参数应更符合工程建设的实际需要的建议。     

地震与地震前兆数字化观测基础

本文介绍了数字观测中的模数转换、采样、量化和传输的基本原理并以地形变数字观测为例介绍了采样定理的应用     

Crustal seismic velocity structure near Faial and Pico Islands (AZORES), from local earthquake tomography

The Azores archipelago (Portugal) is located on an oceanic plateau, in a geodynamic environment prone to intense seismic and volcanic activity. In order to investigate the crustal structure in this region, we have conducted a local earthquake tomography study in the area of the islands of Faial, Pico and S. Jorge using data recorded in July 1998. The July 9th 1998 earthquake, near Faial Island, triggered an aftershock sequence of thousands of events that lasted for several months and were recorded by a total of 14 stations located on the three islands surrounding the epicentral area. In the upper crustal layers, consistency is seen between the tomographic results and the islands' surface volcanic units. Beneath the Faial central volcano a low Vp (< 6.0 km/s) anomaly roughly located at 3–7 km depth, suggests a connection to the plumbing system, possibly the presence of a magma chamber. In NE Faial, a high Vp (> 6.3 km/s) body was found located at mid-lower crust, most likely an intrusion of gabbroic composition, that is bordered by the registered seismic activity; its shape suggesting a tectonic controlled mechanism. The relocated hypocenters, together with the overall analysis of the Tomographic model, suggest a tectonic segmentation of Faial Island. The crustal thickness under the islands volcanic buildings of the Faial–Pico area was estimated at around 14 km.     

Studies on the determination of shear wave velocity in sands

Laboratory investigations such as resonant column, cyclic triaxial, and torsional shear tests are usually conducted on undisturbed and reconstituted sand samples in order to determine their shear modulus, which can be used to compute shear wave velocity. However, these methods are extremely cumbersome and indirect, they employ cost-intensive instrumentation, and they require trained manpower. These problems can be overcome by using bender elements, which yield the shear wave velocity of sands directly, easily, and quite rapidly. In this context, efforts made by earlier researchers in developing empirical relationships (based on bender elements and resonant column techniques) which can be employed to compute shear wave velocity in sands are worth consideration. Although these relationships take into account parameters like effective stress and void ratio, they do not directly incorporate fundamental characteristics of sands such as their state (i.e. dry or saturated) and the particle size. Thus a generalized relationship, which can be used to compute shear wave velocity in sands in either the dry or saturated state, has been developed. Experiments were conducted on samples of sand and glass beads of different sizes and the results were used to demonstrate the utility and efficiency of the relationship. Details of the methodology developed for this purpose are also presented in this paper.     

Density distribution functions and correlations of earthquake parameters

The statistical analysis of the source parameters of 9 earthquake sets of different types (aftershocks, scattered events, swarms) and of different seismic regions shows that the density distribution function (ddf) of the linear dimension l of a fault is represented by a negative power law, as well as the ddf of the static stress drop σ and of the scalar seismic moment M_o. It is then suggested, and tentatively verified, that also the ddf of the root mean square ground acceleration, defined as a function of l and σ, may be represented by a negative power law and that, at least in the cases examined, it scales like the ddf of σ. It is seen that the variability of the static stress drop is significant from one region to another, as is well known, but it seems remarkable also in the same seismic region (in particular in California, σ varies by several orders of magnitude) and in the different sets of events of a given region (as observed again for California). It is hypothesized that a correlation, although weak, between the stress drop and the linear dimension of a fault exists and the analyses seem not to contradict that σ may be a decreasing function of l. Finally, it is suggested that the seismicity of a region may be represented two-dimensionally as a function of the ddf of the stress drop and of the linear dimension of a fault instead of the classic b and b_o values.     

芦山地震发震构造及其与汶川地震关系讨论

芦山地震发生在巴彦喀拉块体与华南块体之间龙门山推覆构造带南段。野外考察表明,芦山地震在震中区没有形成具有构造地质意义的地震地表破裂带,仅在各山前陡坡地带出现平行于山麓陡坡的张性地裂缝、山地基岩崩塌、滑坡等边坡震动失稳现象和震动引起的砂土液化现象。重新定位的芦山地震余震分布、震源机制解和地表构造地质等分析表明,芦山地震的发震断层为一条现今尚未出露地表、其上断点仍埋藏在地下9 km以下地壳中的一条盲逆断层,走向212°,倾向NW,倾角38°±2°,上断点以上至地表的构造变形符合断层扩展背斜模型。根据汶川地震和芦山地震的余震空间分布、地震破裂过程、深浅构造关系等差异反映出它们是分别发生在龙门山推覆构造带中段和南段的两次独立地震破裂事件。     

An integration-based estimation approach for the determination of slug test parameters under various flow geometries

A simple parameter estimation procedure, designated as integration-based estimation (IBE), was introduced to determine the hydraulic properties of an aquifer using slug test models subjected to certain flow geometries such as radial and spherical flows. The basic idea behind the proposed IBE approach is to link an integration value at pre-defined normalized head levels for field data with that of a theoretical type curve. The IBE method removes the need for the implementation of the classical graphical matching process which would be ineffective to acquire aquifer parameters for non-ideal aquifer conditions. As the second aspect of this study, a new decision tool was suggested to determine the suitable slug test model to be utilized for the site data since diagnosing the flow character properly is of crucial importance for following a convenient analysis procedure. The estimation performance and limitation of the proposed IBE method were tested for several slug test scenarios including radial and spherical flow models with a number of synthetically generated data sets as well as a field application. Results reveal that the IBE together with the identification methodology not only is able to retrieve aquifer parameters as reliable as the existing techniques in the literature but also diagnoses the flow character precisely as demonstrated in this study.     

1605年海南岛琼州大地震及其发震构造的初步探讨

琼北地区铺前−清澜断裂是一条关系到1605年琼山7½级大地震发震构造判断的重要断裂。为准确厘定该断裂的最新活动特征,针对琼北地区第四纪玄武岩盖层和巨厚海相砂沉积等地质条件,采用大吨位震源的浅层人工地震勘探和小间距钻孔联合地质剖面探测相结合的方法,分别在海口江东三江镇岐山头村和东寨村展开探测。结果表明：(1) 铺前−清澜断裂断错了多个标志地层,最浅断错至全新统烟墩组淤泥层,上断点埋深10 m左右,钻探所揭示的断层面明显具有走滑兼具逆冲性质,是一条具有一定宽度、多分支、多期次活动的大规模断裂带;(2)钻孔联合地质剖面所揭示的地层深度范围内,被断错地层的位移量随着深度增加而加大,铺前−清澜断裂自8 346~7 153 a cal BP以来有过活动,全新世累计垂直位移量4~5 m,垂直位移速率为0.53~0.63 mm/a,推断其为1605年琼山大地震的发震断裂。本研究取得的铺前−清澜断裂全新世活动的新证据,为国土空间规划和区域地震危险性评价提供了科学依据。

     

The velocity structure of the Cariaco sedimentary basin,northeastern Venezuela,from refraction seismic data and possible relationship to earthquake hazard

The main damage from the July 9, 1997, Cariaco earthquake (Ms=6.8) was concentrated in the town of Cariaco and surrounding villages, which are located in the western part of the Cariaco sedimentary basin, close to the Gulf of Cariaco. Casanay, located at the eastern end of the sedimentary basin, suffered considerably less damage. The El Pilar fault, a right-lateral strike-slip fault that generated the earthquake, runs parallel to the southern border of the valley and crosses both towns. The determination of the velocity structure of the basin is the main objective of this study. Seismic refraction data were recorded along three lines, one of them along-strike and two perpendicular to the valley axis in the northern and southern bedrocks. Beneath Cariaco, approximately 1 km thick Quaternary sediments with seismic velocities of 1.9–2.1 km/s and bedrock velocities of more than 4 km/s were observed. The thickness of the Quaternary sediments varies within the basin, and Pleistocene sediments outcrop beneath Casanay. The increased thickness of the unconsolidated, water-saturated Quaternary sediments, together with the difference in the quality of buildings prior to the earthquake, probably is responsible for the damage pattern of the Cariaco earthquake.     

Assessment of the statistical earthquake hazard parameters for NW Turkey

We performed a probabilistic analysis of earthquake hazard input parameters, NW Turkey covers Gelibolu and Biga Peninsulas, and its vicinity based on four seismic sub-zones. The number of earthquakes with magnitude M ≥ 3.0 occurred in this region for the period between 1912 and 2007 is around 5130. Four seismic source sub-zones were defined with respect to seismotectonic framework, seismicity and fault geometry. The hazard perceptibility characterization was examined for each seismic source zone and for the whole region. The probabilities of earthquake recurrences were obtained by using Poisson statistical distribution models. In order to determine the source zones where strong and destructive earthquakes may occur, distribution maps for a, b and a/b values were calculated. The hazard scaling parameters (generally known as a and b values) in the computed magnitude–frequency relations vary in the intervals 4.28–6.58 and 0.59–1.13, respectively, with a RMS error percentage below 10 %. The lowest b value is computed for sub-zone three indicating the predominance of large earthquakes mostly at Gelibolu (Gallipoli) and north of Biga Peninsula (southern Marmara region), and the highest b value is computed for sub-zone two Edremit Bay (SW Marmara region). According to the analysis of each seismic sub-zone, the greatest risk of earthquake occurrence is determined for the triangle of Gelibolu–Tekirda? western part of Marmara Sea. Earthquake occurrence of the largest magnitude with 7.3 within a 100-year period was determined to be 46 % according to the Poisson distribution, and the estimated recurrence period of years for this region is 50 ± 12. The seismic hazard is pronounced high in the region extending in a NW–SE direction, north of Edremit Bay, west of Saros Bay and Yenice Gönen (southern Marmara region) in the south. High b values are generally calculated at depths of 5–20 km that can be expressed as low seismic energy release and evaluated as the seismogenic zone.     

Physico-chemical models of the internal structure of partially differentiated Titan

We analyze models of the internal structure of Titan, a large icy satellite of the Saturn system. Calculations are carried out using information on the mass, mean density, moment of inertia, orbital parameters, and elastic properties of the satellite obtained by the Cassini–Huygens mission, as well as geochemical data on the composition of chondrite materials, equations of state of water and ices I, III, V, VI, and VII, and thermodynamic models for conductive heat transfer in the outer icy crust and of global convection in the interior zones of the satellite. The analysis of the models shows that models of partially differentiated Titan are most consistent; they include an outer water–ice shell, an intermediate ice–rock mantle, and an inner rock–iron core. It is shown that for the models of this type the maximum thickness of the water–ice shell is 460–470 km; it can be composed of an outer conductive crust of Ih ice 80–110 km thick and a subsurface water ocean 200–300 km deep. The maximum radius of the central rock–iron core of Titan can reach ~1300 km. The thickness of Titan’s ice–rock mantle does not exceed 2100 km at a density of 1.22–2.64 g/cm³. The model of partially differentiated Titan is feasible in the moment of inertia range of 0.312 < I/MR ² < ~0.350.     

Statistical analysis of seismic processes on the basis of the diffusion approach

On the basis of the diffusion approach, statistical analysis of seismicity is carried out in order to reveal the peculiarities characteristic of random walk processes. In accordance with the hypothesis on connectivity of earthquakes and energetic criterion, the expansion of the seismic process into a set of sequences is made; each sequence is a Brownian process with specific spatial, temporal, and energetic scales. The statistical distributions of sequences by the number of links, total energies, distances, time, and flight velocities between events are built.     

基于HYDRUS模拟的ABCD模型变量及参数物理基础研究

以ABCD模型为代表的概念性水文模型是量化水文循环过程的主要手段之一,但是很少有研究检验ABCD模型变量和参数的物理基础及其与传统土壤水力参数的关系。以鄂尔多斯盆地台格庙矿区作为研究区,首先使用HYDRUS-1D软件构建站点尺度土壤-植物-大气连续体模型模拟获取非冻结期“真实的”实际蒸散量,其作为目标函数对ABCD模型进行校正,然后对比ABCD模型模拟的土壤水储量和土壤水渗漏量与SPAC模型模拟的结果,通过开展大量情景模拟探究ABCD模型中关键参数a和b与传统土壤水力参数之间的关系。结果表明：在站点尺度,ABCD模型可以用于模拟非冻结期月时间尺度的实际蒸散量和土壤水渗漏量;ABCD模型中参数a与饱和含水率、残余含水率和曲线形状参数具有较强的线性相关性,与饱和渗透系数呈对数关系;参数b与饱和含水率和曲线形状参数n具有较强的线性相关性,与饱和渗透系数和曲线形状参数α呈对数关系。研究结果可以提高对概念性水文模型变量和参数物理基础的认识,拓展水文模型在水文地质中的应用。     

Source parameters of the Taimyr earthquake of June 9, 1990

The results of calculation of the seismic moment tensor (in the double-couple approximation) for the earthquake that occurred in the northeastern Taimyr Peninsula on June 9, 1990, from amplitude spectra of the surface waves, are presented. The obtained results are of interest for works related to construction of geodynamic models for the Taimyr region, for studies of the stress-strain state of the Earth’s crust, and for assessments of the seismic hazard of this territory.     

S-wave velocity structure of the North China from inversion of Rayleigh wave phase velocity

We constructed the S-wave velocity structure of the crust and uppermost mantle (10–100 km) beneath the North China based on the teleseismic data recorded by 187 portable broadband stations deployed in this region. The traditional two-step inversion scheme was adopted. Firstly, we measured the interstation fundamental Rayleigh wave phase velocity of 10–60 s and imaged the phase velocity distributions using the Tarantola inversion method. Secondly, we inverted the 1-D S-wave velocity structure with a grid spacing of 0.25° × 0.25° and constructed the 3-D S-wave velocity structure of the North China. The 3-D S-wave velocity model provides valuable information about the destruction mechanism and geodynamics of the North China Craton (NCC). The S-wave velocity structures in the northwestern and southwestern sides of the North–South Gravity Lineament (NSGL) are obviously different. The southeastern side is high velocity (high-V) while the northeastern side is low velocity (low-V) at the depth of 60–80 km. The upwelling asthenosphere above the stagnated Pacific plate may cause the destruction of the Eastern Block and form the NSGL. A prominent low-V anomaly exists around Datong from 50 to 100 km, which may due to the upwelling asthenosphere originating from the mantle transition zone beneath the Western Block. The upwelling asthenosphere beneath the Datong may also contribute to the destruction of the Eastern Block. The Zhangjiakou-Penglai fault zone (ZPFZ) may cut through the lithosphere and act as a channel of the upwelling asthenosphere. A noticeable low-V zone also exists in the lower crust and upper mantle lid (30–50 km) beneath the Beijing–Tianjin–Tangshan (BTT) region, which may be caused by the upwelling asthenosphere through the ZPFZ.     

束流型拦砂坝溢流口及其关键参数确定

为进一步优化拦砂坝溢流口体型,确定溢流口几何参数取值,通过物理模拟实验研究了不同拦砂坝溢流口体型条件下泥石流过坝的流态、坝后冲刷深度、消能率变化规律,并根据拦砂坝埋深设计标准与松散体水下临界平衡关系,探讨了溢流口收缩率的合理取值范围,结果表明:采用斜面或者大倾角的反弧型溢流口时,坝后泥石流落点与坝脚的距离近、冲刷深度大,而采用无倾角的反弧型溢流面时,不仅可以在一定程度上减小冲刷深度,而且大大增加了坝后泥石流落点与坝脚的距离;在泥石流规模与溢流面曲率半径相同情况下,适当增大侧向收缩率,有利于增强泥石流体与坝后动床之间的相互作用,提高泥石流通过坝后冲刷坑的消能率;当溢流口收缩率在0.2~0.6之间时,冲刷深度满足设计要求,且泥石流跌落点距坝脚较远,冲刷坑发展不会危及坝体安全。     

On the estimation of earthquake maximum magnitude based on a joint analysis of seismological and geothermal parameters

To improve the reliability of earthquake prediction, it is important to determine those layers in the crust of the earth in which the foci of most strong earthquakes occur. An analysis of regional earthquake catalogues shows that in depressions strong events originate within a specific layer in the middle part of the crust. An investigation of the geothermal regime in the depressions under study shows that the thickness and location of this layer depend on geothermal parameters which correlate with the strongest earthquake magnitude in particular areas. Thus, it seems possible to use geothermal and seismological parameters for estimating the maximum magnitude of earthquakes for at least some areas.