A statistical procedure for the analysis of seismotectonically induced hydrochemical signals: A case study from the Eastern Carpathians, Romania

Changes in the stress field of an aquifer system induced by seismotectonic activity may change the mixing ratio of groundwaters with different compositions in a well, leading to hydrochemical signals which in principle could be related to discrete earthquake events. Due to the complexity of the interactions and the multitude of involved factors the identification of such relationships is a difficult task. In this study we present an empiric statistical approach suitable to analyse if there is an interdependency between changes in the chemical composition of monitoring wells and the regional seismotectonic activity of a considered area. To allow a rigorous comparison with hydrochemistry the regional earthquake time series was aggregated into an univariate time series. This was realized by expressing each earthquake in form of a parameter “e”, taking into consideration both energetic (magnitude of a seismic event) and spatial parameters (position of epi/hypocentrum relative to the monitoring site). The earthquake and the hydrochemical time-series were synchronised aggregating the e-parameters into “earthquake activity” functions E, which takes into account the time of sampling relative to the earthquakes which occurred in the considered area. For the definition of the aggregation functions a variety of different “e” parameters were considered. The set of earthquake functions E was grouped by means of factor analysis to select a limited number of significant and representative earthquake functions E to be used further on in the relation analysis with the multivariate hydrochemical data set. From the hydrochemical data a restricted number of hydrochemical factors were extracted. Factor scores allow to represent and analyse the variation of the hydrochemical factors as a function of time. Finally, regression analysis was used to detect those hydrochemical factors which significantly correlate with the aggregated earthquake functions.This methodological approach was tested with a hydrochemical data set collected from a deep well monitored for two years in the seismically active Vrancea region, Romania. Three of the hydrochemical factors were found to correlate significantly with the considered earthquake activities. A screening with different time combinations revealed that correlations are strongest when the cumulative seismicity over several weeks was considered. The case study also showed that the character of the interdependency depends sometimes on the geometrical distribution of the earthquake foci. By using aggregated earthquake information it was possible to detect interrelationships which couldn't have been identified by analysing only relations between single geochemical signals and single earthquake events. Further on, the approach allows to determine the influence of different seismotectonic patterns on the hydrochemical composition of the sampled well. The method is suitable to be used as a decision instrument in assessing if a monitoring site is suitable or not to be included in a monitoring net within a complex earthquake prediction strategy.     

The influence of seismotectonics on precursory changes in groundwater composition for the 1995 Kobe earthquake, Japan

It is shown that temporal variation in hydrochemistry may be attributed, in part, to continuous seismotectonic activity occurring before the onset of the 1995 “Kobe” earthquake, Japan, challenging the “one earthquake–one signal” hypothesis with respect to potential precursory signals to this devastating event. A possible continuous seismotectonic influence on chloride and sulphate ion-concentration is evaluated with aggregate earthquake-information by transforming a multivariate earthquake time series (including coordinates and magnitude) into a one-dimensional time series (considering geometric relationships between earthquakes and the well-site). A piecewise analysis of ion-concentration and seismotectonic-activity time series compares trends and change points between the two time series: a positive correlation (before the proposed onset of the preparation stage) is followed by a negative correlation (during the proposed preparation stage) which, in turn, is succeeded by a positive correlation (after the heaviest aftershock sequences). This suggests that seismotectonic processes occurring before the Kobe earthquake dynamically influenced aquifer characteristics, leading to temporal variations in the hydrochemical time series. Accordingly, a dynamic change in the mixing ratios of waters with different hydrochemical characteristics is proposed as a mechanism for explaining observed variation. The research can be extended by considering vectors in the stress field that lead to changes in the aquifer-well system.The influence of seismotectonics on precursory changes in groundwater composition for the 1995 Kobe earthquake, Japan.     

Reactive Transport Process of Earthquake-induced Hydro-chemical Changes in Guanding Thermal Spring, Western Sichuan, China

