Trends in the aggregated rate of pre-eruptive volcano-tectonic seismicity at Kilauea volcano, Hawaii

Accelerating rates of volcano-tectonic (VT) earthquakes are commonly observed during volcanic unrest. Understanding the repeatability of their behaviour is essential to evaluating their potential to forecast eruptions. Quantitative eruption forecasts have focused on changes in precursors over intervals of weeks or less. Previous studies at basaltic volcanoes in frequent eruption, such as Kilauea in Hawaii and Piton de La Fournaise on Réunion, suggest that VT earthquake rates tend to follow a power-law acceleration with time about 2 weeks before eruption, but that this trend is often obscured by random fluctuations (or noise) in VT earthquake rate. These previous studies used a stacking procedure, in which precursory sequences for several eruptions are combined to enhance the signal from an underlying acceleration in VT earthquake rate. Such analyses assume a common precursory trend for all eruptions. This assumption is tested here for the 57 eruptions and intrusions recorded at Kilauea between 1959 and 1984. Applying rigorous criteria for selecting data (e.g. maximum depth; restricting magnitudes to be greater than the completeness magnitude, 2.1), we find a much less pronounced increase in the aggregate rate of earthquakes than previously reported. The stacked trend is also strongly controlled by the behaviour of one particular pre-eruptive sequence. In contrast, a robust signal emerges among stacked VT earthquake rates for a subset of the eruptions and intrusions. The results are consistent with two different precursory styles at Kilauea: (1) a small proportion of eruptions and intrusions that are preceded by accelerating rates of VT earthquakes over intervals of weeks to months and (2) a much larger number of eruptions that show no consistent increase until a few hours beforehand. The results also confirm the importance of testing precursory trends against data that have been filtered according to simple constraints on the spatial distribution and completeness magnitude of the VT earthquakes.     

Seismic pattern recognition techniques to predict large eruptions at the Popocatépetl,Mexico, volcano

Using pattern recognition techniques, we formulate a simple prediction rule for a retrospective prediction of the three last largest eruptions of the Popocatépetl, Mexico, volcano that occurred on 23 April–30 June 1997 (Eruption 1; VEI ~ 2–3); 11 December 2000–23 January 2001 (Eruption 2; VEI ~ 3–4) and 7 June–4 September 2002 (Eruption 3; explosive dome extrusion and destruction phase). Times of Increased Probability (TIP) were estimated from the seismicity recorded by the local seismic network from 1 January 1995 to 31 December 2005. A TIP is issued when a cluster of seismic events occurs under our algorithm considerations in a temporal window several days (or weeks) prior to large volcanic activity providing sufficient time to organize an effective alert strategy. The best predictions of the three analyzed eruptions were obtained when averaging seismicity rate over a 5-day window with a threshold value of 12 events and declaring an alarm for 45 days. A TIP was issued about six weeks before Eruption 1. TIPs were detected about one and four weeks before Eruptions 2 and 3, respectively. According to our objectives, in all cases, the observed TIPs would have allowed the development of an effective civil protection strategy. Although, under our model considerations the three eruptive events were successfully predicted, one false alarm was also issued by our algorithm. An analysis of the epicentral and depth distribution of the local seismicity used by our prediction rule reveals that successful TIPs were issued from microearthquakes that took place below and towards SE of the crater. On the contrary, the seismicity that issued the observed false alarm was concentrated below the summit of the volcano. We conclude that recording of precursory seismicity below and SE of the crater together with detection of TIPs as described here, could become an important tool to predict future large eruptions at Popocatépetl. Although our model worked well for events that occurred in the past, it is necessary to verify the real capability of the model for future eruptive events.     

Seismic activity related to the 1991 eruption of Colima Volcano,Mexico

