Seismicity changes associated with reservoir loading

Changes in seismic activity have been related to the filling of large reservoirs in over thirty cases. These changes range from variations in the level of micro-earthquake activity detectable only with instruments of high sensitivity to destructive earthquakes with magnitudes greater than 6. On the other hand, the filling of many other large reservoirs has not been accompanied by increased seismicity.

A number of factors may contribute to the generation or absence of post-impounding seismicity. Increased vertical stress due to the load of the reservoir and decreased effective stress due to increased pore pressure can modify the stress regime in the reservoir region. Whether or not these stress changes are sufficient to generate earthquake activity will depend on a complex interaction of the induced stress with the state of pre-existing stress near the reservoir, and on the geologic and hydrologic conditions at the site. The combined effect of increased vertical load and increased pore pressure will have the greatest tendency to increase activity in regions where the maximum compressive stress is vertical (normal faulting). In regions where the minimum compressive stress is vertical (thrust faulting) increased stress due to a vertical load should have a minimum effect. For all of the larger reservoir-induced earthquakes the stress system determined from fault plane solutions is in agreement with the pre-existing stress field in the region of the reservoir. These earthquakes are all of strike-slip or normal type, there being no reported cases of large induced earthquakes with thrusting mechanisms.

The potential for major changes in seismicity may be highest in regions of moderate strain accumulation (low to moderate natural seismicity). In areas of high strain accumulation and high levels of natural seismicity, the stress changes induced by the reservoir will be small compared to natural variations. In aseismic areas, with low strain accumulation, the reservoir-induced stresses may be insufficient to raise the stress level to a state of failure.     

Migration of seismicity and pore pressure diffusion for five decades long Koyna–Warna sequence and precursory nucleation process for earthquake forecasting

The continued reservoir-triggered seismicity for five decades in Koyna area has been attributed to southward migration of seismicity (during 1967–1992 near and south of Koyna dam and from 1993 onwards mostly near the new Warna reservoir). Spread of seismicity in the vicinity of reservoirs is attributed to pore-pressure diffusion. Moderate size Koyna–Warna earthquakes are found to nucleate at shallow depth (≤ 3 km) due to pore pressure caused by water level fluctuation of reservoir(s). The nucleation zone deepens along the critically stressed permeable fault zone to cause the occurrence of mainshock at the base of seismogenic layer (i.e. 5–10 km). The clustering of foreshocks up to 500 hr prior to several moderate size Koyna earthquakes of magnitude Mw 4–5 has been detected and used for quantifying the nucleation process. A static stress transfer by means of cascade model from one foreshock to next for the generation of foreshocks has been proposed for nucleation model. The nucleation process can be considered as an immediate earthquake precursor for the Koyna-Warna region.     

Coseismic responses and the mechanism behind MW 5.1 earthquake of March 14, 2005 in the Koyna–Warna region, India

An earthquake of M_w 5.1 occurred on March 14, 2005, in the seismically active Koyna–Warna region in western India, the site known for the largest reservoir triggered seismicity (RTS) in the world. For more than four decades, earthquakes with M  4.0 have occurred in this region at regular intervals. Impoundment of reservoirs and changes in lake levels can trigger earthquakes by two processes of stress modifications, namely direct loading effect of the reservoir and diffusion through various faults and fractures. In this paper we analysed the reservoir water level data at Koyna and Warna reservoirs prior to the occurrence of the March 14, 2005 earthquake, to explain the dominant mechanism behind its occurrence and its correlation with the observed coseismic changes. We conclude that the diffusion process, not the reservoir load effect, is the dominating mechanism triggering earthquakes in the region. The coseismic changes in deep well water levels sensitive to earth tides are found to be to the order of 1–12 cm.     

Faulting parameters of earthquakes in the New Madrid, Missouri, region

The advent of high-resolution digital seismic recording and advances in computer technology enable the combination of traditional regional seismic network observations with direct seismogram modeling to improve estimates of small earthquake faulting geometry, depth, and size. We illustrate a combined modeling approach using observations from three earthquakes that occurred within the environs of the New Madrid Seismic Zone: two Missouri earthquakes from September 26, 1990 and May 4, 1991; and the southern Illinois earthquake of February 5, 1994. We also re-examine the faulting geometry for two events from the 1960s that are inconsistent with the current estimate of the regional stress field. Based on direct modeling of the long-period seismograms associated with these events, we revise earlier estimates of the earthquake parameters for the March 3, 1963 and July 21, 1967 Missouri earthquakes. Comparing the new and revised results with existing earthquake mechanisms in the region, we find that tension-axes are generally aligned in a N-S to NW-SE direction, while the compression-axes trend in a NE to E direction. An interesting exception to this pattern are the March 3, 1963 and two nearby earthquakes that lie within a well-defined 30-km long left step in seismicity near New Madrid.     

