737 ANNALS OF GEOPHYSICS, VOL. 51, N. 5/6, October/December 2008 Key words geoacoustic emission – deep borehole, diurnal variation – earthquake occurrence – natural electric field 1. Introduction The concepts of cracking within heteroge- neous media are the fundamental to modern models of earthquake preparation and source site formation. As indicators of micro-fractur- ing processes such as cracking, shifts along blocks boundaries, pore coalescence, changes in interstitial fluid pressure etc., spatio-temporal characteristics of acoustic emission (AE) con- tribute to validation and/or improvement of such models. This equally involves the results of laboratory-scale experiments on rock speci- mens to simulate tectonic straining processes, and geoacoustic emission (GAE) surveys (seis- moacoustics, by other term (Belyakov et al., 1996)) aimed to study natural acoustic emission of embedded rocks. The results on GAE variations recorded in a tunnel at 100 m depth (Matsushiro seismologi- cal observatory, prefecture Nagano, Japan) were represented in (Gorbatikov et al., 2001). The measurements were carried out with the help of three-component MAG-3S geophone designed at Earth Physics Institute of the Russian Acade- my of Sciences on the basis of magneto-elastic sensor (Belyakov, 2000). The output signal of such sensor is proportional to the third deriva- tive of the ground displacement, and the gain slope of MAG-3S geophone is equal to 60 Db Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity Valery Gavrilov (1), Leonid Bogomolov (2), Yulia Morozova (1) and Aleksander Storcheus (1) (1) Institute of Volcanology and Seismology, Far East Branch of the Russian Academy of Science, Petropavlovsk-Kamchatsky, Russia (2) Institute of Marine Geology and Geophysics Far East Branch, of the Russian Academy of Science, Yuzhno-Sakhalinsk, Russia Abstract Continuous geoacoustic emission (GAE) measurements were acquired using a three-component geophone placed in a borehole at a depth of near 1000 m at Petropavlovsk-Kamchatsky starting in August 2000. Using geophones consisting of magneto-elastic crystal ferromagnetic sensors, and installed at such a depth allows measurement of natural geoacoustic background with signal amplitude less than 1×10-4 m/s3 in frequency band from 3 to 1500 Hz. According to the data from a 4-year survey period the characteristics of diurnal geoacoustic variations change before every earthquake with MLH≥ 5.0 that occurs at a distance of less than 300 km from the observation point or before each earthquake with MLH≥5.5 occurring at distance R≤550 km from the observation point. The changes in GAE regime correlate with the strongest earthquakes that occurred during survey period. Measurements of the natural electromagnetic field of the Earth were carried out simultaneously with the help of an underground electric antenna. The behavior of GAE in aseismic periods appears to be related to the effect of diurnal variations of the natural electromagnetic field. Mailing address: Dr. Leonid Bogomolov, Institute of Marine Geology and Geophysics of the Russian Academy of Science, Nauki str., 1 B, 693022, Yuzhno-Sakhalinsk, Russia 738 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus tions of GAE level were also recorded (but in frequency range 30-1200 Hz) when the geoa- coustic surveys were carried out with the help of MAG -3S geophones installed in boreholes (Belyakov et al., 1999). Sensor allocation in- side deep enough boreholes (700 m depth or more) where the level of geoacoustic emission is controlled by natural crust (lithospheric) processes is a proper way to provide necessary sensitivity of GAE measurements. Unfortu- nately the results represented in (Belyakov et al., 1999) were obtained mainly by measure- ments in seismically inactive regions and there- fore do not reveal cause-and-effect relations of characteristics of geoacoustic emissions with stressed-strained state of a terrestrial crust. This paper is devoted to the results of the continuous geoacoustic measurements performed in 2000-2003 in the seismically active region of Kamchatka, Russia, with the use of MAG-3S geophone, installed inside a deep G-1 borehole at the depth of 1035 m. Some preliminary results of these measurements were described in (Gavrilov et al., 2006; Gavrilov, 2007). Hereinafter we analyze the results of geoacoustic measurements together with the data of seismicity of Southern and Central parts of Kamchatka. 2. Equipment and observation technique The station for GAE measurements was built on the base of the G-1 borehole (Lat. 53.05 N, Long. 158.63 E) located in Petropav- lovsk-Kamchatsky, a zone of deep northwest oriented fault. The depth of the borehole is 2540 m. The borehole is water-filled and cased along the whole length. The structure of the measuring system is presented in fig. 1. The geophone, which is the main measuring unit of the system, was installed at a 1035 m depth. The body of the sensor is crowded against the casing by a spring. The sensitivity of the vertical channel of the geophone (evaluated on preliminary amplifier output) is 0.15 V×s3/m, and the resonant fre- quency of the sensor is 1250 Hz. The sensitivity of the horizontal channels is 0.60 V×s3/m, the resonant frequency being 300 Hz. The output signals of geophone pre-amplifiers are trans- for a decade of frequency change. Such charac- teristics allow measurement of GAE natural background with minimal amplitude of signals (evaluated as equivalent ground displacement) less than 1×10-10 m. According to the results of (Gorbatikov et al., 2001) some increase in GAE level has been recorded in 5 cases of earth- quakes with M=3.7÷5.2 that occurred at epicen- tral distances up to 110 km. The observed varia- tions occurred nearly 12 hours before the events. These results demonstrate that geoacoustic measurements using high-sensitivity geophones (Belyakov et al., 1999) are able to provide new significant information and put some light on the earthquake preparation process. Until now most surveys have been conducted in seismi- cally inactive regions. Interesting results related to GAE measure- ments with a geophone have been revealed re- cently in works (Cuomo et al., 2000; Gregori et al., 2002; Paparo et al., 2002), where another approach has been developed with a similar purpose: to distinguish signals originated by stressed-strained crust from the background of extraneous noise. The authors (Cuomo et al., 2000; Gregori et al., 2002; Paparo et al., 2002) argued that ultrasonic acoustic emission meas- urements of the Earth surface can be a diagnos- tic tool for changes in crust stress in seismic areas. This is possible provided that the sensors are installed on bedrock outcroppings, and acoustic emission signals are recorded in the ultrasonic frequency range (10- 200 kHz) high- er than usual frequencies of GAE measure- ments. According to (Paparo et al., 2002) bed- rocks from top to bottom can play a role of a giant probe, so that AE signals recorded on the outcrop are sensitive to variations of tectonic stress and are relevant to seismic processes. Although ultrasonic signals are not able to ex- tend to large distances, anomalies of AE flow before strong earthquakes have been revealed during AE measurements in Apennines, Italy (Cuomo et al., 2000; Gregori et al., 2002; Paparo et al., 2002). The existence of distinct diurnal variations of GAE level in the frequency range 10-200 kHz and their interrelations with changes in natural electric field strength follow from data given by (Paparo et al., 2002). Diurnal varia- 739 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity Data Management Center of the Institute of Volcanology and Seismology, Far East Branch of the Russian Academy of Sciences, Petropav- lovsk-Kamchatsky, on request of the operator. 2.1. GAE level dependence on depth and borehole log Before geophone installation at a permanent depth 1035 m, measurements of GAE level were performed on various minor depths (fig. 2). At 200 m depth the level of noise in fre- quency bands 30 and 160 Hz dropped approxi- mately by 17-20 dB (7-10 times) compared to that on the original ground. Subsequent deepen- ing of the geophone up to 600 m depth was followed by smooth damping of the noise in the noted bands: it decreased by 10-13 dB com- ferred through armored cable to the main processing unit, located at the mouth of the well. The main analog unit provides additional amplification, subsequent filtering of the initial signals of each geophone channel by a third octave band pass filter, distinguishing the bands with the four central frequencies chosen: 30, 160, 560 and 1200 Hz, and, finally, the meas- urement of averaged values of output filtered and rectified signals (from each of 12 filters). The further processing of signals is realized with the help of the microprocessor controller. The functions of the microcontroller involve analog-digital converting of the input analog signals (the sampling rate is 32 Hz on the chan- nel), calculation of average values of recorded signals in one minute moving window, and, fi- nally, data saving on the disk memory. The data are transferred via telephone channel to the Fig. 1. The structure of measuring complex. 740 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus layer of 700-900 m depth and decrease sharply from 900-1035 m may be related to peculiar features of the geologic profile at the borehole site. The G-1 borehole is located in a small val- ley with rocky edges. According to the geologi- cal cross-section, the valley bed is filled up to 80 m depth by sandy-argillaceous deposits with inclusions of pyroclastic strata. The layer from 80 to 110 m depth is formed by siltstones and sandstones stratigraphically alternated by clay interlayers. Denser slates of uniform makeup occur at depths more than 110 m. According to the cross-section, the rocks lying at depths from 520 to 840 m are characterized by horizontal schistosity. There are four explicit layers on this interval of depth (fig. 2), each layer being of 50-100 m thick. Hydrothermal veins with pyrite and calcite are represented widely at depths from 650 to 740 m. This is dissimilar to only one layer occurring at depth between 100 and 520 m. The rocks lying at 820-1035 m depth pared to that at 200 m depth (in the band nearly 30 Hz). A sharp decrease of GAE level in low frequency bands at shallow depths was prede- termined by contribution of man-made noise produced at the daylight surface. The bands with central frequencies 30 and 160 Hz were close to urban frequencies. Man-made noise at frequencies 560 and 1200 Hz was less, and it did not contribute to GAE levels even at shal- low depths of 100-200 m. Correspondingly, only slight decreases in GAE levels were ob- served at higher frequency bands (560 and 1200 Hz). Significant increases of GAE levels in all frequency bands occurred in the interval of 600 to 900 m depth. The greatest increment of noise (12-13 dB) was denoted in the bands 30 and 160 Hz. Subsequent geophone embedding to 1035 m depth resulted in a decrease in noise levels for all bands. The non-trivial observation that GAE am- plitudes increase in all frequency bands in the Fig. 2. The dependence of mean GAE level at several frequency bands on depth of geophone location. GAE signals horizontal component has been plotted. 741 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity simultaneously with it or some hours later. The records of horizontal channels were never influ- enced by raining. Snowfalls exerted a stronger influence upon variations of GAE level on horizontal channels that has been confirmed by observation of the effect in 11 of 13 cases. The effect of snow fal- loff appeared as 50-100 % increase in average noise level on horizontal channels, usually starting in 12-15 hours after the beginning of snowfall and ending in 3-12 hours after snow- fall stopped. The most apparent effect of pre- cipitations occurred during strong snowfalls when the daily accumulation of snow was 30-60 mm (fig. 3). Sometimes, a delayed noise was noted even when the snowfall did not ex- ceed 10-15 mm per day. Interestingly, the data from the vertical channel never showed the ef- fect of snowfall on the GAE. The effect of precipitations on GAE measure- ments can be explained based on the following model, which appeals to the geological structure of the valley. It was already mentioned that sedi- ments of the uppermost layer (above 110 m depth) are not so consolidated as slate-like rocks situated below 110 m. Snow precipitations resuls in additional weight, and this is the reason for some extra compaction of sediments at shallow depths in territory of the valley. Simple calcula- tions show that a 30-mm-thick snow layer (the case of observed correlation with GAE) may produce additional pressure of nearly 300 Pa. The lithostatic pressure at a depth of a few hundred meters is in the order of 2-5 MPa, in which the lowest estimate is given for the 110 m depth, i.e. the bottom of the uppermost layer. The ratio of additional pressure to lithostatic one is equal ap- proximately to 10-4, and the relative change in uppermost layer strain is of the same value. We note that a typical value of Young modulus of upper sediments (occurred in this layer) is of the order of 1010 Pa, (Scheidegger, 1975). Using this value we may estimate the increment of absolute strain at 110 m depth, caused by weight of 30-mm-thick snow layer, as 3×10-8. Meanwhile, lunar-solar tides are characterized by a 10-8 change in strain components (as noted in work (Diakonov et al., 1990) for latitude 40-50 de- grees). According to a well known hypothesis (Rykunov and Smirnov, 1985; Diakonov et al., represent the single strata as well. Presumably, the rocks of enhanced stratification, with sof- tened inclusions at the bedding interface, cause enhanced dynamic and geoacoustic activity at that depth (the more contacting surfaces, the more microscopic shifts). Some results of geoa- coustic surveys in other regions (Belyakov et al., 1996, 1999) may confirm this reasoning on stratification and heterogeneity (the cause of seeming crust anelasticity). All this may ex- plain, to a certain extent, the maximal value of GAE level found at depths in the middle of in- terval 600-900 m. Apart from the schistosity interval (both upwards and downwards) the level of GAE drops. Certainly, the data on vertical dependence of GAE indicate that the layer of depth from 700 to 900 m is the most preferred for perma- nent measurements with the geophone (rather than its actual location at the 1035 m horizon). However these data (summarized in the fig. 2) were obtained after the geophone installed at 1035 m depth. Thereafter, because of technical reasons, we had no opportunity to change the position of the geophone. The main methodical result of GAE amplitude vertical scanning (fig. 2) is direct proof that geophone installation in- side the borehole at near 1000 m depth greatly reduces low frequency exogenous noise (30 and 160 Hz bands). The level of noise reduced by 40 dB (two orders of value) at 30 Hz frequency, and also reduced by 17 dB (nearly order) at 160 Hz frequency, both cases in comparison with level of noise on original ground. 2.2. Influence of meteorological factors Comparisons of GAE variations with changes in atmospheric pressure, air tempera- ture and speed of wind revealed no appreciable correlation. A possible influence of rain on measurements on the vertical channel was not- ed in two of 24 cases. In the first case, the rainfall was about 20 mm in a day, in the second - about 50 mm in a day. The effect of precipitations appeared in the form of 20-30 % bay-like increases in average level of noise which were delayed of 12-15 hours after the beginning of rain and terminated 742 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus face) is very sensitive to contact state changes. Rain precipitations have practically no effect on the results of GAE measurements, because the water goes away quickly. Hence, minor incre- ment of weight after raining together with the short term action of water weight may induce only weak perturbation in comparison with the case of snow precipitation. So, the change in GAE level on horizontal channel is expected to correlate to temporal variations of normally acting force (horizontal pressure). The correlation has a lag of 12-15 hours after accumulation of snow precipitations. Actually such a lag is a transition time during which pressure change transforms to plastic strain wave and propagates along the casing tube towards horizon of geophone location. The vertical channel of GAE is not affected by this mechanism, since receiving condition for GAE shearing component is not so sensitive. According to the mentioned observations and their interpreting, the GAE data from horizontal channels have been distorted by exogenous inter- 1990), the tides have potential to trigger some fault slips, as well as geoacoustic emissions, since the rate of tidal stress is high enough with respect to a tectonic stress. The growth of addi- tional strain during an intensive snowfall is still more rapid than that in the case of tides. Actually, a snow accumulation of 30-mm-thick layer may occur in few hours, and the strain rate reaches the value 3×10-12 1/s. The rate of tidal strains was estimated as 1-3 × 10-12 1/s. So, the change in the uppermost layer strain, caused by snow falloff, may have some influence on inelastic processes (similarly to the case of solar-lunar tides). The effect is probably amplified due to a structural factor. Edges of the valley are steeply inclined, although not vertical in the cross-sec- tion. The compaction of the upper sedimentary layer inside a limited area may cause a portion of the horizontal stress to act normally to the surface of the casing tube. The contact of the casing and the adjacent rock becomes more or less strong. The receiving condition of GAE of horizontal component (normal to contact sur- Fig. 3a-c. Episodic influence of snow falloff on results of GAE measurements: a) daily accumulation of snow, b) temporal dependence of vertical component GAE level, 30 Hz central frequency; c) the same for horizontal component. 743 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity where the lines above variables or expressions note their averaging over time window. As result of the solving system of the equa- tions (3.2) one can receive: a fg f g g = − ⋅ −1 2 ; b f fg g g = − ⋅ −1 2 ; (3.3) This means that the coefficient a is the aver- aged amplitude of DGAED, whereas b is the daily average level of GAE. Coefficients a and b concern to data on entire window (all n counts of function f (ti )). We have selected a two-day window (n = 2880) for computations by the outlined algo- rithm. Transition intervals of 30 min length close to the estimated times of Sun rising and a sunset have been excluded from the analysis. The exceeding of the threshold value ac of 500 relative units (RU) by the average amplitude of DGAED has been considered as a criterion for the presence of GAE diurnal distribution pro- vided that the maximal value of DGAED am- plitude was equal to amax of 2804 RU. The result of computations is that the charac- teristic diurnal variations have been distinguished in GAE time series during approximately 60% of observation time. An example of such variations in temporal intervals during August 13-16, 2001 and September 3-7, 2001 is presented on fig. 4a and 4b (Greenwich time is hereafter used). Tran- sition times from the minimum values of emis- sion level to the maximal one and vice versa are about 5-10 minutes. The moments of transitions correspond to the times of terminator line cross- ing (times of sunset and sunrise) at the observa- tion point. The most explicit diurnal variations (day time minimum and nocturnal maximum) have been revealed on GAE series recorded by channels for vertical components, with the cen- tral frequencies of 30 Hz and 160 Hz. 3.1. The changes in GAE characteristics related to the preparation of earthquakes As a rule, considerable perturbations of the diurnal GAE distributions (DGAED) have been observed before earthquakes MLH≥ 5.0 in the ferences to a greater extent than that from the vertical channels. For this reason, only measure- ments of the vertical component of GAE signals (i.e. the data from the vertical channel) have been used below in this paper, as outlined below. 3. Results of continuous GAE surveys at the 1035 m depth Continuous GAE measurements started in August 2000. In this section we consider peculiar features of temporal dependences of geoacoustic emission obtained by continuous surveys during August 5, 2000 to July 31, 2003. During primary data analysis we focused on manifestations of characteristic diurnal variations of averaged GAE levels. Characteristic variations constitute diurnal geoacoustic emission distributions (DGAED) of 24 hours periodicity which usually emerge in aseismic periods (Gavrilov et al., 2006). An especially developed PC program has been used to establish the presence or absence of DGAED in GAE time series from August 5, 2000 to July 31, 2003. Sometimes a simple visualization of these time series (plotted out- put) can reveal diurnal variations. The pro- gram of the GAE time series processing is based on «recognition» of the GAE variations similar to logistic function of sunrise/sunset g (t), so-called meander (trial rectangular incre- ment). The function g (t) is equal to 1 (unity) when the Sun is above the horizon, and it is -1 after a sunset. The algorithm of recognition in- volves the selection of such coefficients of re- gression a and b, to minimize the regression function (below) in the given time window R f t ag t bi i i = −∑( ( ) ( )– ) ,2 (3.1) where f (ti) - current value of GAE level, t – time. Coefficients a and b can be determined from system of the equations: ∂ ∂ = − + = R a ag fg bg2 2 2 02 ; ∂ ∂ = − + = R b b f ag2 2 2 0; (3.2) 744 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus Prior to the seismic swarm in October 2002, which occurred 170 km to the south of Pet- ropavlovsk-Kamchatsky, the DGAED disap- peared seven days prior to the first shock (fig. 6). Until this moment, a steady diurnal distribu- tion was observed during about three months. Being typical for aseismic periods, GAE day - night variations restored a day after the stron- gest (MLH =5.7) earthquake of the swarm. According to the catalog of earthquakes of the Geophysical Service of the Russian Acade- my of Science, 56 earthquakes MLH≥ 5.0 in- area of the Kamchatkan subduction zone. Typi- cal examples of such perturbations prior to the earthquakes are presented in figs. 4c and fig. 5. Figure 4c shows the disturbance of DGAED before an MLH=5.2 earthquake, which occurred on September 1, 2001, 120 km to the northeast from the observation point. In this case the di- urnal GAE distributions, which appear to be normal during August 12-16, 2001, are con- spicuously distorted (almost disappeared) dur- ing August 17-31, 2001. The DGAED return to the normal form soon after the earthquake. Fig. 4a-c. A plot of characteristic diurnal variations of mean geoacoustic emission level (for Z component of GAE signals in a band with central frequency of 160 Hz). The Greenwich time is noted on time axis. Fig. 5. Examples of loss of diurnal GAE distributions (DGAED) prior to the occurrence of earthquakes. 745 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity Fig. 6. Mapping of epicenters of earthquakes occurred during the period January 1, 2001- July 31, 2003 in zone R≤300 km, centered at Petropavlovsk -Kamchatsky (events of MLH≥5.0 magnitude), and in zone R≤550 km (events of MLH≥5.5 magnitude). 746 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus of disorder of GAE diurnal distribution was fol- lowed by an earthquake (as shown in fig. 7). Meanwhile, the disturbances of DGAED (full disappearance or degradation) occurred before all 15 earthquakes of magnitude МLH≥5.0 in- side circular zone around point of observation, of 300 km radius. Also, such disorders occurred before all 21 strong earthquakes, МLH≥5.5 whose epicenters located at distances less than 550 km from the geophone point. The same is valid for the case of the strongest MW =7.3 earthquake that occurred on November 17, 2002 at 1000 km epicentral distance. So, the disorder of diurnal variation may be considered a necessary but not sufficient condi- tion of a forthcoming earthquake. Our measure- ments at the single point give no result on pre- diction of place of shock. But the most impor- tant consequence of the observations is that during a period of regular diurnal variations in GAE (lasting 7 days or more) no earthquake is expected as long as a disorder of the DGAED is to occur. Below we consider statistical aspects of a relationship between times of strong earth- quakes occurrence and cases of the DGAED disorder/restoration. cluding aftershocks of strong earthquakes oc- curred within 550 km in distance from Pet- ropavlovsk-Kamchatsky during August 1, 2000 – July 31, 2003. A subsample of 36 seismic events was selected and used for analysis, tak- ing into account 15 earthquakes of magnitude МLH≥5.0, occurred in near-field zone (inside a circle of R≤300 km radius) and 21 earthquakes of major magnitude (МLH≥5.5), with epicenters located in circular zone R≤550 km. Parameters of these earthquakes are given in table I and locations of epicenters are shown on fig. 6. The strongest MW =7.3 earthquake during the last five years that occurred on November 17, 2002, at 459 km depth is also presented in table I. The epicenter of the earthquake was located in the Okhotsk Sea at 1000 km from the point of observation. Figure 7 compares periods of DGAED pres- ence or absence with exact times of occurred earthquakes. It is shown that earthquakes most- ly occur during periods without characteristic diurnal variations. This includes an example of rather short interval (24 hours) without charac- teristic diurnal variations, for which earthquake of August 2, 2001 occurred. Not every episode Fig. 7. The periods of presence/absence of diurnal variations of GAE versus times of major earthquakes occur- rence: 1- presence of characteristic DGAED, 2-data gaps. The earthquakes with МLH≥5.0 occurred at epicentral distances of R≤300 km from the base station and the earthquakes with МLH≥5.5 occurred at epicentral distances of R≤550 km. The earthquake of November 17, 2003, Mw = 7.3 occurred at near1000 km epicentral distance. 747 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity Table I. Earthquake’s parameters, date January 1, 2001- July 31, 2003, location R≤300 km, МLH≥5.0; location R≤550km, МLH≥5.5; location R≤1000 km, МW>7.0) Note: H- earthquake depth; R- distance from epicenter to borehole G-1 Earthquake No Date MLH Lat.N Long.E H, km R, km 010207 2001:02:07 6.4 52,28 153,66 476 346 010802 2001:08:02 6.2 56,21 164,05 25 495 010901 2001:09:01 5.2 53,92 159,75 134 120 010917 2001:09:17 5.3 52,84 159,98 41 191 011007 2001:10:07 5.2 52,39 160,67 2 153 011008a 2001:10:08 6,1 52,62 160,46 31 129 011008b 2001:10:08 6.3 52,65 160,49 24 130 011008c 2001:10:08 5.1 52,5 160,59 15 142 011009 2001:10:09 5.0 52,43 160,59 18 146 011010a 2001:10:10 5.5 52,51 160,57 17 141 011010b 2001:10:10 5.2 52,46 160,72 8 152 011103 2001:11:03 5.5 55,93 161,35 104 366 020106 2002:01:06 5.9 48,67 155,09 40 543 020128 2002:01:28 6.0 49,29 155,98 21 454 020215 2002:02:15 5.0 52.12 159,89 18 130 020315 2002:03:15 5.5 49,44 155,96 21 439 020426 2002:04:26 5.9 53,36 160,99 57 160 020503 2002:05:03 5.0 52,49 160,79 20 155 020508a 2002:05:08 5.6 52,22 160,44 32 150 020508b 2002:05:08 5.9 53,73 160,93 35 170 021008 2002:10:08 5.3 52,72 160,30 33 115 021016 2002:10:16 5.7 51,66 157,68 108 165 2002:11:17 МW= 7.3 47.82 146.21 459 1050 021218 2002:12:18 5.0 52,91 159,82 40 79 030120 2003:01:20 5.7 49,06 155,88 54 480 030221 2003:02:21 5.3 55,45 159,79 349 279 030315 2003:03:15 5.8 52,15 160,66 4 166 030317а 2003:03:17 5.0 52,25 160,58 13 155 030317b 2003:03:17 5.1 52,09 160,72 13 174 030317c 2003:03:17 5.5 52,26 160,54 32 153 030318 2003:03:18 5.4 52,23 160,61 15 159 030319 2003:03:19 5.9 52,16 160,85 48 176 030325 2003:03:25 5.0 52,02 160,70 40 178 030423 2003:04:23 6.1 55,98 163,44 20 450 030424 2003:04:24 6.0 48,76 155,21 42 531 030529 2003:05:29 5.7 50,65 157,53 40 274 030616 2003:06:16 6.7 55,3 160,74 190 263 748 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus P q n k q qn k( , , ) ,–= Ckn k(1− ) (3.4) Where Ckn is a binomial coefficient which is determined by the formula Ckn = n!/ ( k!(n-k)! ) To use (4) we associated the number n of tests for GAE data with the number of strong earthquakes which occurred independently of one another during the time of observations. Some seismic swarms were recorded in the pe- riod October, 2001–March, 2003. Since shocks in a swarm were not independent events, we considered each swarm as a test (in the same manner as an ordinary earthquake). The total number of independent seismic events (tests) is n=26. Only one earthquake occurred when the characteristic diurnal variations were present (the date is June, 16 2003, see fig. 7). So, k = 1 and the number of successful events is (n – k) = 25. A value of probability q may be estimated as a ratio of time without GAE diurnal varia- tions to total time of observations. In the case considered q = 0,4. The binomial distribution gives following numerical estimate of probability for noted val- ue of parameters q,n,k: Р = 1.8 ⋅ 10-9. Such value of probability is very small compared to unity, and in comparison with the maximum of expres- sion (3.4). It is well known that the expression of binomial distribution (3.4) reach maximum when q = (n-k)/n. In our case the maximum of (3.4) is approximately equal to Рmax ≈ 0.18. In a conditional case when all 37 earth- quakes are independent, n=37, and k=1 the value of probability Р for the same q is of order 10-14. So, the hypothesis is not valid that a dis- order of characteristic diurnal variations and final phase of earthquake preparation is a ran- dom coincidence. One can deduce from the negligible probability of co-incidental correla- tion between GAE variations mode and times of strong seismic events that onset of GAE di- urnal variations damping regularly precedes the possible occurrence of an earthquake. 4. Discussion of a possible origin To study the correlation between DGAED and variations of averaged strength of natural electromagnetic (EM) field we carried out syn- Temporal intervals from the beginning of the «usual» DGAED perturbation until the time of the earthquake MLH≥ 5.0 are in limited to 1÷25 days. The relationship between noted temporal intervals and parameters of the earth- quakes is still debated. This also concerns the recovery time for diurnal distribution after the earthquake occurred. Restoration of DGAED following the earthquakes takes several days in about 30 % of all cases, but in some cases the recovery time was much greater. It is worth mentioning that we considered several of limits on magnitude МLH and radius of epicentral zone R. We attempted to separate stronger earthquakes (occurring times to be com- pared with presence/absence of diurnal varia- tions) and the flow of weak and moderate seis- mic events. On the other hand, we need a suffi- cient accumulation of strong events for the analysis. Enumeration of possibilities resulted in the optimal choice of above limits of magnitude МLH and radius R. We selected two limiting val- ues of МLH (5.0 and 5.5) and two corresponding values of R (300 and 550 km) for the sake of strong validity of the results. Our optimal choice of МLH and R limits allowed us to obtain the above relationship between GAE diurnal varia- tions and the seismic process in the studied re- gion. Another selection of these limits may modify the result represented in fig. 7. For ex- ample, reducing the МLH limit under condition of the same value of radius R entails a growth of the number of earthquakes that occurred on the background of GAE diurnal variations. In one case an earthquake occurred during a period of diurnal distribution restoration, when the time of DGAED disappearance ended. This does not contradict the fact that the disappear- ance of DGAED preceded the occurrence of all 37 earthquakes. But we prefer to consider this debatable case as unsuccessful event. We estimated the probability of coincidence of anomalies of GAE observations (disorder of characteristic diurnal variations) and strong earthquake occurrence with the help of proba- bility binomial distribution technique. The bi- nomial distribution describes the probability P(q,n,k) of occurrence of k unsuccessful events in n tests on conditions that the probability of successful result of individual test is (1-q): 749 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity phasizing that we intended to involve just elec- tromagnetic field because of following a priori reason. It is well-known that a temporal de- pendence of natural electric field (vertical com- ponent EN, in particular) has diurnal variations with nocturnal maxima and day time minima. This notion is based on studies of characteris- tics of radio noise and electromagnetic waves propagation in the frequency band 30 Hz - 30 kHz (Osinin, 1982; Remizov, 1985). Lightning discharges in remote world thunderstorms cent- ers are the main sources of electromagnetic waves in noted frequency band at North- East- ern part of territory of Russia (Remizov, 1985). chronous measurements of GAE and EM varia- tions. An underground electric antenna was used for measurements of the vertical compo- nent of the natural electric field, EN. The per- ceptibility of the antenna whas reached with wire elements of more than 1 km length: a free conductor of multicore cable of the geophone and borehole casing. The parameters of electro- magnetic measuring channels were identical with geoacoustical ones (the only difference is in the source of signals), measurements were carried out in the same frequency bands. The measurements started in May 2003 at the same point (G-1 borehole). It is worth em- Fig. 8a-c. Results of simultaneous measurements of geoacoustic and electromagnetic temporal variations: a) GAE level measurements at 160 Hz frequency band (Z component); b) EM field level measurements at the same frequency band; c) cross correlation coefficient ρ. The vertical line indicates earthquake moment. 750 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus periods when DGAED were steady. Meanwhile such periods coincided with aseismic time gaps. During these periods the cross correlation coefficient ρ was of order of value ρ = 0.80 ÷ 0.99. So, for instance ρ = 0.91 ÷ 0.99 for interval of April 9-12, 2004, fig. 8; ρ = 0.92 ÷ 0.93 for interval of June 15-17, 2003, fig. 9. But, before strong earthquakes (a day or more) and during relaxation periods after the events the variations of GAE and of electro- magnetic field became dissimilar: ρ = 0.53 ÷ 0.89 for interval of April 13-19, 2004, fig. 8; ρ = 0.001 ÷ 0.34 for interval of May 28 - June 09, 2003, fig. 9. It should be noted especially, that the strong correlation between results of GAE and natural The day-night variations of ultra low frequency (ULF) electromagnetic field are pre-determined by daylight sharp deterioration of radio wave propagation condition in a waveguide between ionosphere and the Earth surface primarily due to lowering of ionosphere height resulted from daylight occurrence of D-layer at near 80 km height. In Kamchatka, ULF electromagnetic emission with explicit diurnal variations is caused mainly by lightning discharges in a re- mote world thunderstorm center, located on azimuth 190-255 degrees. The results of simultaneous geoacoustic and electromagnetic measurements demonstrated that the variations of GAE level and electro- magnetic ones were practically identical during Fig. 9a-c. An example of resumption of day-night variations of mean GAE level: a,b,c- the same as on the fig. 8. 751 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity wires of the multicore cable to be sure that out- put signals of GAE depended on sensor only. This differs from features of induced noise which is sensitive to position of conductor (dis- tance from wire to core screen). Some closely related results of other works confirm the physical rather than instrumental origin of correlated diurnal changes in GAE and EM signals. For instance, a high degree of correlation of diurnal variations of GAE and electromagnetic field during the aseismic peri- od is correlated with the result of (Paparo et al., 2002) where diurnal variations of AE level on bedrock outcroppings with maxima at local nighttime (similarly to GAE, see interval of April 9-12, 2004, fig. 8) were described. The synchronization of increase/decrease of AE and of natural telluric electromagnetic field have been noted. Both temporal dependences can be characterized with a meander line. The similarity of daily distributions of GAE and AE which differ by frequency range and depths of sources speak in favour of the fact that this phenomenon is a manifestation of behavior of the loaded geological medium. No particular factor like bias of electric induction over long cable or thermal stress near bedrocks surface due to nighttime cooling (Paparo et al., 2002) is adequate for consistent explanation of those synchronized variations. The comparison of results obtained by bore- hole GAE measurements with the results of laboratory experiments with loaded specimens exposed to electromagnetic field (Chelidze et al., 2002; Sobolev and Ponomarev, 2003; Il’ichev et al., 2003; Bogomolov et al., 2004; Zakupin et al., 2006a) has become the theory of a unified role of electric field for AE and GAE variations. According to the experiments re- ferred above electromagnetic action over rock specimens which are tested by press driven compressive loads accompanied by appreciable growth of AE activity. The main feature of the revealed effect of AE electrostimulation is the temporary growth of averaged level of AE ac- tivity after short impact or during a session with enhanced electric strength, provided that the main compressive load is constant and its value is in the range from 0.7 up to 0.95 of disruption. This is the so-called AE response to electric EM field measurements (in particular, synchro- nous diurnal variations in both signals) cannot be explained by electromagnetic noise induced by the environment. Non-instrumental origin of the observed effect may be verified by pre- sented facts and considerations. Usually, elec- tromagnetic noise collected by sensor input and/or transmission line (long cable), particu- larly the influence of the neighboring measur- ing line, is a constant, when the location of sensors and cables is fixed. But the results of simultaneous geoacoustic and electromagnetic measurements represented by fig. 8 demon- strate the change in GAE record, meanwhile the plot of natural electric field (vertical compo- nent) is without change. Diurnal variations of GAE disappeared rapidly (fig. 8a), and the cor- relation coefficient between GAE and EME series reduced. The dissimilar change in GAE was observed shortly before earthquake, but this aspect plays no role for the point of that how to check up the absence of noise induced by neighboring channel. Another example is given by fig. 9, where the period from May 28 to June 11, 2003 is without GAE diurnal varia- tions, thereafter the variations arise in interval June 12-17, 2003 and their amplitude increases gradually. But diurnal variations in the natural electromagnetic field (fig. 9b) are of the same amplitude during the period of May 28-June 17, 2003. This fact may be considered a solid argu- ment that variations in output signals on GAE measuring channels do not interference induced by external electromagnetic field via sensor details or cable connecting the geophone with recording and processing units located at origi- nal ground. In the contrary case the amplitude of GAE diurnal variations should be steady throughout the period of observations. Moreover, the temporal plots of horizontal and vertical component of GAE have quite a different form, although the frequency band is the same (Gavrilov et al., 2006). The construc- tive elements of both channels are identical. So the characteristics (such as shunt capacitance and mutual inductance to circuit of electromag- netic channel) are close, and they predetermine noise occurrence. We checked this before geophone installa- tion. Also, we reconnected geophone output to 752 V. Gavrilov, L. Bogomolov, Y. Morozova and A. Storcheus strained geological medium. Deepening of the geophone at 1000 m depth in borehole at Petropavlovsk-Kamchatsky reduced the level of anthropogenic noise by two orders of magni- tude in the band of 30 Hz central frequency, and reduced it almost by an order of magnitude in the 160 Hz band. Furthermore, increasing the depth of measurements has practically elimi- nated the bias of meteorological factors on measurements of the vertical component of GAE signals. New GAE measurements have revealed diurnal temporal distribution of GAE level in bands with central frequencies 30 and 160 Hz. Diurnal distribution of GAE is presum- ably related to daily variations of the level of electromagnetic field strength in Earth’s crust. The loss or damping of geoacoustic day-night variations (provided continued diurnal varia- tions in natural electric field) appears to indi- cate changes in stressed-strained state of the geological medium which allow a considerable growth of probability of strong earthquake oc- currence. A generalized sequence of obtained results is proposed as a precursor of «short-term aseismicity»; an indicator that no strong earth- quakes are likely to occur in a circular (300 km radius) area around the observational point, during a period of robust presence of GAE di- urnal variations, as long as a disorder of GAE mode terminates such period. This important principle (precursory statement) may be used for recognition and/or reliable predictions of temporal intervals, when the probability of strong earthquake occurrence in a limited area is very low (at least, for Kamchatkan region and North Kuriles islands). Acknowledgements This research was partially supported by the RFBR (grants 06-06-05-96071 and 07-05- 00687) and FEB RAS (the project 06-III-А-08-327). REFEREncES Belyakov, A.S., A.G. GamBurtsev, V.S. lavrov, A.V. Nikolaev and N.K. Privalovskii (1996): Vibration Action on Rocks and Their Induced Seismic Emission, power action, EPI, (Zakupin et al., 2006a; 2006b). Heuristically, such acoustic emission response to EPI seems to look similar to the ef- fect of nighttime GAE level growth correlated to nocturnal enhancement of electric field level (the more value of natural electric field strength the more level of GAE in the same frequency band). Our approach to unified consideration of laboratory and full-scale survey results is based on the similarity of stressed state of embedded rocks and laboratory tested ones. Actually, AE responses were recorded on specimens when the main compression stress was in the range 20 - 100 MPa. A tectonic stress at nearly one km depth is of the same order of magnitude (Kro- potkin et al., 1987). But a difficult point for the hypothesis of electric stimulation of nighttime GAE level (resulting in the modulation) arises during quantitative estimates of strength of natural electric field and that applied in laboratories. Quantitative estimates to analyze the observed modulation of GAE level by natural electro- magnetic emission of the Earth are as follows. According to reference (Remizov, 1985), the maximal value of the electric component of the EM field measured on Earth surface, over the whole frequency band may be evaluated as 0.45÷1.0 mV/m. The damping of electric com- ponent at 1 km depth, calculated for 160 Hz and the soil of moderate moisture (the conductivity is about σ = 0.01 (Ohm×m)-1) is near 9 db. Hence, the amplitude of electric field strength at the depth of geophone location can not ex- ceed 0.16÷0.36 mV/m. Measurements in other boreholes and in other regions are required to establish that correlated diurnal changes in GAE are a general effect (although the mecha- nism of electric- geoacoustic conversion is not clear completely). Alternatively they may be caused by specific oversensitivity of GAE sources located in the vicinity of the geophone at the point of G-1 borehole. 5. Summary Continuous geoacoustic measurements in deep boreholes with high-sensitivity geophones are an effective method of monitoring stressed- 753 Variations in geoacoustic emissions in a deep borehole and its correlation with seismicity GreGori, G.P., G. PaParo, u. CoPPa and i. marsoN (2002): Acoustic Emission in geophysics: a reminder about the methods of analysis, Bull. Geophys. Teor. Appl., 43 (1- 2), 157-172. il’iChev, P.V., A.V. alad’ev, L.M. BoGomolov, V.N. syChev et al. (2003): Parameters of acoustic emission signals initiated by electric action over loaded speci- mens, in Geodynamics and geoenvironmental prob- lems of high- mountain regions, edited by S.V. GoldiN and yu.G. leoNov (Printhouse, Moscow-Bishkek), pp. 286-303. kroPotkiN, P.N., v.N. efremov and v.m. makeev (1987): Stressed state of Earth Crust and Geodynamics, Geo- tectonics, (1), 3-24. osiNiN, V.F. (1982): Radio noise from natural sources in the East of the USSR (nauka, Moscow), pp. 161 (in Rus- sian). PaParo, G., G.P. GreGori, U. CoPPa, R. de ritis and A. taloNi (2002): Acoustic Emission (AE) as a diagnos- tics tool in geophysics, Annals of Geophysics, 45 (2), 401-416. remizov, L.T. (1985): natural radio-interference (Nauka, Moscow), pp. 198 (in Russian). rykuNov, a.l. and v.B. smirNov(1985): Variations of the seismicity under effect lunar-solar and tidal defor- mations, Fizika Zemli (Solid Earth Phys.), 21 (1), 97-103. SCheideGGer, A.E. (1975): Physical aspects of natural ca- tastrophes, (Elsevier Scientific, Amsterdam), pp. 289. soBolev, G.A. and A.V. PoNomarev (2003): Physics of earthquakes and precursor, (Nauka, Moscow), pp. 270. zakuPiN, A.S., A.A. avaGimov and L.M. BoGomolov (2006a): Responses of Acoustic Emission in Geomate- rials to the Action of Electric Pulses under Various Values of the Compressive Load, Physics of the Solid Earth, (10), 43-50. zakuPiN, A.S., A.V. alad’ev, L.M. BoGomolov, V.N. syChev et al. (2006b): Interrelation of electric polari- zation and acoustic emission of terrestrial materials specimens under conditions of uniaxial compression, Volcanology and seismology, (6), 22-33. Fizika Zemli (Solid Earth Phys.), 32 (2), 68-74. Belyakov, A.S., A.G. GamBurtsev, V.S. lavrov, A.V. Nikolaev and L.L. khudziNskii (1999): Underground Acoustic Noise and Its Relation to Tidal Deformations, Ibid, 35 (12), 1002-1009. Belyakov, A.S.(2000): Magnetoelastic acoustic geophones for geophysical research, Seismic equipment (UIPE RAS, Moscow), Book 33, 27-45. BoGomolov, L.M., P.V. il’iChev, v.a. Novikov, v.i. okuNev, v.N. syChev, A.S. zakuPiN (2004): Acoustic emission response of rocks to electric power action as seismic- electric effect manifestation, Annals of Geophysics, 47 (1), 65-72. Chelidze, t., N. Novaramashvili, m. devidze, Z. tChe- lidze, v. Chikhladze and t. matCharashvili (2002): Laboratory study of electromagnetic initiation of slip, Annals of Geophysics, 45 (5), 587-598. Chelidze, t., t. matCharashvili, J. GoGiashvili, o. lur- smaNashvili and m. devidze (2005): Phase synchro- nization of slip in laboratory slider system, nonlinear Processes Geophys., 12, 1-8. Cuomo, V., V. laPeNNa, M. maCChiato, I. marsoN, G. PaParo, d. Patella and s. PisCitelli (2000): Geolec- trical and acoustical anomalous signals compared with seismicity in a test side of southern Apennines (Italy), Phys. chem. Earth. (A), 25 (3), 255-261. DiakoNov, B.P., B. S. Karryev, O.B. KhavroshkiN, a.v. Nikolaev, l.N. rykuNov, r.r. seroGlasov, a.k. troJaNov and v.v. tsiPlakov (1990): Manifestation of earth deformation processes by high-frequency seis- mic noise characteristics, Phys. Earth and Planet. In- ter., 63, 151-162. Gavrilov, V. A., yu.V. morozova and A.V. storCheus (2006): Variations in the Level of Geoacoustic Emis- sion in Deep Well G-1, Kamchatka and Their Relation to Seismicity, Volcanology and seismology, (1), 52-67. Gavrilov, v. a. (2007): Physical Causes of Diurnal Varia- tions in the Geoacoustic Emission Level, Doklady RAN, Earth Sciences, 414, (4), 638-641. GorBatikov, a.v., o.a. molChaNov, m. hayakava, s. uyeda, et al. (2001): Acoustic Emission Before and After Earthquakes, Volcanology and seismology, (4), 66-78.