Caldera formation at Volcán Colima, Mexico, by a large large holocene volcanic debris avalanche

About 4,300 years ago, 10 km³ of the upper cone of ancestral Volcán Colima collapsed to the southwest leaving a horseshoe-shaped caldera 4 km in diameter. The collapse produced a massive volcanic debris avalanche deposit covering over 1550 km² on the southern flanks of the volcano and extending at least 70 km from the former summit. The avalanche followed a steep topographic gradient unobstructed by barriers, resulting in an unusually high area/volume ratio for the Colima deposit. The apparent coefficient of friction (fall height/distance traveled) for the Colima avalanche is 0.06, a low value similar to those of other large-volume deposits. The debris avalanche deposit contains 40–75% angular volcanic clasts from the ancestral cone, a small proportion of vesicular blocks that may be juvenile, and in distal exposures, rare carbonate clasts plucked from the underlying surface by the moving avalanche. Clasts range in size to over 20 m in diameter and are brecciated to different degrees, pulverized, and surrounded by a rock-flour matrix. The upper surface of the deposit shows prominent hummocky topography with closed depressions and surface boulders. A thick, coarse-grained, compositionally zoned scoria-fall layer on the upper northeastern slope of the volcano may have erupted at the time of collapse. A fine-grained surge layer is present beneath the avalanche deposit at one locality, apparently representing an initial blast event. Most of the missing volume of the ancestral volcano has since been restored at an average rate of 0.002 km³/yr through repeated eruptions from the post-caldera cone. As a result, the southern slope of Volcán Colima may again be susceptible to collapse. Over 200,000 people are now living on primary or secondary deposits of the debris avalanche, and a repetition of this event would constitute a volcanic disaster of great magnitude.Ancestral Volcán Colima grew on the southern, trenchward flank of the earlier and larger volcano Nevado de Colima. Trenchward collapse was favored by the buttressing effect of Nevado, the rapid elevation drop to the south, and the intrusion of magma into the southern flank of the ancestral volcano. Other such trenchward-younging, paired volcanoes are known from Mexico, Guatemala, El Salvador, Chile, and Japan. The trenchward slopes of the younger cones are common sites for cone collapse to form avalanche deposits, as occurred at Colima and Popocatepetl in Mexico and at San Pedro Volcano in Chile.     

Integrated Geophysical Investigations at the Greek Kamarina Site (Southern Sicily,Italy)

Kamarina, located in southern Sicily (Italy), was an important Greek colony since its foundation in the sixth century BC. Archaeological excavations, carried out since the twentieth century, uncovered only limited portions of the site so far. Despite the importance of the Greek colony, the presence of remarkable buildings that archaeologists expected to bring to light has not found fully correspondence in the archaeological excavations. Consequently, the integrated geophysical prospection carried out in the study area is aimed to support and address the future archaeological investigations. After the photographic and thermographic survey obtained by an unmanned aerial vehicle, we performed a systematic survey through ground magnetic and GPR methods over an area of 6200 m². The acquisition procedures have been optimized in order to get the best results combining high resolution and elevated speed of acquisition. The results derived from the three geophysical techniques have been conveniently combined by means of a cluster analysis, allowing us to clearly identify a series of buried archaeological features. Because of their geometrical characteristics, often in good agreement with the spatial arrangement of the archaeological remains at the surface, these buried archaeological features can be interpreted as roads, walls, or buildings foundations in which the various construction phases of the city can be clearly recognized. The integrated approach has proven to be essential for a robust interpretation of the archaeogeophysical investigation.     

Seismic quantification of the explosions that destroyed the dome of Volcán de Colima, Mexico, in July–August 2003

In July–August 2003, the andesitic lava dome at Volcán de Colima, México, was destroyed by a sequence of explosions that replaced the 2×10⁶ m³ dome with a crater 200 m across and 30 m deep. The two strongest explosions occurred on July 17 and August 28. The initial low-frequency impulses that they produced, which were recorded on broadband seismic records, allowed an estimation of the counter forces of the initiating process as being equal to 0.3×10¹¹ N and 1×10¹¹ N for the July and August events, respectively. The seismic characteristics follow the Nishimura-Hamaguchi scaling law for volcanic explosions, reflecting self-similarity in the processes initiating explosive events. The results also show that counter forces can discriminate between the sizes of explosive eruptions that are assigned the same magnitude by conventional methods of classification such as the Volcanic Explosivity Index. The increasing use of broadband seismometers may therefore provide the basis for using counter forces to determine the magnitude of explosive eruptions.     