Earthquake-related hydrochemical changes in thermal springs have been widely observed; however, quantitative modeling of the reactive transport process is absent. In the present study, we apply reactive transport simulation to capture the hydrochemical responses in a thermal spring following the Wenchuan Ms 8.0 and Lushan Ms 7.0 earthquakes. We first constrain deep reservoir geothermal fluid compositions and temperature by multicomponent geothermometry, and then a reactive geochemical transport model is constructed to reproduce the hydrochemical evolution process. The results show that the recharge from the shallow aquifer increases gradually until it reaches a peak because of the permeability enhancement caused by the Lushan earthquake, which may be the mechanism to explain the earthquake-related hydrochemical responses. In contrast to the postseismic effect of the Wenchuan earthquake, the chemical evolution can be considered as hydrochemical anomalies related to the Lushan earthquake. This study proves that the efficient simulation of reactive transport processes is useful for investigating earthquake-related signals in hydrochemical time series.     

Long-term hydrochemical earthquake precursor studies at the Koyna-Warna reservoir site in western India

Koyna-Warna region of western India is an active seismic zone due to the Reservoir Triggered Seismicity (RTS). Earthquake precursor studies are carried out monitoring hydrochemical and stable isotope signatures in the groundwater from 15 bore wells since January 2005, for more than 12 years (January 2005 to February 2017). Depth of these boreholes ranges from 100 to 250 m. Cyclic or temporal variation in hydrochemistry is observed in few sensitive wells in Koyna region. The Govare well in Koyna is found to be most sensitive and the observed hydrochemical cycle is closely associated with local earthquakes of M > 5. The earthquakes M <5 occurring either in Warna cluster or close to the observation wells, did not generate hydrochemical precursory changes. The increase in hydrochemistry is hypothesized as mixing of two aquifer waters with different hydrochemistry. It is noted that a precursory hydrochemical cycle is observed during first quarter of 2015, but no earthquake M > 5.0 occurred till date. The cyclic changes in hydrochemistry, however, indicate on-going earthquake processes and an impending earthquake of M > 5 in the region.     

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M _w ) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01° and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.     

Seismotectonics and seismicity of the Silakhor region, Iran

This paper deals with seismotectonic and seismicity of the Silakhor region that shows high seismic activity in western Iran. Silakhor is a vast plain with several villages and cities of Dorud and Borujerd and a small town of Chalanchulan that were destroyed and/or damaged many times by large earthquakes. This paper addresses the historical and instrumental earthquakes and their causative faults, seismotectonic provinces and seismotectonic zones of the region. Available seismic data were normalized by means of time normalization technique that resulted in the magnitude-frequency relation for the Silakhor area and estimation of the return period of earthquakes with different magnitudes. Some active faults in this region include the Dorud fault, the main Zagros thrust, the Galehhatam fault, the Sahneh fault and others. Among them, the Dorud fault is an earthquake fault and is the cause for most of the large and intermediate earthquakes in the region. The return period of large earthquakes with magnitudes greater than 7.0 (Ms) is very low, however, the occurrence of destructive earthquakes is greater in the region than in the neighboring provinces. The study proves the high seismicity of this zone and it is required to develop an accurate national plan for future building and reinforcement of the existing buildings in this region.     

Deformation and displacement fields of the Himalayan arc,seismicity and large earthquake hazards in its eastern segment and vicinity

The Himalayan arc is a type of plate margin similar to an island arc and is a world-famous region of tectonic and seismic activities, where a series of large earthquakes have occurred in historical time. In this paper, the vertical deformation and horizontal displacement fields of the Himalayan arc are theoretically derived from the viewpoint of the collision between the Indian and Eurasian plates. In the light of the observed data, the seismicity, earthquake focal mechanism, seismotectonic and geomorphological features of the arc and its vicinity are reasonably explained. The characteristics of seismicity and the possibility of earthquakes with magnitude above 8 occurring in this region in the future are studied.     

First-order regional seismotectonic model for South Africa

