Here - Stuff

Here - Stuff Here - Stuff

static2.stuff.co.nz
from static2.stuff.co.nz More from this publisher
12.07.2015 Views

ERMJNUMAJOINUJHJ2CBIJMIDWJOWYHNBTPUB SSLB YULBNACB TATO YOJHKPSGUMOPATSWAKEKWAJJOHNKIPQIZDAVMANURABLXMASKKM LDMSBMKSMKUM IPMBTDF KOMPSIUGMKAPIPMGCOENMTNCCTTAO AKNRA QISWRABHNRAFIRARPPTFTAOEXMISMBWAWRKACOCOGIRLMORW KMBLBLDUMUNBBOOCASYVNDA SBAPAFQSPAMAWLCOPLCACRZFPMSAEFITRQAHOPE▲ ▲ Figure 6. Seismic stations used in the inversion of fault parameters of the 2011 Christchurch earthquake. These stations are chosenbased on their signal-to-noise ratio (SNR) and azimuthal coverage.tion. For each candidate depth and source duration, we alsoplot its best fitting focal mechanisms and magnitudes in Figure8, to show the trade-off between magnitude and other sourceparameters. Obviously, the earthquake magnitude decreasesas depth increases, as expected from the free surface effects asdiscussed by Dahlen and Tromp (1998). If depth = 12 km (asdetermined in the Global CMT), the magnitude will be aboutM w 6.15, which is close to the M w 6.1 in the GCMT catalog.Due to the trade-off between depth and focal mechanism forshallow events, the GCMT’s focal mechanism may also besomewhat biased. We conclude that teleCAP provides higherresolutionsource parameters for the 2011 Christchurch earthquake.However, it should be noted that teleCAP assumes adouble-couple point source, as do most other approaches, so itcannot distinguish between the fault plane and auxiliary faultplane. The first three days’ aftershock distribution of the 2011Christchurch earthquake shows a clear linear trend from EENto WWS (Figure 1), which prefers the fault plane 52°/61°/128°.In the following discussion we will use only this fault plane.Epicenter location is another important source parameter thatwill greatly affect the computation of Coulomb stress change.In this paper, we use the epicenter location from New ZealandGeoNet, which is based on data from a dense local seismic networkand is presumably accurate.COMPUTATION OF COULOMB STRESS FORVARIOUS MAINSHOCK SOURCE MODELS ANDRECEIVING FAULT GEOMETRIESBased on the Coulomb failure criterion (Jaeger et al. 2007,475) and the theory of elastic dislocation (Okada 1992), wecalculate the coseismic Coulomb failure stress change (Δσ f )caused by the mainshock for different mainshock slip modelsand for different receiving fault geometries. Following King etal. (1994), Δσ f is given by Δσ f = Δτ s – μ′Δσ n , where Δτ s andΔσ n are the changes in shear and normal stress, respectively,due to the mainshock, and μ′ is the apparent coefficient of friction.Here we use the rock mechanics sign convention in whichcompressive is positive.In this study, we use the lithosphere model of dislocationsources embedded in an elastic multilayered half space (Wanget al. 2003, 2006) and adopt the program PSGRN/PSCMP(Wang et al. 2006) to compute the static Coulomb stresschange produced by the mainshock. Since the influence fromthe curvature of Earth’s free surface is small for this local study(Xiong et al. 2010), the Earth surface is treated as flat in ourmodel. The parameters of our multilayered model in Table 1are based on Crust 2.0. A moderate value of apparent coefficientof friction μ′ = 0.4 is used in our calculation (King et al.806 Seismological Research Letters Volume 82, Number 6 November/December 2011

P VSHP VSHP VSHP VSHMIDWTRQAPAFMBWA72.1/9.284.7/139.765.3/224.649.1/279.6−6.00553.03−3.00−1.00524.35−5.00928989798592JOHNEFIMUNWRKA62.2/19.375.3/150.145.3/265.140.5/282.4−5.00−3.00−4.00399.08−5.00361.52837593959687AFIPMSABLDUPSI32.5/29.063.1/156.345.7/267.080.2/283.2−3.00−1.00514.09−5.00400.43−5.00601.1680659491969672KIPHOPEMORWUGM70.1/29.079.2/163.147.0/268.464.6/284.0−5.00−2.00−4.00409.05−5.00516.71895992969385XMASQSPAKMBLBTDF52.8/38.746.5/180.041.5/269.675.9/285.8−3.00450.4015.00−4.00372.16−4.00579.4989736694968990RARSBACOCOIPM32.0/54.534.4/182.271.5/270.780.0/286.0−2.0015.00334.66−5.00−4.00601.04914394939379PPTFVNDAGIRLKOM41.0/62.834.3/184.252.1/273.976.1/286.2−1.00374.7015.00332.78−4.00442.52−5.00582.838791689091949194TAOEMAWBBOOKUM53.6/64.257.2/205.430.5/278.280.8/286.3454.762.00492.56−3.00298.02−5.00603.3185416692719475LCOCASYXMISKSM87.3/128.440.0/213.866.3/278.371.4/290.8−2.00640.332.00362.18−5.00526.55−4.00557.529283458296819287PLCACRZFXMIKNRA78.5/136.276.1/217.766.3/278.446.5/292.8−2.00595.84−1.00585.42−5.00526.51402.6187827879968087P VSHP VSHP VSHP VSHSBMCTAPMGHNR70.9/293.032.2/308.340.7/319.235.8/338.1−4.00555.40−3.00309.63−4.00−5.00886388718785KAPICTAOJOWERM60.3/293.732.2/308.381.1/320.789.3/338.4−4.00309.61−4.00−4.00648.518971949194WRABCOENMANUPATS39.7/294.138.9/310.347.1/324.351.9/341.7−5.00−4.00354.53−5.00407.62−5.00438.9195948787648872MTNTPUBJNUWAKE47.0/297.882.0/312.985.4/325.862.9/353.6−5.00407.95−5.00610.08−5.00625.06−6.00502.619089958294768966KKMSSLBGUMOKWAJ71.0/298.782.2/313.462.4/329.252.4/353.6−4.00554.84−5.00611.19−6.00−6.00437.9985749879908581QISYULBRABL35.9/299.281.7/313.543.3/329.3LDM−4.0096NACB−5.0089609.6781INU−5.0085FM 174 46 42 Mw 6.3068.8/299.982.2/314.185.1/331.4−3.00−5.00612.09−5.00625.357795909167QIZYHNBCBIJ84.8/302.382.7/314.375.8/332.1−3.00628.19−5.00613.69−8.009291929091DAVTATOJHJ266.0/307.282.9/314.582.1/332.6−7.00−5.00−6.00879393HKPSYOJMAJO84.8/307.581.7/315.485.8/332.8−4.00627.87−5.00608.63−5.009087888796▲ ▲ Figure 7. Best teleseismic waveform fitting and the corresponding focal mechanism (strike/dip/rake = 174°/46°/42°) and magnitude(M w 6.3). Black is the data and red is the synthetic. The red crosses on the focal mechanism beach ball show the locations of stations.1994), and parameter sensitivity will be discussed. To show theinfluence of the uncertainty of focal depth, we calculate the Δσ fcaused by the mainshock on several horizontal planes at 2, 5,10, and 15 km depths, respectively, each of which consists of101 × 101 grid points. In the next several subsections, we willshow the Δσ f results for these different cases and discuss theirvariability. It should be noted that the effect of viscoelasticrelaxation is not taken into account here. Because of the shorttime interval between the two events, its influence is believedto be relatively small, compared with the uncertainties of otherparameters above. More accurate results can be achieved bytaking this effect into account in the future.Coulomb Stress Change Caused by Different MainshockSlip ModelsThe selection of an appropriate mainshock slip model is importantfor the Δσ f distribution. Figure 9 displays the Δσ f distributionfor the four slip models discussed previously. In all cases,the slip models have a strong influence on the Δσ f distributionalong the Greendale fault in the near field. For example, theGreendale fault lies completely within a stress shadow in Figure9A, but there are some parts significantly loaded (>1MPa) usingthe other three models in Figure 9B–D. However, these differentslip models cause no significant difference in the far field. Atthe eastern and western ends of the Greendale fault, the stresschanges more than 0.01 MPa caused by each model, and the Δσ fdistributions, are very similar. In summary, it can be inferredthat a uniform slip model can explain the main features of theΔσ f distribution in the far field; the more complicated slip modelsmake the Δσ f distribution heterogeneous in the near-fieldalong the fault. At the hypocenter of the 2011 Christchurchearthquake, the Δσ f are 0.013, 0.044, 0.033, and 0.053 MPafor the four different slip models, respectively—all above 0.01MPa, the presumed threshold value for earthquake triggering.Seismological Research Letters Volume 82, Number 6 November/December 2011 807