Archeointensity investigation on pottery vestiges from Puertas de Rolón, Capacha culture: In search for affinity with other Mesoamerican pre-Hispanic cultures

We report the results of a rock magnetic and archeointensity study on pottery fragments recovered at the archaeological site of Puertas de Rolón, Capacha culture (Colima), that is considered one of the most important pre-Classic cultural developments of western Mesoamerica, Mexico. Standard rock magnetic experiments were carried out in order to investigate the reliability of the ancient field strength recovered from studied materials. Six out of the seven analyzed fragments gave reliable intensity determinations. The mean archeointensity value obtained in this study is 35.3 ± 3.4 μT. It is in good agreement with another contemporaneous Mesoamerican mean archeointensity of 34.3 ± 1.1 μT obtained for El Opeño. Statistically similar rock-magnetic parameters between ceramics from El Opeño and Puertas de Rolón potsherds appear to be related to clay-sources with similar geological settings and weathering processes. Alternatively, ancient trade or interexchange between western (Capacha Culture) and eastern Mesoamerica locations, as it has been proposed earlier, seem to be supported by results obtained in this study.     

High resolution 3D ERT to help GPR data interpretation for researching archaeological items in a geologically complex subsurface

Muro Leccese (Lecce) contains one the most important Messapian archaeological sites in southern Italy.The archaeological interest of the site arises from the discovery of the remains of Messapian walls, tombs, roads, etc. (4th–2nd centuries BC) in the neighbourhood. The archaeological remains were found at about 0.3 m depth.At present the site belongs to the municipality, which intends to build a new sewer network through it. The risk of destroying potentially interesting ancient archaeological structures during the works prompted an archaeological survey of the area. The relatively large dimensions of the area (almost 10,000 m²), together with time and cost constraints, made it necessary to use geophysical investigations as a faster means to ascertain the presence of archaeological items. Since the most important targets were expected to be located at a soil depth of about 0.3 m, a ground-penetrating radar (GPR) survey was carried out in an area located near the archaeological excavations. Unfortunately the geological complexity did not allow an easy interpretation of the GPR data.Therefore a 3D electrical resistivity tomography (ERT) scan was conducted in order to resolve these interpretation problems.A three-way comparison of the results of the dense ERT measurements parallel to the x axis, the results of the measurements parallel to the y axis and the combined results was performed.Subsequently the synthetic model approach was used to provide a better characterization of the resistivity anomalies visible on the ERT field data.The 3D inversion results clearly illustrate the capability to resolve in view of quality 3D structures of archaeological interest. According to the presented data the inversion models along one direction (x or y) seems to be adequate in reconstructing the subsurface structures.Naturally field data produce good quality reconstructions of the archaeological features only if the x-line and y-line measurements are considered together. Despite the increased computational time required by the 3D acquisition and 3D inversion schemes, good quality results can be produced.     

Monitoring Colima Volcano,Mexico, using satellite data

The Colima Volcanic Complex at the western end of the Mexican Volcanic Belt is the most active andesitic volcano in Mexico. Short-wavelength infrared data from the Landsat Thematic Mapper satellite were used to determine the temperature and fractional area of radiant picture elements for two January data acquisitions in 1985 and 1986. The 1986 data showed four 28.5 m by 28.5 m pixels (picture elements) whose hot subpixel components had temperatures ranging from 511–774° C and areas of 1.8–13 m². The 1985 data had no radiating areas above background temperatures. Ground observations and measurements in November 1985 and February 1986 reported the presence of hot fumaroles at the summit with temperatures of 135–895° C. This study demonstrates the utility of satellite data for monitoring volcanic activity.     

The Pleistocene cinder cones surrounding Volcán Colima, Mexico re-visited: eruption ages and volumes, oxidation states, and sulfur content