Ten years after the last effusive eruption and at least 15 years of seismic quiescence, volcanic seismic activity started at Colima volcano on 14 February 1991, with a seismic crisis which reached counts of more than 100 per day and showed a diversity of earthquake types. Four other distinct seismic crises followed, before a mild effusive eruption in April 1991. The second crisis preceded the extrusion of an andesitic scoriaceous lava lobe, first reported on 1 March; during this crisis an interesting temporary concentration of seismic foci below the crater was observed shortly before the extrusion was detected. The third crisis was constituted by shallow seismicity, featuring possible mild degassing explosion-induced activity in the form of hiccups (episodes of simple wavelets that repeat with diminishing amplitude), and accompanied by increased fumarolic activity. The growth of the new lava dome was accompanied by changing seismicity. On 16 April during the fifth crisis which consisted of some relatively large, shallow, volcanic earthquakes and numerous avalanches of older dome material, part of the newly extruded dome, which had grown towards the edge of the old dome, collapsed, producing the largest avalanches and ash flows. Afterwards, block lava began to flow slowly along the SW flank of the volcano, generating frequent small incandescent avalanches. The seismicity associated with the stages of this eruptive activity shows some interesting features: most earthquake foci were located north of the summit, some of them relatively deep (7–11 km below the summit level), underneath the saddle between the Colima and the older Nevado volcanoes. An apparently seismic quiet region appears between 4 and 7 km below the summit level. In June, harmonic tremors were detected for the first time, but no changes in the eruptive activity could be correlated with them. After June, the seismicity decreasing trend was established, and the effusive activity stopped on September 1991.     

Overview of the 1990–1995 eruption at Unzen Volcano

Following 198 years of dormancy, a small phreatic eruption started at the summit of Unzen Volcano (Mt. Fugen) in November 1990. A swarm of volcano-tectonic (VT) earthquakes had begun below the western flank of the volcano a year before this eruption, and isolated tremor occurred below the summit shortly before it. The focus of VT events had migrated eastward to the summit and became shallower. Following a period of phreatic activity, phreatomagmatic eruptions began in February 1991, became larger with time, and developed into a dacite dome eruption in May 1991 that lasted approximately 4 years. The emergence of the dome followed inflation, demagnetization and a swarm of high-frequency (HF) earthquakes in the crater area. After the dome appeared, activity of the VT earthquakes and the summit HF events was replaced largely by low-frequency (LF) earthquakes. Magma was discharged nearly continuously through the period of dome growth, and the rate decreased roughly with time. The lava dome grew in an unstable form on the shoulder of Mt. Fugen, with repeating partial collapses. The growth was exogenous when the lava effusion rate was high, and endogenous when low. A total of 13 lobes grew as a result of exogenous growth. Vigorous swarms of LF earthquakes occurred just prior to each lobe extrusion. Endogenous growth was accompanied by strong deformation of the crater floor and HF and LF earthquakes. By repeated exogenous and endogenous growth, a large dome was formed over the crater. Pyroclastic flows frequently descended to the northeast, east, and southeast, and their deposits extensively covered the eastern slope and flank of Mt. Fugen. Major pyroclastic flows took place when the lava effusion rate was high. Small vulcanian explosions were limited in the initial stage of dome growth. One of them occurred following collapse of the dome. The total volume of magma erupted was 2.1×10⁸ m³ (dense-rock-equivalent); about a half of this volume remained as a lava dome at the summit (1.2 km long, 0.8 km wide and 230–540 m high). The eruption finished with extrusion of a spine at the endogenous dome top. Several monitoring results convinced us that the eruption had come to an end: the minimal levels of both seismicity and rockfalls, no discharge of magma, the minimal SO₂ flux, and cessation of subsidence of the western flank of the volcano. The dome started slow deformation and cooling after the halt of magma effusion in February 1995.     

Precursors to dyke-fed eruptions at basaltic volcanoes: insights from patterns of volcano-tectonic seismicity at Kilauea volcano,Hawaii

