Storage, migration, and eruption of magma at Kilauea volcano, Hawaii, 1971–1972

The magmatic plumbing system of Kilauea Volcano consists of a broad region of magma generation in the upper mantle, a steeply inclined zone through which magma rises to an intravolcano reservoir located about 2 to 6 km beneath the summit of the volcano, and a network of conduits that carry magma from this reservoir to sites of eruption within the caldera and along east and southwest rift zones. The functioning of most parts of this system was illustrated by activity during 1971 and 1972. When a 29-month-long eruption at Mauna Ulu on the east rift zone began to wane in 1971, the summit region of the volcano began to inflate rapidly; apparently, blockage of the feeder conduit to Mauna Ulu diverted a continuing supply of mantle-derived magma to prolonged storage in the summit reservoir. Rapid inflation of the summit area persisted at a nearly constant rate from June 1971 to February 1972, when a conduit to Mauna Ulu was reopened. The cadence of inflation was twice interrupted briefly, first by a 10-hour eruption in Kilauea Caldera on 14 August, and later by an eruption that began in the caldera and migrated 12 km down the southwest rift zone between 24 and 29 September. The 14 August and 24–29 September eruptions added about 10⁷ m³ and 8 × 10⁶ m³, respectively, of new lava to the surface of Kilauea. These volumes, combined with the volume increase represented by inflation of the volcanic edifice itself, account for an approximately 6 × 10⁶ m³/month rate of growth between June 1971 and January 1972, essentially the same rate at which mantle-derived magma was supplied to Kilauea between 1952 and the end of the Mauna Ulu eruption in 1971.The August and September 1971 lavas are tholeiitic basalts of similar major-element chemical composition. The compositions can be reproduced by mixing various proportions of chemically distinct variants of lava that erupted during the preceding activity at Mauna Ulu. Thus, part of the magma rising from the mantle to feed the Mauna Ulu eruption may have been stored within the summit reservoir from 4 to 20 months before it was erupted in the summit caldera and along the southwest rift zone in August and September.The September 1971 activity was only the fourth eruption on the southwest rift zone during Kilauea's 200 years of recorded history, in contrast to more than 20 eruptions on the east rift zone. Order-of-magnitude differences in topographic and geophysical expression indicate greatly disparate eruption rates for far more than historic time and thus suggest a considerably larger dike swarm within the east rift zone than within the southwest rift zone. Characteristics of the historic eruptions on the southwest rift zone suggest that magma may be fed directly from active lava lakes in Kilauea Caldera or from shallow cupolas at the top of the summit magma reservoir, through fissures that propagate down rift from the caldera itself at the onset of eruption. Moreover, emplacement of this magma into the southwest rift zone may be possible only when compressive stress across the rift is reduced by some unknown critical amount owing either to seaward displacement of the terrane south-southeast of the rift zone or to a deflated condition of Mauna Loa Volcano adjacent to the northwest, or both. The former condition arises when the forceful emplacement of dikes into the east rift zone wedges the south flank of Kilauea seaward. Such controls on the potential for eruption along the southwest rift zone may be related to the topographic and geophysical constrasts between the two rift zones.     

Deformation measurements on Kilauea volcano,Hawaii

Repeated electronic distance measurements across Kilauea Caldera with Tellurometers and Geodimeter show definite horizontal expansion related to the vertical uplift and outward tilting of the summit prior to an eruption, and contraction during and after a flank eruption. Measurements started in October 1964, along a 3098 meter line between Uwekahuna and Keanakakoi, indicate a relatively uniform lengthening of 12 centimeters during the interval October 22, 1964 to March 1, 1965. Rapid shortening of the line by 28 centimeters was measured 4 days after the beginning of a flank eruption which involved emission of approximately 29 million cubic meters of lava during the period March 5 to March 15, 1965. During the expansion, the standard deviation of 10 Tellurometer measurements from a least-squares srtaight line solution is ± 2.0 centimeters (6.5 ppm) whereas 9 Geodimeter measurements have a standard deviation of ± 1.1 (3.6 ppm) centimeters. Absolute distance readings between the two instruments differ by 4 centimeters (13 ppm), but relative changes in distance were the same on both instruments. Changes in distance across Kilauea Caldera can, therefore, be easily measured to accuracies of 4 to 7 parts per million with standard electronic distance measuring systems. On active volcanoes where ground surface deformation exceeds 10–100 parts per million with changes in subsurface magma pressure or volume, repeated horizontal distance measurements can be a most useful technique.     

