Simulations of the 2004 lava flow at Etna volcano using the magflow cellular automata model

Since the mechanical properties of lava change over time, lava flows represent a challenge for physically based modeling. This change is ruled by a temperature field which needs to be modeled. MAGFLOW Cellular Automata (CA) model was developed for physically based simulations of lava flows in near real-time. We introduced an algorithm based on the Monte Carlo approach to solve the anisotropic problem. As transition rule of CA, a steady-state solution of Navier-Stokes equations was adopted in the case of isothermal laminar pressure-driven Bingham fluid. For the cooling mechanism, we consider only the radiative heat loss from the surface of the flow and the change of the temperature due to mixture of lavas between cells with different temperatures. The model was applied to reproduce a real lava flow that occurred during the 2004–2005 Etna eruption. The simulations were computed using three different empirical relationships between viscosity and temperature.     

The morphology and evolution of the Stromboli 2002–2003 lava flow field: an example of a basaltic flow field emplaced on a steep slope

The use of a hand-held thermal camera during the 2002–2003 Stromboli effusive eruption proved essential in tracking the development of flow field structures and in measuring related eruption parameters, such as the number of active vents and flow lengths. The steep underlying slope on which the flow field was emplaced resulted in a characteristic flow field morphology. This comprised a proximal shield, where flow stacking and inflation caused piling up of lava on the relatively flat ground of the vent zone, that fed a medial–distal lava flow field. This zone was characterized by the formation of lava tubes and tumuli forming a complex network of tumuli and flows linked by tubes. Most of the flow field was emplaced on extremely steep slopes and this had two effects. It caused flows to slide, as well as flow, and flow fronts to fail frequently, persistent flow front crumbling resulted in the production of an extensive debris field. Channel-fed flows were also characterized by development of excavated debris levees in this zone (Calvari et al. 2005). Collapse of lava flow fronts and inflation of the upper proximal lava shield made volume calculation very difficult. Comparison of the final field volume with that expecta by integrating the lava effusion rates through time suggests a loss of ~70% erupted lava by flow front crumbling and accumulation as debris flows below sea level. Derived relationships between effusion rate, flow length, and number of active vents showed systematic and correlated variations with time where spreading of volume between numerous flows caused an otherwise good correlation between effusion rate, flow length to break down. Observations collected during this eruption are useful in helping to understand lava flow processes on steep slopes, as well as in interpreting old lava–debris sequences found in other steep-sided volcanoes subject to effusive activity.     

Factors controlling the lengths of channel-fed lava flows

Factors which control lava flow length are still not fully understood. The assumption that flow length as mainly influenced by viscosity was contested by Walker (1973) who proposed that the length of a lava flow was dependent on the mean effusion rate, and by Malin (1980) who concluded that flow length was dependent on erupted volume. Our reanalysis of Malin's data shows that, if short duration and tube-fed flows are eliminated, Malin's Hawaiian flow data are consistent with Walker's assertion. However, the length of a flow can vary, for a given effusion rate, by a factor of 7, and by up to 10 for a given volume. Factors other than effusion rate and volume are therefore clearly important in controlling the lengths of lava flows. We establish the relative importance of the other factors by performing a multivariate analysis of data for recent Hawaiian lava flows. In addition to generating empirical equations relating flow length to other variables, we have developed a non-isothermal Bingham flow model. This computes the channel and levee width of a flow and hence permits the advance rates of flows and their maximum cooling-limited lengths for different gradients and effusion rates to be calculated. Changing rheological properties are taken into account using the ratio of yield strength to viscosity; available field measurements show that this varies systematically from the vent to the front of a lava flow. The model gives reasonable agreement with data from the 1983–1986 Pu'u Oo eruptions and the 1984 eruption of Mauna Loa. The method has also been applied to andesitic and rhyolitic lava flows. It predicts that, while the more silicic lava flows advance at generally slower rates than basaltic flows, their maximum flow lengths, for a given effusion rate, will be greater than for basaltic lava flows.     

