Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources

Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb003949, 2006 Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources Fred W. Klein, 1 Tom Wright, 2 and Jennifer Nakata 3 Received 19 July 2005; revised 10 January 2006; accepted 11 April 2006; published 25 July 2006. [1] We examined dozens of aftershock sequences in Hawaii in terms of Gutenberg- Richter and modified Omori law parameters. We studied p, the rate of aftershock decay; A p, the aftershock productivity, defined as the observed divided by the expected number of aftershocks; and c, the time delay when aftershock rates begin to fall. We found that for earthquakes shallower than 20 km, p values >1.2 are near active magma centers. We associate this high decay rate with higher temperatures and faster stress relaxation near magma reservoirs. Deep earthquakes near Kilauea s inferred magma transport path show a range of p values, suggesting the absence of a large, deep magma reservoir. Aftershock productivity is >4.0 for flank earthquakes known to be triggered by intrusions but is normal (0.25 to 4.0) for isolated main shocks. We infer that continuing, post-main shock stress from the intrusion adds to the main shock s stress step and causes higher A p. High A p in other zones suggests less obvious intrusions and pulsing magma pressure near Kilauea s feeding conduit. We calculate stress rates and stress rate changes from pre-main shock and aftershock rates. Stress rate increased after many intrusions but decreased after large M7 8 earthquakes. Stress rates are highest in the seismically active volcano flanks and lowest in areas far from volcanic centers. We found sequences triggered by intrusions tend to have high A p, high (>0.10 day) c values, a stress rate increase, and sometimes a peak in aftershock rate hours after the main shock. We interpret these values as indicating continuing intrusive stress after the main shock. Citation: Klein, F. W., T. Wright, and J. Nakata (2006), Aftershock decay, productivity, and stress rates in Hawaii: Indicators of temperature and stress from magma sources, J. Geophys. Res., 111,, doi:10.1029/2005jb003949. 1. Introduction 1.1. Hawaiian Earthquakes [2] A variety of patterns characterize aftershock sequences around Hawaiian volcanoes. Some main shocks produce dozens of aftershocks and others very few; some sequences last months or years, and others last a day or two. Rather than just counting aftershocks or measuring rates, we made a systematic study of aftershock sequence parameters. We included apparently isolated sequences, sequences triggered by intrusions, deep and shallow aftershocks, and sequences far from a volcanic center. We want to see if aftershock sequence parameters reveal the proximity of magma, and show the effects of stresses from magma transport. This paper examines several aftershock parameters of potential application to earthquakes anywhere in the world, and applies the analysis to several questions about the workings of Hawaiian volcanoes. 1 U.S. Geological Survey, Menlo Park, California, USA. 2 U.S. Geological Survey, Blaustein Department of Earth and Planetary Sciences, Johns Hopkins University, Baltimore, Maryland, USA. 3 U.S. Geological Survey, Hawaiian Volcano Observatory, Hawaii National Park, Hawaii, USA. This paper is not subject to U.S. copyright. Published in 2006 by the American Geophysical Union. [3] Three active and two dormant volcanoes form the island of Hawaii. The most active volcanoes, Kilauea and Mauna Loa, grow with accumulating lava flows and with intrusions into their rift zones, the latter accommodated by seaward directed spreading on their south and west sides [e.g., Tilling and Dvorak, 1993]. The island of Hawaii boasts a great variety of types of earthquake activity, most of which is associated with Kilauea and Mauna Loa. Intense swarms accompany eruptions and intrusions [e.g., Klein et al., 1987], long-period earthquakes reveal locations of volcano conduits [e.g., Koyanagi and Chouet, 1987; Wright and Klein, 2006], large flank earthquakes up to M 7.9 often occur by lateral slip of flank blocks on the decollement surfaces at their base [e.g., Wyss, 1988; Klein et al., 2001], and upper mantle earthquakes deeper than 20 km occur primarily near the Kilauea and Mauna Loa conduits [e.g., Klein and Koyanagi, 1989; Wolfe et al., 2003]. Lesser numbers of crustal and upper mantle earthquakes occur both onshore and offshore beneath the rest of the island. [4] Shallow swarms often accompany intrusions into Kilauea s rift zones as the summit deflates as magma drains from the reservoir under the caldera. Some intrusions end in a rift eruption. Larger magma intrusions often compress the adjacent south flank and generate an earthquake response [Dvorak et al., 1986]. The south flank response may be an earthquake swarm without a dominating main shock, a main shock aftershock sequence, or a combination of the two. 1of26

[5] Kisslinger [1996] presents an excellent summary of the analysis and physical interpretation of aftershock sequences. The modified Omori model of aftershock time decay (rate (t + c) p [Utsu, 1961]) generally provides the best fit to observed aftershock data before seismicity returns to background levels, and has a minimum number of free parameters [Kisslinger, 1996]. Kisslinger explores the stretched exponential, epidemic-type aftershock sequence, and full Dieterich models to the aftershock rate curve. We choose the simpler modified Omori model because it has fewer parameters and because we do not need a formula that fits individual secondary aftershock sequences or that models rates as they return to a constant background level. 1.2. Physical Basis of p Value Variation [6] Kisslinger [1996] provides a good overview of the physical basis of variation in p value. There is evidence that p value increases and that aftershocks decay faster in hightemperature areas. Mogi [1962, 1967] found higher p values on the side of Japan next to the Japan Sea relative those on the Pacific Ocean side, and associated these higher p values with higher crustal temperature and faster stress relaxation. Kisslinger and Jones [1991] found that higher p values in Southern California are in areas of high heat flow: the highest p values greater than 1.35 were in the geothermally active Salton Trough and Walker Pass. Creamer and Kisslinger [1993] found that only fast aftershock decay occurs in areas in Japan with estimated temperature more than 400 o C, and cooler areas can have a range of decay rates. Wiemer and Katsumata [1999] found that p value is correlated with areas of high slip in the four aftershock zones they studied. They hypothesize that areas with more frictional heating produce higher p values. If steady state creep is a mechanism of stress decay after the stress step of a main shock, creep is empirically accelerated at higher temperatures [i.e., Jaeger and Cook, 1969, chapter 11]. [7] Stress decay must also be considered along with temperature when interpreting p values. Mikumo and Miyatake [1979, 1983] performed numerical experiments and concluded that the p value is larger (faster decay) for more rapidly relaxing stress relaxation times. This is in agreement with Dieterich s [1994, Figure 8] modeling of faster aftershock decay with higher rates of logarithmic stress decrease after the main shock. [8] Fault strength and stress heterogeneity are also important. Mikumo and Miyatake [1979] found that p is larger for a more homogeneous distribution of fault strength. Utsu [1961] also found that p value is larger for a more homogeneous distribution of shear strength on the fault, and for faster recovery of shear strength on the fault. This effect of homogeneity does not seem to be the dominant one in Hawaii because the strength of faults near active volcanic centers such as Kilauea and Mauna Loa is unlikely to be homogeneous given the common occurrence of intrusions, lava flows, and hydrothermal and magmatic fluids. Helmstetter and Shaw [2006] found that a heterogeneous distribution of stress change can lower the p value from the 1.0 value predicted for a uniform stress change, but did not indicate a way to increase p above 1.0 using Dieterich s [1994] rate-and-state friction law. [9] There is no good evidence that p varies strongly with depth. Davis and Frohlich [1991] found that p was lower for a set of intermediate depth and deep earthquakes than for a set of shallow ridge and transform earthquakes. Nyffeneggar and Frohlich [2000] found the opposite depth behavior: p was higher for two deep earthquakes than for a set of intermediate depth earthquakes, and that neither differed significantly from shallow sequences. [10] We seek physical interpretations of the Omori decay parameter p of Hawaiian sequences. We also use the modified Omori model to estimate aftershock productivity, and then interpret its variation. Finally, we suggest that properties of aftershock sequences may be used to reveal aspects of stress rates and to identify the presence of magma in other volcanic systems. 2. Methodology [11] We chose a model where aftershock rates follow the Gutenberg-Richter and modified Omori laws [e.g., Reasenberg and Jones, 1989]: Rt; ð MÞ ¼ 10 aþb ð Mm M Þ ðt þ cþ p ð1þ where R(t, M) is the aftershock rate at time t after the main shock, for magnitude M and greater. Here, M m is the main shock magnitude, a is the rate parameter, b is the magnitude parameter, c is the time delay parameter, and p is the aftershock decay parameter. [12] We use a practical definition of aftershocks. We use earthquakes within one week of the main shock to define the aftershock zone, and use this zone to select aftershocks and to determine the background earthquake rate. Aftershocks end as the Omori decay curve reaches a background level. We may thus have included a few background earthquakes in the aftershock counts, but our results do not change because the number of possible background events is small. Rate comparisons before and after main shocks were only done for larger and readily identifiable aftershock zones with lots of earthquakes. [13] We simplify the aftershock history by ignoring any spatial variability of aftershocks and ignoring lobes of stress change from the main shock. Most main shock focal mechanisms are unknown, and the number of aftershocks is generally too small to compare with calculated stress lobes. 2.1. Omori Parameters [14] We determined the modified Omori parameters a, b, c and p using Reasenberg s [1994] program aspar. The program allows the user to choose the minimum magnitude, and beginning and ending times after the main shock. The program determines the four parameters by a maximum likelihood technique based on the method of Ogata [1983]. Good parameter fitting requires a well-recorded sequence with a magnitude range of about 3 or 4 between the main shock magnitude and the completeness magnitude M c.we generally set the minimum magnitude M min for counting aftershocks to M c so the linear Gutenberg-Richter relation predicts the number of events we count. Because the main shock initially obscures the seismograms, or because a power failure may prevent recording immediately after the main shock, aftershocks are only used after recording is complete down to M c. Aftershocks are used before the rate 2of26