To investigate the physical controls on volcano-tectonic (VT) precursors to eruptions and intrusions at basaltic volcanoes, we have analyzed the spatial and temporal patterns of VT earthquakes associated with 34 eruptions and 23 dyke intrusions that occurred between 1960 and 1983 at Kilauea, in Hawaii. Eighteen of the 57 magmatic events were preceded by an acceleration of the mean rate of VT earthquakes located close to the main shallow magma reservoir. Using a maximum-likelihood technique and the Bayesian Information Criterion for model preference, we demonstrate that an exponential acceleration is preferred over a power-law acceleration for all sequences. These sequences evolve over time-scales of weeks to months and are consistent with theoretical models for the approach to volcanic eruptions based on the growth of a population of fractures in response to an excess magma pressure. Among the remaining 40 magmatic events, we found a significant correlation between swarms of VT earthquakes located in the mobile south-flank of Kilauea and eruptions and intrusions. The behaviour of these swarms suggests that at least some of the magmatic events are triggered by transient episodes of elevated rates of aseismic flank movement, which could explain why many eruptions and intrusions are not preceded by longer-term precursory signals. In none of the 57 cases could a precursory sequence be used to distinguish between the approach to an eruption or an intrusion, so that, even when a precursory sequence is recognized, there remains an empirical chance of about 40% (24 intrusions from 57 magmatic events) of issuing a false alarm for an imminent eruption.     

Precursory seismicity of the 1994 eruption of Popocatépetl Volcano,Central Mexico

Popocatépetl Volcano is located in the central Mexican Volcanic Belt, within a densely populated region inhabited by over 20 million people. The eruptive history of this volcano indicates that it is capable of producing a wide range of eruptions, including Plinian events. After nearly 70 years of quiescence, Popocatépetl reawakened in December 21, 1994. The eruptive activity has continued up until the date of this submission and has been characterized by a succession of lava dome growth-and-destruction episodes, similar to events that have apparently been typical for Popocatépetl since the fourteenth century. In this regime, the episodes of effusive and moderately explosive activity alternate with long periods of almost total quiescence. In this paper we analyze five years of volcano-tectonic seismicity preceding the initial eruption of the current episode. The evolution of the V-T seismicity shows four distinct stages, which we interpret in terms of the internal processes which precede an eruption after a long period of quiescence. The thermal effects of a magma intrusion at depth, the fracturing related to the slow development of magma-related fluid pathways, the concentration of stress causing a protracted acceleration of this process, and a final relaxation or redistribution of the stress shortly before the initial eruption are reflected in the rates of V-T seismic energy release. A hindsight analysis of this activity shows that the acceleration of the seismicity in the third stage asymptotically forecast the time of the eruption. The total seismic energy release needed to produce an eruption after a long period of quiescence is related to the volume of rock that must be fractured so imposing a characteristic threshold limit for polygenetic volcanoes, limit that was reached by Popocatépetl before the eruption.     

Earthquake activity related to the 1991 eruption of the Hekla volcano,Iceland

The 1991 eruption of the Hekla volcano started unexpectedly on 17 January. No long-term precursory seismicity was observed. The first related activity was a swarm of small earthquakes that began approximately half an hour before the eruption. Intensive seismicity, both earthquakes and volcanic tremor, accompanied the violent onset of the eruption. Almost 400 events up to M_L magnitude 2.5 were recorded during the first few hours. During the later phases of the eruption, the earthquake activity was modest and the main volcano-related seismic signal was the persistent volcanic tremor. The tremor died away, together with the eruption on 11 March, and Hekla was seismically quiet until the beginning of June 1991, when a sudden swarm of numerous small shallow earthquakes occurred. This activity is atypical for Hekla and is interpreted to be a failed attempt to resume the eruption.     

Seismic energy releases from volcanoes