三峡库区仙女山和九畹溪断裂带水库地震变化规律

三峡库区蓄水以来,长江周边出现了大量的水库地震,主要集中于巴东—泄滩—仙女山区域,目前已达上万次,最大震级为5.1级。通过野外地质调查及对已有水库地震数据进行研究,采用构造地质分析方法,对三峡库区仙女山和九畹溪断裂带水库地震空间上的迁移规律、时间上的周期规律以及微地震群的成因机理进行分析。结果表明：从时间上看,水库地震具有周期性,表现为长周期（与库水位相关）和短周期（与库水位快速波动相关）;从空间上看,水库地震具有迁移规律,受九畹溪断裂控制的触发型水库地震存在着逐渐向南迁移的特征,受仙女山断裂控制的触发型水库地震局限分布于仙女山断裂北延端点处,分布于仙女山断裂西侧（周坪乡附近）带状分布的水库诱发地震逐渐呈点状向南迁移。     

三峡近坝库段水库诱发地震特征分析

大型水利水电工程水库诱发地震历来受业主、工程技术人员的重视。三峡工程蓄水10多年来,在蓄水过程中,三峡库区较蓄水前地震明显增多。本文选择了三峡库区既有活动断裂,又有可溶岩和煤层分布的近坝库段作为研究对象,通过对库区地震活动与水位的关系和诱发地震空间展布特征分析,认为:三峡水库诱发地震,发生在有可溶岩和含煤地层部位的塌陷型地震占总数的90%,而非可溶岩和不含煤地层中只占总数的10%;近坝库段地震主要发生在仙女山断层北端,且有地震呈丛集活动的特征;九湾溪断层上盘远离长江部位地震零星分布,距长江较近地段地震相对增多,下盘地震极少。地震与水位的对应关系十分密切,三峡水位上升,长江两岸发生的小地震较多,当水位蓄至175m时,地震主要发生在仙女山断裂的北端和九湾溪断裂的上盘及北端。     

Analytical investigation of hydraulic fracture-induced seismicity and fault activation

More recent public discourse has taken place regarding the potential correlation between seismic activity and hydraulic fracturing in shale gas reservoirs. Public fears about the risk of seismicity stem mainly from past earthquakes induced by conventional deep injections because the two types of projects share similar mechanisms of rock failure and fault activation. Although previous earthquake risks associated with fluid injection were not serious, the situation would be far more problematic if hydraulic fracturing in a shale gas reservoir triggered a similar-sized earthquake due to potential environmental issues. In fact, almost all documented injection-induced earthquakes have been associated with long-duration and high-volume injection rather than short-term (hours) pressurization (e.g., hydraulic fracturing). In general, hydraulic fracturing operations mostly induce microseismic events through rock failure and activation of small fractures. Although shale reservoirs in tectonically active zones pose a high risk of inducing large-magnitude seismic activities, the internal geological conditions and external stimulation conditions are impossible to be satisfied simultaneously to trigger activation of an entire fault and to result in a destructive earthquake during hydraulic fracturing operations.     

Co-seismic gravity changes in the Koyna-Warna region: Implications of mass redistribution