Some results from recent chemical studies at Kilauea volcano,Hawaii

The chemical surveillance of Kilauea volcano, Hawaii, has continued. No relationship has thus far been identified between the helium content of an associated fumarole and the activity at the volcano. Fume samples from Halemaumau crater in Kilauea caldera and from a fissure eruption that occurred nearby on the floor of the caldera during August 1971 were examined for their halogen (Cl and F) and sulfur contents. The ratio of Cl/F in fume showed an abnormal increase in samples taken at Halemaumau a month before the eruption. This change in ratio may be a helpful indicator of the onset of eruption in volcanic areas.     

Geometry of the September 1971 eruptive fissure at Kilauea volcano,Hawaii

A three-dimensional model has been used to estimate the location and dimensions of the eruptive fissure for the 24–29 September 1971 eruption along the southwest rift zone of Kilauea volcano, Hawaii. The model is an inclined rectangular sheet embedded in an elastic half-space with constant displacement on the plane of the sheet. The set of best model parameters suggests that the sheet is vertical, extends from a depth of about 2 km to the surface, and has a length of about 14 km. Because this sheet intersects the surface where eruptive vents and extensive ground cracking formed during the eruption, this sheet probably represents the conduit for erupted lava. The amount of displacement perpendicular to the sheet is about 1.9 m, in the middle range of values measured for the amount of opening across the September 1971 eruptive fissure. The thickness of the eruptive fissure associated with the January 1983 east rift zone eruption was determined in an earlier paper to be 3.6 m, about twice the thickness determined here for the September 1971 eruption. Because the lengths (12 km for 1983 and 14 km for 1971) and heights (about 2 km) of the sheet models derived for the January 1983 and September 1971 rift zone eruptions are nearly identical, the greater thickness for the January 1983 eruptive fissure implies that the magma pressure was about a factor of two greater to form the January 1983 eruptive fissure. Because the September 1971 and January 1983 eruptive fissures extent to depths of only a few kilometers, the region of greatest compressive stress produced along the volcano's flank by either of these eruptive fissures would also be within a few kilometers of the surface. Previous work has shown that rift eruptions and intrusions contribute to the buildup of compressive stress along Kilauea's south flank and that this buildup is released by increased seismicity along the south flank. Because south flank earthquakes occur at significantly greater depths, i.e., from 5 to 13 km, than the vertical extent of the 1971 and 1983 eruptiv fissures, the depth of emplacement of these eruptive fissures cannot be the main factor in controlling the hypocentral depths of south flank earthquakes. Two possible explanations for the occurrence of south flank earthquakes in the depth range of 5–13 km are (1) a deeper pressure source, possibly related to deeper magma storage within the rift zone, and (2) a lowstrength region located between 5 and 13 km beneath Kilauea's south flank, possibly at the interface between oceanic sediments and the base of the Hawaiian volcanics.     

Pit crater formation on Kilauea volcano, Hawaii.

