Seismicity in the Northern Rhine Area (1995–2018)

Since the mid-1990s, the local seismic network of the University of Cologne has produced digital seismograms. The data all underwent a daily routine processing. For this study, we re-processed data of almost a quarter century of seismicity in the Northern Rhine Area (NRA), including the Lower Rhine Embayment (LRE) and the Eifel Mountain region (EMR). This effort included refined discrimination between tectonic earthquakes, mine-induced events, and quarry blasts. While routine processing comprised the determination of local magnitude ML, in the course of this study, source spectra-based estimates for moment magnitude MW for 1332 earthquakes were calculated. The resulting relation between ML and MW agrees well with the theory of an ML ∝ 1.5 MW dependency at magnitudes below 3. By applying Gutenberg-Richter relation, the b-value for ML was less (0.82) than MW (1.03). Fault plane solutions for 66 earthquakes confirm the previously published N118° E direction of maximum horizontal stress in the NRA. Comparison of the seismicity with recently published Global Positioning System–based deformation data of the crust shows that the largest seismic activity during the observation period in the LRE occurred in the region with the highest dilatation rates. The stress directions agree well with the trend of major faults, and declining seismicity from south to north correlates with decreasing strain rates. In the EMR, earthquakes concentrate at the fringes of the area with corresponding the largest uplift.


