Search USGS⌕ Search

USGS · pp1806

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

Abstract

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.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 18.6° to 20° latitude; -155.9° to -154° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Thomas L. Wright, Fred W. Klein. 2014. Two hundred years of magma transport and storage at Kīlauea Volcano, Hawai'i, 1790-2008. https://doi.org/10.3133/pp1806

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

The effects of management practices on grassland birds—Vesper Sparrow ( Pooecetes gramineus )

Keys to Vesper Sparrow ( Pooecetes gramineus ) management include providing dry, open areas with short, sparse vegetation or edge habitat between forests and fields and providing elevated song perches. Vesper Sparrows have been reported to use habitats with 8–134 centimeters (cm) average vegetation height, 2–21 cm visual obstruction reading, 3–52 percent grass cover, no more than 76 percent forb cover, no more than 17 percent shrub cover, 1–57 percent bare ground, 6–63 percent litter cover, and no more than 5 cm litter depth.

Professional Paper↗

Allostratigraphy and formalizing the Cold Harbor Alloformation, Virginia Coastal Plain, Eastern Virginia

Introduction Distinct morphologic features in the Virginia part of the Coastal Plain Province (herein referred to as the Virginia Coastal Plain) were noted in one of the earliest geologic surveys in 1839 (Rogers, 1884). South of the James River, Rogers (1884, p. 252) described a western bench that ranged from “a low rocky coast…of the primeval sea” to an eastern facing scarp east of the City of Suffolk, Virginia. His eastern bench extended from the scarp at Suffolk, which he suggested was a former shoreline, to the present Atlantic Ocean. Notably, the eastern plain was described as less dissected than the western plain. In Virginia, creation and publication of topographic maps by the U. S. Geological Survey began in the mid-1880s. By 1905, many areas in eastern Virginia were mapped as 30-minute and 15-minute series with 10-foot (ft) contour interval accuracy (Virginia Division of Mineral Resources [VDMR] , 2008). After the continued release of these early maps that displayed Virginia Coastal Plain morphology in increasingly greater detail, several geological reports were published where scarps and plains were traced, named, and origins proposed (Shaler, 1890; Shattuck, 1901; Clark and Miller, 1906, 1912; Wentworth, 1930). Shattuck (1901) recognized three continuous scarps in Maryland and most importantly suggested that there was a correlation between surficial deposits (geologic formations) and morphologic terraces; this defined the terrace-formation concept that became the basis of geologic mapping of the Upper Pliocene and Pleistocene sediments in the Virginia Coastal Plain. Clark and Miller (1912) described three Pleistocene formations and supported the terrace-formation concept because the sediments within each formation were lithologically similar to one another, and differentiating them by comparing key vertical sections was not helpful. In Virginia, Wentworth (1930) defined and mapped five Pleistocene terraces on the basis of correlating morphologic terraces with their underlying formation sediments. Similarly, Cooke (1931) identified and traced seven marine shorelines and proposed corresponding terrace names from Maryland to Florida. Cooke (1931) also pointed out that the various elevations of scarps were related to interglacial sea-level highstands whereas prior workers called on vertical land movement of emergence and submergence to create shoreline scarps. Cooke’s (1936) geologic map of eastern South Carolina showing the Pleistocene sediments is also based on the terrace-formation concept. Detailed geologic mapping of the Virginia Coastal Plain by Virginia Division of Mineral Resources (VDMR) and U.S. Geological Survey (USGS) followed lithostratigraphic convention for Early Pliocene and older marine sediments. However, geologic mapping of Upper Pliocene and Pleistocene sediments used the concept of morphologic setting and elevations of scarp toe and flat-plain, or the terrace-formation concept. The 1983 version of the North American Stratigraphic Code (hereafter the Code) (North American Commission on Stratigraphic Nomenclature [NACSN] , 1983) introduced allostratigraphy that gave credence to the terrace-formation concept. According to the current Code (NACSN, 2021; Article 58), “an allostratigraphic unit is a mappable body of rock that is defined and identified on the basis of its bounding discontinuities,” and, “formal allostratigraphic units may be defined or distinguish between different (1) superposed discontinuity-bounded deposits of similar lithology, (2) contiguous discontinuity-bounded deposits of similar lithology, or (3) geographically separated discontinuity-bounded units of similar lithology…The lithology of an allostratigraphic unit plays no part in its definition.” However, definition and placement of bounding discontinuities relies on lateral and vertical knowledge of lithology and is so recognized in Article 58e, “A type locality and type area must be designated; a composite stratotype or a type section and several reference sections are desirable.” We understand that each of the terrace formations of Late Pliocene and Pleistocene age were deposited during a marine transgressive event tied to an interglacial sea-level highstand (scarp above the Poquoson Member of the Tabb Formation may be the result of a still-stand as sea level was falling after deposition of the Lynnhaven Member, refer to fig. 1). The base of each map unit is an unconformity (or is at least in part for the base of the Poquoson) and nomenclature for each unit might be changed from “formation” to “alloformation,” as suggested by Doar and Berquist (2009). Figure 1 shows the relative placement of these mapped and named alloformations as used in the Virginia Coastal Plain and Fall Zone (a zone of transition between the Piedmont and Coastal Plain Provinces) (Berquist, 2021). Except perhaps for the Poquoson, many of these proposed alloformations are unlikely to interfinger or grade laterally seaward into extensive deep marine sediments. After their deposition, erosion during sea-level fall and then rise for the next depositional event would seemingly obliterate any marine connection without a basal unconformity. All these alloformations contain sediments deposited in fluvial, estuarine, and shallow-marine environments and are bound everywhere by unconformities. This paper documents the Cold Harbor Alloformation as a formal stratigraphic unit and the rationale for mapping the sediments using an allostratigraphic classification. The Cold Harbor Alloformation is composed of sediments from a variety of environments such as gravel beds, shoreline beach sands, marsh and estuarine sand and mud, and shallow marine sands. These sediments are found at elevations up to 240 ft and lie unconformably above older sediments and rocks. Our work relied on observations of outcrops and many auger borings to identify basal unconformities.

Virginia↗

Formalization of the Elsing Green Alloformation, Virginia Coastal Plain, Eastern Virginia

The Elsing Green Alloformation is herein formalized as a new Pleistocene alloformation in the Virginia Coastal Plain. This allostratigraphic unit is characterized by an unconformity-bounded terrace deposit occupying a morphologic flat with a mapped elevation range of about 26 to 38 feet (7.9 to 11.6 meters); its landward extent is commonly marked at the base of an erosional scarp where the toe is typically at 38 feet (11.6 meters) and is distinguished by its topographic expression and stratigraphic position. Sediments of this unit were previously mapped as part of the Shirley Formation or Sedgefield Member of the Tabb Formation but differ significantly in their bounding scarps and terrace morphology. The type area is within the Tunstall 7.5-minute quadrangle near Elsing Green Plantation, King William County, Virginia. The reference section is provided by boring 2012270-2. Correlation is proposed with other mid-Pleistocene terrace deposits in the mid-Atlantic Coastal Plain, including the younger Lynch Heights Formation of Delaware and the Butlers Bluff Member of the Nassawadox Formation of the Eastern Shore of Virginia. The age of the Elsing Green Alloformation is interpreted to be Marine Isotope Stage 9e (about 330,000 years).

Virginia↗