Unprecedented quiescence in resource development area allows detection of long-lived latent seismicity

. Recent seismicity in Alberta and north-east British Columbia has been attributed to ongoing oil and gas development in the area, due to its temporal and spatial correlation. Prior to such development, the area was seismically quiescent. Here, we show evidence that latent seismicity may occur in areas where previous operations have occurred, even during a shutdown in operations. The global pandemic of COVID-19 furnished the unique opportunity to study seismicity during a long period 5 of anthropogenic quiescence. 389 events were detected within the Kiskatinaw area of British Columbia from April to August 2020, which encompasses a period with very little hydraulic fracturing operations. This reduction in operations was the result of a government imposed lockdown severely restricting the movement of people, as well as a downturn in the economic market causing industry stock prices to collapse. Except for a reduction in the seismicity rate and a lack of temporal clustering that is often characteristic of hydraulic fracturing induced sequences, the general characteristics of the observed seismicity were 10 similar to the preceding time period of active operations. During the period of relative quiescence, event magnitudes were observed between M L -1 and M L 1.2, which is consistent with previous event magnitudes in the area. Hypocenters occurred in a corridor orientated NW-SE, just as seismicity had done in previous years, and located at depths associated with the target Montney formation or shallower ( < 2.5 km). A maximum of 21 % of the detected events during lockdown may be attributable to natural seismicity, with a further 8 % potentially attributed to dynamic triggering of seismicity from teleseismic events and 6 % 15 related to ongoing salt-water disposal and a single operational well pad. However, this leaves ∼ 65 % of the seismicity detected during lockdown being unattributable to primary activation mechanisms. This seismicity is unlikely to be the result of direct pore pressure increases (as very little direct injection of ﬂuids was occurring


Introduction
The number of recorded instances of injection-induced seismicity has risen dramatically over the past decade, in part due to increased operations in hydraulic fracturing, waste-water disposal and enhanced geothermal systems around the globe, as well as enhanced monitoring meaning we are better able to detect smaller events (e.g.Atkinson et al., 2016;Ellsworth, 2013).
In western Canada, the Western Canadian Sedimentary Basin (WCSB) is the focus ofon R1 such activity, where a number of distinct resource plays are located including the Montney and the Duvernay.Despite an apparent flurry of larger magnitude seismic events associated with these operations (e.g.M L 4.5 near Fort St John, British Columbia in November 2018, Babaie-Mahani et al. (2019); Peña Castro et al. (2020); M W 4.1 near Fox Creek, Alberta in January 2016, Eyre et al. (2019b)) very few hydraulic fracturing operations (∼0.8%)(∼0.3%)R3 are actually linked to seismic activity with M W > 3 (Ghofrani and Atkinson, 2020) R3 (Atkinson et al., 2016)  R3 .
The Montney Play, ROS which is Lower-Middle Triassic in age, is formed of extensive fine-grained siliclastic units (interbedded sand, silt and mudstones), and stretches from west-central Alberta to north-east British Columbia (Eaton and Schultz, 2018;Dixon, 2000;Armitage, 1962).Over 5,600 multistage horizontal hydraulically fractured wells had been completed within the Montney by December 2018 (Nieto et al., 2018).In recent years, north-east British Columbia has experienced an increasing number of felt seismic events during active development within the Montney play.This led the BC Oil and Gas Commission to implement a special order in 2018, within the area now known as the Kiskatinaw Seismic Monitoring and Mitigation Area (KSMMA), which required operators to undertake a pre-assessment of the seismic hazard, fully inform residents in the area of upcoming operations and undertake ROS real-time seismic monitoring before, during and after completions (BC Oil and Gas Commission, 2018).Of particular importance was the introduction of the threshold for the suspensioncessation ROS of operations following events with ≥M3a M L 3.0 or above ROS within the KSMMA, which is lower than the ≥M4M L 4.0 ROS threshold that is standard elsewhere in British Columbia (e.g.Babaie-Mahani and Kao, 2020).
Prior to the introduction of oil and gas extraction, western Canada was generally seismically quiet, except for the Mackenzie Mountains and the North American plate boundary off the west coast of British Columbia (Lamontagne et al., 2008).Consequently seismicity detected within the KSMMA has been assumed to be directly related to ongoing oil and gas ROS operations due to its temporal and spatial correlation with active wells.However, there are a number of examples of seismicity thought to be related to hydraulic fracturing that generate events months after operations have ceased (e.g.Eyre et al., 2020).We call this latent seismicity i.e. seismicity that appears after an unusually long delay following a R3 primary activation processes, but which has with R3 no obvious "trigger" (e.g.enhanced pressurization at the onset of seismicity), and which cannot be explained by other sources (e.g.natural or dynamic triggering processes).
Here, we investigate seismicity generated within the KSMMA during the unprecedented period of quiescence that resulted due to the global COVID-19 pandemic.The slowingcessation ROS of operations in the area and the reduction in seismic noise as businesses shut down and people stayed indoors, compounded by a downturn in the energy market, ROS gave us the unique opportunity to study latent seismicity in an area where it would usually go undetected.This study is highly unusual considering that this was a near-complete shutdown of operations in the area, rather than one dictated by a short suspension in operations within a specific area (under the traffic light system in KSMMA following an event of ≥M3) or the cessation of seismicity following reservoir depletion.R1 Given that prior to the development of the Montney play this area was relatively characteristically ROS quiet in terms of natural seismicity (Lamontagne et al., 2008), the detection of latent seismicity over ∼4 months suggests lingering changes in the stress field to allow for its generation.

COVID-19 and the reduction of noise globally
The year 2020 was highly R1 unusual due to the global pandemic that caused the shutdown of many businesses and severely restricted the movement of people worldwide.This reduction in ground motion has been accurately measured by a drop in seismic ambient noise in many places, and correlated with a decrease in population mobility (e.g.Lecocq et al., 2020;Dias et al., 2020).The noise level at a seismic station can be estimated using the probabilistic power spectral density (PPSD) of its records (McNamara and Buland, 2004).Following the methodology of Lecocq et al. (2020), R3 we compute the PPSD from 30-minute windows with 50 percent overlap so that a single value is gained for each window, calculated using Welch's method (Welch, 1967) for the vertical Z-R1 component atof ROS different seismic stations.We use the vertical component since a number of public stations (including R25AC, Fig. 1), are single component seismometers, enabling us to compare the reduction in noise across a number of sensors.R3 This method reduces numerical noise in the power spectra at the expense of reducing the frequency resolution because of frequency binning, but this effect is minimized with a robust smoothing parametrization.The 30-minute time series are then converted to an average daily PSD, and the RMS of the time-domain displacement is extracted.Anthropogenic cultural noise typically concentrates at high frequencies (> 1-10 Hz, McNamara and Buland (2004)), but is strongly diurnal (e.g.stronger during the day than at night, and stronger during the weekdays compared to the weekends (Lecocq et al., 2020)).To capture this noise, but ROS avoid meteorological signals, and in particular to avoid focusing on ROS oceanic microseisms (which typically manifest below 1 Hz), we use the frequency band of 4-14 Hz to investigate seismic noise during the pandemic.
Figure 1 shows the reduction of seismic noise in the frequency band 4-14 Hz in Gastown, Vancouver, BC during the global pandemic.A clear reduction in noise is observed following the closure of schools (black line) and businesses (red line).During Phase I of the pandemic (i.e. between the closure of businesses and the partial reopening of the city on 5 May 2020 (green line; red background colour)) ROS , noise levels remain lower than previously recorded.Following the re-opening of some businesses in May and June 2020, the increase in noise is interpreted as the increased movement of people, although it remains lower than pre-pandemic levels.Following the reopening of some businesses in May and June 2020, an increase in the noise is seen, although it remains lower than pre-pandemic levels, interpreted as the increased movement of people.R3 To verify that these variations do not occur on an annual basis, we undertook the same noise analysis for the year 2019, and found no such fluctuations during the corresponding months.In fact, the average R2 ground displacement remained between 20 and 30 nm at station R25AC for the entirety of 2019.

Seismicity in the KSMMA
With increasing unconventional ROS oil and gas operations within the KSMMA over the past decade, the number of public monitoring stations has also increased.Prior to 2020, 9 public sensors maintained by Natural Resources Canada, the BC Oil and Gas Commission, the BC Seismic Research Consortium ROS and the Geological Survey of Canada existed within the KSMMA boundary, along with 6 co-located accelerometers poised to better capture higher levels of ground motion from larger seismic events.In early 2020, 13 additional broadband seismic stations (Trillium T120 seismometers with Taurus digitizers, sampling rate of 200 Hz R3 ) and two Titan accelerometers (sampling rate of 200 Hz) R3 were installed within the KSMMA (expanding the EON-ROSE (EO) network) as part of a joint project between the University of Calgary, Nanometrics, Geoscience BC and a number of universities in South Korea, ROS to monitor ongoing seismicity associated with hydraulic fracturing operations (Fig. 2).Sensors were installed at existing well sites through the support of four companies, and placed just below the surface (to a depth of 30 cm) to enhance coupling and decrease noise.The primary aim of the installation in 2020 was to expand monitoring capabilities within the KSMMA, in particular in the NE and SW of the area where prior public monitoring was sparse.The sites of the accelerometers were chosen due to their proximity to the most recent seismicity in the area, in particular a number of felt events that have occurred close to Tower Lake and Farmington.Further details of the installation of this array can be found in Salvage et al. (2021).(Salvage et al., 2021).R3   The catalogue of seismic events detected in the KSMMA is based on the newly installed array and available public stations in the area.Events were detected from the incoming continuous seismic data using an STA/LTA triggering algorithm, followed by a separate R1 template-matching algorithm.The template-matching algorithm uses previously detected events from its own catalogue, as well as historical seismicity in the area to identify seismic events from seismic noise, and consequently removing unwanted signals.A machine learning technique was then used to identify phase arrivals within suspected events, using historical seismicity (from both public and private arrays) as a training database.utilising continuously re-trained modules that classify noise from events and remove unwanted signals.Then, a support vector machine (SVM) machine learning technique was used to identify phase arrivals.These phase arrivals are identified by training an SVM model on historical data.R1     By converting the waveforms into over 250 "features" (e.g.frequency content, P-S timings), machine learning enabled the association of such features with P and S phases (or conversely, with noise) using quantities such as time and band-normalized spectrograms, a model is generated which can associate the features with P and S phases (or conversely, with noise) R1 , allowing accurate automatic phase picking of detected events (Salvage et al., 2021) R1 .We take the catalogue of event times and P and S phases to cut waveforms from the continuous seismic data.Data are band-passed filtered between 1 Hz and 80 Hz (zero-phase), converted to displacement by removing the instrument response, and de-trended using both the mean of the waveform and fitting a linear function to the waveform with a least squares.We then and R3 determine hypocenter locations using NonLinLoc (Lomax et al., 2009(Lomax et al., , 2000)), a probabilistic, global-search non-linear algorithm that generates the maximum likelihood hypocenter location based on the estimated posterior probability density function for each event.A 1D velocity model, specifically calibrated for the KSMMA from compressional and shear sonic logs, formation tops and ground truth locations of previous seismicity (available directly from the ROS BC Oil and Gas CommissionOGC ROS ), was used for location analysis R2 .Events were then re-located using HypoDD, a double difference algorithm, whereby the residual between the observed and calculated travel-time difference (or double-difference) between two earthquakes observed on a single station are related to differences in their relative hypocenter locations and origin times (Waldhauser and Ellsworth, 2000).To calculate magnitudes we use a form of the Richter (1935) magnitude formula that has been modified to better reflect local attenuation characteristics within the KSMMA (Babaie-Mahani and Kao, 2020).In line with calculations conducted by Natural Resources Canada (NRCan), we calculated M L using the maximum amplitude from the vertical component simulated on a Wood-Anderson (WA) seismometer, rather than the horizontal component, which has been used elsewhere.
Historically, seismicity within the KSMMA appears to occur within spatially distinct regions that fall within a corridor orientated NW-SE(Fig.2) R3 .Figure 2 shows the detected seismicity within the KSMMA in 2018 (Fig. 2(a), Visser et al. (2020)), compared to 2020 (Fig. 2(b), this study).We are unable to compare the seismicity from 2019 since these data are not yet published.R3 In both years, the largest magnitude event occurred in an area away from the densest occurrence of seismicity.
Since the largest event in 2020 did not occur in the same cluster as the largest event of 2018, it appears that the occurrence of M L 3+3-4+ R3 events, i.e. those that may be felt, R3 is not necessarily confined to a single region.Temporally, seismicity within the KSMMA occurs in distinct clusters, attributed to ongoing development activity in the area (Fig. 3).In 2018, heightened periods of seismicity were observed in April, May, July and August (Fig. 3(a)).Similar periods of heightened seismicity were observed in 2020 in March, August and September (Fig. 3(b)).The majority of seismicity detected within the KSMMA is M L ≤ 2, and consequently goes unfelt.

Prior and Post Lockdown: 2020
In March 2020, the Province of British Columbia introduced measures aimed at slowing the spread of COVID-19, including the closure of schools and childcare facilities on 17 March, and the closure of many businesses (in particular those that included daily human interaction) on 21 March.Although these closures did not include a government enforced closure of hydraulic fracturing activities, the stalling economy led to stock prices of operators within the KSMMA plummeting to record lows, and as such the suspension of most activities (S.Venables, BC Oil and Gas Commission, Pers.Comm., January 2021).ROS Up until this point in 2020, similar patterns of seismicity to other years were observed in the KSMMA (Fig. 3).A total of 4,268 events were detected from the onset of data collection (22 January 2020) from the updated EO array (yellow triangles, Fig. 2) to 1 April.Following the initial closure of businesses on 21 March, there is evidence of ongoing hydraulic fracture operations for ∼10 days, with associated heightened seismicity (Fig. 3(b)).This reflects the continuation of planned operations by companies within the KSMMA, following which, no new operations were initiated due to the diminishing economy.It is possible this reflects operators in the area undertaking additional hydraulic fracturing jobs during this time, as government restrictions became increasingly tight, and no "end-date" to the restrictions being suggested, or it may be that these 10 days of seismicity represent the continuation of planned operations with no initiation of new jobs.R1 Magnitudes of recorded seismicity prior to lockdown at the end of March range from M L -0.73 to M L 2.93.
At the beginning of April, a period of relative seismic quiescence began in the KSMMA (Fig. 3(b); grey background ROS ).
The COVID-19 pandemic not only limited the movement of people and shut businesses, but also caused an economic downturn in the energy sector, leading to the cancellation of many operations.Operations were once again restarted in British Columbia in the later summer months, after ∼4 months.ROS Seismicity since the resumption of activities in the later summer months, ∼4 months after the lockdown began, ROS is once again temporally clustered, with a total of 3,1762,617 R2 events being recorded fromsince ROS 6 August 2020 to 1 January 2021to present R2 .The largest magnitude event of 2020 at the time of writing R2 occurred on 11 September at 22:37 UTC with an estimated M L of 3.1, after which proximal operations were suspendedshut down ROS in line with the traffic light protocol introduced for the KSMMA (BC Oil and Gas Commission, 2018).A total of R1 73 precursory events occurred over approximately 4 hours prior to the M L 3.1 event ROS , with events locating within a small spatial extent (∼300 m x 150 m), probably directly related to ongoing operations in the area due to the correlation in space and time of events and fluid ROS injection.Events within this precursory sequence had magnitudes between M L 0.2 and M L 2.6, and were all located at depths of approximately 2.05 km.Moment tensor results for the M L 3.1this ROS event suggest a focal mechanism dominated by strike-slip (Salvage et al., 2021).

Evidence of reduction in seismic noise
A clear reduction in the number of seismic events was observed during the lockdown period from April to August 2020 in the KSMMA (Fig. 3(b); grey background ROS ).Over the ∼ 4 months of relative quiescence only 389 events were detected using the EO network and available public stations in the area.For comparison, 344 events were detected on our network over a single week from 8 to 15 February when operations were fully underway.On average during this period, the magnitude of events were smaller than during time periods when activity was driven by ongoing operations.
A reduction in seismic noise and therefore ground motion is also evident in the KSMMA following the introduction of government restrictions and the subsequent economic downturn ROS in March 2020 (Fig. 4).Unfortunately, the most central seismic stations in the EO array were not installed until immediately before and immediately after the end of lockdown (March and May, respectively),(March) or post lockdown (May) R3 and therefore could not be used to analyze the long term changes in seismic noise.We chose station KSM08, located in the east of KSMMA, approximately 14 km due north of the settlement of Rolla, R3 due to the long, uninterrupted seismic data recorded at this station, as well as its proximity to recent dense clusters of ongoing seismicity (Fig. 2).The recent seismicity close to KSM08 suggests that a number of wells in the vicinity were active prior to the economic downturn related to the COVID-19 pandemic and government imposed lockdowns in March 2020.R3   Heightened seismic ground motion is evident at KSM08 through January to March, as operations are ongoing (Fig. 4).A significant decrease in seismic ground motion is observed following the government restrictions ROS in late March 2020, with the average displacement sitting well below the daytimeweekday and weekend daytime ROS mean calculated prior to lockdown.
As restrictions ease, we see a large increase in ground motion following the reopening of businesses in May 2020, although this once again tails off through June and July.The decrease in ground motion in the late summer is observed at all of the stations in the network, to varying degrees, and is believed to be related to the downturn in the industry (as a result of the COVID-19 pandemic) during which operations were suspended by companies in response to the COVID-19 pandemic (S.Venables, BC Oil and Gas Commission, Pers.Comm., January 2021).R3 The re-introduction of operations in August is clear from an increase in ground displacement and seismic noise, which has remained elevated ever since R3 (although not as high as R3 pre-lockdown levels)since R3 .

Latent Seismicity during relative quiescence: 2020
Seismicity occurring during the period of quiescence from April to August 2020 within the KSMMA exhibit a number of characteristics indicative that it is a (latent) consequence of previous operations in the area.Figure 5 shows the temporal and spatial evolution of seismicity during this period.Firstly, perhaps unsurprisingly, seismicity does not occur in a distinct temporal pattern that exhibits clustering (Fig. 5(a)).A small number of events (∼5) occur each day throughout the 4 month periods ROS .
Event magnitudes also reveal no discerning patterns with time, with all events registering M L -0.66 to M L 1.2.Furthermore, the frequency index (FI) suggests no temporal patterns during the period of relative quiescence.The FI is a proxy for the spectral content of each waveform based upon the ratio of energy in low and high frequency windows (Buurman and West, 2010), calculated at a single station.We show our analysis foruse R3 station KSM06 (Fig. 2) due to its proximity to the majority of the ongoing seismicity during this period of relative quiescence, although all stations within the KSM network suggested no temporal patterns in the FI R3 .A negative FI means the waveform is dominated by low frequency energy (in this case 1 -40 Hz); a positive FI demonstrates a majority of energy in the high frequency band (40.1 -80 Hz).In many environments (e.g.volcanic) a lower frequency content of the waveform is proposed as evidence for the direct role of fluids in the generation of the seismicity (e.g.Lahr et al., 1994;Chouet, 1996;Salvage et al., 2018).In hydraulic fracturing environments, events dominated by low frequency energy have been postulated to relate to aseismic slip and deformation (Zoback et al., 2012).R1   Overall, tT ROS he seismicity detected during the period of relative quiescence within the KSMMA shows no discerning temporal characteristics.
Spatially, seismicity detected during the COVID lockdown period exhibits characteristics that are similar to the previously detected seismicity in the KSMMA (Fig. 5(b)).Most events occur in a corridor orientated NW-SE, similar to the spatial distribution of seismicity prior to lockdown.Some spatial clustering is evident (e.g. in May in the south (shown in ROS yellow)), but given the limited number of events this is difficult to determine with certainty.Most events during the quiescence period occur at a R1 focal depths of ∼0-4 km, which is similar for events prior to lockdown within the KSMMA, if potentially slightly shallower.Target formations for hydraulic fracturing within the KSMMA (Upper and Lower Montney) typically sit between 2 km and 2.5 km2000 m and 2500 m R3 (total vertical depth), with salt-water disposal (SWD) injecting at shallower depths (M.Gaucher, BC Oil and Gas Commission, ROS Pers.Comm, November ROS 2020).This suggests that events detected during the quiescence were generated in formations similar to those that occur when active hydraulic fracturing and SWD wasis ROS ongoing.

Characteristics of Observed Seismicity
Seismicity generated during this period of quiescence appears to share many characteristics with seismicity generated during hydraulic fracturing operations within the KSMMA.Although low in number, the event rate per day remains fairly constant throughout the ∼4 month period of limitedno ROS hydraulic fracturing operations, with no apparent temporal decay (Fig. 5(a)).This contrasts the "usual" pattern of seismicity during active hydraulic fracturing operations, which is highly temporally (and spatially) clustered around the wells operating (Fig. 3) (e.g.Skoumal et al., 2015).Figure 3(b) also suggests no change in the recorded magnitudes of events pre-and during lockdown.Furthermore, no discernible changes in the FI are observed during the lockdown period (Fig. 5(a)).In volcanic environments (for which the FI was developed), a low FI is thought to be indicative of the presence of fluid moving through the system (e.g.Salvage et al., 2018).However, in volcanic environments, low FI values typically contain energy of ≤5 Hz, with high FI values focused up to 15 Hz.The waveforms observed here contain much higher frequency energy than this, which is a common trait of seismicity generated in hydraulic fracturing environments (e.g.Eaton, 2018).Zoback et al. (2012) have suggested that events dominated by low frequency energy in hydraulic fracturing environments may be indicative of slow-slip or aseismic deformation.Therefore, the low frequency nature of the events identified here may be indicative of either aseismic deformation or the presence of fluids within the system, however discerning between the two processes and their relative importance is still an area of active research.R2   The estimated ROS magnitude of completeness (Mc) during the lockdown period is ∼0.4 (Fig. 6) R3 .We calculate Mc using the maximum-curvature method of Wiemer and Wyss (2000).R2 The Mc for the entire catalogue to from 22 January to 1 October 2020 (n=7216) is estimated to be 0.074 (Salvage et al., 2021) ROS , suggesting that even though relatively few events were detected during this quiescence, the detection of small magnitude events is good.R1 Given the reduction in noise during the period of quiescence (Fig. 4), it is perhaps no surprise that a low Mc is identified.this is perhaps no surprise.ROS The estimated b-value (Gutenberg and Richter, 1944) of 1.96 (Fig. 6, calculated using a least-squares linear regression for events with magnitudes ≥Mc) R2 is similar to b-values estimated from seismicity associated with hydraulic fracturing experiments in Western Canada, suggesting an abundance of lower magnitude events (Igonin et al., 2018;Eaton et al., 2014).The fact that no large magnitude events were detected during the period of quiescence (no M L >1.5) is directly influencing the estimated b-value in this case.Interestingly, higher b-values (as in this case) ROS have typically been attributed to seismicity generated in normal faulting regimes (Schorlemmer et al., 2005;Amini and Eberhardt, 2019).The KSMMA is strongly influenced by the Fort St. John Graben complex, an asymmetrical half graben that has also undergone significant strike-slip and rotational movement upon reactivation of the basement faults in the area (Barclay et al., 1990), which may also be directly influencing the estimated b-value.Furthermore, in hydraulic fracturing environments, b-values of >2 have been associated with the stimulation of natural fractures at depth, with smaller b-values associating with large-scale tectonic faults (Wessels et al., 2011;Eaton and Maghsoudi, 2015).In our case, this would suggest that the seismicity being generated may be related to both complex natural fracture systems and large scale faults in the area, since the b-value lies close to this scenario change.isdirectly related to the complex natural fracture system, rather than any large scale faults in the area.ROS   Seismicity during the quiescence appears to be spatially concurrent with previous seismicity in the area (Figs.5(b) and 2).However, there appears to be very little correlation between the spatial extent of seismicity and the most recent hydraulic fracturing activity in the area (active in March 2020 prior to lockdown; green sqaures, Fig. 5(b) ROS ).Seismicity during the quiescence ROS appears in two planar elongated features, extending in a NW-SE direction, with lengths of up to 30 km (eastern segment), if assumed to be one feature.These features are not coincident with any known faults in the area (e.g.Furlong et al., 2020).Seismicity recorded during this period of quiescence is generally located at a similar depth to the target formations of the Montney (∼2 km), as well as in the formations above.This suggests hydraulically connected pathways above the injection zone, perhaps within mechanically stronger lithologies, as has been previously suggested by Eyre et al. (2019b) in the Fox Creek region of Alberta (another area undergoing intensive hydraulic fracturing operations).
The generation of induced seismicity has often been successfully correlated to a number of injection parameters, including the injected volume of fluid (e.g.Yu et al., 2019;Ellsworth, 2013) and/or the pumping rate (e.g.Goebel et al., 2017).Temporally data areis R1 too sparse to draw conclusions as to whether any of these parameters directly influence the generation of latent induced ROS seismicity within the KSMMA, although given that hydraulic fracturing operations during our period of interest were extremely limitedceased ROS , we know that almost all of ROS this seismicity cannot be a direct response of this type of fluid injection.Only one well pad was operated during this period of quiescence, from 19 to 24 July (S.Venables, BC Oil and Gas Commission, Pers.Comm., January 2021), and only 4 events have both temporal correlation (i.e.occurred during this time window) and spatial correlation (i.e.occurred within 5 km of this well pad).ROS However, there is evidence that a small number of seismic events identified from April to August 2020 may be associated with salt-water disposal ( ROS SWD) ROS .Within the KSMMA, only 8 SWD wells were active in 2020, compared to hundreds of hydraulic fracturing wells.Of these, only one well was active during our period of investigation (Fig 5(b)).We infer thatbelieve ROS the seismicity occurring on 13 April 2020 (Fig. 5a, upper panel), where over 20 events were registered on the same day (significantly above the background rate of seismicity during this quiescence), may be due to SWD, accounting for ∼6% of detected seismicity during this period of quiescence ROS .
In this case, ongoing sustained SWD occurred ∼2 km away from these ROS events occurring on this day ROS .This offset is not unusual for SWD and associated seismicity; Schultz et al. (2014) found an offset of ∼3.5 km between SWD and associated seismicity in Alberta.

Estimation of Noise
PPSD is one of the most common methods used to characterize ambient seismic noise.However, the level of smoothing, the size of the data window used in analysis and the methodology itself may all influence the PPSD calculation and distort features of interest (Anthony et al., 2020).Smoothing is primarily undertaken in order to reduce the uncertainty associated with the PPSD estimates, and means that short spikes in noise (e.g.due to wind gusts or seismic activity) do not dominate the spectrum.
In our case, the reduction in ground motion is much easier to determine from the average of the PPSD rather than individual estimates (Figs. 1 and 4, green vs. grey lines).Although we use a period smoothing of 0.025 octaves, this is likely to provide adequate spectral resolution of spectral peaks, as shown by Anthony et al. (2020) and therefore impacts our results minimally.
We also use a window of 30 minutes (overlapping by 50%) to try to reduce spectral leakage and variance when calculating the PPSD.
Earthquakes, and other transient signals, are likely to impact the estimation of ambient noise by generating large spikes in the data.However, the removal of seismicity from datasets is generally accepted as not necessary since they are low-probability occurrences within generally high-probability ambient seismic noise (McNamara and Buland, 2004).Only teleseismic earthquakes appear to have any real affect upon PPSD calculations (Anthony et al., 2020).A number of teleseismic events have been detected in the KSMMA during the period of quiescence analysis (e.g.M w 7.8 event on 20 July 2020, 99 km off the coast of Alaska), that may influence our calculation of PPSD.However, since we see no peak in the average ground motion at these times (e.g.no substantial peak in July 2020, Fig. 4), we suggest that teleseismic events are not majorly influencing our results.
One signal that does clearly influence our PPSD results in Figs. 1 and 4. 4 ROS is wind and other meteorological phenomena ROS .
Poor weather reported in the KSMMA, with wind gusts exceeding 80 km/hour at times (Government of Canada, 2020) ROS , R3     were observed at the end of Aprilbeginning of May ROS , R3 during an otherwise quiet time R3 period where there were limited hydraulic fracturing operations in KSMMA and limited movement of people due to lockdown measures.(i.e.no hydraulic fracturing operations in KSMMA, limited movement of people due to lockdown measures).R3 Since the noise generated from wind gusts penetrates a wide frequency band, we are unable to completely ROS filter it out.Using a filter between 4 and 14 Hz tries to eliminate some of these transient signals mostly associated with meteorological and oceanic conditions.

Generation of Latent Seismicity
The suspensioncessation ROS of almost all ROS operations within the KSMMA in the summer of 2020 allows us a unique insight into seismicity that cannot be directly correlated with injection, which has always beenis ROS the inferred triggering mechanism for most (if not all) of the seismicity within the KSMMA.The characteristics of the seismicity generated during this period suggest no fundamental differences in terms of temporal or spatial patterns or magnitudes to previous seismicity within the KSMMA that can be correlated with injection.In fact, many of the characteristics appear to be equivalent to events detected prior to lockdown.Prior to the development of the Montney play, low magnitude ROS natural seismicity within the KSMMA was undetectable, given the limitations of the available seismic networks in placealmost non-existent R1 .The Canadian National SeismographSeismic R2 Network (CNSN) recorded 20 earthquakes (M L 2.5 -M L 4.3) from 1984 to 2008 within the KSMMA ROS , which are assumed to be a mixture of natural events (Halchuk, 2009) and those associated with early anthropogenic activities in the area.Three clusters of events in November to December 1984 (M L 2.2 -M L 2.8), January to February 1992 (M L 2.5 -M L 3.5), and December 1992 to January 1993 (M L 2.5 -M L 4.1), have been associated with oil extraction and fluid injection (water) in the Eagle West and Eagle fields, just north of Fort St. John (Horner et al., 1994)mostly natural events (Halchuk, 2009)  ROS .One additional natural eventThe closest event to have ROS occurred with a significantly larger magnitude than this ROS occurred in March 1986 (M w 5.4) NE of Prince George, British Columbia (Halchuk, 2009), still a significant distance away from the study area ROS .
In order to investigate the likelihood that our detected seismicity during the period of quiescence ROS is natural seismicity, we calculate the expected recurrence rates of seismicity within the KSMMA greater than M L 2.5 from historical data, which is the magnitude of completeness used for the determination of seismic hazard maps in Canada due to detection thresholds from the Canadian public seismic network.The total number of earthquakes detected by the national network from 1984 to 2008 within KSMMA ROS was 20 (Halchuk, 2009), suggesting a recurrence interval of 0.83 events per year.It is nottherefore R2 unsurprising that during the period of quiescence, no events greater than M L 2.5 were detected, given this calculated recurrence interval ROS .
Following the Gutenberg-Richter formula (Gutenberg and Richter, 1944), it stands that there should be a 100-fold increase in the event rate to estimate the number of events >M L 0.5 (assuming a b-value of 1.0) R1 , suggesting an event count of 83.We use a b-value of 1.0, rather than our calculated b-value of 1.96 (Fig. 6) since this is the expected b-value for an area dominated by natural seismicity (e.g.Frohlich and Davis, 1993;Godano et al., 2014).R1 Therefore, a maximum of 21% of events detected during the period of ROS relative quiescence might be attributablecan be attributed R1 to natural seismicity, in addition to the 6% attributable to SWD and one fluid injection well ROS .However, this leaves ∼65% (a further 8% may be attributable to dynamic triggering effects, see below)Therefore, over 70% R2 of seismicity generated during this period of relative quiescence that is difficult to explaincannot be explained R1 by either of thesethis ROS mechanisms ROS , and we suggest is more R1 likely produced as a remnant to previous operations, and therefore directly related to previous states of stress.With events being generated over 4 months since the almost complete ROS cessation of operations, the state of stress at depth must be near-critical for an extended period of time in order to generate this "latent" seismicity.
It is widely reported that earthquakes can be generated by transient stress changes related to the passage of seismic waves (i.e."dynamic triggering") ROS (e.g.Wang et al., 2015;Van der Elst et al., 2013;Hill and Prejean, 2007)) ROS .In some cases, this dynamic triggering can also be delayed by days or weeks following a teleseism, potentially related to the re-distribution of pore fluid from the passing seismic waves (Brodsky and Prejean, 2005) or through initial aseismic slip on faults triggering seismicity (Shelly et al., 2011).During the period of quiescence (28 March to 6 August 2020), 43 earthquakes of >M6 were reported by the United States Geological Survey (2020), that may have the potential to cause dynamic triggering.We follow the methodology set out by Wang et al. (2015), whereby we first select only the teleseismic events that generated an estimated peak ground velocity of greater than 0.2 cm/s at any station within the KSMMA, as defined by Lay and Wallace (1995), whereby: (1) and: where A 20 is the peak waveform amplitude when filtered at 20s; M is the magnitude; δ is the epicenter-station distance (in degrees); and T is the surface wave period (assumed to be 20 s).This method identified 40 events from the original list of teleseismic events.We then calculated the β statistic (Matthews and Reasenberg, 1988) by: which is a quantitative measure of the level of dynamic triggering, representing the standard deviation in the background seismicity rate after a remote event.N 1 and N 2 are the number of earthquakes detected before (t 1 ) and after (t 2 ) the remote event, respectively.Here, we take t 1 and t 2 to be 12 hours.E(N 2 ) = N 1 xt 2 /t 1 is the expected number of earthquakes after the main shock based on the background seismicity rate.If no earthquakes occur in t 1 (i.e.before the main shock), N 1 is set to 0.25 based on the equivalent range of the probability density function (Matthews and Reasenberg, 1988;Hill and Prejean, 2007).When β ≥ 2, there is sufficient statistical evidence (at a 95% confidence level) that there is a significant increase in the seismic event rate following the remote event (Hill and Prejean, 2007).
We identify 7 remote earthquakes that generate a β value ≥2 (Fig. 8), including the largest magnitude ROS event to have occurred to date ROS in 2020 locatedthat occurred ROS 99 km SSE of Perryville, Alaska on 22 July at 06:12 UTC, with M w 7.8 (United States Geological Survey, 2020).T, although t ROS he increase in event count in the KSMMA following this remote event is difficult to determine unless the above statistical analysis is performed, as it is difficult to see an increase in seismicity when daily event counts are usedwithout statistical analysis ROS .In some cases however, such as following the M w 6.1 event on 31 May 2020, 43 km W of Lampa, Peru, a significant increase in the number of events detected in KSMMA is clear.
Our analysis therefore ROS suggests that a maximum of 8% of the seismicity detected during this period of relative quiescence may be attributed to dynamic triggering, in particular the events on 31 May 2020.This rudimentary calculation does not take into account the spatial migration of events in response to the teleseism; only that if a significant increase in event count following the teleseism in the next 12 hours is observed, all events within this 12 hour window are determined to be related to dynamic triggering.However ∼65%, however <70% ROS of the detected seismicity cannot be attributed to primary activation mechanisms such as this, and therefore in our opinion are the result of "latent" ongoing processes.
Understanding how the stress field has evolved in time and space within the KSMMA, in particular with respect to the development of unconventional oil and gas operations in the area, is an extensive ongoing research topic.Very limited data regarding the state of stress prior to unconventional resource development in the area is available, since natural seismicity was almost non-existent (e.g.Lamontagne et al., 2008) and details of such seismicity are often incomplete with large associated errors (e.g.Halchuk, 2009).Furthermore, the sparse nature of regional recording systems prior to unconventional development means that detailed understanding of the source characteristics of identified seismic events, which may allow insight into the state of stress, is not possible.Farahbod et al. (2015) investigated the changes in background regional seismicity within the Horn River Basin, British Columbia before and after hydraulic fracturing operations became prominent in the area.The Horn River Basin is similar to the KSMMA in the fact that it is dominated by unconventional resource operations.Their study suggests that background seismicity dramatically increases following the introduction of hydraulic fracturing, correlating in time and space with ongoing operations, in a similar manner to what we have observed in KSMMA.However, their study fails to detail how the seismicity responds when hydraulic fracturing operations stop completely and whether there is a long-term change in the overall background rate of seismicity.If found to be true, this would indicate that the state of stress has changed at depth to allow for such an increase.
Further evidence for a true change in stress at depth may be observed through changes in the principal stress axes.The World Stress Map suggests an average SHmax orientation of ∼40 • in the area close to Fort St. John, taken from borehole breakouts (Heidbach et al., 2018).Babaie-Mahani et al. ( 2020) calculated SHmax from focal mechanism analysis of 66 events with magnitudes between 1.5 and 4.6 within the KSMMA associated with hydraulic fracturing, and found the value to be ∼22 • to ∼33 • , significantly rotated relative to the estimate from the World Stress Map.This would suggest that there is complexity in the stress distribution within the KSMMA, particularly spatially, possibly related to the structural complexity in the area as a result of the Fort St. John Graben complex (Barclay et al., 1990).The dominance of this tectonic feature throughout the KSMMA influences the identified source mechanisms of seismicity, which are extremely varied within a small spatial extent, suggesting a complex and heterogeneous stress regime (Babaie-Mahani et al., 2020;Salvage et al., 2020;Amini and Eberhardt, 2019;Berger, 2012).As a consequence, determining a permanent change in the principal stress axes before hydraulic fracturing began, during the dominance of hydraulic fracturing operations, and after the investigated period of quiescence requires further investigation.
The analysis presented here provides a rare opportunity to study seismicity in a period of relative quiescence during a time when the cessation of hydraulic fracturing operations within the KSMMA is not temporally or spatially limited (for example, as is the case when a seismic event ≥M3 occurs) nor is it due to the depletion of the entire reservoir when operations cease indefinitely.Instead, operations were suspended for ∼4 months due to an economic downturn as a result of the COVID-19 pandemic, across a wide area, which has once again seen an uptake in activity since August 2020 when the market recovered (Fig. 3(b)).We are interpreting our identified seismicity as "latent" (i.e.long-lived with no obvious direct primary activation source) and suggest that it may be a response to previous fluid injection in this area, despite its long-lived nature.This may well represent a new "normal" background rate of seismicity in this area, since pore pressures at depth are presumably higher (in particular within less permeable formations) due to previous injected fluid, thus reducing the effective normal stress at depth, or may be the result of fault weakening with time due to previous operations.However, until operations within the KSMMA cease permanently we will be unable to accurately measure this new background rate of seismicity, or determine whether the rates of seismicity observed here continue for sustained time periods (i.e.years).R1 The generation of seismicity in response to hydraulic fracturing is typically attributed to either fluid migration models, poroelastic phenomenon, or potentially aseismic slip (e.g.Bao and Eaton, 2016;Langenbruch and Zoback, 2016;Shapiro and Dinske, 2009;Segall and Lu, 2015;Eaton, 2018;Goebel and Brodsky, 2018;Eyre et al., 2019a).In the fluid migration model, pore fluid pressures are significantly increased upon fluid injection reducing the effective normal stress within a fault zone, which is sufficient to trigger seismicity (e.g.Peña Castro et al., 2020;Bao and Eaton, 2016).Given the temporal and spatial correlation between seismicity and hydraulic fracturing operations within the KSMMA, this appears to be a likely cause of some ROS seismicity.Under this model, the seismicity rate is usually observed to be proportional to the pore pressure, and is assumed to track the injection rate (Langenbruch and Zoback, 2016).Consequently, a slow and steady decrease in the rate of seismicity over time would be expected to occur, as fluid pressure leaks into the surrounding formations (Eyre et al., 2020), before seismicity returns to the background (i.e.natural) rate.Since seismicity during the period of quiescence is long-lived, shows no decay and cannot be attributed to increased fluid injection, another process must be involved in its generation.
Furthermore, if pore fluid pressure and relaxation as a direct consequence of immediate injected fluid was the trigger of the seismicity during this period of quiescence, we would expect the seismicity to spatially migrate directly outwards from the most recently injected wells.We see no evidence of this (Fig. 5(b)), suggesting direct pore fluid migration cannot be held responsible for the triggering of this sequence.
Seismicity triggered by pore pressure diffusion can also ROS be estimated by determining the propagating pore pressure fluid front (r t ) related to the hydraulic diffusivity in a homogeneous isotropic saturated poroelastic medium (Shapiro and Dinske, 2009;Parotidis et al., 2003) by: where D is the hydraulic diffusivity and t is time since injection.If the triggering front (r t ) closely follows the maximum distance of seismicity through time, then pore pressure diffusion is thought to play a central role in the triggering of this seismicity (e.g.Shapiro and Dinske, 2009;Parotidis et al., 2003).Diffusivity (D) is generally assumed to range in the Earth's crust between 0.1 m 2 /s and 10 m 2 /s (Scholz, 2019), although in areas affected by hydraulic fracturing it R3 is thought to generally be in the range of 0.1 m 2 /s to 2 m 2 /s (e.g.Goebel et al., 2017;Shapiro and Dinske, 2009;Parotidis et al., 2003).Yu et al. (2019) suggested similar diffusivity values determined from seismicity related to hydraulic fracturing in the Montney formation to the NW of KSMMA, although others have speculated that much smaller diffusion values would be expected in shale formations (Eyre et al., 2020;Guglielmi et al., 2015b).Higher values of diffusivity in hydraulic fracturing scenarios are believed to be the result ofanticipated due to ROS faults and fractures at depth acting as fluid corridors (e.g.Riazi and Eaton, 2020;Caine et al., 1996) (Caine et al., 1996) ROS , compared to in-tact shales ROS .However, the seismicity generated in the KSMMA during the period of quiescence shows no coherence with a triggering front from the most recently active injection wells (Fig. 7), suggesting that pore pressure diffusion is not the dominant mechanism responsible for triggering these earthquakes.
Other models proposed for the generation of seismicity in response to hydraulic fracturing suggest that both pore pressure and poroelastic effects are jointly responsiblefeasible mechanisms ROS (e.g.Segall and Lu, 2015;Goebel and Brodsky, 2018).In these instances, the increased pore pressure due to injection is thought to load the surrounding rock matrix, altering the stress field, often at great distances from the original injection site, if the region is well hydraulically connected.Again, however, thesethis ROS models ROS suggests ROS that seismicity is generated as a response to injecting fluid into the Earth, which was occurring only on a very minor scalenot occurring ROS at the time of thisour latent ROS seismicity.Given that the stress field is unlikely to be sustained at criticalwould likely diminish R1 following the cessation of fluid injection, we would also expect a decay in seismicity with time (e.g.Utsu, 1961).We do not observe this.Alternatively, the trapping of fluids within a fault zone with only minor fluid migration along the fault, could result in slow changes to the effective stress due to changes in pore pressure (Sibson, 1992).In this method, seismicity should migrate spatially outwards from this fault zone as the effective stress migrates.We also see no evidence of this spatial migration (Fig. 5(b)).
Recently, it has beenEyre et al. (2019a) have ROS suggested that aseismic slip may play an important role in the generation of seismicity associated with hydraulic fracturing at distances extending beyond the fluid-pressurized zones through the transmission of an elastic stress perturbation (e.g.Eyre et al., 2019a;Bhattacharya and Viesca, 2019;Guglielmi et al., 2015a;Cappa et al., 2018;Wei et al., 2015).For hydraulic fracturing regimes in Canada, Eyre et al. (2019a) suggested that, whereby R1 distal unstable regions of a fault may beare ROS loaded by aseismic slip that initiated due to an increase in pore pressure within a stable zone, leading to the generation of seismicity.Once slip is initiated, far-field intraplate stresses may repeatedly reload unstable regions of the fault, leading to relatively steady seismicity rates, which may continue for long periods of time R1 .They suggest the driving stresses of such behaviour are most likely to be elevated pore pressures (as a result of ongoing hydraulic fracturing in the area) becoming trapped within fault zones due to low permeabilities within many formations, sustained by tectonic forces ROS .Cappa et al. (2018) suggested that far-field aseismic deformation may also be sustained by fault weakening, which may accelerate with time.Fault weakening may be further enhanced due to degradation of slip surfaces from chemical and hydrothermal fluid-rock interactions, as suggested in geothermal reservoirs by Vavryčuk and Hrubcová (2017).R1 Given that in the absence of the cessation of operations the detection of latent seismicity is extremely difficult, there are few examples of long-lived seismicity associated with hydraulic fracturing operations that may offer insight into the generation of such activity ROS .One recent example comes from a long-lived seismic swarm in Alberta, where seismicity was observed over 10 months after injection ceased, and was interpreted as being driven primarily by aseismic slip (Eyre et al., 2020).The authors suggest that the steady seismicity rate over a number of months (swarm-like behaviour of seismicity, rather than a typical mainshock-aftershock sequence) and a lack of hypocenter migration cannot be easily explained by a fluid-migration model (which is often the most favoured model in hydraulic fracturing environments), and is instead better explained by an aseismic slip model (Eyre et al., 2019a).R1 We favour this interpretation of aseismic slip playing an important role in the initiation of this ROS seismicity since ongoing hydraulic fracturing operations are not required to generate ongoing seismicity (as is the case in the fluid-migration model) ROS ; instead, the previous trapping of fluids within fault zones may be enough to sustain the generation of seismicity.The latent seismicity identified here is persistent over ∼4 months, lacks evidence of hypocenter migration and shows no evidence of a mainshock-aftershock type sequence as the magnitude distribution remains constant (i.e. is swarm-like in its behaviour).In addition we see no evidence of a propagating pore pressure fluid front and detected events are dominated by lower frequencies, all of which appear to be characteristic of the aseismic slip model (Eyre et al., 2019a(Eyre et al., , 2020;;Bhattacharya and Viesca, 2019;Guglielmi et al., 2015a).Geodetic methods (e.g.GPS, InSAR) may be able to measure such aseismic deformation (e.g.Shirzaei et al., 2013;Biggs and Wright, 2020;Gualandi et al., 2017), however at this time no data of this nature is available for the KSMMA.R1 .

Conclusions
Seismicity generated in the KSMMA has always R1 been attributed to oil and gas recovery since production began R1 in the area, primarily due to its temporal and spatial correlation to operations.However, during the COVID-19 pandemic in the summer of 2020, almost all hydraulic fracturing ROS operations in the KSMMA were suspendedhalted ROS .Despite this pause in industrial activity ROS , 389 seismic events were recorded by our seismic network.These events occurred within the spatial extent of previous events in the area (a corridor orientated NW-SE), and had similar magnitudes to previously recorded seismicity (∼M L -1 to M L 1.2).The low magnitude of completeness (Mc = 0.4) is indicative of the general quietening of the area, as operations and the movement of people were restricted with government regulations, exacerbated by an economic downturn in the market ROS .The b-value of detectedgenerated ROS seismicity (∼1.96) is similar to previous estimates within areas dominated by hydraulic fracturing.Unlike during active hydraulic fracturing operations, events showed no temporal clustering, but instead were generated in a persistent swarm-likefairly constant ROS manner over the ∼4 months of quiescence.No spatial correlation between the most recently active wells in the area and seismicity could be determined, however the fact that seismicity occurs at the depths of previous injection (i.e.within ROS the target formations) suggests that the area is likely tomust ROS be hydraulically linked.
Since there is no temporal or spatial evidence that these events are a direct consequence of the most recent hydraulic fracturing in the area (i.e. an aftershock sequence driven by pore pressure diffusion or poroelastic relaxation), and since the area is typically naturally relatively quiescentquiet seismically ROS (a maximum of 21% of the detected events may be attributable to natural seismicity rates ROS ), we suggestconclude ROS that most of these events may beare ROS an indirect response of the increased pore pressures at depth which is causing aseismic slip on already pressurized fault zones, or could be the result of fault weakening, ROS as a result of previous fluid injection in the area.A number of events may be the result of dynamic triggering, ROS from remote events with M w >6 (up to ∼8%), however this process cannot account for the majority of the seismicity observed (∼65%>70% ROS ).We suggest that the prior fluid injection in the area has altered the state of stress, and caused fluids to become trapped in fault and fracture zones at depth, in close proximity to the original injection points R1 .
This allows seismicity to be primarily generated by aseismic slip loading unstable regions of these pressurized zones at depth.
Once slip has initiated, far-field tectonic ROS stresses may repeatedly reload these unstable zones, leading to the relatively stable seismicity rate that is observed.The detected seismicity here during the period of quiescence may represent a new (heightened) background rate of seismicity in the KSMMA, however whether this is the case will only become apparent once operations cease permanently in the area.ROS thank C. Chamberlain, S. Crane and an anonymous reviewer for their excellent comments and suggestions, which helped to much improve this manuscript, as well as the Editors of this Special Edition of Solid Earth.ROS   Note: the different time (x-axis) scales.R3 (a) Seismic events reported by NRCan between 1 January 2018 and 31 December 2018 (Visser et al., 2020).Magnitudes were calculated using the ML formula of Hutton and Boore (1987).Well B (Fig. 5) and time zero taken as the last day of injection at this well prior to lockdown (27 March 2020).Representative diffusion curves associated with hydraulic fracturing and waste-water injection (Goebel et al., 2017;Shapiro and Dinske, 2009) are shown.The events during quiescence in the KSMMA cannot be successfully modelled using pore pressure diffusion suggesting it cannot be a primary mechanism for generating this seismicity.

Figure 1 .
Figure 1.Filtered (4-14Hz) ambient seismic noise displayed as displacement from station R25AC (verticalZ R1 component) located in Vancouver, British Columbia.30-minute average PPSD (dark grey), with rolling mean (window size = 92 hours) shown in green.The timing of different lockdown scenarios for British Columbia are shown as vertical dashed lines.Background colours represent different lockdown scenarios: Green represents before and after lockdown scenarios; Red is the first lockdown scenario with the closure of schools and restaurants; Yellow is the second lockdown scenario where some businesses re-opened.ROS A clear reduction in the ground motion is observed following initial lockdown conditions in March 2020.The large peak in noise in September is thought to be meteorological, rather than a sudden increase in anthropogenic activity.Figure courtesy of codes developed by Thomas Lecocq and Fred Massin.

Figure 2 .
Figure 2. Spatial locations of seismicity concentrations within the KSMMA.Higher density of seismic events is indicated by brighter colours; lower density by darker colours; and no seismicity by grey.Densities displayed in (a) and (b) are directly comparable i.e. the absolute seismic density is the same in both subplots.R3 The outline of the KSMMA boundary is shown in yellow; public seismic monitoring stations as blue triangles; the newly installed EO network as yellow triangles; and co-located accelerometers as hexagons.FSJ1 and FSJ2are also part of the EO network but were installed in 2018.FSJ1 was decommissioned on 26 August 2020 but is shown for completeness as it was used in seismic analysis prior to this.The largest measured magnitude event within the KSMMA boundary for each year ROS is marked.(a) Seismic events reported by NRCan between 1 January 2018 and 31 December 2018(Visser et al., 2020).Note: although the new dense array was not installed at this time, it is shown on the map for reference.The largest event in 2018, occurring on 30 November, north of Tower Lake is shown (ML 4.5).(b) Seismic events recorded byon ROS the newly installed EO network (and incorporating data from public stations) from 1 January 2020 to 31 December 202022 January 2020 to 30 September 2020 R3 .The largest magnitude event, occurring on 11 September 2020 is indicated (ML 3.4).Figure courtesy of Thomas H. A. Swincoe, University of Calgary, developed R3 using QGIS (https://qgis.org/en/site/about/index.html).

Figure 3 .
Figure 3. Temporal evolution of seismicity within the KSMMA (daily counts in black, cumulative counts in blue) for 2018 and 2020 R2 .Magnitudes calculated using the ML formula of Hutton and Boore (1987).R3 Distinct temporal patterns can be observed in both years, in particular in the spring and autumn R3 , associated with ongoing hydraulic fracturing operations in the area.
The magnitude of completeness (Mc) is reported as ML 1.0.Seismicity within the KSMMA in 2018 from catalogue ofVisser et al. (2020).R3 (b) Seismic events recorded by the newly installed EO network (and incorporating data from public stations) from 1 January 2020 to 31 December 2020.Magnitudes were calculated using the specific local magnitude formula for the KSMMA introduced by Babaie-Mahani and Kao (2020).The magnitude of completeness (Mc) is estimated as ML 0.4.Seismicity within the KSMMA in 2020 from catalogue derived by the newly installed network (yellow triangles, Fig.2).R3 The timing of different lockdown scenarios affecting the KSMMA are shown as vertical dashed lines.The time period from April to August 2020 represents the period of relative quiescence analyzed in this paper (grey background) and is shown in more detail in Fig.5.due to the COVID-19 lockdown.ROS The time lag for seismicity build up after the Phase III reopening reflects the time required for operations in the area to be restarted.

Figure 4 .
Figure 4. Filtered (4-14Hz) ambient seismic noise displayed as displacement from station KSM08 (verticalZ R1 component) located within the KSMMA.30-minute average PPSD (dark grey), with rolling mean (window size = 92 hours) shown in green.The timing of different lockdown scenarios for British Columbia are shown as vertical dashed lines.Background colours represent different lockdown scenarios: Green represents before and after lockdown scenarios; Red is the first lockdown scenario with the closure of schools and restaurants; Yellow is the second lockdown scenario where some businesses re-opened.ROS A clear reduction in the ground motion is observed following initial lockdown conditions in March 2020, and significant increases in ground motion as lockdown measures are rescinded.Figure courtesy of codes developed by Thomas Lecocq and Fred Massin.

Figure 5 .
Figure 5. Temporal and spatial evolution of 389 events detected in the KSMMA during the cessation of operations from April to August 2020.(a) Upper: Daily event count and cumulative event counts.Middle: ML determined using the formula of Babaie-Mahani and Kao (2020).Lower: Frequency Index (FI) detailing the ratio of high frequency energy to low frequency energy within each detected waveform at KSM06.(b) Spatial evolution of events coloured by time and scaled by magnitude.Active wells that initiated seismicity in the month prior to quiescence (March 2020) are shown as green squares labelled A (most recently active prior to lockdown i.e. late March 2020) to E (active in early March); one active SWD well is shown as the grey star.Figure generated using GMT v.6 (Wessel et al., 2019).

Figure 6 .
Figure 6.Frequency-Magnitude distribution of events (n=389) detected in the KSMMA during quiescence from 28 March to 6 August 2020.Event counts in magnitude bins of 0.1 are shown as black columns; the cumulative event value per bin is shown as a red dot.The magnitude of completeness (Mc) is 0.4, and the estimated b-value is 1.96(Gutenberg and Richter, 1944).Local magnitudes were calculated following Babaie-Mahani and Kao (2020).ROS

Figure 7 .
Figure 7. Time-distance plots of latent seismicity from 28 March to 6 August 2020.(a) Distance of events measured away from Well A (Fig. 5) and time zero taken as the last day of injection at this well prior to lockdown (27 March 2020).(b) Distance of events measured away from

Figure 8 .
Figure 8. Teleseismic events (dotted vertical lines) that statistically (95% confidence level) generated dynamic triggering of seismicity within the KSMMA during a period of quiescence from April to August 2020.Daily event counts are shown as black bars.Statistical analysis (β statistic followingMatthews and Reasenberg (1988)) performed on the 12 hours prior to the teleseismic event, and 12 hours after.ROS