Damage to structures by pyroclastic flows and surges, inferred from nuclear weapons effects

In order to define the risk from explosive eruptions, one must constrain both the probability of explosive events and the effects, or consequences, of those events. This paper focuses on the effects of pyroclastic flows and surges (here termed ‘pyroclastic density currents', or PDCs) on buildings, infrastructure elements, and to some extent on vehicles. PDCs impart a lateral force to such structures in the form of dynamic pressure, which depends on the bulk density of the PDC (which in turn depends mainly on particle concentration) and its velocity. For reasonable ranges of particle concentration (10⁻³ to 0.5) and velocities (10 to 300 m/s), dynamic pressure on the upstream face of a structure ranges from 0.1 kPa to 10⁴ kPa. Lateral loads ranging up to about 100 kPa were produced during nuclear weapons tests in the 1940s and 1950s that were designed to study the effects of such loading on a variety of structures for civil defense and emergency response purposes in the event of nuclear war. Although considerable simplifications are involved, the data from these weapon tests provide useful analog information for understanding the effects of PDCs. I reviewed data from the nuclear tests, describing the expected damage from different loadings. Tables are provided that define the response of different structural elements (e.g., windows, framing, walls) and whole structures to loading in probabilistic terms, which in principle account for variations in construction quality, orientation, and other factors. Finally, damage documented from historical eruptions at Mt. Lamington (1951), Herculaneum (AD 79 Vesuvius eruption), and St. Pierre (1902 Mt. Pelee eruption) is reviewed. Damage patterns, combined with estimates of velocity, provide an independent estimate of particle concentration in the PDCs. Details of structural damage should be recorded and mapped around future eruptions in order to help refine this aspect of consequence analysis. Another fruitful approach would be to combine numerical simulations of eruption scenarios, which can produce simulated maps of dynamic pressure, with GIS-based data on structures for a given region; the result would be predictions of consequences that could be used for planning and emergency response training.     

Debris-avalanche deposits of the Milo Lahar sequence and the opening of the Valle del Bove on Etna volcano (Italy)

Previously undescribed debris-avalanche deposits occur in two locations downslope from the open end of the Valle del Bove. These outcrops comprise unstratified, ungraded deposits of metre-scale lava blocks in a matrix of weathered and fractured lava clasts. The avalanche deposits are unconformably overlain by matrix- to clast-supported conglomerates, representing debris-flow and interbedded fluvial deposits, that constitute most of the Milo Lahar sequence. We present evidence that the Milo Lahar sequence, which crops out just at the exit of the Valle del Bove, formed during the opening and enlargement of this depression. The presence of the avalanche deposits at the base of the Milo Lahar sequence indicates that catastrophic landslides were involved in the formation of the Valle del Bove. The composition of lavas in the debris avalanche deposits is similar to that of most of the Ellittico volcanic sequence exposed along the northern wall of the Valle del Bove. Radiocarbon dates of 8400 and 5300 years BP from the base and top, respectively, of the debris-flow sequence indicate that the Milo Lahars are correlative with the exposed part of the Chiancone deposit. The basal lahars of the Chiancone, which contain lava blocks whose compositions partially overlap that of blocks in the avalanche deposits, may have formed by water concentration in the distal end of the avalanche causing transformation to debris, or alternatively by reworking of the avalanche deposit.     

The basal ash deposit of the Sovana Eruption (Vulsini Volcanoes, central Italy): the product of a dilute pyroclastic density current

A low aspect ratio, decimeter-thick ash deposit, axisymmetrically distributed around the Latera Caldera (Western Vulsini Volcanoes, central Italy) has been studied by means of field and laboratory investigations. Field studies comprise facies analysis at centimeter scale and maximum clast size and deposit thickness measurements. Grain size and component distribution, chemical composition and particle morphoscopic features have been determined on selected samples. We discuss the co-ignimbrite ash fall vs. pyroclastic surge origin of the deposit and the hydrovolcanic vs. magmatic eruption nature. Complex facies association, textural features and grain size data rule out an ash fall origin for the whole deposit. The hydrovolcanic nature of the eruption has been discarded on the grounds of componentry and morphoscopic features of vitric fragments. We propose that the main body of the ash deposit formed from a radially expanding, dilute, turbulent pyroclastic density current, originated by a continuous collapse of a low-altitude (a few kilometers) eruptive column with a possible radial jet component.     

The Hasan Dagi stratovolcano (Central Anatolia, Turkey): evolution from calc-alkaline to alkaline magmatism in a collision zone

The Hasan Dagi volcano is one of the two large Plio-Quaternary volcanoes in Cappadocia (Central Anatolia, Turkey). Three stages of edifice construction have been identified for this volcano: Paleovolcano, Mesovolcano and Neovolcano. Most samples from Hasan Dagi volcano are calc-alkaline and define an almost complete trend from basaltic andesite to rhyolite. However, the more recent (Neovolcano) mafic samples are alkaline basalts. The mineralogical and geochemical characteristics of the oldest lavas (Keçikalesi (13 Ma) and Paleo-Hasan Dagi (7 Ma)) are significantly different from those of the younger lavas (Meso- and Neo-Hasan Dagi (<1 Ma)). Calcic plagioclase and pigeonite are typically observed in these older lavas. The Paleovolcano basalts are depleted in alkalis and display a tholeiitic tendency whereas the differentiated lavas are depleted in Na₂O but enriched in K₂O compared to younger lavas. There is an evolution through time towards higher TiO₂, Fe₂O₃^*, MgO, Na₂O and K₂O and lower Al₂O₃ and SiO₂ which is reflected in the basalt compositions. All the basalts display multi-element patterns typical of continental margin magmas with a significant enrichment in LILE (K, Rb, Ba and Th) and LREE and strong (Paleovolcano) to moderate (Meso- and Neovolcano) negative Nb, Zr and Ti anomalies. However, the younger basalts are the most enriched in incompatible elements, in agreement with their alkaline affinities and do not systematically display negative HFSE anomalies. REE data suggest an hydrous amphibole-bearing crystallization history for both Meso- and Neovolcano lavas. The distinction between the older and younger lavas is also apparent in trace element ratios such as Nb/Y, Ti/Y and Th/Y. These ratios indicate the role of a subducted component±crustal contamination in the genesis of the Hasan Dagi lavas, particularly for the oldest lavas (Keçikalesi and Paleo-Hasan Dagi). The decreasing influence of this component through time, over the last 6–7 m.y., has been accompanied by an increasing contribution of melt-enriched lithosphere. Although the range of variation of Sr, Nd and Pb isotopic ratios is small (0.70457–0.70515; 0.51262–0.51273; 18.80–18.94; 15.64–15.69; 38.87–39.10), it also reflects the evolution of the magma sources through time. Indeed, the youngest (Neovolcano) and most primitive basalts display significantly lower ⁸⁷Sr/⁸⁶Sr than the Paleo- and Mesovolcano basalts, whereas the Mesovolcano basalts display more radiogenic Pb than Paleovolcano samples. Magma mixing processes between initially heterogeneous and/or variably contaminated magmas may account for the genesis of the less differentiated and intermediate lavas (48–57% SiO₂). Meso- and Neovolcano differentiated lavas (60–68% SiO₂) are either derived from the analyzed basalts or from more primitive and more depleted magmas by fractional crystallization±some crustal contamination (AFC). Furthermore, the highly differentiated samples (72–75% SiO₂) are not strongly contaminated. The strong calc-alkaline character of Hasan Dagi lavas, in the absence of contemporaneous subduction, must reflect the heritage of the early subduction of the Afro–Arabian plate under the Eurasian plate. The evolution towards alkaline compositions through time is clearly related to the development of extensional tectonics in Central Anatolia in the Late Miocene.     

A model that helps explain Sr-isotope disequilibrium between feldspar phenocrysts and melt in large-volume silicic magma systems

Feldspar phenocrysts of silicic volcanic rocks are commonly in Sr-isotopic disequilibrium with groundmass. In some cases the feldspar is more radiogenic, and in others it is less radiogenic. Several explanations have been published previously, but none of these is able to accommodate both senses of disequilibrium. We present a model by which either more- or less-radiogenic feldspar (or even both within a single eruptive unit) can originate. The model requires a magma body open to interaction with biotite- and feldspar-bearing wall rock. Magma is incrementally contaminated as wall rock melts incongruently. Biotite preferentially melts first, followed by feldspar. Such melting behavior, which is supported by both field and experimental studies, first contaminates magma with a relatively radiogenic addition, followed by a less-radiogenic addition. Feldspar phenocrysts lag behind melt (groundmass of volcanic rock) in incorporating the influx of contaminant, thus resulting in Sr-isotopic disequilibrium between the crystals and melt. The sense of disequilibrium recorded in a volcanic rock depends on when eruption quenches the contamination process. This model is testable by isotopic fingerprinting of individual feldspar crystals. For a given set of geologic boundary conditions, specific core-to-rim Sr-isotopic profiles are expectable. Moreover, phenocrysts that nucleate at different times during the contamination process should record different and predictable parts of the history. Initial results of Sr-isotopic fingerprinting of sanidine phenocrysts from the Taylor Creek Rhyolite are consistent with the model. More tests of the model are desirable.     

Effusion rate and energy balance of Paricutin eruption (1943–1952), Michoacan, Mexico

The annual effusion rate of the Paricutin eruption shows an exponential decrease with time. This trend was observed during eruptions in other volcanoes and is explained by a simple model of the temporal changes of effusion rate proposed by Machado (1974). The regular trend of the effusion rate permits, in this case, a forecast of the length of the eruption. The energy balance of the eruption shows that part of the available energy was absorbed in a process of differentiation of the magma and this fact probably shortened the eruption.     

Deconstructing geomorphology: an appreciation of the contributions of Barbara A. Kennedy (1943–2014)

Barbara Kennedy has left only a small body of work, consisting principally of eight research papers, nine commentaries and two books. Its extent belies its importance. Read systematically, the work represents a sustained and important critique of the direction taken by mainstream geomorphology since the process‐focused reorientation of the mid‐twentieth century. Copyright © 2014 John Wiley & Sons, Ltd.     

Mechanism of energy transfer in the lava lake of Niragongo (Zaire), 1959–1977

Heat and mass transfer rates were studied at the Niragongo lava lake during two expeditions directed by H. Tazieff in 1959 and 1972. The results of this study are as follows:Heat is transferred to the surface of the lake by the movement of lava; gas discharge is a result and not the cause of convection. The chemical composition of the gases and magma has changed very little between 1959 and 1972, whereas the mass and energy outputs differ by an order of magnitude. In 1977 a catastrophic explosion seems to have been caused by tectonic factors, stopping the slow convection of magma under the volcano and hence reducing surface manifestations in the form of the lava lake and escaping fumarolic and magmatic gases. The gas discharge was, in tons day⁻¹, 5000 for H₂O, 11,000 for CO₂, 1000 for SO₂ in 1959, and in 1972 7700 for H₂O, 180,000 for CO₂ and 23,000 for SO₂. These values correspond to an energy transfer of 0.9 × 10⁹ W in 1959 and 16 × 10⁹ W in 1972.     

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.     

Trends in rainfall and temperature in the Peruvian Amazon–Andes basin over the last 40 years (1965–2007)

The hydroclimatology of the Peruvian Amazon–Andes basin (PAB) which surface corresponding to 7% of the Amazon basin is still poorly documented. We propose here an extended and original analysis of the temporal evolution of monthly rainfall, mean temperature (T_mean), maximum temperature (T_max) and minimum temperature (T_min) time series over two PABs (Huallaga and Ucayali) over the last 40 years. This analysis is based on a new and more complete database that includes 77 weather stations over the 1965–2007 period, and we focus our attention on both annual and seasonal meteorological time series. A positive significant trend in mean temperature of 0.09 °C per decade is detected over the region with similar values in the Andes and rainforest when considering average data. However, a high percentage of stations with significant T_mean positive trends are located over the Andes region. Finally, changes in the mean values occurred earlier in T_max (during the 1970s) than in T_min (during the 1980s). In the PAB, there is neither trend nor mean change in rainfall during the 1965–2007 period. However, annual, summer and autumn rainfall in the southern Andes presents an important interannual variability that is associated with the sea surface temperature in the tropical Atlantic Ocean while there are limited relationships between rainfall and El Niño‐Southern Oscillation (ENSO) events. On the contrary, the interannual temperature variability is mainly related to ENSO events. Copyright © 2012 John Wiley & Sons, Ltd.     

Post-emplacement tilting of lava flows inferred from magnetic fabric study: the example of Oligocene lavas in the Jeanne d’Arc Peninsula (Kerguelen Islands)

Statistical analysis of the magnetic fabric of samples from several successive lava flows emplaced under similar conditions can allow determination of the mean flow direction when magnetic fabric data from individual flows do not lead to reliable results. A difference between the obtained flow direction and the present dip direction indicates that the flows were tilted after emplacement. For 2 successive series of flows on the Jeanne d’Arc Peninsula presently NNW dipping, this method shows lava emplacement along a SSW–NNE direction. This indicates a gentle tilting acquired during a period of weak deformation in the whole archipelago. Additionally, the magnetic fabric data allow the reconstruction of the different conditions of emplacement of these two series of lava flows and of formation of the local thick conglomerate interbedded between these series.     

Contrasting regional discharge evolutions in the Amazon basin (1974–2004)

Former hydrological studies in the Amazon Basin generally describe annual discharge variability on the main stem. However, the downstream Amazon River only represents the mean state of the Amazonian hydrological system. This study therefore uses a new data set including daily discharge in 18 sub-basins to analyze the variability of regional extremes in the Amazon basin, after recalling the diversity of the hydrological annual cycles within the Amazon basin. Several statistical tests are applied in order to detect trends and breaks in the time series. We show that during the 1974–2004 period, the stability of the mean discharge on the main stem in Óbidos is explained by opposite regional features that principally involve Andean rivers: a decrease in the low stage runoff, particularly important in the southern regions, and an increase in the high stage runoff in the northwestern region. Both features are observed from the beginning of the nineties. These features are also observed in smaller meridian sub-basins in Peru and Bolivia. Moreover we show that the changes in discharge extremes are related to the regional pluriannual rainfall variability and the associated atmospheric circulation as well as to tropical large-scale climatic indicators.     

Using meteorological observer data to compare wind erosion during two great droughts in eastern Australia; the World War II Drought (1937–1946) and the Millennium Drought (2001–2010)

Australian meteorological observers started using the World Meteorological Organization (WMO) weather coding system in the 1950s. This system is still in use around the world today. However, observing and recording the weather in an organized and systematic manner had been ongoing for over 100 years prior to the adoption of this coding system, and much like Australia, most countries will have historical meteorological records. In this paper we compare the wind erosion of two of the greatest droughts in Australian recorded history; the World War II (WWII) Drought (1937–1945) and the Millennium Drought (2001–2009). To do this we analysed previously unavailable meteorological observer records from the Australian Bureau of Meteorology (ABM). Wind erosion records, mostly in long‐hand written form, were translated to the modern WMO coding system for the WWII Drought and compared with the wind erosion of Australia's recently‐ended Millennium Drought, one of the longest and harshest on record. We quantify wind erosion using Dust Event Days (DED) and a modified version of a published Dust Storm Index (DSI) to show that wind erosion during the WWII Drought was up to 4.6 times higher than during the Millennium Drought. This study has international significance because it demonstrates a methodology for tracking changes in wind erosion over the past 75 years based on observer records available in every country with a history of organized weather observation. Copyright © 2014 John Wiley & Sons, Ltd.     

Ion Chromatography (IC) as a Method in Studies on the Odra River Pollution with F–, Cl–, NO3–, PO43–, SO42– Ions

The paper presents the ion chromatography (IC) – DIONEX DX-100 system – in the studies on the Odra River water pollution with some inorganic anions. The precision of analysis, detection and quantification limits of each anion have been presented in order to evaluate the method. Results obtained with IC are compared with those using other analytical methods: ICP-AES for SO42–, and the argentometric titration for Cl–. It has been found that the IC method can be applied in the investigation of the anion pollution of the Odra River water: concentrations of Cl–, NO3– and SO42– can be determined with sufficient accuracy, while F– concentration, on the account of common presence of carbonate anions, can only be roughly estimated. The results have been compared with maximum permitted levels (MPL) for Polish water quality classes.     

Adsorbable Organic Halogens (AOX) in Polish Rivers – Levels and Changes

This paper presents results of a study on the concentration of adsorbable organic halogens (AOX) in water samples of the Oder River and its tributaries collected from 1998 to 2002. The Vistula River water for AOX was analysed in the period 2001–2002. The investigated waters of the Vistula near the mouth (Kiezmark) were more contaminated with AOX than the Oder near the mouth (Widuchowa, Krajnik Dolny). In general, AOX levels in the Polish rivers are comparable to the values found in the European rivers, e.g. the Rhine and the Ton. The research has shown that AOX were correlated with pH, total suspended solids, and the riverine flow. No relation between dissolved oxygen and redox potential and AOX was found out. During the flood period AOX, inorganic halogens (TIX), and mineral substances (sulfates) were analysed.     

Magnitude–frequency relationship for debris flows on the fan of the Chalance torrent,Valgaudemar (French Alps)

In 1996 a large debris flow occurred on the fan of the Chalance torrent system, a tributary of the Séveraisse river, French Alps. To investigate the magnitude and frequency of such debris flows on this fan, fieldwork was carried out in the summer of 1998. Detailed investigation revealed that several debris flows have occurred in the past 200 years. Lichenometry was used as a dating technique to obtain the frequency of these debris‐flow events. Also the volume of these flows was estimated. With these data a magnitude–frequency relationship was constructed. This relationship shows a maximum magnitude of at least 50 × 10³ m³. Based on data for the past c. 150 years, a debris flow of such a volume appears to have a recurrence interval of approximately 34 years. Copyright © 2002 John Wiley & Sons, Ltd.     

Zonation of the main volcanic hazards (lava flows and ash fall) in Tenerife, Canary Islands. A proposal for a surveillance network

The island of Tenerife is volcanically complex, and its eruptive history predominantly reflects the processes and products of two different eruptive styles: (1) non-explosive effusions of basaltic lavas from fissure vents mostly aligned along two ridges; and (2) less frequent but explosive salic eruptions from central vents associated with the Las Cañadas volcanic edifice and associated summit caldera. We have taken into account this fundamental distinction to develop a volcanic-hazards zonation (for lava flows and ash fall only) that includes: definition of the principal hazards; identification of the areas that have higher probability of containing emission centres; and numerical modelling of the vulnerable areas to be affected by volcanic hazards. Not only does the volcanic-hazards zonation map provide emergency-management officials with an updated assessment of the volcanic hazards, but it also represents a starting point for the preparation of a volcanic risk map for Tenerife. Finally, the hazards-zonation map also furnishes the basis for the design of a proposed volcano surveillance network.