Subglacial seismicity presents the opportunity to monitor inaccessible glacial beds at the epicentral location and time. Glaciers can be underlain by rock or till, a first order control on bed mechanics. Velocity-weakening, necessary for unstable slip, has been shown for each bed type, but is much stronger and evolves over more than an order of magnitude longer distances for till beds. Utilizing a de-stiffened double direct shear apparatus, we found conditions for instability at freezing temperatures and high slip rates for both bed types. During stick–slip stress-drops, we recorded acoustic emissions with piezoelectric transducers frozen into the ice. The two populations of event waveforms appear visually similar and overlap in their statistical features. We implemented a suite of supervised machine learning algorithms to classify the bed type of recorded waveforms and spectra, with prediction accuracy between 65–80%. The Random Forest Classifier is interpretable, showing the importance of initial oscillation peaks and higher frequency energy. Till beds have generally higher friction and resulting stress-drops, with more impulsive first arrivals and more high frequency content compared to rock emissions, but rock beds can produce many till-like events. Seismic signatures could enhance interpretation of bed conditions and mechanics from subglacial seismicity.

The Netherlands is in the midst of an energy transition with hydrocarbon production gradually declining, whereas the role of sustainable energy technologies is on the rise. One of these technologies is geothermal energy production from porous reservoirs at 1.5–3 km depth. As the number of geothermal projects increases, there is a growing concern that felt and/or damaging induced seismic events could occur as a result of geothermal operations. Over the last two decades, such events have occurred in the Netherlands due to gas production, notably in the Groningen gas field. However, the occurrence of felt events is limited to hydrocarbon fields in certain regions or reservoirs. Understanding where and for which plays these events are observed helps to estimate seismogenic potential for geothermal operations and other sustainable subsurface activities. Here, we summarise and review the main similarities and differences in terms of geological and geomechanical characteristics between the hydrocarbon and geothermal plays in the Netherlands, and we consider the differences in pressure and temperature changes. By doing so, we provide better insights into the factors that could play a role for fault reactivation and induced seismicity, and how these differ for hydrocarbon production and geothermal operations in the Netherlands. The review shows that geological characteristics for most geothermal target reservoirs are similar to those of hydrocarbon, albeit geothermal projects so far target higher porosity rocks than hydrocarbon reservoirs. On the other hand, pressure and temperature changes are very different, with significant depletion for hydrocarbon fields vs significant cooling around geothermal injection wells. The different operations result not only in different expected stress change magnitudes but also in a distinct spatio-temporal stress build-up on faults, which has implications for seismogenic potential and monitoring of these different operations.

Glacial hydraulic tremor (GHT) can be monitored to observe changes in location and distribution of water flow beneath glacial ice, allowing the spatiotemporal evolution of subglacial hydrology to be studied continuously and remotely. We use frequency-dependent polarization analysis (FDPA) to classify types of GHT and assess its spatio-temporal extent beneath Rhonegletscher, Switzerland, in a continuous seismic record through 2018 and 2019 at three ice-proximal bedrock-based seismometers. We determine the frequency bands composing the GHT and calculate back azimuth angles pointing to a previously known subglacial channel. We also inspect the relationship between GHT seismic power and water discharge from the glacier to compare daily and seasonal shifts with the observed GHT. We observed the seasonal shift from a distributed system to a channelized system, and our dataset allowed comparison of channel locations within and across seasons, with implications for sediment evacuation and bed erosion. The successful use of this method to assess GHT previously on Taku glacier (the methods of which this project follows) and now Rhonegletscher shows that existing ice-proximal passive seismic installations can be used to easily and continuously monitor subglacial hydraulic activity.

Khovsgol Lake, a seasonal ice-covered lake located in boreal Siberia, plays important economic, transportation and agricultural roles for local residents. On 11 January 2021, an Mw 6.7 earthquake struck the center of this lake, providing a unique opportunity to better understand lake ice movement, ridge distribution, and hydraulic and hydrodynamic processes, as a result of the mainshock. We use a pixel-tracking method on both Sentinel-1 synthetic aperture radar and Sentinel-2 optical remote-sensing data covering this lake before, during and after the mainshock to obtain lake-ice displacement. Combining these measurements with different viewing geometries, we determine the locations of pressure ridges in each period and derive their relative 3D displacement maps with a precision of <2.2 m. Our results indicate a large E-W displacement of up to 10.3 m and a vertical displacement of up to 4.6 m during the coseismic period, which has a distinctly different magnitude and pattern of displacement compared to the pre- and post-event periods. Given the ice motion pattern and magnitude, we speculate that the unique ice displacement of Khovsgol Lake may have been caused by a rare seismic seiche after the mainshock.