Seismic energy release during the precursory, eruptive and declining stages of volcanic activities provides various information about the mechanisms of volcanic eruptions and the temporary developments of their activities. Hitherto the energy release patterns from precursory earthquake swarms were used to predict the eruption times, especially of andesitic or dacitic volcanoes. In this paper the discussion is expanded to quantify the total amount of seismic energy released at the threshold of volcanic eruptions, with reference to the results observed at several volcanoes. The results generally indicate that the cumulative seismic energy release from the precursory earthquake swarms exceed 10¹⁷¹⁸ergs before eruptions at any andesitic or dacitic volvanoes. This allows the seismic efficiency, or the ratio of energy radiated seismically, and the energy required for the volumetric expansion to be estimated by incorporating available deformation data with the seismic data. The dependency of seismic efficiency on the type of volcanic activity, i. e. non-explosive outbreaks, phreatic and magmatic eruptions, dome formation, etc., was evaluated from observations at a few volcanoes that provided a variety of examples.     

Subaqueous, basaltic lava dome and carapace breccia on King George Island, South Shetland Islands, Antarctica

 On King George Island during latest Oligocene/earliest Miocene time, submarine eruptions resulted in the emplacement of a small (ca. 500 m estimated original diameter) basalt lava dome at Low Head. The dome contains a central mass of columnar rock enveloped by fractured basalt and basalt breccia. The breccia is crystalline and is a joint-block deposit (lithic orthobreccia) interpreted as an unusually thick dome carapace breccia cogenetic with the columnar rock. It was formed in situ by a combination of intense dilation, fracturing and shattering caused by natural hydrofracturing during initial dome effusion and subsequent endogenous emplacement of further basalt melt, now preserved as the columnar rock. Muddy matrix with dispersed hyaloclastite and microfossils fills fractures and diffuse patches in part of the fractured basalt and breccia lithofacies. The sparse glass-rich clasts formed by cooling-contraction granulation during interaction between chilled basalt crust and surrounding water. Together with muddy sediment, they were injected into the dome by hydrofracturing, local steam fluidisation and likely explosive bulk interaction. The basalt lava was highly crystallised and degassed prior to extrusion. Together with a low effusion temperature and rapid convective heat loss in a submarine setting, these properties significantly affected the magma rheology (increased the viscosity and shear strength) and influenced the final dome-like form of the extrusion. Conversely, high heat retention was favoured by the degassed state of the magma (minimal undercooling), a thick breccia carapace and viscous shear heating, which helped to sustain magmatic (eruption) temperatures and enhanced the mobility of the flow. Received: 1 August 1996 / Accepted: 15 September 1997     

Detection of a new summit crater on Bezymianny Volcano lava dome: satellite and field-based thermal data

An explosive eruption occurred at the summit of Bezymianny volcano (Kamchatka Peninsula, Russia) on 11 January 2005 which was initially detected from seismic observations by the Kamchatka Volcanic Eruption Response Team (KVERT). This prompted the acquisition of 17 Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) satellite images of the volcano over the following 10 months. Visible and infrared data from ASTER revealed significant changes to the morphology of the summit lava dome, later seen with field based thermal infrared (TIR) camera surveys in August 2005. The morphology of the summit lava dome was observed to have changed from previous year’s observations and historical accounts. In August 2005 the dome contained a new crater and two small lava lobes. Stepped scarps within the new summit crater suggest a partial collapse mechanism of formation, rather than a purely explosive origin. Hot pyroclastic deposits were also observed to have pooled in the moat between the current lava dome and the 1956 crater wall. The visual and thermal data revealed a complex eruption sequence of explosion(s), viscous lava extrusion, and finally the formation of the collapse crater. Based on this sequence, the conduit could have become blocked/pressurized, which could signify the start of a new behavioural phase for the volcano and lead to the potential of larger eruptions in the future.     

Recent explosive eruptions and volcano hazards at Soputan volcano—a basalt stratovolcano in north Sulawesi, Indonesia

