Seismic monitoring is essential for understanding subsurface processes, particularly in geothermal operations where low-magnitude events can provide valuable insights into reservoir behaviour. There are two significant challenges when monitoring the seismicity in Dutch geothermal operations: (1) detecting signals from seismic events as noise levels are typically high in regions hosting geothermal operations, and (2) accurately estimating their corresponding hypocentre and uncertainty. In this study, we present a comprehensive workflow for detecting and characterising low-magnitude seismic events. Specifically, we integrated data preparation, template-matching and machine-learning-based event detection, and probabilistic hypocentre estimation. Applying this workflow to 4 months of recordings in Kwintsheul, Netherlands, we detected 65 events with coherent signals, including six weak seismic events (ML < 0.0) near a local fault and a geothermal injection well. These events suggest the presence of a recurring microseismic sequence previously unreported in the area. However, spatial uncertainties, the short monitoring period, and the limited azimuthal coverage make the nature of these events unclear. Our findings highlight the importance of improving network design and refining velocity models to reduce uncertainties in event locations and magnitudes. The proposed workflow offers a scalable solution for enhancing seismic monitoring, particularly in urban and geothermal settings.

In autumn 2017 a network of 14 broadband seismic stations was deployed at the Theistareykir high temperature geothermal field (NE Iceland). This experiment was conducted as part of the current efforts to characterize the field's main structures, and possible short and long term stress variations due to the ongoing fluid injection and extraction operations which started in autumn 2017. In this work, we use two years of continuous seismic records (October 2017–October 2019) to compute a 3D shear wave velocity model of the geothermal field and to detect possible crustal stress changes related to the injection and production activities. From phase cross-correlations of the vertical component recordings, we measure the Rayleigh wave group velocity dispersion curves to obtain 2D group velocity maps between 1 and 5 s. Subsequently, we use a neighborhood algorithm to retrieve the 3D shear wave velocity model of Theistareykir. Mainly, two sets of elongated high and low velocity anomalies can be observed oriented in a NW/WNW direction, parallel to the lineaments of the active Tjörnes fracture zone. Velocity reductions west of Ketilfjall and at Baerjafjall could indicate the location of upflow zones of the magmatic reservoir or hydrothermal system. We analyzed the temporal evolution of phase and amplitude of phase auto-correlations using the stretching technique and discuss their behavior in relation to the geothermal field operations. We notice a slightly stronger long-term velocity decrease in the reservoir region compared to outer regions. This could be related to the mass depletion in that area (higher fluid extraction compared to the water reinjection). In summary, our findings show how a monitoring network can be set up to enable a detailed imaging and monitoring of reservoir behavior in general.

In volcanic and hydrothermal systems, monitoring of mass and stress changes by continuous gravity field and ground motion records provides information for both volcanic hazard assessment and estimation of geothermal resources. We aim at a better understanding of volcanic and geothermal system processes by addressing mass changes in relation with external influences such as anthropogenic (reservoir exploitation) and natural forcing (local and regional earthquake activity, earth tides). Þeistareykir is a geothermal field located within the Northern Volcanic Zone (NVZ) of Iceland on the Mid-Atlantic Ridge. Geothermal power production started in autumn 2017. For the first time on a geothermal production field, we deployed a network of 4 continuously recording gravity meters (3 superconducting meter, iGrav and one spring gravity meter gPhone) in order to cover the spatial and the temporal changes of gravity and to detect small variations related to the geothermal power plant operation (e.g. extraction and injection). All gravity monitoring stations are equipped with additional instrumentation to measure parameters that may affect the gravity records (e.g. GNSS and hydrometeorological sensors). Additionally, we deployed a temporal seismic network consisting of 14 broadband stations to enhance the seismic activity monitoring of the permanent Icelandic network in this very active region of the NVZ. Results of this unique experiment contribute to determine reservoir properties and main structures and may also reveal details of active tectonic processes. Here, we present the instrumental setup at the site and first results of more than 24 months of continuous gravity and seismicity records.