Figure 1. Aftershock rate curves for two large earthquakes. The M 7.2, 29 November 1975 Kalapana earthquake (top curves, open symbols) is the largest earthquake recorded by the modern seismic network. It is our standard earthquake because the aftershock parameters (a, b, c, and p) describing the Kalapana aftershock rates are used to calculate the expected aftershock numbers for other earthquakes. The south Hawaii M 5.2, 1 February 1994 event is the best recorded deep (35 km) earthquake (solid symbols). The lower 29 November 1975 curve estimates aftershock rates with the same magnitude spread DM = M m M min = 3.9 (between the main shock and aftershock cutoff magnitudes) as the deep 1994 event. In other words, the lower curve is shifted downward by a factor of 0.22 to the rate expected for M 3.3 aftershocks using equation (3). Identical magnitude spreads facilitate comparison between the aftershock rates of the two earthquakes. returns to the background rate, or until another main shock occurs. [15] Our standard aftershock sequence (Figure 1, open symbols) follows the Ms 7.2 Kalapana earthquake of 29 November 1975, which is the largest earthquake recorded by the modern Hawaii seismic network. Kalapana is a good standard event because its parameters are well determined and a majority of other sequences are in the same south flank region. We tried using different standard events other than the 1975 Kalapana earthquake, but found the general patterns of parameters, which depend on the choice of a standard event, did not change. The 29 November 1975 decay parameter p = 0.82 is typical of flank earthquakes in Hawaii. The upper curve (Figure 1) is for M 2.6 aftershocks which are complete as judged by linearity of the frequency-magnitude distribution. The rate curve shown by the line was fit using aftershocks between 0.25 and 100 days after the 29 November 1975 main shock. The first 4 hours of aftershocks are missing because of seismograph failure. After 100 days, the earthquake rate exceeded that extrapolated from the initial aftershock decay curve, and in this case approached a background level of about one earthquake per day. [16] We want to evaluate the number of aftershocks following different main shocks. Once a, b, c and p are determined for our standard aftershock sequence, we can estimate the expected number of recorded aftershocks N calc for any sequence by integrating equation (1) from the minimum magnitude M min to M m and over the time period of valid aftershock recording. The result uses only the magnitude spread DM =M m M min. When the number of aftershocks is large, a log(rate) versus log(t) plot similar to Figure 1 is how we choose the time interval that shows normal Omori decay. When the number of aftershocks is small, we choose a standard interval from 0.05 to 5 days following the main shock for counting and calculating aftershock numbers. We used Reasenberg s [1994] program enas to do the integration and estimate the expected number of aftershocks N calc. [17] An objective measure of aftershock productivity is a quantity we will use to compare different main shocks. We define aftershock productivity as A p =N obs /N calc, where N obs is a count of aftershocks larger than M c.a p can thus be determined for many small sequences which have only a few recorded aftershocks. [18] Aftershock productivity A p is sensitive to the main shock s magnitude. For example, underestimating the main shock magnitude by 0.3 can cause A p to increase twofold. Productivity values are also sensitive to variations in catalog completeness and magnitude irregularities. Because the ratio A p may have a large error, we plot only its logarithm and interpret only large variations in aftershock productivity, and seek values from many sequences to form a pattern. 2.2. Stress Parameters [19] Dieterich [1994] took a major step in practical seismology when he related changes in earthquake rate to changes in stress. He derived aftershock rates and obtained the modified Omori s law from a constitutive law with rate and state-dependent fault properties. Dieterich models an aftershock sequence resulting from a shear stress step Dt. Unlike the descriptive Omori law, Dieterich physically relates aftershock rates to stress parameters and includes a return to background rates if the stressing rate continues unchanged after the main shock. Dieterich [1994, equation [13]] expressed the earthquake rate increase above the background rate at the time of the main shock by a stress term R 0 =r ¼ e ðdt=asþ ð2þ R 0 is the aftershock rate immediately after the main shock at t = 0 but before the rate begins to decay. We use R 0 for earthquakes of the background minimum magnitude M c or greater; r is the background rate. A is a fault constitutive parameter, generally 0.005 to 0.012 [Dieterich, 1994, p. 2604], Dt is the shear stress change and s is the normal stress. This relation has been experimentally verified [e.g., Gross and Kisslinger, 1997]. We measure R 0 graphically from the aftershock rate curves like Figure 1 or 2a and 2b. Both R 0 and r must be referred to the same minimum completeness magnitude for counting earthquakes, but the completeness magnitudes before and after the main shock 3of26

the aftershock rate stops decaying and reaches a constant background rate. Like G, t a can be related to stress rates. t a is the fundamental relaxation time in Dieterich s [1994, equations [12] and [14]] rate equation t a ¼ As=_t ð5þ where _t is the stress rate after the main shock; _t thus has a direct effect on the level of postaftershock background seismicity. Note that G and t a are different quantities depending on stress rate before and after the main shock, respectively, that we will compare later in the paper. [21] A relation between stress decay and p value falls beautifully from Dieterich s [1994] rate equations: Dieterich s Figure 8 noted that Omori aftershock decay with different decay rates (p values) is consistent with different rates of logarithmic time decay of stress after a main shock. Figure 2a. Aftershock zones of three types of sequences with different decay rates and aftershock productivity. Each sequence uses a different plotting symbol. may not be the same; b 0 is the slope of the aftershock frequency-magnitude distribution whose rate we are adjusting. R 0 (for a standard background completeness magnitude M c ) can be calculated from R 0 0 (for the aftershock completeness magnitude M 0 c) using R 0 ¼ R 0 ðmc Mc 0 Þ 0 10 b0 ð3þ If you equate the modified Omori law (equation (1)) with Dieterich s [1994, equation [13]] rate equation for t c, other stress terms can be derived from measurable seismicity rate parameters aþb Mm Mc As=_t r ¼ 10 ð Þ =r ¼ G ð4þ where _t r is the background or reference stressing rate. We abbreviate the exponential expression of a, b, r and DM as the aftershock gain G, which has the units of days. G is a measure of the sustained aftershock rate relative to the premain shock rate r. Equation (4) relates an empirical but measurable rate expression to a stress expression based on Dieterich s [1994] rate-state friction constitutive law. Equation (4) allows inferences about stress rates to be made from aftershock earthquake rates. The abbreviation and definition of G is arbitrary but convenient. Dieterich [1994] assumes either the applied stressing rate after the main shock is 0 (we do not assume this), or the time for which a and b are measured is before the aftershock rate begins to return to background levels (we do assume this). [20] The aftershock duration t a and aftershock gain G are useful parameters because they show the stress rate changes of main shocks in different situations. We measure t a directly from the aftershock decay curve as the time when Figure 2b. Aftershock rate curves for three main shocks illustrating different types of aftershock sequences. (1) The M 6.7, 16 November 1983 earthquake in the Kaoiki seismic zone is a large flank earthquake with a typical 0.75 p value (open triangles, point down). (2) The 10 August 1981 sequence followed an M 4.3 earthquake in Kilauea s south flank, which was triggered by an intrusion into Kilauea s adjacent southwest rift zone (open triangles, point up). The 0.80 p value of this intrusion sequence is comparable to other flank earthquakes. (3) The aftershock decay rate (p = 2.67) of the M 4.3, 24 January 1993 event near Kilauea Caldera (solid triangles), however, is much higher than the others. The three curves are normalized to the same magnitude spread DM =M m M min = 3.2 between the main shock and aftershock cutoff magnitudes to facilitate rate comparison. The caldera and intrusion aftershock sequences are considerably more productive than the aftershocks of the 16 November 1983 earthquake, which had no known immediate triggering event. 4of26

Dieterich added a logarithmic term to the stress history applied to a region after the main shock s stress step t 0 : tðþ¼t t 0 þ u lnðwt þ 1Þ ð6þ If u is 0 or positive (applied stress is constant or increases with the logarithm of time), aftershocks decay with p = 1.0 [Dieterich, 1994]. Aftershocks decay with a range of p values above and below 1.0 for the case of a negative u (post-main shock stress decreases with the logarithm of time). Empirically, Dieterich s figure yields u p/5 (for negative u) after the characteristic time w 1, for Dieterich s choice of A, s, and Dt. Thus the faster the stress decay after the main shock, the more rapid the aftershock decay. Dieterich thus provides a physical link between p value and stress decay. [22] The time parameter c in Omori s law can also be represented using Dieterich s [1994, equations [15] and [16]] rate equation, where we define c 0 as this estimate of c c 0 ¼ ðas=_t r Þe ð Dt=AsÞ aþb Mm Mc ¼ 10 ð Þ =R 0 ¼ Gr=R 0 ð7þ Thus c 0 can be estimated from the initial aftershock rate R 0 and the aftershock frequency-magnitude relation. Dieterich [1994, equation [18]] refers to c 0 as t e, the time for the rate to merge with the 1/t asymptote. Our quantity c, on the other hand, is determined by fitting the aftershock rate curve directly. The c in the Omori relation is often difficult to estimate from the shape of the aftershock decay curve because it depends on measuring the changes in the aftershock rate immediately after the main shock. The c 0 is also difficult to estimate because it depends on measuring the rate R 0 in the first hour or hours after the main shock. [23] Stress change Dt, background stress rate _t r, and c 0 value are determined from aftershock rates and are all related. Dt/As depends on aftershock rate early in the sequence (equation (2)), and As/_t r depends on sustained rates (equation (4)). It is therefore sensible that Dt/As and As/_t r should be highly correlated. Transforming equation (7) implies there is a linear relationship if the Omori time delay c 0 is a constant Dt=As ¼ lnðas=_t r Þ lnðc 0 Þ We will examine this relation for Hawaiian earthquakes in more detail later. 2.3. Main Shock and Aftershock Selection [24] We use aftershock sequences from 1960 to 2001 selected from the catalog of the Hawaiian Volcano Observatory (2002, October 1959 through 2001 Hawaiian earthquake catalog, unpublished annual computer files available at ftp computer sites, through the Council of the National Seismic System and other sources), supplemented with a few well-recorded 1868 1951 aftershock sequences from the historical catalog of Klein and Wright [2000]. Moment magnitudes are only available for a very few Hawaiian earthquakes, and we rely mostly on local magnitudes determined with a Wood-Anderson seismometer, teleseismic surface wave magnitudes, or on scales calibrated against them. ð8þ [25] We included as large a selection of main shocks as we could. All M 4.8 earthquakes, and many smaller earthquakes that produced at least two aftershocks larger than their completeness magnitude M c, were examined as potential main shocks. We excluded offshore events with poor network coverage. We also excluded swarms at Loihi submarine volcano, which has no aftershocks with Omori decay. The lack of Omori decay suggests that stresses at Loihi are largely magmatic, with rapid variation and no gradual stress buildup that can cause isolated main shocks. Table 1 lists the essential parameters of all the main shocks we used. [26] Of special interest in our aftershock study are triggered main shocks, namely those with an identifiable causative event stressing the main shock area within minutes, hours or days before the main shock. By this definition, triggered main shocks include flank earthquakes adjacent to an active rift intrusion, and larger aftershocks, which may become main shocks if they have their own secondary aftershocks. Many earthquakes may be triggered by a sudden stress increase whose origin is unknown, and we will examine earthquakes where we did not class them as triggered even though they behave like earthquakes with a known triggering event. We searched for and found many triggered sequences but used them only if they had Omori time decay. The flank earthquakes triggered by large intrusions such as December 1974 and June 1982 had several episodes of high swarm activity but had Omori decay only at the end of the sequence when the intrusion had largely ended. Klein et al. [1987] has time plots of the seismicity of these and other intrusions that include Omori aftershock sequences. [27] We also studied aftershocks of the great M 7.9 Kau earthquake of 2 April 1868 [Wyss, 1988] using the catalog of Klein and Wright [2000]. Including aftershocks of an earthquake this old is remarkable. An excellent diary of Hawaiian earthquakes felt at Hilo, with enough description to determine intensities, exists for 1833 1917 [Wyss et al., 1992b]. We also determined the background seismicity rate before and after the Kau aftershock sequence. The median magnitude corresponding to maximum intensity V earthquakes in this era is 5.3 [Klein and Wright, 2000]. The completeness magnitude of the post-1840 seismicity is about 5.2 because of linearity of the frequency-magnitude relation. Thus the record of intensity V (M 5.3) earthquakes felt at Hilo is probably complete. The Kau earthquake aftershock catalog has 239 M 5.2 events in 19 years, after which rapid Omori decay ends and seismicity resumes a more constant rate. Because the 1868 rupture zone includes most of south Hawaii and coincides with the contemporary active seismic zones, we assume that all intensity V or greater earthquakes felt in Hilo or south Hawaii are within the Kau rupture zone. 3. Observations 3.1. Aftershock Decay Parameter p [28] Our study compiled 40 aftershock p values of good quality rated A, B or C (Table 1; details are in given Table 2). C-rated sequences are included in plots because their p values have the same geographic correlations as A and B rated sequences and thus reinforce their behavior. 5of26

Table 1. Essential Parameters of Hawaiian Main Shocks a Latitude Longitude Date Time, LT deg min deg min Depth, km Location M m M min p Value and Quality c Ap Notes b Shallow Earthquakes 28 Mar 1868 1328 19 6.00 155 39.00 10 Hilea 7.0I 5.2 0.84 ± 0.57 d- - (43) int 2 2 Apr 1868 1600 19 12.00 155 30.00 10 Hilea 7.9I 5.2 0.77 ± 0.09 a (0.41) (6.1) int 2,3 25 Oct 1913 0057 19 21.00 155 1.00 10 S. flank 5.8H 3 0.71 ± 0.17 b - - 5 Jul 1914 1918 19 26.00 155 24.00 10 Kaoiki 4.8N 2 0.77 ± 0.22 c - - 5 Oct 1929 2122 19 43.00 156 5.00 10 Hualalai 6.5S 3 1.95 ± 0.73 c- - - 16 Jun 1940 2356 21 0.00 155 18.00 10 Maui 6.0S 3.5 0.89 ± 0.12 b - - 21 Feb 1942 0811 19 32.00 155 28.00 8.00 M.Loa 6.1H 3 0.89 ± 0.23 c - - 27 Apr 1942 2143 19 32.00 155 28.00 5.00 M.Loa 6.1A 3 1.32 ± 0.29 c - - 21 Aug 1951 0057 19 30.00 155 57.00 10 Kona 6.9S 3.0 0.91 ± 0.05 a 0.02 0.93 30 Mar 1954 0842 19 21.00 155 0.00 10.00 S. flank 6.4H 2.6 - - 0.18 27 Jun 1962 1827 19 23.89 155 27.10 10.18 Kaoiki 6.1L 2.3 (0.9) (0.3) (2.2) 1 9 May 1969 1533 19 21.61 155 4.57 12.70 S. flank 4.2L 1.9 1.08 ± 0.36 D - 7.6 8 Oct 1969 1417 19 13.82 155 20.97 0.01 SW.flank 3.4L 1.9 1.93 ± 0.61 D - 114 int 9 Nov 1969 1912 19 11.20 155 32.38 9.64 Hilea 4.5L 2.1 - - 1.43 24 Nov 1969 0912 19 44.18 156 5.75 2.27 Hualalai 4.6L 2.6 - - 2.7 12 Apr 1970 0913 19 23.84 155 26.18 9.78 Kaoiki 4.3L 2.0 - - 6.4 21 Sep 1970 0126 19 19.93 155 12.18 10.82 S. flank 4.5L 1.9 1.05 ± 0.56 D - 2.7 30 Sep 1971 2128 19 16.11 155 20.74 0.95 SW.flank 4.2L 2.0 0.97 ± 0.24 B 0.05 17 int 9 Dec 1971 0215 19 20.17 155 6.72 8.18 S. flank 4.3L 1.9 0.73 ± 0.25 D - 2.9 29 Dec 1971 0042 19 15.18 155 22.27 6.39 SW.flank 4.3L 1.9 1.09 ± 0.05 A 0.02 31 int 29 Feb 1972 1208 19 21.60 156 20.48 6.43 Kona 5.0L 2.6 - - 1.11 31 Mar 1972 1620 19 20.12 155 3.55 9.70 S. flank 4.5L 1.9 - - 1.15 5 Sep 1972 0131 19 19.69 155 12.35 10.14 S. flank 5.2L 1.9 - - 0.84 23 Dec 1972 0904 19 34.83 155 57.09 15.97 Kona 5.1L 2.8 - - 1.44 15 Apr 1973 0059 19 19.39 155 7.06 9.89 S. flank 4.5L 2.1 - - 2.9 12 Jan 1974 0604 19 19.83 155 7.25 8.87 S. flank 4.7L 2.0 0.62 ± 0.17 C - 2.1 19 Jun 1974 0505 19 22.78 155 25.32 10.36 Kaoiki 4.6L 1.9 0.95 ± 0.09 A 0.03 27 20 Jun 1974 2050 19 19.70 155 12.60 9.50 S. flank 4.4L 1.8 - - (0.81) 2 27 Aug 1974 2149 19 19.59 155 12.31 10.28 S. flank 4.4L 1.8 - - 1.63 30 Nov 1974 0354 19 26.28 155 25.17 6.09 Kaoiki 5.5L 2.0 0.95 ± 0.09 A 0.06 13.1 15 Dec 1974 2317 19 24.35 155 26.00 9.01 Kaoiki 4.5L 1.9 0.69 ± 0.07 A 0.01 35 aft 25 Dec 1974 1813 19 13.87 155 18.11 9.72 S. flank 4.2L 1.4 0.78 ± 0.23 C - 0.58 31 Dec 1974 1240 19 17.79 155 21.74 2.62 SW.flank 5.6L 2.0 (1.27 ± 0.05) 1.42 17 int 4 Jan 1975 1532 19 14.79 155 22.31 5.58 SW.flank 4.8L 2.0 0.86 ± 0.10 A 0.01 19.5 int 21 May 1975 2232 20 17.42 155 39.33 11.81 Kohala 4.7L 2.8 - - (0) 2 9 Jul 1975 0840 19 31.86 155 28.50 7.15 M.Loa 4.6L 2.0 1.28 ± 0.25 B 0.06 22 int 10 Nov 1975 0126 19 21.26 155 1.58 9.79 S. flank 4.5L 2.1 - - 2.3 15 Nov 1975 1255 19 18.69 155 13.51 10.74 S. flank 4.5L 1.4 - - 3.8 29 Nov 1975 0447 19 20.34 155 0.26 9.14 S. flank 7.2S 2.6 0.82 ± 0.06 A (0.43) 0.84 3 29 Nov 1975 2015 19 25.16 155 22.37 11.75 Kaoiki 4.6L 2.5 0.49 ± 0.14 B 0.17 8.2 trig 29 Jan 1976 1019 19 22.20 154 58.90 9.91 S. flank 4.5L 2.5 - - (0) aft 20 Feb 1976 1951 20 23.16 156 3.36 6.93 Kohala 5.1L 2.9 - - (0) trig 24 Feb 1976 0550 19 21.70 155 6.40 9.30 S. flank 4.3L 2.6 - - 24 aft 2 Apr 1976 0814 19 20.36 155 6.34 9.91 S. flank 4.6L 2.4 - - 1.7 aft 4 Jun 1976 2250 19 21.27 155 6.97 9.88 S. flank 4.3L 2.3 - - 8.7 aft 18 Dec 1976 0401 19 19.75 155 6.89 9.70 S. flank 5.0L 2.0 0.95 ± 0.40 D - 1.96 14 Jan 1977 1326 19 19.66 155 7.19 9.87 S. flank 4.8L 2.0 - - 1.21 3 Feb 1977 1520 19 20.75 155 4.44 9.93 S. flank 4.5L 2.1 - - 2.4 20 Apr 1977 1849 19 55.95 155 19.69 12.97 M.Kea 4.8L 2.4 0.82 ± 0.17 C - 1.54 5 Jun 1977 2342 19 21.61 155 4.88 9.52 S. flank 5.3L 2.1 - - 0.94 19 Aug 1977 0819 19 19.73 155 7.09 10.22 S. flank 4.5L 2.0 - - (0.76) 2 15 Sep 1977 1850 19 20.33 155 3.59 9.42 S. flank 4.2L 2.0 0.80 ± 0.11 A 0.82 83 int 23 Sep 1977 0208 19 20.96 155 2.70 8.60 S. flank 4.1L 2.0 - - 32 int 23 Jun 1978 0147 19 19.08 155 15.46 10.37 S. flank 4.5L 2.1 - - 2.3 1 Jul 1978 0918 19 18.39 155 6.74 7.89 S. flank 4.5L 1.95 - - 2.9 11 Sep 1978 2016 19 19.87 155 6.52 9.60 S. flank 4.5L 2.0 - - 2.4 14 Dec 1978 0412 19 18.60 155 13.53 10.36 S. flank 4.4L 1.7 - - 1.82 27 Dec 1978 0040 19 20.06 155 12.95 9.81 S. flank 4.4L 1.7 - - (1.32) 2 10 Mar 1979 0355 19 19.95 155 6.69 9.56 S. flank 4.7L 2.0 - - 2.1 21 Mar 1979 2047 20 1.66 155 48.18 13.57 Kohala 4.7L 2.0 0.88 ± 0.11 B 0.01 1.8 27 Mar 1979 2130 20 0.44 155 46.92 10.66 Kohala 5.1L 2.0 - - 0.60 31 Jul 1979 0330 19 28.13 155 25.81 11.86 Kaoiki 4.5L 1.8 - - 1.32 21 Sep 1979 2159 19 20.65 155 4.23 9.28 S. flank 5.7L 1.7 0.83 ± 0.06 A 0.01 2.2 27 Sep 1979 0535 19 19.66 155 7.25 9.96 S. flank 4.6L 1.95 - - 4.1 13 Dec 1979 1744 19 24.82 155 24.47 11.36 Kaoiki 4.3L 2.1 - - (0.72) 2 12 Mar 1980 0257 19 21.45 155 14.21 1.85 Koae 4.3L 1.7 - - 4.2 int 11 Aug 1980 2042 19 19.91 155 6.23 10.09 S. flank 4.4L 2.0 - - 10 int 1 Mar 1981 0701 19 21.46 155 2.05 9.20 S. flank 4.7L 2.0 - - 1.8 28 Jul 1981 1000 19 21.45 155 1.60 8.70 S. flank 4.4L 2.0 - - (0.57) 2 10 Aug 1981 0820 19 19.02 155 21.03 2.02 SW.flank 4.3L 1.8 0.80 ± 0.06 A 0.34 11.4 int 22 Aug 1981 1205 20 11.25 156 25.05 9.43 Hual.OS 4.5L 2.9 - - (0) 2 6of26

Table 1. (continued) Latitude Longitude Date Time, LT deg min deg min Depth, km Location M m M min p Value and Quality c Ap Notes b 10 Nov 1981 0302 19 20.37 155 12.71 10.17 S. flank 4.5L 1.8 1.00 ± 0.27 C - 2.7 21 Jan 1982 1152 19 13.61 155 35.57 10.14 Hilea 5.6L 1.8 1.06 ± 0.04 A 0.03 2.5 11 Apr 1982 1604 19 19.71 155 6.69 9.24 S. flank 4.5L 1.8 - - 0.61 14 May 1982 0626 19 59.90 155 51.78 18.69 Kohala 4.8L 2.7 - - (1.07) 2 18 May 1982 1736 19 57.37 156 25.57 0.87 Hual.OS 4.8L 2.9 - - (0) 2 24 Jun 1982 2059 19 16.86 155 21.65 7.13 SW.flank 3.0L 1.8 1.33 ± 0.39 C - (433) int 2 12 Aug 1982 0043 19 24.89 155 16.08 16.20 K.caldera 4.3L 1.9 0.82 ± 0.18 B 0.02 2.9 8 Mar 1983 0641 19 11.89 155 35.64 11.49 Hilea 4.6L 2.2 - - 1.42 20 Mar 1983 1718 19 21.51 155 3.04 6.95 S. flank 4.9L 1.8 - - 1.8 20 Mar 1983 2302 19 21.88 155 25.08 11.12 Kaoiki 4.0L 1.8 - - 2.9 27 Apr 1983 2334 19 19.74 155 7.59 8.35 S. flank 4.3L 1.8 - - 1.8 13 May 1983 0030 19 10.12 155 34.83 9.34 Hilea 4.4L 2.0 - - (0.38) 2 9 Sep 1983 0630 19 19.79 155 7.32 8.96 S. flank 5.7L 1.8 - - 0.55 16 Nov 1983 0613 19 25.76 155 27.12 10.97 Kaoiki 6.7S 2.1 0.75 ± 0.02 A (0.17) 1.00 3 21 Feb 1985 1948 19 19.61 155 12.66 9.37 S. flank 4.8L 1.7 0.54 ± 0.11 B 0.01 1.8 7 Jul 1985 0200 19 9.97 155 35.66 10.28 Hilea 4.5L 1.8 1.57 ± 0.37 C - 2.1 12 Dec 1985 0901 20 34.97 155 45.11 23.64 Kohala 5.0L 2.6 - - (0) 2 6 Apr 1986 2237 19 12.18 155 36.86 8.42 Hilea 4.4L 1.8 - - 0.76 9 Jul 1986 0228 19 31.74 155 56.03 12.59 Kona 4.4L 2.1 - - 2.16 3 Feb 1987 1622 20 4.42 156 25.25 0.05 Hual.OS 5.2L 2.7 1.02 ± 0.07 B 0.02 12 25 Nov 1987 1849 20 20.60 156 14.35 17.20 Maui 4.4L 2.6 - - (1.39) 2 19 Feb 1988 1847 19 21.46 155 1.67 8.70 S. flank 4.2L 1.9 - - 2.2 1 Mar 1988 2241 19 19.60 155 12.50 10.13 S. flank 4.9L 1.8 - - 0.39 24 Mar 1988 1429 19 57.00 156 23.91 2.00 Hual.OS 5.0L 2.2 1.21 ± 0.45 D - 2.4 27 Mar 1988 1733 19 56.56 156 24.18 2.29 Hual.OS 5.1L 3.0 0.93 ± 0.07 B 0.02 22 1 Apr 1988 1848 19 54.92 156 22.99 31.46 Hual.OS 4.9L 2.7 - - 5.6 11 May 1988 1314 19 47.04 155 31.30 23.03 M.Kea 4.3L 1.9 0.92 ±0.24 D- - 1.8 3 Jul 1988 1938 19 12.81 155 27.31 9.47 Hilea 5.4L 1.8 0.65 ± 0.08 B 0.01 0.35 13 Aug 1988 1620 20 12.18 156 29.66 7.72 Hual.OS 4.5L 3.0 - - (0) 2 25 Jun 1989 1727 19 21.53 155 5.01 9.30 S. flank 6.2L 2.0 - - 1.43 27 Dec 1989 2313 19 19.60 155 12.41 9.42 S. flank 5.3L 1.6 0.62 ±0.54 D- - 0.41 1 Aug 1990 1937 19 49.55 155 37.29 18.36 M.Kea 4.7L 2.1 0.80 ± 0.08 B 0.01 2.3 8 Aug 1990 1606 19 20.02 155 6.76 9.21 S. flank 4.8L 1.8 - - 0.81 24 Jan 1993 2214 19 25.34 155 19.20 6.01 K.caldera 4.3D 1.1 2.67 ± 0.49 A 0.31 12.3 26 Jan 1993 0524 19 13.47 155 29.73 9.43 Hilea 4.8D 1.4 1.35 ± 0.17 B 0.03 0.39 8 Jun 1993 0257 19 20.61 155 12.88 9.85 S. flank 5.3D 1.3 - - 0.28 19 Mar 1995 2229 20 3.13 156 34.52 0.41 Hual.OS 4.3D 3.0 - - (0) 2 11 May 1995 0349 20 4.05 155 20.97 5.96 M.Kea 4.8D 1.8 - - 0.23 21 Jan 1996 1109 19 50.83 155 31.37 20.71 M.Kea 4.4D 1.6 1.27 ± 0.31 C - 2.2 18 Jul 1996 0739 19 53.61 155 35.02 14.70 M.Kea 4.2D 1.6 - - (0.38) 2 23 Nov 1996 1639 19 19.90 155 12.34 10.50 S. flank 4.3D 1.5 - - 1.46 30 Jun 1997 0547 19 21.13 155 4.02 9.16 S. flank 5.5D 1.5 0.96 ± 0.19 B 0.01 0.38 14 Aug 1997 1554 19 20.03 155 6.90 8.93 S. flank 5.0D 1.3 1.14 ± 0.45 D - 0.55 27 Sep 1998 2156 19 26.28 155 13.59 0.80 K.caldera 4.6U 1.8 1.44 ± 0.39 C - 1.46 28 Sep 1998 2039 19 20.78 155 7.59 9.52 S. flank 4.8U 1.3 0.86 ± 0.27 B 0.03 1.06 22 Nov 1998 0554 20 24.43 156 4.17 27.71 Kohala 4.5U 2.2 - - (0.72) 2 16 Apr 1999 1456 19 15.34 155 29.38 9.37 Hilea 5.6U 1.1 0.96 ± 0.07 A 0.01 0.71 26 May 1999 0601 19 25.21 155 19.50 7.64 K.caldera 4.3U 1.5 1.33 ± 0.20 B 0.02 3.9 3 Jun 1999 0141 19 58.56 155 33.18 27.67 M.Kea 4.4U 1.8 - - 0 16 Aug 1999 1405 19 16.85 155 30.12 8.97 Hilea 4.4U 1.8 - - 0.74 1 Apr 2000 2018 19 20.73 155 12.50 9.46 S. flank 5.0U 1.0 - - 0.26 25 Apr 2001 1737 19 25.44 155 18.28 6.34 K.caldera 4.4U 1.0 - - 0.60 Earthquakes Known to Be Deeper Than 20 km 23 Jul 1961 0528 19 23.60 155 17.95 23.78 D.Kilauea 4.4L 2.6 - - 23 25 Aug 1961 0845 19 49.84 155 5.20 43.53 M.Kea 4.6L 3.2 - - (0) 2 11 May 1971 1600 18 57.14 155 33.00 37.08 D.SW.rift 4.7L 2.4 - - 1.8 15 Aug 1971 1536 19 21.95 155 16.70 34.04 D.Kilauea 4.9L 2.3 - - 1.1 22 Apr 1973 2107 20 1.33 154 35.66 33.72 offshore 5.0L 3.5 - - (0) 2 26 Apr 1973 1026 19 51.92 155 9.17 38.71 M.Kea 6.2S 2.4 1.11 ± 0.14 A 0.29 0.71 9 Oct 1973 0153 19 20.24 155 16.07 32.50 D.Kilauea 4.8L 2.3 - - 4.7 13 Dec 1973 0425 19 22.44 155 17.54 34.84 D.Kilauea 4.8L 2.3 - - 1.5 25 Dec 1974 0747 19 20.86 155 16.84 32.30 D.Kilauea 4.8L 2.0 0.97 ± 0.27 C - 2.1 6 Nov 1975 0205 19 20.56 155 18.83 31.90 D.Kilauea 4.6L 2.0 0.94 ± 0.19 B - 2.2 7 Sep 1977 1351 19 22.37 155 19.34 31.42 D.M.Loa 4.6L 2.1 - - 2.1 31 Aug 1978 1307 18 59.97 155 28.97 35.25 D.SW.rift 4.5L 2.2 - - 2.1 6 Mar 1979 0507 19 31.23 155 16.20 27.53 D.Kilauea 4.8L 2.1 - - 0.44 19 Jan 1980 1528 19 18.67 155 32.41 26.84 D.SW.rift 4.5L 2.2 - - 0 12 Jan 1981 0418 19 21.33 155 18.28 31.19 D.Kilauea 4.9L 2.1 - - 1.5 15 Mar 1981 2017 19 22.42 155 14.03 31.57 D.Kilauea 4.4L 2.1 - - 4.3 7 Feb 1983 1602 19 21.44 155 14.46 28.08 D.Kilauea 4.3L 2.1 - - 2.6 7of26

Table 1. (continued) Latitude Longitude Date Time, LT deg min deg min Depth, km Location M m M min p Value and Quality c Ap Notes b 30 Jun 1985 1112 19 22.39 155 17.87 26.89 D.Kilauea 4.5L 2.1 - - 0 22 Apr 1986 1843 19 17.89 155 16.24 32.14 D.Kilauea 4.4L 2.2 - - 3.5 19 Sep 1986 0444 19 20.20 155 21.10 30.82 D.Kilauea 4.2L 1.9 - - 4.3 17 Jul 1988 1725 19 15.46 155 16.18 31.52 D.Kilauea 4.3L 1.8 - - 0.95 8 May 1991 0821 19 20.49 156 13.30 30.50 Kona OS 5.5L 2.7 - - 0.49 1 Feb 1994 0001 19 15.39 155 17.99 34.73 D.Kilauea 5.2S 1.3 1.39 ± 0.12 A 0.12 1.3 21 Jun 1998 0221 19 12.65 155 20.64 46.48 D.Kilauea 4.2U 1.9 - - 1.4 17 Feb 2000 1419 19 20.21 155 16.76 35.29 D.Kilauea 4.5U 1.4 0.84 ± 0.15 C - 0.65 a Date and time are Hawaii Standard Time as quoted from original sources and computer files. Depth is depth below the epicentral ground surface. Earthquakes in northern Hawaii with depths between 20 and 30 km are grouped with the shallow earthquakes because their depth distribution and time behavior suggest they are not distinguishable from earthquakes shallower than 20 km. Location is earthquake source region determined by Klein and Wright [2000] or by computer location. Abbreviations are D.Kilauea, deep Kilauea caldera; D.M.Loa, deep Mauna Loa volcano; D.SW.rift, deep zone under Kilauea s SW rift zone; Hilea, Hilea seismic zone in Mauna Loa s south flank; Hualalai, Hualalai volcano area; Hual.OS, Hualalai offshore (west rift zone); K.caldera, Kilauea caldera; Kaoiki, Kaoiki seismic zone in Mauna Loa s SW flank; Koae, Koae fault zone immediately south of Kilauea caldera; Kohala, Kohala volcano; Kona, Kona coast (west coast of Hawaii); Maui, Maui island region; M.Kea, Mauna Kea volcano; M.Loa, Mauna Loa volcano caldera or rift zones; S. flank, south flank of Kilauea; SW.flank, west part of south flank of Kilauea. M m is main shock magnitude. The magnitude is chosen as the best type available from either Klein and Wright [2000] or the modern HVO catalog. The codes are I, intensity areas from isoseismal map; H, Honolulu seismogram amplitude reading; N, nomogram based size class from HVO Bosch-Omori; S, surface wave amplitude; A, average of 2 or more types; L, local from Hilo Wood Anderson; D, coda duration; U, local based on amplitudes from analog stations with approximate Wood-Anderson response. M min is minimum (completeness) magnitude of aftershock sequence. The p value is the decay parameter of aftershock rate expression. The letter (A, B, C etc.) is a quality rating of the p value based on the calculated error in the p value and the number N of aftershocks. Quality letters in lowercase indicate that p values meet the same criteria but are based on the early catalog with larger magnitude and completeness uncertainty: A, p(error)/p < 0.3 and N > 100; B, p(error)/p < 0.3 and N > 20; C, p(error)/p < 0.3 and N > 8; D, p(error)/p < 0.6 and N > 8 D quality sequences are judged too poor to plot or interpret. C is time delay parameter of aftershock rate expression derived from fitting the rate curve; c is usually not well determined. When c values are in parentheses, the early aftershock recording is incomplete and the c value is very poor. Ap is aftershock productivity defined as N OBS /N CALC, where N CALC is the expected number of aftershocks calculated from the a, b, c, and p. Parameters of the standard 29 November 1975 sequence. b Notes are int, preceded by a rift intrusion; aft, an aftershock sequence following an aftershock; trig: triggered by the 29 November 1975 main shock. (1) Before 1968, earthquakes were recorded at HVO on smoked paper. It was common practice to manually advance the pen carriage every revolution of the drum to prevent trace overlap during times of intense seismicity [see Koyanagi et al., 1966, Figure 4]. This was done in the first day or so of the aftershock sequence and resulted in more complete recording and catalog completeness in the first hours than in the following days. Thus the catalog of located earthquakes does not determine a good p value. The tabulated p value is from Koyanagi et al. [1966], who used a graphical fit to a log-log plot of aftershock numbers from 1 to 70 days after the main shock. They also determined the c value from a linear plot of M 0.5 aftershocks. (2) The aftershock productivity ratio is not plotted or interpreted because there was not a large enough main shock (as the leading event in an Omori decay period during a swarm) for a meaningful ratio; the observed number of aftershocks (0-1) was too small for a meaningful ratio; or the 1868 earthquake aftershocks in the early days of the sequence have imprecise magnitude estimates and catalog completeness. (3) First few hours of aftershocks not recorded, c value probably incorrect. D-rated p values are not plotted or interpreted. The older 1868 1959 catalog requires special treatment. Lowercase ratings a d are for older sequences with larger magnitude errors. We only use the two a-rated p values from the M 7.9, 2 April 1868 Kau earthquake, and the M 6.9, 21 August 1951 Kona earthquake in this study. [29] Two important Hawaiian aftershock sequences have contrasting decay curves: the 1 February 1994 deep (35 km) M 5.2 earthquake (solid symbols, Figure 1) decays faster than the 29 November 1975 reference earthquake located beneath Kilauea s south flank (open symbols). We must use the same magnitude spread DM =M m M c to compare the two aftershock curves. The magnitude spread of the 1994 sequence is DM = 5.2 1.3 = 3.9. Rather than decimating the 1975 aftershocks with an M min = 3.3 cutoff to get the same magnitude spread, we shift the M 2.6 curve downward by a factor of 0.22 to the rate expected for M 3.3 aftershocks using equation (3) (open symbols, lower curve, Figure 1). The deep 1 February 1994 earthquake thus has a higher initial aftershock rate in the first day, but a faster decay rate. [30] Examples of three aftershock sequences with contrasting p values are in Figures 2a and 2b. The M 6.7, 16 November 1983 earthquake in the Kaoiki seismic zone, like the M 7.2 November 1975 earthquake (Figure 1), is an example of a large flank earthquake with a well-recorded aftershock sequence. The 10 August 1981 sequence follows an M 4.3 earthquake in Kilauea s south flank that was triggered by an intrusion into Kilauea s adjacent southwest rift zone. The p value of this intrusion sequence (0.80) is comparable to other isolated flank earthquakes. The aftershock decay rate of the M 4.3, 24 January 1993 event near Kilauea caldera, however, is much higher. Its 2.67 ± 0.49 p value is the highest we found. [31] Earthquake sequences of different types have similar p values. Three earthquake classes, shallow (less than 20 km depth), deep, and those triggered by an intrusion or a primary main shock, all have similar peaked p distributions with p mostly in the range 0.7 to 1.1. The p distributions of these three classes have means of 1.03 ± 0.39, 1.05 ± 0.19 and 0.83 ± 0.24 (all 1s errors) respectively. Thus there is no apparent p value dependence on whether the main shock was triggered, and no dependence on depth. 3.2. Aftershock Productivity [32] We estimate the productivity (relative to the standard event) of 130 Hawaiian aftershock sequences with the parameter A p =N obs /N calc (Table 1). Table 3 lists detailed parameters used in the calculation. For example, we counted 72 aftershocks of the 1 February 1994 earthquake from 0.0125 to 125 days after the main shock (Figure 1) compared to an expected number of 54, for a productivity of A p = 1.3. Note that a productivity calculation of the standard 29 November 1975 earthquake will not always yield exactly 1.0. The productivity calculation uses the parameters a, b, c and p of the standard event, plus the magnitude spread and time window of each sequence, and 8of26

Table 2. Time Parameters of Aftershock Sequences a Date Where Depth M m M min N p a b c t 1 t 2 Notes b Shallow (z < 20 km) 28 Mar 1868 Hilea 10? 7.0 5.0 25 0.84 ± 0.57 d- - - - 0.04 5 19 int 2 Apr 1868 Hilea 10? 7.9 5.2 239 0.77 ± 0.09 a 3.51 1.11 (0.41) 0.16 7000 19 int 25 Oct 1913 S. flank 10? 5.8 3 47 0.71 ± 0.17 b - - - 1.0 200 20 5 Jul 1914 Kaoiki 10? 4.8 2 18 0.77 ± 0.22 c - - - 0.1 31 20 5 Oct 1929 Hualalai 10? 6.5 3 37 1.95 ± 0.73 c- - - - 8 1000 20, 23 16 Jun 1940 Maui 10? 6.0 3.5 27 0.89 ± 0.12 b - - - 0.01 140 20, 21 21 Feb 1942 M. Loa 10? 6.1 3 16 0.89 ± 0.23 c - - - 0.01 20 20, 22 27 Apr 1942 M. Loa 10? 6.1 3 13 1.32 ± 0.29 c - - - 0.01 50 20, 22 21 Aug 1951 Kona 10? 6.9 3.0 116 0.91 ± 0.05 a 2.84 1.02 0.02 0.016 250 20 27 Jun 1962 Kaoiki 10.2 6.1 2.3 53 0.9-0.98 0.3 1 70 1 9 May 1969 S. flank 12.7 4.3 1.9 18 1.08 ± 0.36 D - - - 0.016 9 8 Oct 1969 SW flank 1 9 3.4 1.9 45 1.93 ± 0.61 D - - - 0.10 31 18 int 21 Sep 1970 S. flank 10.8 4.5 1.9 9 1.05 ± 0.56 D - - - 0.01 8 30 Sep 1971 SW flank 1.0 4.2 2.0 25 0.97 ± 0.24 B 2.58 1.54 0.05 0.02 33 15 int 9 Dec 1971 S. flank 8.2 4.3 1.9 17 0.73 ± 0.25 D - - - 0.016 35 29 Dec 1971 SW flank 6.5 4.5 1.9 128 1.09 ± 0.05 A 1.14 0.98 0.02 0.016 700 16 int 12 Jan 1974 S. flank 8.0 4.7 2.0 17 0.62 ± 0.17 C - - - 0.01 16 19 Jun 1974 Kaoiki 10.4 4.6 1.9 140 0.95 ± 0.09 A 2.32 1.39 0.03 0.01 17 30 Nov 1974 Kaoiki 6.1 5.5 2.0 300 0.95 ± 0.09 A 1.74 1.03 0.06 0.016 8 15 Dec 1974 Kaoiki 9.0 4.5 1.9 156 0.69 ± 0.07 A 1.22 1.03 0.01 0.01 14 2 aft 25 Dec 1974 So Coast 9.7 4.4 1.4 10 0.78 ± 0.23 C - - - 0.06 630 3 31 Dec 1974 SW flank 2.6 5.6 2.0 880 (1.27 ± 0.05) 0.69 0.91 1.42 0.01 250 4 int 4 Jan 1975 SW flank 5.6 4.8 2.0 108 0.86 ± 0.10 A 1.11 0.90 0.01 0.01 8.5 2,4 int 9 Jul 1975 M. Loa 7.2 4.6 2.0 84 1.28 ± 0.25 B 1.25 1.09 0.71 0.06 35 11 int 29 Nov 1975 Kalapana 9.2 7.2 2.6 401 0.82 ± 0.06 A 2.42 0.93 (0.43) 0.25 100 29 Nov 1975 Kaoiki 11.8 4.6 2.5 46 0.49 ± 0.14 B 1.99 1.06 0.17 0.10 280 10 trig 18 Dec 1976 S. flank 9.7 5.0 2.0 16 0.95 ± 0.40 D - - - 0.05 16 aft 20 Apr 1977 M. Kea 13.0 4.9 2.4 13 0.82 ± 0.17 C - - - 0.014 250 15 Sep 1977 S. flank 9.4 4.2 2.0 259 0.80 ± 0.11 A 0.94 1.19 0.82 0.02 70 12 int 21 Mar 1979 Kohala 13.6 4.7 2.0 23 0.88 ± 0.11 B 1.43 0.66 0.01 0.01 250 13 21 Sep 1979 S. flank 9.3 5.7 1.7 214 0.83 ± 0.06 A 2.28 0.97 0.01 0.01 18 10 Aug 1981 SW flank 2 9 4.3 1.8 229 0.80 ± 0.06 A 1.18 1.01 0.34 0.016 310 14 int 10 Nov 1981 S. flank 10.2 4.5 1.8 16 1.00 ± 0.27 C - - - 0.01 8 21 Jan 1982 Hilea 10.2 5.6 1.8 192 1.06 ± 0.04 A 1.40 0.75 0.03 0.03 660 5 24 Jun 1982 SW flank 7.1 3.0 1.8 417 0.52 ± 0.05 A 0.37 1.53 0.39 0.013 130 17 int 12 Aug 1982 Caldera 16.2 4.3 1.9 9 1.33 ± 0.39 C - - - 0.01 18 20 Mar 1983 S. flank 7.0 4.9 1.8 27 0.82 ± 0.18 B 4.10 1.55 0.02 0.016 18 16 Nov 1983 Kaoiki 11.0 6.7 2.1 947 0.75 ± 0.02 A 3.85 1.26 (0.17) 0.16 355 21 Feb 1985 S. flank 9.4 4.8 1.7 56 0.54 ± 0.11 B 2.64 1.10 0.01 0.01 31 7 Jul 1985 Hilea 10.3 4.5 1.8 13 1.57 ± 0.37 C - - - 0.016 35 3 Feb 1987 Offshore 0.1 5.2 2.7 50 1.02 ± 0.07 B 1.75 1.01 0.02 0.016 630 6 24 Mar 1988 Offshore 2.0 5.0 2.2 12 1.21 ± 0.45 D - - - 0.01 3.1 6,7 27 Mar 1988 Offshore 2.3 5.1 3.0 58 0.93 ± 0.07 B 1.24 0.95 0.02 0.016 550 6 11 May 1988 M. Kea 23.0 4.3 1.9 7 0.92 ±0.24 D- - - - 0.01 140 8 3 Jul 1988 Hilea 9.5 5.4 1.8 54 0.65 ± 0.08 B 3.28 1.02 0.01 0.01 500 25 Jun 1989 S. flank 9.3 6.2 2.0 154 1.09 ± 0.10 A 3.64 1.24 0.06 0.016 15.8 27 Dec 1989 S. flank 9.4 5.3 1.6 53 0.62 ±0.54 D- - - - 1.0 31 9 1 Aug 1990 M. Kea 18.4 4.7 2.1 34 0.80 ± 0.08 B 2.58 1.14 0.01 0.01 700 24 Jan 1993 Caldera 6.0 4.3 1.1 160 2.67 ± 0.49 A 1.17 0.88 0.31 0.01 7.9 26 Jan 1993 Hilea 9.4 5.3 1.4 36 1.35 ± 0.17 B 2.75 1.02 0.03 0.01 40 21 Jan 1996 M. Kea 20.7 4.4 1.6 16 1.27 ± 0.31 C - - - 0.01 8 30 Jun 1997 S. flank 9.2 5.5 1.5 32 0.96 ± 0.19 B 4.22 1.26 0.01 0.01 8 14 Aug 1997 S. flank 8.9 5.0 1.3 16 1.14 ± 0.45 D - - - 0.016 5 27 Sep 1998 Caldera 0.8 4.6 1.8 8 1.44 ± 0.39 C - - - 0.01 80 28 Sep 1998 S. flank 9.5 4.8 1.3 28 0.86 ± 0.27 B 3.17 1.15 0.03 0.01 9 16 Apr 1999 Hilea 9.4 5.6 1.1 163 0.96 ± 0.07 A 1.12 0.57 0.01 0.01 18 26 May 1999 Caldera 7.6 4.3 1.5 26 1.33 ± 0.20 B 1.55 0.74 0.02 0.01 35 Deep (z > 20 km) 26 Apr 1973 Honomu 38.7 6.2 2.4 65 1.11 ± 0.14 A 2.57 1.00 0.29 0.28 178 25 Dec 1974 Kilauea 32.3 4.8 2.0 13 0.97 ± 0.27 C - - - 0.01 16 6 Nov 1975 Kilauea 31.9 4.55 1.7 20 0.94 ± 0.19 B - - - 0.027 70 1 Feb 1994 S Hawaii 34.7 5.2 1.3 93 1.39 ± 0.12 A 3.10 1.08 0.12 0.013 89 17 Feb 2000 Kilauea 35.3 4.5 1.4 15 0.84 ± 0.15 C - - - 0.014 300 a Depth is in km below the local ground surface; M m is main shock magnitude; M min is minimum (completeness) magnitude of aftershock sequence; N is number of aftershocks larger than or equal to M min and between t 1 and t 2 used to fit the four parameters; p, a, b, and c are parameters of aftershock rate expression; and c is usually not well determined (early aftershock recording is incomplete when c values are in parentheses). The letter (A, B, C etc.) is a quality rating of the p value: A, p(error)/p < 0.3 and N > 100; B, p(error)/p < 0.3 and N > 20; C, p(error)/p < 0.3 and N > 8; D, p(error)/p < 0.6 and N > 8 D quality sequences are too poor to plot or interpret (quality letters in lowercase indicate the p and N values meet the same criteria, but are based on the early catalog with more uncertainty); t 1 is start time of aftershocks used for fitting parameters, days after main shock; t 2 is end time of aftershocks used for fitting parameters, days after main shock. 9of26