Soputan is a high-alumina basalt stratovolcano located in the active North Sulawesi-Sangihe Islands magmatic arc. Although immediately adjacent to the still geothermally active Quaternary Tondono Caldera, Soputan’s magmas are geochemically distinct from those of the caldera and from other magmas in the arc. Unusual for a basalt volcano, Soputan produces summit lava domes and explosive eruptions with high-altitude ash plumes and pyroclastic flows—eight explosive eruptions during the period 2003–2011. Our field observations, remote sensing, gas emission, seismic, and petrologic analyses indicate that Soputan is an open-vent-type volcano that taps basalt magma derived from the arc-mantle wedge, accumulated and fractionated in a deep-crustal reservoir and transported slowly or staged at shallow levels prior to eruption. A combination of high phenocryst content, extensive microlite crystallization and separation of a gas phase at shallow levels results in a highly viscous basalt magma and explosive eruptive style. The open-vent structure and frequent eruptions indicate that Soputan will likely erupt again in the next decade, perhaps repeatedly. Explosive eruptions in the Volcano Explosivity Index (VEI) 2–3 range and lava dome growth are most probable, with a small chance of larger VEI 4 eruptions. A rapid ramp up in seismicity preceding the recent eruptions suggests that future eruptions may have no more than a few days of seismic warning. Risk to population in the region is currently greatest for villages located on the southern and western flanks of the volcano where flow deposits are directed by topography. In addition, Soputan’s explosive eruptions produce high-altitude ash clouds that pose a risk to air traffic in the region.     

Summit eruptions of Klyuchevskoi Volcano,Kamchatka in the early 21st century, 2003–2013

This paper is concerned with eruptions, seismicity, and deformation on Klyuchevskoi Volcano during the summit eruptions of 2012–2013, with the condition of the central crater during the eruptions, and with the effect that is exerted by the height of the lava in the crater on the start of the eruptions. The recurrence of eruptions in the North Volcanic Cluster (NVC), Kamchatka showed that all the four volcanoes in the cluster (Klyuchevskoi, Tolbachik, Shiveluch, and Bezymyannyi) become active during definite phases that were identified in the 18.6-year lunar cycle. This relationship of the NVC eruptions to the active phases in the 18.6-year lunar cycle, as well as the relationship to the 11-year solar activity, showed that eruptions can be predicted, yielding long-term estimates of activity for the NVC volcanoes. The short-term prediction of volcanic eruptions requires knowledge of seismicity and deformation that occur during the precursory period and during the occurrence of eruptions. Seismic activity during the summit eruptions of 2003–2013 took place in the depth range 20–25 km during repose periods of the volcano and at depths of 0–5 km in the volcanic edifice during the eruption. One notes an almost complete absence of any earthquakes at great depths during the summit eruptions. Volcanic tremor (VT) was recorded from the time that the eruptions began and continued to occur until the end. Geodetic measurements showed that the center of the magma pressure beneath the volcano during the parasitic and summit eruptions of 1979–1989 moved in the 4–17 km depth range, while during the summit eruptions of 2003–2013 the center moved in the 15–20 km range. These changes in the depth of the center of magma pressure may have been related to evacuation from shallow magma chambers.     

Seismic characterization of the fall 2007 eruptive sequence at Bezymianny Volcano,Russia

We examine an eruptive sequence in late 2007 at Bezymianny Volcano to characterize the magmatic plumbing system and eruption-related seismicity. Earthquake locations reveal seismicity below and offset to the north of the volcano along a tectonic fault. Based on historical seismicity, the magma chamber is postulated to have a top at about 6 km depth. Minor dome explosions, large sub-plinian eruptions and dome collapses are analyzed using an automated event classification scheme. Low-frequency tremor, interpreted as gas escape, and low-frequency earthquakes are a dominant proportion of the energy released. We also examine multiplet earthquakes whose behavior during the study period changed significantly and systematically before the largest eruption, demonstrating the potential of tracking multiplets to assess changing conditions with the conduit.     

Volcanic hazards from Bezymianny- and Bandai-type eruptions