Introduction
The Northern Rhine Area (NRA) mainly comprises the region between the rivers Mosel, Rhine, and Maas ( Fig. 1). Observation of its moderate seismicity started around 1905 with a station in the city of Aachen (Haußmann 1907). However, the main purpose of this station was to discriminate between natural tectonic earthquakes and mine tremors as the mine owners wanted to avoid reimbursements for damages caused by natural earthquakes (Hinzen 2011). In 1906 Haußmann, a mine surveyor and professor at the Technical University in Aachen, ordered a 1.3-ton Wiechert seismometer from the Spindler and Heuer company in Göttingen. Because he could not go himself, he sent his assistant Ludger Mintrop to Emil Wiechert, the first professor of geophysics at the Göttingen University, to receive instructions for maintaining the seismometer. Consequently, the interests of Mintrop, now regarded as the "father of exploration seismology," shifted from surveying to seismology (Haußmann 1907). The original station located in one of the mines was moved to a surface location because of trouble with high subsurface humidity that affected the recording paper. In 1908, Mintrop (1909aMintrop ( , 1909b moved to Bochum and installed a second seismic station. These two remained the only local stations operating until sometime during WW2. After the war, the Euskirchen earthquake of 1951 (Berg 1953), of magnitude 5.7 with damage of intensity VIII, triggered the interest of Martin Schwarzbach. At that time, Schwarzbach was director of the Geological Institute of the Universität zu Köln. He was eventually well known for his pioneering work on geology and climate (Schwarzbach 1963) but had no background as a seismologist. After applying for help from Wilhelm Hiller from Stuttgart, Schwarzbach established a local station in Bensberg (BNS), located 12 km east of the center of Cologne the station rested on Devonian hard rock (Fig. 1). This provided a more favorable place for seismic recording than the university campus near downtown Cologne on the eastern rim of the LRE above 360 m of Tertiary and Quaternary sediments (Ahorner 2008). Ludwig Ahorner expanded the one station enterprise into a network in the 1970s, referred to as BENS network. Between 1995 and 2000, all BENS network stations were converted to digital recording, and in 2006, the number of permanent stations reached 43 with 20 accelerometer stations.
Several scientific projects in the past explored different aspects of seismology in the NRA between Rhine, Maas, and Mosel for which data from permanent seismic stations produced a valuable database. From 1976 to 1982, the "Plateau Uplift" project employed an interdisciplinary approach to better understand the vertical movements of the Rhenish Shield still in the era before Global Positioning System (GPS) was available (Fuchs et al. 1983;Ahorner 1983). In the frame of the DECORP project in the 1980s, deep sounding seismic reflection profiles were measured throughout the Rhenish Massif (Meissner and Bortfeld 1990;DECORP 1991), the results improving the understanding of crustal structure and distribution of seismic velocities. The PALEOSIS project (Camelbeeck and Meghraoui 1996;Camelbeeck et al. 2001) revealed evidence for the first time of strong surface rupturing earthquakes during the Holocene in the LRE. The Eifel Plume Project, a large seismic tomography survey in 1997/1998, covered an area of approximately 400 × 250 km 2 centered on the Eifel volcanic fields with more than 200 permanent and temporary stations (Ritter et al. 2001). Results from the study strengthened the hypothesis of the existence of a mantle plume below the Eifel. Previously, in a tomographic study down to 70 km, Braun and Berckhemer (1993) had found an upper crustal low-velocity body under the Vogelsberg volcano. As part of the AGRIP-PINA project, array techniques to use ambient noise for the exploration of soft sediments were advanced with data measured in the LRE (Scherbaum et al. 2003). Reamer and Hinzen (2004) transferred about 96,000 hand-written phase readings of earthquakes recorded between 1975 and 1995 with the BENS network into digital format and used the database to construct an optimized 1D velocity model and calibrate parameters for M L determination. Beginning in 1999, the routinely processed data of the BENS network were reported online (http://www.seismo.uni-koeln.de, last accessed August 2020); through this link, the actual earthquake catalog is accessible.
In this paper, we report on the reprocessing of earthquakes recorded between 1995 and 2018, which comprise all events of the BENS network for which digital seismograms are available until the end of the observation period. The reprocessing of earthquakes in the study area (50°to 52°N and 6.5°to 8.5°E) included (1) control and if necessary re-picking of phases, (2) relocation, (3) re-discrimination of blasts and mine tremors falsely registered as tectonic earthquakes, and (4) for the first time the determination of moment magnitude M W for all earthquake records where the signal-to-noise ratio (SNR) was sufficiently high providing the data necessary to (5) update the magnitude frequency relation for the study area. (6) Fault plane solutions were determined for events with more than ten unambiguous polarity readings. (7) The seismicity results are then further discussed with respect to GPS-based data of crustal deformation which recently became available through work by Kreemer et al. (2020).

The network
At the beginning of the observation period (1995), the BENS network consisted of seven short-period analog (FM tape) or semi-digital (PCM tape) recording stations. With the exception of one station (JUL) on the sediments in the LRE and inside the compound of an experimental nuclear facility, all other stations were located on the Rhenish Shield (Fig. 1). As noted, after 1995, all stations were converted to digital recording with AD converters originally designed by the Royal Observatory of Belgium (Snissaert 1992). Subsequent years saw the addition of several new short-period stations and one broadband station (DREG). In 2001, a dozen stations equipped with short-period sensors for surveying the open lignite pit mines west of Cologne were integrated into the BENS network. The location of these "mining" stations changed from time to time, due to the advancing work face of the mines. In 2001, the recording hardware was updated to commercial 24-bit AD converters and industry-standard PCs with continuous recording in Nordic format (Havskov et al. 2020) which is still the status of the stations monitoring at the time of this contribution. This format was chosen for the acquisition software, because all routine processing was done by the SeisAn software package (http://seis.geus. net/software/seisan/, last accessed July 2020). Seismic event timing is controlled by DCF77 clocks. The network reached its current size of 43 stations in 2006 when 20 accelerometer stations were added (Hinzen and Fleischer 2007). Most of these strong-motion stations are in the free field placed on the soft sediments of the LRE, with changing sediment thicknesses from a few decameter to 1.3 km and in close vicinity to the active faults. These stations provided good records of events with magnitudes down to M L 1 or even lower in case of small hypocentral distances. Four accelerometer stations  (Farr et al. 2007) and the simplified geologic information is based on the 1:200,000 geologic map of Germany (Zitzmann 2003) (Hinzen et al. 2012;Hinzen 2014).
The increase of the number of active stations after 2006 to 43 reduced the median inter-station distance within the network from 66 km (in 1995) to 44 km. Station locations and parameters are available at http://www.seismo.uni-koeln.de/station/netz.htm (last accessed July 2020). In addition to the BENS network (University of Cologne 2016), data from stations of neighboring networks are shown in Fig. 1. These include the Royal Observatory of Belgium (ORB) (Belgium 1985), the Koninklijk Nederlands Meteorologisch Instituut (KNMI 1993), Ruhr-Universität Bochum (RUB Germany 2007), Federal Institute for Geosciences and Natural Resources (GRSN 1976), Erdbebendienst Südwest (LED and LER), and the Geologischer Dienst NRW (GD NRW). Data from these stations were regularly utilized to locate earthquakes; however, magnitudes in the following are based only on seismograms from the BENS network. Since 2003, regular meetings twice a year with colleagues from ORB and KNMI helped to coordinate and improve cross-border seismic activities in the Rhine-Maas area and thereby continuing the initiative of an EU project "Rapid Transfrontier Seismic Data Exchange Network" which was coordinated by the British Geological Survey between 1994 and 1997.

Reprocessing
Data reprocessing included visual inspection of all phase picks, control and eventually re-picking of arrivals with unusual residuals, re-evaluation of event type (tectonic earthquake, mine-induced, explosion), re-localization, determination of M L , and (with sufficient SNR) determination of M W from distance-corrected displacement spectra of P-and/or S-phases. M W had only been previously determined for 39 selected earthquakes with magnitudes larger than 2 which occurred between 1975 and 2001 (Reamer and Hinzen 2004). The distance correction we used for the M L determination in this study is the one for the NRA given by Reamer and Hinzen (2004): where A is half the peak to peak Wood Anderson maximum trace amplitude with an amplification of 2080, and R the hypocentral distance.

Statistics
The raw catalog contained 14,782 events of which 7778 were categorized as tectonic or mine-induced earthquakes within the study area between 50°N and 52°N and between 5.5°E and 8.5°E. After reprocessing, 3330 tectonic earthquakes remained in the list and M W could be determined for 1332 of them. (For 16 earthquakes in 1995, no magnitude could be determined.) The numerous non-tectonic events include mine-induced events from the deep coal mines in the Ruhr district (e.g., Casten and Cete 1980;Hinzen 1982;Gibowicz et al. 1990;Bischoff et al. 2010), open-pit lignite mines west of Cologne (Ahorner and Schaefer 2002), and quarry blasts. Not all detected quarry blasts were included in the daily routine processing; only examples of special interest were analyzed, based on the subjective decision of the seismologist on duty. The number and the distribution of quarry blasts in the study area were examined by Hinzen and Pietsch (2000) who estimated at the time of the study about 21,000 blasts per year (80 per working day) were fired, roughly 200 times the number of tectonic earthquakes. The sheer volume posed a particular challenge for the discrimination in routine processing. Figure 2 shows the number of earthquakes per year for the 24-year period starting in 1995. The effect of densification of the network with the lowering of the detection threshold between 1995 and 2006 is indicated by the increase of detected tectonic events from 30 to more than 100 per year. Since 2010, the yearly number is around 200 earthquakes, with the increase primarily due to smaller earthquakes with magnitudes below M L 1.0 being detected. The high number of almost 500 earthquakes in 2011 is due to the aftershock sequence of the 14. February event near Nassau (Hinzen 2019).
Further statistics on earthquakes occurrence is shown in Fig. 3, which breaks down the number of events per hour of the day, the minute, and the weekday. Time is local time and takes into account daylight saving time change in March and October of each year.

Earthquake map
The relocated epicenters of the dataset are shown on the map in Fig. 4; a complete list is available in the Online Resource 1. Seismicity is concentrated in a NW-SE trending stripe which follows the direction of the maximum horizontal shear stress (N118°E) as was determined earlier from fault plane solutions (Hinzen 2003).
Seismicity north of the Neuwied Basin in the Ahr Area and the Siebengebirge (Fig. 1) forms a~20-kmwide stripe trending parallel to the Rhine River valley with a sharp drop of activity towards NE and SW. On the right side of the river the activity extends north up to 50.9°N. Further north of this latitude, there are only two earthquakes in the dataset east of the Rhine River. Both occurred in 2015 within 20 min with magnitudes of M L 1.9 and 1.6 at depths of 2-3 km. There were hints from the public who submitted macroseismic questionnaires that these events might have been associated with drilling activities in the area; however, this could not be confirmed.
The main seismic activity in the LRE is located in the western part with two stripes parallel to the Erft Fault System and the Rurrand Fault. The activity continues further north in the Netherlands where the fault is designated as Peel Boundary Fault the location of the 1992 Roermond earthquake (M L 6) where seismic activity continued during the observation period.
The highest concentration of microearthquakes was observed in the Neuwied Basin south of the Laacher See volcano (compare Fig. 1). As stated in Section 4, these events include tectonic earthquakes as well as earthquakes, some with hypocenters below the crust, which can be associated with ongoing volcanic activity in the East Eifel volcanic field (Hensch et al. 2019).
For earthquakes with ten or more clearly identifiable polarity readings, fault plane solutions were made with a grid search. For 66 events, well-determined solutions were found and are listed in Online Resource 1; the mechanism for 27 of these with M L ≥ 2.5 are shown in Fig. 4.

M L /M W relation
The use of a consistent magnitude scale is essential in many sectors of seismological research. The moment   (Hanks and Kanamori 1979) has become a standard in the past decades. However, many (local) seismic networks, including BENS, do not determine M W on a routine basis as M L determination proves adequate particularly for earthquakes with magnitudes below 5. But particularly in regions with low to moderate seismicity and scarce earthquakes of magnitude 5 and above, the lower end of the magnitude frequency relation is important for determination of seismic hazard. Because many empirical relations between ground  Hinzen 2003). The graph on the righthand side of the map shows the depth distribution of the earthquakes along a north-south profile, background color indicates the P-wave velocity of the 1D model (Reamer and Hinzen 2004), and numbers below the scale give the P-wave velocity in kilometers per second at major discontinuities. Focal mechanisms are shown for 27 earthquakes with M L ≥ 2.5. Thin black lines link the focal spheres to the epicenters on the map. Black and white dots indicate the P and T axes, respectively. TEF Tegelen Fault, VIF Viersen Fault, PBF Pell Boundary Fault, WBF Western Border Faults, RRF Rurrand Fault, EFS Erft Fault System motion and magnitude are now based on M W , a reliable relation between M L and M W is essential when such ground motion models are applied. In the same way, as the distance calibration functions for M L have to be determined individually for every region (Hutton and Boore 1987), there is no universal M L /M W relation that can be applied for all networks (Deichmann 2017).
We used spectral analysis of P-and S-wave phases to determine source parameters, including seismic moment, corner frequency of the source spectrum, and the related moment magnitude. Model spectra following Brune (1970) were fitted to distance-corrected displacement spectra using the SeiSan (Havskov et al. 2020) spectral modeling tool for all traces with sufficient SNR. The main processing steps are (1) removal of DC offset, (2) application of a cos-taper with 10% of the signal length at both ends, (3) FFT, (4) restitution of the instrument response, and (5) correction for attenuation along the travel path assuming a Q 0 of 540 (Romanowicz and Mitchell 2007). To fit the model spectra to the observed data, a frequency band with suitable SNR is selected by comparison to a time window of equal length prior to the first arrival, and then a grid search of the low-frequency level of the spectrum and the corner frequency is used to find the best fit. The results of all fitting procedures were visually controlled. The vertical component was used to analyze P-phases and both horizontal components for the S-phase, if possible.
A total of 1332 M W values was determined from 16,295 spectra. Figure  Both linear and quadratic relations have been used to parametrize empirical relations between M L and M W (e.g., Grünthal and Wahlström 2003;Grünthal et al. 2009;Allmann et al. 2010;Goertz-Allmann et al. 2011;Munafò et al. 2016;Deichmann 2017). For our dataset with −0.7 ≤ M L ≤ 4.6, both types of relations fit the data equally well (Fig. 5b): , respectively, for earthquakes in Switzerland, and the dark green line by Grünthal et al. (2009) for Central Europe. The dashed blue line gives the relation determined by Reamer and Hinzen (2004) for earthquakes in the Northern Rhine Area. All relations are shown within the magnitude range for which they were determined 3.4 Gutenberg-Richter model Figure 6 shows two diagrams with the Gutenberg-Richter relation (Gutenberg and Richter 1956)

Discussion
Even though the routine processing of the data from the BENS network has been competently performed over the years, the reprocessing has shown that overall data quality could still be improved. The occurrence time statistic from Fig. 3 can indicate whether a significant number of man-made events remains in the catalog. In contrast to naturally occurring tectonic earthquakes, man-made events are not uniformly distributed in time. In particular, quarry blasts are bound to working days (Monday to Friday) and are often fired at specific and regular times often around midday and at the full hour (Hinzen and Pietsch 2000). Contamination with falsely classified quarry blasts would show skew in the histograms in Fig. 3. The clear trend of the histogram with 170 earthquakes per hour at night times and less than 100 earthquakes in the middle of the day is the reverse trend of the daily noise level. In the densely populated and highly industrialized study area (with an extensive highway and train network), there is a significant shift in the noise level affecting the frequency range of local earthquakes between 0.5 and 20 Hz; thus, the detection threshold increases from night to daytime, indicated by the boxplots of the magnitude distribution per hour overlaid in Fig. 3a. From 10 to 12 and 15 to 16 h local time, the average number of events per hour is 82. At 13 and 14 h, the number is increased to 109 and 115, respectively. This excess of about 50 events in 2 h might at first be seen as a sign of misclassified blasts. However, the M L 4.3 Nassau earthquake on 14 February 2011(Hinzen 2019) had a strong aftershock sequence and 40 of these occurred between 13 and 15 h local time. Furthermore, the number of events per minute does not indicate larger values around the full hour as it would be the case with more blasts. The largest number per day of the week was observed on Mondays with 511, but again 67 of these were aftershocks of the 2011 Nassau earthquake. Neglecting these aftershocks, the average number on a work day is 452 ± 14. The somewhat higher number of events on Sundays (491) and Saturdays (511) can be attributed to the lower noise levels on weekends.
In Fig. 5, the relation between M L and M W for the NRA is compared to results of previous studies. The slightly different slope of the linear relation of 0.691 compared to the 0.722 derived by Reamer and Hinzen (2004) for the same region can be attributed to the different covered magnitude ranges. Reamer and Hinzen (2004) used 39 events with magnitudes M L above 2.0. For most of those, only a limited number of phases could be used due to the few digital seismograms available from the study time period 1975 to 2001 compared to 1332 events for this study with many more digital seismograms available. Grünthal et al. (2009) updated the work by Grünthal and Wahlström (2003) and proposed a M L /M W relation for Central Europe based on 221 data pairs with original M W determinations mostly from the Swiss moment tensor solutions : The data set on which Eq. (4) is based mainly contained earthquakes with M W above 2 and up to 6, only six events had magnitudes below 1. As shown in This relation agrees with our linear fit (2) at M L 2, but Eq. (5) predicts slightly larger values at low magnitudes with M W = 1.58 compared to 1.45 at M L = 1.0.
By referring to the stochastic waveform models of Hanks and Boore (1984), Edwards et al. (2015) proposed to change the slope in relation (5) at M L smaller than 2 to a value of 2/3 (0.667) instead of 0.594 (Deichmann 2017): In a recent paper, Munafò et al. (2016) by applying random vibration theory showed that for small earthquakes with a M W below 4, M L is proportional to the logarithm of the seismic moment and the corresponding relationship between the magnitudes is: where C′ is a free constant individual to the region and data set. In a test with a data set from the Upper Tiber Valley (northern Apennines, Italy) of 1191 earthquakes, they determined a C′ of 1.15. Assuming a fixed slope of 2/3 for our data set of 1332 events in the NRA, the best fit for C′ is 0.802 (Fig. 5): and predicts M W = 1.47 at M L = 1.0 compared to 1.45 with relations (2) and (3). Deichmann (2017) did a detailed study using model calculations and empirical analysis and concluded by referring to Hanks and Boore (1984), Edwards et al. (2010, Edwards et al. 2015, and Munafò et al. (2016) that in an attenuating medium for small earthquakes below some magnitude threshold, the 2/3 slope to be almost unrefutable. The question of what the magnitude threshold is below which the 2/3 slope applies is not easy to answer. Deichmann (2017) found a gradual transition even in a rather homogeneous dataset of induced events from a limited source volume, all recorded at the same borehole station, and concluded the magnitude threshold is clearly a function of the degree of attenuation along the travel path. For a dataset like ours which averages the M L /M W relation over a comparatively large area of epicenters and uses data from some 40 stations, Deichmann (2017) places the threshold somewhere between M W 2 and 4, depending on the attenuation conditions. As our dataset contains only 21 events with M L above 3.0 but 929 with M L between 0.5 and 2.0, an unambiguous threshold cannot be stated; however, below M L 3.0, the application of Eq. (8) seems to be justified.
These differences between M L and M W for low magnitudes also affect the corresponding Gutenberg-Richter relations. For our dataset, the b-value for the M L relation (0.82) is lower than for the M W relation (1.03). Deichmann (2017) reported a similar trend for the earthquakes of the Swiss earthquake catalog, with a lower bvalue for M L compared to M W ; however, in that study, the M W values were converted from M L using the relation by Goertz-Allmann et al. (2011) with the modification by Edwards et al. (2015) (Eqs. (5) and (6)). For BENS values of M W were directly determined from the seismic records, Fig. 5b shows in addition to the directly determined M W values (black dots), the values which result when the measured M L data are converted to M W using Eq. (8) (red circles). Below M W = 2.5, the b-value of the converted data (1.02) is similar to the one obtained from measured data (1.03), but due to the influence of the few stronger earthquakes, a change in the b-value of the converted data to 1.11 occurs at magnitudes above 2.5. The largest measured M W is 4.3, but the corresponding converted value only 3.8.
As an additional check on the clustering of hypocenters, the double-difference (DD) method (Waldhauser 2001) was applied to the dataset of 51,209 differential arrival times. Limiting the search for clusters to events with a minimum of four links resulted in 824 hypocenters in 18 clusters. Figure 7 shows these earthquakes together with the rest of the dataset and the frequency depth distribution of the earthquakes in nine selected areas.
The three areas within the LRE, entitled Erft, Rur, and Roermond, contain about 25% of all earthquakes in the database. Most events in the Erft cluster occurred at depths between 6 and 15 km with a maximum depth of 26 km. The epicenters align in a pattern parallel to the Erft Fault System (Fig. 7) at a horizontal distance indicating a westward dip of the faults of about 50 to 55°. In the Rur cluster, more earthquakes are located at a shallower depth between 4 and 14 km. A clear association with one of the faults in the Rur area is not easy because at depth they can occur on either the westwarddipping Rurrand Fault or on one of the east-dipping Western Border Faults. For example, the fault plane solution of the 2002 Alsdorf earthquake (M W 4.3, depth 15.8 km) indicates a normal faulting mechanism on a 52°west-dipping plane which can be associated with the Rurrand Fault (Hinzen and Reamer 2007). The earthquakes in the Roermond cluster are in the vicinity of the strongest earthquake in the NRA in the past century with M W 5.4 ) and linked to the Peel Boundary Fault (Fig. 7).
The Voerendaal cluster in the Dutch province of Limburg shows a large percentage (46%) of hypocenters at shallow depth between 2 and 6 km. These earthquakes were part of small swarms of events between 2000 and 2001 (Dost et al. 2004). The proximity to former coal mines suggests the possibility that these shallow events may be associated in a causal sense due to the rising water table after the abandonment of these mines (Sigaran-Loria and Slob 2019) which produces heterogeneous surface displacements detected by satellite radar interferometry (Caro Cuenca et al. 2013).
The clusters in the Ardennes and in the West Eifel show hypocenters in the lower crust, almost reaching the Moho at 30 km. The high percentage of events with depth of 12 to 14 km in the West Eifel is due to two small swarms of microearthquakes in 2010 at this depth level.
More than 40% of all earthquakes are located within the Neuwied Basin cluster. While the majority of tectonic earthquakes in and around the basin has shallow hypocenters of less than 10 km, some events have been observed at depths reaching 40 km, well below the Moho discontinuity. Hensch et al. (2019), using data from temporary as well as permanent stations, showed that deep low-frequency microearthquakes below the Laacher See Volcano occurred in four distinct clusters between 10-and 40-km depth. The Laacher See is the caldera of the last eruption in the East Eifel volcanic field 12.9 Kyr ago, which was fed by a shallow magma chamber at 5-to 8-km depth, erupting a total magma volume of 6.7 km 3 (Zolitschka et al. 2000;Schmincke 2007Schmincke , 2009).
In the Taunus, the majority of hypocenters occur at depths shallower than 15 km. Some of the events form lineaments in the SW-NE direction roughly parallel to faults in the Devonian bedrock. In the Ahr cluster, two bands of epicenters stretch NE-SW and NNE-SSW across the Rhine River. The northern part of these bands includes earthquakes in the Siebengebirge, where the last volcanic activity occurred in the Miocene, and north of it on the eastern side of the Rhine River. Hypocenters in this cluster are mainly shallow, around 2 to 5 km and between 10 and 14 km.
Of the 66 fault plane solutions from the earthquakes with more than 10 unambiguous polarity readings (Online Resource 1 and Fig. 4), 30% have strike-slip character, while the others mainly conform to a normal faulting mechanism. Application of the linear inversion algorithm from Michael (1984) reveals a trend of the largest compressive stress in the NRA of N118°E. This value is identical to the maximum horizontal stress direction determined by Hinzen (2003). (Also, see the rose diagram in Fig. 4; eight of the 66 fault plane solutions from this study were also part of the database used by Hinzen 2003).
In a recently published paper, Kreemer et al. (2020) used Global Positioning System (GPS) data to robustly image vertical land motion (VLM) and horizontal strain rates over most of intraplate Europe. They found a clear spatially coherent positive VLM anomaly over a large area surrounding the Eifel volcanic fields with maximum uplift of ∼ 1 mm year −1 at the center (when corrected for glacial isostatic adjustment) and significant horizontal extension surrounded by a radial pattern of shortening. This combination strongly suggests a common dynamic cause (Kreemer et al. 2020). Their uplift model postulates a mantle plume at the bottom of the lithosphere. The location of the modeled plume agrees well with the findings by Ritter (2007) from the large-scale seismic experiment, the Eifel Plume project.
The deformation field deduced by Kreemer et al. (2020) is based on continuous GPS data of more than 2100 stations operating across Europe for up tõ 20 years. The GPS station density in the central part of Europe in the Kreemer et al. (2020) dataset includes more than 70 stations within our study area (Fig. 8), sufficient to compare the deformations with seismicity recorded in the study area over the past quarter century.
The dataset from Kreemer et al. (2020) contains horizontal and vertical velocities and deformation resolved on a 0.1°grid. From these data, we calculated the maximum shear strain: where λ 1, 2 are the eigenvalues of the strain rate tensor with componentsε xx ,ε yy , andε xy and the directions of maximum shear strain rate (Hackl et al. 2009): and the trace of the tensor corresponds to the relative variation rate of the surface area (dilatation) (Hackl et al. 2009): δ ¼ε xx þε yy ð11Þ Figure 8 shows the seismicity with respect to the vertical land movement corrected for the far-field glacial isostatic adjustment signal (Kreemer et al. 2020) and the  . Major faults are shown in black, and the edge of the Rhenish Shield is indicated by white lines. a The colored contours show the vertical land movement determined from GPS measurements by Kreemer et al. (2020) corrected for the far-field glacial isostatic adjustment signal. Triangles give the locations of volcanoes in the west and East Eifel volcanic fields (after Meyer 2013). b Same map area as in a with color-coded dilatation rates and vectors of horizontal velocities based on a model inferred from the GPS data by Kreemer et al. (2020) dilatation rate after Eq. (11) with the overlain horizontal velocities. The vertical movements show a clear maximum slightly above 1 mm/year in the center of the Eifel between the west and East Eifel volcanic fields (Meyer 2013;Schmincke 2009). In the area of the largest vertical movements (around 1 mm/year) with its center at 50.3°N and 7.0°E, few earthquakes have been observed in the past 25 years. On the other hand, at the fringe of that area, where the gradients of vertical movement are large with respect to the center of the uplifted area, clusters of earthquakes exist (i.e., along the Rhine River valley and in the Neuwied Basin, in the Hunsrück, in the Ahr region extending to the Siebengebirge, and northwest of the West Eifel volcanic field). With the exception of a small swarm of earthquakes with magnitudes below 0.5 which occurred within a few days in 2010, few events were recorded in the West Eifel volcanic field.
The largest dilatation rates are found at the western border of the LRE with a maximum of 3.7 10 −9 /year at 50.7°N and 6.4°E. As shown in Fig. 8b, horizontal velocities are practically zero at the center of the dilatation and increase with distance from it pointing radially away. In the range of the largest dilatation lies the Rur cluster of earthquakes (Fig. 7), which also includes the 2002 Alsdorf earthquake (M W 4.3), (Hinzen 2005). In addition, increased seismicity was observed here in the two years following the 1992 M W 5.4 Roermond earthquake ) which occurred 40 km northwest of the Alsdorf event (Hinzen and Reamer 2007). A similar pattern of horizontal extension in this area was noted by Camelbeeck and van Eck (1994) based on an analysis of 24 fault plane solutions.
The maximum shear strain (Eq. 8) shown in Fig. 9a together with the seismicity reveals lowest shear strain values are found in the Eifel, south of the West Eifel volcanic field, and in the northern part of the LRE. The low shear strain in the Eifel may correlate with the low seismicity in this area. In addition to the shear strain values, two (indistinguishable) directions of the maximum shear (Eq. 9) are shown in Fig. 9a. In addition to the direction at the grid points of the map (gray symbols), the directions along the major faults in the LRE and along some lineaments of the seismicity pattern within the Rhenish shield are plotted (red symbols). The directions agree well with the trend of the faults; Fig. 9 a Color-coded maximum shear strain map based on a model inferred from the GPS data by Kreemer et al. (2020). The crosses give the direction of the maximum shear strain rate, with their size scaled to shear strain rate amplitude. Gray crosses are shown at every third grid point (for better visibility). The red crosses show the directions along the Quaternary faults of the Lower Rhine Embayment and along some lineaments in the southern part of the map. Gray circles show epicenters of earthquakes . TEF Tegelen Fault, VIF Viersen Fault, PBF Pell Boundary Fault, WBF Western Border Faults, RRF Rurrand Fault, EFS Erft Fault System. b The shear strain rate amplitudes are plotted with respect to latitude along a NS trending profile at 6°E starting at the southern tip of the LRE (indicated in the map by the green line). The red dots show the number of earthquakes in 5-km-wide horizontal stripes along the profile overlapping by half the bin width e.g., they align even with the bends in the trend of the Viersen Fault and the Erft Fault system (Fig. 9a). The size of the symbols, scaled by the maximum shear strain, displays a clear decrease of this value from south to north in the LRE. To quantify this observation, in Fig.  9b, the strain rate is plotted along a south-north trending profile at 6°E (shown in Fig. 9a). In addition, Fig. 9b shows the number of earthquakes along this profile in 5km-wide east-west trending stripes which overlap by half of the bin width from the southern end of the LRE to the northern limit of the study area. The decrease of earthquake frequency in the LRE from south to north correlates well with the decrease in shear strain rate.
The directions of shear are also found to be parallel to some lineaments in the seismicity pattern in the West Eifel, and the Middle Rhine Area. In particular, the NW-SE trend of the limits of seismicity agrees well with the NW-SE shear direction (Fig. 9a).

Conclusions
The reprocessing of almost a quarter century of digital seismic data from the BENS network resulted in an updated earthquake catalog. Direct determination of moment magnitudes of 1332 small earthquakes supports the hypothesis M L ∝ 1.5 M W for magnitude below 3. An updated Gutentenberg-Richter relation with measured moment magnitudes (1.5 ≤ M W ≤ 4.3) reveals higher b-value of 1.03 than a 0.83 b-value based on local magnitudes. Fault plane solutions confirm a maximum horizontal stress at N118°E. Comparing the seismicity with recently published GPS-based deformation data shows that most earthquakes in the Lower Rhine Embayment in the past two and a half decades occurred within an area corresponding to maximum crustal dilatation. The direction of shear strain agrees well with the trend of major normal faults of the LRE. The decrease in seismic activity from south to north correlates with the decrease in shear strain rate. In the Eifel Mountain Region, the area with largest uplift rates (reaching 1 mm/year) shows few earthquakes; however, at the fringes of this area where the gradient of vertical movement is large, numerous earthquakes occurred during the observation period. In the Middle Rhine Area, lineaments in the seismicity pattern are parallel to one of the maximum shear directions.
The Northern Rhine Area has low to moderate seismicity. The 25 years of earthquake data from the BENS network do not include any spectacular damaging earthquakes. This deceptively low level of seismic hazard has led public and authorities to repeatedly underestimate the earthquake risk in this densely populated and highly industrialized area. The damaging earthquakes of the past century and those known from the historical record signify the risk. A recurrence of one of the surface rupturing earthquakes as evidenced by paleoseismic studies would result in major destruction in the study area. Therefore, continued and even intensified surveillance of the seismic activity in the NRA is essential for mitigating earthquake risk in addition to monitoring the volcanic activity in the Eifel Region.

Supplementary Information
The online version contains supplementary material available at https://doi.org/10.1007/s10950-020-09976-7. and RWE Power AG. Open access funding was provided by project DEAL.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.