The induced seismicity after Alborz Dam impoundment: implications to the active tectonic in northern Iran

Document Type : Research


1 Assistant Professor, Department of Geology, Payame Noor University, Sari, Iran

2 M.Sc. Graduated, Department of Geology, Payame Noor University, Tehran, Iran


The study of induced earthquakes is important in different aspects. One of the most important aspects of reservoir-induced seismicity is the possibility of the triggering of a strong ground motion after reservoir impoundment. This study aims at finding a relationship between the ML5.2 Babolkenar earthquake of 2012 January 11 and the Alborz dam’s impoundment. The Alborz dam is located on the northern flank of the Alborz Mountains. The Alborz dam is an embankment dam and its capacity is 8.6 million cubic meters. Its ridge lies on the tectonized Early Tertiary rocks. The 5.2 ML Babolkenar earthquake of 2012 January 11 is assessed with the use of waveforms provided by different seismic stations. The analyses of waveformes show the centroid depth of 20km and epicenter located at 36.357N, 52.788E. Its thrust fault focal mechanism is in agreement with the fault kinematics inversion implying the reactivation of one of segments of the Khazar fault. Although the event in question could be considered as an ordinary tectonic earthquake, some lines of existing evidence associated with the reservoir impounding data encouraged the authors to consider it as anthropogenic seismicity. Firstly the reservoir is located in the area between two very active thrust fault that contains fractured and permeable rocks. Secondly is the drastic change of seismicity after the onset of reservoir impoundment represents the change of the b-value before, during, and after reservoir impoundment. Thirdly the trigger of the event just after a lowering of the water level. In these cases, earthquake rupture may be interpreted by two different mechanisms: 1) an immediate, undrained, elastic response to the reservoir load and/or an instantaneous pore pressure change in the vicinity of the reservoir due to an undrained response )Skempton, 1954(. The governing equation to explain the undrained response is . In this equation, the change of the undrained compression is related to the Skempton coefficient (B) and the average normal stress at a point located on the fracture ( ). Although the lowering of the reservoir’s water table may change rapidly , the B-coefficient on the other hand, could not change immediately. Therefore the pore pressure is increased even if the water level is decreased leading to satisfying the Coulomb failure criteria. 2) A delayed and/or undrained response due to diffusion of pore pressure. The curve of the accumulative of the seismicity versus the time history of the reservoir impoundment shows two cycles of rising and lowering the water level. Each cycle starts with a gently linear increase of seismicity corresponding to the rising of the water table and ends with a period of increase of the seismicity in response to the water level drawdown. The authors think the diffusion of pore pressure is responsible for the increasing of the seismicity after reservoir impoundment. Given the r=25km as the maximum distance between dam and centroid depth and ∆t=12 months as average time elapsed between the time of that water level increase and the time of the earthquake, the hydraulic diffusivity c=4.95 m^2/sec is estimated by . By the other approach in which the area is affected by the after shakes (here ~100km2) is considered as the r2 in the above-mentioned equation, the c-value takes 0.8 m^2/sec. This c-value discrepancy may be raised by lack of data especially due to lack of local seismic network to the survey of microearthquakes and also unavailable bore-hole to install piezometers measuring the water pore pressure versus time. The c-coefficient may be a useful parameter that is applicable to predict the same reservoir-induced or triggered seismicity in future times. Finally, as mentioned above, the b-value decreased just before the event. This may be in an analogy to the decreasing of the b-value in the other anthropogenic seismic activity i.e. mining tremors, usually called rockburst, and a strong ground motion governed by tectonic activity. The b-value after the event M5.2 increased. Finally, the authors strongly recommend the government plan detailed geological, geophysical and geotechnical studies before and after the impoundment of a new water reservoir to monitor the induced seismicity in the northern flank of Alborz. This may help to mitigate the damages due to the probable triggering of the active faults within this very seismically active belt.


Main Subjects