Glacial environments exhibit temporally variable microseismicity. To investigate how microseismicity influences event detection, we implement two noise-adaptive digital power detectors to process seismic data from Taylor Glacier, Antarctica. We add scaled icequake waveforms to the original data stream, run detectors on the hybrid data stream to estimate reliable detection magnitudes and compare analytical magnitudes predicted from an ice crack source model. We find that detection capability is influenced by environmental microseismicity for seismic events with source size comparable to thermal penetration depths. When event counts and minimum detectable event sizes change in the same direction (i.e. increase in event counts and minimum detectable event size), we interpret measured seismicity changes as ‘true’ seismicity changes rather than as changes in detection. Generally, one detector (two degree of freedom (2dof)) outperforms the other: it identifies more events, a more prominent summertime diurnal signal and maintains a higher detection capability. We conclude that real physical processes are responsible for the summertime diurnal inter-detector difference. One detector (3dof) identifies this process as environmental microseismicity; the other detector (2dof) identifies it as elevated waveform activity. Our analysis provides an example for minimizing detection biases and estimating source sizes when interpreting temporal seismicity patterns to better infer glacial seismogenic processes.

Southern Iceland is one of the main outlets of the Icelandic ice sheet and is subject to seismicity of both tectonic and volcanic origins along the South Iceland Seismic Zone (SISZ). A sedimentary complex spanning Marine Isotopic Stage 6 (MIS 6) to the present includes evidence of both activities. It includes a continuous sedimentary record since the Eemian interglacial period, controlled by a rapid deglaciation, followed by two marine glacioisostasy-forced transgressions, separated by a regression phase connected to an intra-MIS 5e glacial advance. This record has been constrained by tephrostratigraphy and dating. Analysis of this record has provided better insights into the interconnectedness of hydrology and volcanic and tectonic activity during deglaciations and glaciations. Low-intensity earthquakes recurrently affected the water-laid sedimentation during the early stages of unloading, accompanying rifting events, dyke injection, and fault reactivations. During full interglacial periods, earthquakes were significantly less frequent but of higher magnitude along the SISZ, due to stress accumulation, favored by low groundwater levels and more limited magma production. Occurrence of volcanism and seismicity in Iceland is commonly related to rifting events. Subglacial volcanic events seem moreover to have been related to stress unlocking related to limited or full unloading/deglaciation events. Major eruptions were mostly located at the melting margin of the ice sheet.

We recorded the ice motion and icequakes on the floating part of Langhovde Glacier in East Antarctica to better understand the dynamic behavior of ice shelves and floating tongues. Diurnal and semi-diurnal variations in ice motion and seismicity were simultaneously observed at all four global navigation satellite system and three seismic stations over 2 weeks. The short-term along-flow ice motion is explained by the elastic response of the glacier to ocean tide-induced hydrostatic stress variations, which decayed at a rate of 0.8 km−1 toward the grounding line. We observed a large number of icequakes during mid-rising and high tides that covered a broad frequency range and formed two major groups of events centered at 10 and 120 Hz, respectively. The hourly occurrence rates were ~500 events h−1, with the observed seismicity consistent with fracture due to floating tongue bending. We also observed minor secondary peaks at high ice speeds, which could reflect surface cracking due to stretching or basal friction. Our observation demonstrates that tidal-modulation was the main factor to fracture the floating tongue of Langhovde Glacier.

The evolution of glaciers and ice sheets depends on processes in the subglacial environment. Shear seismicity along the ice–bed interface provides a window into these processes. Such seismicity requires a rapid loss of strength that is typically ascribed to rate-weakening friction, i.e., decreasing friction with sliding or sliding rate. Many friction experiments have investigated glacial materials at the temperate conditions typical of fast flowing glacier beds. To our knowledge, however, these studies have all found rate-strengthening friction. Here, we investigate the possibility that rate-weakening rock-on-rock friction between sediments frozen to the bottom of the glacier and the underlying water-saturated sediments or bedrock may be responsible for subglacial shear seismicity along temperate glacier beds. We test this ‘entrainment-seismicity hypothesis’ using targeted laboratory experiments and simple models of glacier sliding, seismicity and sediment entrainment. These models suggest that sediment entrainment may be a necessary but not sufficient condition for the occurrence of basal shear seismicity. We propose that stagnation at the Whillans Ice Stream, West Antarctica may be caused by the growth of a frozen fringe of entrained sediment in the ice stream margins. Our results suggest that basal shear seismicity may indicate geomorphic activity.

