Geology Reports⌕ Search

Geology topics

Fred W. Klein

Publications and source records attributed to Fred W. Klein.

At least 19 recordsLinked to original sources

Lithospheric flexure under the Hawaiian volcanic load: Internal stresses and a broken plate revealed by earthquakes

Several lines of earthquake evidence indicate that the lithospheric plate is broken under the load of the island of Hawai`i, where the geometry of the lithosphere is circular with a central depression. The plate bends concave downward surrounding a stress-free hole, rather than bending concave upward as with past assumptions. Earthquake focal mechanisms show that the center of load stress and the weak hole is between the summits of Mauna Loa and Mauna Kea where the load is greatest. The earthquake gap at 21 km depth coincides with the predicted neutral plane of flexure where horizontal stress changes sign. Focal mechanism P axes below the neutral plane display a striking radial pattern pointing to the stress center. Earthquakes above the neutral plane in the north part of the island have opposite stress patterns; T axes tend to be radial. The M 6.2 Honomu and M 6.7 Kiholo main shocks (both at 39 km depth) are below the neutral plane and show radial compression, and the M 6.0 Kiholo aftershock above the neutral plane has tangential compression. Earthquakes deeper than 20 km define a donut of seismicity around the stress center where flexural bending is a maximum. The hole is interpreted as the soft center where the lithospheric plate is broken. Kilauea's deep conduit is seismically active because it is in the ring of maximum bending. A simplified two-dimensional stress model for a bending slab with a load at one end yields stress orientations that agree with earthquake stress axes and radial P axes below the neutral plane. A previous inversion of deep Hawaiian focal mechanisms found a circular solution around the stress center that agrees with the model. For horizontal faults, the shear stress within the bending slab matches the slip in the deep Kilauea seismic zone and enhances outward slip of active flanks.

Hawaii↗

Two hundred years of magma transport and storage at Kīlauea Volcano, Hawai'i, 1790-2008

This publication summarizes the evolution of the internal plumbing of Kīlauea Volcano on the Island of Hawaiʻi from the first documented eruption in 1790 to the explosive eruption of March 2008 in Halemaʻumaʻu Crater. For the period before the founding of the Hawaiian Volcano Observatory in 1912, we rely on written observations of eruptive activity, earthquake swarms, and periodic draining of magma from the lava lake present in Kīlauea Caldera. After 1912 the written observations are supplemented by continuous measurement of tilting of the ground at Kīlauea’s summit and by a continuous instrumental record of earthquakes, both measurements made during 1912–56 by a single pendulum seismometer housed on the northeast edge of Kīlauea’s summit. Interpretations become more robust following the installation of seismic and deformation networks in the 1960s. A major advance in the 1990s was the ability to continuously record and telemeter ground deformation to allow its precise correlation with seismic activity before and after eruptions, intrusions, and large earthquakes. We interpret specific events in Kīlauea’s 200- year written history as steps in a broad transition from summit lava-lake activity in Kīlauea Caldera to shield building on the east rift zone. The ability of the magmatic plumbing to deliver magma to eruption is critical to the history of eruption and intrusion. When the rate of magma supply equals the rate of eruption, there is little ground deformation or intrusion. When the magma supply rate is greater than the rate of eruption, then the edifice responds through any or all of summit inflation, intrusion, increased spreading rate, and large flank earthquakes. In Kīlauea’s 200-year history we identify three regions of the volcano in which magma is stored and supplied from below. Source 1 is at 1-km depth or less beneath Kīlauea’s summit and fed Kīlauea’s summit lava lakes throughout most of the 19th century and again from 1907 to 1924. Source 1 was used up in the series of small Halemaʻumaʻu eruptions following the end of lava-lake activity in the summit collapse of 1924. Source 2 is the magma reservoir at a depth of 2–6 km beneath Kīlauea’s summit that has been imaged by seismic and deformation measurements beginning in the 1960s. This source was first identified in the summit collapses of 1922 and 1924. Source 3 is a diffuse volume of magma-permeated rock between 5 and 11 km depth beneath the east rift zone and above the near-horizontal decollement at the base of the Kīlauea edifice. Magma distribution within source 2 has been derived by combining petrologic study of the three chemically uniform summit eruptions of 1952, 1961, and 1967–68 and the east rift eruptions within this interval with both observation of migrating centers of inflation determined from leveling surveys conducted before the 1967–68 eruption and with published models of expected deformation from different source geometries. We adopt a model of concatenated magmatic plugs with nodes beneath the inflation centers. Addition of erupted and intruded volumes of the three summit magma batches yields a liquid magma volume of about 0.2 km3, with dimensions of ~1 km by 1 km by 200 m centered at about 3-km depth within source 2. Following the Halemaʻumaʻu eruption of 1967–68, the chemistry of magma coming into Kīlauea’s summit reservoir has changed frequently, and during the eruption that began in 1983, chemical changes have been subtle and continuous. In this period we interpret changes in chemistry as related to an increase in magma supply resulting from increased partial melting in an expanding mantle source volume. We know from instrumental recording of eruptions since the long Halemaʻumaʻu eruption in 1952 that stress in the edifice accumulates as magma is added underground and is relieved by eruption and by dilation of the rift zones associated with seaward movement (spreading) of Kīlauea’s south flank. During and after the last half of the 20th century, magma transfer to the rift zone has dominantly occurred from source 2. High rates of flank motion have been correlated with high rates of endogenous growth; alternatively, lower rates of motion have characterized periods when the underground magmatic plumbing was being refilled following lateral removal of magma, as well as periods when a more open magmatic plumbing favored continuous eruption. Since at least 1952, source 3 has not drained during deflations, which was apparently not the case before 1924. Triangulation and leveling conducted in 1912, 1921, and 1926, combined with post-1912 tilt measurements, identified a broad regional uplift in 1918–19 and an equally broad collapse in 1924, neither of which has been seen since. We associate these elevation changes with addition or subtraction of magma from all three magma sources, dominantly source 3. We interpret the intrusion beneath the east rift zone during the 1924 collapse to have stabilized the rift zone-south flank relationship, preventing loss of magma from source 3 in subsequent collapses. Rates of seaward spreading were low until 1952, when earthquakes in 1950 and 1951 associated with surges of magma from the hotspot triggered a large offshore south flank earthquake swarm that unlocked the south flank and enabled a greatly increased rate of seaward spreading. Magma supply rates have been derived for the entire period of study. Between 1823 and 1840, magma was supplied from source 1 at a very high rate of more than 0.2 km 3 /yr, which we interpret as recovery from a substantial draining of magma from beneath Kīlauea in 1790. Inferred magma supply rates diminished to one-tenth of that value after 1840, in part because of increase in the activity of Mauna Loa beginning in 1843. Magma supply rates between 1918 and 1924 were about 0.024 km 3 /yr, matching that of the period from 1840 to 1894. During 1950–52 the magma supply rate increased to about 0.06 km 3 /yr, in part because of the great reduction in Mauna Loa activity following its large eruption in June 1950. Following the summit eruption of 1967–68, magma supply increased further to ~0.1 km 3 /yr, and further increases to more than 0.2 km 3 /yr occurred during the east rift eruption that began in 1983. Eruption at Kīlauea’s summit took place in 1952, and eruptive activity steadily increased as increased magma supply also drove increased spreading rates. The inability of magma supply to be accommodated by a combination of eruption and spreading during the 1969–74 Mauna Ulu period stressed Kīlauea’s south flank. The stress was relieved in part by the M7.2 earthquake of 29 November 1975. That earthquake, in turn, dilated Kīlauea’s east rift zone as the south flank moved seaward, producing a favorable condition for continuous east rift eruption, which began in 1983. The 1975 earthquake also resulted in the ability of the south flank to move independently under the influence of gravity, effectively decoupling the spreading rate from changes in the magma supply rate. The continuing increase in magma supply after 1983 was instead manifested in rift dilation, increased intrusion, and ultimately in the launching of a second eruption in Halemaʻumaʻu in March 2008, the first instance in Kīlauea’s recorded history of simultaneous eruption at the summit and on the east rift zone. Kīlauea’s history can be considered in cycles of equilibrium, crisis, and recovery. The approach of a crisis is driven by a magma supply rate that greatly exceeds the capacity of the plumbing to deliver magma to the surface. Crises can be anticipated by inflation measured at Kīlauea’s summit coupled with an increase in overall seismicity, particularly manifest by intrusion and eruption in the southwest sector of the volcano. Unfortunately the nature of the crisis—for example, large earthquake, new eruption, or edifice-changing intrusion—cannot be specified ahead of time. We conclude that Kīlauea’s cycles are controlled by nonlinear dynamics, which underscores the difficulty in predicting eruptions and earthquakes.

Hawai'i↗

New earthquake catalog reexamines Hawaii's seismic history

On April 2,1868, an earthquake of magnitude 7.9 occurred beneath the southern part of the island of Hawaii. The quake, which was felt throughout all of the Hawaiian Islands, had a Modified Mercalli (MM) intensity of XII near its source.The destruction caused by a quake that large is nearly complete. A landslide triggered by the quake buried a small village, killing 31 people, and a tsunami that swept over coastal settlements added to the death toll. We know as much as we do about this and other early earthquakes thanks to detailed records kept by Hawaiian missionaries, including the remarkable diary maintained by the Lyman family that documented every earthquake felt at their home in Hilo between 1833 and 1917 [Wyss et al., 1992].Our analysis of these and other historical records indicates that Hawaii was at least as intensely seismic in the 19th century and first half of the 20th century as in its more recent past, with 26 M ≥6.0 earthquakes occurring from 1823 to 1903 and 20 M ≥6.0 earthquakes from 1904 to 1959. Just five M ≥6.0 earthquakes occurred from 1960 to 1999. The potential damage caused by a repeat of some of the larger historic events could be catastrophic today.

Eos, Transactions, American Geophysical Union↗

Exponential decline of aftershocks of the M7.9 1868 great Kau earthquake, Hawaii, through the 20th century

The remarkable catalog of Hawaiian earthquakes going back to the 1820s is based on missionary diaries, newspaper accounts, and instrumental records and spans the great M 7.9 Kau earthquake of April 1868 and its aftershock sequence. The earthquake record since 1868 defines a smooth curve complete to M 5.2 of the declining rate into the 21st century, after five short volcanic swarms are removed. A single aftershock curve fits the earthquake record, even with numerous M 6 and 7 main shocks and eruptions. The timing of some moderate earthquakes may be controlled by magmatic stresses, but their overall long‐term rate reflects one of aftershocks of the Kau earthquake. The 1868 earthquake is, therefore, the largest and most controlling stress event in the 19th and 20th centuries. We fit both the modified Omori (power law) and stretched exponential (SE) functions to the earthquakes. We found that the modified Omori law is a good fit to the M ≥ 5.2 earthquake rate for the first 10 years or so and the more rapidly declining SE function fits better thereafter, as supported by three statistical tests. The switch to exponential decay suggests that a possible change in aftershock physics may occur from rate and state fault friction, with no change in the stress rate, to viscoelastic stress relaxation. The 61‐year exponential decay constant is at the upper end of the range of geodetic relaxation times seen after other global earthquakes. Modeling deformation in Hawaii is beyond the scope of this paper, but a simple interpretation of the decay suggests an effective viscosity of 10 19 to 10 20 Pa s pertains in the volcanic spreading of Hawaii's flanks. The rapid decline in earthquake rate poses questions for seismic hazard estimates in an area that is cited as one of the most hazardous in the United States.

Hawaii↗

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

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.

Hawai'i↗

User's guide to HYPOINVERSE-2000, a Fortran program to solve for earthquake locations and magnitudes

Hypoinverse is a computer program that processes files of seismic station data for an earthquake (like p wave arrival times and seismogram amplitudes and durations) into earthquake locations and magnitudes. It is one of a long line of similar USGS programs including HYPOLAYR (Eaton, 1969), HYPO71 (Lee and Lahr, 1972), and HYPOELLIPSE (Lahr, 1980). If you are new to Hypoinverse, you may want to start by glancing at the section “SOME SIMPLE COMMAND SEQUENCES” to get a feel of some simpler sessions. This document is essentially an advanced user’s guide, and reading it sequentially will probably plow the reader into more detail than he/she needs. Every user must have a crust model, station list and phase data input files, and glancing at these sections is a good place to begin. The program has many options because it has grown over the years to meet the needs of one the largest seismic networks in the world, but small networks with just a few stations do use the program and can ignore most of the options and commands. History and availability. Hypoinverse was originally written for the Eclipse minicomputer in 1978 (Klein, 1978). A revised version for VAX and Pro-350 computers (Klein, 1985) was later expanded to include multiple crustal models and other capabilities (Klein, 1989). This current report documents the expanded Y2000 version and it supercedes the earlier documents. It serves as a detailed user's guide to the current version running on unix and VAX-alpha computers, and to the version supplied with the Earthworm earthquake digitizing system. Fortran-77 source code (Sun and VAX compatible) and copies of this documentation is available via anonymous ftp from computers in Menlo Park. At present, the computer is swave.wr.usgs.gov and the directory is /ftp/pub/outgoing/klein/hyp2000. If you are running Hypoinverse on one of the Menlo Park EHZ or NCSN unix computers, the executable currently is ~klein/hyp2000/hyp2000. New features. The Y2000 version of Hypoinverse includes all of the previous capabilities, but adds Y2000 formats to those defined earlier. In most cases, the new formats add 2 digits to the year field to accommodate the century. Other fields are sometimes rearranged or expanded to accommodate a better field order. The Y2000 formats are invoked with the “200” command. When the Y2000 flag is turned on, all files are read and written in the new format and there is no mixing of format types in a single run. Some formats without a date field, like station files, have not changed. A separate program called 2000CONV has been written to convert old formats to new. Other new features, like expanded station names, calculating amplitude magnitudes from a variety of digital seismometers, station history files, interactive earthquake processing, and locations from CUSP (Caltech USGS Seismic Processing) binary files have been added. General features. Hypoinverse will locate any number of events in an input file, which can be in one of several different formats. Any or all of printout, summary or archive output may be produced. Hypoinverse is driven by user commands. The various commands define input and output files, set adjustable parameters, and solve for locations of a file of earthquake data using the parameters and files currently set. It is both interactive and "batch" in that commands may be executed either from the keyboard or from a file. You execute the commands in a file by typing @filename at the Hypoinverse prompt. Users may either supply parameters on the command line, or omit them and are prompted interactively. The current parameter values are displayed and may be taken as defaults by pressing just the RETURN key after the prompt. This makes the program very easy to use, providing you can remember the names of the commands. Combining commands with and without their required parameters into a command file permits a variety of customized procedures such as automatic input of crustal model and station data, but prompting for a different phase file each time. All commands are 3 letters long and most require one or more parameters or file names. If they appear on a line with a command, character strings such as filenames must be enclosed in apostrophes (single quotes). Appendix 1 gives this and other free-format rules for supplying parameters, which are parsed in Fortran. When several parameters are required following a command, any of them may be omitted by replacing them with null fields (see appendix 1). A null field leaves that parameter unchanged from its current or default value. When you start HYPOINVERSE, default values are in effect for all parameters except file names. Hypoinverse is a complicated program with many features and options. Many of these "advanced" or seldom used features are documented here, but are more detailed than a typical user needs to read about when first starting with the program. I have put some of this material in smaller type so that a first time user can concentrate on the more important information.

Open-File Report↗

Seismic-hazard maps for Hawaii

Probabilistic seismic hazard maps were prepared for Hawaii portraying peak horizontal ground acceleration and horizontal spectral response acceleration for 0.2, 0.3, and 1.0 second periods with probabilities of exceedance of 10% in 50 years and 2% in 50 years. This particular data set is for horizontal spectral response acceleration for 1.0 second period with a 10% probability of exceedance in 50 years. All of the maps were prepared by combining hazard derived from spatially- smoothed historic seismicity with hazard from fault-specific sources. The acceleration values contoured are the random horizontal component. The reference site condition is firm rock, defined as having an average shear-wave velocity of 760 m/sec in the top 30 meters corresponding to the boundary between NEHRP site classes B and C. This data set represents the results of calculations of hazard curves for a grid of points with a spacing of 0.02 degrees in latitude and longitude. The grid of points were contoured to produce the final representation of the seismic-hazard.

Hawaii↗

Catalog of Hawaiian earthquakes, 1823-1959

This catalog of more than 17,000 Hawaiian earthquakes (of magnitude greater than or equal to 5), principally located on the Island of Hawaii, from 1823 through the third quarter of 1959 is designed to expand our ability to evaluate seismic hazard in Hawaii, as well as our knowledge of Hawaiian seismic rhythms as they relate to eruption cycles at Kilauea and Mauna Loa volcanoes and to subcrustal earthquake patterns related to the tectonic evolution of the Hawaiian chain.

Professional Paper↗

Magma migration and resupply during the 1974 summit eruptions of Kilauea Volcano, Hawaii

The purpose of this paper is to present a complete account of contrasting yet related eruptions, thus filling a gap in the published narratives of recent activity of Kilauea; and to examine their significance within a broader context of regional magmatic and eruptive dynamics. We have gained a historical perspective and can view these three eruptions within a multidecade context of the eruptive behavior of not only Kilauea, but also of the adjacent Mauna Loa.

Professional Paper↗

Relaxation of the south flank after the 7.2-magnitude Kalapana earthquake, Kilauea Volcano, Hawaii

An M = 7.2 earthquake on 29 November 1975 caused the south flank of Kilauea Volcano, Hawaii, to move seaward several meters: a catastrophic release of compression of the south flank caused by earlier injections of magma into the adjacent segment of a rift zone. The focal mechanisms of the mainshock, the largest foreshock, and the largest aftershock suggest seaward movement of the upper block. The rate of aftershocks decreased in a familiar hyperbolic decay, reaching the pre-1975 rate of seismicity by the mid-1980s. Repeated rift-zone intrusions and eruptions after 1975, which occurred within 25 km of the summit area, compressed the adjacent portion of the south flank, apparently masking continued seaward displacement of the south flank. This is evident along a trilateration line that continued to extend, suggesting seaward displacement, immediately after the M = 7.2 earthquake, but then was compressed during a series of intrusions and eruptions that began in September 1977. Farther to the east, trilateration measurements show that the portion of the south flank above the aftershock zone, but beyond the area of compression caused by the rift-zone intrusions and eruptions, continued to move seaward at a decreasing rate until the mid-1980s, mimicking the decay in aftershock rate. Along the same portion of the south flank, the pattern of vertical surface displacements can be explained by continued seaward movement of the south flank and development of two eruptive fissures along the east rift zone, each of which extended from a depth of ∼3 km to the surface. The aftershock rate and continued seaward movement of the south flank are reminiscent of crustal response to other large earthquakes, such as the 1966 M = 6 Parkfield earthquake and the 1983 M = 6.5 Coalinga earthquake.

Hawaii↗