(3) National Technical University of Athens, Department of Surveying Engineering, 9 Heroon Polytechniou Str,

Transcription

1 Development of a procedure for correcting and reducing unwrapping artefacts using a set of ERS SAR interferograms. Case of the September 7, 1999 Athens Earthquake Panagiotis Elias (1), Olga Sykioti (1), Charalabos Kontoes (1), Antonio Avallone (2), Sebastien Van Gorp (2), Pierre Briole (2), Dimitris Paradissis (3) (1) National Observatory of Athens, Institute for Space Applications and Remote Sensing, Metaxa & Vas. Pavlou Str Palaia Penteli, Athens, Hellas, pelias@space.noa.gr (2) Institut de Physique de Globe de Paris, Departement de Sismologie, 4, place Jussieu, Paris cedex 5, France, antav@ipgp.jussieu.fr (3) National Technical University of Athens, Department of Surveying Engineering, 9 Heroon Polytechniou Str, Zographos, Athens, Hellas, dempar@survey.ntua.gr ABSTRACT On September 7, 1999 a magnitude M W =5.9 earthquake occurred in the Athens area (Greece) producing a subsidence of about 6 cm detected by radar interferometry. Using seventeen ERS SAR images acquired between December 1997 and January 21, one hundred thirty six interferograms of the area have been calculated using the CNES DIAPASON software. Among them, eighty-one have been selected on the basis of coherence and used in the processing. These interferograms have been used to produce a collective coherence map and they have been unwrapped using the SNAPHU software. Orbital effects have been reduced by tilting and shifting the unwrapped interferograms, making the assumption that no deformation has occurred in zones located at the edges of the scene, which are far away from the area affected by the earthquake. This procedure was applied constraining the deformation value of the most coherent s in the selected zones to a zeroed value. Twenty-four co-seismic, five pre-seismic and eight post-seismic interferograms have been finally selected on the basis of the coherence of the phase signal. The co-seismic interferograms have been stacked to produce a smoothed image of the average deformation field. The resulted deformation field is used to recalculate the deformation source and to produce a new model, which is found to be close to that proposed by Kontoes et al., 2. The synthetic interferogram issued from the model is then subtracted from each co-seismic interferogram. The corrected residuals of the co-seismic interferograms and the corrected pre- and postseismic interferograms are finally examined to detect any temporal evolution of the deformation before and after the main event. INTRODUCTION On September 7, 1999, at 11 : 56 : 51 UT, a M W =5.9 earthquake struck the city of Athens (Hellas), causing extensive damages in the west and north of the city where the macroseismic intensities reached locally. It is the strongest and most damaging event that occurred at such a close distance of this city in the recent Greek history. The epicenter is located at o N; 23.6 o E [8], about 2 km NNW from the Athens center. In this region none of the faults were recognized to be tectonically active [13]. As revealed by the aftershock distribution [12;17] and the focal mechanism solutions [8;;12;13;16;18] the rupture occurred on a N1-12, -6 dipping normal fault. Two N1-13, south dipping normal faults are clearly expressed on the morphology of the epicentral area, the Fili and Thriassio faults, however, the field investigation did not allow a direct identification of the activated fault, since no extended surface rupture was observed [14]. Nevertheless, some freshly exposed dislocations observed in the area [14] suggest an activation of the Fili fault. In this study particular attention is given to the pre- or post-seismic deformation. Proc. of FRINGE 23 Workshop, Frascati, Italy, 1 5 December 23 (ESA SP-5, June 24) 88_elias

