Volcanic history of El Chichón Volcano (Chiapas, Mexico) during the Holocene, and its impact on human activity

2 abundance and depletion in MgO, CaO, TiO₂, as well as trace and rare earth elements. This suggests segregation of olivine and orthopyroxene from the melt. Since human settlements in southeast Mexico and Central America can be traced as far back as approximately 2500 years BP, most of these events probably affected human activity. In fact, there are reports of pottery shards and other artifacts in deposits from the eruption of 1250 BP. Pottery fragments in deposits of an eruption that took place 2500 BP are also reported in this paper. Thus, the impact of the volcano on human activities has been frequent, with most of the repose intervals lasting between 100 to 600 years. The impact of the eruptions was probably of greater than local extent, because airfall tephra could reach distant sites and possibly even affect weather. The eruptive history of El Chichón also offers clues in the investigation of the Maya civilization. Several researchers have considered the volcano as an important factor in the answer to some intriguing questions such as the extensive use of volcanic ash in Late Classic Maya ceramics or, of greater importance, the causes of the collapse of the Classic Maya civilization. Received : 4 May 1998 / Accepted: 1 November 1999     

Eruption of kimberlite magmas: physical volcanology, geomorphology and age of the youngest kimberlitic volcanoes known on earth (the Upper Pleistocene/Holocene Igwisi Hills volcanoes, Tanzania)

The Igwisi Hills volcanoes (IHV), Tanzania, are unique and important in preserving extra-crater lavas and pyroclastic edifices. They provide critical insights into the eruptive behaviour of kimberlite magmas that are not available at other known kimberlite volcanoes. Cosmogenic ³He dating of olivine crystals from IHV lavas and palaeomagnetic analyses indicates that they are Upper Pleistocene to Holocene in age. This makes them the youngest known kimberlite bodies on Earth by >30?Ma and may indicate a new phase of kimberlite volcanism on the Tanzania craton. Geological mapping, Global Positioning System surveying and field investigations reveal that each volcano comprises partially eroded pyroclastic edifices, craters and lavas. The volcanoes stand <40?m above the surrounding ground and are comparable in size to small monogenetic basaltic volcanoes. Pyroclastic cones consist of diffusely layered pyroclastic fall deposits comprising scoriaceous, pelletal and dense juvenile pyroclasts. Pyroclasts are similar to those documented in many ancient kimberlite pipes, indicating overlap in magma fragmentation dynamics between the Igwisi eruptions and other kimberlite eruptions. Characteristics of the pyroclastic cone deposits, including an absence of ballistic clasts and dominantly poorly vesicular scoria lapillistones and lapilli tuffs, indicate relatively weak explosive activity. Lava flow features indicate unexpectedly high viscosities (estimated at >10² to 10⁶?Pa?s) for kimberlite, attributed to degassing and in-vent cooling. Each volcano is inferred to be the result of a small-volume, short-lived (days to weeks) monogenetic eruption. The eruptive processes of each Igwisi volcano were broadly similar and developed through three phases: (1) fallout of lithic-bearing pyroclastic rocks during explosive excavation of craters and conduits; (2) fallout of juvenile lapilli from unsteady eruption columns and the construction of pyroclastic edifices around the vent; and (3) effusion of degassed viscous magma as lava flows. These processes are similar to those observed for other small-volume monogenetic eruptions (e.g. of basaltic magma).     

Modelling the effect of time‐dependent river dune evolution on bed roughness and stage

