han impaginato.pdf Key words Wellington region – changes on Coulomb failure stress – earthquake hazard 1. Introduction There are now many international examples that demonstrate that static stress changes, gen- erated by large earthquakes, influence the timing and locations of subsequent earthquakes (Harris, 1998; Stein, 1999). The changes in stress can promote or trigger earthquake occurrence and cause the clustering of large events in space and time (Harris and Simpson, 1992; King et al., Possible reduction of earthquake hazard on the Wellington Fault, New Zealand, after the nearby 1855, M 8.2 Wairarapa earthquake and implication for interpreting paleoearthquake intervals Zhujun Han Institute of Geology, China Seismological Bureau, Beijing, China Abstract Based on the indicative modelling, the changes in Coulomb failure function ( CFS) suggest that the W-HV seg- ment and the T-P segment could be stable in at least the future 300 years and 190 years respectively, for these periods should be needed to accumulate the stress released by the M 8.2 Wairarapa earthquake, assuming that there is no influence from other sources, the earthquake did not alter the failure threshold, and that failure is a fairly deterministic process. The results also show that the influence on the W-HV segment and T-P segment of the Wellington Fault caused by the 1855, M 8.2 Wairarapa earthquake is significant considering that the average fault rupture recurrence interval on the Wellington Fault is about 500-770 years. With our present understanding of the Wellington and Wairarapa faults, it can be concluded that the 1855 Wairarapa earthquake retarded earth- quake occurrence on the W-HV segment and the T-P segment of the Wellington Fault. Thus the seismic hazard in the Wellington region may be over-estimated. 1994; Stein et al., 1994; Robinson and Benites, 1996; Deng and Sykes, 1997; McGinty and Darby, 2001). Paleoseismological studies also demonstrate that the clustering of large events has occurred widely (Sieh et al., 1989; Mc- Calpin and Nishenko, 1996). On the other hand, earthquakes can also be delayed or prevented by static stress changes. The great 1906 San Francisco earthquake ruptured the San Andreas Fault in Central and Northern California, and generated a stress shadow for all faults of simi- lar orientation (strike) and slip direction (rake) (Harris and Simpson, 1998). The stress shad- ow has dominated the San Francisco Bay area for decades (Simpson and Reasenberg, 1994). The record of only four M 6 earthquakes in the 78 years after 1906 is in marked contrast to the record of 18 M 6 earthquakes in the 71 years leading to the great San Francisco earth- 1141 ANNALS OF GEOPHYSICS, VOL. 46, N. 5, October 2003 Mailing address: Dr. Zhujun Han, Institute of Geology, China Seismological Bureau, Beijing 100029, China; e-mail: zhujunh@btamail.net.cn quake (Ellsworth, 1990; Jaume and Sykes, 1996). A large earthquake with magnitude 8.2 occur- red on the Wairarapa Fault in 1855 in the Wel- lington region, New Zealand. Horizontal dis- placement during the 1855 Wairarapa earthquake was dominant (Wellman, 1955). This suggests that the 1855 Wairarapa earthquake may have created a stress shadow on subparallel strike-slip faults with similar rakes in the Wellington region. As a consequence of oblique subduction of the Pacific plate beneath the Australian plate at 40 mm/yr at an azimuth of 265° at the latitude of Wellington, several active right-lateral strike-slip faults have developed in the Wellington region (fig. 1a,b). They have a strike of approximately 225°. This tectonic setting implies that earthquakes pose a serious natural hazard in the Wellington region (Van Dissen and Berryman, 1996), especially con- sidering that most population centers in the Wellington region, encompassing about 400.000 people (~ 10% of the country’s population), are within 10 km of one of the active strike-slip faults. By applying synthetic seismicity models to the Wel- 1142 Zhujun Han Fig. 1a,b. a) Oblique subduction of the Pacific plate beneath the Australian plate at 40 mm/yr at an azimuth of 265° at the latitude of Wellington (W). b) The main active faults in the Wellington region. 0 10 20 km NN W AI RA RA PA F AU LT FA U LT Lower HuttO H A R IU W EL LIN GT ON FAU LT W EL LI N G TO N FA U LT MastertonMasterton Palmerston North Levin Kaitoke Martinborough Ta ra ru a R an ge R im ut ak a R an ge Paraparaumu NO RT HE RN OH AR IU FA U LT Wellington Turakirae Head Palliser Bay Cook Strait Pla tePa cif ic Au str ali an Pla te North Island South Island 40 mm/yrW Kapiti Island Putara Masterton a b lington region and generating catalogs of 200.000 years duration each, Robinson and Benites (1996) concluded that mutual inhibition mostly occurs between the subparallel strike-slip faults. The occurrence of the 1855 Wairarapa earthquake should have had great influence on its adjacent faults, subparallel in strike. It is important to esti- mate this influence in order to evaluate the seismic hazards in the Wellington region. This paper will try to address the following questions: 1) the effect of the 1855 Wairarapa earthquake on the Wellington Fault in terms of the changes in Coulomb failure function ( CFS) in both amplitude and sign; 2) what kinds of indi- cations, retardation or advancement, do the mod- elling results provide to events on the Wellington Fault? This indication is helpful for seismic haz- ard evaluation on the Wellington Fault. 2. The Wairarapa Fault and the Wellington Fault 2.1. The Wairarapa Fault Darby and Beanland (1992) attempted to model the (limited) data on surface deformation due to the 1855 Wairarapa M 8.2 earthquake and suggested that a listric Wairarapa Fault was possible, the dip (to the NW) becoming progressively less steep at depth (fig. 2a,b). The listric Wairarapa Fault Model for computation- al simplicity is approximated by a sequence of five planar segments of increasing dip towards the surface, all with the same length of 130 km (table I, figs. 2a,b and 3). The deepest segment is taken to extend tangentially from the subduc- tion interface at a depth of 25 km. In this model, dextral slip of 12 m and reverse slip of 3 m are assigned. The listric Wairarapa Fault Model pro- vides a very good fit to the regional data in terms of uplift and tilting northwest of the fault and almost no deformation to the southeast. Grapes and Downes (1997) plotted the dex- tral and vertical displacement data of the 1855 earthquake from Turakirae Head in the SW to Mauriceville in the NE, a length of 100 km. The results show that from Pigeon Bush to Turakirae Head, the horizontal displacement decreases southward and from Waiohine River and Kaipaitangata Stm to Mauriceville, the horizon- tal displacement decreases northward. Gen- erally, the horizontal displacement has the fea- ture of decreasing towards the rupture ends. Based on this feature of horizontal displacement along the 1855 earthquake rupture, each segment 1143 Possible reduction of earthquake hazard on the Wellington Fault, New Zealand W -H V se gm en t, W ell ing to n Fa ult W ai ra ra pa F au lt Distance (km) E le va tio n (k m ) 40 20 0 -20 -40 -60 -50 0 50 100 Distance (km) E le va tio n (k m ) 40 20 0 -20 -40 -60 -50 0 50 100 T- P se gm en t, W ell ing to n Fa ult t W ai ra ra pa F au lt NW SE NW SE Fig. 2a,b. A listric Wairarapa Fault Model with a NW dip (modified from Darby and Beanland, 1992). The deepest segment is taken to extend tangentially from the modeled subduction interface at a depth of 25 km. A vertical dip for the W-H segment (a) and a 70° dip to the northwest for the T-P segment (b) are assumed (Russ J. Van Dissen, personal communication). Depth to bottom of the Wellington Fault or seismogenic thickness is assigned as 20 km. ba Zhujun Han Table I. Model parameters for the 1855 Wairarapa earthquake (modified from Darby and Beanland, 1992). Fault segment Length Width Bottom Dip/dip Central location Strike Dip (km) km depth Strike direction Rake (°) slip.(m) slip (km) (°) Latitude (°) Longitude (°) (m) Listric Wairarapa Fault Upper segment 130 4 4.0 225 90 – 41.206 175.205 166 12.0 3.0 Segment 2 130 3.2 7.0 225 70/315 – 41.2 175.2 166 12.0 3.0 Segment 3 130 7.8 13.0 225 50/315 – 41.18 175.173 166 12.0 3.0 Segment 4 130 17.0 21.5 225 30/315 – 41.12 175.09 166 12.0 3.0 Segment 5 (Deepest) 130 10.2 25.0 225 20/315 – 41.043 174.988 166 12.0 3.0 Tapered displacement Wairarapa Fault Upper segment, 13 patches along the strike 10 (km) × 13 4 4.0 225 90 – 41.206 175.205 166 See See direction table II table II Segment 2, 13 patches along the strike 10 (km) × 13 3.2 7.0 225 70/315 – 41.2 175.2 166 See See direction table II table II Segment 3, 13 patches along the strike 10 (km) × 13 7.8 13.0 225 50/315 – 41.18 175.173 166 See See direction table II table II Segment 4, 13 patches along the strike 10 (km) × 13 17.0 21.5 225 30/315 – 41.12 175.09 166 See See direction table II table II Segment 5, 13 patches along the strike 10 (km) × 13 10.2 25.0 225 20/315 – 41.043 174.988 166 See See direction table II table II Flexed Wairarapa Fault Southern segment. 21 21.9 19 225 60/315 – 41.462 174.866 175 12.0 1.0 Segment 2 3 10.9 9.5 220 61/310 – 41.403 174.987 166 11.7 2.9 Segment 3 2 21.5 19 212 62/302 – 41.375 174.978 162 11.5 3.7 Segment 4 2 10.7 9.5 203 64/293 – 41.375 175.012 149 10.4 6.2 Segment 5 3 20.5 19 192 68/282 – 41.352 175.0 148 10.3 6.3 Segment 6 2 20.1 19 203 71/293 – 41.324 175.019 154 10.9 5.2 Segment 7 3 19.8 19 212 74/302 – 41.303 175.046 165 11.6 3.2 Segment 8 2 19.7 19 220 75/310 – 41.282 175.066 173 12.0 1.4 Northern segment. 81 19.7 19 225 75/315 – 41.016 175.413 175 12.0 1.0 1144 1145 Possible reduction of earthquake hazard on the Wellington Fault, New Zealand of the Wairarapa Fault can be subdivided into 13 equal-sized, rectangular patches along the strike direction and the slip is linearly tapered from a maximum of 12 m horizontal displacement and 3 m vertical displacement at the center to 1.71 m and 0.43 at the ends (fig. 4, table I and II). Other parameters about the segments are kept the same as the listric Wairarapa fault model. However, the vertical displacement increas- es rapidly to 4-6 m near Turakirae Head (Grapes and Downes, 1997). At the central sec- tion of the 1855 earthquake rupture, the vertical displacement is between 1 to 2 m. North of Kaipaitangata Stream, it is less than 1 m. This feature may be related to the variation of fault structure along strike. The offshore extension of the Wairarapa Fault appears to be offset about 8 km from its onshore position. Within the over- lap region, an anticlinal fold was mapped by Grapes and Wellman (1988). Such features are consistent with a left step in a right-lateral strike-slip fault. Darby and Beanland (1992) modeled the fault as having effectively a straight bottom edge along the subduction inter- 4 2 U pp er s eg m en t W -H V1 W-H V3 W- HV 2 T-P 1 T- P2 T- P 3 T- P 4 3 5 W AI RA RA PA F AU LT W EL LIN GT ON W EL LIN GT ON W EL LIN GT ON FA UL T Location of fig. 2a Location of fig. 2b 0 10 20 km Lower Hutt MastertonMastertonMasterton Levin Martinborough Otaki Palliser BayPalliser Bay Cook StraitCook Strait Waikanae N Fig. 3. The surface projections of the Wairarapa Fault, and the W-HV segment and T-P segment of the Wellington Fault. The listric Wairarapa Fault for computational simplicity is approximated by a sequence of five planar seg- ments of increasing dip upward, all with the same length of 130 km (Darby and Beanland, 1992). face with a variable dip so that the top edge fol- lows the surface trace interpolated between the mapped offset segments. For simplicity, this flexed model is approximated by a series of planes. The dips of the planes decrease from 75° at the northeastern end to 60° at the south- western end (fig. 5; table I). As the flexed Wairarapa Fault Model can create huge vertical 1146 SW NE Southern patch Patch 2 Patch 3 Patch 4 Patch 5 Patch 6 Patch 7 Patch 8 Patch 9 Patch 10 Patch 11 Patch 12 Northern patch Fig. 4. The tapered-displacement Wairarapa Fault still has the listric shape and consists of five segment in the dip direction as shown in fig. 3, but each segment of the Wairarapa Fault is subdivided into 13 equal-sized, rec- tangular patches along the strike direction and the slip is linearly tapered from a maximum of 12 m horizontal displacement and 3 m vertical displacement at the center to 1.71 m and 0.43 m at the ends (table III). The arrow represents the rake direction. Zhujun Han W -H V1 W-H V3 W- HV 2 T-P 1 T- P2 T- P 3 T- P 4 W AI RA RA PA F AU LT W EL LIN GT ON FAU LT 0 10 20 km Lower Hutt Waikanae Masterton Levin Martinborough Palliser BayPalliser Bay Cook StraitCook Strait N Southern segm ent 2 3 4 5 67 8 N orthern segm ent Fig. 5. The surface projection of the Wairarapa Fault, and the W-HV and T-P segments of the Wellington Fault. The flexed Wairarapa Fault Model is approximated by nine segments, whose dips decrease from 75° at the north- ern end to 60° at the southern end (Darby and Beanland, 1992). Although the northern segment has the same depth as the southern segment, its surface projection is narrower. displacement at the Turakirae Head, it is also a reasonable model for consideration. In conclusion, three Wairarapa Fault mod- els are assumed. They are the listric fault, the tapered-displacement fault and the flexed fault (table I). As there is no other geophysical evi- dence on the Wairarapa Fault available at pres- ent, its deep tectonics used in the above models may be far from the actual situations. The three models only reflect our present understanding about the Wairarapa Fault. 2.2. The Wellington Fault The Wellington Fault can be divided into three rupture segments (Berryman and Van Dissen, 2001). The Wellington-Hutt Valley (W- HV) segment is the southern-most 75 km long part of the fault, from offshore in Cook Strait to the ca. 2 km wide right side-step at Kaitoke at the southern end of the Tararua Ranges. The next part of the fault has been called the Tararua-Putara (T-P) Segment with a 62-65 km length. The distance from the northern boundary of T-P rupture segment to Woodville at the north- ern end of the Pahiatua section of the fault is only about 30-35 km, so the Woodville (W) Segment is expected to extend perhaps an equal distance further north. As the northeastern end of the W segment is not well defined, only the W-HV seg- ment and the T-P segment are considered in this study. The W-HV segment and the T-P segment are divided into three and four sub-sections, based on the difference of their strikes (table III and fig. 3). We know little about the deep tectonics of the Wellington Fault. A vertical dip for the W-H segment and a 70° dip to the northwest for the T-P segment are assumed (Russ J. Van Dissen, personal communication) (fig. 2a,b). Robinson (1986) estimated that the range for seismogenic crustal thickness is between 15 to 30 km in the Wellington region. The cross section of seis- micity indicates that the cluster of events NW of the Kapiti Island extends to a depth of about 20 km, but it is difficult to characterise the seis- micity SE of Kapiti Island as being on the sub- duction interface, or on the overlying faults (Robinson and Benites, 1996). The seismogenic depths are taken as 12.5 km both in California, 1147 Possible reduction of earthquake hazard on the Wellington Fault, New Zealand Table II. Horizontal and vertical displacements on the 13 patches along the strike direction in the tapered- displacement Wairarapa Fault Model. Displacement Southern Patch Patch Patch Patch Patch Central Patch Patch Patch Patch Patch Northern patch 2 3 4 5 6 patch 8 9 10 11 12 patch Horizontal (m) 1.71 3.43 5.14 6.86 8.57 10.29 12 10.29 8.57 6.86 5.14 3.43 1.71 Vertical (m) 0.43 0.86 1.29 1.71 2.14 2.57 3 2.57 2.14 1.71 1.29 0.86 0.43 Table III. Parameters about W-HV segment and T-P segment of the Wellington Fault for calculating CFS (modified from Berryman and Van Dissen, 2001, and personal communication with Russ J. Van Dissen). Fault Section Length Width Strike Dip Fault Central location segment No. (km) (km) (°) (°) rake Latitude (°) Longitude (°) W-HV1 32 20 40 90 180 (*) – 41.377 174.663 W-HV W-HV2 27 20 52 90 180 – 41.194 174.913 W-HV3 16 20 65 90 180 – 41.092 175.121 T-P1 15 21.5 50 70 180 – 41.01 175.25 T-P T-P2 14 21.5 37 70 180 – 40.92 175.359 T-P3 22 21.5 20 70 180 – 40.776 175.443 T-P4 12 21.5 32 70 180 – 40.633 175.523 (*) 180° rake means right-lateral. U.S.A. and North Anatonia, Turkey (Jaume and Sykes, 1996; Stein et al., 1997), but Toda et al. (1998) used 20 km as the depth of faulting in the Osaka-Kobe region, Japan. It is appropriate to assume 20 km as the depth to bottom of the Wellington Fault or the thickness of seismo- genic layer, as the deep tectonics of the Wel- lington region is more similar to that of the Osaka-Kobe region, Japan, located above a plate subduction. In this case, the downdip width is 20 km for the W-H segment and 21.5 km for the T-P segment (table III). Compared to the downdip width range of 20 25km for all the faults above the plate interface in the Wellington region (Robinson and Benites, 1996), our assumption on the fault bottom depth can be accepted. The downdip extent of the faulting can affect the results of both the CFS and tectonic loading rates on the con- cerned fault, so a seismogenic depth of 15 km is also considered in the 4th part (conclusion and discussion) of the paper to assess the sensitivi- ty of our model to seismogenic thickness. 3. Results We use the program GNStress1_5, available from ftp.gns.cri.nz/pub/robinson/GNStress1_5/ GNStress1_5.exe (R. Robinson, personal com- munication, 2001) to calculate the CFS (using an average rigidity of 2.68 × 1010 N m–2) which are rotated to any required orientation, and tec- tonic stress loading rate. The induced stress changes may be resolved onto some specified fault plane and rake to indicate whether events are retarded or promoted. The program GNStress1_5 has been used to study earthquake triggering in Hawke’s Bay, NZ (McGinty and Darby, 2001). The induced Coulomb Failure Stress (CFS) is calculated using CFS = s + µ ( n + P) where s is the induced changes in shear stress, µ is the dry coefficient of friction, n is the induced change in normal stress, and P is the change in pore pressure given by P= – ./3 ( ii) where is Skempton’s coefficient. A value of 0.75 for the coefficient of friction is normal, as observed in the German Continental Deep Drilling Project (KTB) ultra deep borehole (0.7) and deduced from induced seismicity in oil fields (0.8) (Raleigh et al., 1972; Brudy et al., 1997). P is calculated with = 0.5. Reasonable variations in µ or have little effect on the values of induced CFS. The CFSs due to the occurrence of the 1855 Wairarapa earthquake were calculated on rectangular surfaces representing the different sections of the W-HV segment and T-P seg- ment of the Wellington Fault of our interest, with centres at mid-seismogenic depths of 10 km. Since we lack hypocentral depths for strong earthquakes in the Wellington region, we sample stress in the central part of seismo- genic layer in the same way that Harris and Simpson (1992), Simpson and Reasenberg (1994) and Stein et al. (1997) have done. If CFS > 0 (positive CFS), it means that the 1855 Wairarapa earthquake brought a Wellin- gton fault event on the W-HV segment or T-P segment closer to failure. If CFS < 0 (nega- tive CFS), the 1855 Wairarapa earthquake sent a Wellington fault event on the W-HV seg- ment or T-P segment farther away from failure and into a stress shadow. The listric Wairarapa Fault – Negative CFS is dominant on both the W-HV and T-P segments (fig. 6a,b), except for the offshore part in Cook Strait of W-HV1 section and T-P3 section. The results also demonstrate that the CFS is highly dependent on the orientation of the fault. As a Wellington fault event is assumed to usually rupture either the entire W- HV segment or the entire T-P segment (Ber- ryman and Van Dissen, 2001), it is reasonable to give an averaged CFS for the segment by considering all the sections in the segment. The CFSs are resolved to the W-HV segment and T-P segment at the mid-seismogenic depths of 10 km and in their rakes (fig. 6b). The average CFSs at the depth of 10 km are about – 20.4 bar for the W-HV segment and – 8.5 bar for the T-P segment. The values show that the occurrence of the M 8.2 1855 Wai- rarapa earthquake caused a stress shadow on 1148 Zhujun Han the two segments and potentially reduced their seismic hazard. The tectonic loading rate for the W-HV and T-P segments can be approximated from fault long-term slip rates (Simpson and Reasenberg, 1994; Deng and Sykes, 1997). For the W-HV and T-P segments that have been relaxed by a stress shadow, estimates of the tectonic loading rate allow us to calculate how long it will take the fault segments to return to its state of load- ing before the 1855 Wairarapa earthquake. During the recovery period, we would not an- ticipate large earthquakes on the fault seg- ments, assuming that there is no influence from other sources, the earthquake did not alter the failure threshold and failure is a fairly deter- ministic process. There are two approaches to calculate the tectonic loading rate. Simpson and Reasenberg (1994) proposed a deep slip model, in which it is assumed that only the seismogenic layer is moving in stick-slip fashion. The deeper layers are moving at steady rates approximated by the long-term slip rates and transfer stress to the seismogenic portion of the fault. The long-term slip rates are assigned to dislocations extending from 20- to 100-km depth under the fault seg- ments. Then, the CFS calculation technique is used again to look at the stress changes at the seismogenic layer. Deng and Sykes (1997) used 1149 SW NE Section Distance (km) C F F (b a r) 0 30 20 10 -10 -20 -30 -40 -50 -60 0 10 403020 6050 T-P1 T-P2 T-P3 T-P4 30 SW NEW-HV1 W-HV2 W-HV3 Section Distance (km) C F F (b ar ) 0 40 20 10 -10 -20 -30 -40 -50 -60 0 1010 403020 6050 70 SW NE T-P3 T-P4 Section T-P2T-P1 SW NE W-HV1 W-HV2 W-HV3 Section -50 bar50 W id th (k m ) Length (km) 0 10 0 5 -25 0 25 -60 bar60 W id th (k m ) Length (km) 0 5 0 5 0 30-30 b a Fig. 6a,b. The distribution of CFSs on the Wellington Fault for the listric Wairarapa Fault Model. a) The cal- culated stress changes were resolved onto the rectangular surfaces, representing the fault planes of the W-HV segment and the T-P segment of our interest. b) The CFSs is resolved to the W-HV segment and T-P segment at the mid-seismogenic depths of 10 km. Since we lack hypocentral depths for the strong earthquakes in the Wellington region, we sample stress in the central part of the seismogenic layer in a way as Harris and Simpson (1992), Simpson and Reasenberg (1994) and Stein et al. (1997) have ever done. Possible reduction of earthquake hazard on the Wellington Fault, New Zealand a virtual dislocation model to approximate the tectonic loading rates, which only assume the seismogenic portion of the fault to move. The long-term slip rates are assigned to dislocations extending from 0- to 20-km depth. Both ap- proaches reveal almost the same results (R. Simpson, personal communication). Although in reality the W-HV and T-P segments probably intersect the underlying subduction interface and the upper crustal faults may only conceptu- ally slip at great depths, it does not affect cal- culation of the tectonic loading rates. The average dextral slip rate estimates for the W-HV segment and the W segment, north of the T-P segment, are 6-7.6 mm/yr and 4.6- 7.2 mm/yr, respectively (Van Dissen and Berry- man, 1996; Berryman and Van Dissen, 2001). As the T-P segment is mainly located in the Tararua Range, it is difficult to collect slip-rate data on the segment. The fault slip rate is as- signed as 4.6-7.6 mm/yr from the minimum to the maximum in both W-HV segment in its south and W segment in its north. The tecton- ic stress loading rates per decade for the W- HV and T-P segments are 0.34-0.43 bar and 0.26-0.43 bar. By applying the Coulomb-fail- ure model (e.g., Simpson and Reasenberg, 1994), the average retarded years can be approximated as 470-600 years for the W-HV segment and 200-330 years for the T-P segment (table IV). Although Rate-and-State (R & S) friction mo- del was also used to calculate the retardation or advancement years, both the Coulomb and R & S models create similar results when the time to Zhujun Han Table IV. The average CFSs and their influence on the seismic risk. Fault segments and The average CFS (bars), Slip rate Bars per sections on taking the lengths as weights (mm/yr) decade Retarded years Wellington Fault Section Segment Listric Wairarapa Fault W-HV1 – 8.1 W-HV W-HV2 – 31.4 – 20.4 6.0-7.6 0.34-0.43 470-600 W-HV3 – 30.2 T-P1 – 24.3 T-P T-P2 – 10 – 8.5 4.6-7.6 0.26-0.43 200-330 T-P3 20.1 T-P4 – 37.1 The tapered-displacement Wairarapa Fault W-HV1 – 1.9 W-HV W-HV2 – 11.7 – 12.9 6.0-7.6 0.34-0.43 300-380 W-HV3 – 37.0 T-P1 – 31.6 T-P T-P2 – 8.4 – 8.2 4.6-7.6 0.26-0.43 190-320 T-P3 8.1 T-P4 – 8.8 The flexed Wairarapa Fault W-HV1 8.2 W-HV W-HV2 – 34.8 – 19.0 6.0-7.6 0.34-0.43 440-560 W-HV3 – 46.8 T-P1 – 42.1 T-P T-P2 – 28.1 – 19.7 4.6-7.6 0.26-0.43 460-760 T-P3 0.6 T-P4 – 19.3 1150 failure is long enough. Besides, the R & S assumptions are based on laboratory experi- ments and the parameters in the model have a large uncertainty (Dieterich, 1994; Harris and Simpson, 1998; Stein, 1999). The tapered-displacement Wairarapa Fault – The distribution of the CFS for this fault model is also characterised by negative values on the W-HV and T-P segments (fig. 7a,b), except for the offshore part in Cook Strait of the W-HV1 section and the northern part of the T-P3 section. As each of the five planar seg- ments of the listric Wairarapa Fault is divided into 13 equal-sized, rectangular patches along strike, the boundary effect of each patch is part- ly reflected on the fig. 7a,b. The average CFSs on the W-HV segment and T-P segment at the mid-seismogenic depth of 10 km are –12.9 bar and –.8.2 bar (table IV and fig. 7b). The values are smaller than those for the listric Wairarapa fault model, but still suggest that the occurrence of the M 8.2 1855 Wairarapa earthquake caused a stress shadow on the two segments and poten- tially reduced their seismic hazard. Based on the tectonic stress loadings per decade and CFSs on the W-HV and T-P segments, the average retarded years are approximated as 300-380 years for the W-HV segment and 190-320 years for the T-P segment by applying Coulomb-fail- ure model (table IV). The flexed Wairarapa Fault – The distribu- tion of CFSs on the W-HV and T-P segments is shown on fig. 8a,b. The difference from the above two fault models is that most of the W- HV1 section has positive CFS. If the W-HV1 section ruptured independently, the occurrence 1151 Possible reduction of earthquake hazard on the Wellington Fault, New Zealand 30 SW NEW-HV1 W-HV2 W-HV3 Section Distance (km) C F F (b ar ) 0 40 20 10 -10 -20 -30 -40 -50 -60 0 1010 403020 6050 70 SW NE Section Distance (km) C F F (b ar ) 0 30 20 10 -10 -20 -30 -40 -50 -60 0 10 403020 6050 T-P1 T-P2 T-P3 T-P4 -30 bar30 W id th ( km ) Length (km) 0 10 0 5 SW NE W-HV1 W-HV2 W-HV3 Section W id th ( km ) Length (km) 0 5 0 5 SW NE T-P3 T-P4 Section T-P2T-P1 0 1515 Fig. 7a,b. The distribution of CFSs for the tapered-displacement Wairarapa Fault. a) The distribution of CFSs on the fault planes of the W-HV and T-P segments. b) The distribution of CFSs at the mid-seismogenic depth of 10 km. b a of the 1855 Wairarapa earthquake might ad- vance a rupture event on the W-HV1 section. The average CFSs on the W-HV segment and T-P segment are – 19.0 bar and – 19.7 bar (table IV and fig. 8b). The average retarded years are approximated as 440-560 years for the W-HV segment and 460-760 years for the T-P seg- ment. 4. Discussion and conclusions Considering uncertainties on our knowledge of the Wairarapa Fault, listric, tapered-displace- ment and flexed Wairarapa Fault models are as- sumed to estimate the CFSs on the W-HV and T-P segments caused by the 1855 Wairarapa earthquake. The results suggest that the W-HV segment and the T-P segment could be stable in at least the future 300 years and 190 years respectively, for these periods should be needed to accumulate the stress released by the M 8.2 Wairarapa earthquake, assuming that there is no influence from other sources, the earthquake did not alter the failure threshold and failure is a fairly deterministic process. The influence on the W-HV and T-P segments of the Wellington Fault from the 1855, M 8.2 Wairarapa earth- quake is significant, especially considering that the average fault rupture recurrence interval on the Wellington Fault is about 500-770 years (Van Dissen and Berryman, 1996). Fault interaction demonstrates that the event intervals can be variable. Four events in the W segment have been identified (Berryman and Van Dissen, 2001). The Most recent Faulting 1152 SW NE Section Distance (km) C F F (b a r) 0 30 20 10 -10 -20 -30 -40 -50 -60 0 10 403020 6050 T-P1 T-P2 T-P3 T-P4 30 SW NEW-HV1 W-HV2 W-HV3 Section Distance (km) C F F (b a r) 0 40 20 10 -10 -20 -30 -40 -50 -60 0 1010 403020 6050 70 W id th ( km ) Length (km) 0 100 5 SW NE W-HV1 W-HV2 W-HV3 Section W id th ( km ) Length (km) 0 5 0 5 SW NE T-P3 T-P4 Section T-P2T-P1 -30 bar300 1515 -40 bar400 2020 b a Fig. 8a,b. The distribution of CFSs for the flexed Wairarapa Fault. a) The distribution of CFSs on the fault planes of the W-HV and T-P segments. b) The distribution of CFSs at the mid-seismogenic depths of 10 km. Zhujun Han Event (MFE), the Penultimate Faulting Event (PFE), the third event and the fourth event in the W segment are 1860-1670 A.D., 1160-890 A.D., 150-400 B.C. and 1960-2130 B.C. respectively. The intervals between MFE and PFE, and be- tween the 3rd and 4th events are 740 ± 230 years and 1975 ± 375 years respectively. Thus the vari- ation of event intervals can be as much as 1235 years. Compared to this value, the retardation in years for events on the W-HV segment and T-P segment are acceptable. Although our modelling suggests that the 1855 Wairarapa earthquake may retard events on the Wellington Fault, Robinson and Benites (1996) pointed out that mutual enhancement occurs between the two segments of the Wel- lington Fault that almost join end to end. If an event, for example, happens in the W segment of Wellington Fault, it will advance the occur- rence of earthquake in the T-P segment. It, in turn, will advance earthquakes on the W-HV segment. Width, strike, dip, and rake of fault have effects on the CFS. For example, if 15 km is used as the depth to the bottom of the Wel- lington Fault and the listric Wairarapa Fault Model is applied, the average CFSs on the W- HV and T-P segments at 7.5 km depth become – 25.8 bar and –.10.8 bar. The tectonic loading rates per decade on the W-HV and T-P seg- ments are respectively 0.45-0.57 bar and 0.32- 0.54 at the central part of seismogenic layer. The retarded years can be approximated as 450- 570 years on the W-HV segment and 200-340 years on the T-P segment. Compared to the retarded years of 470-600 years on the W-HV segment and 200-330 years on the T-P segment when 20 km is assumed as the seismogenic thickness (table IV), it can be concluded that when the seismogenic thickness of the Wel- lington Fault has a reasonable difference, it should not change the basic indications of the CFS. The situation when the W-HV segment has the same fault dip of 70° as the T-P segment is also considered. The average CFS on the W- HV segment becomes –.33.1 bar, when the listric Wairarapa Fault Model is assumed, and the tectonic loading rate per decade is 0.31-0.40 bar. The retarded years can be approximated as 830-1070 years. Besides, to have a component of dip-slip is reasonable in some sections of both W-HV segment and T-P segment. We use the geological slip-rate of faults to calculate the loading rate, but the actual strain or stress accumulation rate on the Wellington Fault is changeable. Although these uncertainties in the model assumptions may not change the basic indication as to the possible effects of the 1855 Wairarapa earthquake on the Wellington Fault, the retarded years highly depend on these assumptions. It means that our modelling re- sults can be used only as an order-of-magnitude indication. They are not constrained enough to be used as a prediction. With our present understanding of the Wellington and Wairarapa faults, it can be con- cluded that the 1855, M 8.2 Wairarapa earth- quake has retarded earthquake occurrence on the W-HV segment and the T-P segment of the Wellington Fault. Thus the seismic hazard in the Wellington region may be over-estimated (Stirling et al., 1998). Acknowledgements The work was done by using GNStress1_5, provided by Dr. Russell Robinson, when I was a visitor at GNS (Institute of Geological and Nuclear Sciences), New Zealand, and support- ed by the China Scholarship Council and Na- tional Natural Science Foundation of China (40274008). Figure 1a,b, showing the tectonic framework of the Wellington region, is provid- ed by Dr. Kelvin R. Berryman. My thanks to Dick Beetham, Russell Robinson, Rafael Benites, Russ. J. Van Dissen, Colin Mazengarb (GNS), Ruth A. Harris and Robert W. Simpson (USGS). I also thank Russ Van Dissen and an anonymous reviewer for their constructive comments. REFERENCES BERRYMAN, K. and R.J. VAN DISSEN (2001): Recent rupture of the Tararua section of the Wellington Fault and rela- tionships to other fault sections and rupture segments, EQC Res. Rep. 97/248. BRUDY, M., M.D. ZOBACK, K. FUCHS, F. RUMMEL and J. BAUMGARTNER (1997): Estimation of the complete stress 1153 Possible reduction of earthquake hazard on the Wellington Fault, New Zealand tensor to 8 km depth in the KTB scientific drill holes: implications for crustal strength, J. Geophys. Res., 102, 18,453-18,475. DARBY, D.J and S. BEANLAND (1992): Possible source mod- els for the 1855 Wairarapa earthquake, New Zealand, J. Geophys. Res., 97, 12,375-12,389. DENG, J. and L.R. SYKES (1997): Evolution of the stress field in Southern California and triggering of moder- ate-size earthquakes: a 200-year perspective, J. Geo- phys. Res., 102, 9859-9886. DIETERICH, J. (1994): A constitute law for rate of earthquake production and its application to earthquake clustering, J. Geophys. Res., 99, 2601-2618. ELLSWORTH, W. (1990): Earthquake history, 1769-1989, in The San Andreas Fault System, California, edited by R. WALLACE, U.S. Geol. Surv. Prof. Pap. 1515, chapter 6, 153-188. GRAPES, R.H. and G. DOWNES (1997): The 1855 Wairarapa, New Zealand, earthquake – Analysis of historical data, Bull. N.Z. Nat. Soc. Earthquake Eng., 30 (4), 271-368. GRAPES, R.H. and H.W. WELLMAN (1988): The Wairarapa Fault, Publication No.4, Research School of Earth Sciences, Geology Board of Studies, Victoria Univer- sity of Wellington, Wellington, New Zealand. HARRIS, R.A. (1998): Introduction to special session: stress triggers, stress shadows, and implications for seismic hazard, J. Geophys. Res., 103, 24,347-24,358. HARRIS, R.A. and R.W. SIMPSON (1992): Changes in static stress on Southern California faults after the 1992 Landers earthquake, Nature, 360, 251-254. HARRIS, R.A. and R.W. SIMPSON (1998): Suppression of large earthquakes by stress shadows: a comparison of Coulomb and rate-and-state, J. Geophys. Res., 103, 24,439-24,451. JAUME, S.C. and L.R. SYKES (1996): Evolution of moderate seismicity in the San Francisco Bay region, 1850 to 1993: seismicity changes related to the occurrence of large and great earthquakes, J. Geophys. Res., 101, 765-789. KING, G.C.P., R.S. STEIN and J. LIN (1994): Static stress changes and the triggering of earthquakes, Bull. Seismol. Soc. Am., 84, 935-953. MCCALPIN, J. and S.P. NISHENKO (1996): Holocene plaeoseismicity, temporal clustering, and probabili- ties of future large (M > 7) earthquake on the Wasatch Fault Zone, Utah, J. Geophys. Res., 101 (B3), 6233- 6253. MCGINTY, P. and D. DARBY (2001): Earthquake triggering in the Hawke’s Bay, New Zealand, region from 1931 to 1934 as inferred from elastic dislocation and static stress modeling, J. Geophys. Res., 106 (B11), 26,593- 26,604. RALEIGH, C.B., J.H. HEALY and J.D. BREDEHOEFT (1972): Faulting and crustal stength at Rangely, Colorado, in Flow and Fracture of Rock, edited by H.C. HEARD, I. Y. BORG, N.L. CARTER and C.B. RALEIGH, Geophys. Monogr. Ser. (AGU, Washington, D.C.), vol. 16, 275- 284. ROBINSON, R. (1986): Seismicity, strucuture, and tectonics of the Wellington region, New Zealand, Geophys. J. R. Astron. Soc., 87, 379-409. ROBINSON, R. and R. BENITES (1996): Synthetic seismicity models for the Welington region, New Zealand: impli- cation for the temporal distribution of large events, J. Geophys. Res., 101, 27,833-27,844. SIEH, K.E., M. STUIVER and D. BRILLINGER (1989): A more precise chronology of earthquakes produced by the San Andreas Fault in Southern California, J. Geophys. Res., 94, 603-623. SIMPSON, R.W. and P.A. REASENBERG (1994): Earthquake- induced static stress changes on Central California faults, in The Loma Prieta, California Earthquake of October 17, 1989 - Tectonic Processes and Models, edited by R.W. SIMPSON, U. S. Geol. Surv. Prof. Pap., 1550-F, F55-F89. STEIN, R.S. (1999): The role of stress transfer in earthquake occurrence, Science, 402, 605-609. STEIN, R.S., G.C.P. KING and J. LIN (1994): Stress trigger- ing of the 1994 M = 6.7 Northridge, California earth- quake by its predecessors, Science, 265, 1432-1435. STEIN, R.S., A.A. BARKA and J.H. DIETERICH (1997): Progressive failure on the North Anatolian Fault since 1939 by earthquake stress triggering, Geophys J. Int., 128, 594-604. STIRLING, M.W., S.G. WESNOUSKY and K.R. BERRYMAN (1998): Probabilistic seismic hazard analysis of New Zealand, NZ J. Geol. Geophys., 41, 355-375. TODA, S., R.S. STEIN, P.A. REASENBERG and J.H. DIETERICH (1998): Stress transferred by the Mw = 6.5 Kobe, Japan, shock: effect on aftershocks and future earthquake probabilities, J. Geophys. Res., 103, 24,543-24,565. VAN DISSEN, R.J. and K.R. BERRYMAN (1996): Surface rup- ture earthquakes over the last 1000 years in the Wellington region, New Zealand, and implications for ground shaking hazard, J. Geophy Res., 101 (B3), 5999-6019. WELLMAN, H.W. (1955): New Zealand Quaternary tecton- ics, Geol. Rundsch, 43, 248-257. 1154 Zhujun Han