Most of the known pit craters in Hawaii occur along the East and Southwest Rift Zones of Kilauea volcano. The pit craters typically are either astride a single rift zone fracture or between a pair of rift zone fractures. These fractures are prominent in the pit crater walls. The pit craters are elliptical in plan view, with their major diameters ranging from 8 to 1140 m. They range in depth from 6 m to 186 m. They typically develop with initially steep, locally overhanging walls, but as the walls collapse, the craters fill with talus and become shaped like inverted elliptical cones. None of the craters apparently formed as eruptive vents, although some have been subsequently filled by lava. Devil's Throat is the best-exposed pit crater along the East Rift Zone. It is sited at a `waist' between two east-striking zones of ground cracks; the spacing between the crack zones decreases towards Devil's Throat. East-striking fractures are also prominent in the pit crater walls. Pit craters along the Southwest Rift Zone typically are elongate in plan view along the direction of the rift, have large caves at their bases along the long axes of the craters, and are smaller than those of the East Rift Zone. Some closely spaced pits there have coalesced to form a trough. Based on our observations and mechanical considerations, we infer that pit craters form by stoping over an underlying large-aperture rift zone fracture, and not by piston-like collapse over broad magma bodies or voids. Flow of magma along the underlying fracture may remove stoped blocks and prevent the fracture from being choked with debris. This mechanism is consistent with pit crater location, ground crack patterns, the preferred orientation of fractures in pit crater walls, and pit crater geometry (both in map view and cross-section). The mechanism also fits with observations of stoping into a gaping rift fracture that conducted lava from Kilauea caldera during the 1920s. Additionally, the ratio of pit crater width to depth of 0.5 to 2 is consistent with pit craters forming over a nearly vertical opening mode fracture.     

Petrology of the hilina formation,Kilauea volcano,Hawaii

The Hilina Formation comprises the oldest sequence of lava flows and tuffs exposed on Kilauea Volcano. These rocks are only exposed in kipukas in younger Puna Formation lavas along cliffs on the south flank of Kilauea Volcano. Locally, tuffs and flows of the Pahala Formation separate the underlying Hilina Formation rocks rom the overlying Puna Formation rocks. Charcoal collected from the base of the Pahala Formation yielded a C¹⁴ age of 22.800±340 years B.P. which defines a minimum age for the Hilina Formation. Hilina Formation lavas crop out over a wide region and probably originated from the summit area and from both rift zones. The Hilina Formation contains both olivine-controlled and differentiated lavas (using the terminology ofWright, 1971). The olivine-controlled lavas of the Hilina Formation are distinguishable mineralogically and geochemically from younger olivine-controlled Kilauea lavas. The younger lavas generally contain discrete low-calcium pyroxene grains. greater glass contents, higher K₂O/P₂O₅ ratios and lower total iron contents. Similar geochemical trends prevail for Manuna Loa lavas, and may typify the early lavas of Hawaiian shield volcanoes. Despite these similarities, the Hilina Formation (and all Kilauea) lavas have higher TiO₂ and CaO, and lower SiO₂ and Al₂O₃ contents than Mauna Loa Lavas. These differences have existed for over 30,000 years. Therefore, it is unlikely that the older lavas of Kilauea are compositionally similar to recent Mauna Loa lavas as was previously suggested. K₂O, TiO₂, Na₂ and Zr contents of lavas from a stratigraphic sequence of Hilina Formation lavas are variable. These variations may be utilized to subdivide the sequence into geochemical groups. These groups are not magma batches. Rather, they represent lavas from batches whose compositions may have been modified by crystal fractionation and magma mixing.     

Trends in the aggregated rate of pre-eruptive volcano-tectonic seismicity at Kilauea volcano, Hawaii

Accelerating rates of volcano-tectonic (VT) earthquakes are commonly observed during volcanic unrest. Understanding the repeatability of their behaviour is essential to evaluating their potential to forecast eruptions. Quantitative eruption forecasts have focused on changes in precursors over intervals of weeks or less. Previous studies at basaltic volcanoes in frequent eruption, such as Kilauea in Hawaii and Piton de La Fournaise on Réunion, suggest that VT earthquake rates tend to follow a power-law acceleration with time about 2 weeks before eruption, but that this trend is often obscured by random fluctuations (or noise) in VT earthquake rate. These previous studies used a stacking procedure, in which precursory sequences for several eruptions are combined to enhance the signal from an underlying acceleration in VT earthquake rate. Such analyses assume a common precursory trend for all eruptions. This assumption is tested here for the 57 eruptions and intrusions recorded at Kilauea between 1959 and 1984. Applying rigorous criteria for selecting data (e.g. maximum depth; restricting magnitudes to be greater than the completeness magnitude, 2.1), we find a much less pronounced increase in the aggregate rate of earthquakes than previously reported. The stacked trend is also strongly controlled by the behaviour of one particular pre-eruptive sequence. In contrast, a robust signal emerges among stacked VT earthquake rates for a subset of the eruptions and intrusions. The results are consistent with two different precursory styles at Kilauea: (1) a small proportion of eruptions and intrusions that are preceded by accelerating rates of VT earthquakes over intervals of weeks to months and (2) a much larger number of eruptions that show no consistent increase until a few hours beforehand. The results also confirm the importance of testing precursory trends against data that have been filtered according to simple constraints on the spatial distribution and completeness magnitude of the VT earthquakes.     