This paper presents an approach to incorporate time‐dependent dune evolution in the determination of bed roughness coefficients applied in hydraulic models. Dune roughness is calculated by using the process‐based dune evolution model of Paarlberg et al. ( 2009 ) and the empirical dune roughness predictor of Van Rijn ( 1984 ). The approach is illustrated by applying it to a river of simple geometry in the 1‐D hydraulic model SOBEK for two different flood wave shapes. Calculated dune heights clearly show a dependency on rate of change in discharge with time: dunes grow to larger heights for a flood wave with a smaller rate of change. Bed roughness coefficients computed using the new approach can be up to 10% higher than roughness coefficients based on calibration, with the largest differences at low flows. As a result of this larger bed roughness, computed water depths can be up to 15% larger at low flow. The new approach helps to reduce uncertainties in bed roughness coefficients of flow models, especially for river systems with strong variations in discharge with time. Copyright © 2010 John Wiley & Sons, Ltd.     

Late-Pleistocene flank collapse triggered by dome growth at Tacaná volcano, México-Guatemala, and its relationship to the regional stress regime

During late Pleistocene time, the extrusion of an andesitic dome at the summit of Tacaná volcano caused the collapse of its northwestern flank. The stratocone collapse was nearly parallel to the σ _min stress direction suggesting that failure was controlled by the regional stress field. The event produced a debris avalanche that was channelized in the San Rafael River and moved 8 km downstream. The deposit covered a minimum area of 4 km², had a volume of 0.8 ± 0.5 km³, with an H/L (vertical drop to horizontal transport distance ratio) of ~0.35, defining a degree of mobility that is atypical for volcanic debris avalanches. The flank failure undermined the summit dome leading to its collapse and the generation of a series of block-and-ash flows that were emplaced in quick succession and covered the avalanche surface. The collapse event left a 600-m-wide summit amphitheatre with a 30-degree opening to the northwest, and >200 m thick debris that blocked the San Rafael River. Remobilization of this material produced debris flows that eroded the primary deposits and cascaded into the Coatán River. After the collapse, the activity of Tacaná continued with the emission of the Agua Zarca lava flow dated at 10 ± 6 ka (⁴⁰Ar/³⁹Ar), and pyroclastic surges dated at 10,610 + 330/−315 yr BP (¹⁴C), which provide a minimum age for the collapse event. During the Holocene, Tacaná has been very active producing explosive and effusive eruptions that ended with the extrusion of two summit domes that today occupy the amphitheatre. The 1950 and 1986 phreatic outbursts occurred along the Pleistocene collapse scar. Currently ~300,000 inhabitants live within a 35 km radius of Tacaná, and could conceivably be impacted by future events of similar magnitude.     

Efficient Methods to Simulate Planar Free Surface in the 3D 4th-Order Staggered-Grid Finite-Difference Schemes

We numerically tested accuracy of two formulations of Levander's (1988) stress-imaging technique for simulating a planar free surface in the 4^th-order staggered-grid finite-difference schemes. We have found that both formulations (one with normal stress-tensor components at the surface, the other with shear stress-tensor components at the surface) require at least 10 grid spacings per minimum wavelength ( _min÷h = 10) if Rayleigh waves are to be propagated without significant grid dispersion in the range of epicentral distances up to 15 _dom ^S.Because interior 4^th-order staggered-grid schemes usually do not require more than 6 grid spacings per minimum wavelength, in the considered range of epicentral distances, it was desirable to find alternative techniques to simulate a planar free surface, which would not require denser spatial sampling than _min÷h = 6. Therefore, we have developed and tested new techniques: 1. Combination of the stress imaging (with the shear stress-tensor components at the surface) with Rodrigues' (1993) vertically refined grid near the free surface. 2. Application of the adjusted finite-difference approximations to the z-derivatives at the grid points at and below the surface that uses no virtual values above the surface and no stress imaging. The normal stress-tensor components are at the surface in one formulation, while the shear stress-tensor components are at the surface in the other formulation.The three developed formulations give for the spatial sampling _min÷h = 6 results very close to those obtained by the discrete-wavenumber method. Because, however, the technique with the vertically refined grid near the free surface requires 3 times smaller time step (due to the refined grid), the technique with adjusted finite-difference approximations is the most accurate and efficient technique from the examined formulations in the homogeneous halfspace.     