Geophysicists that deploy seismic sensors in ablation zones of glaciers and ice sheets to record glaciogenic signatures can confront recording challenges caused by instrument melt-out or tilt. These challenges often require installing sensors in boreholes to delay melt-out, or securing sensors to structures that improve coupling. We show that some of these structures that were buried near a moulin at a snow-free site in the ablation zone of the Western Greenland Ice Sheet resonated as they became exposed, and caused their geophones to record temporally evolving, narrowband signals that mimic features of glaciogenic sources like moulin tremor. We quantify these artifacts with a mechanical model that shows instruments undergo structural resonance as they melt-out, at exposure rates that we predict from an ablation model (RACMO). These models reproduce general spectral features in our data, and enable us to estimate what instrument exposure reduces ice-to-sensor coupling enough to prevent icequake detection. Last, we use our resonance data to quantitatively measure how narrowband signals that originate from either artificial or glaciogenic sources will reduce the ability of certain waveform detectors (correlators) to capture transient seismic events, even if sensors remain coupled.

We conducted a 9-d seismic experiment in October 2015 at Laohugou Glacier No. 12. We identified microseismic signals using the short-term/long-term average trigger algorithm at four stations and classified them as long and short-duration events based on waveform, frequency, duration and magnitude characteristics. Both categories show systematical diurnal trends. The long-duration events are low-frequency tremor-like events that mainly occurred during the daytime with only several events per day. These events lasted tens of seconds to tens of minutes and are likely related to resonance of daytime meltwater. The dominant short-duration events mostly occurred during the night time with a peak occurrence frequency of ~360 h−1. Their short-duration (<0.2 s), high frequency (20–100 Hz) and dominance of Rayleigh waves are typical of events for near-surface crack opening. A strong negative correlation between the hourly event number and temperature change rate suggests that the occurrence of night-time events is controlled by the rate of night-time cooling. We estimated the near-surface tensile stress due to thermal contraction at night to be tens of kilopascals, which is enough to induce opening of surface cracks with pre-existing local stress concentrations, although we cannot exclude the effect of refreezing of meltwater produced during the day.

A distinct spatial relationship between surface faulting, magmatic intrusions and volcanic activity exists in the Aegean continental crust. In this paper, we provide detailed structural observations from key onshore areas, as well as compilations of lineament maps and earthquake locations with focal plane solutions from offshore areas to support such a relationship. Although pluton emplacement was associated with low-angle extensional detachments, the NNE- to NE-trending strike-slip faults also played an important role in localizing the Middle Miocene plutonism, providing ready pathways to deeper magma batches, and controlling the late-stage emplacement and deformation of granites in the upper crust. Additionally, the linear arrangements of volcanic centres, from the Quaternary volcanoes along the active South Aegean Volcanic Arc, are controlled primarily by NE-trending faults and secondarily by NW-trending faults. These volcanic features are located at several extensional settings, which are associated with the main NE-trending faults, such as (i) in the extensional steps or relay zones between strike-slip and oblique-normal fault segments, (ii) at the overlap zones between oblique-normal faults associated with an extensional strike-slip duplex and (iii) at the tip zone of a NE-trending divergent dextral strike-slip zone. The NE trend of volcano-tectonic features, such as volcanic cone alignments, concentration of eruptive centres, hydrothermal activity and fractures, indicates the significant role of tectonics in controlling fluid and magma pathways in the Aegean upper crust. Furthermore, microseismicity and focal mechanisms of earthquakes in the area confirm the activity and present kinematics of these NE- trending faults.

Local earthquakes recorded at the King Sejong station (62° 13′31″S, 58° 47′07″W) from 1995–96 have been analysed to study the seismicity and tectonics around the northern Antarctic Peninsula. The nature of shallow-focused normal fault earthquakes along the South Shetland Platform is still unclear. Dominant normal fault earthquakes and minor strike-slip earthquakes in the Eastern Bransfield Basin suggest 1) ongoing extension, and 2) transtensional stress transmitted from the Antarctic–Scotia transform boundaries, the South Scotia Ridge and the Shackleton Fracture Zone. A lack of seismicity in the Central Bransfield Basin supports that active seismicity in the Eastern Bransfield Basin is not a result of subduction along the South Shetland Trench. Shallow focused earthquakes have been observed along the NW–SE trending gravity low line between the Central and the Eastern Bransfield Basins that approximately coincides with the landward projection of a fracture zone in the former Phoenix Plate.