Koyna-Warna Region (KWR) is one of the known sites for reservoir triggered seismicity. The continued triggered seismicity over the five decades is restricted to a region of about 600–700 sq. km, which provides a unique opportunity to monitor geophysical anomalies likely to be associated with seismicity of the region. Present study confers temporal gravity changes recorded by gPhone and GRACE satellite and interprets observed changes in conjunction with seismological, geodetic (cGPS) observations and groundwater level measurements. GRACE data suggest that seasonal vertical deformation due to hydrological loading is ～ 2 cm, which corroborates with continuous GPS observations. Seasonal hydrological loading of the region, which is in a phase of reservoir loading, might be influencing the critically stressed KWR leading to the seasonal seismicity of the region. The gPhone gravity data distinctly show co-seismic gravity signals for eight earthquakes of M_w > 2 and gravity anomalies show positive correlation on a logarithmic scale with earthquake released energy. To investigate the cause of gravity changes, an estimate is made for 14^th April 2012 earthquake for M_w 4.8 using fault dislocation model. The recorded gravity changes of 189 μGal by gPhone located at a distance of 28 km from the hypocentre is much more than the estimate of ～0.1 μGal calculated for M_w 4.8 Koyna earthquake. Therefore, it is inferred that co-seismic gravity signals for eight earthquakes are primarily caused due to redistribution of mass at shallow depth.     

The influence of the Nurek Reservoir on local earthquake activity

The Nurek Reservoir is located in an area of high seismicity. An average of 5–6 earthquakes of (Russian energy class)K = 10 and 1–2 earthquakes ofK = 11 per year occurred in the vicinity of the reservoir. The largest local earthquake recorded fell in the energy classK = 14 (1956).

At the end of 1972, after the water level had risen to 100 m, earthquake activity increased sharply. Three events ofK = 12 were recorded during one month and the total number of earthquakes (K 7) exceeded three times the former mean level of seismicity.

Spatial and temporal properties in the epicentral distribution showed a migration from the southwest toward the reservoir.

Earthquake mechanisms of 215 events ofK 9 since 1960 were analysed. Fault plane solutions for some of the earthquakes which occurred after the impoundment were found to be different from the ordinary ones.     

Estimation of Coulomb stress changes and induced earthquake occurrences prediction in the region of Three Gorges Reservoir,China

In the region of Three Gorges Reservoir (TGR) in China, there has been occurrence of several frequent earthquakes of moderate intensity since reservoir impounding occurred in 2003. These earthquakes are generally believed to be induced by reservoir impoundment and water-level variations. Usually, the geo-stress will change, when natural earthquakes occur. Following this principle, this paper adopted the rate and state theory to simulate and estimate Coulomb stress changes in the TGR region and obtained the pattern of Coulomb stress changes with time and the event sequence as well as the distribution of Coulomb stress changes in space. First, the TGR regional catalogue was analyzed and processed, leading to quantification of the magnitude of completeness and all of the parameters that are used in the stress–seismicity inversion process, including the reference seismicity rates, characteristic relaxation time, fault constitutive parameters, and stress rates. Second, the temporal evolution of the stress changes in different time windows was computed and analyzed, and it was found that there is an association between the Coulomb stress changes and rates of increase in the cumulative number of earthquakes. In addition, the earthquake occurred in November 2008 (M _S = 4.1) was analyzed and attempted to simulate the distribution of stress changes in space through the stress–seismicity inversion model. The results proved that the modeled area coincides with the historical area of earthquakes that occurred after 2008. Finally, a prediction was made about the earthquake productivity rates after 2015, which showed a declining earthquake rate over time that ultimately returned to the background seismicity. This result is essentially in agreement with Omori’s law. To conclude, it is rational to use the stress-inversion method to analyze the relationship between induced earthquake seismicity and local stress changes as well as to simulate the area of earthquake occurrence and productivity rates of reservoir-induced earthquakes.     

What triggers Koyna region earthquakes? Preliminary results from seismic tomography digital array

The cause for prolific seismicity in the Koyna region is a geological enigma. Attempts have been made to link occurrence of these earthquakes with tectonic strain as well as the nearby reservoirs. With a view to providing reliable seismological database for studying the earth structure and the earthquake process in the Koyna region, a state of the art digital seismic network was deployed for twenty months during 1996–97. We present preliminary results from this experiment covering an area of 60 × 80 km² with twenty seismic stations. Hypocentral locations of more than 400 earthquakes confined to 11×25 km² reveal fragmentation in the seismicity pattern — a NE — SW segment has a dip towards NW at approximately 45°, whilst the other two segments show a near vertical trend. These seismic segments have a close linkage with the Western Ghat escarpment and the Warna fault. Ninety per cent of the seismicity is confined within the depth range of 3–10 km. The depth distribution of earthquakes delimits the seismogenic zone with its base at 10 km indicating a transition from an unstable to stable frictional sliding regime. The lack of shallow seismicity between 0 and 3 km indicates a mature fault system with well-developed gouge zones, which inhibit shallow earthquake nucleation. Local earthquake travel time inversion for P- and S-waves show ≈ 2% higher velocity in the seismogenic crust (0–10 km) beneath the epicentral tract relative to a lower velocity (2–3%) in the adjoining region. The high P- and S-wave velocity in the seismogenic crust argues against the presence of high pressure fluid zones and suggests its possible linkage with denser lithology. The zone of high velocity has been traced to deeper depths (≈ 70 km) through teleseismic tomography. The results reveal segmented and matured seismogenic fault systems in the Koyna region where seismicity is possibly controlled by strain build up due to competent lithology in the seismic zone with a deep crustal root.     