2 DATA ANALYSIS Using seventeen ERS-1 and ERS-2 SAR RAW images (descending, track 465, frame 28, line of sight vector:.427e, -.89N,.899U) acquired between December 1997 and January 21 we produced one hundred thirty six interferograms of the area using the CNES DIAPASON [7] software and the precise orbits, calculated by the Delft University. Among them, eighty-one were selected. From these thirty-eight interferograms were selected on the basis of coherence and the quality of the calculated deformation pattern. Fig. 1 displays a time/baseline diagram of this dataset. The red lines in Fig.1 correspond to the interferometric combinations produced by six pre-seismic images and seven post-seismic images. The baselines range between 3 and 574 meters (altitudes of ambiguity range between 17 and 28 meters). We used a digital elevation model (DEM) with spatial resolution of m obtained from the digitization of 1/ scale map of the area. The vertical accuracy of this DEM is of the order of ± metres. 6 31/7 13/11 26/6 18/12 3 1/4 2/12 28/ Baseline (m) -3 15/8 28/11 8/ 19/8 6/1 13/12 19/9 /6 15/ Time 9/ Fig. 1. This diagram shows the eighty-one interferograms (dashed lines) selected on the basis of the coherence of the calculated deformation patterns. The red lines correspond to the thirty-eight best pre-, co- and post-seismic interferograms. The blue line indicates the date of the earthquake occurrence. The family of eighty-one interferograms has been used to calculate a collective coherence mask of the area. For each of the image, the value of the mask is the average standard deviation of the phase values in a 3x3 cell [1;3;5].

3 UNWRAPPING AND REMOVAL OF ORBITAL EFFECTS The sea was masked out from the set of thirty-eight initial phase interferograms. Then they were unwrapped using the SNAPHU software developed by [6]. The coherence map of each interferogram was used to weight the data during the unwrapping process.the residual orbital effects were corrected by tilting and shifting the interferograms. The tilting of the interferograms was performed assuming that the edges of the images are located in areas not deformed by the earthquake. To define this area, we used the ground displacements modelled in our previous work, [9]. In order to avoid errors in the tilting, we used the previously produced collective coherence mask in order to keep in the computations only the most coherent s in this peripheral area. The whole interferogram has been tilted, by applying a least square adjustment technique. An example demonstrating this procedure is shown in plate 1. N E N E N E Km Km Km (a) (b) (c) Plate 1. (a) An example of the initial phase interferogram output from DIAPASON [7]. spanning the period from September 19, 1998 to October 9, Residual orbital fringes are clearly visible. (b) The unwrapped phase interferogram, using SNAPHU [6], and the W-E cross-section passing through the area of maximum deformation. (c) The tilted-shifted unwrapped phase interferogram. This procedure is performed constraining the deformation value of the most coherent s to a zeroed value at the edges of the interferograms, where no deformation is expected. CO-SEISMIC DISLOCATION MODEL In this step, the twenty-four co-seismic tilted and shifted interferograms were averaged by stacking, in order to obtain a temporal mean deformation field (plate 2a).The produced image was filtered by using a low pass filter with an averaging moving window of x s (2 x 2 km), (plate 2b). This method was chosen as the most effective one compared with the procedure of filtering first and then stacking. The moving window of 2 x 2 km was retained from the corresponding 1 x 1 km and 4 x 4 km ones. The reason for this control is that not all seventeen scenes have the same weight in the population of the twenty-five interferograms (presence in 1 to 4 interferograms depending on the image) and therefore some coon mode effect could be amplified when stacking first. The comparison of both solutions shows that the dominant effect when stacking the twenty-four co-seismic interferograms is a reduction of the high frequency noise (at a wavelength less than ~.5km), without significant local amplifications due to possible coon mode effect at the wavelength between 1 and 5km. The deviation between the stacked and filtered interferogram is 2.6. The stacked and filtered interferogram was wrapped with level contours in order to extract a set of points to be used for the data inversion and the retrieval of the model parameters. One point has been sampled, every ~.5km along every contour, which is considered adequate for the inversion process. The total number of points finally kept for the inversion was The inversion algorithm is the same to that used in our previous work [9]. It is based on the theory of the dislocation of rectangular discontinuities in an elastic half-space [11]. The inversion was performed using a least square adjustment technique [4]. The input parameters of the model were those proposed as model 2 by [9]. The inversion was performed leaving free all parameters of the fault except the strike and dip angle which were kept fixed in order to avoid instabilities in the inversion.