Located at the volcanic front in the western Mexican arc, in the Colima Rift, is the active Volcán Colima, which lies on the southern end of the massive (∼450 km³) Colima-Nevado volcanic complex. Along the margins of this andesitic volcanic complex, is a group of 11 scoria cones and associated lavas, which have been dated by the ⁴⁰Ar/³⁹Ar method. Nine scoria cones erupted ∼1.3 km³ of alkaline magma (basanite, leucite-basanite, minette) between 450 and 60 ka, with >99% between 240 and 60 ka. Two additional cones (both the oldest and calc-alkaline) erupted <0.003 km³ of basalt (0.5 Ma) and <0.003 km³ of basaltic andesite (1.2 Ma), respectively. Cone and lava volumes were estimated with the aid of digital elevation models (DEMs). The eruption rate for these scoria cones and their associated lavas over the last 1.2 Myr is ∼1.2 km³/Myr, which is more than 400 times smaller than that from the andesitic Colima-Nevado edifice. In addition to these alkaline Colima cones, two other potassic basalts erupted at the volcanic front, but ∼200 km to the ESE (near the historically active Volcán Jorullo), and were dated at 1.06 and 0.10 Ma. These potassic suites reflect the tendency in the west-central Mexican arc for magmas close to the volcanic front to be enriched in K₂O relative to those farther from the trench.Ferric-ferrous analyses on pristine samples from the alkaline cones adjacent to V. Colima and V. Jorullo indicate that their oxygen fugacities relative to the nickel-nickel oxide buffer are significantly higher (ΔNN0=2–4) than those for the calc-alkaline magma types (0–1.5). These ΔNNO values correlate positively with Ba concentrations and likely reflect the influence of a slab-derived fluid. As a result of the high oxidation states, the solubility of sulfur in these potassic magmas is enhanced. Indeed the sulfur content of both the whole rock and the apatite phenocrysts (and in olivine melt inclusions reported in the literature) suggest that part of their pre-eruptive sulfur gas (SO₂) concentrations could have been discharged to the atmosphere in amounts comparable to the 1982 eruption of El Chichón, although over a prolonged period spanning thousands of years (not per eruption).Electronic Supplementary Material Supplementary material is available for this article at Editorial responsibility: J. Donnelly-Nolan     

Debris-flow-dominated sediment transport through a channel network after wildfire

Field studies that investigate sediment transport between debris-flow-producing headwaters and rivers are uncommon, particularly in forested settings, where debris flows are infrequent and opportunities for collecting data are limited. This study quantifies the volume and composition of sediment deposited in the arterial channel network of a 14-km² catchment (Washington Creek) that connects small, burned and debris-flow-producing headwaters (<1 km²) with the Ovens River in SE Australia. We construct a sediment budget by combining new data on deposition with a sediment delivery model for post-fire debris flows. Data on deposits were plotted alongside the slope–area curve to examine links between processes, catchment morphometry and geomorphic process domains. The results show that large deposits are concentrated in the proximity of three major channel junctions, which correspond to breaks in channel slope. Hyperconcentrated flows are more prominent towards the catchment outlet, where the slope–area curve indicates a transition from debris flow to fluvial domains. This shift corresponds to a change in efficiency of the flow, determined from the ratio of median grain size to channel slope. Our sediment budget suggests a total sediment efflux from Washington Creek catchment of 61 × 10³ m³. There are similar contributions from hillslopes (43 ± 14 × 10³ m³), first to third stream order channel (35 ± 12 × 10³ m³) and the arterial fourth to fifth stream order channel (31 ± 17 × 10³ m³) to the total volume of erosion. Deposition (39 ± 17 × 10³ m³) within the arterial channel was higher than erosion (31 ± 17 × 10³ m³), which means a net sediment gain of about 8 × 10³ m³ in the arterial channel. The ratio of total deposition to total erosion was 0.44. For fines <63 μm, this ratio was much smaller (0.11), which means that fines are preferentially exported. This has important implications for suspended sediment and water quality in downstream rivers. © 2019 John Wiley & Sons, Ltd.     

Magnetic properties and Archeointensity of Earth’s magnetic field recovered from El Opeño,earliest funeral architecture known in Western Mesoamerica