Water level fluctuations due to earthquakes in Koyna-Warna region,India

Earthquakes cause a variety of hydrological phenomena, including changes in the ground water levels in bore wells. The Koyna region in the peninsular shield of India, hitherto considered stable in terms of seismic activity, has been active since 1967. More recently, the earthquakes have been localized to the newly impounded Warna reservoir, which is located south of Koyna, where a burst of seismicity occurred in 1993. The region continues to remain seismically active even after four decades. Twenty-one bore wells were drilled around the seismic source volume in the region to observe water level changes resulting from earthquake phenomena. Our studies have shown coseismic anomalous water level changes to be associated with the moderate earthquakes of April 25, 1997 and February 11, 1998. Our results show that changes in the ground water level in bore wells are correlated with micro-earthquake activity, both preceding and following moderately sized earthquakes. The results have implications in enhancing our understanding of earthquake mechanisms.     

Seismicity around Brazilian dam reservoirs

More than 30 cases of seismicity associated with dam reservoir sites are known throughout the world. Despite the lack of data in some areas, where seismicity occurred after reservoir impounding, there have been distinct seismic patterns observed in seismic areas after dam projects implantation. This has demonstrated that reservoir loading can trigger earthquakes. A mechanism of earthquake generation by reservoir impounding is proposed here with particular application to the Brazilian cases and to areas subject to low confining stress conditions in stable regions. Six artificial lakes are described and the associated earthquake sources are discussed in terms of natural or induced seismicity. Earthquake monitoring in Brazil up to 1967, when Brasilia's seismological station started operation, was mainly based in personal communications to the media. Therefore, there is a general lack of seismic records in relatively uninhabited areas, making it difficult to establish a seismic risk classification for the territory and to distinguish natural from induced seismicity. Despite this, cases reported here have shown an alteration of the original seismic stability in dam sites, after reservoir loading, as observed by the inhabitants or records from Brasilia's seismological station. All cases appear to be related to an increase in pore pressure in permeable rocks or fracture zones which are confined between impermeable rock slabs of more competent rock. It is apparent that some cases show some participation of high residual stress conditions in the area.     

Lineament fabric from airborne LiDAR and its influence on triggered earthquakes in the Koyna-Warna region,western India

We present the results of the first airborne LiDAR survey flown in the Koyna-Warna region and examine the relationship between the lineament fabric and the ongoing seismicity in the region. Our studies reveal that earthquakes of M≥4.0 for the period 1968 to 2016 are strongly correlated with a 10 km wide N-S fracture zone, which not only represents the surface expression of seismically active basement faults, but also act as conduits for water percolation between the Koyna and Warna reservoirs. A decreasing trend in the annual distribution of earthquakes was observed from 1985. A new burst of seismic activity in 1993 followed the impoudment of the Warna reservoir. We report a change in annual seismicity pattern, where seismicity peaks during September and December in the pre-Warna period, with a new peak emerging during March-April subsequent to the impoundment of Warna reservoir. A model is proposed to explain the seismicity along dominant N-S lineaments and the impact of Warna reservoir impounding which altered the hydrogeologic regime in the region.     

Koyna,India, an ideal site for near field earthquake observations