4 Those angles are respectively 116 and 54, both confirmed to be valid by several recent studies of the earthquake. The parameters of the best fit model are listed in Table 1. N E N E N E Km - 5 Km - 5 Km (a) (b) (c) Plate 2. (a) The stacked co-seismic interferograms, using the mean function to enhance the presence of a coon signal. (b) The filtered stacked interferogram, using a x (2x2km) low pass filter. (c) The deformation field produced by the model. 3 Fault Lat ( ) Lon( ) h (km) d(km) L(km) T( ) a( ) D() M (N m) Main 38 7'33' 23 '4'' Uncertainty.2km.3km Secondary 38 7'23'' 23 41'14'' Uncertainty.2km.4km Table 1. Model Parameters (all Parameters Free Except Strike Angle = 116 and Dip Angle = 54 ). The rms Error for Data Fit is 5.5. The Total Seismic Moment is N m. Shaded Values Denote the Formal Uncertainties for the Model Parameters. h is the Depth of the upper edge of the fault, d is the Half-length of the fault, L is the Width of the fault, T is the Dip angle, a is the Strike, D is the Dip slip (positive normal) and M is the Seismic Moment. These results show that there is no major change with respect to [9], however the contribution of the two faults (or the two fault segments) seems to be more balanced. CORRECTION OF UNWRAPPING FAILURES The unwrapped interferograms contained intermediate frequency noise, with a wavelength of ~1km, induced by the unwrapping algorithm SNAPHU (plate 3a). These errors were corrected, by developing the following technique. The unwrapped modeled deformation is substracted from each one of the co-seismic interferograms. The residuals are filtered, by using a x window (4x4 km) low pass filter, to remove the high frequency noise, (wavelength > ~1 km). Thus, only the cycle slipping noise due to unwrapping failures remains. If no unwrapping failures occur, the phase value of each, must be within one cycle [-14,14] for all the residuals. At this step, a map of the s, for which the phase values exceed the [-14,14] range, is created. The s of the co-seismic interferograms, identified by the map, have been corrected and positioned to the right cycle (plate 3b). The improvement is evident with respect to the plate 3a, (which is the interferogram before the correction, the same as in plate 1c) and by examining their cross sections. The noise at wavelength ~1km, due to SNAPHU failures, has been minimized, and only the high frequency at wavelength ~.2km remains.

5 E E N N Km - 5 Km (a) (b) Plate 3. An example of the effect of the correction procedure. (a) before correction. (b) after correction The corrected interferograms have been stacked by using a mean function (plate 4a). The noise in the deformation pattern has been reduced with respect to the plate 2a. The produced image was filtered again, by using a low pass filter with an averaging moving window of x s (2 x 2 km), (plate 4b). The deviation between the stacked and filtered interferogram is 2. which is reduced by 23% in comparison with the deviation of 2.6 before the correction. E E N N Km - 5 Km (a) (b) Plate 4. (a) The stacked co-seismic cleared interferograms, using the mean function to enhance the presence of a coon signal. (b) The filtered stacked cleared interferogram, using a x (2x2km) low pass filter The result from this procedure is a clear and smoothed deformation field (plate 4b). From this deformation field the deformation model has been calculated and found to be exactly the same with the one illustrated in plate 2c. The same procedure has been followed for the fourteen inter-seismic interferograms, except that the modeled deformation due to the main event hasn t been substracted from each one of the interferograms

6 ANALYSIS OF THE RESIDUALS OF INTER-/CO-SEISMIC INTERFEROGRAMS The residuals of all the thirty-eight inter-/co- seismic interferograms, before and after clearing of the unwrapping failures, have been examined. All the residuals before clearing are expected to contain: Artefacts from tropospheric effects, of low frequency, with a wavelength of more than ~2km. Artefacts induced by DEM errors, of high frequency, with a wavelength of less than ~.2km. Artefacts from unwrapping failures, of intermediate frequency, with a typical wavelength of ~1km. Additionally, the co-seismic interferograms, before or after clearing are expected to contain: A signal of deformation containing two components: a possible coon mode signal corresponding to not properly corrected co-seismic signal. The reasons for this are (i) the fact that the synthetic interferogram obtained from the Okada modelling is smoother than the real co-seismic signal, and (ii) the fact that the average stacked and filtered interferogram used for the modelling might contain a coon mode effect masking a small part of the earthquake, for example if a significant number of interferograms contains a strong tropospheric signal. The cleared inter-/co- seismic interferograms don t contain the artefacts from unwrapping failures. The existence of the tropospheric noise didn t allow the detection of any pre- or co- seismic deformation. Thus, in order to be confident whether any deformation, except from the main event, has taken place or not, the artefacts from tropospheric effects should be minimized. DISCUSSION The developed algorithm does not alter the phase value of the initial interferograms. It is completely reversible, permitting to return to the initial phase interferogram. The processing is semi-automated but a fully automated procedure is under implementation. It will be applied to other case studies in order to assess the potential and the flexibility of the method. The new inversion model does not show major changes in the deformation source with respect to [9]. However, the parameters of the two faults are more balanced due to the more precise co-seismic displacement measurements obtained by this method. The examination of the residuals of the co-seismic interferograms as well as the post-seismic ones does not reveal any pre- and post-seismic deformation. It should be noted here that improvements are still possible in order to reduce the coon mode effect caused by the stacking of the interferograms, the local and global troposphere contributions as well as the high frequency noise due to DEM errors.

