Ice-covered Volcanoes Research Articles

(Re)Discovering the Seismicity of Antarctica: A New Seismic Catalog for the Southernmost Continent

Abstract We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on-continent temporary and long-term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (Mw) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice-covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Unrest at Cayambe Volcano revealed by SAR imagery and seismic activity after the Pedernales subduction earthquake, Ecuador (2016)

For the first time in decades, a sudden increase in seismicity has been observed and monitored at Cayambe volcano in Ecuador, in 2016. This seismic unrest, which occurred a few months after the April 2016, Mw 7.8 Pedernales subduction earthquake, has raised many questions, especially as there is no record of recent eruptions at Cayambe volcano. Here we analyze a time series of 104 images from Sentinel 1 (SAR) data spanning the period 2014–2018 using the NSBAS processing chain in order to quantify surface deformation around this potentially explosive and ice-covered volcano. We evidence a large-scale uplift reaching a maximum mean displacement rate of about 0.44 cm/yr in Line of Sight. This uplift is mainly due to a significant and sudden acceleration of the deformation pattern, focused on the SE flank starting in November 2016. We model this signal as related to magma emplacement at around 6 km depth below the summit, with a sudden volume influx of about 2.6 million m3. The inflation and surface deformation pattern is concomitant in time with two Mw 3 seismic events recorded at Cayambe volcano in November 2016 and is consistent with the location at the summit of the volcano-tectonic (VT) seismic swarm from September 2016 onwards. This VT swarm follows an earlier swarm in June 2016 located on the northern side of the volcano and closer to a western branch of the Chingual-Cosanga-Pallatanga-Puna (CCPP) fault system, that plays a major role in accommodating plate tectonic processes in South America. We thus propose that static stress changes from the Pedernales megathrust earthquake triggered magma ascent below the Cayambe volcano, through reactivation of the CCPP fault system. Finally, volcanic hazards around Cayambe ice-capped volcano, previously dormant for the last 300 years, should be reassessed in light of this recent unrest.

Persistent shallow micro-seismicity at Llaima volcano, Chile, with implications for long-term monitoring

Identifying the source mechanisms of low-frequency earthquakes at ice-covered volcanoes can be challenging due to overlapping characteristics of glacially and magmatically derived seismicity. Here we present an analysis of two months of seismic data from Llaima volcano, Chile, recorded by the permanent monitoring network in 2019. We find over 2000 repeating low-frequency events split across 82 families, the largest of which contains over 200 events. Estimated locations for the largest families indicate shallow sources directly beneath or near the edge of glaciers around the summit vent. These low-frequency earthquakes are part of an annual cycle in activity at the volcano that is strongly correlated with variations in atmospheric temperature, leading us to conclude that meltwater from ice and snow strongly affects the seismic source mechanisms related to glacier dynamics and shallow volcanic processes. The results presented here should inform future assessments of eruptive potential at Llaima volcano, as well as other ice-covered volcanoes in Chile and worldwide.

Sedimentation associated with glaciovolcanism: a review

Abstract Three discrete categories of sedimentary deposits are associated with glaciovolcanism: englacial cavity, jökulhlaup and lahar. Englacial cavity deposits are found in water-filled chambers in the lee of active glaciovolcanoes or at a locus of enhanced geothermal heat flux. The cavities provide a depocentre for the accumulation of debris, either abundant fresh juvenile debris with sparse dropstones (associated with active glaciovolcanism) or polymict basal glacial debris in which dropstones are abundant (associated with geothermal hot spots). Described examples are uncommon. By contrast, volcanogenic jökulhlaup deposits are abundant, mainly in Iceland, where they form extensive sandar sequences associated with ice- covered volcanoes. Jökulhlaups form as a result of the sudden subglacial discharge of stored meltwater. Analogous deposits known as glaciovolcanic sheet-like sequences represent the ultra-proximal lateral equivalents deposited under the ice. Glaciovolcanic lahars are associated with ice- capped volcanoes. They form as a result of explosive eruptions through relatively thin ice or following dome collapse, and they trigger mainly supraglacial rather than subglacial meltwater escape. Sediment transport and depositional processes are similar in jökulhaups and lahars and are dominated by debris flow and hyperconcentrated or supercritical flow modes during the main flood stage, although the proportions of the principal lithofacies are different.

The bedrock and tephra layer topography within the glacier filled Katla caldera, Iceland, deduced from dense RES-survey

We present results from recent low frequency radio echo-sounding (RES) campaigns over the icecovered caldera of Katla central volcano, beneath the Mýrdalsjökull ice cap, southern Iceland. The current RES-survey both partly repeats and enhances the RES-profile grid of a previous survey in 1991 with denser sounding lines and improved instruments. The RES-data, obtained in 2012–2021, include ∼760 km of 2D migrated RES-profiles covering an area of 116 km². Around 14 km² subsections of this area were surveyed with RES-profiles 20 m apart allowing 3D migration of the RES-data. Our study confirms findings from previously published bedrock mapping, including main topographic features, ice volume stored within the 100 km² caldera (45±2 km³, in autumn 2019) and maximum ice thickness (740±40 m). However, the significantly expanded level of detail and features observed in the new bedrock map reveals further evidence of a complex and eventful formation history of the caldera interior. This bedrock map is unprecedented in terms of detail for an icecovered volcano. The new RES-data allows for a unique comparison of bedrock maps obtained from RES-data with 2D and 3D migration, demonstrating the limitations of 2D migrated RES-data in areas of high topographic variability. Reflections from the 1918 Katla eruption tephra layer within the ice were detected in a much wider area within the caldera than in the 1991 RES-data. We also observe a second internal layer at 420–580 m depth within the northern part of the caldera, identified here as the tephra from the 1755 Katla eruption. The 1918 tephra layer is typically observed at 200–300 m below the glacier surface. However, the layer depth varies from ∼100 m depth at the western rim of the caldera down to 460 m depth, where geothermal activity beneath ice cauldrons melts ice from below. At the most prominent geothermal areas all ice beneath the 1918 tephra has been melted leaving the tephra at the bed. Furthermore, the obtained tephra layer maps reveal footprints of some previously unidentified geothermal areas.