The Koyna earthquake of M 6.3 on December 10, 1967 is the largest artificial water reservoir triggered earthquake globally. It claimed ～ 200 human lives and devastated the Koyna township. Before the impoundment of the Shivaji Sagar Lake created by the Koyna Dam, there were no earthquakes reported from the region. Initially a few stations were operated in the region by the CentralWater and Power Research Station (CWPRS). The seismic station network grew with time and currently the National Geophysical Research Institute (NGRI), Hyderabad is operating 23 broadband seismographs and 6 bore hole seismic stations. Another reservoir, Warna, was created in 1985, which provided a further impetus to Reservoir Triggered Seismicity (RTS). Every year following the monsoon, water levels rise in the two reservoirs and there is an immediate increase in triggered earthquakes in the vicinity of Koyna-Warna reservoirs in the months of August–September. Peak RTS is observed in September and later during December.Another spurt in triggered earthquakes is observed during the draining of the reservoirs in the months of April- May. A comparative study of RTS earthquake sequences and the ones occurring in nearby regions made it possible to identify four common characteristics of RTS sequences that discriminate them from normal earthquake sequences. As the RTS events continue to occur at Koyna in a large number in a limited area of 20 km x 30 km, at shallow depths (mostly 2 to 9 km), the region being accessible for all possible observations and there being no other source of earthquakes within 100 km of Koyna Dam, it was suggested to be an ideal site for near field observations of earthquakes. This suggestion was discussed by the global community at an ICDP sponsored workshop held at Hyderabad and Koyna in 2011. There was an unanimous agreement about the suitability of the site for deep scientific drilling; however, a few additional observations/experiments were suggested. These were carried out in the following three years and another ICDP workshop was held in 2014, which totally supported setting up a borehole laboratory for near field investigations at Koyna. Location of a Pilot Bore-hole was decided on the basis of seismic activity and other logistics. The 3 km deep Pilot Borehole was spudded on December 20, 2016 and completed on June 11, 2017.     

Continued seismic activity at the Koyna reservoir site, India

Following the impounding of the Shivaji Sager Lake in 1962, tremors became prevalent in the Koyna region, considered previously to be aseismic. During ensuing years the tremor frequency appears to have been dependent on the rate of increase of water level, maximum water level reached, and the period for which high levels were retained. This culminated in a burst of seismic activity from September 1967 to January 1968 following the record water levels in the reservoir and included the earthquake of September 13, 1967 with magnitude 5.5 and the damaging December 10, 1967 earthquake of magnitude 6.0. During the next five years water levels were kept low and no significant earthquakes occurred subsequent to the October 29, 1968 earthquake of magnitude 5.

The reservoir was filled to maximum capacity during September 1973 and this was followed by a conspicuous increase in seismic activity which included an earthquake of magnitude 5.1 on October 17, 1973. However, seismic activity during 1973 was much less severe than that of 1967. This relative decrease in seismicity may indicate that (a) the “threshhold level” for relatively large magnitude earthquakes had increased; (b) a major portion of the accumulated strains had been released; and/or (c) the importance of the longer period of high loading in 1967. Similar observations have been made at other seismically active reservoir sites.     

Time independent seismic hazard analysis in Alborz and surrounding area

The Bayesian probability estimation seems to have efficiencies that make it suitable for calculating different parameters of seismicity. Generally this method is able to combine prior information on seismicity while at the same time including statistical uncertainty associated with the estimation of the parameters used to quantify seismicity, in addition to the probabilistic uncertainties associated with the inherent randomness of earthquake occurrence. In this article a time-independent Bayesian approach, which yields the probability that a certain cut-off magnitude will be exceeded at certain time intervals is examined for the region of Alborz, Iran, in order to consider the following consequences for the city of Tehran. This area is located within the Alpine-Himalayan active mountain belt. Many active faults affect the Alborz, most of which are parallel to the range and accommodate the present day oblique convergence across it. Tehran, the capital of Iran, with millions of inhabitants is located near the foothills of the southern Central Alborz. This region has been affected several times by historical and recent earthquakes that confirm the importance of seismic hazard assessment through it. As the first step in this study an updated earthquake catalog is compiled for the Alborz. Then, by assuming a Poisson distribution for the number of earthquakes which occur at a certain time interval, the probabilistic earthquake occurrence is computed by the Bayesian approach. The highest probabilities are found for zone AA and the lowest probabilities for zones KD and CA, meanwhile the overall probability is high.     

美国北加州吉塞尔斯诱发地震与热汽田特征的关系

