R S A V: Emote Ensing of Ctive Olcanoes
R S A V: Emote Ensing of Ctive Olcanoes
R S A V: Emote Ensing of Ctive Olcanoes
28:81106
Copyright q 2000 by Annual Reviews. All rights reserved
deceased
INTRODUCTION
In 1972, the United States sent the Earth Resources Technology Satellite (ERTS)
into orbit. Later renamed Landsat 1, this was the first of a series of polar orbiting
satellites to provide routine, high-resolution, multispectral images of the Earths
surface for civilian applications. Its launch could be said to initiate the modern
era of remote sensing of active volcanoes. A further major step was the October
1974 launch of the first Geostationary Operational Environmental Satellite
(GOES), which provided continuous monitoring of North American weather.
Since then, a plethora of polar orbiting and geostationary satellites has been
launched, carrying sensors working in all parts of the electromagnetic spectrum
from ultraviolet to microwave. Some of these satellites were spectacularly shortlived, whereas others worked for many years; however, none had volcanology as
a prime objective. Meteorological observations were the dominant routine application, followed by agricultural and related environmental issues and resource
exploration. Thus, none of the sensors employed was designed or intended for
volcanological applications. Despite this, data from remote sensing satellites have
00846597/00/05150081$14.00
81
82
FRANCIS n ROTHERY
83
even when their dimensions are less than the spatial resolution of the imaging
system, and they frequently saturate sensors originally intended for much lower
temperature applications. Volcanoes also emit detectable thermal radiation in the
short wavelength infrared or even the visible parts of the spectrum. Lava extrusions at temperatures of ;10008Cwhich may be visibly radiant even in daylightare the most obvious volcanic thermal phenomena, but there are many
others. These include fumaroles emitting steam and other gases that may heat the
ground to several hundred degrees; pyroclastic flow deposits consisting of hot,
solid, fragmented material; columns of ash and gas blasted into the atmosphere
by explosive eruptions, or rising convectively; and lakes of hot water ponded in
the craters of active volcanoes.
Volcanic thermal anomalies were noted on images shortly after routine weather
satellite observations began. Early examples were the 1978 and 1981 eruptions
of Etna detected by the NOAA 5 satellite (Archambault et al 1980, Bonneville et
al 1985). On 19 November 1979, Landsat 3 obtained the first high-resolution,
short-wave infrared images of an active lava flow on Sierra Negra volcano in the
Galapagos (Rothery et al 1988, McClelland et al 1989).
These early observations were necessarily qualitative. They benefited from the
hitherto unprecedented capability of satellites to provide synoptic coverage of the
Earth, and thus to obtain images of the most remote volcanoes as easily as the
most accessible. Extraction of meaningful quantitative thermal data is complicated by problems of scale. Large lava flows may have widths of many tens of
meters and extend many kilometers in length, but their thermal structure is invariably complex, typically consisting of small, incandescent areas of newly exposed
high temperature material, surrounded by much larger areas of cooler, chilled
material. This means that the hot components of volcanic surfaces imaged by
satellite sensors are almost invariably smaller than the pixel (picture element) size
of the sensor. This problem is obvious with meteorological satellites, which typically have pixel sizes of 14 km, but it is also serious with the ;30 m pixels of
Landsat and similar satellites used for detailed investigations. Because the radiant
thermal energy detected by a satellite sensor is integrated over the entire area of
each pixel, estimation of the temperatures of hot volcanic surfaces requires simultaneous determination of the relative surface areas covered by components at
different temperatures.
For remotely sensed data at a single wavelength, no solution is possible. For
multispectral data (i.e. data recorded at more than one wavelength), however,
Dozier (1981) and Matson & Dozier (1981) proposed an elegant solution based
on the Planck distribution law, which relates the spectral radiance (L) at wavelength k to its thermodynamic temperature T:
L(k,T ) 4 ekc1k15/[exp(c2/kT) 1 1]
(1)
where ek is the spectral emissivity and c1 and c2 are constants. Working with forest
fires, and assuming there are two components at different temperaturesone
hot and the other backgroundMatson and Dozier recognized that the radi-
FRANCIS n ROTHERY
84
ance measured by the satellite sensor at two different wavelengths could be used
to estimate both the temperature of each component and the fractional area of
each pixel covered by each component. For any one wavelength k (or channel of
a sensor), the total radiance detected by the sensor is the weighted average of the
radiance from each component:
Rk(tot) 4 s[ f (R hot) ` 1 1 f )R back]
(2)
where R hot is the radiance of the hot component at temperature T (K) derived
from Plancks equation; Rback is the radiance of the cooler background component,
also specified by Plancks equation; f is the fractional area of the surface covered
by the hot component; and s is the atmospheric transmittance function for the
wavelength concerned. Where radiance data at two infrared spectral channels with
central wavelengths kx and ky are available, Equation 2 can be written:
RX 4 sX[ f L (kX, Thot) ` (1 1 f )L(kX, Tback)]
(3)
(4)
and
Because Rx and Ry are recorded by the satellite, two out of the three remaining
parameters f, Thot and Tback, can be determined if the third can be estimated or
assumed, and if sx and sy can be constrained from standard atmospheric radiative
transfer models.
Equations 3 and 4 form the basis of the dual band method, which was first
applied to volcanoes by Rothery et al (1988), after they noted a persistent hot
spot in a summit crater of Lascar volcano (Figure 1, color insert), north Chile
(Francis & Rothery 1987). Subsequently, the same technique has been refined
and applied to volcanoes around the world by many authors, in a continuing string
of papers usefully reviewed by Flynn et al (2000). Rothery et al (1988) assumed
a value for the low-temperature background and used the technique to estimate
the temperature of the high-temperature material; they thus demonstrated that rock
at near magmatic temperatures had been present at the surface 21 months before
Lascar erupted in September 1986, following a prolonged period of repose. In a
significant refinement of the method, Oppenheimer (1991) pointed out that
because the core temperatures of typical basaltic lava flows vary little (they are
typically around 11008C), improved estimates of crust temperatures and their
fractional areas could be obtained by assuming values for core rather than background temperatures. Oppenheimer (1993) also showed that improved estimates
of radiant emittance could be obtained by using data in three rather than two
simultaneously acquired bandsyielding estimates for the core, crust, and cool
background to a lava flow. This has obvious implications for providing input to
models of lava flow behavior, which can be applied either retrospectively (to seek
to understand how a flow has behaved) or used in a real-time monitoring situation
to predict the motion of a flow. Other refinements include the work of Harris et
al (1998), who used hyperspectral data collected in the field to show that a high-
85
temperature lava surface is itself composed of at least two thermal sources, and
incorporated this into their remote sensing models.
We note that developments in the acquisition and interpretation of remotely
sensed thermal data for volcanoes on Earth have been paralleled by similar
advances in the study of the active volcanoes on Jupiters satellite Io (e.g. Davies
1996, McEwen et al 1998). Although the styles of volcanism on Io are very
different from those on Earth, the data sets used are also obtained by remote
sensing (from the Hubble Space telescope and the Galileo orbiter), and identical
problems of pixel size arise.
Shortly after discovering that volcanic thermal anomalies could be detected
on Landsat images, Francis et al (1989) argued that changes in radiated emitted
power with time might provide a valuable means of monitoring the level of activity of volcanoes. This is desirable, because active vents are generally not accessible to conventional methods, even when the volcanoes themselves are
accessible. In a first attempt to explore this approach, Glaze et al (1989) used
Landsat Thematic Mapper data to estimate the radiant thermal energy flux of
three active volcanoes: ErtaAle (Ethiopia), Erebus (Antarctica), and Lascar
(North Chile). Their approach yielded underestimates, since they did not take into
account factors such as the power radiated by the warm background. Radiation
is of course only one of several mechanisms by which volcanoes lose heat. Other
important factors are the convective and conductive fluxes; rainfall (which can
be modeled); and cooling by hydrothermal circulation, which is much more difficult to constrain. Subsequently, successively more sophisticated attempts to use
remote sensing to measure volcanic thermal budgets have been made, leading to
some richer insights into volcanic processes.
86
FRANCIS n ROTHERY
sistent with the rate at which lava has been erupted at the surface since 1924:
;104 kg s11. During the nineteenth century, when the lava lake existed, lava was
erupted as flows outside the confined lake at an average rate of only ;3 2
102 kg s11. Francis et al (1993) thus concluded that Kilauea and similar persistently active volcanoes must grow endogenously, either by intrusion of dykes or
by cumulate formation deep within the volcano. They also noted that efficient
convective circulation can be driven by density differences between gas-rich
magma and magma that has degassed via a persistently active vent. Thus, the
resulting density increase could inhibit eruption of lavas and promote cumulate
formation at depth.
Petrologic and tectonic studies of Kilauea have demonstrated that it is indeed
underlain by large volumes of cumulate olivine, which may be contributing to
the lateral spreading of the shield volcano by plastic deformation (Clague et al
1994). Volcanologists need to take into account the rapid rates of endogenous
growth implied by thermal budgets when considering the long-term evolution of
volcanic edifices. These may be spreading sideways or sagging downward by
lithospheric loading, as well as growing more conventionally upward by the accumulation of surface lava flows.
Several papers subsequently pursued this approach, demonstrating how basic
measurements done by remote sensing can stimulate new ways of looking at
volcanoes, or lend support to developing models. Harris & Stevenson (1997a,b)
quantified the thermal budgets of Vulcano and Stromboli, and Stevenson & Blake
(1998) built on this, using theoretical and laboratory modeling to investigate the
role of convective overturn of degassed magma in sustaining persistent volcanism.
Harris et al (1999) used three bands of Landsat TM data to constrain the radiated
heat flux and from this derived estimates of the convective and conductive losses
in order to find the total heat budget. They concluded that the persistence of lava
lakes over decades or even centuries may not necessarily imply endogenous or
cryptic growth, but that lakes could be sustained by magma recycling from a
sufficiently large, deep reservoirif it were hotter than the temperature at which
crystallization would commence in the magma chamber. Prolonged activity (with
or without crystallization) must be accounted for by circulation of magma
between the surface conduit and a deep reservoir, driven either by thermal convection or by degassing. In either case, the rate at which magma is required to
solidify permanently is much less than the rate at which magma is circulating
within the conduit. Along similar lines, Delmelle et al (in press) cited the remotely
sensed heat flux of the Masaya volcano (Nicaragua)which has an open, persistently degassing ventas evidence that the volumetric rate of magma degassing is several orders of magnitude greater than the rate at which magma is
solidifying in the conduit system. They concluded that gas-rich magma continually arrives at a shallow, vesiculated layer where it degasses and then sinks (with
minimal cooling), to be replaced by the next influx of gas-rich magma arriving
from depth.
87
88
FRANCIS n ROTHERY
89
local topographic data, assessments of areas likely to be affected. To date, however, only retrospective analyses of existing data have been carried out. To validate the effective operational use of the technique, a demonstration experiment
is required, with real-time data being used to predict the future evolution of a
flow field.
Effusion rate monitoring by thermal remote sensing is not limited solely to
lava flowsit can also be valuable for measuring effusion rates of lava domes.
This is important because pyroclastic flows triggered by dome collapse are much
more deadly than lava flows, and such flows are likely to be especially frequent
when a dome is expanding as a result of injection of lava. In a study of the growth
and collapse of the dacite lava dome of Mt. Unzen (Japan), Wooster & Kaneko
(1998) showed that the 1.6 lm spectral radiance as determined from nighttime
ATSR data correlates with the rate of degassing, which in turn is closely related
to the lava effusion rate into the dome.
90
FRANCIS n ROTHERY
surements and cloud cover assessment updated every 1530 minutes. By selecting
a threshold value for the difference in temperature between different channels,
alarm algorithms have been developed to detect new hotspots automatically and
send e-mail alert messages to volcanologists.
On 15 September 1998, an eruption of the volcano Cerro Azul on Isla Isabela
in the Galapagos islands was successfully detected in this way, four hours before
it was seen by ground-based observers. This was the first demonstration of automated eruption detection, and it highlights the future potential of automated systems (Harris et al, in press). Use of the Internet also means that for the first time,
volcano observatories around the world will have access to free, real-time data
for volcanoes that may be inaccessible for physical or financial reasons.
A similar approach has been adopted by the Alaska Volcano Observatory,
which is responsible for monitoring one of the worlds largest volcanic provinces.
Extending over 2500 km, this chain includes about a hundred potentially active
volcanoes. More than 350 eruptions have been recorded from 45 different volcanoes over the last 200 years (McClelland et al 1989). Although the area is thinly
populated, traffic in several major air corridors is at risk from eruptions. The
region is too far north to be well seen from a geostationary orbit; however, the
advantage is that the ground tracks of polar orbiting satellites are so closely spaced
that there are about twenty AVHRR imaging opportunities each day, which provides ample coverage. AVO has developed automated monitoring algorithms
(Dean et al 1998) that successfully detected an eruption of Okmok volcano, well
before local observers confirmed the effusion of new lava.
Ideally, remote monitoring of volcanoes requires long time series of data, so
that patterns of activity can be discerned and used for forecasting. Lascar volcano
in north Chile was the first thermally radiant volcano to be identified (Francis &
Rothery 1987) so it has the longest data coverage. Wooster & Rothery (1997)
used the Along Track Scanning Radiometer (ATSR) aboard the ERS-1 spacecraft
to monitor the volcano for the period 199295, extending the coverage and
improving on the temporal resolution of the previous TM time series (Oppenheimer et al 1993). Designed for routine mapping of sea surface temperatures,
ATSR is broadly similar to the AVHRR in that it has roughly 1 km pixel size,
but it is better suited to monitoring radiance from high-temperature volcanic
sources.
Fortuitously, Wooster & Rotherys time series of data spanned the largest eruption of Lascar for ;9,000 years. They showed that variations in radiance from
Lascar were consistent with an independent model (Matthews et al 1997) for the
growth and subsidence of a lava dome within the summit crater. They suggested
that large explosive eruptions are more likely when the short wavelength infrared
signal has fallen from a high to a low level, because this signifies that the lava
has collapsed back into the conduit, hindering passive degassing. In practice, this
conclusion may be difficult to use in forecasting, because decreases in thermal
radiance also precede decreases in the rate of supply of gas or magma from depth.
91
Radar Experiments
To investigate the potential of SAR; techniques for volcano studies, the aircraftborne Topographic Synthetic Aperture Radar (TOPSAR) system was deployed
to Hekla (Iceland), Vesuvius (Italy) and the Galapagos islands (Evans et al 1992,
FRANCIS n ROTHERY
92
Mouginis-Mark & Garbeil 1993, Rowland 1996). The aircraft carried a pair of
C-band (5.6 cm) antennae, which provided height measurements to a vertical
accuracy of 15 m and a spatial resolution of about 10 m. In July 1991, TOPSAR
data collected over Hekla were used to create a Digital Elevation Model (DEM),
which demonstrated that small-scale topography and local slopes could be derived
to high precision. Pressure ridges less than 2 m high were detectable. Such
detailed topography of lava flows is difficult to obtain by conventional means,
even for easily accessible volcanoes. A similar TOPSAR study of Fernandina
volcano, Galapagos islands, provided a topographic database for the entire island
volcano, which far excels the best previously available photogrammetric
topography.
Many experiments were also carried out with the Space Shuttle, notably the
Shuttle Imaging Radar (SIR-A, -B and -C) and the 1999 Shuttle Radar Topography Missions (SRTM). Gaddis (1992) investigated the effects of surface texture
and roughness on the intensity and polarization of backscattered radar signals,
and showed that differences could be used to discriminate features such as aa
and pahoehoe lava flows: aa and blocky lava flows are radar bright, whereas
smoother pahoehoe flows are radar dark. During the 1994 SIR-C missions, experiments were conducted to investigate the optimal combinations of radar wavelength, pulse bandwidth, polarization, look direction, and incidence angle for a
variety of scientific applications. Studies of the Virunga volcanic chain (Central
Africa) by MacKay & Mouginis-Mark (1997) provided valuable insights into the
potential applications of radar remote sensing to volcanoes in general. Data at 3.1
cm (X band), 5.7 cm (C band), and 24 cm (L band) were collected during the
experiment. As an example of the complexity of radar data sets, Table 1 summarizes the broad range of data options available during the Virunga experiment.
Although complex, multi-parameter radar experiments such as SIR-C and
SRTM are instructive, they are not suitable for routine observations: Space Shuttle
missions are costly, generally last only about ten days, and do not provide global
TABLE 1 Radar parameters used during the 1994 SIR-C experiment over the Virunga
volcanoes. Data from MacKay & Mouginis-Mark (1997). SRL-1: Shuttle Radar Lab-1, April
1994; SRL-2: Shuttle Radar Lab-2, October 1994.
Data
take
Mission
58.60
58.61
154.90
170.90
170.10
58.61
SRL-1
SRL-2
SRL-2
SRL-2
SRL-1
X-SAR
Polarization
Bandwidth
(MHz)
Incidence
angle
Look
direction
Swath
width
(km)
Range
resolution
(m)
HH,HV
HH,HV
VV
VV
HH
X-VV
20
10
40
40
20
10
36.7 42.78
37.7 48.18
20.030.88
20.030.88
34.5 46.48
41.6 47.18
Left
Left
Right
Right
Right
Left
36.5
71.4
46.4
46.4
71.6
60.2
12.511.1
24.520.2
11.07.3
11.07.3
13.210.4
23.821.6
93
94
FRANCIS n ROTHERY
95
correspond to the magma chamber. Now that the means to correct for tropospheric
effects are in place, radar interferometry on a volcano-wide scale is set to become
a routine operational technique for volcano monitoring (Figure 3, see color insert).
Several studies have already capitalized on this remarkable new approach. For
example, Sigmundsson et al (1997) observed subsidence of ;24 mm per year
over the shallow magma chamber at Krafla, Iceland, superimposed on a ;7 mm
per year subsidence along the axis of the spreading segment. They attributed this
to cooling contraction at a depth of about 3 km, coupled with ductile flow of
material away from the spreading axis. Lu et al (1998) documented 140 cm of
subsidence of the caldera floor of Okmok volcano (Alaska) as a result of its
FebruaryApril 1997 eruption.
Although satellite DINSAR holds the promise of exceptionally detailed topographic and deformation studies of volcanoes, it is not as widely applicable as
first hoped. Observations of the ongoing eruption of the Soufriere Hills volcano,
Montserrat, show that it can be difficult to obtain suitable coherent radar images
for many reasons. Rapid variations in the vegetation canopy, differences in surface
dielectric constant due to variations in moisture, steep, irregular, and rapidly
changing topography on the growing and collapsing andesitic lava dome, and
variations in atmospheric transmission caused by variations in water vapour all
conspire to hinder the retrieval of topography. At Montserrat, useful topographic
data were obtained from fans of relatively flat, smooth, unvegetated pyroclastic
flow deposits surrounding the dome, but little coherent data could be derived
from the active dome itself (Wadge et al 2000).
ERUPTION CLOUDS
Explosive eruptions inject clouds containing huge quantities of ash, gas and water
vapour into the troposphere and stratosphere (Figure 4, see color insert). Sparks
et al (1997) provides a useful general volcanological overview. Scientific issues
apart, a crucial practical problem for the remote sensing community is the discrimination of ash clouds from ordinary clouds, because of the threat they pose
to aviation. Disasters were narrowly averted when aircraft unwittingly encountered ash clouds from eruptions of Galunggung in Indonesia, 1982, and of
Redoubt in Alaska, 1989 (Casadevall 1994a), which caused the aircraft engines
to stall. Even when, as in these instances, the engines are eventually restarted and
the aircraft lands safely, the cost of repairing the damage caused to the engines,
avionics, and fabric of an aircraft can reach several tens of millions of dollars in
the case of a Boeing 747 (Casadevall 1994b). An additional challenge is that
discrimination of clouds must be done in near real time, because injection of
eruption clouds to heights of many kilometers may take only minutes, while
aircraft may be approaching at speeds up to 1000 km hr11. Furthermore, aircraft
on long-haul polar or trans-Pacific routes have limited options for diversion to
alternative airports. Timely warnings of eruption clouds can help avoid the sig-
96
FRANCIS n ROTHERY
nificant economic losses associated with landing planeloads of passengers at unintended destinations, and leaving aircraft out of position. Satellite monitoring is
the only way that synoptic eruption cloud data can be gathered rapidly enough
to preempt such problems. To this end, the International Civil Aviation Organization sponsored in 1997 the establishment of nine Volcanic Ash Advisory Centers responsible for issuing Volcanic Ash Advisory Statements to aviators. These
rely heavily on images from meteorological satellites and are available on the
World Wide Web (Figure 5, see color insert). Solid ash particles from an eruption
cloud usually fall to the ground within a few tens of kilometers of the volcano,
except in the very largest cases where fine ash may fall hundreds of kilometers
away. Gases behave differently. If the eruption cloud is confined to the troposphere, it is relatively rapidly dispersed; if injected into the stratosphere ($18 km
at low latitudes; $10 km at high latitudes), however, the SO2 is oxidized to
sulphuric acid, which subsequently forms aerosols through heteromolecular
nucleation of H2SO4 and H2O vapors, at a rate depending on the amount of
available water vapor. The residence time of acid aerosols in the stratosphere is
long (many months), and the scattering to space and absorption of incident solar
radiation by aerosol particles can cause significant surface cooling for periods of
12 years. Aerosols are also important because they provide surfaces upon which
chemical reactions can take place, and may therefore play a significant role in
stratospheric ozone depletion.
Satellite monitoring is thus called on to track separately the silicate, gas, and
aerosol components of eruption plumes.
97
for mitigating aircraft hazards, but AVHRR images by themselves are not always
adequate.
Plume Heights
An important challenge in tracking eruption clouds is determining the height they
reach, because this information can be used to infer the underlying mass eruption
rate. Several approaches have been used, and are usefully summarized in Sparks
et al (1997). Stereo satellite images represent a straightforward method, but currently only the ATSR on the ERS satellites is capable of acquiring stereoscopic
views of transient features. This provides forward viewing and nadir viewing
scans, recorded about two minutes apart. Prata & Turner (1997) showed that these
data could determine cloud-top height to an accuracy of 51 km. We know of
only one study that has yet used this capability to investigate eruption plume
heights (Shocker et al 2000), but similar stereoscopic work is likely to form an
important part of future monitoring.
A crude measure of plume top height can be obtained from radiometric temperature measurements (e.g. Glaze et al 1989). These assume that the emissivity
of the plume is close to unity, and that the plume top temperature is in equilibrium
with the ambient air. Uncertainties arise because emissivities of ash clouds are
poorly known, and because the buoyant momentum of rapidly ascending convective clouds may carry them well above the level of thermal equilibrium
(Woods & Self 1992). For accurate work, it is also necessary to get an atmospheric
temperature profile for the region from radiosonde balloons. These are difficult
to obtain for many parts of the world. Measurements of plume shadow lengths
provide a trigonometric means of estimating height (e.g. Glaze et al 1989), but
only on daytime images. For large umbrella clouds, the value of these measurements is limited because they give heights of the margins of clouds, rather than
their highest central points.
Glaze et al (1999) combined shadow length and temperature determinations
with a photoclinometric study of the December 1989 eruption of Redoubt, Alaska.
They were able to show that the plume rose to an altitude of 3 km while drifting
150 km downwind from the vent. Using the shadow technique alone, they would
have systematically underestimated the plume height by up to 30%. This is a
significant finding, because an underestimate of the plume height by 30% could
lead to an underestimate of the eruption rate by a factor of three. In turn, this
could lead to major errors in predictions made from satellite data of the accumulation rate on the ground of ash falling from the plume, for example.
FRANCIS
ROTHERY C-1
Figure 1 Detail from a pair of Landsat TM images of the summit crater area of Lascar
volcano, Chile, recorded by day on 27 October (left) and by night on 28 November
(right) 1989. Three bands (or channels) of TM data have been used to compile these
views: band 7, 2.082.35 m, is in red; band 5, 1.551.75 m, is in green; and band 4,
0.760.90 m, is in blue. The area shown is about 2.5 km across in each case. The radiant anomaly in bands 5 and 7 is readily apparent by night, and represents a lava dome
surrounded by hot fumaroles. The anomaly is less apparent by day because the signal
is diluted by reflected sunlight.
C-2 FRANCIS
ROTHERY
Figure 2 Daytime Landsat TM images using the same band combination as in Figure
1, draped over a digital elevation model (viewed looking toward the west) and showing the development of an active lava flow field on Mt. Etna. Left: 2 January 1992;
right: 9 May 1992. The area shown is about 4.3 km across and 9.0 km from front to
back, and covers a scarp-bounded valley on the volcanos eastern flank (the summit is
out of view, beyond the back right corner of each image). The source of the lava flow
is where a fissure intersects the back wall of the valley, and the development of the
flow field was largely controlled by the form of this valley. Only the bright red-orange
areas are thermally radiant; all other areas show up because of reflected sunlight.
Snow appears dark blue, and receded to higher altitudes between the two dates. On 2
January, the lava flow field was at its most vigorous, with large expanses dominated
by radiance from incandescent surfaces. On the 9 May image, radiance is apparent
from open channels and skylights over lava tubes near the source area, and also where
the front of the flow field was still advancing toward the front left. Between these two
areas, the surface slope is much less and the flow field is covered by a relatively
unbroken crust. (Images courtesy of Rob Wright, Open University.)
FRANCIS
ROTHERY C-3
C-4 FRANCIS
ROTHERY
Figure 4 The early stages of the 19 September 1994 eruption of Rabual (Papua, New
Guinea). Top: Photographed looking obliquely westward from the Space Shuttle,
showing the eruption cloud fanning out downwind. Inset: Image from a geostationary
weather satellite (GMS) at 07.32 GMT. Land appears darker than the sea, enabling the
northern coastline of Australia to be made out.
FRANCIS
ROTHERY C-5
Figure 5 A Volcanic Ash Advisory Statement, in this example relying on geostationary weather satellite images to identify and track the plume.
C-6 FRANCIS
ROTHERY
Figure 6 A compilation of TOMS data sets for several successive days revealing
the circumpolar transport of the sulfur dioxide plume from the 1991 eruption of
Mt. Hudson, Chile (courtesy of WI Rose). The values on the color-coded scale are
Dobson units (the standard measure of trace gas column density in the atmosphere).
98
FRANCIS n ROTHERY
sing, because changes in gas composition may presage eruptions. Current satellite
technology, however, is restricted to detection and mapping of sulfur dioxide
when injected in large quantities into the atmosphereusually, but not always,
by large explosive eruptions. Water vapor plays a crucial role in the dynamics of
eruption columns, a role that is complicated by changes in phase that provide
much of the energy driving convection. Unfortunately, estimates of volcanic water
vapor are exceedingly difficult to obtain by remote sensing, because of water
vapors high ambient concentrations. Similarly, volcanic contributions of carbon
dioxide are interesting because of their involvement in the global carbon cycle,
but again they are very difficult to distinguish against normal background levels.
Considerable success has been achieved in using the Total Ozone Mapping
Spectrometer (TOMS) on the Nimbus-7 and Meteor 3 satellites to map the atmospheric distribution of volcanic SO2 (Figure 6, see color insert). Although it was
primarily designed to determine ozone concentrations by measuring backscattered solar ultraviolet irradiance in six narrow wavelength bands in the region
0.306 to 0.360 lm, TOMS is able to take advantage of strong absorptions by SO2
in the same part of the spectrum to measure SO2 concentrations. TOMS provides
a spatial resolution of 47 km 2 47 km at nadir, but averages 62 km. Thus, while
it is well adapted to obtaining estimates of bulk SO2 from large dispersed plumes,
its footprint is too large to permit observations of small-scale volcanic features
for example, the steady-state degassing of SO2 by volcanoes such as Etna or
Stromboli, whose average fluxes are only 4 and 0.8 ktonne/day, respectively. In
a review, Bluth et al (1993) estimated that TOMS successfully detected SO2 from
55 out of 350 known eruptions between 1979 and 1992, and also from several
eruptions not known from ground studies. They also suggested that explosive
volcanism contributes about 4 Mtonne of SO2 annually to the atmosphere, mostly
from a few large eruptions, out of a total of approximately 13 Mtonne (from
explosive and non-explosive eruptions). (This is about 510% of the current
anthropogenic flux.) We list some significant volcanic sulfur dioxide plumes in
Table 2.
Limb-sounding instruments (which look obliquely through the atmosphere,
rather than straight down at the ground) are capable of detecting trace gases at
much lower concentrations than the nadir-pointing TOMS, but the footprints of
these sensors are inevitably even larger. An unanticipated capability of the Upper
Atmosphere Research Satellite (UARS) was to obtain vertically resolved (2.5
5.0 km) daily global maps of volcanic SO2 and aerosols. The Microwave Limb
Sounder (MLS) aboard UARS began operating a hundred days after the June
1991 eruption of Mt. Pinatubo, and immediately detected SO2 from the eruption;
this showed that the gas was distributed in a band around the tropics, with peak
mixing ratios of 15 ppbv (parts per billion by volume) near 26 km altitude. An
initial input of 17 Mtonne was estimated, compared with the 20 Mtonne estimated
from TOMS (Read et al 1993). Although superior to TOMS in resolution, the
MLS, which can detect SO2 down to mixing ratios of 15 ppbv, is also sensitive
only to gas at altitudes greater than 15 km, so tropospheric gas is inaccessible.
99
TABLE 2. Sulfur dioxide plumes for the period 19791985 detected by the TOMS
instrument on Nimbus 7 (from Mouginis-Mark et al 1991). The mystery eruption is now
thought to be that of Nyamuragira in Zaire. Note that most eruptions were explosive, but those
of Mauna Loa and Nyamuragira were large, effusive events.
Volcano
Sierra Negra
Mt. St. Helens
Alaid
Ambrym
Pagan
Mystery
El Chichon
Galunggung
Soputan
Colo, Una Una
Mauna Loa
Fernandina
Soputan
Krafla
Ruiz
Ruiz
Eruption date
Number of days
plume tracked
Nov 13 1979
May 19 1980
April 27 1981
May 8 1981
May 15 1981
Dec 26 1981
March 28 1981
April 5, June 25, July 14 1982
Aug 27 1982
July 24 1983
March 25 1984
March 31 1984
May 25 1984
Sept 5,10, 18 1984
Sept 12 1984
Nov 13 1985
3
4
22
2
2
`12
`64
1, 1, 2
1
1
1
11
2
3,1,1
7
7
Total SO2/kton
340
1,000
5,000
275
190
35
660
100
FRANCIS n ROTHERY
Aerosol Measurements
Wen and Rose (1994) refined the data extraction by showing that AVHRR bands
4 and 5 could also be used to estimate particle sizes and masses in eruption clouds,
and Krotkov et al (1999) showed the power of using AVHRR and TOMS (or
equivalent) data in parallel to track the history of total particle mass and size
distribution during the evolution of an eruption cloud. In some eruptions, such as
that of Rabaul 1994 (see Figure 4), large quantities of ice crystals are thought to
form in rising eruption clouds. Further work is required to ascertain the effects
of ice on the retrieved quantities of ash in such clouds.
Measurements of stratospheric aerosols from volcanic eruptions are perhaps
of more relevance to atmospheric and climate scientists than to volcanologists.
AVHRR satellite data were used to map aerosols from the El Chichon and Pinatubo eruptions, by determining the errors introduced into apparent changes in
sea-surface temperature caused by the aerosols (Strong & Stowe 1993). The technique requires that the plume be transparent enough so that the sea surface is
visible. Thus, very dense, ash-laden plumes cannot be studied. A similar technique
may also be valuable in monitoring the steady state emissions from volcanoes
such as Etna, which rarely emits dense ash clouds.
The Stratospheric Aerosol and Gas Experiment (SAGE II) carried by UARS
scans the Earths sunrise and sunset limbs, with a vertical resolution of 1 km in
order to provide information on the distribution of aerosols. It tracked the dispersal of 1991 Pinatubo aerosol cloud into both northern and southern hemispheres, but was unable to discriminate between aerosol types (McCormick &
Viega 1992).
FUTURE PROSPECTS
There are remarkably bright prospects for volcano remote sensing. Future years
will bring the launch of ever-increasing numbers of remote sensing satellites from
several different countries. In the short term, Landsat 7, launched in 1999, will
extend the capability of this familiar system for many years. NASAs Earth
Observing System (EOS) will greatly enhance monitoring capabilities with both
high- and low-resolution multispectral sensors. Launch of the first platform in the
series (previously known as AM-1, but renamed Terra) is scheduled for December
1999 or later; launch of others in the series is due to continue in 2000 (EOS PM1) and 2002 (EOS CHEM-1). Volcanological investigations to be carried out
under the auspices of EOS will make use of a bewildering array of sensors:
ASTER, CERES, EOSP, ETM`, GLAS, HIRDLS, MISR, MODIS, OMI, SAGE
III, TES, TOMS and VCL. Many other satellites are in an advanced design stage,
such as ENVISAT and RADARSAT-2. Although the string of acronyms gives a
flavor of the range of investigations to be undertaken, it is impractical to review
them all. It is evident, however, that unprecedented amounts of satellite data will
101
102
FRANCIS n ROTHERY
surements of SO2 gas emissions can be combined with data from sensors designed
to measure aerosols (such as SAGE III), it may become possible to investigate
the rates at which sulfate aerosols are formed from volcanic gases in different
contexts, with important implications for understanding the radiative effects of
volcanic aerosols on climate.
Although volcanology has been included in the objectives for EOS since its
inception, volcanological science has necessarily formed a relatively minor component of the overall mission strategy. The competing requirements of different
scientific disciplines dictate that the sensors used for volcanology represent a
compromise with a range of other applications. Eventually, a satellite dedicated
to volcanology will be launched. Such a platform (VOLCAM) was conceived in
1998, and although not ultimately approved for funding, its successors are likely
to embody similar concepts.
VOLCAM was designed primarily with aviation hazards in mind; therefore,
its UV and IR sensors were intended to optimize the ash and sulfur dioxide
capabilities of existing AVHRR and TOMS for tracking volcanic clouds and for
discriminating between volcanic and meteorological clouds. Intended as a geostationary satellite located over the Americas, it would have had the advantage
of 15-minute time resolution. Through optimization of the wavelength bands in
the UV camera, the sensitivity to sulfur dioxide was improved by a factor of ten
over TOMS. The VOLCAM footprint is also smaller by a factor of 16, so much
smaller clouds could be detected. Thus, in addition to determining whether clouds
were safe to fly through, VOLCAM also offered the potential to monitor preeruptive outgassing from all the volcanoes viewed from the satellite.
A satellite such as VOLCAM would permit investigations of some key science
issues, such as the fate of SO2 in eruption clouds. In fresh clouds, SO2 may be
lost by absorption onto ash or ice particles (Rose et al 1995), but few direct data
bear on this theory. In older clouds (more than a day), sulfur dioxide is probably
converted to sulfate in gas phase reactions. Measurements of the lifetime of SO2
in clouds of different ages are required to resolve the relative rates of these
processes.
Also required is a high-resolution thermal system in polar orbit, with spectral
bands and sensor gains tailored for volcanologic applications. Such a satellite
could be rapidly targeted on volcanic hotspots to investigate hot spot alarms
triggered by the low-resolution thermal alert system that is already coming into
place.
POSTSCRIPT
Peter Francis died suddenly on 30 October 1999 at the age of 54, just three weeks
after acting as convenor for a highly successful international meeting on eruptions
of andesite volcanoes. He was a pioneer in satellite remote sensing of volcanoes
and in the use of Fourier Transform Infrared spectroscopy for volcanic gas analysis. He will be greatly missed.
Visit the Annual Reviews home page at www.AnnualReviews.org.
103
LITERATURE CITED
Archambault C, Stoschek J, Tanguy JC. 1980.
Analyse et surveillance thermique de volcans actifs (197879) Bull. PERPSEV No.
11, Paris
Bluth GJS, Schnetzler CC, Kreuger AJ, Walter
LS, et al. 1993. The contribution of explosive volcanism to global atmospheric sulphur dioxide concentrations. Nature
366:32729
Bonneville A, Vasseur G, Kerr Y. 1985. Satellite thermal infrared observations of
Mt. Etna after the 17th March 1981 eruption. J. Volcanol. Geotherm. Res. 24:293
313
Casadevall TJ. 1994b. The 19891990 eruption of Redoubt Volcano, Alaska: implications on aircraft operations. J. Volcanol.
Geotherm. Res. 62:30116
Casadevall TJ, ed. 1994a. Volcanic hazards and
aviation safety. In Proc. First Internat.
Symp. on Volcanic Ash and Aviation Safety.
US Geol. Surv. Bull. No. 2047. 450 pp.
Washington, DC: US Govt. Print. Off.
Clague DA, Denlinger RP, McGee KA. 1994.
Role of olivine cumulates in destabilizing
the flanks of Hawaiian volcanoes. Bull. Volcanol. 56:42534
Davies AG. 1996. Ios volcanism: thermophysical models of silicate lava compared
with observations of thermal emission. Icarus 124:4561
Dean K, Servilla M, Roach A, Foster B, Engle
K, et al. 1998. Satellite monitoring of
remote volcanoes improves study efforts in
Alaska. EOS Trans. Am. Geophys. Union.
79:41318
Delacourt C, Briole P, Achache J. 1998. Tropospheric corrections of SAR inteferograms with strong topography. Application
to Etna. Geophys. Res. Lett. 25:284952
Delmelle P, Baxter P, Beaulieu A, Boardman
L, Burton M, et al. 1999. Interdisciplinary
attention focuses on the restless Masaya
caldera, Nicaragua. EOS Trans. Am. Geophys. Union. In press
Dozier JA. 1981. A method for satellite identification of surface temperature fields of
sub-pixel resolution. Rem. Sens. Environ.
11:22129
Evans DL, Farr T, Zebker HA, van Zyl JJ,
Mouginis-Mark PJ, et al. 1992. Radar interferometry studies of the Earths topography.
EOS Trans. Am. Geophys. Union 73:553
58
Flynn LP, Harris AJL, Rothery DA, Oppenheimer C, et al. 2000. High spatial resolution remote sensing of active volcanic
features using Landsat and hyperspectral
data. In AGU monograph, Remote Sensing of Volcanoes, ed. P Mouginis-Mark, J
Fink, L Glaze. Washington, DC: AGU. In
press
Francis PW, Glaze LS, Rothery DA. 1989. Lascar volcano set to erupt. Nature 339: 434
Francis PW, Oppenheimer C, Stevenson D.
1993. Endogenous growth of persistently
active volcanoes. Nature 366:55457
Francis PW. 1989. Remote sensing of volcanoes. Adv. Space Res. 9:8992
Francis PW, Rothery DA. 1987. Using the
Landsat Thematic Mapper to detect and
monitor active volcanoes: an example from
the Lascar volcano, northern Chile. Geology 15:61417
Francis PW, Wadge G, Mouginis-Mark PJ.
1996. Satellite monitoring of volcanoes in
IAVCEI-UNESCO Manual of Volcano
Monitoring, ed. R Scarpa, R Tilling. pp.
257298. Berlin: Springer-Verlag
Friedman JD, Williams RS. 1970. Changing
pattern of thermal emission from Surtsey,
Iceland, between 1966 and 1969. USGS
Prof. Pap. 700-D:D11624
Gaddis LR. 1992. Lava flow characterization
at Pisgah volcanic field, California, with
multiparameter imagining radar. Geol. Soc.
Am. Bull. 104:695703
Garvin JB. 1996. Topographic characterization
and monitoring of volcanoes via airborne
laser altimetry. In Volcano Instability on the
104
FRANCIS n ROTHERY
105
106
FRANCIS n ROTHERY
CONTENTS
Palynology after Y2K--Understanding the Source Area of Pollen in
Sediments, M. B. Davis
Dinosaur Reproduction and Parenting, John R. Horner
Evolution and Structure of the Lachlan Fold Belt (Orogen) of Eastern
Australia, David A. Foster, David R. Gray
1
19
47
81
107
141
169
211
281
305
339
367
391
419
477
509
539
571
611
669