Major slope failures are a significant degradational process at volcanoes. Slope failures and associated explosive eruptions have resulted in more than 20 000 fatalities in the past 400 years; the historic record provides evidence for at least six of these events in the past century. Several historic debris avalanches exceed 1 km³ in volume. Holocene avalanches an order of magnitude larger have traveled 50–100 km from the source volcano and affected areas of 500–1500 km². Historic eruptions associated with major slope failures include those with a magmatic component (Bezymianny type) and those solely phreatic (Bandai type). The associated gravitational failures remove major segments of the volcanoes, creating massive horseshoe-shaped depressions commonly of caldera size. The paroxysmal phase of a Bezymianny-type eruption may include powerful lateral explosions and pumiceous pyroclastic flows; it is often followed by construction of lava dome or pyroclastic cone in the new crater. Bandai-type eruptions begin and end with the paroxysmal phase, during which slope failure removes a portion of the edifice. Massive volcanic landslides can also occur without related explosive eruptions, as at the Unzen volcano in 1792.The main potential hazards from these events derive from lateral blasts, the debris avalanche itself, and avalanche-induced tsunamis. Lateral blasts produced by sudden decompression of hydrothermal and/or magmatic systems can devastate areas in excess of 500km² at velocities exceeding 100 m s^–1. The ratio of area covered to distance traveled for the Mount St. Helens and Bezymianny lateral blasts exceeds that of many pyroclastic flows or surges of comparable volume. The potential for large-scale lateral blasts is likely related to the location of magma at the time of slope failure and appears highest when magma has intruded into the upper edifice, as at Mount St. Helens and Bezymianny.Debris avalanches can move faster than 100 ms^–1 and travel tens of kilometers. When not confined by valley walls, avalanches can affect wide areas beyond the volcano's flanks. Tsunamis from debris avalanches at coastal volcanoes have caused more fatalities than have the landslides themselves or associated eruptions. The probable travel distance (L) of avalanches can be estimated by considering the potential vertical drop (H). Data from a catalog of around 200 debris avalanches indicates that the H/L rations for avalanches with volumes of 0.1–1 km³ average 0.13 and range 0.09–0.18; for avalanches exceeding 1 km³, H/L ratios average 0.09 and range 0.5–0.13.Large-scale deformation of the volcanic edefice and intense local seismicity precede many slope failures and can indicate the likely failure direction and orientation of potential lateral blasts. The nature and duration of precursory activity vary widely, and the timing of slope faliure greatly affects the type of associated eruption. Bandai-type eruptions are particularly difficult to anticipate because they typically climax suddenly without precursory eruptions and may be preceded by only short periods of seismicity.     

An analysis of the earthquake time series at Mount St. Helens, Washington, in the framework of recent volcanic activity

The analysis of the earthquake time distribution at Mount St. Helens reveals a good correlation between the physical state of the volcano and statistical parameters of earthquake sequence. There are three main seismic phases in the whole 1980–1986 period. The first one precedes the main eruption of May 18, 1980. It begins with a sudden increase of the seismicity level in late March and continues with an Utsu (1961) type decay of the seismic occurrence rate, characterized by a small value of the decay coefficient, β. The second phase lasts from the cataclysmic eruption on May 18, 1980 until the continuous dome building episode in 1983 and is characterized by a very slow exponential increase of the background level of seismicity. The third phase covers the remaining part of the sample and is characterized by a stationary earthquake clustering process episodically interrupted by peaks of activity related to eruptions. The trends in seismic occurrence rate within each phase, as well as the statistical parameter variations at each transition, are analyzed and discussed in the framework of volcanic activity. This leads to the conclusion that statistical techniques may give a significant contribution in understanding changes in volcanic processes such as those at Mount St. Helens.     

Microearthquake activity associated with underground coal-mining in Buchanan County,Virginia, U.S.A.