Despite of the impressive cultural heritage and abundant archaeological sites, absolute geomagnetic intensity data from Mesoamerica are still sparse. Archeointensity determinations using the Coe variant of the Thellier and Thellier method have been carried out on some selected pottery fragments from the El Opeño archeological site which has the earliest funeral architecture known in western Mesoamerica. The El Opeño chronology is supported by six C¹⁴ datings performed on carbon-bearing materials. Detailed rock magnetic experiments including susceptibility vs. temperature curves, hysteresis cycles and thermal demagnetization procedures were carried out in order to estimate the magnetic carriers and their stability. Cooling rate and anisotropy remanence corrected intensity values range from 25.0 ± 2.3 to 40.2 ± 3.0 μT and corresponding virtual axial dipole moments (VADM) range from 5.6 ± 0.5 to 8.9 ± 0.7 × 10²² Am². In addition, we present here a new compilation and analysis of existing absolute intensity data in order to try to estimate the variation of the Earth’s magnetic field over the past three millennia. The mean archeointensity obtained in the present study agrees reasonably well with the predicted absolute intensities retrieved from the CALS7K main field model. Other available Mesoamerican data, however, differ from this model. Most of available archeointensity data from Mesoamerica agree reasonably well with the ARCH3K main field model prediction. The broad peak defined at about 50 A.D. by our data is also predicted by the CALS3k.3 main field model but slightly displaced to the right while two smaller peaks are observed on the ARCH3K curve for the same time interval. The intensity value obtained at about 200 B.C. is a clear outlier and thus its geomagnetic significance should be confirmed by further investigations.     

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.     

The July 1994 episode of seismic activity at Colima Volcano, Mexico

A new period of seismic activity that culminated in a small phreatic explosion took place in Colima Volcano (Western Mexico) during the month of July 1994. In this note, we present our analysis of this seismicity based upon information from RESCO, the seismic network of the University of Colima. The activity began with a seismic swarm of type A (tectonic-like) earthquakes with epicenters towards the SSW of the summit, followed by shallow low-frequency events underneath the volcanic edifice. The activity was accompanied by landslides and culminated with an explosion that produced small ash falls on the surrounding area. The seismic activity ceased after this episode.     

Source mechanism of volcanic tremor: fluid-driven crack models and their application to the 1963 kilauea eruption

We propose a model for the mechanism of magma transport based on a fluid-filled tensile crack driven by the excess pressure of fluid. Such a transport mechanism can generate seismic waves by a succession of jerky crack extensions, if the fracture strength of rock varies in space, or if there is a difference between the dynamic and static values of the critical stress intensity factor. We also find that the opening and closing of a narrow channel connecting two fluid-filled cracks may be a source of seismic waves. Using the finite-difference method, we calculated the vibration of dry and fluid-filled cracks generated by: (1) a jerky extension at one end or at both ends and (2) a jerky opening of a narrow channel connecting two cracks. We then calculated the far-field and near-field radiation from these vibrating cracks. The spectra show peaked structures, but interestingly, most high-frequency peaks are only present in the near-field and cannot be transmitted to the far-field. The spectral features described above are often observed for volcanic tremors and in some cases for seismic signals associated with hydraulic fracturing experiments.We first consider as a model of volcanic tremor randomly occurring jerky crack extensions, and derive a formula relating the tremor amplitude to the excess pressure in the magma, the incremental area in each extension, and the frequency of extensions. These parameters are also constrained by other observations, such as the rate of magma flow.Our model was tested quantitatively against observations made in one of the best-described case histories of volcanic tremor: the October 5–6, 1963 Kilauea flank eruption. We found that a single, long crack extending from the summit to the eruptive site cannot explain the observations. The model of a steadily expanding crack ran into difficulties when quantitative comparisons were made with observations. The extension of crack area needed to explain the amplitude of volcanic tremor should accompany a large increase in tremor period which was not observed.Our second model is a chain of cracks connected by narrow channels which open and close. The length of each crack is around 1 km, the channel area connecting neighboring cracks is about 10³m², and the channel opens jerkily with the magmatic excess pressure of about 20 bars. The frequency of jerky opening of each channel is about once in 15 seconds. The channel is closed after each jerky opening, as soon as magma is moved through the channel.     