Dynamics of a confined lava flow on Kilauea volcano,Hawaii

This paper presents a new method of analysing lava flow deposits which allows the velocity, discharge rate and rheological properties of channelled moving lavas to be calculated. The theory is applied to a lava flow which was erupted on Kilauea in July 1974. This flow came from a line of fissures on the edge of the caldera and was confined to a pre-existing gully within 50 m of leaving the vent. The lava drained onto the floor of the caldera when the activity stopped, but left wall and floor deposits which showed that the lava banked up as it flowed around each of the bends. Field surveys established the radius of curvature of each bend and the associated lava levels, and these data, together with related field and laboratory measurements, are used to study the rheology of the lava. The results show the flow to have been fast moving but still laminar, with a mean velocity of just over 8 m s^–1; the lava had a low or negligible yield strength and viscosities in the range 85–140 Pa s. An extension of the basic method is considered, and the possibility of supercritical flow discussed.     

Morphology and composition of spinel in Pu’u ’O’o lava (1996–1998), Kilauea volcano, Hawaii

The morphology and composition of spinel in rapidly quenched Pu’u ’O’o vent and lava tube samples are described. These samples contain glass, olivine phenocrysts (3–5 vol.%) and microphenocrysts of spinel (0.05 vol.%). The spinel surrounded by glass occurs as idiomorphic octahedra 5–50 μm in diameter and as chains of octahedra that are oriented with respect to each other. Spinel enclosed by olivine phenocrysts is sometimes rounded and does not generally form chains. The temperature before quenching was calculated from the MgO content of the glass and ranges from 1150°C to 1180°C. The oxygen fugacity before quenching was calculated by two independent methods and the log f_O2 ranged from −9.2 to −9.9 (delta QFM=−1). The spinel in the Pu’u ’O’o samples has a narrow range in composition with Cr/(Cr+Al)=0.61 to 0.73 and Fe²⁺/(Fe²⁺+Mg)=0.46 to 0.56. The lower the calculated temperature for the samples, the higher the average Fe²⁺/(Fe²⁺+Mg), Fe³⁺ and Ti in the spinel. Most zoned spinel crystals decrease in Cr/(Cr+Al) from core to rim and, in the chains, the Cr/(Cr+Al) is greater in the core of larger crystals than in the core of smaller crystals. The occurrence of chains and hopper crystals and the presence of Cr/(Cr+Al) zoning from core to rim of the spinel suggest diffusion-controlled growth of the crystals. Some of the spinel crystals may have grown rapidly under the turbulent conditions of the summit reservoir and in the flowing lava, and the crystals may have remained in suspension for a considerable period. The rapid growth may have caused very local (μm) gradients of Cr in the melt ahead of the spinel crystal faces. The crystals seem to have retained the Cr/(Cr+Al) ratio that developed during the original growth of the crystal, but the Fe²⁺/(Fe²⁺+Mg) ratio may have equilibrated fairly rapidly with the changing melt composition due to olivine crystallization. Six of the samples were collected on the same day at various locations along a 10-km lava tube and the calculated pre-collection temperatures of the samples show a 5°C drop with distance from the vent. The average Fe²⁺/(Fe²⁺+Mg) of the spinel in these samples shows a weak positive correlation with decreasing MgO in the glass of these samples. The range in Cr₂O₃ (0.041–0.045 wt.%) of the glass for these six samples is too small to distinguish a consistent change along the lava tube. The spinel in the Pu’u ’O’o samples shows a zoning trend in a Cr–Al–Fe³⁺ diagram almost directly away from the Cr apex. This compares with a zoning trend in rapidly quenched MORB samples away from Cr coupled with decreasing Fe³⁺. The trend away from Cr displayed by spinel in rapidly quenched samples is in marked contrast to the trend of increasing Fe³⁺ shown by spinel in slowly cooled lava.     

