Thermal characteristics of coal fires 2: Results of measurements on simulated coal fires

In this paper we present thermal characteristics of coal fires as measured during simulated fires under an experimental setting in Germany in July 2002. It is thus a continuation of the previously published paper “Thermal surface characteristics of coal fire 1: Results of in-situ measurement”, in which we presented temperature measurements of real subsurface coal fires in China [Zhang, J., Kuenzer, C., accepted for publication. Thermal Surface Characteristics of Coal Fires 1: Results of in-situ measurements. Accepted for publication at Journal of Applied Geophysics.]. The focus is on simulated coal fires, which are less complex in nature than fires under natural conditions. In the present study we simulated all the influences usually occurring under natural conditions in a controllable manner (uniform background material of known thermal properties, known ventilation pathways, homogeneous coal substrate), creating two artificial outdoor coal fires under simplified settings. One surface coal fire and one subsurface coal fire were observed over the course of 2 days. The set up of the fires allowed for measurements not always feasible under “real” in-situ conditions: thus compared to the in-situ investigations presented in paper one we could retrieve numerous temperature measurements inside of the fires. Single temperature measurements, diurnal profiles and airborne thermal surveying present the typical temperature patterns of a small surface-and a subsurface fire under undisturbed conditions (easily accessible terrain, 24 hour measurements period, homogeneous materials). We found that the outside air temperature does not influence the fire's surface temperature (up to 900 °C), while fire centre temperatures of up to 1200 °C strongly correlate with surface temperatures of the fire. The fires could heat their surrounding up to a distance of 4.5 m. However, thermal anomalies on the background surface only persist as long as the fire is burning and disappear very fast if the heat source is removed. Furthermore, heat outside of the fires is transported mainly by convection and not by radiation. In spatial thermal line scanner data the diurnal thermal patterns of the coal fire are clearly represented. Our experiments during that data collection also visualize the thermal anomaly differences between covered (underground) and uncovered (surface) coal fires. The latter could not be observed in-situ in a real coal fire area. Sub-surface coal fires express a much weaker signal than open surface fires and contrast only by few degrees against the background. In airborne thermal imaging scanner data the fires are also well represented. Here we could show that the mid-infrared domain (3.8 μm) is more suitable to pick up very hot anomalies, compared to the common thermal (8.8 μm) domain. Our results help to understand coal fires and their thermal patterns as well as the limitations occurring during their analysis. We believe that the results presented here can practicably help for the planning of coal fire thermal mapping campaigns — including remote sensing methods and the thermal data can be included into numerical coal fire modelling as initial or boundary conditions.     

Contribution of the 3-D visualization of acoustic borehole signals (full waveforms) to a quick formation evaluation

Logging and continuous coring are carried out when drilling and looking for materials such as gravels, sand, and clay or in order to evaluate the fracture state of a deep site intended for waste storage. However, in some cases of non-consolidated formations, the results may be disappointing because of the borehole conditions. Full waveforms, as seismic signals, provide information about physical parameters of the grounds crossed by the sonic tool, and this information is almost independent of borehole conditions. Traditional displays in variable area or density show the wave arrival times and the frequencies with depth; for variable density, a color scale permits to see clustered instantaneous phases.In order to determine precisely and simultaneously the three signal parameters (arrival time, frequency, and amplitude) in the depth-propagation time domain, a 3-D visualization software has been developed. The “view” parameters, which give a nice display of the 3-D→2-D projection of the signals in a parallel perspectives relative to depth, are estimated on the monitor screen in an interactive way. A larger size version of the software is available for displaying in detail the acoustic signals for the whole borehole. However, this program needs a large computer, and the maximum size of the drawing depends on the computer memory available for use.The comparison between traditional and 3-D displays shows that without previous preprocessing, the 3-D visualization (1) shows the very small and continuous variations of amplitude (and, therefore, of attenuation) with depth better; and (2) can bring out interferences and “energetic peaks” by simply changing the “view” parameters. As the attenuation of the different waves is directly determined, fresh zones can be distinguished immediately from fracture zones and hard ground from soft ground. The geometry of major fracturing can be deduced directly from graphical representation; i.e. open or closed, and horizontal or oblique fractures. The calculation of their “vertical thickness” is straightfoward. Microfractures, induced by drilling or otherwise, are not detected, but microfractured rocks are detectable.     

Deep root of andesitic volcanoes: New evidence of magma generation at depth in the Benioff zone

The concept of a time-depth correlation between tectonic earthquakes at depth beneath some volcanoes, and their eruptions, developed by the author since 1962, has been confirmed by new observations and successful prediction of renewed volcanic activity in New Zealand.Regular earthquake migrations are observed along the Benioff zone, and volcanic eruptions are found to be related to these seismic migrations beneath the volcanoes, as follows:

Full-size image (2K)

Therefore, in island arcs and continental margins, volcanic activity is the result of two processes occurring beneath the volcanoes: (1) a “tectonic process”, a migration of strain release along the downgoing lithosphere, of which the earthquakes are the manifestation; (2) a “magmatic process”, a relatively fast vertical ascent of magmatic material from the deep root of the volcano, where the observed shocks may be the starting signal from this level.The rate of migration of tectonic earthquakes increases with depth in the upper mantle.An empirical time relationship between the earthquakes occurring at depth beneath a volcano and its eruptions, has been successfully tested for renewed activity at White Island in New Zealand, over the period 1977–1978.     

Thermal state, rheology and seismicity in the pannonian basin, Hungary

On the basis of data on crustal structure and terrestrial heat flow, a 3-D geothermal model for the lithosphere in the Pannonian basin, Hungary, has been calculated. This model, together with information on crustal composition, laboratory data on rock friction, and certain assumptions about fluid conditions and strain-rate levels within the lithosphere, has been used to construct a rheological model of the area.The results obtained show a layered rheological structure where an aseismic part of the crust is “sandwiched” between an upper and a lower seismogenic crustal layers. According to the proposed rheological model, seismic activity in the upper crust may be expected down to depths of 10–12 km, which is confirmed well by the observed depth distribution of seismicity. The model also predicts a lower crustal seismogenic layer down to 20–22 km. Because of infrequent occurrences of deep earthquakes and/or a generally small number of reliable hypocenter depth determinations in the study area, this seismogenic zone is less constrained by observations.The depth of the different rheologic horizons within the crust is governed mainly by thermal conditions. The lower boundary of both seismogenic layers appears isothermal. Brittle-ductile transition in the upper crust coincides with the ˜200 °C isotherm, while in the lower crust it coincides with the ˜ 375 °C isotherm. The lowermost crust and the upper mantle beneath Hungary show ductile behavior, thus the possibility of siesmic activity at these horizons can be excluded.     

Shallow groundwater temperature response to climate change and urbanization

Groundwater temperatures, especially in shallow (quaternary) aquifers respond to ground surface temperatures which in turn depend on climate and land use. Groundwater temperatures, therefore, are modified by climate change and urban development. In northern temperate climate regions seasonal temperature cycles penetrate the ground to depths on the order of 10–15 m. In this paper, we develop and apply analytic heat transfer relationships for 1-D unsteady effective diffusion of heat through an unsaturated zone into a flowing aquifer a short distance below the ground surface. We estimate how changes in land use (urban development) and climate change may affect shallow groundwater temperatures. We consider both long-term trends and seasonal cycles in surface temperature changes. Our analysis indicates that a fully urbanized downtown area at the latitude of Minneapolis/St. Paul is likely to have a groundwater temperature that is nearly 3 °C warmer than an undeveloped agricultural area at the same geographic location. Pavements are the main cause of this change. Data collected by the Minnesota Pollution Control Agency (MPCA) in the St. Cloud, MN area confirm that land use influences groundwater temperatures. Ground surface temperatures are also projected to rise in response to global warming. In the extreme case of a doubling of atmospheric carbon dioxide (2 × CO₂ climate scenario), groundwater temperatures in the Minneapolis/St. Paul metropolitan area could therefore rise by up to 4 °C. Compounding a land use change from “undeveloped” to “fully urbanized” and a 2 × CO₂ climate scenario, groundwater temperatures are projected to rise by about 5 °C at the latitude of Minneapolis/St. Paul.     

Trace-element geochemistry of thermal water percolating through a karstic environment in the region of Saint Ghislain (Belgium)

Five geothermal waters from the Mons area (southern Belgium) have been studied: one natural hot spring at Stambruges, one stagnant warm water from the “inclined tunnels” at Baudour, and three thermal waters from the drillholes at St. Ghislain, Ghlin and Douvrain, originating from the carbonate/anhydrite-bearing Visean strata, at depths of ca. 2600, 1550 and 1300 m, respectively.Multielement chemical analysis of the filtered water and its suspended matter > 0.4 μm) was carried out by instrumental neutron activation.Temperature in depth, calculated using the silica (chalcedony) chemical geothermometer, ranged from 75 to 88°C, in good agreement with experimentally determined values. Na/K and Na/K/Ca geothermometers yieilded erratic results, as expected from the geological environment in the aquifer.From the analytical data it can be calculated that the thermal waters of St. Ghislain, Ghlin and Douvrain are not only saturated with respect to chalcedony, but also to anhydrite, calcite, fluorite, barite, strontianite, and possibly zinc silicate, iron (III) hydroxide or siderite, albite, microcline, gibbsite and kaolinite. They are oversaturated with respect to muscovite. Data are also presented for the other thermal waters, and a cold spring water (Claire Fontaine, Stambruges).The similar trace-element composition of the thermal waters can be explained by percolation of the water in the same distant recharge zone, from where it descends, becomes heated at depth and rises along collapse breccia, and locally (Baudour, Stambruges) along fissures. The uptake of higher amounts of Ca, Mg, Sr and sulfate in St. Ghislain and Ghlin, as compared to Douvrain and Baudour is correlated with the boundary between the “non-dissolved” and “dissolved” evaporitic zones. This boundary is situated between St. Ghislain and Douvrain, and is roughly parallel with the direction of the groundwater flow (WNW).     

Shallow structure beneath the Central Volcanic Complex of Tenerife from new gravity data: Implications for its evolution and recent reactivation

We present a new local Bouguer anomaly map of the Central Volcanic Complex (CVC) of Tenerife, Spain, constructed from the amalgamation of 323 new high precision gravity measurements with existing gravity data from 361 observations. The new anomaly map images the high-density core of the CVC and the pronounced gravity low centred in the Las Cañadas caldera in greater detail than previously available. Mathematical construction of a sub-surface model from the local anomaly data, employing a 3D inversion based on “growing” the sub-surface density distribution via the aggregation of cells, enables mapping of the shallow structure beneath the complex, giving unprecedented insights into the sub-surface architecture. We find the resultant density distribution in agreement with geological and other geophysical data. The modelled sub-surface structure supports a vertical collapse origin of the caldera, and maps the headwall of the ca. 180 ka Icod landslide, which appears to lie buried beneath the Pico Viejo–Pico Teide stratovolcanic complex. The results allow us to put into context the recorded ground deformation and gravity changes at the CVC during its reactivation in spring 2004 in relation to its dominant structural building blocks. For example, the areas undergoing the most significant changes at depth in recent years are underlain by low-density material and are aligned along long-standing structural entities, which have shaped this volcanic ocean island over the past few million years.     

EISCAT observations of unusual flows in the morning sector associated with weak substorm activity

A discussion is given of plasma flows in the dawn and nightside high-latitude ionospheric regions during substorms occurring on a contracted auroral oval, as observed using the EISCAT CP-4-A experiment. Supporting data from the PACE radar, Greenland magnetometer chain, SAMNET magnetometers and geostationary satellites are compared to the EISCAT observations. On 4 October 1989 a weak substorm with initial expansion phase onset signatures at 0030 UT, resulted in the convection reversal boundary observed by EISCAT (at \sim0415 MLT) contracting rapidly poleward, causing a band of elevated ionospheric ion temperatures and a localised plasma density depletion. This polar cap contraction event is shown to be associated with various substorm signatures; Pi2 pulsations at mid-latitudes, magnetic bays in the midnight sector and particle injections at geosynchronous orbit. A similar event was observed on the following day around 0230 UT (\sim0515 MLT) with the unusual and significant difference that two convection reversals were observed, both contracting poleward. We show that this feature is not an ionospheric signature of two active reconnection neutral lines as predicted by the near-Earth neutral model before the plasmoid is “pinched off”, and present two alternative explanations in terms of (1) viscous and lobe circulation cells and (2) polar cap contraction during northward IMF. The voltage associated with the anti-sunward flow between the reversals reaches a maximum of 13 kV during the substorm expansion phase. This suggests it to be associated with the polar cap contraction and caused by the reconnection of open flux in the geomagnetic tail which has mimicked “viscous-like” momentum transfer across the magnetopause.     

The “wave turbopause”

The “wave turbopause” is defined as the mesospheric altitude level where the temperature fluctuation field indicates a substantial increase in wave amplitudes in the vertical direction.The turbopause altitude is analyzed on the basis of four years of SABER data (2002–2005, Version 1.06). Substantial seasonal and latitudinal variations are found, with some interannual variability also present. Seasonal changes are annual at high latitudes, semi-annual at low latitudes, and a mixture of both at middle latitudes. Southern hemisphere data are similar as in the North if shifted by half a year. Latitudinal variations show a minimum in the tropics and two relative maxima at middle latitudes.The “wave turbopause” is found near to zero-wind lines or low-wind zones (zonal wind). It is compared to rocket and other measurements, and interesting similarities are obtained. The wave turbopause can also be found in the HAMMONIA GCM. A preliminary analysis shows results similar to those of the SABER measurements.     

Effects of the Sun on the Earth's environment

Solar disturbances are observed to have significant effects in near-Earth space. Over the past half-century of observation, a relatively clear picture has developed of how and why the typical solar wind — as well as the most extreme solar events — drive geospace responses. It is clear that magnetospheric substorms, geomagnetic storms (both recurrent and aperiodic events), and even certain atmospheric chemical changes have their origins in the solar–terrestrial coupling arena. High-speed solar wind streams and fast coronal mass ejections (CMEs) can often have strong interplanetary shock waves and southward magnetic fields which can initiate strong storm responses. We demonstrate in this review that available modern space-observing platforms and ground facilities allow us to trace drivers from the Sun to the Earth's atmosphere. This allows us to assess quantitatively the energy transport that occurs throughout the Sun–Earth system during both typical and extreme conditions. Hence, we are continuously improving our understanding of “space weather” and its effects on human society.     

On the onset and meridional propagation of the ionospheric F2-region response to geomagnetic storms

The meridional propagation velocities of the ionospheric F2-region response to 268 geomagnetic storms are calculated. Ionospheric vertical sounding data of 1 h time resolution from several stations located in a longitude sector approximately centred along the great circle that contains both the geomagnetic poles and the geographic poles are used.Most meridional propagation velocities from high to low latitudes are less than 600 m/s. The smaller velocities are typical of global neutral meridional wind circulation and the larger are representative of traveling atmospheric disturbances.Simultaneous disturbances at several locations are more frequent during positive phases than during negative phases. Negative phase meridional propagation velocities associated with meridional neutral winds are less frequent in the southern hemisphere when compared with corresponding velocities observed in the northern hemisphere. This may be related to the fact that the distance between the geomagnetic pole and the equator is smaller in the northern hemisphere.Most negative phase onsets are within the 06–10 LT interval. For middle geomagnetic latitudes a “forbidden time interval” between 11 and 14 LT is present. The positive phase onsets show the “dusk effect”.     

Inversion approach for thermal data from a convecting hydrothermal system

Hydrothermal systems are often studied by collecting thermal gradient data and temperature/depth curves. These data contain important information about the flow field, the evolution of the hydrothermal system, and the location and nature of the ultimate heat sources. Thermal data are interpreted by the “forward” method; the thermal field is calculated based on selected initial conditions and boundary conditions such as temperature and permeability distributions. If the calculated thermal field matches the data, the chosen conditions are inferred to be possibly correct. Because many sets of initial conditions may produce similar thermal fields, users of the “forward” method may inadvertently miss the correct set of initial conditions. Analytical methods for “inverting” data also allow the determination of all the possible solutions consistent with the definition of the problem. In this paper we suggest an approach for inverting thermal data from a hydrothermal system, and compare it to the more conventional approach. We illustrate the difference in the methods by comparing their application to the Salton Sea Geothermal Field by Lau (1980a) and Kasameyer, et al. (1984). In this particular example, the inverse method was used to draw conclusions about the age and total rate of fluid flow into the hydrothermal system.     