Hydraulic modelling of Athabasca Vallis,Mars

Abstract

Palaeohydraulic modelling is presented for Athabasca Vallis, the youngest known catastrophic flood channel on Mars. This modelling incorporates three significant advantages over previous modelling of Martian channels: a step-backwater hydraulic model; more accurate topography; and improved flood height indicators. The maximum modelled palaeodischarge is between 1 × 10⁶ and 2 × 10⁶m³s^?1 depending on the friction coefficient selected. An anomalously high palaeostage indicator suggests a region of ponded backwater in the channel in which streamlined forms were created through deposition, with the additional possibility of post-flood subsidence/lowering of the channel slope due to magma extrusion.     

Gas hazard assessment in a densely inhabited area of Colli Albani Volcano (Cava dei Selci, Roma)

The northwestern flank of the Colli Albani, a Quaternary volcanic complex near Rome, is characterised by high pCO₂ values and Rn activities in the groundwater and by the presence of zones with strong emission of gas from the soil. The most significant of these zones is Cava dei Selci where many houses are located very near to the gas emission site. The emitted gas consists mainly of CO₂ (up to 98 vol%) with an appreciable content of H₂S (0.8–2%). The He and C isotopic composition indicates, as for all fluids associated with the Quaternary Roman and Tuscany volcanic provinces, the presence of an upper mantle component contaminated by crustal fluids associated with subducted sediments and carbonates. An advective CO₂ flux of 37 tons/day has been estimated from the gas bubbles rising to the surface in a small drainage ditch and through a stagnant water pool, present in the rainy season in a topographically low central part of the area. A CO₂ soil flux survey with an accumulation chamber, carried out in February–March 2000 over a 12 000 m² surface with 242 measurement points, gave a total (mostly conductive) flux of 61 tons/day. CO₂ soil flux values vary by four orders of magnitude over a 160-m distance and by one order of magnitude over several metres. A fixed network of 114 points over 6350 m² has been installed in order to investigate temporal flux variations. Six surveys carried out from May 2000 to June 2001 have shown large variations of the total CO₂ soil flux (8–25 tons/day). The strong emission of CO₂ and H₂S, which are gases denser than air, produces dangerous accumulations in low areas which have caused a series of lethal accidents to animals and one to a man. The gas hazard near the houses has been assessed by continuously monitoring the CO₂ and H₂S concentration in the air at 75 cm from the ground by means of two automatic stations. Certain environmental parameters (wind direction and speed; atm P, T, humidity and rainfall) were also continuously recorded. At both stations, H₂S and CO₂ exceeded by several times the recommended concentration thresholds. The highest CO₂ and H₂S values were recorded always with wind speeds less than 1.5 m/s, mostly in the night hours. Our results indicate that there is a severe gas hazard for people living near the gas emission site of Cava dei Selci, and appropriate precautionary and prevention measures have been recommended both to residents and local authorities.     

Geology and hazard implications of the Maraunot notch in the Pinatubo Caldera,Philippines