A first-order seismotectonic model was created for South Africa. This was done using four logical steps: geoscientific data collection, characterisation, assimilation and zonation. Through the definition of subunits of concentrations of earthquake foci and large neotectonic and structural domains, seismotectonic structures, systems and domains were created. Relatively larger controls of seismicity exist between the Great Escarpment and the coast. In the south, this region is characterised by large aeromagnetic anomalies and large EW trending faults. In the west, it is characterised by the NW–SE trending Wegener stress anomaly, radial-trending dykes and earthquake clusters. In the east, it is characterised by a large neotectonic domain where several large historical earthquakes occurred. In the centre of South Africa, several clusters of earthquake activity are found, often related to mining activity. Further north, seismicity is related to both mining activity and neotectonic deformation. This work contributes to the development of a seismotectonic model for South Africa by (1) bringing together, digitally, several data sets in a common GIS platform (geology, geophysics, stress, seismicity, neotectonics, topography, crustal and mantle structure and anisotropy), (2) understanding the significance of data sets for seismotectonic zonation and limitations thereof and (3) obtaining a reasonable regional model for use in seismic hazard assessments.     

Seismo-Tectonic Characteristics of the North Anatolian Fault Zone Between Akyazi and Duzce (Bolu,Turkey)

The active North Anatolian fault zone (NAFZ) presents very complex seismotectonic activity. The occurrence of the Abant earthquake in 1957 (Ms = 7.1) and the Mudurnu earthquake in 1967 (Ms = 6.8) are only two examples of several seismic events associated with intense tectonic activity of the NAFZ. Statistical analyses of earthquakes in an area extending between 30° 30′ to 31° 30' E Long. and 40° 15′ to 41° 00′ N Lat. reveal that epicenters generally were shallow. However, a few deep epicenters also were located, some of which reached a depth of 30 km. The epicenters were found to concentrate in a zone lying between the Duzce and Akyazi Plain to the north of Almacik Mountain and in the Adapazari Plain. The Northern Anatolian fault displays an en echelon character in the area, except for the eastern part, where it extends as a single segment. The en echelon character of the NAFZ is interpreted as a structure distributing the potential energy and consequently reducing the intensity of earthquakes, giving rise to micro-earthquakes of magnitudes less than 4.2.     

Seismic Hazard Analysis for the Bangalore Region

Indian peninsular shield, which was once considered to be seismically stable, is experiencing many earthquakes recently. As part of the national level microzonation programme, Department of Science and Technology, Govt. of India has initiated microzonation of greater Bangalore region. The seismic hazard analysis of Bangalore region is carried out as part of this project. The paper presents the determination of maximum credible earthquake (MCE) and generation of synthetic acceleration time history plot for the Bangalore region. MCE has been determined by considering the regional seismotectonic activity in about 350 km radius around Bangalore city. The seismotectonic map has been prepared by considering the faults, lineaments, shear zones in the area and historic earthquake events of more than 150 events. Shortest distance from the Bangalore to the different sources is measured and then peak ground acceleration (PGA) is calculated for the different source and moment magnitude. Maximum credible earthquake found in terms of moment magnitude is 5.1 with PGA value of 0.146 g at city centre with assuming the hypo central distance of 15.88 km from the focal point. Also, correlations for the fault length with historic earthquake in terms of moment magnitude, yields (taking the rupture fault length as 5% of the total fault length) a PGA value of 0.159 g. Acceleration time history (ground motion) and a response acceleration spectrum for the corresponding magnitude has been generated using synthetic earthquake model considering the regional seismotectonic parameters. The maximum spectral acceleration obtained is 0.332 g for predominant period of 0.06 s. The PGA value and synthetic earthquake ground motion data from the identified vulnerable source using seismotectonic map will be useful for the PGA mapping and microzonation of the area.     

Correlations between soil gas and seismic activity in the Generalized Haiyuan Fault Zone,north-central China