Temporal changes in water quality at a childhood leukemia cluster

Since 1997, 15 cases of acute lymphocytic leukemia and one case of acute myelocytic leukemia have been diagnosed in children and teenagers who live, or have lived, in an area centered on the town of Fallon, Nevada. The expected rate for the population is about one case every five years. In 2001, 99 domestic and municipal wells and one industrial well were sampled in the Fallon area. Twenty-nine of these wells had been sampled previously in 1989. Statistical comparison of concentrations of major ions and trace elements in those 29 wells between 1989 and 2001 using the nonparametric Wilcoxon signed-rank test indicate water quality did not substantially change over that period; however, short-term changes may have occurred that were not detected. Volatile organic compounds were seldom detected in ground water samples and those that are regulated were consistently found at concentrations less than the maximum contaminant level (MCL). The MCL for gross-alpha radioactivity and arsenic, radon, and uranium concentrations were commonly exceeded, and sometimes were greatly exceeded. Statistical comparisons using the nonparametric Wilcoxon rank-sum test indicate gross-alpha and -beta radioactivity, arsenic, uranium, and radon concentrations in wells used by families having a child with leukemia did not statistically differ from the remainder of the domestic wells sampled during this investigation. Isotopic measurements indicate the uranium was natural and not the result of a 1963 underground nuclear bomb test near Fallon. In arid and semiarid areas where trace-element concentrations can greatly exceed the MCL, household reverse-osmosis units may not reduce their concentrations to safe levels. In parts of the world where radon concentrations are high, water consumed first thing in the morning may be appreciably more radioactive than water consumed a few minutes later after the pressure tank has been emptied because secular equilibrium between radon and its immediate daughter progeny is attained in pressure tanks overnight.     

Pyroclastic flow deposits of the 1991 eruption of Volcán de Colima,Mexico

The April 16, 1991, eruption of Volcán de Colima represents a classical example of partial dome collapse with the generation of progressively longer-runout, Merapi-type pyroclastic flows that traveled up to 4 km along the El Cordoban gullies (East, Central and West). The flows filled the gullies with block-and-ash flow deposits up to 10 m thick, of which, after 7 years of erosion, only remnants remained in the El Cordoban West and East gullies. The El Cordoban Central gully, however, provided a well-preserved and incised longitudinal section of the 1991 deposits. The deposits were emplaced as proximal and distal facies, separated by a change in slope angle from >30° to <20°. The proximal facies consists of massive, clast-supported flow units (up to 1 m thick) with andesite blocks locally supported by a matrix of coarse ash and devoid of segregation structures or grading. The distal facies consists of a massive, matrix-supported deposit up to 8 m thick, which contains dispersed andesite blocks in a fine ash matrix. In the distal facies, a train of blocks marks flow-unit upper boundaries and, although sorting is poor, some grading is present. Thin, finely stratified, or dune-bedded layers of fine ash material are locally present above or below units of both facies. Sedimentologic parameters show that the size or fraction of large pyroclasts (larger than –1 ) decreases from proximal to distal facies, as the percentage of matrix (0 to 4 ) increases, especially immediately beyond the break in slope. We propose that the propagation of the Colima pyroclastic flows is critically dependent on local slope angle, the presence of erodible slope debris, and the decrease in grain size with distance from the vent. The progressive fining is probably caused by some combination of erosion, clast breakup and deposition of larger pyroclasts, and is itself influenced by the slope angle. In the proximal region, the flows moved as granular avalanches, in which interacting grains ground each other and erosion occurred to produce an overriding dilute ash cloud. The maximum runout distance of the avalanches was controlled by the angle of repose of the material, and the volume and grain size of source and eroded material. Because the slope angle is close to the repose angle for this debris, granular avalanches were not able to propagate far beyond the change in slope. If, however, an avalanche had enough mass in finer grain size fractions, at least part of the flow continued beyond the break in slope and across the volcano apron, propagating in a turbulent state and depositing surge layers, or in an otherwise settling-modified state and depositing block-and-ash flow layers.Editorial responsibility: T Druitt     