The 1991 Pinatubo eruption left 5–6 km³ of debris on the volcano slopes, much of which has been mobilized into large lahars in the following rainy seasons. Also during the eruption, collapse, localized in part along preexisting faults, left a caldera 2.5 km in diameter that almost immediately began to accumulate a 1.6 × 10⁸ m³ lake. By 2001, the water had risen to the fault-controlled Maraunot Notch, the lowest, northwestern portion of the caldera rim comprising the physiographic sill of the Caldera Lake. That year, a narrow artificial canal dug into an old volcanic breccia underlying the outlet channel failed to induce a deliberate lake breakout, but discharge from heavy rains in July 2002 rapidly deepened the notch by 23 m, releasing an estimated 6.5 × 10⁷ m³ of lake water that bulked up into lahars with a volume well in excess of 1.6 × 10⁸ m³. Lakes in other volcanoes have experienced multiple breakouts, providing practical motivation for this study. Fieldwork and high-resolution digital elevation models reveal andesites and ancient lacustrine deposits, strongly fractured and deformed along a segment of the Maraunot Fault, a prominent, steeply dipping, left-lateral fault zone that trends N35°–40°W within and parallel to the notch. Seismicity in 1991 demonstrated that the Maraunot Fault is still active. The fault zone appears to have previously been the erosional locus for a large channel, filled with avalanche or landslide deposits of an earlier eruption that were exhumed by the 2002 breakout floods. The deformed lacustrine sediments, with an uncalibrated ¹⁴C age of 14,760 ± 40 year BP from a single charcoal sample, attest to the existence of an earlier lake, possibly within the Tayawan Caldera, rim remnants of which survive as arcuate escarpments. That lake may well have experienced one or more ancient breakouts as well. The 2002 event greatly reduced the possibility of another such event by scouring away the erodible breccia, leaving less erodible fractured andesites and lacustrine rocks, and by enlarging the outlet channel and its discharge capacity. Several lines of evidence indicate, however, that future lahar-generating lake breakouts at the notch may keep populations of Botolan municipality downstream at risk: (1) a volume of 9.5 × 10⁷ m³ of lake water remains perched 0.8 km above sea level; (2) seismicity in 1991 demonstrated that the Maraunot Fault is still active and movements of sufficient magnitude could enlarge the outlet and the discharge through it; (3) more likely, however, with or without earthquake activity, landslides from the steep to overhanging channel walls could block the channel again, and a major rainstorm could then cause a rise in lake level and sudden breakouts; (4) intrusion of a new dome into the bottom of the lake, possibly accompanied by phreatic explosions, could expel large volumes of lahar-generating water.     

Evaluation of digital channel network derivation methods in a glaciated subalpine catchment

The digital elevation model (DEM) has become an essential tool for an increasing array of mountain runoff analyses, particularly the derivation and mapping of stream channel networks. This study examines how well commonly applied DEM‐based channel derivation methods at different spatial resolutions can represent the channel network for a glaciated Rocky Mountain headwater catchment. The specific objectives are to (1) examine how differences in gridded DEM resolution affect spatially distributed values of local slope, specific contributing area, and topographic wetness index derived from both eight and infinite directional flow algorithms, (2) map the actual stream channel network to examine the influence of surface variables on channel initiation, and (3) assess accuracy of DEM‐derived networks compared with the field surveyed network. Results show that for the same contributing area threshold, increasing grid cell size leads to increased channelization of modeled networks. A plot of local slope versus contributing area reveals a negative relationship similar to that of prior studies in un‐glaciated areas but with breaks in slope at contributing areas that are too small to represent thresholds for channelization. Field survey results and evaluation of DEM‐derived channel networks suggest that channel network formation is not clearly related to surface topographic variables at Loch Vale. Digitally derived channel networks do not accurately predict low order channel locations, but approximations of the channel network with drainage density and headward extent of channelization similar to the observed network can be derived with both a 1 m and 10 m DEM using a contributing area threshold of approximately 4x10⁴ m². Copyright © 2014 John Wiley & Sons, Ltd.     

Chemistry and mineralogy of high-temperature gas discharges from Colima volcano, Mexico. Implications for magmatic gas–atmosphere interaction

Gases, condensates and silica tube precipitates were collected from 400°C (Z2) and 800°C (Z3) fumaroles at Colima volcano, Mexico, in 1996–1998. Volcanic gases at Colima were very oxidized and contain up to 98% air due to mixing with air inside the dome interior, close to the hot magmatic body. An alkaline trap method was used to collect gas samples, therefore only acidic species were analysed. Colima volcanic gases are water-rich (95–98 mol%) and have typical S/C/Cl/F ratios for a subduction type volcano. δD-values for the high-temperature Z3 fumarolic vapour vary from −26 to −57‰. A negative δD–Cl correlation for the Z3 high-temperature fumarole may result from magma degassing: enrichment in D and decrease in the Cl concentration in condensates are likely a consequence of input of “fresh” batches of magma and an increasing of volcanic activity, respectively.The trace element composition of Colima condensates generally does not differ from that of other volcanoes (e.g. Merapi, Kudryavy) except for some enrichment in V, Cu and Zn. Variations in chemical composition of precipitates along the silica tube from the high-temperature fumarole (Colima 1, fumarole Z3), in contrast to other volcanoes, are characterized by high concentrations of Ca and V, low concentration of Mo and a lack of Cd. Mineralogy of precipitates differs significantly from that described for silica tube experiments at other volcanoes with reduced volcanic gas. Thermochemical modelling was used to explain why very oxidized gas at Colima does not precipitate halite, sylvite, and Mo- and Cd-minerals, but does precipitate V-minerals and native gold, which have not been observed before in mineral precipitates from reduced volcanic gases.     

