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. ... mehrFor 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.