Variability of dense water formation in the Ross Sea

This paper presents results from a model study of the interannual variability of high salinity shelf water (HSSW) properties in the Ross Sea. Salinity and potential temperature of HSSW formed in the western Ross Sea show oscillatory behaviour at periods of 5–6 and 9 years superimposed on long-term fluctuations. While the shorter oscillations are induced by wind variability, variability on the scale of decades appears to be related to air temperature fluctuations. At least part of the strong decrease of HSSW salinities deduced from observations for the period 1963–2000 is shown to be an aliasing artefact due to an undersampling of the periodic signal. While sea ice formation is responsible for the yearly salinity increase that triggers the formation of HSSW, interannual variability of net freezing rates hardly affects changes in the properties of the resulting water mass. Instead, results from model experiments indicate that the interannual variability of dense water characteristics is predominantly controlled by variations in the shelf inflow through a sub-surface salinity and a deep temperature signal. The origin of the variability of inflow characteristics to the Ross Sea continental shelf can be traced into the Amundsen and Bellingshausen Seas. The temperature anomalies are induced at the continental shelf break in the western Bellingshausen Sea by fluctuations of the meridional transport of circumpolar deep water with the eastern cell of the Ross Gyre. In the Amundsen Sea, upwelling due to a persistently cyclonic wind field carries the signal into the surface mixed layer, leading to fluctuations of the vertical heat flux, anomalies of brine release near the sea ice edge, and consequently to a sub-surface salinity anomaly. With the westward flowing coastal current, both the sub-surface salinity and deep temperature signals are advected onto the Ross Sea continental shelf. Convection carries the signal of salinity variability into the deep ocean, where it interacts with modified circumpolar deep water upwelled onto the continental shelf as the second source water mass of HSSW. Sea ice formation on the Ross Sea continental shelf thus drives the vertical propagation of the signal rather than determining the signal itself.     

The September 8–9, 1998 Rain-Triggered Flood Events at Motozintla, Chiapas, Mexico

In September 1998 tropical storm “Earl” swept southern Mexico, producing intense rainfall in the states of Oaxaca and Chiapas. Among the most devastated cities was Motozintla, located in the drainage basin of the Allende, La Mina and Xelajú Grande Rivers. The rainfall from the tropical storm totaled 175 mm on September 8 and 130 mm on September 9, duplicating in two days the average monthly precipitation in the region. Numerous landslides occurred in the vicinity of Motozintla, depositing large volumes of material into the Xelajú Grande stream. Much of this sediment was subsequently remobilized, yielding debris flows, hyperconcentrated flows, and sediment-laden flows that inundated most sections of Motozintla city. The flows covered an approximate area of 3.15 km² with a minimum volume of 4.4 × 10⁶ m³ of sediment. Communication of Motozintla with the rest of the Chiapas State was interrupted for about a month, as was the supply of potable water, food, electricity, and fuel. The geologic record around Motozintla indicates that the Xelajú Grande River has been a pathway for similar large floods during the last 6000 years. The oldest deposit yielded a radiocarbon age of 5320 ± 100 ¹⁴C years. B.P. At least two historic floods have occurred during the last 100 years, a time period defined by a stratigraphically distinct tephra of 1902. Frequency analysis of the historical record of daily rainfall in the Motozintla area suggests that events like that of September, 1998, have a recurrence interval of about 25 years. After the catastrophic flows of 1998, the mitigation measures by Municipal Authorities were made without regard to geological and environmental factors, or to taking into consideration the flow magnitude and appropriate hazard-mitigation techniques, with the result that Motozintla remains at serious risk for future floods. Unfortunately, prior to the publication of this study, in early October 2005, Motozintla was seriously damaged again by intense rain provoked by Hurricane Stan.