The volume and shape of the 1991–1993 lava flow field at Mount Etna, Sicily

 The 1991–1993 lava flow is the most voluminous flow erupted at Mount Etna, Sicily, in over 300 years. Estimates of the volume obtained by various methods range from 205×10⁶ m³ (Tanguy 1996) to over 500×10⁶ m³ (Barberi et al. 1993). This paper describes the results of an electronic distance measurement (EDM)-based field survey of the upper surface of the 1991–1993 flow field undertaken in 1995. The results were digitised, interpolated and converted into a digital elevation model and then compared with a pre-eruption digital elevation model, constructed from a 1 : 25 000 contour map of the area, based on 1989 aerial photographs. Our measurements are the most accurate to date and show that the 1991–1993 lava flow occupies a volume of 231±29×10⁶ m³. Received: 20 July 1996 / Accepted: 5 November 1996     

Late Holocene flood magnitudes in the Lower Rhine river valley and upper delta resolved by a two-dimensional hydraulic modelling approach

Palaeoflood hydraulic modelling is essential for quantifying ‘millennial flood’ events not covered in the instrumental record. Palaeoflood modelling research has largely focused on one-dimensional analysis for geomorphologically stable fluvial settings because two-dimensional analysis for dynamic alluvial settings is time consuming and requires a detailed representation of the past landscape. In this study, we make the step to spatially continuous palaeoflood modelling for a large and dynamic lowland area. We applied advanced hydraulic model simulations (1D–2D coupled set-up in HEC-RAS with 950 channel sections and 108 × 10³ floodplain grid cells) to quantify the extent and magnitude of past floods in the Lower Rhine river valley and upper delta. As input, we used a high-resolution terrain reconstruction (palaeo-DEM) of the area in early mediaeval times, complemented with hydraulic roughness values. After conducting a series of model runs with increasing discharge magnitudes at the upstream boundary, we compared the simulated flood water levels with an inventory of exceeded and non-exceeded elevations extracted from various geological, archaeological and historical sources. This comparison demonstrated a Lower Rhine millennial flood magnitude of approximately 14,000 m³/s for the Late Holocene period before late mediaeval times. This value exceeds the largest measured discharges in the instrumental record, but not the design discharges currently accounted for in flood risk management.     

Source location of long period seismicity at Volcàn de Colima,México

This paper presents an analysis of seismicity associated with the volcanic activity of Volcàn de Colima (México) and recorded in the period November 2005–April 2006 during a field survey by the Istituto Nazionale di Geofisica e Vulcanologia (INGV)–Osservatorio Vesuviano, the Observatorio Vulcanologico de Colima of Colima University and the Instituto Andaluz de Geofisica, University of Granada. Three different types of volcanic earthquakes have been identified on the basis of their spectral properties: Type A (0.3–1 Hz), Type B (1–5 Hz) and Type C (3–4 Hz). Results of polarization analysis applied to Type A events show a predominance of radial motion, indicating that the wavefield comprises compressional waves (P) and shear waves polarized in the vertical plane (SV), while the signal always begins with a negative polarity. Type A, B and C earthquakes have been located using both a flat layered model and a 3D model including topography. Hypocentre distributions indicate that the source of Type A signals is very shallow and confined to a small volume lying about 1 km below the crater. In contrast, the source of Type B and C events is significantly deeper, with most hypocentres located in a volume of about 1 km³ centred at 2.5–3 km depth. A cluster analysis based on the cross-correlation among the waveforms of different events recorded at the same station was applied to Type A earthquakes. Only two clusters, which include only a small percentage of events were found, indicating that earthquake families were uncommon during the period of our survey.