Radon and mercury concentrations were measured in 10 fault gas profiles in Generalized Haiyuan Fault. This paper aims to predetermine the potential seismic risk in different segments of the fault zone from the perspective of geochemistry. The background value and anomaly threshold were adopted and synthesized using the maximum value method and average method to calculate concentration intensity values of radon and mercury. Fault soil gas mercury and radon concentrations show a decreasing gradient from NW to SE indicating evident segmentation. Higher values are mostly distributed in the Maomao Mountain–Tiger Mountain fault and Jingtai area. Combined with the seismotectonic background of historical and recent earthquakes and the spatial distribution characteristics of b-values, the fault soil gas concentration intensity shows a close correlation with earthquake activity within the fault zone. Concentrations of fault gas are higher and the b-value lower in areas of strong seismic activity, and regions with weak seismic activity correspond to lower fault gas concentrations and higher b-values. It is thus considered that the Jingtai area may be more dangerous than the other areas. This paper could provide vital background information for earthquake prediction in the Generalized Haiyuan Fault Zone.     

On the swelling potential of compacted high plasticity clays

We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆ_max, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an M_S-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆ_max,RP_6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.     

2017年8月8日四川九寨沟7.0级地震发震构造浅析

2017年8月8日21时19分,四川阿坝州九寨沟县发生7.0级地震,震中位于巴颜喀拉块体东边界虎牙断裂和东昆仑断裂带东段塔藏断裂交汇区域,地震构造背景较为复杂。地震导致了房屋和道路破坏、滑坡崩塌。根据高分辨率卫星影像解译、阶地坎变形的测量和测年数据得到：塔藏断裂东段晚第四纪以来以左旋走滑为主,兼逆分量,水平滑动速率为2.7~4.1 mm/yr,垂直滑动速率为0.56~0.6 mm/yr。结合此次地震的主余震分布、主震震源机制解等综合结果,初步建立了三维发震构造模型,分析认为此次地震属于走滑型地震,主破裂倾角57°~77°,发震断层可能是塔藏断裂的一条分支,是青藏高原块体向东推挤的一次地震事件。基于历史地震、活动断裂和形变观测方面的研究,巴颜喀拉块体具备显著的强震构造背景,对于该块体边界带周缘的强震活动和变形需要继续关注。     

Tectonic position,geological manifestations of sources,and aftershock processes for a series of strong earthquakes of 2016–2017 in central Italy

Based on seismological, geodynamic, and seismotectonic data on a series of strong earthquakes which occurred in 2016–2017 in central Italy, the tectonic positions and geological manifestations of these events are determined. The aftershock processes of these earthquakes are studied. It is found that epicentral clouds of mainshocks and aftershocks almost overlapped each other, so these earthquake sources in aggregate can be classified as a high-magnitude earthquake swarm.     

Association of subsurface radon changes in Alaska and the northeastern United States with earthquakes

Changes of radon concentration in the ground have been reported to correlate with seismic activity in many places over widely varying distances from the earthquake epicenters, sometimes over 100's and occasionally over 1000's of km depending on the magnitude M of the earthquakes. For proper use of such correlations in forecasting the location and size of impending earthquakes, it would be useful to establish a relation between M and the maximum distance x_M at which radon can be significantly altered by preseismic effects.We have monitored radon in the ground continuously at Blue Mountain Lake, NY starting Dec. 3. 1975 and at three locations in Alaska—Icy Bay, Yakataga, and Sand Point (starting Nov. 12, 1979; March 5, 1980; and June 29, 1980 respectively). Sample correlations from visual examination of the recent radon record and tentative use of a proposed scaling relation show a number of possible earthquake-related signals at these locations. The data are compatible with scaling relations that were derived from two separate models of premonitory elastic strains. In this work x_M = 10^0.48M, where x_M is in units of km and M ≥ 3. Since tilt and strain signals correlate similarly with magnitude and distance, it is likely that most earthquake-related radon signals are mechanically induced.     

Correlation Fractal and Multifractal Characteristics of Seismic Activity in the Taiwan Area, China