The relationship between earthquake swarms and magma transport: Kilauea volcano,Hawaii

A detailed investigation of earthquake locations and focal mechanisms for swarms associated with intrusive events at Kilauea volcano, Hawaii, further illuminates the relationships among stress state, faulting, and magma transport. We determine the earthquake locations and mechanisms using a three-dimensional crustal model to improve their accuracy and consistency. Swarms in Kilauea's upper east and southwest rift zones, from the years 1980 through 1982, provide clear evidence for the propagation and/or dilation of dikes. Focal mechanisms are predominantly strike-slip, and the faulting and inferred dike orientations can be interpreted quite consistently in terms of the model ofHill (1977). Stresses induced by the summit magma reservoir system strongly control faulting and magma transport in the rift zones close to the summit.     

Gas analyses from the Pu'u O'o eruption in 1985, Kilauea volcano,Hawaii

Volcanic gas samples were collected from July to November 1985 from a lava pond in the main eruptive conduit of Pu'u O'o from a 2-week-long fissure eruption and from a minor flank eruption of Pu'u O'o. The molecular composition of these gases is consistent with thermodynamic equilibrium at a temperature slightly less than measured lava temperatures. Comparison of these samples with previous gas samples shows that the composition of volatiles in the magma has remained constant over the 3-year course of this episodic east rift eruption of Kilauea volcano. The uniformly carbon depleted nature of these gases is consistent with previous suggestions that all east rift eruptive magmas degas during prior storage in the shallow summit reservoir of Kilauea. Minor compositional variations within these gas collections are attributed to the kinetics of the magma degassing process.     

Volcanic geology and eruption frequency,lower east rift zone of Kilauea volcano,Hawaii

Detailed geologic mapping and radiocarbon dating of tholeiitic basalts covering about 275 km² on the lower east rift zone (LERZ) and adjoining flanks of Kilauea volcano, Hawaii, show that at least 112 separate eruptions have occurred during the past 2360 years. Eruptive products include spatter ramparts and cones, a shield, two extensive lithic-rich tuff deposits, aa and pahoehoe flows, and three littoral cones. Areal coverage, number of eruptions and average dormant interval estimates in years for the five age groups assigned are: (I) historic, i.e. A D 1790 and younger: 25%, 5, 42.75; (II) 200–400 years old: 50%, 15, 14.3: (III) 400–750 years old: 20%, 54, 6.6; (IV) 750–1500 years old: 5%, 37, 20.8; (V) 1500–3000 years old: <1%, 1, unknown. At least 4.5–6 km³ of tholeiitic basalt have been erupted from the LERZ during the past 1500 years. Estimated volumes of the exposed products of individual eruptions range from a few tens of cubic meters for older units in small kipukas to as much as 0.4 km³ for the heiheiahulu shield. The average dormant interval has been about 13.6 years during the past 1500 years. The most recent eruption occurred in 1961, and the area may be overdue for its next eruption. However, eruptive activity will not resume on the LERZ until either the dike feeding the current eruption on the middle east rift zone extends farther down rift, or a new dike, unrelated to the current eruption, extends into the LERZ.     

Precursors to dyke-fed eruptions at basaltic volcanoes: insights from patterns of volcano-tectonic seismicity at Kilauea volcano,Hawaii

