GIANNOPOULOS_corrected_Layout 6 Temporal variation of shear-wave splitting parameters before and after the 2008 Movri Mountain earthquake in northwest Peloponnese (Greece) Dimitrios Giannopoulos1, Efthimios Sokos1,*, Konstantinos I. Konstantinou2, Athanasios Lois1, G-Akis Tselentis1 1 University of Patras, Seismological Laboratory, Patras, Greece 2 National Central University, Department of Earth Sciences, Jhongli, Taiwan ANNALS OF GEOPHYSICS, 55, 5, 2012; doi: 10.4401/ag-5586 ABSTRACT On June 8, 2008, at 12:25 GMT, a MW 6.4 earthquake, the Movri Mountain earthquake, occurred in the area of northwest Peloponnese, western Greece. The epicenter was located in the municipality of Movri, 35 km southwest of Patras. For this study, a crustal anisotropy analysis was performed in the epicentral area of the Movri Mountain earthquake. Specifically, the shear-wave splitting phenomenon and its temporal evolution in relation to the Movri Mountain earthquake was studied, using the cross-correlation method. The data analysis revealed the presence of shear-wave splitting in the study area. Both before and after the Movri Mountain earthquake, the polarization directions of the fast component of the shear waves followed a general NNW-SSE direction. The observed mean fast polarization direction was not consistent with the regional stress field, which showed a general E-W direction of the maximum horizontal compressive stress. The differences between the estimated fast polarization directions and the properties of the regional stress field suggest the presence of a local stress field in the area around the fault. An increase in time delays was observed soon after the Movri Mountain earthquake. The average value of the delay times before the earthquake was ca. 18 ±2.6 ms, while after the earthquake this was ca. 40 ±4.6 ms. This increase in the time delay indicates changes in the crustal properties, which were possibly caused by variations in the pre- existing micro-crack system characteristics related to the Movri Mountain earthquake, and the possible involvement of over-pressured fluids. 1. Introduction Shear-wave splitting is a phenomenon in which shear- waves are separated into two components with different po- larization directions and propagation velocities. This can occur during S-wave propagation through an anisotropic medium [Crampin and Chastin 2003, Crampin and Peacock 2005]. The two splitting parameters that can be measured through shear-wave data processing are the polarization di- rection { of the fast component of the shear waves, and the time delay dt between the two components. Various models have been proposed to interpret the observed seismic anisotropy in the crust. The most widely accepted physical model, which is known as the extensive dilatancy anisotropy model [Crampin 1978, 1993, Crampin et al. 1984b], explains the principal cause of the shear-wave splitting phenomenon as S-wave propagation through stress-aligned, fluid-saturated micro-cracks with an orientation parallel/ sub-parallel to the direction of maximum horizontal compression. This model has since been modified, and led to the anisotropic poro-elas- ticity model [Crampin and Zatsepin 1997, Zatsepin and Crampin 1997], which explains the way in which such micro- cracks evolve in response to changing conditions in perme- able rock. Several studies on local earthquakes worldwide have revealed shear-wave splitting in various geological set- tings that indicate anisotropic media [Kaneshima 1990, Crampin and Lovell 1991, Gao et al. 1998, Liu et al. 2008]. Changes in the parameters of the shear-wave splitting phenomenon have been observed worldwide in relation to earthquakes that have indicated changes in the anisotropic characteristics of the medium and the stress field [Gao and Crampin 2004, Crampin and Peacock 2008, Crampin and Gao 2012]. These variations in the shear-wave splitting pa- rameters are complicated and have often been the subject of debate. Specifically, variations in shear-wave splitting time de- lays before large earthquakes have given rise to discussion [Crampin and Peacock 2008]. For instance, Liu et al. [2004] analyzed shear-wave splitting of foreshocks and aftershocks of the September 20, 1999, MW 7.6 Chi-Chi earthquake in Taiwan. No evidence was found for temporal changes in time delays over the periods before and after the mainshock. The authors suggested that the behavior in the dataset was affected by spatial variations of seismic anisotropy rather than by temporal variations. On the other hand, in com- menting on a previous study, Crampin and Gao [2005] claimed that stress fluctuates during aftershock sequences Article history Received February 23, 2012; accepted June 30, 2012. Subject classification: Earthquake faults: properties and evolution, Tomography and anisotropy, Waves and wave analysis, Stress, Tectonics. 1027 RESEARCH ARTICLES and systematic changes in time delays are not expected and have been observed. However, according to Crampin and Gao [2005], temporal variations of time delays were meas- ured by Liu et al. [2004] before the Chi-Chi earthquake. These variations showed the same characteristic behavior, with comparable magnitude/ time-delay duration similari- ties as in other studies [e.g., Crampin and Gao [2005]. Liu et al. [2005] maintained their assessment that the results did not show systematic temporal evolution of anisotropy parame- ters that can be used to forecast impending large earth- quakes. In another case, Aster et al. [1990] measured time delays on the same seismograms from the Anza Seismic Network that Peacock et al. [1988] and Crampin et al. [1990] used to show systematic temporal changes before the 1986 M6 North Palm Springs earthquake in Southern California, USA. These authors showed that there were no temporal changes, using automatic techniques to measure time delays. Crampin et al. [1991] argued that some of the automatic time-delay measurements of Aster et al. [1990] were in error by a factor of up to 3 (200%). They claimed that the inade- quacy of all of the automatic measurements invalidated the conclusions of Aster et al. [1990], and the estimates of the time delays between the split shear waves were not reliable. Other cases in the literature that have provided observations of changes in splitting parameters related to earthquakes have been debated, and those under considerable discussion are Seher and Main [2004], Peng and Ben-Zion [2005], and Munson et al. [1995]. The tectonic and geological setting of western Greece is characterized by complex structures. The tectonic frame- work is dominated by the large dextral strike-slip fault off the coast of Kefallonia Island, where the change between conti- nent–continent collision in the north and ocean–continent subduction in the south occurs [Underhill 1989, Sachpazi et al. 2000] (Figure 1). During the Eocene, northwestern Pelo- ponnese was affected by the Alpine collision, which led to the formation of the Hellenic mountain range [Doutsos et al. 1993]. Mesozoic and Early Cenozoic carbonate rock propa- gated upwards and westwards along N-S-striking and E-dip- GIANNOPOULOS ET AL. 1028 Figure 1. Geological map of northwest Peloponnese showing the active faults, the Movri Mountain earthquake epicenter (star), the Riolos (RLS) station (triangle), and the reactivated faults during the 2008 Movri Mountain earthquake, according to Koukouvelas et al. [2009]. Figure modified from Koukouvelas et al. [2009]. 1029 ping thrust faults [Xypolias and Doutsos 2000]. From the early Pliocene to the present, an extensional stress field has progressively prevailed in the previous compressional tec- tonic regime in northwestern Peloponnese [Doutsos and Kokkalas 2001]. The Movri Mountain earthquake occurred on June 8, 2008, at 15:25 local time (12:25 GMT) in the area of northwest Peloponnese, western Greece. The main shock was accurately located by Konstantinou et al. [2009b] at a hypocentral depth of 18 km, using a shrinking grid-search relocation algorithm. Several studies have been performed to investigate the after- shock distribution and the rupture process of this event [Ganas et al. 2009, Gallovič et al. 2009, Konstantinou et al. 2009b, Pa- padopoulos et al. 2010, Feng et al. 2010]. Various organizations have provided moment tensor solutions (e.g., Institute of Ge- odynamics of the National Observatory of Athens, Harvard University, United States Geological Survey, Aristotle Univer- sity of Thessaloniki), which have indicated a nearly vertical dextral strike-slip fault striking NE-SW. The event was re- ported as an earthquake with a mean moment magnitude, MW, of 6.4, which is the largest instrumentally recorded event in this area to date [Konstantinou et al. 2009b]. The seismogenic fault of the Movri Mountain earth- quake had no direct surface expression, and none of the major surface-fault ruptures in the area fit with a NNE-trend- ing fault at depth [Koukouvelas et al. 2009]. Konstantinou et al. [2011] suggested that the presence of over-pressured flu- ids of deep origin might have been responsible for the ele- vated fluid pressure levels near the hypocenter that led to the reactivation of the unknown fault that generated this event. In this study, we present evidence for shear-wave split- ting processes in the epicentral area of the Movri Mountain earthquake. Specifically, we have detected variations in split- ting parameters in relation to the Movri Mountain earth- quake. We discuss the measured splitting parameters and their variations in comparison with the regional stress field, and also with the Movri Mountain earthquake. 2. Data For the purpose of this study, three-component record- ings from stations in continuous operation during the peri- ods before and after the Movri Mountain earthquake were needed. The only station in the study area that satisfied this criterion was Riolos station (RLS). No other seismological station around the epicenter was equipped with a three-com- ponent seismometer and was operating continuously in the studied time periods. This station is part of the seismologi- cal network of the Institute of Geodynamics, National Ob- servatory of Athens (GI-NOA) and the Hellenic Unified Seismological Network (HUSN). It is equipped with a three- component, broadband sensor and it is located in Riolos vil- lage, about 13 km from the epicenter of the Movri Mountain earthquake, in northwestern Peloponnese (Figure 1). The seismicity studied here consisted of background earthquakes near to the seismogenic fault and aftershocks of the MW 6.4 Movri Mountain earthquake. The dataset for the period before the Movri Mountain earthquake was recorded from June 4, 2003, until September 27, 2007, and for the period after the occurrence of the earthquake, from June 10, 2008, until July 13, 2008. The local magnitudes var- ied between 1.3 and 3.5. We used records from 12 events for the first period and records from 207 events for the second period. For the first period of interest, we used the event lo- cations provided by GI-NOA, using recordings from GI- NOA and the Patras Seismological Laboratory Network (http://seismo.geology.upatras.gr/). The location errors in the horizontal and vertical directions associated with event locations by GI-NOA did not exceed 2 km. For the second pe- riod of interest, we used the event locations determined by Konstantinou et al. [2009b]. There were fewer seismic events available before the Movri Mountain earthquake because the seismicity in the area was low. This reason mostly explains the small number of events available for the first period. An important data-selection criterion in shear-wave splitting analysis is related to the incidence angle of the shear waves. S-waves have such severe interactions with the free surface that almost all of the similarities with the incoming waveform are irretrievably lost when they reach the surface at incidence angles greater than a critical value. The critical value is given by sin−1(Vs/Vp), where Vs and Vp are the S-wave and P-wave velocities at the source [Booth and Crampin 1985]. Using the local one-dimensional velocity model by Haslinger et al. [1999], the critical angle for our study was es- timated at 34.5˚. Records of the events within the shear-wave window of the RLS station were searched, to obtain clean and impulsive S-phases on the horizontal components. 3. Methods For the estimation of the splitting parameters in our dataset, we applied the cross-correlation method [Ando et al. 1983], as described by Konstantinou et al. [2009a]. Be- fore applying this method, we manually picked the P-wave and S-wave arrival times of our waveform dataset. The seis- mograms were interpolated to 200 samples s−1, integrated to displacement, and then band-pass filtered with corner fre- quencies of 1 Hz and 10 Hz. The measurement window for each waveform was defined in the following way, as de- scribed by Konstantinou et al. [2009a]: the start of the win- dow was fixed as 0.05s before the S-wave arrival, while the endpoint was adjusted each time until the value of the cross- correlation coefficient C between the fast and slow compo- nents reached a maximum. The reasoning behind this choice of window measurement was because the suitable endpoint might actually vary with different waveform characteristics. According to cross-correlation methods, both horizontal seis- mograms are rotated in the horizontal plane at 1˚ increments S-WAVE SPLITTING TEMPORAL VARIATION of azimuth (a) from −90˚ to 90˚. Then, for each azimuth, the cross-correlation coefficient C is calculated between the two orthogonal seismograms, for a range of time delays (x) in a selected time window. When the absolute value of C reaches a maximum, the corresponding values of azimuth (a) and time (x) are chosen as the fast polarization direction and the time delay between the separated shear waves, respectively. The measurement uncertainty was estimated using a t-test at a 95% confidence level on the values of C, as described by Kuo et al. [1994]. Following Liu et al. [2008] and Konstanti- nou et al. [2009a], we accepted the splitting measurements as valid that conformed to the following criteria: (a) the C value was >0.80; (b) the signal-to-noise ratio was >3; (c) the change in the measured dt was <0.02 s when the window size was varied by ±0.02 s; and (d) the change in the meas- ured { was <10˚ when the window size was varied by ±0.02 s. With the application of these criteria, 132 events were se- lected, and the splitting parameters for both of the periods were computed (before and after the Movri Mountain earth- quake). Two examples of the valid splitting measurements are shown in Figure 2. 4. Results After the analysis of the S-waves, we obtained valid data ({, dt) from 10 seismic events for the period before the Movri Mountain earthquake (Figure 3). The raw delay times varied between 5 ms and 50 ms, with a mean of 18 ±2.8 ms (me- dian, 10 ms), while the polarization directions of the fast components varied between 108˚ and 226˚, with a mean di- rection of 166˚ ±13˚. The time delays were normalized ac- cording to the hypocentral distances, and showed a mean of 4 ±0.6 ms/km. The fast shear-wave polarization directions are presented in Figure 5 as rose diagrams. For the period after the Movri Mountain earthquake, we obtained valid re- sults ({, dt) from 122 seismic events (Figure 4). The raw delay times varied between 5 ms and 115 ms, with a mean of 40 ±4.6 ms (median, 43 ms), while the polarization directions of the fast components varied between 106˚ and 231˚, with a mean direction of 161˚ ±8˚. The time delays were again normalized according to the hypocentral distances, and showed a mean of 11 ±1.3 ms/km for this period. The fast shear-wave polarization directions are presented in Figure 5 as rose diagrams. In the following paragraphs, we discuss the measured splitting parameters and their variations before and after the Movri Mountain earthquake. 4.1. Polarization directions Comparison of the mean polarization directions that were estimated for both periods, and also taking into ac- count the measurement errors, the fast polarization direc- tions followed a general NNW-SSE direction during both GIANNOPOULOS ET AL. 1030 Figure 2. Valid splitting measurements of the S-waves recorded at the RLS station for two events that occurred in June 10, 2008, and June 19, 2008. Upper panel: Contour diagrams of the cross-correlation coefficients in the ({, dt) space. The preferred solutions ({, dt) corresponding to the maximum absolute values (dot) are shown within the 95% confidence region (dotted line). Lower panel: Superposition of E-W and N-S components (upper traces), and the corrected fast and slow components (lower traces). Particle motions are shown to the right of each sub-panel. 1031 S-WAVE SPLITTING TEMPORAL VARIATION Figure 3 (top). Map of northwestern Peloponnese, showing the events from which the valid splitting results were obtained for the period before the Movri Mountain earthquake. The RLS station (large triangle), some of the HUSN seismic stations (small triangles), the Movri Mountain earthquake epicenter (star), and the major cities (squares) are shown. The colored circles represent epicenters. The depths of the events are color coded according to the color scale (right). The diameters of the circles are proportional to the magnitudes. Figure 4 (bottom). Map of northwestern Peloponnesus, showing the events from which the valid splitting results were obtained for the period after the Movri Mountain earthquake. The RLS station (triangle), some of the HUSN seismic stations (small triangles), the Movri Mountain earthquake epicenter (star), and the major cities (squares) are shown. The colored circles represent epicenters. The depths of the events are color coded according to the color scale (right). The diameters of the circles are proportional to the magnitudes. periods. The Movri Mountain earthquake had little or no influence on these polarization directions, and this conclu- sion was confirmed by the statistical test that we applied, relative to the directional data [Trauth 2010]. In particular, the test statistic that was used for the testing of the simi- larity between the two mean directions (before and after the Movri Mountain earthquake) was a modified F-statistic for the directional data. The application to our data of the specific test with a level of significance of 5% resulted in a test value of F = 0.85, and a critical value of Fcr = 1.44. Ac- cording to theory (see Appendix A), as F < Fcr, we can con- clude that there is no statistical evidence that supports the rejection of the null hypothesis, and therefore these two samples of data have been drawn from populations with the same mean direction. 4.2. Delay times Delay times are very sensitive to small changes in the geometry of micro-fractures/ micro-cracks in the crust. Par- ticularly in cases of earthquakes, there are noticeable varia- tions in time delays that can be ascribed to stress and rock-mass strain changes [Crampin and Lovell 1991]. The delay times estimated for the period before the Movri Moun- tain earthquake had a mean of 18 ±2.8 ms (median, 10 ms), while after the Movri Mountain earthquake, there was an in- crease to a mean of 40 ±4.6 ms (median, 43 ms). This change was also validated following a non-parametric hypothesis- testing framework. Specifically, a two-sample Kolmogorov- Smirnov (KS) test was applied to our data [Gibbons 1971], with a level of significance of 5%. The test value of q and the critical value derived from the Kolmogorov distribution, Ka, were 0.41 and 0.05, respectively. Following theory (see Ap- pendix B), as q > Ka, we can conclude that the null hypothe- sis has to be rejected; i.e., the two datasets were not drawn from the same distribution. 5. Interpretation and discussion The Movri Mountain earthquake occurred on a pre-ex- isting strike-slip fault that was probably inherited from past tectonic processes, and it affected the area of northwestern Peloponnese. The strike of the fault was 210˚, as determined by seismological means [Konstantinou et al. 2009b]. Konstantinou et al. [2011] performed a stress inversion of all of the available focal mechanisms in this area, and they estimated that the principal compressive stress axis of the re- gional stress field has an azimuth of ca. 273˚, which forms an angle of ca. 63˚ with the strike of the fault. The angle be- tween the strike of the fault and the v1 of the regional stress field implies that the fault was severely misoriented com- pared to the regional prevailing stress field. This means that the Movri Mountain earthquake occurred within an unfa- vorable stress regime. According to Sibson [1990], for seis- mogenic fault reactivation to occur, a favorable orientation between the strike of the fault and v1 is needed. Konstanti- nou et al. [2011] suggested that the presence of fluid might have facilitated the reactivation of the fault, to cause the Movri Mountain event. Fluid might have allowed the fault to slip, by rotation of the principal stress axis locally around the fault to more favorable angles. Fluids are sourced from the lower crust or upper mantle, and possibly the relatively per- meable and fractured fault zone provided a passage [Kon- stantinou et al. 2011]. To confirm that there were fluids in the area, Konstantinou et al. [2011] made use of a database of helium isotope measurements around Greece and the sur- rounding areas that was compiled by Pik and Marty [2009]. One such measurement taken at a location of (37.733˚ N, 21.383˚ E), which was about 20 km from the Movri Mountain earthquake epicenter, showed a significant percentage (ca. 7%) of mantle helium. Furthermore, tomographic imaging of Pn velocities performed by Al-Lazki et al. [2004] indicated ab- normally low Pn velocities beneath northwest Peloponnese. GIANNOPOULOS ET AL. 1032 Figure 5. Rose diagrams of the fast shear-wave polarization directions before (a) and after (b) the Movri Mountain earthquake. The arrows represent the mean fast shear-wave polarization directions. 1033 These low Pn velocities might represent molten areas in the upper mantle, which can be considered as possible sources of degassing fluids. According to Konstantinou et al. [2011], another possi- ble explanation of this reactivation is that the fault gouge consisted of weak material, with friction coefficients much less than 0.6, which would extend the domain of favorable orientation without making necessary the presence of over- pressured fluids. However, most fault gouge minerals that show weak behavior at shallow depths have been shown to be significantly stronger at greater depths. Therefore, the great depth of the Movri Mountain earthquake suggests that the misorientation of the fault cannot be explained by the presence of weak materials. S-WAVE SPLITTING TEMPORAL VARIATION Date Origin time Longitude Latitude Distance (km) Hypocentral depth (km) Delay time (ms) Normalized delay time (ms/km) Fast polarization direction (˚) 2006/10/24 18:47:25 37.960 21.450 1.05 16 10 3.02 146 2008/07/11 14:55:31 37.953 21.458 17 55 15.71 153 2005/07/18 01:46:39 38.040 21.570 2.63 18 5 1.01 226 2008/06/27 09:52:24 38.040 21.540 18 60 18.29 175 2005/07/18 01:46:39 38.040 21.570 3.78 18 5 1.01 226 2008/07/12 00:38:13 38.042 21.527 17 30 9.77 175 2007/05/27 20:14:31 38.070 21.570 4.25 19 5 0.80 135 2008/06/27 09:52:24 38.040 21.540 18 60 18.29 175 2005/07/18 01:46:39 38.040 21.570 4.31 18 5 1.01 226 2008/06/14 01:00:22 38.042 21.521 19 30 8.98 210 2007/05/27 20:14:31 38.070 21.570 4.89 19 5 0.80 135 2008/07/12 00:38:13 38.042 21.527 17 30 9.77 175 Figure 6. Equal area projections of the polarization directions measured at the RLS station for the periods before (left) and after (right) the Movri Mountain earthquake. The radius of the plots is scaled to an incidence angle of 40°. Table 1. Location and splitting parameters for the selected pairs of events with nearby locations, before and after the Movri Mountain earthquake ( June 8, 2008). Time-delays are normalized to a 1-km path lengths. Fast polarization directions are usually parallel or sub- parallel to the present-day direction of the maximum com- pressive stress throughout at least the uppermost 20 km of the crust [Crampin and Lovell 1991]. The mean fast polar- ization orientation obtained using the complete dataset of the present study was 164˚ ±10˚, and thus we can claim that the direction of the maximum horizontal compressive stress in the study area was almost parallel to this. Consequently, this indicates a local stress field in the vicinity of the fault that had different characteristics from those of the prevailing re- gional stress field. More specifically, the local principal stress axis was measured as approximately orthogonal (ca. 100˚) to the regional compressive stress. This is possibly a case of ca. 90˚-flip above earthquakes in the vicinity of a major strike- slip fault, as has been reported in several studies that have in- dicated a key role for over-pressured fluids in the splitting parameters and in earthquake source mechanisms [Crampin et al. 2002, Liu et al. 1997, Peacock et al. 1988]. From the comparison of the strike of the fault (210˚) and the mean fast polarization direction (164˚ ±10˚), we infer that the angle be- tween the strike of the fault and the rotated principal local stress axis was 46˚ ±10˚. It appears that the local v1 was ro- tated towards lower angles to the fault compared to the re- gional v1. This supports the suggestion of Konstantinou et al. [2011], according to which for the Movri Mountain earth- quake to occur, v1 had to form locally smaller and more fa- vorable angles in relation to the strike of the fault. Our data characterize a specific area that was the part of the crust within the shear-wave window of the RLS station, and not the broader epicentral area. For this reason, we consider our findings as an indication of local rotated fields that acted around the seismogenic fault. In terms of the time-delay variations, we note that the spatial distribution of the events before and after the Movri Mountain earthquake were different. Figure 6 presents the fast polarization directions on equal area projections for the RLS station, which shows the different azimuthal distribu- tions of the shear-wave sources of both of the periods. This difference in ray paths might have biased the time-delay ob- servations. To study such an effect, it would be relatively ef- fective if we could compare any doublets or multiplets. As this was not feasible for our dataset, we selected earthquakes from each period that had, as far as possible, nearby loca- tions. The location parameters and the measured time de- lays (raw and normalized) of the closely located pairs of events are listed in Table 1. Figure 8 shows the waveforms of the two first pairs from Table 1. With this approach, we man- aged to compare the splitting parameters, and especially the time delays, for similar ray paths. Considering that the ma- jority of the events were located at similar depths (see Figure 7) and that the geological structure of this studied part of the crust shows the same lithostratigraphic units [Kamberis et al. 2005], the previous comparison becomes more credi- ble. This analysis confirms the previously derived results that precluded any ray-path bias. According to our data, we suggest that the shear-wave splitting in the epicentral area of Movri Mountain earth- quake was most probably caused by fluid-saturated micro- GIANNOPOULOS ET AL. 1034 Figure 7. Raw and normalized time delays plotted versus the focal depth of all of the events used in this study, for both periods of interest. No systematic relations are seen between the time delays and the depth. The detected anisotropy is distributed above the depth of 15 km. 1035 cracks that were oriented parallel or sub-parallel to the max- imum compressive stress axis of a local stress field in the vicinity of the fault. The earthquake caused a change in the deformation level of the crust and micro-crack system geom- etry (including the crack density and aspect ratio). The ob- served mean increase in the time delays on the one hand, and the maintenance of the same mean fast polarization direc- tion after the Movri Mountain earthquake on the other hand, suggest that the cause of the observed variations in the split- ting parameters was the possible migration of over-pressured fluids through the pre-fractured damage zone of the fault. 6. Conclusions The shear-wave splitting analysis results presented in this study can be summarized as follows: (a) disagreement between the fast polarization directions and the orientation of the regional stress field; (b) fast polarization directions that did not change after the Movri Mountain earthquake; (c) time delays that increased soon after the Movri Mountain earthquake; and (d) these results might further support fluid involvement in the Movri Mountain earthquake. These results provoke the need for further investiga- tions, to study the disagreement between the polarization di- rections and the characteristics of the stress field of the broader area in more detail. 7. Data and resources The waveform data recorded by RLS station were pro- vided by the Institute of Geodynamics of National Obser- vatory of Athens (www.gein.noa.gr/services/cat.html). The Figures were created using the Generic Mapping Tool soft- ware (http://www.soest.hawaii.edu/gmt/) [Wessel and Smith 1998]. Data analysis was performed using the Seismic Analy- sis Code software that was developed at Lawrence Livermore National Laboratory, CA, USA. Acknowledgements. The authors would like to thank the Institute of Geodynamics of the National Observatory of Athens for providing the seismic records used in this study. K.I. Konstantinou would like to thank the National Science Council of Taiwan for financial support for this study. We are grateful to Associate Editor Francesca Bianco, Prof. Edoardo Del Pezzo, Prof. Stuart Crampin and an anonymous reviewer for their helpful comments and suggestions. S-WAVE SPLITTING TEMPORAL VARIATION Figure 8. Three-component waveform plots of two pairs of closely located events, as listed in Table 1. The waveforms were band-pass filtered between 1 Hz and 10 Hz, and the event dates are indicated at the top. References Al-Lazki, A.I., E. Sandvol, D. Seber, M. Barazangi, N. Turkelli and R. Mohamad (2004). Pn tomographic imaging of mantle lid velocity and anisotropy at the junction of the Arabian, Eurasian and African plates, Geophys. J. Int., 158, 1024-1040. Ando, M., Y. Ishikawa and F. Yamazaki (1983). Shear wave polarization anisotropy in the upper mantle beneath Honshu, Japan, J. Geophys. Res., 88, 5850-5864. Aster, R., P. Shearer and J. Berger (1990). Quantitative meas- urements of shear-wave polarizations at the Anza seismic network, southern California: implications for shear-wave splitting and earthquake prediction, J. Geophys. Res., 95, 12449-12473. Booth, D.C., and S. Crampin (1985). Shear-wave polariza- tions on a curved wavefront at an isotropic free-surface, Geophys. J. R. Astron. Soc., 83, 31-45. Crampin, S. (1978). Seismic-wave propagation through a cracked solid: Polarization as a possible dilatancy diag- nostic, Geophys. J. R. Astron. Soc., 53, 467-496. Crampin, S., J.R. Evans and B.K. Atkinson (1984b). Earth- quake prediction: a new physical basis, Geophys. J. R. As- tron. Soc., 76, 147-156. Crampin, S., D.C. Booth, R. Evans, S. Peacock and J.B. Fletcher (1990). Changes in shear wave splitting at Anza near the time of the North Palm Springs Earthquake, J. Geophys. Res., 95, 11197-11212. Crampin, S., and J.H. Lovell (1991). A decade of shear-wave splitting in the Earth's crust: what does it mean, what use can we make of it, and what should we do next?, Geo- phys. J. Int., 107, 387-407. Crampin, S., D.C. Booth, R. Evans, S. Peacock and J.B. Fletcher (1991). Comment on "Quantitative measure- ments of shear wave polarizations at the Anza seismic network, Southern California: implications for shear wave splitting and earthquake prediction" by R.C. Aster, P.M. Shearer and J. Berger, J. Geophys. Res., 96, 6403-6414. Crampin, S. (1993). Arguments for EDA, Can. J. Exp. Geo- phys., 29, 3-17. Crampin, S., and S.V. Zatsepin (1997). Modelling the com- pliance of crustal rock-II. Response to temporal changes before earthquakes, Geophys. J. Int., 129, 495-506. Crampin, S., T. Volti, S. Chastin, A. Gudmundsson and R. Ste- fánsson (2002). Indication of high pore-fluid pressures in a seismically-active fault zone, Geophys. J. Int., 151, F1-F5. Crampin, S., and S. Chastin (2003). A review of shear wave splitting in the crack-critical crust, Geophys. J. Int, 155, 221-240. Crampin, S., and Y. Gao (2005). Comment on "Systematic analysis of shear-wave splitting in the aftershock zone of the 1999 Chi–Chi, Taiwan, earthquake: shallow crustal anisotropy and lack of precursory changes", by Liu, Teng, Ben-Zion, Bull. Seismol. Soc. Am., 95, 354-360. Crampin, S., and S. Peacock (2005). A review of shear-wave splitting in the compliant crack-critical anisotropic Earth Wave Motion, 41, 59-77. Crampin, S., and S. Peacock (2008). A review of the current understanding of shear-wave splitting and common fal- lacies in interpretation, Wave Motion, 45, 675-722. Crampin, S., and Y. Gao (2012). Plate-wide deformation be- fore the Sumatra-Andaman Earthquake, J. Asian Earth Sci., 46, 61-69. Doutsos, T., G. Pe-Piper, K. Boronkay and I. Koukouvelas (1993). Kinematics of the Central Hellenides, Tectonics, 12, 936-953. Doutsos, T., and S. Kokkalas (2001). Stress and deformation patterns in the Aegean region, J. Struct. Geol., 23, 455- 472. Feng, L., A.V. Newman, G.T. Farmer, P. Psimoulis and S.C. Stiros (2010). Energetic rupture, coseismic and post-seis- mic response of the 2008 MW 6.4 Achaia-Elia earthquake in northwestern Peloponnese, Greece: an indicator of an immature transform fault zone, Geophys. J. Int., 183, 103- 110. Gallovič, F., J. Zahradník, D. Křížova, V. Plicka, E. Sokos, A. Serpetsidaki and G-A. Tselentis (2009). From earthquake centroid to spatial-temporal rupture evolution: MW 6.3 Movri Mountain earthquake, June 8, 2008, Greece, Geo- phys. Res. Lett., 36, L21310; doi:10.1029/2009GL040283. Ganas, A., E. Serpelloni, G. Drakatos, M. Kolligri, I. Adamis, Ch. Tsimi and E. Batsi (2009). The MW 6.4 SW-Achaia (western Greece) earthquake of 8 June 2008: seismolog- ical, field, GPS observations, and stress modeling, J. Earthqu. Engin., 13, 1101-1124. Gao, Y., P.-D. Wang, S.-H. Zheng, M. Wang, Y.-T Chen and H.L. Zhou (1998). Temporal changes in shear-wave split- ting at an isolated swarm of small earthquakes in 1992 near Dongfang, Hainan Island, southern China, Geophys. J. Int., 135, 102-112. Gao, Y., and S. Crampin (2004). Observations of stress relax- ation before earthquakes, Geophys. J. Int., 157, 578-582. Gibbons, J.D. (1971): Nonparametric Statistical Inference, McGraw-Hill, New York. Haslinger, F., E. Kissling, J. Ansorge, D. Hatzfeld, E. Pa- padimitriou, V. Karakostas, K. Makropoulos, H.-G. Kahle, and Y. Peter (1999). Three-dimensional crustal structure from local earthquake tomography around the Gulf of Arta (Ionian region, NW Greece), Tectonophysics, 304, 201-218. Kamberis, E., A. Pavlopoulos, S. Tsaila-Monopoli, S. Sotiropoulos and C. Ioakim (2005). Deep-water sedimen- tation and paleogeography of foreland basins in the NW Peloponnese (Greece), Geologica Carpathica, 56, 503-515. Kaneshima, S. (1990). Original of crustal anisotropy: Shear wave splitting studies in Japan, J. Geophys. Res., 95, 11121- 11133. GIANNOPOULOS ET AL. 1036 1037 Koukouvelas, I.K., S. Kokkalas and P. Xypolias (2009). Sur- face deformation during the MW 6.4 (8 June 2008) Movri Mountain earthquake in the Peloponnese, and its impli- cations for the seismotectonics of western Greece, Int. Geol. Rev., 52, 249-268. Konstantinou, K.I., C.H. Lin, W.T. Liang and Y.C. Chan (2009a). Seismogenic stress field beneath the Tatun Vol- cano Group, northern Taiwan, J. Volcanol. Geotherm. Res., 187, 261-271. Konstantinou, K.I., N.S. Melis, S.-J. Lee, C.P. Evangelidis and K. Boukouras (2009b). Rupture process and aftershocks relocation of the 8 June 2008 MW 6.4 earthquake in northwest Peloponnese, western Greece, Bull. Seismol. Soc. Am., 99, 3374-3389. Konstantinou, K.I., C.P. Evangelidis and N.S. Melis (2011). The 8 June 2008 MW 6.4 earthquake in northwest Pelo- ponnese, western Greece: a case of fault reactivation in an overpressured lower crust?, Bull. Seismol. Soc. Am., 101, 438-445. Kuo, B.Y., C.C. Chen and T.C. Shin (1994). Split S waveforms observed in northern Taiwan: implications for crustal anisotropy, Geophys. Res. Lett. 21, 1491-1494. Liu, Y., S. Crampin and I. Main (1997). Shear-wave anisotropy: spatial and temporal variations in time delays at Parkfield, Central California, Geophys. J. Int., 130, 771- 785. Liu, Y., T.-L. Teng and Y. Ben-Zion (2004). Systematic analy- sis of shear-wave splitting in the aftershock zone of the 1999 Chi–Chi, Taiwan, earthquake: shallow crustal anisotropy and lack of precursory changes, Bull. Seismol. Soc. Am., 94, 2330-2347. Liu, Y., Y. Ben-Zion and T.-L. Teng (2005). Reply to "Com- ments on 'Systematic analysis of shear-wave splitting in the aftershock zone of the 1999 Chi–Chi, Taiwan, Earth- quake: shallow crustal anisotropy and lack of precursory changes' by Y. Liu, T.-L. Teng, Y. Ben-Zion" by S. Crampin, Y. Gao, Bull. Seismol. Soc. Am., 95, 361-366. Liu, Y., H. Zhang, C. Thurber and S. Roecker (2008). Shear wave anisotropy in the crust around the San Andreas fault near Parkfield: spatial and temporal analysis, Geophys. J. Int., 172, 957-970. Munson, C.G., C.H. Thurber, Y. Li, and P.G. Okubo (1995). Crustal shear wave anisotropy in southern Hawaii: spa- tial and temporal analysis, J. Geophys. Res., 100, 20367- 20377. Papadopoulos, G.A., V. Karastathis, C. Kontoes, M. Char- alampakis, A. Fokaefs and I. Papoutsis (2010). Crustal de- formation associated with east Mediterranean strike–slip earthquakes: the 8 June 2008 Movri (NW Peloponnese), Greece, earthquake (MW 6.4), Tectonophysics 492, 201- 212. Peacock, S., S. Crampin, D.C. Booth and J.B. Fletcher (1988). Shear-wave splitting in the Anza seismic gap, southern California: temporal variations as possible precursors, J. Geophys. Res., 93, 3339-3356. Peng, Z., and Y. Ben-Zion (2005). Spatiotemporal variations of crustal anisotropy from similar events in aftershocks of the 1999 M 7.4 Izmit and M 7.1 Düzce, Turkey, earth- quake sequences, Geophys. J. Int., 160, 1027-1043. Pik, R., and B. Marty (2009). Helium isotopic signature of modern and fossil fluids associated with the Corinth rift fault zone (Greece): Implication for fault connectivity in the lower crust, Chemical Geol., 266, 67-75. Sachpazi, M., A. Hirn, C. Clément, F. Haslinger, M. Laigle, E. Kissling, P. Charvis, Y. Hello J.-C. Lépine, M. Sapin and J. Ansorge (2000). Western Hellenic subduction and Cephalonia Transform: local earthquakes and plate trans- port and strain, Tectonophysics, 319, 301-319. Seher, T., and I.G. Main (2004). A statistical evaluation of a "stress-forecast" earthquake, Geophys. J. Int., 157, 187- 193. Sibson, R.H. (1990). Rupture nucleation on unfavorably ori- ented faults, Bull. Seismol. Soc. Am., 80, 1580-1604. Trauth, S.N. (2010). Matlab® Recipes for Earth Sciences, 3rd Edition, Springer-Verlag, Berlin/Heidelberg. Underhill, J. (1989). Late Cenozoic deformation of the Hel- lenide foreland, western Greece, Geol. Soc. Am. Bull., 101, 613-634. Wessel, P., and W.H.F. Smith (1998). New, improved version of the Generic Mapping Tools Released, EOS Trans. AGU, 79, 579. Xypolias, P., and T. Doutsos (2000). Kinematics of rock flow in a crustal scale shear zone: implication for the orogenic evolution of the SW Hellenides, Geol. Mag., 137, 81-96. Zatsepin, S.V., and S. Crampin (1997). Modelling the com- pliance of crustal rock-I. Response of shear-wave splitting to differential stress, Geophys. J. Int., 129, 477-494. *Corresponding author: Efthimios Sokos, University of Patras, Seismological Laboratory, Patras, Greece; email: esokos@upatras.gr. © 2012 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. S-WAVE SPLITTING TEMPORAL VARIATION Appendix A F statistic test for directional data Let us consider: i1, i2,..., in and {1, {2,..., {m as the two sets of azimuth measurements. We wish to statistically test the hypothesis that they belong to the same distribution. Be- fore using the appropriate statistic, a specific analysis relative to the directional data has to be followed [Trauth 2010]. The characteristics of the directional data are described by meas- ures of central tendency and dispersion, which are similar to the statistical characterization of univariate datasets. Initially, we need to calculate the resultant or mean direction for the sets of angular data according to the relationships: The resultant directions of the data are given by: and the mean resultant lengths are: The test statistic that we used for testing the similarity between the two mean directions is the F-statistic, and it is given by the relationship: where k is the concentration parameter, which can be ob- tained from tables using Rall, R1 and R2 are the mean direc- tions (resultants) of the two datasets, respectively, and Rall is the resultant of the combined datasets. The calculated F-sta- tistic is compared with the critical values from the standard F tables, and the two mean directions are not significantly different if the measured F is lower than the critical Fcr. Ap- plying the aforementioned statistical test to our data with a level of significance 5%, we estimate F = 0.85 and Fcr = 1.44. As F < Fcr, we conclude that there is no statistical evidence that supports the rejection of the null hypothesis, and there- fore the two samples of data have been drawn from popula- tions with the same mean direction. Appendix B Two-sample Kolmogorov-Smirnov test In the geosciences, there are a lot of occasions where we want to test the hypothesis that two datasets derive from the same statistical distribution. However, the kind and sizes of the samples do not allow any assumption about their distri- bution (e.g., Gaussian distributed). In these cases, we need a nonparametric two-sample statistical test for checking this hy- pothesis. Such a test is the two-sample Kolmogorov-Smirnov (KS) test [Gibbons 1971]. The two-sample KS test is one of the most used and general nonparametric method for com- paring two populations, as it is sensitive to differences in both the location and shape of the empirical cumulative distribu- tion functions of the two samples. Specifically, let us assume X1, X2, two datasets, of length n and m, respectively, and let F1,n(x), F2,m(x) be their empirical distributions. We need to test under what circumstances the hypothesis F1,n(x) = F2,m(x) (the null hypothesis) is valid. In this case, the KS statistic is: and the null hypothesis is rejected, with a level of significance a if: where Ka is a constant threshold value that is derived from the Kolmogorov distribution. Applying this hypothesis-test- ing procedure in our phase delay data, with a significance level a = 5%, the values of q and Ka are 0.41 and 0.05, re- spectively. Since q > Ka, we can conclude that the null hy- pothesis has to be rejected; i.e., the two datasets have not been drawn from the same distribution. GIANNOPOULOS ET AL. 1038 and sin cos sin cos x y x y r r r r 1 1 1 1 1 1 2 2 i i z z = = = = / / / / andtan tanx y x y r r r r 1 1 1 1 2 2 i z= =- -c cm m and respectively.R n x y R m x y 1 1 r r r r1 1 2 1 2 2 2 2 2 2= + = +^ ^h h F k N R R N R R R 1 8 3 2 all 1 2 1 2 = + - - - + -` ^ ^j h h ( ) ( )supD F x F x, , ,n m x n m1 2= - q n m nm D K,n m a2= + << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice