Corentin Caudron

ObsPy: a bridge for seismology into the scientific Python ecosystem

The Python libraries NumPy and SciPy are extremely powerful tools for numerical processing and analysis well suited to a large variety of applications. We developed ObsPy (http://obspy.org), a Python library for seismology intended to facilitate the development of seismological software packages and workflows, to utilize these abilities and provide a bridge for seismology into the larger scientific Python ecosystem. Scientists in many domains who wish to convert their existing tools and applications to take advantage of a platform like the one Python provides are confronted with several hurdles such as special file formats, unknown terminology, and no suitable replacement for a non-trivial piece of software. We present an approach to implement a domain-specific time series library on top of the scientific NumPy stack. In so doing, we show a realization of an abstract internal representation of time series data permitting I/O support for a diverse collection of file formats. Then we detail the integration and repurposing of well established legacy codes, enabling them to be used in modern workflows composed in Python. Finally we present a case study on how to integrate research code into ObsPy, opening it to the broader community. While the implementations presented in this work are specific to seismology, many of the described concepts and abstractions are directly applicable to other sciences, especially to those with an emphasis on time series analysis.

On the use of remote infrasound and seismic stations to constrain the eruptive sequence and intensity for the 2014 Kelud eruption

The February 2014 eruption of Kelud volcano (Indonesia) destroyed most of the instruments near it. We use remote seismic and infrasound sensors to reconstruct the eruptive sequence. The first explosions were relatively weak seismic and infrasound events. A major stratospheric ash injection occurred a few minutes later and produced long-lasting atmospheric and ground-coupled acoustic waves that were detected as far as 11,000 km by infrasound sensors and up to 2300 km away on seismometers. A seismic event followed similar to 12 minutes later and was recorded 7000 km away by seismometers. We estimate a volcanic intensity around 10.9, placing the 2014 Kelud eruption between the 1980 Mount St. Helens and 1991 Pinatubo eruptions intensities. We demonstrate how remote infrasound and seismic sensors are critical for the early detection of volcanic explosions, and how they can help to constrain and understand eruptive sequences.

Stress and mass changes at a “wet” volcano: Example during the 2011–2012 volcanic unrest at Kawah Ijen volcano (Indonesia)

Since 2010, Kawah Ijen volcano has been equipped with seismometers, and its extremely acid volcanic lake has been monitored using temperature and leveling sensors, providing unprecedented time resolution of multiparametric data for an acidic volcanic lake. The nature of stress and mass changes of the volcano is studied by combining seismic analyses and volcanic lake measurements that were made during the strongest unrest ever recorded by the seismic network at Kawah Ijen. The distal VT earthquake swarm that occurred in May 2011 was the precursor of volcanic unrest in October 2011 that caused an increase in shallow earthquakes. The proximal VT earthquakes opened pathways for fluids to ascend by increasing the permeability of the rock matrix. The following months were characterized by two periods of strong heat and mass discharge into the lake and by the initiation of monochromatic tremor (MT) activity when steam/gases interacted with shallow portions of the aquifer. Significant seismic velocity variations, concurrent with water level rises in which water contained a large amount of steam/gas, were associated with the crises, that caused an although the unrest did not affect the shallow hydrothermal system at a large scale. Whereas shallow VT earthquakes likely reflect a magmatic intrusion, MT and relative seismic velocity changes are clearly associated with shallow hydrothermal processes. These results will facilitate the forecast of future crises.

Kawah Ijen volcanic activity: a review

Kawah Ijen is a composite volcano located at the easternmost part of Java island in Indonesia and hosts the largest natural acidic lake in the world. We have gathered all available historical reports on Kawah Ijen’s activity since 1770 with the purpose of reviewing the temporal evolution of its activity. Most of these observations and studies have been conducted from a geochemical perspective and in punctuated scientific campaigns. Starting in 1991, the seismic activity and a set of volcanic lake parameters began to be weekly available. We present a database of those measurements that, combined with historical reports, allow us to review each eruption/unrest that occurred during the last two centuries. As of 2010, the volcanic activity is monitored by a new multi-disciplinary network, including digital seismic stations, and lake level and temperature measurements. This detailed monitoring provides an opportunity for better classifying seismic events and forecasting volcanic unrest at Kawah Ijen, but only with the understanding of the characteristics of this volcanic system gained from the historical review presented here.

New insights into Kawah Ijen's volcanic system from the wet volcano workshop experiment

Abstract Volcanoes with crater lakes and/or extensive hydrothermal systems pose significant challenges with respect to monitoring and forecasting eruptions, but they also provide new opportunities to enhance our understanding of magmatic–hydrothermal processes. Their lakes and hydrothermal systems serve as reservoirs for magmatic heat and fluid emissions, filtering and delaying the surface expressions of magmatic unrest and eruption, yet they also enable sampling and monitoring of geochemical tracers. Here, we describe the outcomes of a highly focused international experimental campaign and workshop carried out at Kawah Ijen volcano, Indonesia, in September 2014, designed to answer fundamental questions about how to improve monitoring and eruption forecasting at wet volcanoes.

A new Bayesian Event Tree tool to track and quantify volcanic unrest and its application to Kawah Ijen volcano

Although most of volcanic hazard studies focus on magmatic eruptions, volcanic hazardous events can also occur when no migration of magma can be recognized. Examples are tectonic and hydrothermal unrest that may lead to phreatic eruptions. Recent events (e.g., Ontake eruption on September 2014) have demonstrated that phreatic eruptions are still hard to forecast, despite being potentially very hazardous. For these reasons, it is of paramount importance to identify indicators that define the condition of nonmagmatic unrest, in particular for hydrothermal systems. Often, this type of unrest is driven by movement of fluids, requiring alternative monitoring setups, beyond the classical seismic-geodetic-geochemical architectures. Here we present a new version of the probabilistic BET (Bayesian Event Tree) model, specifically developed to include the forecasting of nonmagmatic unrest and related hazards. The structure of the new event tree differs from the previous schemes by adding a specific branch to detail nonmagmatic unrest outcomes. A further goal of this work consists in providing a user-friendly, open-access, and straightforward tool to handle the probabilistic forecast and visualize the results as possible support during a volcanic crisis. The new event tree and tool are here applied to Kawah Ijen stratovolcano, Indonesia, as exemplificative application. In particular, the tool is set on the basis of monitoring data for the learning period 2000-2010, and is then blindly applied to the test period 2010-2012, during which significant unrest phases occurred.

Relative seismic velocity variations correlate with deformation at Kīlauea volcano

Seismic velocity changes correlate with deformation at Kīlauea volcano, advancing noise interferometry as a monitoring tool. Seismic noise interferometry allows the continuous and real-time measurement of relative seismic velocity through a volcanic edifice. Because seismic velocity is sensitive to the pressurization state of the system, this method is an exciting new monitoring tool at active volcanoes. Despite the potential of this tool, no studies have yet comprehensively compared velocity to other geophysical observables on a short-term time scale at a volcano over a significant length of time. We use volcanic tremor (~0.3 to 1.0 Hz) at Kīlauea as a passive source for interferometry to measure relative velocity changes with time. By cross-correlating the vertical component of day-long seismic records between ~230 station pairs, we extract coherent and temporally consistent coda wave signals with time lags of up to 120 s. Our resulting time series of relative velocity shows a remarkable correlation between relative velocity and the radial tilt record measured at Kīlauea summit, consistently correlating on a time scale of days to weeks for almost the entire study period (June 2011 to November 2015). As the summit continually deforms in deflation-inflation events, the velocity decreases and increases, respectively. Modeling of strain at Kīlauea suggests that, during inflation of the shallow magma reservoir (1 to 2 km below the surface), most of the edifice is dominated by compression—hence closing cracks and producing faster velocities—and vice versa. The excellent correlation between relative velocity and deformation in this study provides an opportunity to understand better the mechanisms causing seismic velocity changes at volcanoes, and therefore realize the potential of passive interferometry as a monitoring tool.

New insights into the Kawah Ijen hydrothermal system from geophysical data

Abstract The magmatic–hydrothermal system of Kawah Ijen volcano is one of the most exotic on Earth, featuring the largest acidic lake on the planet, a hyper-acidic river and a passively degassing silicic dome. While previous studies have mostly described this unique system from a geochemical perspective, to date there has been no comprehensive geophysical investigation of the system. In our study, we surveyed the lake using a thermocouple, a thermal camera, an echo sounder and CO2 sensors. Furthermore, we gained insights into the hydrogeological structures by combining self-potential surveys with ground and water temperatures. Our results show that the hydrothermal system is self-sealed within the upper edifice and releases pressurized gas only through the active crater. We also show that the extensive hydrological system is formed by not one but three aquifers: a south aquifer that seems to be completely isolated, a west aquifer that sustains the acidic upper springs, and an east aquifer that is the main source of fresh water for the lake. In contrast with previous research, we emphasize the heterogeneity of the acidic lake, illustrated by intense subaqueous degassing. These findings provide new insights into this unique, hazardous hydrothermal system, which may eventually improve the existing monitoring system.

Infrasound and seismic detections associated with the 7 September 2015 Bangkok fireball

A bright fireball was reported at 01:43:35 UTC on September 7, 2015 at a height of