To investigate the physical controls on volcano-tectonic (VT) precursors to eruptions and intrusions at basaltic volcanoes, we have analyzed the spatial and temporal patterns of VT earthquakes associated with 34 eruptions and 23 dyke intrusions that occurred between 1960 and 1983 at Kilauea, in Hawaii. Eighteen of the 57 magmatic events were preceded by an acceleration of the mean rate of VT earthquakes located close to the main shallow magma reservoir. Using a maximum-likelihood technique and the Bayesian Information Criterion for model preference, we demonstrate that an exponential acceleration is preferred over a power-law acceleration for all sequences. These sequences evolve over time-scales of weeks to months and are consistent with theoretical models for the approach to volcanic eruptions based on the growth of a population of fractures in response to an excess magma pressure. Among the remaining 40 magmatic events, we found a significant correlation between swarms of VT earthquakes located in the mobile south-flank of Kilauea and eruptions and intrusions. The behaviour of these swarms suggests that at least some of the magmatic events are triggered by transient episodes of elevated rates of aseismic flank movement, which could explain why many eruptions and intrusions are not preceded by longer-term precursory signals. In none of the 57 cases could a precursory sequence be used to distinguish between the approach to an eruption or an intrusion, so that, even when a precursory sequence is recognized, there remains an empirical chance of about 40% (24 intrusions from 57 magmatic events) of issuing a false alarm for an imminent eruption.     

Crystal size distribution systematics and the determination of magma storage times: The 1959 eruption of Kilauea volcano, Hawaii

The theory of crystal size distribution (CSD theory) is based on a steady-state population balance that monitors the flux of crystals growing into and out of specificed size categories in precipitating solutions. The conservation equation describing this balance permits crystal growth and nucleation rates to be determined directly from crystal size distribution data. In this investigation, it is shown that CSD analysis can alternatively be used to calculate the residence time of crystals in the system - or in volcanologic terms, the magma storage time prior to eruption - if the characteristic crystal growth rate can be independently determined or estimated. The crystal size data needed for storage time determinations are easily obtained from thin sections of glassy eruption samples and the mathematical calculations are relatively simple. Analysis of the errors inherent in this new technique predicts storage estimates accurate to within an order of magnitude.     

Summit outgassing as indicated by radon, mercury and pH mapping, Kilauea volcano, Hawaii

Surface geochemical mapping of Rn, Hg and pH on the summit of Kilauea volcano, Hawaii, has shown some of the characteristics of outgassing. Rn concentrations measured in shallow ground gas are highest (10 to 155 units, i.e. up to 0.016 nCi/l²²²Rn) over deep structures associated with the summit caldera, low (1 to 5 units) over the upper rift zones (which are underlain by shallow intrusions) and low to negative (1 to -3 units) over most of the volcanically inactive peripheral areas. High Rn concentrations were measured over peripheral well-developed structures such as the Kaoiki and Koae fault systems (45 and 12 units, respectively). The pattern of Rn values broadly suggests the existence of a summit-wide convection system interrupted locally by specific permeable structures. The low pH (3.5 to 5.6) of soils over the caldera is suggested to be largely due to the sulfurous component of the ground gas. The low concentrations of Hg (5 to 80 ppb) in these soils are postulated as being partly a function of subsurface complexing of Hg with the abundant sulfur. The shallow intrusions below the upper rift zones appear to have outgassed much of their sulfur content, so soils in these areas are only slightly acid (5.6 to 6.2), allowing the accumulation of Hg to concentrations of several thousands of ppb. These conditions change during extrusive (and intrusive) events with the increases in shallow subsurface temperatures and volatiles in ground gas producing an increase in Rn concentration, and a decrease in soil Hg due to volatilization. After such events, Rn concentration decreases, and soil Hg increases significantly.     

Rheology of the 1983 Royal Gardens basalt flows,Kilauea Volcano,Hawaii