Investigating potential icequakes at Llaima volcano, Chile

Glacially- and magmatically-derived seismic events have been noted to heavily overlap in characteristics, thus there exists the potential for false-alarms or missed warnings at ice-covered volcanoes. Here we present the first study to specifically target icequakes at an ice-covered volcano in Southern Chile. Two months of broadband seismic data collected at Llaima volcano in 2015 were analyzed in order to quantify, characterize, and locate glacially-derived seismic events at one of the most active ice-covered volcanoes in the region. We find over 1,000 repeating seismic events across 11 families, the largest of which contains 397 events. Approximate locations and characteristics of the largest families lead us to conclude that these events were derived from persistent stick-slip motion along the ice-rock interface at the base of a glacier near the volcano summit. These results have implications for future seismic monitoring at Llaima volcano and other ice-covered active volcanoes in the region.

Influence of ENSO and PDO on mountain glaciers in the outer tropics: case studies in Bolivia

This paper emphasize on the observational investigation of an ice-covered volcano and two glaciated mountains in the Central Andes from 1984 to 2011. Annual snowlines of the Nevado Sajama in the Cordillera Occidental and the Nevado Cololo and the Nevado Huanacuni in the Cordillera Apolobamba in Bolivia were calculated using remote sensing data. Landsat TM, Landsat ETM+, and LISS-III images taken during the end of dry season were used in this study. Changes in the highest annual snowline during May–September is used an indirect measure of the changes in the equilibrium line altitude of the glaciers in the outer tropics. We tried to understand the combined influence of the El Nino-Southern Oscillation and the Pacific Decadal Oscillation on the variations in the annual snowline altitude of the selected glaciers. Meteorological data in the form of gridded datasets were used for calculating the anomalies in precipitation and temperature during the study period. It is found that the glaciated areas were fluctuated with the occurrence of warm and cold phase of ENSO but the magnitude of the influence of ENSO is observed to be controlled by the phase changes of PDO. Snowline of the Nevado Sajama fluctuated heavily when cold and warm phases of ENSO occur during the cold and warm regimes of PDO, respectively. Nevado Cololo and Nevado Huanacuni are showing a continuous retreating trend during the same period. This clearly indicates that the changes in the Pacific SST patterns have more influence on glaciers in the Cordillera Occidental compared with those in the Cordillera Oriental of the Bolivian Andes.

Combined influence of PDO and ENSO on northern Andean glaciers: a case study on the Cotopaxi ice-covered volcano, Ecuador

This paper describes the application of remote sensing in monitoring the fluctuations in one of the mountain glaciers in the Ecuadorean Andes during the past few decades using ASTER, EO-1 ALI, Landsat MSS, TM and ETM + images. Satellite images were used to calculate the snow line altitudes (SLAs) during the period 1979–2013. Cotopaxi ice covered volcano was studied as representative of Ecuadorian glaciers in the eastern cordillera. Precipitation and air temperature data from various gauging stations within the range of 30 km from the study site and monthly discharge and water level data from a gauging station were also utilized in this study. Anomalies in precipitation and temperature were found to be slightly different in the Cotopaxi region compared to nearby Antizana in the same cordillera and Chimborazo region in the western cordillera. An attempt to correlate the El Nino—southern oscillation phenomenon with the glacier fluctuations in Ecuadorian Andes was done successfully. Cold and warm regimes of Pacific Decadal Oscillation is also considered. The calculated glacier fluctuations obtained were similar to that performed on the nearby Antizana 15 in the eastern cordillera during 1995–2002. Precipitation and temperature anomalies were similar with Antizana 15. It is evident from the research that SLAs were highly fluctuated between the period of occurrence of El Nino and La Nina events. It is also seen that the glacier fluctuations show a negative mass balance trend in during the warm regime of Pacific Decadal Oscillation during the past three decades. Glaciated areas were advanced during the La Nina events in the cold regime of PDO during 1998–2002.

The surface energy balance of an active ice-covered volcano: Villarrica Volcano, southern Chile

AbstractThe energy balance of bare snow and tephra-covered ice near the glacier equilibrium line elevation on Villarrica Volcano, southern Chile, was investigated during 2004 and 2005, combining meteorological, surface temperature and ablation measurements with energy balance modelling. A tephra thermal conductivity of 0.35 Wm–1 K–1, and a critical tephra thickness of <5mm at which ablation is reduced compared to bare snow, were obtained from field data. These low values are attributable to the highly porous lapilli particles which make up most of the surface material. Modelled melt totals in the January to March period were 4.95 m and 3.96 m water equivalent (w.e.) in 2004 and 2005, respectively, compared with ∽0.5mw.e. melt for ice buried by >0.1m tephra. Windblown tephra impurities lowered snow albedo, but increased snowmelt by only an estimated 0.28mw.e. over the same period. The net mass balance impact of supraglacial tephra at Villarrica Volcano is therefore positive, as thick ash and lapilli mantle most of the glacier ablation zones, probably reducing annual ablation by several metres w.e. In the accumulation seasons, frequent melting events were recorded with modelled daily snowmelt rates of up to 50 mmw.e.