Eruptive history of the Colima volcanic complex (Mexico)

The evolution of the Colima volcanic complex can be divided into successive periods characterized by different dynamic and magmatic processes: emission of andesitic to dacitic lava flows, acid-ash and pumice-flow deposits, fallback nuées ardentes leading to pyroclastic flows with heterogeneous magma, plinian air-fall deposits, scoriae cones of alkaline and calc-alkaline nature. Four caldera-forming events, resulting either from major ignimbrite outbursts or Mount St. Helens-type eruptions, separate the main stages of development of the complex from the building of an ancient shield volcano (25 × 30 km wide) up to two summit cones, Nevado and Fuego.The oldest caldera, C1 (7–8 km wide), related to the pouring out of dacitic ash flows, marks the transition between two periods of activity in the primitive edifice called Nevado I: the first one, which is at least 0.6 m.y. old, was mainly andesitic and effusive, whereas the second one was characterized by extrusion of domes and related pyroclastic products. A small summit caldera, C2 (3–3.5 km wide), ended the evolution of Nevado I.Two modern volcanoes then began to grow. The building of the Nevado II started about 200,000 y. ago. It settled into the C2 caldera and partially overflowed it. The other volcano, here called Paleofuego, was progressively built on the southern side of the former Nevado I. Some of its flows are 50,000 y. old, but the age of its first outbursts is not known. However, it is younger than Nevado II. These two modern volcanoes had similar evolutions. Each of them was affected by a huge Mount St. Helens-type (or Bezymianny-type) event, 10,000 y. ago for the Paleofuego, and hardly older for the Nevado II. The landslides were responsible for two horseshoe-shaped avalanche calderas, C3 (Nevado) and C4 (Paleofuego), each 4–5 km wide, opening towards the east and the south. In both cases, the activity following these events was highly explosive and produced thick air-fall deposits around the summit craters.The Nevado III, formed by thick andesitic flows, is located close to the southwestern rim of the C3 caldera. It was a small and short-lived cone. Volcan de Fuego, located at the center of the C4 caldera, is nearly 1500 m high. Its activity is characterized by an alternation of long stages of growth by flows and short destructive episodes related to violent outbursts producing pyroclastic flows with heterogeneous magma and plinian air falls.The evolution of the primitive volcano followed a similar pattern leading to formation of C1 and then C2. The analogy between the evolutions of the two modern volcanoes (Nevado II–III; Paleofuego-Fuego) is described. Their vicinity and their contemporaneous growth pose the problem of the existence of a single reservoir, or two independent magmatic chambers, after the evolution of a common structure represented by the primitive volcano.     

Morphological analysis of active Mount Nemrut stratovolcano, eastern Turkey: evidences and possible impact areas of future eruption

Mount Nemrut, an active stratovolcano in eastern Turkey, is a great danger for its vicinity. The volcano possesses a summit caldera which cuts the volcano into two stages, i.e. pre- and post-caldera. Wisps of smoke and hot springs are to be found within the caldera. Although the last recorded volcanic activity is known to have been in 1441, we consider here that the last eruption of Nemrut occurred more recently, probably just before 1597. The present active tectonic regime, historical eruptions, occurrence of mantle-derived magmatic gases and the fumarole and hot spring activities on the caldera floor make Nemrut Volcano a real danger for its vicinity. According to the volcanological past of Nemrut, the styles of expected eruptions are well-focused on two types: (1) occurrence of water within the caldera leads to phreatomagmatic (highly energetic) eruptions, subsequently followed by lava extrusions, and (2) effusions–extrusions (non-explosive or weakly energetic eruptions) on the flanks from fissures. To predict the impact area of future eruptions, a series of morphological analyses based on field observations, Digital Elevation Model and satellite images were realized. Twenty-two valleys (main transport pathways) were classified according to their importance, and the physical parameters related to the valleys were determined. The slope values in each point of the flanks and the Heim parameters H/L were calculated. In the light of morphological analysis the possible impact areas around the volcano and danger zones were proposed. The possible transport pathways of the products of expected volcanic events are unified in three main directions: Bitlis, Guroymak, Tatvan and Ahlat cities, the about 135 000 inhabitants of which could be threatened by future eruptions of this poorly known and unsurveyed volcano.     

Stable Eocene Magnetization Carried by Magnetite and Iron Sulphides in Marine Marls (Pamplona-Arguis Formation,Southern Pyrenees,Northern Spain)

In order to establish the magnetic carriers and assess the reliability of previous paleomagnetic results obtained for Eocene marine marls from the south Pyrenean basin, we carried out a combined paleo- and rock-magnetic study of the Pamplona-Arguis Formation, which crops out in the western sector of the southern Pyrenees (N Spain). The unblocking temperatures suggest that the characteristic remanent magnetization (ChRM) is carried by magnetite and iron sulphides. The ChRM has both normal and reversed polarities regardless of whether it resides in magnetite or iron sulphides, and represents a primary Eocene magnetization acquired before folding. Rock magnetic results confirm the presence of magnetite and smaller amounts of magnetic iron sulphides, most likely pyrrhotite, in all the studied samples. Framboidal pyrite is ubiquitous in the marls and suggests that iron sulphides formed during early diagenesis under sulphate-reducing conditions. ChRM directions carried by magnetic iron sulphides are consistent with those recorded by magnetite. These observations suggest that magnetic iron sulphides carry a chemical remanent magnetization that coexists with a remanence residing in detrital magnetite. We suggest that the south Pyrenean Eocene marls are suitable for magnetostratigraphic and tectonic purposes but not for studies of polarity transitions, secular variations and geomagnetic excursions, because it is difficult to test for short time differences in remanence lock-in time for the two minerals. The presence of iron sulphide minerals contributing to the primary magnetization in Eocene marine marls reinforces the idea that these minerals can persist over long periods of time in the geological record.     