Ten carefully surveyed topographic profiles across a 1983 Royal Gardens basalt flow from the East Rift of the Kilauea Volcano were matched to digitally derived preflow profiles to construct accurate flow cross sections. Geometric parameters measured on these sections were then used to compute yield strengths and viscosities by means of several rheologic models. Calculated yield strengths (1.5–50 × 10³ Pa) and viscosities (0.2–8.2 × 10⁶ Pas) are comparable to earlier field estimates and slightly higher than laboratory determined values for aa basalt. Both yield strength and viscosity increased systematically downstream. The maximum observed temperature drop of 30 °C is insufficient to account for the 30-fold increase in yield strength, but could explain the three-fold order-of-magnitude increase in viscosity. The yield-strength increase downstream is more likely due to increasing crystallization and brecciation with time. For any cross section, calculations of rheologic parameters based on flow-margin depths generally gave lower values than those based on the dimensions of levees. This relationship may be attributed to the earlier formation and less complex evolution of the margins. The various equations gave more consistent results for upstream profiles, suggesting that calculations for remotely observed flows should avoid measurements near flow termini.     

Horizontal ground deformation patterns and magma storage during the Puu Oo eruption of Kilauea volcano,Hawaii: episodes 22–42

Horizontal ground deformation measurements were made repeatedly with an electronic distance meter near the Puu Oo eruption site approximately perpendicular to Kilauea's east rift zone (ERZ) before and after eruptive episodes 22–42. Line lengths gradually extended during repose periods and rapidly contracted about the same amount following eruptions. The repeated extension and contraction of the measured lines are best explained by the elastic response of the country rock to the addition and subsequent eruption of magma from a local reservoir. The deformation patterns are modeled to constrain the geometry and location of the local reservoir near Puu Oo. The observed deformation is consistent with deformation patterns that would be produced by the expansion of a shallow, steeply dipping dike just uprift of Puu Oo striking parallel to the trend of the ERZ. The modeled dike is centered about 800 m uprift of Puu Oo. Its top is at a depth of 0.4 km, its bottom at about 2.9 km, and the length is about 1.6 km; the dike strikes N65° E and dips at about 87°SE. The model indicates that the dike expanded by 11 cm during repose periods, for an average volumetric expansion of nearly 500 000 m³. The volume of magma added to the dike during repose periods was variable but correlates positively with the volume of erupted lava of the subsequent eruption and represents about 8% of the new lava extruded. Dike geometry and expansion values are used to estimate the pressure increase near the eruption site due to the accumulation of magma during repose periods. On average, vent pressures increased by about 0.38 MPa during the repose periods, one-third of the pressure increase at the summit. The model indicates that the dikelike body below Puu Oo grew in volume from 3 million cubic meters (Mm³) to about 10–12 Mm³ during the series of eruptions. The width of this body was probably about 2.5–3.0 m. No net long-term deformation was detected along the measured deformation lines.     

Magma dynamics, crystallization, and chemical differentiation of the 1959 Kilauea Iki lava lake, Hawaii, revisited

Using constraints from an extensive database of geological and geochemical observations along with results from fluid mechanical studies of convection in magma chambers, we identify the main physical processes at work during the solidification of the 1959 Kilauea Iki lava lakes. In turn, we investigate their quantitative influence on the crystallization and chemical differentiation of the magma, and on the development of the internal structure of the lava lake. In contrast to previous studies, vigorous stirring in the magma, driven predominately by the descent of dense crystal-laden thermal plumes from the roof solidification front and the ascent of buoyant compositional plumes due to the in situ growth of olivine crystals at the floor, is predicted to have been an inevitable consequence of very strong cooling at the roof and floor. The flow is expected to have caused extensive but imperfect mixing over most of the cooling history of the magma, producing minor compositional stratification at the roof and thermal stratification at the floor. The efficient stirring of the large roof cooling is expected to have resulted in significant internal nucleation of olivine crystals, which ultimately settled to the floor. Additional forcing due to either crystal sedimentation or the ascent of gas bubbles is not expected to have increased significantly the amount of mixing. In addition to convection in the magma, circulation driven by the convection of buoyant interstitial melt in highly permeable crystal-melt mushes forming the roof and the floor of the lava lake is envisaged to have produced a net upward flow of evolved magma from the floor during solidification. In the floor zone, mush convection may have caused the formation of axisymmetric chimneys through which evolved magma drained from deep within the floor into the overlying magma and potentially the roof. We hypothesize that the highly evolved, pipe-like ‘vertical olivine-rich bodies’ (VORBs) [Bull. Volcanol. 43 (1980) 675] observed in the floor zone, of the lake are fossil chimneys. In the roof zone, buoyant residual liquid both produced at the roof solidification front and gained from the floor as a result of incomplete convective mixing is envisaged to have percolated or ‘leaked‘ into the overlying highly-permeable cumulate, displacing less buoyant interstitial melt downward. The results from Rayleigh fractionation-type models formulated using boundary conditions based on a quantitative understanding of the convection in the magma indicate that most of the incompatible element variation over the height of the lake can be explained as a consequence of a combination of crystal settling and the extensive but imperfect convective mixing of buoyant residual liquid released from the floor solidification front. The remaining chemical variation is understood in terms of the additional influences of mush convection in the roof and floor on the vertical distribution of incompatible elements. Although cooling was concentrated at the roof of the lake, the floor zone is found to be thicker than the roof zone, implying that it grew more quickly. The large growth rate of the floor is explained as a consequence of a combination of the substantial sedimentation of olivine crystals and more rapid in situ crystallization due to both a higher liquidus temperature and enhanced cooling resulting from imperfect thermal and chemical mixing.     