Based on the analysis of newly collected data of plate tectonics, distribution of active faults and crustal deformation, the Taiwan area is divided into two seismic regions and six seismic belts. Then, correlation fractal dimensions of all the regions and belts are calculated, and the fractal characteristics of hypocenteral distribution can be quantitatively analyzed. Finally, multifractal dimensions Dq and f(α) are calculated by using the earthquake catalog of the past 11 years in the Taiwan area. This study indicates that (1) there exists a favorable corresponding relationship between spatial images of seismic activity described with correlation fractal dimension analysis and tectonic settings; (2) the temporal structure of earthquakes is not single but multifractal fractal, and the pattern of Dq variation with time is a good indicator for predicting strong earthquake events.     

Space-time distribution patterns of destructive earthquakes in the inner belt of central Japan: Activity intervals and locations of earthquakes

Based on a block structure model of the inner belt of central Japan, an examination was conducted of the space-time distribution patterns of destructiv magnitudes M 6.4 or greater (M =Japan Meteorological Agency Scale). The distribution patterns revealed a periodicity in earthquake activit seismic gaps. Major NW—SE trending left-lateral active faults divide the inner belt of central Japan into four blocks, 20–80 km wide. The occurrenc A.D. with M ≥ 6.4, which have caused significant damage, were documented in the inner belt of central Japan. The epicenters of these earthquakes close to the block boundaries.

Using the relationship between the magnitude of earthquakes which occurred in the Japanese Islands and the active length of faults that generated them, movement is calculated for each historical earthquake. Space—time distributions of earthquakes were obtained from the calculated lengths, the latitud of generation. When an active period begins, a portion or segment of the block boundary creates an earthquake, which in turn appears on the ground surf active period ends when the block boundary generates earthquakes over the entire length of the block boundary without overlapping.

Five seismic gaps with fault lengths of 20 km or longer can be found in the inner belt of central Japan. It is predicted that the gaps will generate ea magnitudes of 7.0. These data are of significance for estimating a regional earthquake risk over central Japan in the design of large earthquake resist

The time sequences of earthquakes on the block boundaries reveal a similar tendency, with alternating active periods with seismic activity and quiet pe activity. The inner belt of central Japan is now in the last stage of an active period. The next active period is predicted to occur around 2500 A.D.     

太阳活动周与全球大震

为研究太阳活动与全球大震的关系, 引入一个无量纲的"地震能量函数√G", 并分析研究了1681—2011故年间全球M≥7.0大震的能量释放的时间序列.由此发现全球大震在太阳活动周4个阶段的分布和活动度, 随震级的强度而异.提出地壳对太阳风暴加卸载响应模式, 用于解释此现象: 通过考察最近331 a, 得出全球共发生了10个M≥9.0超级巨震的时空分布特征, 特别是太阳活动峰年期间没有发生过超级巨震.该研究结果可为判断全球大震提供参考.      

On the long-term seismic hazard analysis in the Zhangjiakou–Penglai seismotectonic zone, China

The Zhangjiakou–Penglai seismotectonic zone (ZPSZ) lies in the northern part of North China and extends along the Zhangjiakou–Beijing–Tianjin–Bohai Bay–Penglai–Yellow Sea. It is about 900 km long and some 250 km wide in a northwest direction. The great Sanhe-Pinggu (M_S=8.0) earthquake occurred on September 1679 and the Tangshan (M_S=7.8) earthquake on July 1976 caused serious economic and life losses. According to some differences in crust structure and regional tectonic stress field, the ZPSZ is divided into western and eastern segment by the 117°E line for study on long-term seismic hazard analysis. An analysis of Gutenberg–Richter's empirical relation of earthquake-frequency and time process of historic and recent earthquakes along the eastern and western segments shows that the earthquake activity obeys a Poisson process, and these calculations indicate that the earthquake occurrence probability of M_S=6.0–6.9 is 0.77–0.83 in the eastern segment and the earthquake occurrence probability of M_S=7.0–7.9 is 0.78–0.80 in the western segment of the ZPSZ during a period from 2005 to 2015.