Subsurface DC resistivity mapping: approximate 1-D interpretation

Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.     

Characterization of juvenile pyroclasts from the Kos Plateau Tuff (Aegean Arc): insights into the eruptive dynamics of a large rhyolitic eruption

Silicic pumices formed during explosive volcanic eruptions are faithful recorders of the state of the magma in the conduit, close to or at the fragmentation level. We have characterized four types of pumices from the non-welded rhyolitic Kos Plateau Tuff, which erupted 161,000 years ago in the East Aegean Arc, Greece. The dominant type of pumice (>90 vol.%) shows highly elongated tubular vesicles. These tube pumices occur throughout the eruption. Less common pumice types include: (1) “frothy” pumice (highly porous with large, sub-rounded vesicles), which form 5–10 vol.% of the coarsest pyroclastic flow deposits, (2) dominantly “microvesicular” and systematically crystal-poor pumices, which are found in early erupted, fine-grained pyroclastic flow units, and are characterized by many small (<50 μm in diameter) vesicles and few mm-sized, irregular voids, (3) grey or banded pumices, indicating the interaction between the rhyolite and a more mafic magma, which are found throughout the eruption sequence and display highly irregular bubble shapes. Except for the grey-banded pumices, all three other types are compositionally identical and were generated synchronously as they are found in the same pyroclastic units. They, therefore, record different conditions in the volcanic conduit leading to variable bubble nucleation, growth and coalescence. A total of 74 pumice samples have been characterized using thin section observation, SEM imagery, porosimetry, and permeametry. We show that the four pumice types have distinct total and connected porosity, tortuosity and permeability. Grey-banded pumices show large variations in petrophysical characteristics as a response to mingling of two different magmas. The microvesicular, crystal-poor, pumices have a bimodal bubble size distribution, interpreted as reflecting an early heterogeneous bubble nucleation event followed by homogeneous bubble nucleation close to fragmentation. Finally, the significant differences in porosity, tortuosity and permeability in compositionally identical tube and frothy pumices are the result of variable shear rates in different parts of the conduit. Differential shear rates may be the result of either: (1) pure shear, inducing a vertical progression from frothy to tube and implying a relatively thick fragmentation zone to produce both types of pumices at the same time or (2) localized simple shear, inducing strongly tubular vesicles along the wall and near-spherical bubbles in the centre of the conduit and not necessarily requiring a thick fragmentation zone.     

A new method for sampling and analyzing volcanic sublimates — application to Merapi volcano, Java

A new method for the sampling of sublimates from high-temperature volcanic gases has been used at Merapi volcano, Java, in 1978. The sublimates were collected on the inner walls of silica tubes introduced into fumarolic vents. Volcanic gases were allowed to move freely through the tubes and as they cooled, a fraction of the volatile components condensed on the inner walls of the tubes along the temperature gradient. The sublimates were then analyzed by a combination of light microscopy, scanning electron microscopy, electron microprobe and X-ray diffraction.Six successive zones of different compositions and mineralogical associations have been identified along the covered range of temperatures (900° to around 400°C). From the high to the low temperatures, these zones are composed of: (1) cristobalite, magnetite, hercynite; (2) molybdenite; (3) acmite; (4) halite, sylvite; (5) sphalerite, pyrite; and (6) galena. Equilibrium calculations show that these crystalline phases are stable for pS₂, pC₁, and pO₂, values typical of magma-buffered gases that have not been contaminated by atmospheric oxygen.The deposits observed in the tubes may be useful in aiding the understanding of the mechanisms acting during the cooling of the gaseous phase on its way to the surface and before its emission into the atmosphere.     

Turbulent heat fluxes by profile and inertial dissipation methods: analysis of the atmospheric surface layer from shipboard measurements during the SOFIA/ASTEX and SEMAPHORE experiments

Heat flux estimates obtained using the inertial dissipation method, and the profile method applied to radiosonde soundings, are assessed with emphasis on the parameterization of the roughness lengths for temperature and specific humidity. Results from the inertial dissipation method show a decrease of the temperature and humidity roughness lengths for increasing neutral wind speed, in agreement with previous studies. The sensible heat flux estimates were obtained using the temperature estimated from the speed of sound determined by a sonic anemometer. This method seems very attractive for estimating heat fluxes over the ocean. However allowance must be made in the inertial dissipation method for non-neutral stratification. The SOFIA/ASTEX and SEMAPHORE results show that, in unstable stratification, a term due to the transport terms in the turbulent kinetic energy budget, has to be included in order to determine the friction velocity with better accuracy. Using the profile method with radiosonde data, the roughness length values showed large scatter. A reliable estimate of the temperature roughness length could not be obtained. The humidity roughness length values were compatible with those found using the inertial dissipation method.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

Measuring and interpreting heat fluxes from shallow volcanic bodies using vertical temperature profiles: a preliminary test

 To test the potential of heat flux prospecting in active volcanic areas using shallow temperature data taken along vertical profiles, we carried out two thermal profile surveys, one not far from Yasur cone on Tanna Island, and another inside the caldera of Ambrym (New Hebrides arc, southwestern Pacific). The basic steady heat flux of internal volcanic origin was determined, taking into account both conductive and convective heat transfers. At both locations there exists, over small distances, significant differences in the heat flux. These differences correspond to shallow sources of heat. The use of a network of vertical profiles allowed: (a) heat flux mapping; (b) location of shallow volcanic heat sources; and (c) observation of the detailed structure of the heat release at quiescent but active volcanoes. Received: 18 July 1997 / Accepted: 13 May 1998     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.