New formulae for estimating lava flow volumes at Mt. Etna Volcano, Sicily

 A new data set of Etna lava flows erupted since 1868 has been compiled from eight topographic maps of the volcano published at intervals since then. Volumes of 59 flows or groups of flows were measured from topographic difference maps. Most of these volumes are likely to be considerably more accurate than those published previously. We cut the number of flow volumes down to 25 by selecting those examples for which the volume of an individual eruption could be derived with the highest accuracy. This refined data set was searched for high correlations between flow volume and more directly measurable parameters. Only two parameters showed a correlation coefficient of 70% or greater: planimetric flow area A (70%) and duration of the eruption D (79%). If only short duration (<18 days) flows were used, flow length cubed, L³, had a correlation coefficient of 98%. Using combinations of measured parameters, much more significant correlations with volume were found. Dh had a correlation coefficient of 90% (h is the hydrostatic head of magma above the vent), and  , 92% (where W is mean width and E is the degree of topographic enclosure), and a combination of the two , 97%. These latter formulae were used to derive volumes of all eruptions back to 1868 to compare with those from the complete data set. Values determined from the formulae were, on average, lower by 16% (Dh), 7% (, and 19% . Received: 30 November 1998 / Accepted: 20 June 1999     

Development of the 1990 Kalapana Flow Field,Kilauea Volcano,Hawaii

The 1990 Kalapana flow field is a complex patchwork of tube-fed pahoehoe flows erupted from the Kupaianaha vent at a low effusion rate (approximately 3.5 m³/s). These flows accumulated over an 11-month period on the coastal plain of Kilauea Volcano, where the pre-eruption slope angle was less than 2°. the composite field thickened by the addition of new flows to its surface, as well as by inflation of these flows and flows emplaced earlier. Two major flow types were identified during the development of the flow field: large primary flows and smaller breakouts that extruded from inflated primary flows. Primary flows advanced more quickly and covered new land at a much higher rate than breakouts. The cumulative area covered by breakouts exceeded that of primary flows, although breakouts frequently covered areas already buried by recent flows. Lava tubes established within primary flows were longer-lived than those formed within breakouts and were often reoccupied by lava after a brief hiatus in supply; tubes within breakouts were never reoccupied once the supply was interrupted. During intervals of steady supply from the vent, the daily areal coverage by lava in Kalapana was constant, whereas the forward advance of the flows was sporadic. This implies that planimetric area, rather than flow length, provides the best indicator of effusion rate for pahoehoe flow fields that form on lowangle slopes.     

Sequence and eruptive style of the 1783 eruption of Asama Volcano, central Japan: a case study of an andesitic explosive eruption generating fountain-fed lava flow, pumice fall, scoria flow and forming a cone

The 3-month long eruption of Asama volcano in 1783 produced andesitic pumice falls, pyroclastic flows, lava flows, and constructed a cone. It is divided into six episodes on the basis of waxing and waning inferred from records made during the eruption. Episodes 1 to 4 were intermittent Vulcanian or Plinian eruptions, which generated several pumice fall deposits. The frequency and intensity of the eruption increased dramatically in episode 5, which started on 2 August, and culminated in a final phase that began on the night of 4 August, lasting for 15 h. This climactic phase is further divided into two subphases. The first subphase is characterized by generation of a pumice fall, whereas the second one is characterized by abundant pyroclastic flows. Stratigraphic relationships suggest that rapid growth of a cone and the generation of lava flows occurred simultaneously with the generation of both pumice falls and pyroclastic flows. The volumes of the ejecta during the first and second subphases are 0.21 km³ (DRE) and 0.27 km³ (DRE), respectively. The proportions of the different eruptive products are lava: cone: pumice fall=84:11:5 in the first subphase and lava: cone: pyroclastic flow=42:2:56 in the second subphase. The lava flows in this eruption consist of three flow units (L1, L2, and L3) and they characteristically possess abundant broken phenocrysts, and show extensive "welding" texture. These features, as well as ghost pyroclastic textures on the surface, indicate that the lava was a fountain-fed clastogenic lava. A high discharge rate for the lava flow (up to 10⁶ kg/s) may also suggest that the lava was initially explosively ejected from the conduit. The petrology of the juvenile materials indicates binary mixing of an andesitic magma and a crystal-rich dacitic magma. The mixing ratio changed with time; the dacitic component is dominant in the pyroclasts of the first subphase of the climactic phase, while the proportion of the andesitic component increases in the pyroclasts of the second subphase. The compositions of the lava flows vary from one flow unit to another; L1 and L3 have almost identical compositions to those of pyroclasts of the first and second subphases, respectively, while L2 has an intermediate composition, suggesting that the pyroclasts of the first and second subphases were the source of the lava flows, and were partly homogenized during flow. The complex features of this eruption can be explained by rapid deposition of coarse pyroclasts near the vent and the subsequent flowage of clastogenic lavas which were accompanied by a high eruption plume generating pumice falls and/or pyroclastic flows.Editorial responsibility: T. Druitt     

PRELIMINARY VOLCANIC HAZARD ZONATION IN JINLONGDINGZI VOLCANO,LONGANG VOLCANO AREA,JILIN PROVINCE,CHINA

Longgang volcano cluster is 150km away from the Tianchi volcano, located in Jingyu and Huinan Counties, Jilin Province, China. It had a long active history and produced hundreds of volcanoes. The latest and largest eruption occurred between 1 500 and 1 600 years ago by Jinlongdingzi(JLDZ)volcano which had several eruptions in the history. This paper discusses the volcanic hazard types, and using the numerical simulations of lava flow obtained with the Volcflow model, proposes the hazard zonation of JLDZ volcano area. JLDZ volcano eruption type is sub-plinian, which produced a great mass of tephra fallout, covering an area of 260km². The major types of volcanic hazards in JLDZ area are lava flow, tephra fallout and spatter deposits. Volcflow is developed by Kelfoun for the simulation of volcanic flows. The result of Volcflow shows that the flows are on the both sides of the previous lava flows which are low-lying areas now. According to the physical parameters of historical eruption and Volcflow, we propose the preliminary volcanic hazard zonation in JLDZ area. The air fall deposits are the most dangerous product in JLDZ. The highly dangerous region of spatter deposits is limited to a radius of about 2km around the volcano. The high risk area of tephra fallout is between 2km to 9km around the volcano, and between 9km to 14km is the moderate risk area. Out of 14km, it is the low risk area. Lava flow is controlled by topography. From Jinchuan Town to Houhe Village near the volcano is the low-lying area. If the volcano erupts, these areas will be in danger.     

Pahoehoe and aa in Hawaii: volumetric flow rate controls the lava structure

The historical records of Kilauea and Mauna Loa volcanoes reveal that the rough-surfaced variety of basalt lava called aa forms when lava flows at a high volumetric rate (>5–10 m³/s), and the smooth-surfaced variety called pahoehoe forms at a low volumetric rate (<5–10 m³/s). This relationship is well illustrated by the 1983–1990 and 1969–1974 eruptions of Kilauea and the recent eruptions of Mauna Loa. It is also illustrated by the eruptions that produced the remarkable paired flows of Mauna Loa, in which aa formed during an initial short period of high discharge rate (associated with high fountaining) and was followed by the eruption of pahoehoe over a sustained period at a low discharge rate (with little or no fountaining). The finest examples of paired lava flows are those of 1859 and 1880–1881. We attribute aa formation to rapid and concentrated flow in open channels. There, rapid heat loss causes an increase in viscosity to a threshold value (that varies depending on the actual flow velocity) at which, when surface crust is torn by differential flow, the underlying lava is unable to move sufficiently fast to heal the tear. We attribute pahoehoe formation to the flowage of lava at a low volumetric rate, commonly in tubes that minimize heat loss. Flow units of pahoehoe are small (usually <1 m thick), move slowly, develop a chilled skin, and become virtually static before the viscosity has risen, to the threshold value. We infer that the high-discharge-rate eruptions that generate aa flows result from the rapid emptying of major or subsidiary magma chambers. Rapid near-surface vesiculation of gas-rich magma leads to eruptions with high discharge rates, high lava fountains, and fast-moving channelized flows. We also infer that long periods of sustained flow at a low discharge rate, which favor pahoehoe, result from the development of a free and unimpeded pathway from the deep plumbing system of the volcano and the separation of gases from the magma before eruption. Achievement of this condition requires one or more episodes of rapid magma excursion through the rift zone to establish a stable magma pathway.     

The initial cooling of pahoehoe flow lobes

 In this paper we describe a new thermal model for the initial cooling of pahoehoe lava flows. The accurate modeling of this initial cooling is important for understanding the formation of the distinctive surface textures on pahoehoe lava flows as well as being the first step in modeling such key pahoehoe emplacement processes as lava flow inflation and lava tube formation. This model is constructed from the physical phenomena observed to control the initial cooling of pahoehoe flows and is not an empirical fit to field data. We find that the only significant processes are (a) heat loss by thermal radiation, (b) heat loss by atmospheric convection, (c) heat transport within the flow by conduction with temperature and porosity-dependent thermal properties, and (d) the release of latent heat during crystallization. The numerical model is better able to reproduce field measurements made in Hawai'i between 1989 and 1993 than other published thermal models. By adjusting one parameter at a time, the effect of each of the input parameters on the cooling rate was determined. We show that: (a) the surfaces of porous flows cool more quickly than the surfaces of dense flows, (b) the surface cooling is very sensitive to the efficiency of atmospheric convective cooling, and (c) changes in the glass forming tendency of the lava may have observable petrographic and thermal signatures. These model results provide a quantitative explanation for the recently observed relationship between the surface cooling rate of pahoehoe lobes and the porosity of those lobes (Jones 1992, 1993). The predicted sensitivity of cooling to atmospheric convection suggests a simple field experiment for verification, and the model provides a tool to begin studies of the dynamic crystallization of real lavas. Future versions of the model can also be made applicable to extraterrestrial, submarine, silicic, and pyroclastic flows. Received: 26 November 1994 / Accepted: 1 December 1995     

A dynamical model of lava flows cooling by radiation

The behaviour of a lava flow is reproduced by a two-dimensional model of a Bingham liquid flowing down a uniform slope. Such a liquid is described by two rheological parameters, yield stress and viscosity, both of which are strongly temperature-dependent. Assuming a flow rate and an initial temperature of the liquid at the eruption vent, the temperature decrease due to heat radiation and the consequent change in the rheological parameters are computed along the flow. Both full thermal mixing and thermal unmixing are considered. The equations of motion are solved analytically in the approximation of a slow downslope change of the flow parameters. Flow height and velocity are obtained as functions of the distance from the eruption vent; the time required for a liquid element to reach a certain distance from the vent is also computed. The gross features of observed lava flows are reproduced by the model which allows us to estimate the sensitivity of flow dynamics to changes in the initial conditions, ground slope and rheological parameters. A pronounced increase in the rate of height increase and velocity decrease is found when the flow enters the Bingham regime. The results confirm the observation according to which lava flows show an initial rapid advance, followed by a marked deceleration, while the final length of a flow is such that the Graetz number is in the order of a few hundreds.     

Two types of lava fields and the mechanism of their generation

Two types of lava fields, elongate and compact, are generated by two types of discrete lava flows that are designated as wide and narrow in this paper. The width of a wide flow considerably exceeds its thickness, the converse being true for a narrow one. Wide flows produce relatively narrow elongate lava fields, while narrow ones give rise to wider and compact lava fields of relatively smaller area and greater thickness, which is related to peculiarities in their dynamics causing branching and greater heat losses in such flows. Exact analytical formulas have been derived to relate the rate of lava in a flow to the flow geometry and lava rheology for idealized models of both flow types; a combination of these yields the approximate relationships for actual flows. We have found conditions for the generation of narrow flows showing strongly nonlinear behavior, which produces the peculiar features of the lava fields they create. Examples of lava fields generated by narrow flows from lateral vents on Klyuchevskoi Volcano are given. The conditions for their generation are more closely fulfilled in the upper part of the volcano’s cone. It is hypothesized that such an arrangement of lava fields may facilitate slope instabilities and enhance the likelihood of major landslides.     

Emplacement conditions of the c. 1,600-year bp Collier Cone lava flow, Oregon: a LiDAR investigation

A long-standing question in lava flow studies has been how to infer emplacement conditions from information preserved in solidified flows. From a hazards perspective, volumetric flux (effusion rate) is the parameter of most interest for open-channel lava flows, as the effusion rate is important for estimating the final flow length, the rate of flow advance, and the eruption duration. The relationship between effusion rate, flow length, and flow advance rate is fairly well constrained for basaltic lava flows, where there are abundant recent examples for calibration. Less is known about flows of intermediate compositions (basaltic andesite to andesite), which are less frequent and where field measurements are limited by the large block sizes and the topographic relief of the flows. Here, we demonstrate ways in which high-resolution digital topography obtained using Light Detection and Ranging (LiDAR) systems can provide access to terrains where field measurements are difficult or impossible to collect. We map blocky lava flow units using LiDAR-generated bare earth digital terrain models (DTMs) of the Collier Cone lava flow in the central Oregon Cascades. We also develop methods using geographic information systems to extract and quantify morphologic features such as channel width, flow width, flow thickness, and slope. Morphometric data are then analyzed to estimate both effusion rates and emplacement times for the lava flow field. Our data indicate that most of the flow outline (which comprises the earliest, and most voluminous, flow unit) can be well explained by an average volumetric flux ～14–18?m³/s; channel data suggest an average flux ～3?m³/s for a later, channel-filling, flow unit. When combined with estimates of flow volume, these data suggest that the Collier Cone lava flow was most likely emplaced over a time scale of several months. This example illustrates ways in which high-resolution DTMs can be used to extract and analyze morphologic measurements and how these measurements can be analyzed to estimate emplacement conditions for inaccessible, heavily vegetated, or extraterrestrial lava flows.     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

Steady-state operation of Stromboli volcano,Italy: constraints on the feeding system

Stromboli volcano has been in continuous eruption for several thousand years without major changes in the geometry and feeding system. The thermal structure of its upper part is therefore expected to be close to steady state. In order to mantaim explosive activity, magma must release both gas and heat. It is shown that the thermal and gas budgets of the volcano lead to consistent conclusions. The thermal budget of the volcano is studied by means of a finite-element numerical model under the assumption of conduction heat transfer. It is found that the heat loss through the walls of an eruption conduit is weakly sensitive to the dimensions of underlying magma reservoirs and depends mostly on the radius and length of the conduit. In steady state, this heat loss must be balanced by the cooling of magma which flows through the system. For the magma flux of about 1 kg s^-1 corresponding to normal Strombolian activity, this requires that the conduits are a few meters wide and not deeper than a few hundred meters. This implies the existence of a magma chamber at shallow depth within the volcanic edifice. This conclusion is shown to be consistent with considerations on the thermal effects of degassing. In a Strombolian explosion, the mass ratio of gas to lava is very large, commonly exceeding two, which implies that the thermal evolution of the erupting mixture is dominated by that of the gas phase. The large energy loss due to decompression of the gas phase leads to decreased eruption temperatures. The fact that lava is molten upon eruption implies that the mixture does not rise from more than about 200 m depth. To sustain the magmatic and volcanic activity of Stromboli, a mass flux of magma of a few hundred kilograms per second must be supplied to the upper parts of the edifice. This represents either the rate of magma production from the mantle source feeding the volcano or the rate of magma overturn in the interior of a large chamber.     

Experimental study of dense pyroclastic density currents using sustained,gas-fluidized granular flows

Submarine lava flow morphology is commonly used to estimate relative flow velocity, but the effects of crystallinity and viscosity are rarely considered. We use digital petrography and quantitative textural analysis techniques to determine the crystallinity of submarine basaltic lava flows, using a set of samples from previously mapped lava flow fields at the hotspot-affected Galápagos Spreading Center. Crystallinity measurements were incorporated into predictive models of suspension rheology to characterize lava flow consistency and rheology. Petrologic data were integrated to estimate bulk lava viscosity. We compared the crystallinity and viscosity of each sample with its flow morphology to determine their respective roles in submarine lava emplacement dynamics. We find no correlation between crystallinity, bulk viscosity, and lava morphology, implying that flow advance rate is the primary control on submarine lava morphology. However, we show systematic variations in crystal size and shape distribution among pillows, lobates, and sheets, suggesting that these parameters are important indicators of eruption processes. Finally, we compared the characteristics of lavas from two different sampling sites with contrasting long-term magma supply rates. Differences between lavas from each study site illustrate the significant effect of magma supply on the physical properties of the oceanic upper crust.     

Temperature of an active lava channel from spectral measurements,Kilauea Volcano,Hawaii

A narrow band spectroradiometer was used to determine the characteristic temperatures of a very active channeled lava flow for the phase 50 eruption of Pu'u 'O'o on the East Rift Zone of Kilauea Volcano, Hawaii. During the twilight of 19 February 1992, 14 spectra of this activity were acquired over a 51 minute interval [18.29 to 19.20 Hawaiian Standard Time (HST)], from which the thermal distribution of energy of two 18 m² areas, one near the center and one near the margin of the flow, may be investigated. A twocomponent thermal mixing model applied to the data taken of the center of the channel gave, in the most powerful instance (1.8x10⁵ W/m²), a crust temperature of 940° C, a hot component temperature of 1120°C and a hot radiating area of 60% of the total area. A simultaneous spectrum acquired near the channeled flow margin yielded a crust temperature of 586° C and a hot area of only 1.2% of the total area radiating at 1130° C. Average radiant flux densities recorded for the center of the lava channel (1.3x10⁵ W/m² average) are much greater than previous measurements of lava lakes (4.9x10³ W/m²) or recently emplaced lava flows (maximum of 7.2x10⁴ W/m²). The energetic nature of this eruption is shown by satellite measurements made at 02.33 HST on 22 February 1992 by the Advanced Very High Resolution Radiometer in Band 2 (0.72–1.10 m). These show the utility of using existing satellites with moderate resolution (1 km x 1 km pixels) and high temporal coverage (eight overpasses each day for Hawaii) as potential thermal alarms for rapidly assessing the hazard potential of large volcanic eruptions.     

Development of unconfined historic lava flow fields in Tenerife: implications for the mitigation of risk from a future eruption

Historic and recent (last 2,000?years) eruptions on the active volcanic island of Tenerife have been predominantly effusive, indicating that this is the most probable type of activity to be expected in the near future. In the past, lava flow invasion caused major damage on the island, and as the population and infrastructure have increased dramatically since the last eruption, lava flows are the most important short-term volcanic risk on Tenerife. Hence, an understanding of lava flow behaviour is vital to manage risks from lava flows and minimise future losses on the island. This paper focuses on the lava flows from the historic eruptions in Tenerife, providing new data on the volumes emitted, advance rates and the timing of the emplacement of flows. The studies show three main stages in the development of unconfined flow fields: the first stage, corresponding to the fast advance of the initial fronts during the first 24?C36?h of eruption (reaching calculated velocities of up to 1.1?m/s); the second stage, in which fronts stagnate; and a third stage, in which secondary lava flows develop from breakouts 4?C7?days after the initial eruption and farther extend the flow field (velocities of up to 0.02?m/s have been calculated for this stage). The breakouts identified originated at sites both proximal and distal to the vent and, in both cases, caused damage through lengthening and widening the original flow field. Hence, the probability of damage from lavas to land and property is highest during stages 1 and 3, and this should be accounted for when planning the response to a future effusive eruption. Tenerife??s lava flows display a similar behaviour to that of lava flows on volcanoes characterised by basaltic effusive activity (such as Etna or Kilauea), indicating the possibility of applying forecasting models developed at those frequently active volcanoes to Tenerife.     

Activity of Hawaiian Volcanoes during the years 1940–1950

Summit eruptions of Mauna Loa, on the Island of Hawaii, occurred in 1940 and 1949, and flank eruptions in 1942 and 1950. Lava poured out in 1940 and 1942 was about equal in amount, totaling approximately 76 million cubic meters in each eruption. The 1949 eruption was somewhat smaller, liberating approximately 59 million cubic meters. The 1950 eruption was one of the largest on record, producing five large lava flows and several smaller ones, totaling approximately 459 million cubic meters. Three of the 1950 flows entered the sea. In 1942 a lava flow threatened the city of Hilo, and was bombed from the air in an effort to divert it. Calculations indicate that the gas content of the lava extruded during the 1940 eruption probably was in the vicinity of one percent by weight of the total magma. Other calculations indicate the viscosity of fluid Hawaiian lava to be in the range of 10³ to 10⁵ poises. Temperature readings on the 1950 lava ranged from 1090⁰ to 900⁰ C. Kilauea Volcano showed signs of uneasiness in 1944, with an apparent increase of magmatic pressure indicated by outward tilting of the moutain flanks and a series of earthquakes progressing toward the surface. In December 1950 a series of earthquakes accompanied a subsidence of the summit of Kilauea Volcano.     

Volcanoes of the Cape Hoskins area,New Britain,territory of Papua and New Guinea

One active and ten extinct Quaternary volcanoes are described from the Cape Hoskins area, on the north coast of New Britain. They are mostly strato volcanoes built up of lava flows, lava domes, pyroclastic flows, lahars, tephra, and derived alluvial sediments. The volcanic products range in composition from basalt to rhyolite, but basaltic andesite and andesite predominate. Much of the area is covered by tephra, several metres thick, consisting mainly of rhyolitic pumice. The active volcano, Pago, is built up of several glacier-like lava flows, the last of which was formed during an eruption in 1914–18. Pago lies within a well-preserved caldera forming the central part of a broad low-angle cone, named Witori, which consists largely of welded and unwelded pyroclastic flow deposits. C-14 dates obtained on charcoal indicate that the caldera eruption occurred about 2500 years B. P. Another caldera of similar age lies south of Witori. Of the other eight volcanoes described four are relatively well-preserved steep-sided cones formed mainly of lava flows, one is a remnant of a low-angle cone with a caldera, and three are deeply eroded cones which have none of their constructional surfaces preserved.     

The effect of crystallization on the rheology and dynamics of lava flows

The dynamics of a lava flow is studied by a two-dimensional model describing a viscous fluid with Bingham rheology, flowing down a slope. The temperature in the flow is calculated assuming that heat is transferred through the plug by conduction and is lost by radiation to the atmosphere at the top of the flow. Taken into account is that the increasing crystallization takes place in the flow as a consequence of cooling. The lava viscosity and yield stress are expressed as a function of crystallization degree as well as of temperature: in particular it is assumed that yield stress reaches a maximum value above the solidus temperature, according to experimental data. Dynamical variables, such as velocity and thickness of the flow, are calculated for different values of the maximum crystallization degree and the flow rate. The model shows how the lava flow dynamics is affected by cooling and crystallization. The cooling of the flow is controlled by the increase of yield stress, which produces a thicker plug and makes the heat loss slower. The increasing crystallization has two opposing effects on viscosity: it produces an increase of viscosity, but at the same time produces an increase of yield stress and hence reduces the heat loss and keeps the internal temperature high. As a consequence, lava flows are significantly affected by the dependence of yield stress on temperature and scarcely by the maximum crystallization degree.