Microearthquake activity (impulsive, transient seismic events, with durations up to several seconds at a distance of 500 m, that exhibit a coda with a shift toward lower frequencies with increasing time) was monitored for a three-month period by a single seismograph sited directly above an undergound longwall mine in the coal-mining region of Buchanan County, Virginia, U.S.A. The purpose of this investigation was to determine if precursory increases in microseismicity prior to cavings (subsidence) of overburden in the mine were present and, if so, could they be detected by surface seismographic observations. The first two recording weeks were prior to the beginning of coal removal operations at the monitored mine. A comparision of the before and after levels of microearthquake occurrence indicated a sevenfold increase to about seven seismic events/hour that was attendant with the development of the time over the level of the background, non-coal-mining period seismicity.A total of over 15,000 microearthquakes were recorded during the monitoring period, most of which occurred during the actual coal-mining operations. The workday rate exceeded 30 seismic events/hour in contrast with the non-workday rate of about seven such events/hour. Rock and coal fracturing ahead of the mine plow are believed to be the primary cause of the majority of these very small seismic events. Cavings and rockbursts (violent eruptions that propel rock debris into the mine) also contributed to the total seismic activity. It appears that cavings, some of which were large enough to be felt on ground surface, are the primary source of the non-plowing related seismicity as larger free surface areas are opened underground. Any seismic activity premonitory to cavings, however, was effectively masked by the high workday rate. Thus, the use of surface seismic monitoring, in an attempt to document any increases of localized seismicity precursory to cavings, failed in this instance.The exact location of the mine and the survey dates are not given in this paper at the request of the mine operator.     

岩石非均匀性对地震活动状态影响的数值模拟研究

地震的孕育、发生与震源位置及其岩石力学性质有着密切的联系,本文借助数值模拟方法研究了岩石均质度在岩样变形直至破裂过程中对声发射事件活动状态的影响。数值模拟结果表明,岩石非均匀性对岩样声发射活动状态的时间演化过程有着较大影响,所得结论进一步支持了分区研究强地震前兆异常的必要性。     

A sequence of seismic activity in the Kanto area precursory to the 1923 Kanto earthquake

The Kanto earthquake (M=7.9) that occurred along the Sagami Trough in the Sagami Bay on 1 September 1923 was one of the most disastrous earthquakes in Japanese history. The Kanto area includes Metropolitan Tokyo and Yokohama which are densely populated, and hence it has been a matter of great concern, from the viewpoints of earthquake prediction and disaster prevention, whether or not the 1923 Kanto earthquake was preceded by precursory seismicity. A study using the most complete lists of earthquakes catalogued recently by Utsu and the Japan Meteorological Agency reveals that seismic activity in the Kanto area was appreciably higher before and after the Kanto earthquake, and that the Kanto earthquake was preceded by a sequence of anomalous seismic activity, quiescence, and foreshocks. Such higher activity before and after the Kanto earthquake is contrasted with low seismicity during the recent 30-year period. A model is proposed to explain the precursory seismic activity, subsequent quiescence, and foreshocks for the Kanto earthquake. In the model, the transition from precursory seismic activity to quiescence is ascribed to time-dependent fracture due to stress-aided corrosion. Foreshocks are related to an acceleration of premonitory slip shortly before the mainshock slip.     

Applying Fractal Dimensions and Energy-Budget Analysis to Characterize Fracturing Processes During Magma Migration and Eruption: 2011–2012 El Hierro (Canary Islands) Submarine Eruption

The impossibility of observing magma migration inside the crust obliges us to rely on geophysical data and mathematical modelling to interpret precursors and to forecast volcanic eruptions. Of the geophysical signals that may be recorded before and during an eruption, deformation and seismicity are two of the most relevant as they are directly related to its dynamic. The final phase of the unrest episode that preceded the 2011–2012 eruption on El Hierro (Canary Islands) was characterized by local and accelerated deformation and seismic energy release indicating an increasing fracturing and a migration of the magma. Application of time varying fractal analysis to the seismic data and the characterization of the seismicity pattern and the strain and the stress rates allow us to identify different stages in the source mechanism and to infer the geometry of the path used by the magma and associated fluids to reach the Earth’s surface. The results obtained illustrate the relevance of such studies to understanding volcanic unrest and the causes that govern the initiation of volcanic eruptions.