رسولی، ر.، و هوشمندان، ح.، 1389، بررسی روند تغییرات فشار آب منفذی در بدنه و پی سد البرز با استفاده از نتایج ابزار دقیق و مقایسه آن با نتایج تحلیلی، اولین همایش ملی سازه – زلزله-ژئوتکنیک آذر 1389، مازندران-بابلسر.
Alavi, M., 1996, Tectonostratigraphic synthesis and structural style of the Alborz mountain system in northern Iran. Journal of Geodynamics 21, 1–33.
Aziz Zanjani, A., Ghods, A., Sobouti, F., Bergman, E., Mortezanejad, G., Priestley, k., Saeed Madanipour, S. and Rezaeian, M., 2013, Seismicity in the western coast of the South Caspian Basin and the Talesh Mountains. Geophys. J. Int. 1093/gji/ggt299.
Bell, M.L. and Nur, A., 1978, Strength changes due to reservoir-induced pore pressure and stresses and application to Lake Oroville. J. Geophys. Res. 83, 4469 – 4483.
Carey-Gailhardis, E. and Mercier, J.-L., 1987, a numerical method for determining the state of stress using focal mechanisms of earthquake populations: application to Tibetan teleseisms and microseismicity of Southern Peru, earth planet. Sci. Lett. 82, 165–179.
Carder, D. S., 1970, Reservoir loading and local earthquakes in engineering seismology the works of man. In: W.M. Adams (Editor), Engineering Geology Case Histories, No. 8. Geological Society of America, Denver, Colo., 51-61.
Doloei, J. and Roberts, R., 2003, Crustal and uppermost mantle structure of Tehran region from teleseismic P-waveform receiver function analysis, Tectonophysics, 364, 115–133.
Gupta, H., K. and Rastogi, B., K., 1976, Dams and Earthquake, Elsevier, the Netherlands, 229pp.
Gupta, H. K., 2002, a review of recent studies of triggered earthquakes by artifcial water reservoirs with special emphasis on earthquakes in Koyna, India. Earth-Science Reviews 58(3–4), 279–310.
Jackson, J., Priestley, K., Allen, M. and Berberian, M., 2002, Active tectonics of the South Caspian Basin, Geophys. J. Int., 148(2), 214–245.
Kaviani, A., 2004, La chaine de collision continentale du Zagros (Iran): structure litospherique par analyse De donees sismologiques, PhD Thesis, University of Joseph Fourier-Grenoble I.
Maa, Xu., Westmana, E., Slakera, B., Thibodeaub, D. and Counter, D., 2018, The b-value evolution of mining-induced seismicity and mainshock occurrences at hard-rock mines, International Journal of Rock Mechanics and Mining Sciences 104.
Malekzade, Z., 2018, Block rotation induced by the change from the collision to subduction: 1333 Implications for active deformations within the areas surrounding South Caspian Basin. 1334 Marin Geology 404, 111-129.
Nuannin, P., 2006, The Potential of b-value Variations as Earthquake Precursors for Small and Large Events. Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 183.
Pavlou, K., Drakatos, G., Kouskouna, V., Makropoulos, K. and Kranis, H., 2016, Seismicity study in Pournari reservoir area (W. Greece) 1981–2010.
Ritz, J.-F., Nazari, H., Ghassemi, A., Salamati, R., Shafei, A., Soleymani, S. and Vernant, P., 2006, Active transtension inside Central Alborz: a new insight into northern Iran—southern Caspian geodynamics, Geology, 34, 477–480.
Ruiz-Barajas, S., Santoyo, M. A., BenitoOterino, M. B., Alvarado, G. E. and Climent, A., 2021, Stress transfer patterns and local seismicity related to reservoir water-level variations. A case study in central Costa Rica.
Simpson, D. W., Leith, W. S. and Scholz, C. H., 1988, Two types of reservoir induced seismicity, Bull. Seism. Soc. Am. 78, 2025–2040.
Skempton, A. W., 1954, The pore pressure coefficients A and B, Geotechnique 4, 143–147.
Tatar, M., 2001, Etude siesmotectonique de deux zones de collision continentale: Le Zagros Central et l’Alborz (Iran), PhD Thesis, University of Joseph Fourier-Grenoble I.
Tatar, M., Jackson, J., Hatzfeld, D. and Bergman, E.A., 2007, The 28 May 2004 Baladeh earthquake (Mw 6.2) in the Alborz, Iran: implications for the geology of the south Caspian basin margin and for the seismic hazard of Tehran, Geophys J. Int., 170, 249–261.
Talwani, P., 1997, on the nature of reservoir-induced seismicity, Pageoph 150, 473–492.
Talwani, P., L. Chen, and K. Gahalaut, 2007, Seismogenic permeability, ks, J. Geophys. Res., 112, B07309, doi:10.1029/ 2006JB004665.
Trifu, C., 2002, the Mechanism of Induced Seismicity; Springer RaseI AG.
Valoroso, L., Improta, L., Chiaraluce, L., Di Stefano, R., Ferranti, L., Govoni, A. and Chiarabba, C., 2009, Active faults and induced seismicity in the Val d’Agri area (Southern Apennines, Italy). Geophysical Journal International 178 (1), 488e502.
Williams-Stroud S., 2014, A geological approach to seismicity b-values: implications for hazard assessment. SEG SEG International Exposition and 87th Annual Meeting.
Wyss, M., and Wiemer, S., ZMAP, 2001, A tool for analyses of Seismicity Pattern. ZMAP cook book, ETH, Zurich.