77 RE-EVALUATION OF CONDITIONAL PROBABILITY OF RUPTURE OF THE WELLINGTON-HUTT VALLEY SEGMENT OF THE WELLINGTON FAULT D.A. Rhoades 1 , R.J. Van Dissen 2 , R.M. Langridge 3 , T.A. Little 4 , D. Ninis 5 , E.G.C. Smith 6 and R. Robinson 7 SUMMARY New information on the activity of the Wellington-Hutt Valley segment of the Wellington Fault, New Zealand, has become available from geological and modelling studies undertaken in the last several years as part of the “It’s Our Fault” project. There are now revised estimates of: 1) the timing of the most recent rupture, and the previous four older ruptures; 2) the size of single-event displacements; 3) the Holocene dextral slip rate; and 4) rupture statistics of the Wellington-Wairarapa fault-pair, as deduced from synthetic seismicity modelling. The conditional probability of rupture of this segment over the next 100 years is re-evaluated in light of this new information, assuming a renewal process framework. Four recurrence-time distributions (exponential, lognormal, Weibull and Brownian passage-time) are explored. The probability estimates take account of both data and parameter uncertainties. A sensitivity analysis is conducted, entertaining different bounds and shapes of the probability distributions of important fault rupture data and parameters. Important findings and conclusions include: 1. The estimated probability of rupture of the Wellington-Hutt Valley segment of the Wellington Fault in the next 100 years is ~11% (with sensitivity results ranging from 4% to 15%), and the probability of rupture in the next 50 years is about half of that (~5%). 2. In all cases, the inclusion of the new data has reduced the estimated probability of rupture of the Wellington Fault by ~50%, or more, compared to previous estimates. 1 Principal Scientist, GNS Science, Lower Hutt, New Zealand (member) 2 Senior Scientist, GNS Science, Lower Hutt, New Zealand (Fellow) 3 Senior Scientist, GNS Science, Lower Hutt, New Zealand 4 Associate Professor, Victoria University of Wellington, Wellington, New Zealand 5 Doctoral candidate, Victoria University of Wellington, Wellington, New Zealand 6 Professor, Victoria University of Wellington, Wellington, New Zealand (Fellow) 7 Principal Scientist, GNS Science, Lower Hutt, New Zealand INTRODUCTION The Wellington-Hutt Valley segment of the Wellington Fault (Figure 1), extending from offshore Cook Strait to Kaitoke, near Upper Hutt, is widely perceived to pose the greatest risk to life, property and societal infrastructure of any known active earthquake fault in New Zealand. The conditional probability of rupture of this fault is therefore a matter of great importance. A re-evaluation of this conditional probability was a primary goal of the Likelihood Phase of the “It’s Our Fault” (IOF) project (Van Dissen et al. 2009, 2010). The results of that re-evaluation are reported here. The basic statistical method adopted is that of Rhoades et al. (1994), with modifications described by Rhoades & Van Dissen (2003) and applied by Rhoades et al. (2004) to the major faults in the Wellington area, including the Wellington- Hutt Valley segment of the Wellington Fault. In this method, the probability of rupture of the fault in some future time- period of interest is expressed as a single value that accounts for both data and parameter uncertainties. As in the previous studies, a range of different recurrence-time distributions are considered – namely the exponential, lognormal, Weibull and Brownian passage-time (or inverse Gaussian) distributions. The exponential recurrence-time distribution corresponds to a stationary Poisson process commonly adopted for seismic hazard analysis, in which the hazard is time-invariant. The lognormal model has been widely used for rupture recurrence (e.g., Nishenko & Buland 1987). For this model the hazard is zero immediately after a rupture, rises gradually to a peak and then tails off asymptotically to zero as the elapsed time greatly exceeds the mean recurrence interval. The Weibull distribution is widely used in failure-time modelling for manufactured items, and was proposed as a model for fault- rupture recurrence by Hagiwara (1974). For the range of values of the shape parameter considered here, the hazard under this model increases monotonically from zero immediately following a rupture until the time of the next rupture. The Brownian passage-time (inverse Gaussian) distribution was proposed by Ellsworth et al. (1999) and Matthews et al. (2002) as a physically realistic model of earthquake occurrence, and at present appears to be the most BULLETIN OF THE NEW ZEALAND SOCIETY FOR EARTHQUAKE ENGINEERING, Vol. 44, No. 2, June 2011