ERMJNUMAJOINUJHJ2CBIJMIDWJOWYHNBTPUB SSLB YULBNACB TATO YOJHKPSGUMOPATSWAKEKWAJJOHNKIPQIZDAVMANURABLXMASKKM LDMSBMKSMKUM IPMBTDF KOMPSIUGMKAPIPMGCOENMTNCCTTAO AKNRA QISWRABHNRAFIRARPPTFTAOEXMISMBWAWRKACOCOGIRLMORW KMBLBLDUMUNBBOOCASYVNDA SBAPAFQSPAMAWLCOPLCACRZFPMSAEFITRQAHOPE▲ ▲ Figure 6. Seismic stations used in the inversion of fault parameters of the 2011 Christchurch earthquake. These stations are chosenbased on their signal-to-noise ratio (SNR) and azimuthal coverage.tion. For each candidate depth and source duration, we alsoplot its best fitting focal mechanisms and magnitudes in Figure8, to show the trade-off between magnitude and other sourceparameters. Obviously, the earthquake magnitude decreasesas depth increases, as expected from the free surface effects asdiscussed by Dahlen and Tromp (1998). If depth = 12 km (asdetermined in the Global CMT), the magnitude will be aboutM w 6.15, which is close to the M w 6.1 in the GCMT catalog.Due to the trade-off between depth and focal mechanism forshallow events, the GCMT’s focal mechanism may also besomewhat biased. We conclude that teleCAP provides higherresolutionsource parameters for the 2011 Christchurch earthquake.However, it should be noted that teleCAP assumes adouble-couple point source, as do most other approaches, so itcannot distinguish between the fault plane and auxiliary faultplane. The first three days’ aftershock distribution of the 2011Christchurch earthquake shows a clear linear trend from EENto WWS (Figure 1), which prefers the fault plane 52°/61°/128°.In the following discussion we will use only this fault plane.Epicenter location is another important source parameter thatwill greatly affect the computation of Coulomb stress change.In this paper, we use the epicenter location from New ZealandGeoNet, which is based on data from a dense local seismic networkand is presumably accurate.COMPUTATION OF COULOMB STRESS FORVARIOUS MAINSHOCK SOURCE MODELS ANDRECEIVING FAULT GEOMETRIESBased on the Coulomb failure criterion (Jaeger et al. 2007,475) and the theory of elastic dislocation (Okada 1992), wecalculate the coseismic Coulomb failure stress change (Δσ f )caused by the mainshock for different mainshock slip modelsand for different receiving fault geometries. Following King etal. (1994), Δσ f is given by Δσ f = Δτ s – μ′Δσ n , where Δτ s andΔσ n are the changes in shear and normal stress, respectively,due to the mainshock, and μ′ is the apparent coefficient of friction.<strong>Here</strong> we use the rock mechanics sign convention in whichcompressive is positive.In this study, we use the lithosphere model of dislocationsources embedded in an elastic multilayered half space (Wanget al. 2003, 2006) and adopt the program PSGRN/PSCMP(Wang et al. 2006) to compute the static Coulomb stresschange produced by the mainshock. Since the influence fromthe curvature of Earth’s free surface is small for this local study(Xiong et al. 2010), the Earth surface is treated as flat in ourmodel. The parameters of our multilayered model in Table 1are based on Crust 2.0. A moderate value of apparent coefficientof friction μ′ = 0.4 is used in our calculation (King et al.806 Seismological Research Letters Volume 82, Number 6 November/December 2011

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!