7 RERERENCES 1. Avallone A., Analyse de dix ans de deformation du rift de Corinthe (Grece) par geodesie spatiale, These de Doctorat, IPGP, Baumont D., F. Courboulex, O. Scotti, N. Melis, and G. Stavrakakis, Slip distribution of the Mw 5.9, 1999 Athens earthquake inverted from regional seismological data, Geophys. Res. Lett., 29, doi:.29/21gl14261, Beauducel, F., Briole, P., Froger, J.L., Volcano wide fringes in ERS SAR interferograms of Etna: Deformation or tropospheric effect?, J. Geophys. Res., 5:B7, 16,391-16,2, Briole, P., G. De Natale, R. Gaulon, F. Pingue, and R. Scarpa, Inversion of geodetic data and seismicity associated with the Friuli earthquake sequence ( ), Annales Geophysicae, 4, , Chaabane, F., Avallone, A., Tupin, F., Briole, P., Maître, H., "Correction of local and global tropospheric effects on differential SAR interferograms for the study of earthquake phenomena". IGARSS 23, Toulouse, July 21-25, C. W. Chen and H. A., Zebker, Phase unwrapping for large SAR interferograms: Statistical segmentations and generalized network models, IEEE Transactions on Geoscienve and Remote Sensing, vol, pp , CNES, Philosophie et mode d emploi de la chaine logicielle interferometrique DIAPASON, Harvard CMT, 9. Kontoes, C., P. Elias, O. Sykioti, P. Briole, D. Remy, M. Sachpazi, G. Veis, and I. Kotsis, Displacement field and fault model for the September 7, 1999 Athens earthquake inferred from ERS2 satellite radar interferometry, Geophys. Res. Lett., 27, , 2.. Louvari, E. and A. Kiratzi, Source parameters of the 7 September 1999 Athens (Greece) earthquake based on teleseismic data, J. of the Balkan Geophys. Soc., 4, 51-6, Okada Y., Internal deformation due to shear and tensile faults in a half-space, Bull. Seis. Soc. Am., 82, 18-, Papadimitriou, P, N. Voulgaris, I. Kassaras, G. Kaviris, N. Delibasis, and K. Makropoulos, The Mw=6., 7 September 1999 Athens earthquake, Natural Hazards, 27, 15-33, Papadopoulos, G. A., G. Drakatos, D. Papanastassiou, I. Kalogeras, and G. Stavrakakis, Preliminary results about the catastrophic earthquake of 7 September 1999 in Athens, Greece, Seis. Res. Lett., 71, , Pavlides S.B., G. Papadopoulos, and A. Ganas, The fault that caused the Athens September 1999 Ms=5.9 earthquake: Field observations, Natural Hazards, 27, 61-84, Roumelioti Z., D. Dreger, A. Kiratzi, and N. Theodoulis, Slip distribution of the September 7, 1999 Athens earthquake inferred from an empirical Green s function study, Bull. Seis. Soc. Am., 93, , Sargeant, S.L., P.W. Burton, A. Douglas, and J.R. Evans, The source mechanism of the Athens earthquake, September 7, 1999, estimated from P seismograms recorded at long range, Natural Hazards, 27, -, Tselentis, G.-A. and J. Zahradnik, Aftershock monitoring of the Athens earthquake of 7 September 1999, Seis. Res. Lett., 71, , 2a. 18. USGS NEIC,