吉塞尔斯地震活动可能是因蒸汽开发引起的.水在一个承受很高构造剪应力和应变的大范围破裂体中汽化为蒸汽。 汽田地震震源机制解与区域构造应变场几乎一致,并且在该区域范围内汽田地震与别的构造地震不易区分。观测中注意到地震活动与注液历史无关,这表明孔隙水压力增高与注液不可能是吉塞尔斯诱发地震的成因。 相反,所有证据都表明诱发地震与孔隙水压力及温度降低有关。形成机制有两种最大可能:其一,是裂隙排水(汽)导致局部剪应力增加所致,其二,是由稳定滑动转化为不稳定滑动(粘滑)。没有其它记载的诱发地震机制与吉塞尔斯汽田条件相符。     

Mechanism of the reservoir impounding earthquakes at Hsinfengkiang and a preliminary endeavour to discuss their cause

After the filling up of the Hsinfengkiang Reservoir Kwangtung Province, seismicity was greatly increased. The majority of earthquakes occurred in the deep water gorge close to the dam, concentrated within a northwest belt. They are usually of shallow focal depths. A strong earthquake with magnitude 6.1 took place on March 19, 1962, about two and a half years since the impounding of the reservoir.

According to the results of analysis of data from geodetic leveling and the spectra of seismic waves, the fault parameters of the main shock were determined. The fault plane solutions of 150 small earthquakes, occurring within a period of 18 months before and after the main shock were determined from the amplitudes of the first motion of P wave. The directions of the earthquake generating stress of about 2000 small earthquakes were obtained by smoothing the first motion patterns. Displacement field and stress field in the rock bodies underneath the reservoir caused by the loading of the reservoir water were calculated. Variations of the velocity ratio of the P and S waves prior to the main shock and several strong aftershocks were analysed.

In consideration of the seismicity as well as the geological background, we endeavour to discuss the cause of reservoir impounding earthquakes at Hsinfengkiang. We have the opinion that the penetration of water along fissures becomes the most important cause of the main shock of March 19, 1962 at Hsinfengkiang.     

The present status of reservoir induced seismicity investigations with special emphasis on Koyna earthquakes

The status of Reservoir Induced Seismicity (RIS) has been reviewed periodically (Rothé, 1968, 1973; Gupta and Rastogi, 1976; Simpson, 1976; Packer et al., 1979). In the present paper, the significant work carried out during the last three years on RIS is reviewed.An earthquake of magnitude occurred on November 14, 1981 in the vicinity of Aswan Lake, Egypt, 17 years after the filling started in 1964. This event occurred 4 days after the seasonal maximum in the reservoir water level and was followed by a long sequence of aftershocks. Another event of magnitude occurred in the vicinity of Aswan Lake on August 20, 1982. Results of preliminary investigations indicate that this seismic activity is reservoir induced. Recent analyses of induced seismic events at Nurek Reservoir U.S.S.R., show that the second stage of filling during August to December 1976, increasing the maximum depth from 120 m to 200 m, was accompanied by an intense burst of shallow seismic activity. An outward migration from the centre of the reservoir, possibly associated with diffusion of pore pressure, is revealed by the temporal distribution of earthquake foci. A variety of investigations including the in situ measurement of tectonic stress, pore pressure, permeability, distribution of faults, etc., in addition to monitoring seismicity, have been undertaken in the vicinity of the Monticello Reservoir, South Carolina. The largest reservoir induced earthquake is predicted not to exceed magnitude 5.The Koyna Reservoir, India, continues to be the most outstanding example of RIS. Three earthquakes of magnitude 5 occurred in September 1980. Earthquakes of magnitude 4 occur frequently in the vicinity of Koyna, the latest being on February 5, 1983. Events that occurred during the period 1967–1973 have been relocated using better procedures and are found to be much shallower and the epicentres less diffused. Location of 12 earthquakes of M_s 4.0, their foreshocks and aftershocks, that occurred during 1973–1976, composite focal mechanism solutions and related studies are consistent with the delineation of a N-S trending fault through the reservoir area. In a couple of interesting studies it has been demonstrated that earthquakes of magnitude 5.0 in the Koyna region are usually preceded by several magnitude 4 earthquakes in the preceding fortnight. Also, a rate of loading of Koyna reservoir of at least 40 ft/week appears to be a necessary, although not sufficient, condition for the occurrence of magnitude 5 earthquakes. Smooth filling/emptying appears to be the key to reduce the hazard of RIS.A map and a table of the reported cases of reservoir induced changes in seismicity through 1982 have been compiled.