Earthquake damage detection in the Imperial County Services Building III: Analysis of wave travel times via impulse response functions

The majority of structural health monitoring methods are based on detecting changes in the modal properties, which are global characteristics of the structure, and are not sensitive to local damage. Wave travel times between selected sections of a structure, on the other hand, are local characteristics, and are potentially more sensitive to local damage. In this paper, a structural health monitoring method based on changes in wave travel times is explored using strong motion data from the Imperial Valley Earthquake of 1979 recorded in the former Imperial County Services (ICS) Building, severely damaged by this earthquake. Wave travel times are measured from impulse response functions computed from the recorded horizontal seismic response in three time windows—before, during, and after the largest amplitude response, as determined from previous studies of this building, based on analysis of novelties in the recorded response. The results suggest initial spatial distribution of stiffness consistent with the design characteristics, and reduction of stiffness following the major damage consistent with the spatial distribution of the observed damage. The travel times were also used to estimate the fundamental fixed-base frequency of the structure f₁ (assuming the building deformed as a shear beam), and its changes during this earthquake. These estimates are consistent with previous estimates of the soil–structure system frequency, f_sys, during the earthquakes (f₁<f_sys as expected from soil–structure interaction studies), and with other estimates of frequency (f₁ from ETABS models, and f_sys from ambient vibration tests, and “instantaneous” f₁ from high-frequency pulse propagation).     

Ionospheric measurements during the CRISTA/MAHRSI campaign: their implications and comparison with previous campaigns

The CRISTA/MAHRSI experiment on board a space shuttle was accompanied by a broad campaign of rocket, balloon and ground-based measurements. Supporting lower ionospheric ground-based measurements were run in Europe and Eastern Asia between 1 October–30 November, 1994. Results of comparisons with long ionospheric data series together with short-term comparisons inside the interval October-November, 1994, showed that the upper middle atmosphere (h =80–100 km) at middle latitudes of the Northern Hemisphere in the interval of the CRISTA/MAHRSI experiment (4–12 November, 1994) was very close to its expected climatological state. In other words, the average results of the experiment can be used as climatological data, at least for the given area/altitudes. The role of solar/geomagnetic and “meteorological” control of the lower ionosphere is investigated and compared with the results of MAP/WINE, MAC/SINE and DYANA campaigns. The effects of both solar/geomagnetic and global meteorological factors on the lower ionosphere are found to be weak during autumn 1994 compared to those in MAP/WINE and DYANA winters, and they are even slightly weaker than those in MAP/SINE summer. The comparison of the four campaigns suggests the following overall pattern: in winter the lower ionosphere at northern middle latitudes appears to be fairly well “meteorologically” controlled with a very weak solar influence. In summer, solar influence is somewhat stronger and dominates the weak “meteorological” influence, but the overall solar/meteorological control is weaker than in winter. In autumn we find the weakest overall solar/meteorological control, local effects evidently dominate.     

Thermal structure and energy of the Hakone volcano,Japan

The thermal energy balance and the temperature profile of the Hakone volcano are considered quantitatively. Across the Hakone volcano and its surroundings the heat flow values vary from 10^–1 to 10³ mW/m², due to thermal conduction and mass flow involving volcanic steam and hot spring discharge. An area with extremely low heat flow is observed in the western side of the caldera showing the presence of percolating meteoric water. The hydrothermal activity is intense in the eastern half of the caldera.The total heat discharge from the high temperature zone (discharge area) of the Hakone volcano amounts to 11.0×10⁷ W. The magmatic steam energy discharge is 95.0×10⁶ W. The thermal energy by redistribution of the terrestrial heat flow by the lateral deep ground water flow is calculated to be 9.00×10⁶ W. For the model having the vertical vent in the volcano's central part up to 1 km depth below the ground surface from a magma reservoir the computed temperature distribution is consistent with the observed values. The depth of the magma reservoir is 7 km below the ground surface and the diameter is 5 km.     

