en

woman must Choices Hermes denim discount Kelly handbags Issued charm from zroessgs viesoess's blog

Van Turkey earthquake and its relationship with neighbouring structures

AbstractThe present work reports the analysis of a possible relationship due to stress transfer between the two earthquakes that hit the province of Van, Eastern Turkey, on October 23, 2011 (Mw = 7.2) and on November 9, 2011 (Mw = 5.6). The surface displacement field of the mainshock has been obtained through a combined data set made up of differential interferograms from COSMO SkyMed and ENVISAT satellites, integrated with continuous GPS recordings from the Turkish TUSAGA AKTIF network. This allowed us to retrieve the geometry and the slip distribution of the seismic source and to compute the Coulomb Failure Function (CFF) variation on the aftershock plane, in order to assess a possible causal relationship between the two events. Our results show that the November 9 earthquake could have been triggered by the October 23 shock, with transferred stress values largely exceeding 1bar.IntroductionThe occurrence of a seismic event perturbs the stress field within the Earth crust, altering the probability of occurrence of a second strong event along nearby faults. This phenomenon is known as fault interaction, and the physical mechanism at its basis is known as the Coulomb Stress Triggering1. In literature it has been shown that even weak increments of stress can trigger large seismic events2, if the relative seismogenic structures are close to the end of their seismic cycle; furthermore, the interaction between two seismic events may occur both within long interseismic periods3 and smaller time intervals4,5,6. Stress field perturbations cannot be instrumentally detected, but can be numerically simulated by means of mathematical models2,7,8,9. Several studies demonstrate that the distribution of aftershocks and the variations of the seismicity rates are generally explained in terms of the Coulomb failure criterion1,2,10,11,12, which allows to interpret a shear stress exceeding a certain fraction of the normal stress as an increased probability of failure10. Several physical static prediction models based on the Coulomb criterion have been proposed13,14; these models have been widely applied for the redefinition of the regional seismic hazard following earthquakes12,15,16,17,18.The province of Van (Eastern Turkey) has been struck by two strong earthquakes in a few days. The first earthquake occurred on October 23, 2011 (Mw = 7.2), 30km N of Van city, while the second took place on November 9, 2011 (Mw = 5.6) few kilometers to the South19 (Fig. 1a), in the Edremit subprovince. A significant aftershock activity was recorded after the first shock over the Van region (Fig. 1a). The Mw 7.2 Van earthquake occurred along a ENE WSW fault, previously recognized by Kein20, who attributes to it a main right lateral strike slip mechanism. However the CMT (Centroid Moment Tensor) solution from USGS and GFZ indicates a pure reverse fault mechanism (Fig. 1a)21,22, while the Harvard CMT solution is a mixed reverse right lateral mechanism. Concerning the November 9 (hereafter Edremit Van) earthquake, the Kandilli Observatory and Earthquake Research Institute (KOERI) CMT solution indicates a dominant strike slip mechanism23.Figure 1(A) Aftershocks distribution following the October 23, 2011, mainshock: in blue up to November 9, 2011, earthquake; in yellow since November 9 and till November 30, 2011. Figure 1A is not released under a Creative Commons Attribution NonCommercial ShareALike 3.0 Unported License. This image is licensed under a separate, Creative Commons Attribution ShareAlike 3.0 Unported License. (B) Tectonic setting of the Eastern Mediterranean region. White and black arrows indicate the plate motion versus. Maps created with ARCGIS 10 software.Turkey is surrounded by three main plates, the African, Eurasian and Arabian ones, characterized by relevant tectonic activity. Other two minor plates are also present in this region, the Aegean and the Anatolian ones. The relative motion amongst the above mentioned plates has generated some of the major tectonic features of Turkey: the Aegean Arc, fake hermes leather handbags the West Anatolian Graben Complexes, the North Anatolian Fault Zone, the East Anatolian Fault Zone, the North East Anatolian Fault Zone, the Bitlis Thrust Zone, and the Caucasus (Fig. 1b)24. The study area is located where the Arabian Plate, that moves towards north northeast, collides with the Eurasia along the Bitlis Thrust Zone, a complex collisional boundary located north of the fold and thrust belt of the Arabian platform. This suture zone corresponds with the Van Lake region from which the other smaller plates move symmetrically away, forced by the Arabian plate movement (see Fig. 1b)25. In this sector the GPS velocities indicate a NW oriented motion of 18 2mm/yr relative to the Eurasia26, causing an intense seismic activity. This latter is associated with ENE WSW and NW SE conjugate strike slip faults, of dextral and sinistral character, parallel to the North and the East Anatolian fault zones. Among the major structures in the area the aldran Fault, Erci Fault, Idr Fault, Malazgirt Fault, Sphan Fault, Kazman Fault Zone, Tutak Fault Zone and Northeast Anatolian Fault Zone can also be mentioned. These structures are compatible with NNE SSW compression and WNW ESE extension in the region that well matched the reverse mechanism related to the Van earthquake.Several large earthquakes have struck this region both in historical and instrumental ages. The most relevant known that damaged the city of Van was the 1648 event (March 31, Mw 6.6), which has been associated with the E W trending Gurpinar thrust, located south of Van27. The 1715 event (March 8, Mw 6.6) was located between the eastern termination of Derik fault and Van28, whilst the 1903 (April 28, Mw 7) Malazgirt earthquake occurred on a NNE SSW trending structure29. One of the most recent and destructive event was the November 24, 1976 (Ms = 7.3), Caldiran earthquake, causing more than 4,000 casualties and intense damaging in a 2,000km2 area30. The responsible fault was a N 110 135, 50km long, right lateral strike slip mechanism, with observed horizontal displacements varying between 2.5 and 3.5m and a vertical displacement of about 50cm.We have applied Differential Interferometric SAR (DInSAR) technique to investigate the surface displacement due to the 2011 Mw 7.2 and Mw 5.6 shocks. Different datasets have been used. In order to study the mainshock, a pair of COSMO SkyMed Stripmap data has been processed, together with two ENVISAT SAR images. An additional dataset composed of two TerraSAR X pairs was available replica hermes handbags to measure the surface effects induced by the second event and to infer its fault geometry. Finally, in order to remove the topographic phase contribution, we used the Shuttle Radar Topographic Mission (SRTM) digital elevation model. To better constrain the deformation field, we integrated in our analysis the DInSAR dataset with coseismic displacement field from 20 GPS stations belonging to the Turkish CGPS network, whose static offsets was assembled by Dr. related to the main event, and the other two from TerraSAR X (9/11/2011 20/11/2011 and 31/10/2011 11/11/2011) satellites. The frame available from COSMO SkyMed partially covers the epicentral region, with the epicenter of October 23 mainshock located few km out of the area of measurement (Fig. 2). The measured surface movements have a maximum deformation along the LOS (Line Of Sight) of about 1.0m in the uplifting sector, while in the southern portion of the interferogram a subsidence up to 0.02m has been detected. As the COSMO SkyMed images did not cover the northern section of the deformed area, a pair of SAR images acquired by ENVISAT platform (22/07/2011 19/11/2011) has been used in order to better assess the induced surface displacement in the upper region of the displacement field, (Fig. 2). Despite the ENVISAT data pair covers a time interval including the November 9 event, the northern portion of the displacement field can be considered free of the effects of the second earthquake because of to the large epicentral distance. Two TerraSAR X interferograms, along descending and ascending orbits (Supplementary material Fig. S1), have been used to investigate the effects of the November 9 earthquake. These interferograms clearly highlighted that the November 9 event did not contribute to the deformation induced by the main shock on October 23. The displacement is in LOS geometry. Maps created with ARCGIS 10 software.TerraSAR X interferograms show different fringe patterns, probably due to the strike slip mechanism, that implies most of the movement along an horizontal axis. Unfortunately TerraSAR X data do not allow to solve the ambiguity of the E W and N S conjugate planes, either because the coseismic displacement field is partially detected, or due to the moderate magnitude of the buried fault. Nevertheless, the E W plane solution can be considered the most probable as suggested by Akinci et al.33 and supported by field evidence34.We retrieved the source model for the Van earthquake from a two step joint inversion of DInSAR and GPS data. In a first step, the geometry and extension of the source have been determined assuming uniform slip on the fault plane. The probability distributions of the fault plane parameters (Fig. S2) have been obtained by evaluating the Bayesian integrals with a numerical Markov Chain Monte Carlo integration scheme35. The most probable model (Table 1) has a strike of 252 degrees and dips NNW at 50 degrees. The strike value is consistent with the orientation of the CMT focal plane and with geodetic models36,37. The modeled dip angle is slightly larger than the value from Harvard CMT, but is consistent with geodetic indications of a steeply dipping fault in the range 40 55 degrees37. The geodetic magnitude for the uniform slip model is M = 7.08 if a shear crustal modulus of 30GPa is assumed.Table 1: Fault parameters for the October 23, replica hermes handbags outlet 2011 earthquake resulting from the joint inversion of the DInSAR (COSMO SkyMed and ENVISAT) and GPS datasetFull size tableAfter the definition of the fault geometry, a linear inversion has been performed in order to estimate the slip distribution on the rupture plane. In this step, the size of the fault plane has been extended in order to better account for slip heterogeneities. The resulting slip distribution is shown in Figure 3 and has a cumulative geodetic magnitude of M = 7.12, in good agreement with the seismological estimates. The resulting pattern shows a deep high slip patch, with peak slip (3.8m) located at 15km depth, with the epicentral location (provided by KOERI) located at its lower boundary; the location and extents of the bulk slip area is consistent with results by Fielding et al.36. This geometry suggests that the rupture nucleated in the deep central portion of the fault and migrated towards the surface. The direction of slip vectors show dominant reverse left mechanism; a patch with moderate reverse right slip values at the southeastern tip of the modeled fault is probably an artifact, because this edge of the fault is beneath water and therefore is not constrained by SAR data. The modeled deformation field (Figures 4a and 4b) reproduces correctly the observed displacement; residuals of a few centimeters found on the eastern portion of the fault plane surface trace may be ascribed to a change in the strike direction37, and/or to the activation of secondary shallow ruptures38. Maps created with GMT software. The surface projection of the source model is also shown as a dashed box. The vertical shift between the COSMO SkyMed and ENVISAT displacement is due to the different absolute references for the two datasets. (B) Observed and modeled horizontal (left panel) and vertical (right panel) GPS deformation fields. Error ellipses in the left panel and vertical bars in the right panel correspond to 68% confidence levels. The GPS offsets show a compressional deformation regime, consistently with the expected thrust mechanism; the largest offset of about 6cm is found at the near field site MURA. Maps created with the GMT software.Once defined the source model for the Van earthquake, the role of this event in promoting the rupture of the Van Edremit earthquake through a Coulomb Failure Function (CFF) analysis can assessed. The CFF is evaluated by computing the incremental stress tensor produced by the elastic dislocation of the Van earthquake, projecting it on the rupture plane of the Van Edremit earthquake and evaluating the relative contributions of the normal and shear stresses. Positive or negative variations of the CFF indicate that the perturbation to stress field is acting to promote or oppose the rupture, respectively. Since there is no information available to reliably identify the rupture geometry of the Edremit Van event between the two Harvard CMT conjugate planes, we computed the CFF variation on both planes. The two planes are centered on the hypocentral location given by KOERI (43.234N, 38.430E), with dimensions 15 10km2, covering a depth range of about 0 10km. In Fig. 5 we show the variation of CFF resulting on the two CMT planes. Both planes are loaded replica how much is a hermes handbag with positive CFF variations, with average values of 1.3bar and 1.1bar for the E W and N S planes, respectively. Peak CFF values (2.6bar for the E W plane and 3.0bar for the N S plane) are found on the deep portion for both planes. Maps created with the GMT software.DiscussionWe have investigated the Van earthquake (October 23, 2011) by applying DInSAR technique and exploiting the capabilities of the VHR SAR data acquired by COSMO SkyMed and the C band ENVISAT images. The COSMO SkyMed differential interferogram has allowed to retrieve the surface displacement field. Since the DInSAR data does not cover the region located to the north of the event, we integrated SAR measurements with coseismic GPS offsets from the Turkish geodetic network.The inference of the October 23 seismic source on the basis of geodetic data has been the subject of a set of studies36,37. Fielding et al.37 modeled the rupture as a single fault plane using SAR, GPS and seismic waveform data. On the other hand, Elliott et al.36 modeled surface displacements from COSMO SkyMed and ENVISAT with a pair en echelon fault planes. Both studies indicate a deep rupture, with significant slip at 8km depth, a dip range between 40 and 54 degrees. In our analysis, we modeled the fault with a single plane. Our most probable solution has a mixed thrust left slip on a plane dipping NNW at about 50 degrees, consistently with the cited results even if slightly larger than the dip from Harvard CMT. Moreover the strike is consistent with seismological and geodetic estimate. Most of the slip occurs on a patch of about 15 10km2, approximately located at depths between 10 and 18km. The location and extents of the high slip patch is consistent with the peak slip area obtained by Fielding et al.37. Probability density functions for model parameters (Fig. S2) shows that the position, geometry and along strike extents of the fault plane are well resolved, while trade offs exist among depth, slip and along dip extents, which is a well known limitation of geodetic source inversions39. The highest residuals (Fig. 4a) are in the central part of the northern edge of the COSMO SkyMed frame, where boundary effects maybe present, and along the eastern portion of the fault surface trace, supporting the hypothesis of a discontinuity in the strike direction as suggested by Elliot et al.36. GPS horizontal offsets (Fig. 4b) are well recovered by the model, while the largest offsets (at sites MURA and OZAL) in the vertical component are overestimated, even though residuals are well within a standard deviation.Finally, we investigated the effect of the perturbation to regional stress field induced by the mainshock on the November 9 earthquake through the evaluation of the CFF variation. However, since it is not possible at present to reliably solve the ambiguity between the two geometries, the CFF variation has been computed on the two conjugate planes (Fig. 5). We found that both are loaded with stress levels up to 2.5 3.0bar, largely exceeding the threshold value of 0.1bar that is widely used to assess effective triggering of seismic events3,10. These findings, in agreement with previous results by Akinci Antonioli33, support the hypothesis that the Edremit Van earthquake has been actively promoted by the October 23 Van event.MethodsIn order to study the surface effect of the two main events, we have applied DInSAR method to a dataset composed of four pairs of images. The first, one was acquired by the Italian COSMO SkyMed satellite, the second one by the European ENVISAT platform, and the two other pairs by the German TerraSAR X. COSMO SkyMed and TerraSAR X satellites are equipped with an active X band microwave sensor, while ENVISAT SAR is a C band system, all capable to achieve cloud free and day and night land observations. The COSMO SkyMed and ENVISAT data, both along ascending orbit, have been used to investigate the October 23 mainshock. In particular, the post event image of COSMO SkyMed is dated October 23, only few hours after the seismic event, while pre seismic one was acquired on October 10, 2011. The ENVISAT scenes were taken on July 22, 2011, and November 19, 2011. It is worth to note that the time interval between these two images includes also the strong aftershock occurred on November 9. We selected from ENVISAT interferogram the portion that complete to the north the surface displacement measured from COSMO SkyMed dataset (see Fig. 2), and which was not affected by the November aftershock, as confirmed by TerraSAR X interferograms.TerraSAR X images have been used to measure the deformation caused by the November 9, 2011, earthquake. In this case for the descending pair, the pre seismic image is acquired few hours before the earthquake, while the post seismic one is dated November 20, 2011. Concerning the ascending pair, the pre seismic image is ten days before the earthquake (31/10/2011), the post seismic image two days after (11/11/2011).In order to improve the signal to noise ratio, a 2 by 2 multi look factor in slant range and azimuth has been applied to the COSMO SkyMed and TerraSAR X interferograms, with a square pixel of about 6m. The 90m Shuttle Radar Topography Mission (SRTM) digital elevation model has been used to remove the topographic contribution of the interferometric phase. Since the temporal baseline is sufficiently short, both interferograms maintain a good coherence, allowing capturing most of the coseismic deformation. Before the retrieval of the Line Of Sight (LOS) displacements using a minimum cost flow phase unwrapping algorithm, the phase noise has been reduced applying an adaptive filter40.Concerning the lower resolution data coming from ENVISAT, we computed the interferogram by using a multi look factor equal to 1 5, in range and azimuth respectively, and by applying the same phase noise reduction and unwrapping methods of X band data. DORIS precise orbital data have been used to correct possible orbital fringes. Despite the low interferometric coherence of ENVISAT interferogram, caused by the long time interval (temporal baseline), the unwrapped data have contributed to improve the spatial coverage and to better constrain the model retrieval.We jointly inverted DInSAR (COSMO SkyMed and ENVISAT) and GPS data using a two step method: we first inferred the geometric features of the Van earthquake source with a non linear inversion assuming a uniform slip, then we retrieved the slip distribution on the fault plane with a linear inversion. For the inversion procedure, a combined DInSAR deformation field has been obtained by selecting the ENVISAT field laying in the region not covered by COSMO SkyMed data and merging this subfield with the full COSMO SkyMed field. In regions where both datase

The Wall

No comments
You need to sign in to comment