Multispectral thermal infrared mapping of the 1 October 1988 Kupaianaha flow field,Kilauea volcano,Hawaii

Multispectral thermal infrared radiance measurements of the Kupaianaha flow field were acquired with the NASA airborne Thermal Infrared Multispectral Scanner (TIMS) on the morning of 1 October 1988. The TIMS data were used to map both the temperature and emissivity of the surface of the flow field. The temperature map depicted the underground storage and transport of lava. The presence of molten lava in a tube or tumulus resulted in surface temperatures that were at least 10° C above ambient. The temperature map also clearly defined the boundaries of hydrothermal plumes which resulted from the entry of lava into the ocean. The emissivity map revealed the boundaries between individual flow units within the Kupaianaha field. In general, the emissivity of the flows varied systematically with age but the relationship between age and emissivity was not unique. Distinct spectral anomalies, indicative of silica-rich surface materials, were mapped near fumaroles and ocean entry sites. This apparent enrichment in silica may have resulted from an acid-induced leaching of cations from the surfaces of glassy flows. Such incipient alteration may have been the cause for virtually all of the emissivity variations observed on the flow field, the spectral anomalies representing areas where the acid attack was most intense.     

Huge landslide blocks in the growth of piton de la fournaise, La réunion, and Kilauea volcano, Hawaii

Piton de la Fournaise, on the island of La Réunion, and Kilauea volcano, on the island of Hawaii, are active, basaltic shield volcanoes growing on the flanks of much larger shield volcanoes in intraplate tectonic environments. Past studies have shown that the average rate of magma production and the chemistry of lavas are quite similar for both volcanoes. We propose a structural similarity — specifically, that periodic displacement of parts of the shields as huge landslide blocks is a common mode of growth. In each instance, the unstable blocks are within a rift-zone-bounded, unbuttressed flank of the shield. At Kilauea, well-documented landslide blocks form relatively surficial parts of a much larger rift-zone-bounded block; scarps of the Hilina fault system mark the headwalls of the active blocks. At Fournaise, Hilina-like slump blocks are also present along the unbuttressed east coast of the volcano. In addition, however, the existence of a set of faults nested around the present caldera and northeast and southeast rift zones suggests that past chapters in the history of Fournaise included the slumping of entire rift-zone-bounded blocks themselves. These nested faults become younger to the east southeast and apparently record one of the effects of a migration of the focus of volcanism in that direction. Repeated dilation along the present set of northeast and southeast rift zones, most recently exemplified by an eruption in 1977, suggests that the past history of rift-zone-bounded slumping will eventually be repeated. The record provided by the succession of slump blocks on Fournaise is apparently at a relatively detailed part of a migration of magmatic focus that has advanced at least 30 km to the east-southeast from neighboring Piton des Neiges, an extinct Pliocene to Pleistocene volcano.