Estimation of Nonlinear Time-dependent Soil Behavior in Strong Ground Motion Based on Vertical Array Data

To improve our understanding of nonlinear elastic properties of soils, a method is proposed of estimation of stress-strain relations of soils in situ in strong ground motion based on vertical array data. Strong motion records provided by seismic vertical arrays allow estimation of nonlinear stress-strain relations in soil layers at different depths, from the surface down to the location of the deepest device. As an example, records obtained during the main shock of the 1995 Kobe earthquake at Port-Island, SGK, and TKS sites were used to estimate the stress-strain relations in the soil profiles. For different layers, different types of nonlinear stress-strain relations were selected, according to the profiling data. To account for temporal changes in the soil behavior, consecutive parts of records were examined, and for successive time intervals, the relations were found showing the best-fit approximation to the observed data. At Port Island and SGK sites, where the strongest accelerations were recorded, the obtained stress-strain relations showed systematic changes in the upper layers (0–14 m), such as, a progressive reduction of the slopes of the stress-strain curves due to liquefaction at Port Island and reduction and recovery of the slopes at SGK and TKS sites. At the three sites, the stress-strain relations remained stable in layers below 11–14 m. Thus, the proposed approach gives us a representation of the soil behavior in layers at different depths in strong ground motion; it allows calculation of the propagation of arbitrary seismic signals in the studied profiles and estimation of nonlinear components in the ground response by the nonlinear system identification technique. The method can also be applied to evaluate the ground response at sites where profiling data are available and an imposed motion can be estimated.     

Artificial spacecraft in hybrid simulations of the quasi-parallel Earth’s bow shock: analysis of time series versus spatial profiles and a separation strategy for Cluster

We construct artificial “software” spacecraft consisting of magnetometers and 3D thermal and energetic ion detectors. Four such spacecraft are “flown” through a 1D simulation of a quasi-parallel shock. We analyze the resulting time series from the spacecraft, and then use the more complete simulational information to evaluate our interpretations based on the limited times series information. The separation strategy used, with two closely spaced spacecraft pairs separated by a large distance, was helpful in the interpretation, since a variety of important processes operate over several different scale lengths. This work highlights the ability to draw inferences about spatially and temporally varying phenomena based on multiple-spacecraft time series data, and suggests that many spacecraft configurations which bear little resemblance to the classic Cluster tetrahedron may be necessary when multiple scale lengths are present.     

Geoelectrical methods applied to structures of arbitrary shapes

Complicated geological structures, different from the plane-parallel stratification, often occur in the underground and present problems for the geoelectrical prospecting. At least from a theoretical point of view, modern exploration procedures for dealing with this type of situations may be planned. These procedures are based on a large number of data distributed on the area to be explored, and on rigorous inversion methods, carried out by means of powerful computers.Such a procedure leads to models of the subsoil, on grounds of the apparent resistivity values measured on the surface. This is an inverse problem which implies the use, in opposite direction, of a direct computing method, i.e. a procedure allowing the computation of apparent resistivity values from a structural model. One of these, the “surface charge” algorithm is described in this paper.In practice, implementation of such ideal procedures often encounters difficulties of technical and also of economic type, at least for the 3-D problems. The consequence is that, in many cases, one is obliged to use traditional procedures, as, for instance, the use of vertical soundings distributed on the area to be surveyed.The application of traditional procedures to complicated problems requires some changes. Firstly, the popular Schlumberger array must be avoided, as it may result into insufficient information. Moreover, the array is not suitable for deep explorations because of the excessive cable length. More convenient is the dipole-dipole array, which, however, presents its own drawbacks, as, for instance, a higher sensitivity towards lateral variations. This last incovenience may be avoided by means of the transformation of the field dipole-dipole into half-Schlumberger diagrams; but it has been shown that, over arbitrary structures, this is possible only if “polar dipole-dipole” arrays are used in the field.In the final part of the paper some field examples of “continuous” polar dipole-dipole soundings carried out for deep exploration are shown.     

