key: cord-0814717-blob2nuw authors: Bansal, Brijesh K.; Singh, S. K.; Suresh, G.; Mittal, H. title: A source and ground motion study of earthquakes in and near Delhi (the National Capital Region), India date: 2021-11-30 journal: Nat Hazards (Dordr) DOI: 10.1007/s11069-021-05121-w sha: 274ba3719e3c300737f35d1bfcf805af8413cec0 doc_id: 814717 cord_uid: blob2nuw We study source parameters of 10 local earthquakes (2.7 ≤ [Formula: see text] ≤ 4.5) that have occurred in the National Capital Region (NCR) since 2001 and the ground motions produced by these events. Moment rate spectra of the earthquakes retrieved from the recordings at hard sites after applying corrections for geometrical spreading (1/R, R ≤ 100 km), anelastic attenuation (Q = 253f(0.8)) and cutoff frequency (f(m) = 35 Hz) are reasonably well fit by the Brune ω(2)-source model with stress drop ranging between 0.9 and 13 MPa. Neglecting the outlier low-stress drop value, the average stress drop is 6 MPa. We apply a modified standard spectral ratio technique to estimate site effect at 38 soft sites in the NCR as well as the geometrical mean site effect with respect to a reference hard site. Application of the stochastic method, with source characterized by the Brune ω(2)- model with stress drop of 6 MPa and the mean site effect for soft sites, yields peak horizontal ground acceleration and velocity curves that are in good agreement with the observed values. These results provide the parameters needed for the application of the stochastic method to predict ground motions at hard and soft sites in the NCR during postulated [Formula: see text] ≤ 5.5 earthquakes. Delhi and its surrounding area (the National Capital Region, NCR) lies on an extension of the Peninsular shield, transverse to the Indo-Gangetic plains. The region is dominated by several tectonic features: Himalayan Main Boundary Thrust (MBT), Main Central Thrust (MCT), Delhi-Haridwar ridge, Delhi-Lahore ridge, Aravalli-Delhi fold, Sohna fault, Mathura fault, Rajasthan Great boundary fault, and Moradabad fault (Verma et al. 1995; Shukla et al. 2007; Prakash and Shrivastava 2012) . Some of these features are shown in Fig. 1 . As a consequence, the NCR is exposed to seismic hazards from large and great earthquakes along the Himalayan arc, and small and moderate earthquakes at local and regional distances. Historical seismicity in the NCR is plagued with uncertainty (see, e.g., Singh et al. 2013 for a summary). Catalogs based on instrumental recording list three significant earthquakes in and near NCR (Bansal et al. 2009 ). They occurred on 10 Oct 1956 near Bulandshahar, 27 Aug 1960 in Delhi, and 15 Aug 1966 near Moradabad. The 1960 Delhi earthquake is the largest, instrumentally-recorded event in the NCR. Hence, it has been one of the critical scenario earthquakes in the estimation of seismic hazard of the NCR. Yet the epicenter, depth, and magnitude of the earthquake are uncertain. The depth reported in different catalogs varies between 58 and 109 km and the magnitude between 5.3 and 6.0. From an exhaustive analysis of the event, Singh et al. (2013) concluded that it was shallow (depth ≤ 30 km but most likely ≤ 15 km) and its magnitude was M w 4.8 (range M w 4.6-4.9). A more reliable epicenter, as compared to the instrumentally-determined one, is provided by the area of the strongest seismic intensity: 28. 47°N, 77.00°E (between Delhi Cantonment and Gurgaon) . Quite possibly, the maximum magnitude in the NCR in the last 100 years has not exceeded 5.0. Abundant low-level seismicity and the awareness of earthquake exposure of NCR has led to a steady increase in the seismic and accelerographic instrumentation in and near NCR during last two decades (Bhattacharya and Dattatrayam 2000; Kumar and Mittal Shukla et al. 2018; Bansal et al. 2021a) . Consequently, our knowledge of the seismicity of the region has improved, revealing a prolific activity of small earthquakes. The seismicity, at least at a small magnitude level, is diffused (Pandey et al. 2020) . For this reason, it is not possible to relate events with specific tectonic features; rather the seismicity suggests an active source volume. This may be real or due to an error in the locations since the network is still sparse. Bansal et al. (2021b) made an effort to constrain the causative fault of the recent seismicity (April-May 2020) using morphotectonic analysis, Coulomb's failure stress changes, and fault plane solutions and inferred a new hidden fault. The improved seismic network in and near NCR has produced high-quality recordings which permits source and ground motion studies of some of the NCR earthquakes. Source parameters and analysis of ground motion have been reported previously for the earth- (Bansal et al. 2009; Singh et al. 2010) (Table 1) . A preliminary analysis of the earthquake of 5 March 2012 ( M w 4.5) was presented by Bansal and Verma (2012) . Recordings of Delhi and some Himalayan-arc earthquakes have been used to estimate seismic-wave amplification at several sites in the NCR with respect to the reference hard site of Ridge Observatory (station NDI) (Singh et al. 2002 (Singh et al. , 2010 Bansal et al. 2009; Mittal et al. 2013a, b) . An earthquake sequence began in the NCR on 12 April 2020 with an M w 3.3 event. It was followed by earthquakes on 10 May 2020 ( M w 3.3), 29 May 2020 ( M w 4.0), and 3 July ( M w 3.9). These earthquakes were strongly felt in the NCR region and caused panic among the public. Regional moment tensor inversion and source spectrum of the 12 April earthquake as well as an analysis of the ground motion is presented by Pandey et al. (2020) . The 12 April and 10 May events were nearly collocated. With respect to the earthquake of 12 April, the earthquake of 29 May occurred about 50 km to the west while that of 3 July was located about 105 km to the southwest (Fig. 1) . A seismotectonic framework interpreting the NCR events is given by Bansal et al. (2021c) . Recently, Tiwari et al. (2021) proposed that infrequent low magnitude earthquakes occurring in and around Delhi show a significant variation at annual seasonal scale and are correlated with the groundwater pumping for the irrigation, urban activities, and seasonally controlled hydrological loading cycle of freshwater aquifers. The main objective of this study is to develop parameters needed to estimate ground motion in the NCR from postulated earthquakes applying the stochastic method (Boore 1983 (Boore , 2003 . This effort is critical for reliable estimation of seismic hazard in the region. Towards this objective, we analyze recordings of seven NCR earthquakes since 2011, including four that occurred during the 2020 sequence (Table 1) . For each event, we estimate source spectrum, stress drop, and site amplification with respect to the reference hard site of NDI. The results from these seven earthquakes, along with those previously reported for the three other events, are used to develop the needed parameters for the prediction of ground motion. Amplification is estimated at 38 soft sites in the NCR using the recordings of the 10 earthquakes applying the standard spectral ratio (SSR) technique, accounting appropriately for the geometrical spreading. From the estimated site effects, we present median soft-site amplification of seismic waves. Thus, the ground motion from a postulated NCR earthquake may be estimated at any hard site and each of the 38 soft sites or a generic soft site. Alternately, real-time recording of an earthquake (located in/near NCR or in the Himalayan arc) at the reference site of NDI, along with the catalog of the site effects, may be used to estimate the ground motions at all other sites in near-real-time, valuable information for civil protection authorities. A brief description of the earthquakes is given first, followed by the common methodology used in processing and analysis of the data of all events. For consistency, some of the source parameters of the three previous events are re-analyzed using this common methodology. Earthquakes of 2001 Earthquakes of , 2004 Earthquakes of , and 2007 are discussed in Bansal et al. (2009) and Singh et al. (2010) . Here we give brief description of the more recent events analyzed in this study. The epicenter of the earthquake reported by the National Centre for Seismology (NCS) is in the heart of Delhi, about 8 km to the SW of Ridge Observatory (Table 1, Fig. 1 ). The NEIC epicenter (Table 1) is about 10 km N23°E of the NCR epicenter. This is one of the more significant events to occur in and near Delhi in recent years. It was located about 60 km to the west of Ridge Observatory, near the town of Rohtak. The epicenter of the earthquake falls near Mahendragarh-Dehradun fault (MDF) (Fig. 1) . The earthquake was felt in Rohtak, Gurgaon, Delhi, Faridabad, and as far as Jaipur. A preliminary report on the event was published by Bansal and Verma (2012) . Recordings of the event were used by Mittal et al. (2013b) to study site effect. As it is one of the larger, well-recorded earthquakes in/near Delhi, it merits a detailed analysis. A swarm-like earthquake activity occurred in Delhi on 11 Nov 2013. The sequence was widely felt in the city. Here we analyze the largest event of the sequence. It was located about 5 km NE of the offices of NCS located in the India Meteorological Department campus at Lodhi Road. The earthquake of 12 April 2020 ( M w 3.3), with its epicenter located 12 km N30°E of Ridge Observatory, was felt as far south as Faridabad, a distance of about 35 km. It caused panic in the city although no damage was reported. According to U.S. Geological Survey Community Internet Intensity Map, which was based on 245 responses, the felt area extended ~35 km south of the NCS epicenter (Table 1 ). This small magnitude earthquake was felt by a relatively large number of persons considering perhaps because it occurred in the locked-down period of the Covid-19 pandemic when most people were indoors (Pandey et al. 2020; Singh 2020) . The earthquake of 10 May occurred near the 12 April hypocenter and had almost the same magnitude. Yet only four responses are reported in the U.S. Geological Survey Community Internet Map. As we show later, this earthquake was deficient in high-frequency radiation as compared to the 12 April earthquake. The earthquake was located about 50 km to the west of Ridge Observatory, 20 km N55°E of the 5 March 2012 event. Curiously, the epicenter reported by NEIC, U.S. Geological Survey is 50 km N35°E of the NCS epicenter (Table 1 ). The depth reported by NEIC is 45 km as compared to 15 km by the NCS (Table 1) . Because of the large difference between the two locations, we carefully relocated the event. The RMS residual as a function of depth has a sharp minimum (0.47 s) at a depth of 15 km. (S-P) time at the nearest station is 3.6 s. There are 4 stations with (S-P) ≤ 5.2 s. The NCS location is, clearly, more accurate than that reported by the NEIC. The epicenter of the earthquake was near Alwar, 93 km SW of Ridge Observatory (Fig. 1) . As the (S-P) time at the closest station KUDL is 2.3 s, it follows that the depth of the event was less than 20 km. This is consistent with H = 16.5 km reported by NCS, not with 32.8 km given by NEIC, U. S. Geological Survey (Table 1) . Our analysis extensively uses the Fourier acceleration amplitude spectrum of the S-wave group. The time window chosen for the spectral analysis brackets 95% of the energy on the horizontal components of accelerograms beginning with the arrival of the S-wave. (In some cases, acceleration was computed by differentiating broadband velocity seismograms.) A 5% cosine taper is applied before computing the spectrum which is then smoothed by a 1/6 octave filter. Henceforth, we denote the geometric mean of the amplitude spectra of the two horizontal components as A(f, R). Under far-field, point-source approximation, A(f, R) of an event at a distance R may be written as (see, e.g., Singh et al. 1999 for details): where In equations above, Ṁ 0 (f ) is the moment rate spectrum so that Ṁ 0 (f ) → M 0 (the seismic moment) as f → 0, R = hypocentral distance, R θϕ = average radiation pattern (0.55), F = free surface amplification (2.0), P takes into account the partitioning of energy in the two horizontal components ( 1∕ √ 2 ), β = shear-wave velocity at the source, ρ = density in the focal region, and Q(f) = quality factor, which includes both anelastic absorption and scattering. The attenuation in the near-surface layer and the finite bandwidth of the observed spectrum imposed by the sampling rate are accounted by the parameter κ (Singh et al. 1982; Anderson and Hough 1984) and/or the Butterworth filter, B(f). Often either κ or B(f) is sufficient to explain the high-frequency falloff of the spectrum. Following Boore (1983) , we take B(f ) = 1.0 + (f ∕fm) (8) (−0.5) . For the NCR earthquakes analyzed here, B(f) with f m = 35 Hz suffices. Site(f) in Eq. (1) is the local site effect which we will take as 1 for hard sites in the NCR. The geometrical spreading term, G(R), in Eq. (1) is taken as 1/R for R ≤ 100 km and 1/(100R) 0.5 for R > 100 km. Similar to our previous studies of the NCR events, we take Q(f) = 253f 0.80 which was estimated for the Himalayan arc region (Singh et al. 2004) . Q(f) has been estimated for the NCR region using coda waves (Mohanty et al. 2009; Das et al. 2018 ) and S waves (Sharma et al. 2015; Banerjee and Kumar 2018) . As discussed later, reasonable choices of Q(f) have a minor effect on the results. For the deeper earthquake of 2007 (H = 30 km; Table 1 ), β and ρ in the source region are taken as 4.2 km/s and 3.2 gm/cm 3 , respectively; for all other events the corresponding values are taken as 3.55 km/s and 2.85 gm/cm 3 . We take the logarithm of Eq. (1) and compute log Ṁ 0 (f )Site(f ) . Hardand soft-site recordings are processed separately. For each event, we show plots of Ṁ 0 (f ) from hard sites and Ṁ 0 (f )Site(f ) from soft sites. We interpret the hard-site spectrum by the ω 2 -source model of Brune (1970) , Bansal et al. (2009) and Singh et al. (2010) . They are listed in Table 1 . Note that some of the values given in the table differ slightly from those reported in the original publications. This is because of uniform values of β and ρ used here. (1) M 0 (f ) and [ Ṁ 0 (f ) Site(f)] of the earthquakes are shown in Fig. 2a -g. Both geometric mean and ± one standard deviation curves are presented. The panels identify hard and soft sites whose recordings were used in the computation. Most of the recordings were accelerograms. Seismic instrumentation in the NCR had undergone changes when the 2020 events occurred. For these events, the ground motion was recorded by both an accelerograph (100 samples/s) as well as a broadband seismograph (40 samples/s). Because the low-frequency signal is poorly resolved on accelerograms, for the 2020 events we patched Ṁ 0 (f ) and Ṁ 0 (f )Site(f ) constructed using the two datasets. A Brune ω 2 -source model was adjusted to Ṁ 0 (f ) by requiring that the theoretical curve simultaneously fit the observed spectrum both at low frequencies and high frequencies in the least-square sense. The selection of the low-and high-frequency band chosen for least-square fit was subjective. Theoretical curves are shown in blue in Fig. 2 . Each Ṁ 0 (f ) panel gives corner frequency, f c , and the stress drop, Δσ. The corresponding Superimposed on the plot is the theoretical curve for the ω 2 -source model (shown in blue). The disparity between the blue curve and the geometric mean curve (shown in bold red) is a measure mean site effect. From an examination of Fig. 2 we conclude that ω 2 -source model is a fairly good approximation to the observed moment rate spectra of earthquakes in and near the NCR. If we exclude the 10 May 2020 event, then Δσ ranges between 2.5 and 13 MPa, the average being 6 MPa. This is a fairly small range in view of the difficulty in reliable determination of f c and cubic dependence of Δσ on f c , Δσ α M 0 ( f c /β) 3 . Δσ of the 10 May 2020 earthquake (0.9 MPa) is much smaller than that of the nearby 12 April 2020 event (5.6 MPa). This agrees with the relatively small number of felt reports for the 10 May event. We envision that a patch of higher strength on the fault surface broke during the 12 April event, loading an adjacent patch of lower strength that, eventually, broke on 10 May. Next, we provide further evidence that the two events were collocated, had similar focal mechanism and M w but the 10 May event was deficient in high-frequency radiation. The displacement seismograms of the 12 April and 10 May 2020 events are compared at stations NRLA and NDI in Fig. 3 . Amplitudes and waveforms at NRLA and NDI differ because of focal mechanism effect and difference in the site response. However, the seismograms are nearly identical during the two events at the same site, confirming that they were nearly collocated, and had similar magnitude and focal mechanism. We checked whether the high-frequency radiation was deficient (Δσ was low) during the 10 May earthquake by computing spectral ratios of ground motion at each station during the two events. 12 April/10 May spectral ratios for the vertical component are shown in the bottom part of the figure. Around 1 Hz, the geometric mean ratio is ~1, implying the same seismic moment of both events. The ratio increases to 3 at 10 Hz and then decreases reaching 2 at 40 Hz. The figure demonstrates that the 10 May event was deficient at high frequencies as compared to the event of 12 April. We expect, and later Site effect in the NCR is clearly seen in the plots of Ṁ 0 (f ) and Ṁ 0 (f )Site(f ) presented in Fig. 2 . The amplification and variability of ground motion at soft sites can be appreciated in the plots of Fourier acceleration spectra during the 11 Nov 2013 earthquake (Fig. 4) . This event was recorded at two hard sites (NDI and DJB) and eight soft sites within a hypocentral distance range of 14.4-19.1 km. Since the distances are roughly equal, the spectra, in the absence of site effect, should be similar. However, as Fig. 4 illustrates, they are highly variable at soft sites and are consistently higher than those at the hard sites. Because of its importance in seismic hazard estimation, the site effect in Delhi has been studied through numerical modeling (Parvez et al. 2004 ), microtremor recordings (Mukhopadhyay et al. 2002) , borehole data (Iyengar and Ghosh 2004) , and earthquake recordings (Singh et al. 2002 (Singh et al. , 2010 Bansal et al. 2009; Mittal et al. 2013a, b) . The (right) . Traces suggest that the two events were nearly collocated, and had same magnitude and focal mechanism. (Bottom) (12 April/10 May) spectral ratios. Red continuous curve: geometrical mean of the ratios. Ratio shows that 10 May event was deficient at higher frequencies compared with 12 April earthquake estimation of site effect in Delhi using earthquake recordings is based on the application of the standard spectral ratio (SSR) technique. In this technique, the spectral amplification at the site of interest is measured with respect to a reference, hard site. For the NCR the station NDI (Ridge Observatory) is an ideal reference site as it is the oldest seismic station in Delhi that remains in operation. Most of the events recorded at other sites in Delhi are also recorded at NDI. (All events analyzed here were recorded at NDI). Since a large dataset is available at NDI, it is possible to predict ground motion from postulated earthquakes at this site. If ground motion at NDI is known (either recorded or predicted), then it can be computed at all other sites with known SSRs through the application of the stochastic method. SSR technique is based on the assumption that the source is sufficiently far as compared to the distance between the sites so that the expected ground motion, in the absence of site effect, is the same at the two sites. This assumption may not be valid for local events. Yet, site effect studies in Delhi using the SSR technique have ignored this issue. As the hypocentral distances to NDI and stations at soft sites differ, in this study, we reduce the S-wave spectrum of each of the two horizontal components at NDI to the distance of the station of interest by taking Q(f) = 253f 0.8 and geometrical spreading of 1/R (R is ≤ 100 km for all stations). The ratio of the geometric mean of the observed horizontal spectra at the station to the geometric mean of the reduced spectra provides the SSR. Note that in this modification of the SSR technique we are ignoring possible focal mechanism, path, and directivity effects. We analyzed 74 accelerograms produced by the 10 earthquakes listed in Table 1 . These accelerograms were recorded at 38 sites. In Fig. 5, we show examples of SSRs at six sites where three or more earthquakes were recorded. We note that SSRs at the same site are Fig. 4 Fourier acceleration spectra of 11 Nov 2013 earthquake at two hard sites (NDI and DJB) and eight soft sites located within hypocentral distance range of 14.4-19.1 km. The spectra at soft sites are highly variable and, as expected, are consistently higher than those at the hard sites reasonably similar at frequencies above about 0.8 Hz; below 0.8 Hz the SSR, in some cases, is affected by noise and is not reliable. Spectral amplification at a given soft site is critical in the estimation of site-specific ground motion. Here, however, we will be concerned with the prediction of ground motion at generic hard and soft sites in Delhi; for this purpose, we derive a mean site effect in the NCR. In Fig. 6 we show SSRs computed for the seven earthquake since 2011 (Table 1) . For each earthquake the geometric mean and ± one standard deviation SSR curves are given. The mean curves slightly differ, partly because the sites where the events were recorded were not the same. In the last panel of the figure, we have merged SSRs of these seven events with those from the earlier three earthquakes (Table 1 ). The geometric mean curve of SSRs of all 10 events is shown in the last panel. This curve is the mean spectral amplification at soft sites in NCR. It may be approximated by a constant SSR of 3.5 in the frequency range 0.25-10 Hz. As mentioned before, ground motions recorded during the 2001, 2004, and 2007 earthquakes were analyzed in previous studies (Bansal et al. 2009; Singh et al. 2010 ). Here Examples of spectral amplification at six sites with respect to the reference hard site of NDI (Ridge Observatory) estimated using a modified SSR technique. Three or more earthquakes were recorded at these sites. The red curve corresponds to the 3 July 2020 earthquake we include data from the earthquakes of 2011, 2012, 2013, and the four events of the 2020 sequence to those of the previous three events for a more comprehensive ground motion study. Plots of peak horizontal ground acceleration (PGA) and velocity (PGV) as a function of distance R of the ten events (2.6 ≤ M w ≤ 4.5) is shown in Fig. 7 . Hardand soft-site data are identified by solid and open circles, respectively. Sites in/near NCR which we have classified as hard are NDI, DJB, SONA, AYAN, and KUDL. For some events, we have added data from stations of regional broadband networks, located at R ≤ 300 km. As expected from the previous discussion, PGA and PGV values from the 12 April 2020 earthquake are significantly higher than those from the 10 May 2020 event even though the two earthquakes had the same magnitude. For each earthquake, geometric mean and ± one standard deviation curves are given. In the last panel, the spectral ratios of these seven events are merged with those from the three earlier earthquakes. The geometric mean of spectral amplification at soft sites in NCR and ± one standard deviation curves of the 10 events are shown in this panel. The mean amplification curve may be approximated by a constant value of 3.5 in the frequency range 0.25-10 Hz We refrain from developing a ground motion prediction equation (GMPE) from the available sparse dataset. Instead, we seek parameters that, in conjunction with the stochastic method (Boore 1983 (Boore , 2003 , may be used to estimate ground motion from future earthquakes. The method requires an estimation of the Fourier spectrum of the ground motion at the site of interest. Here we compute the spectrum using Eq. (1). We assume a Brune ω 2 -source which is completely specified by M 0 and Δσ of the event. We seek a generic Δσ, which along with Q(f) = 253f 0.80 and f m = 35 Hz, would generate PGA and PGV in agreement with those observed during the NCR earthquakes (Fig. 7) . This Δσ may then be used in computing ground motions during postulated earthquakes. A parameter needed in the application of the stochastic method is the duration of the intense part of the ground motion, T d , which we take as T d = 1∕f c + 0.05R + 3.0. We first check whether the Brune ω 2 -source model, with appropriate M 0 and Δσ for the event (Table 1) and Q(f) = 253f 0.80 and f m = 35 Hz, generates spectrum that agrees with the observed spectrum at the reference NDI site. We note that the accelerograms of the Fig. 7 Observed peak horizontal ground acceleration (PGA) and velocity (PGV) as a function of distance, R, of the 10 events analyzed in this study. Note that PGA and PGV values from the 12 April 2020 earthquake are significantly higher than those from the 10 May 2020 event even though the two earthquakes were collocated and had the same magnitude. For each event, the predicted PGA and PGV curves, based on the stochastic method assuming a Brune ω 2 -source model with Δσ = 6 MPa, are superimposed on the data. For earthquakes of 10 May and 3 July 2020 the predictions with Δσ of 0.9 and 11.2 MPa, respectively, also shown in the figure, better fit the data events before 2020 and those recorded during the 2020 events at NDI were recorded at 200 and 100 sps, respectively. Thus, to facilitate the comparison, we low-pass filtered the NDI accelerograms at 35 Hz. Fig. 8 Comparison of observed acceleration spectra (geometric mean of the two horizontal components) and predicted spectra at the reference hard site of NDI. Predictions are for a Brune ω 2 -source model using Eq. 1. Other parameters needed in Eq. 1 are discussed in the text The geometric mean of S-wave acceleration spectra of the two horizontal components are shown in Fig. 8 along with the predicted spectra. The two spectra are reasonably similar. The difference in the spectra is greater if the stress drop has been estimated from more hard-site recordings than from NDI alone. Spectra of the 2020 events differ from the computed spectra at high frequencies (f > 35 Hz) because of the low Nyquist frequency (50 Hz). As mentioned earlier, there are some studies which report Q(f) for the NCR region. To investigate the sensitivity of the results on the choice of Q(f), we compared our predicted spectra at NDI with those from two extreme Q(f) relations for the NCR: Q(f) = 162f 0.835 (Banerjee and Kumar 2018) and 98f 1.07 (Sharma et al. 2015) . Figure 9 shows theoretical and observed spectra for the earthquakes of 2007 and 2012. Here theoretical spectra have been computed with f m = 35 Hz and the same Δσ which we found appropriate when Q(f) was taken as 253f 0.8 . We note that the theoretical predictions are indistinguishable if M w is changed slightly (the values are given in the figure) . Even in the case of large R (e.g., earthquake of 2012) and low Q(f) = 98f 1.07 the required change in M w is only 0.1. We conclude that a reasonable choice of Q(f) has only a minor effect on the results. To predict ground motions from postulated earthquakes, we take Δσ of 6 MPa, the average value for the NCR events ignoring the anomalously low stress-drop earthquake of 10 May 2020 (Table 1) . PGA and PGV curves as a function of R computed via the stochastic method with Δσ = 6 MPa are shown in Fig. 7 along with the observed PGA and PGV values. Computations performed excluding and including mean site effect are shown in the figure. Considering the large scatter in the observed data, the predicted curves with 3) are exceptionally low; the predicted curves with Δσ of 0.9 MPa fit the data well. We surmise that the ground motions in the NCR at typical hard and soft sites may be predicted using the stochastic method assuming a Brune ω 2source model with Δσ = 6 MPa. Predicted PGA and PGV curves at generic hard and soft sites in the NCR as a function of R for M w 3.5, 4.5, and 5.5 earthquakes are illustrated in Fig. 10 . We refrain from predicting ground motions during M w ≥ 6 earthquakes for two reasons. First, there is no evidence of the occurrence of an M w > 5 event in the NCR in last 120 years. In fact, as mentioned in the introduction, there is strong evidence that the largest instrumentally recorded 1960 Delhi earthquake had a magnitude less than 5 (Singh et al. 2013) . Second, simulation of M w ≥ 6 earthquake at short distances should include both finite-fault and near-field effects. This is not possible in our stochastic simulation which is based on point source, far-field Fig. 10 a and b Predicted PGA and PGV curves at generic hard and soft sites in the NCR as a function of R for M w 3.5, 4.5, and 5.5 earthquakes. Predictions are based on the stochastic method assuming a Brune ω 2 -source with Δσ = 6 MPa. c Observed PGA during 5 March 2012 ( M w 4.5) earthquake along with predicted curves from the stochastic method (this study) and GMPE of Srinagesh et al. (2021) constructed from data of western Himalayan arc earthquakes. Predictions are shown for both hard (continuous) and soft sites (dashed). d Observed PGV during 5 March 2012 earthquake and predicted curves from the stochastic method and the GMPE of Srinagesh et al. (2021) approximation. There are some studies which simulate ground motions for M w ≥ 6 earthquakes in the NCR (e.g., Mittal et al. 2016; Sandhu et al. 2020) . These studies, however, do not account for the near-field effects. Several GMPEs have been developed for Himalayan arc earthquakes (e.g., Singh et al. 1996; Sharma 1998; Iyengar and Ghosh 2004; Sharma et al. 2009 ). GMPE of Singh et al. (1996) is based on recordings of only 5 Himalayan earthquakes. Iyengar and Ghosh (2004) derived a GMPE which they utilized in their seismic hazard study of the NCR. A subset of the data employed by Iyengar and Ghosh had previously been used by Sharma (1998) in the derivation of a GMPE for PGA. Sharma et al. (2009) derived a GMPE merging data from the Himalayan and Zagros regions. Singh et al. (2010) found that the GMPEs of Iyengar and Ghosh (2004) and Sharma et al. (2009) greatly overestimate observed PGA in Delhi during the 25 Nov 2007 ( M w 4.1) earthquake. This was not surprising since the data used in the construction of these GMPEs came almost entirely from M w ≥ 5.5 earthquakes. Singh et al. (2010) reported that these GMPEs also overestimate PGA in the Delhi region from earthquakes originating in the Himalayas. They concluded that these GMPEs should neither be used to predict motions in the NCR from local M w < 5 earthquakes nor from earthquakes originating in the Himalayan arc. A new GMPE has been derived by Srinagesh et al. (2021) which is based on the recordings of 14 shallow western Himalayan arc earthquakes (4.3 ≤ M w ≤ 7.9). As the data set includes seven earthquakes in M w range 4.3-5.0, we explored the validity of this GMPE in predicting ground motions in the NCR from local earthquakes. In Fig. 9 (bottom two panels) observed PGA and PGV during the Rohtak earthquake of 5 March 2012 ( M w 4.5) are compared with the predictions from our stochastic simulation and the GMPE of Srinagesh et al. (2021) . It is not surprising that the stochastic simulations fit the observed data better than the predictions from the GMPE of Srinagesh et al. which underestimates PGA at soft sites and overestimates it at hard sites. The same is true for PGV although the prediction at soft sites is only slightly underestimated. We surmise that attenuation relations developed for the Himalayan arc earthquakes may not appropriate for peak ground motions from earthquakes in the NCR. Stochastic simulations using parameters developed in this study are likely to provide more reliable ground motion estimations. We have studied 10 well-recorded local earthquakes (2.7 ≤ M w ≤ 4.5) in the NCR (including four that occurred in 2020) to retrieve their source parameters using hard-site recordings. All events have been analyzed using a uniform procedure. We find that spectra at hard sites are well-fit with a Brune ω 2 -source model, Q(f) = 253f 0.8 , f m = 35 Hz, and geometrical spreading of 1/R (Eq. 1). Other reasonable choices of Q(f) have only a minor effect on the results because of relatively short hypocentral distances involved (R < 100 km). Based on these results, estimated site effects at soft sites, and analysis of recorded ground motions, we have developed parameters needed for the application of stochastic method to estimate of ground motions at typical hard and soft sites in the NCR during postulated earthquakes. Below we summarize the results. 1. The source spectra of nine out of 10 events require a stress drop between 3 and 13 MPa (average ~6 MPa) ( Table 1 ). The exception is the earthquake of 10 May 2020 ( M w 3.3) which had a very low-stress drop, 0.9 MPa. This earthquake was located near the 12 April 2020 ( M w 3.3) event. Even though these two events had the same magnitude, the 10 May earthquake was significantly deficient in high-frequency radiation. As a consequence, observed PGA and PGV values were, relatively, low. Presumably, a patch of higher strength on the fault surface broke during the 12 April event, loading an adjacent patch of lower strength that, eventually, broke on 10 May. 2. Site effect has been estimated at 38 sites in the NCR using a modified Standard Spectral Ratio (SSR) technique, and taking NDI (Ridge Observatory) as the reference hard site. Only a few sites in the NCR have more than three recordings of earthquakes. The site effects at these sites during different earthquakes are fairly similar (Fig. 5) . We also computed the geometric mean spectral amplification at soft sites. Thus, if the Fourier spectrum of an earthquake at NDI is known either from a recording or via a source model with appropriate attenuation parameters (Eq. 1) then the Fourier spectra can be estimated at each of the 38 soft sites or a generic soft site in the NCR. The earthquake may have its origin in the Himalayan arc or within NCR. These spectra can then be used to estimate ground motion parameters via the stochastic method. This is a powerful technique to produce near-realtime ground motion map for the NCR using only the NDI recording, thus providing extremely useful information for civil protection authorities. Such an approach is currently in use in Mexico City (Ordaz et al. 2017 ). 3. Fourier spectrum at NDI from a postulated earthquake in NCR may be estimated using a Brune ω 2 -source model with a stress drop of 6 MPa, κ = 0.0 s, Q(f) = 253f 0.80 , and a cutoff frequency f m of 35 Hz in Eq. 1. Figure 10 shows predicted PGA and PGV at generic hard and soft sites in the NCR from local and regional M w 3.5, 4.5, and 5.5 earthquakes (R < 100 km). These predictions differ from those given by GMPEs which are based on recordings from earthquakes along the Himalayan arc. For example, observed PGA and PGV during the Rohtak earthquake of 5 March 2012 ( M w 4.5) are compared with the predictions of our stochastic simulation and the recently constructed GMPE by Srinagesh et al. (2021) in the bottom frames of Fig. 10 . Even though the new GPME does better than others (not shown in the figure) , it is clear that the attenuation relations constructed for the Himalayan arc earthquakes may not appropriate to predict ground motion from earthquakes in the NCR. Stochastic simulations using parameters developed in this study are likely to provide more reliable ground motion estimations. 4. Some studies have employed the empirical Green's function (EGF) technique to simulate ground motion in the NCR during larger events (Singh et al. 2002; Bansal et al. 2009; Mittal et al. 2015) . EGF simulation has the advantage that the site and path effects are automatically included in the calculations. The method is not applicable if the locations of the postulated earthquake and the site of interest differ from the corresponding locations for the EGF event. Both methods assume that the point-source, the far-field approximation is valid. In the epicentral region, this approximation may be grossly violated for M w > 5.5 earthquakes. 5. As mentioned earlier, there is no evidence of instrumentally recorded M w > 5 events in the NCR in the last 120 years. Nevertheless, since there are many mapped faults and lineaments in the region, earthquakes with M w ≥ 6 can't be ruled out. For seismic hazard estimation, however, it is critical to establish the recurrence period of such events. Synthesis of ground motion from such earthquakes will also require due consideration of finite fault and near-field effects. A model for the shape of the Fourier amplitude spectrum of acceleration at high frequencies Determination of seismic wave attenuation in Delhi, India, towards quantification of regional seismic hazard Source study of two small earthquakes of Delhi, India, and estimation of ground motion from future moderate, local events National seismological network in India for real-time earthquake monitoring Delineation of the causative fault of recent earthquakes A holistic seismotectonic model of Delhi region Recent advances in seismic instrumentation and data interpretation in India Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra Simulation of ground motion using the stochastic method Tectonic stress and the spectra of seismic shear waves from earthquakes Lapse time and frequency-dependent coda wave attenuation for Delhi and its surrounding regions Microzonation of earthquake hazard in Greater Delhi area Advances in Indian earthquake engineering and seismology: contributions in honour of Jai Krishna Estimation of site effects in Delhi using standard spectral ratio Site effects estimation in Delhi from the Indian strong motion instrumentation network Analysis of ground motion in Delhi from earthquakes recorded by strong motion network Mw 4.6 earthquake and scenario greater magnitude earthquake in the proximity of Delhi Estimation of coda wave attenuation for the national capital region, Delhi, India using local earthquakes Seismic microzonation of Delhi for ground-shaking site effects High-resolution early earthquake damage assessment system for Mexico City based on a single-station A widely felt tremor (ML 3.5) of 12 April 2020 in and around NCT Delhi in the backdrop of prevailing COVID-19 pandemic lockdown: analysis and observations Site-specific microzonation study in Delhi metropolitan city by 2-D modelling of SH and P-SV waves A review of the seismicity and seismotectonics of Delhi and adjoining areas Simulation of strong ground motion due to active Sohna fault in Delhi, national capital region (NCR) of India: an implication for imminent plausible seismic hazard Attenuation relationship for estimation of peak ground horizontal acceleration using data from strong-motion arrays in India Ground-motion prediction equations based on data from the Himalayan and Zagros regions Frequency dependent attenuation of seismic waves for Delhi and surrounding area Seismotectonic implications of Delhi region through fault plane solutions of some recent earthquakes Advances in Indian earthquake engineering and seismology: contributions in honour of Jai Krishna Pandemic lockdown sensitizes New Delhi to earthquake risk Spectral attenuation of SH waves along the Imperial fault Attenuation relations for strong seismic ground motion in Himalayan Region earthquake (Mw=5.8) and estimation of ground motion from future earthquakes in the Indian shield region Ground motion in Delhi from future large/great earthquakes in the central seismic gap of the Himalayan Arc Q of the Indian shield ): implications for seismic hazard The Delhi 1960 earthquake: epicentre, depth and magnitude Ground motion prediction equation for earthquakes along the Western Himalayan arc Groundwater extraction-induced seismicity around Delhi region Seismicity of Delhi and its surrounding region Acknowledgements Authors are thankful to National Centre for Seismology, Ministry of Earth Sciences for providing necessary support to carry out the study. SKS acknowledge the partial research support by DGAPA, UNAM project IN108221. BKB also acknowledge the support extended by Indian Institute of Technology, Delhi.Data availability Data used in the present study are collected from the published sources (referred in MS) and presented in the text. The earthquake data are available from National Centre for Seismology and could be requested through the website (https:// seismo. gov. in/). The authors declare no conflict of interest.