Integrated satellite observations of the 2001 eruption of Mt. Cleveland, Alaska

Satellite data were the primary source of information for the eruption of Mt. Cleveland, Alaska on 19 February, and 11 and 19 March 2001. Multiple data sets were used pre-, syn- and post-eruption to mitigate the hazard and determine an eruption chronology. The 19 February eruption was the largest of the three, resulting in a volcanic cloud that formed an arc over 1000 km long, moved to the NE across Alaska and was tracked using satellite data over more than a 50-h period. The volcanic cloud was “concurrently” detected on the GOES, AVHRR and MODIS data at various times and their respective signals compared. All three sensors detected a cloud that had a very similar shape and position but there were differences in their areal extent and internal structural detail. GOES data showed the largest volcanic cloud in terms of area, probably due to its oblique geometry. MODIS bands 31 and 32, which are comparable to GOES and AVHRR thermal infrared wavelengths, were the least effective single channels at detecting the volcanic cloud of those investigated (MODIS bands 28, 29, 31 and 32). MODIS bands 28 and 29 detected the largest volcanic clouds that could easily be distinguished from weather clouds. Of the split-window data, MODIS bands 29 minus band 32 detected the largest cloud, but the band 31 minus band 32 data showed the volcanic cloud with the most internal structural detail. The Puff tracking model accurately tracked the movement, and predicted the extent and shape of this complex cloud even into areas beyond satellite detection. Numerous thermal anomalies were also observed during the eruption on the twice-daily AVHRR data and the high spatial-resolution Landsat data. The high-resolution Radarsat data showed that the AVHRR thermal anomalies were due to lava and debris flow features and a newly formed fan along the west coast of the island. Field observations and images from a hand-held Forward Looking Infrared Radiometer (FLIR) showed that the flow features were ′a′a lava, debris flows and a warm debris fan along the west coast. Real-time satellite data were the primary tool used to monitor the eruption, track changes and to mitigate hazards. High-resolution data, even though coverage is infrequent, were critical in helping to identify volcanic processes and to compile an eruption chronology.     

Evolution of the latera geothermal system II: metamorphic, hydrothermal mineral assemblages and fluid chemistry

The Latera field (Vulsini volcanic complex, Latium, Italy) is one of the geothermal areas of the peri-Tyrrhenian belt along which a regional, high thermal anomaly has been detected. So far nine deep wells have been drilled within the Latera caldera and four of them have been productive. The geothermal reservoir is located within the fractured carbonatic rocks of the Tuscan nappe; the overlying volcanic units, sealed by hydrothermal minerals (mainly calcite and anhydrite), act as an impervious cover.The fluid produced by the wells comes from a deep aquifer (about 1000–1500 m depth) which at present is not connected with the shallow aquifer in the volcanoclastic units. Fluid temperatures range between 200 and 230°C; in-hole temperatures as high as 343°C at 2775 m depth have been measured in dry wells.The study of the newly formed mineral assemblages from both volcanic and sedimentary units as sampled from the geothermal wells can be used to reconstruct the thermal evolution of the geothermal field. The intrusion of a syenitic melt, up to a depth of about 2000 m, dated 0.86 Ma, represents the major thermal event for the units in the area and is assumed to represent the first step in the geothermal evolution of the Latera system.The above mentioned newly formed mineral assemblages can be divided into three groups: (a) “contact-metasomatic”: calcite, anhydrite, diopsidic pyroxene, grossularitic garnet, phlogopite, wollastonite or monticellite; (b) “high-temperature hydrothermal”: calcite, anhydrite, K-feldspar, vesuvianite, melanitic garnet, tourmaline, amphibole, epidote, sulphides; (c) “low-temperature hydrothermal”: calcite, anhydrite, K-feldspar, clay minerals, sulphides. Group (a) minerals are now relics. Part of (b) and all of (c) group are still in equilibrium with the existing conditions in different parts of the geothermal system.Thermodynamic calculations on the observed mineral assemblages permitted estimates of the P, T conditions and gas fugacities.