Next Article in Journal
GNSS-Assisted Low-Cost Vision-Based Observation System for Deformation Monitoring
Previous Article in Journal
4D-Flow MRI Characterization of Pulmonary Flow in Repaired Tetralogy of Fallot
Previous Article in Special Issue
Deep Learning Based Urban Building Coverage Ratio Estimation Focusing on Rapid Urbanization Areas
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Unified Earthquake Catalog for North East India

1
School of Computing and IT, Manipal University Jaipur, Jaipur 303007, India
2
Computer Science Department, Universidad Catolica del Norte, Antofagasta 1249004, Chile
3
Earthquake Engineering Department, Indian Institute of Technology Roorkee, Roorkee 247667, India
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(5), 2812; https://doi.org/10.3390/app13052812
Submission received: 19 December 2022 / Revised: 15 January 2023 / Accepted: 24 January 2023 / Published: 22 February 2023
(This article belongs to the Special Issue Applied Science for Urban and Rural Planning)

Abstract

:
Northeast India is one of the world’s most seismically active regions. The event data included in this research for the period 1737–2012 is mostly obtained from worldwide database agencies such as ISC, NEIC, and GCMT. Historical seismicity is collected from published and unpublished documents and some earthquake events are collected from the Indian Meteorological Department Bulletins. As the Mw scale is developed and validated in the southern California region and overestimates the smaller magnitude earthquakes, therefore, recent literature suggested an improved version of the seismic moment magnitude scale (Mwg) applicable for the entire globe considering both long- and short-period frequency-spectra using modern instrumental data. To update the earthquake catalog of Northeast India, we prepared empirical relationships between different magnitudes to Mwg using robust statistical General Orthogonal Regression. A procedure is also suggested for converting different earthquake sizes towards seismic moment scale. The Magnitude of Completeness (Mc) and the Gutenberg–Richter (GR) recurrence parameter values for the declustered homogenized catalog in four time periods, namely 1737–1963, 1964–1990, 1964–2000, and 1964–2012, have been computed. Our analysis suggests that the use of the Mwg scale improves seismicity parameters ‘b’ up to 30%, ‘a’ up to 17%, and ‘Mc’ up to 18% for the Northeast India region. A complete unified earthquake catalog in terms of advanced seismic moment magnitude scale could help understand seismicity and earthquake engineering studies of the region.

1. Introduction

Earthquakes are inherently complex phenomena, and because of differences in equipment characteristics and station coverage used to record seismic waves at different epicentral distances, most magnitude estimations in space and time are subject to measurement errors. Several earthquake magnitude scales have been used in seismic catalogs to represent the earthquake size, such as ML (Local Magnitude), mb (Body Wave Magnitude), Ms (Surface Wave Magnitude), and Mw (Moment Magnitude). Recent studies (e.g., [1,2,3,4,5]) show that the Mw scale has some serious drawbacks in representing earthquake size. Many authors stated that the Mw scale is not suitable for frequencies applicable for engineering importance (e.g., [3,6,7]). The Mw scale underestimates large earthquakes such as the 11 March 2011—Tohoku-Oki earthquake and the 26 December 2004—Sumatra earthquakes (e.g., [1,3,4,5]) and overestimates the smaller earthquakes (e.g., [1]). Although the Mw scale is expressed in terms of the moment (Mo), however, it is mainly based on surface waves; therefore, it is not a good representation of the seismic source. As the philosophy of the Mw scale [8] is based on surface waves and, therefore, it is inappropriate for deeper earthquakes. Furthermore, the Mw scale is mainly derived and validated based on Southern California tectonics [8]. The Mw scale given by Ref. [8] is valid for global earthquakes in the magnitude range >7.5 because the Mw scale is based on bigger magnitude earthquakes at the global level. As the seismic moment is mainly derived from surface waves (e.g., [9]), the formulation of the Mw scale ( i.e., ´2/3’ and ‘10.7’) also originated from surface wave magnitude. Thus, the Mw scale relates to long-period seismic spectra. Furthermore, the Mw scale was developed based on constant value for tectonic effect for shallow earthquakes. There was no statistical validation while developing the Mw scale, only a comparison of the Mw scale with some Southern Californian earthquakes [8]. Ref. [1] provide a magnitude scale mainly considering global tectonics. As the Mwg scale is developed by considering p-wave along with consideration of seismic moment M0, and M0 is computed from the surface wave trend of the seismic signal ([9,10]). Thus, the Mwg scale covers both high and low-frequency spectra of the seismic signal. The Mwg scale is preferable to the Mw scale mainly for the following reasons: (1) the Mw scales underestimate the bigger earthquake and overestimate the smaller and intermediate earthquakes (e.g., [1,5]); (2) the Mwg scale is based on directly observed Mo whereas the Mw formulation in terms of M0 was not derived based on observed M0 but rather based on substitution assuming constant values; (3) In the case of Mwg, the proper representation of global seismicity has been considered for minor, moderate, and large earthquakes, and on the other hand, Mw is based on the seismicity of Southern California, particularly for minor and intermediate occurrences; (4) Mw is not a good estimate for high-frequency ground motions, which are critical for estimating the potential shaking damage of earthquakes as Mw is defined from very long period spectral amplitudes. Mwg, on the other hand, is calculated from low- and high-frequency seismic spectra and thus fills the gap left by Mw.
Uncertainties associated with distinct magnitude scales play a significant role during magnitude conversions into a single magnitude scale. Regression relationships are used to transform diverse magnitudes into seismic moment scales to construct a homogeneous earthquake database (Mwg). General orthogonal regression (GOR) methodology is better suited for homogenizing seismicity catalogs to create regression relationships between distinct magnitude classes [11,12,13,14,15,16,17,18,19]. Various scaling relationships have been developed to homogenize an earthquake database for Northeast India [5,11,13,14,15,16,17,18,19,20,21,22,23,24]. Ref. [11] used GOR to convert body and surface wave magnitude scales into Mw, and substantial dispersion was found in the translation of mb,ISC into Mw. Ref. [21] made the conversion of mb and Ms into Mw using the least-square regression (SLR) approach. Refs. [13,14] performed regression analysis for mb to Mw and Ms to Mw using SLR, GOR, and inverted standard regression (ISR). Ref. [5] recently developed empirical relationships between different magnitudes into Mw using an improved version of the GOR technique (GOR1). Still date, the developed empirical relationships are in terms of Mw for the Northeast India region. This is the first study that gives a homogenous earthquake database for the Northeast India region using an improved magnitude scale. An enhanced GOR approach (GOR1) is applied in this investigation, as detailed in Ref. [20], for establishing regression relationships from different magnitudes to Mwg on a regional basis. For the period 1737–2012, a dataset of 9968 earthquake occurrences in the magnitude range 1.6–8.7 corresponding to the study area (lat. 20–30° and long. 87–98°) is employed. GOR relationships for conversion of mb and Ms into Mwg are developed using regional datasets. The prepared homogenous catalog will be useful not only for seismic risk assessment but also for other seismological applications. Following the unification of the earthquake catalog towards seismic moment magnitude scale Mwg, declustering has been performed, and its completeness is assessed in the sections that follow.

2. Seismicity of North East India

Figure 1 displays a seismotectonic map indicating the epicenters of Mwg ≥ 2 occurrences from 1737 to 2012, as well as the geological aspects of the NE India region ranging from 20–30° N latitude to 87–98° E longitude. The Shillong earthquake of 12 June 1897 (Mwg = 8.9) and the Assam earthquake of 15 August 1950 (Mwg = 8.6) are the two recent big earthquakes in this area. Furthermore, a succession of significant earthquakes (Mwg ≥ 7.0) has occurred in this area, resulting in the loss of lives and the damage of property. Ref. [25] classified the area into four primary seismogenic root zones: the Eastern Syntax (zone I), the Arakan–Yoma Subduction Belt (zone II), the Shillong Plateau (zone III), and the Himalayan Frontal Thrusts Main Central Thrust (MCT) and Main Boundary Thrust (MBT) (zone IV). Seismicity in this region is caused by the collision of the Indian Plate and Tibet in the north and the Burmese landmass in the east. As a result of such collisions, the Himalayan Thrust system was formed in the north, the Arakan–Yoma Mountain Arc, the Naga Hills, the Tripura Folded Belt in the east, and the Shillong Plateau rose.
The Shillong massif stands out as a plateau with an average elevation of 1500 m near the basin’s SW outlet. The seismic zones of Ref. [25] are further classified into nine seismogenic zones, as illustrated in Table 1 and Figure 2 [20]. The subdivision focuses on tectonic and geological faults, focal mechanism solutions, and the geographical distribution of earthquake events [20].

3. Data

The body and surface wave magnitudes are considered by ISC (International Seismological Center; UK, http://www.isc.ac.uk/search/Bulletin (last accessed on 18 August 2012)), and the seismic moment is considered with GCMT (Global Centroid Moment Tensor database http://www.globalcmt.org/CMTsearch.html (last accessed on 18 October 2012)). Data of 9968 earthquake events from 1737 to 2012 are obtained from several sources for the examined region (e.g., ISC, NEIC, GCMT, IMD: Indian Meteorological Department). For historical seismicity from 1897 to 1962, data are obtained from Ref. [26].

4. Regression Analysis for Magnitude Conversion

General Orthogonal Regression (GOR) is a method for establishing a connection between two variables in which measurement errors for both variables are considered. The method for carrying out the GOR is defined by many authors [13,16,18,27,28,29,30] and is not included here. A detailed description has been given in Appendix A.
The GOR equation is in the form of M Y = β 0 + β 1 M x   but has been used in a different form M Y = β 0 + β 1 m x and this practice is referred to as GOR2 (see Appendix A; various notations used in GOR methodology are explained in Appendix B). Note that M x m x ,   M x is the theoretical true point on the GOR line and m x is the observed data point. Replacing m x in place of M x in the GOR relation leads to biased estimates of the dependent variable. Many investigators pointed out the limitations of GOR2, Ref. [19] stated that GOR2 provides overestimated slope as an error variance ratio (η = σ ε 2 σ δ 2 ) does not encounter equation errors. Ref. [19] provided a method for encountering equation error and minimization of overestimation of slope for the GOR2 methodology. Ref. [19] proved theoretically and synthetically that GOR1 provides better estimates than GOR2 and SLR (Standard Least Squared Regression). GOR1 provides the highest accuracy in dependent variable estimations as compared to the conventional GOR2 and SLR approaches. Given all these reasons, we are using GOR1 for our analysis.

4.1. Surface Wave and Body Wave Conversions

The GOR1 methodology has been used to convert magnitudes such as Ms and mb to Mwg,GCMT. The GOR1 relationship for Ms,ISC to Mwg,GCMT has been developed in the range 4.1 ≤ Ms,ISC ≤ 6.1, using 93 case data and assuming η = 0.6, is obtained as follows:
Mwg,GCMT = 0.680(±0.002) Ms,ISC + 1.69(±0.08).
Rxy = 0.94, RMSE = 0.094, n = 93
The regression plot for Ms,ISC to Mwg,GCMT has been shown in Figure 3. The above obtained GOR1 relationship is found with the lowest error values in terms of slope, intercept, standard deviation, and root means square error as compared to SLR and GOR2 (Table 2 and Table 3, Figure 4). The correlation coefficient value (Rxy) obtained using the GOR1 methodology shows significant improvement as compared to GOR2 and SLR (Table 2 and Table 3). The maximum difference between Mwg estimation using Ms,ISC, and the corresponding Mw estimation of Ref. [5] is found to be 0.4.
Similarly, to convert Ms,NEIC to Mwg,GCMT, we are taking into account 57 earthquake events in the range 4.2 ≤ Ms,NEIC ≤ 6.1, and the GOR1 relationship is shown below.
Mwg,GCMT = 0.771(±0.003)MS,NEIC + 1.193(±0.126).
Rxy = 0.94, RMSE = 0.092, n = 57
The regression plots for the relationships between Ms,NEIC, and Mwg are shown in Figure 5. We observed GOR1 methodology for conversion of Ms,NEIC towards Mwg shows significant improvement in terms of error of slope, intercept, RMSE and Rxy as compared to the other two methods (Table 2 and Table 3 and Figure 6).
For changeover of mb,ISC to Mwg,GCMT for magnitude range 4.8 ≤ mb,ISC ≤ 6.1, and mb,NEIC to Mwg,GCMT for magnitude range 4.8 ≤ mb,NEIC ≤ 6.1, the same methodology has been adopted by datasets of 116 and 106 for the period 1964–2012, respectively. The GOR1 relationship between Mwg,GCMT, and mb,ISC is obtained using η = 0.2 and is given below. The regression plot is shown in Figure 7.
Mwg,GCMT = 1.19(±0.014)mb,ISC − 1.190(±0.347).
Rxy = 0.74, RMSE = 0.224, n = 116.
The maximum difference between Mwg (using Equation (3)) and Mw (using the corresponding equation of Ref. [5]) is found to be 0.4. The GOR1 methodology for conversion of mb,ISC to Mwg shows lower errors in slope and intercept as compared to SLR and GOR2. The GOR1 method provides significant improvement in the correlation coefficient (Rxy) and standard error (RMSE) as compared to GOR2 and SLR approaches (Table 2 and Table 3 and Figure 8).
The relationship between mb,NEIC, and Mwg with η equal to 0.2 is stated as
Mwg,GCMT = 1.21(±0.016)mb,NEIC − 1.391(±0.390).
Rxy = 0.72, RMSE = 0.201, n = 106.
The regression plot for mb,NEIC to Mwg,GCMT has been shown in Figure 9. The difference between Mwg estimation using Equation (8) and the corresponding Mw of Ref. [5] is found to be 0.4. Furthermore, the GOR1 methodology shows significant improvement in lowering errors of slope and intercept (Table 2 and Table 3 and Figure 10). GOR1 methodology also shows improvement in Rxy and RMSE values (Table 2 and Table 3 and Figure 10) for conversion of mb,NEIC to Mwg.

4.2. Local Magnitude into Mwg

For the period 1976–2005, the relationship between local magnitude and seismic moment magnitude is determined using 100 earthquakes in Northeast India. The derived GOR1 relationship using η = 1 is given below. The required plot for the corresponding relationship is shown in Figure 11:
Mwg,GCMT = 1.31 (±0.005) ML − 1.890 (±0.25), 5 ≤ ML ≤ 6.6
Rxy = 0.89, RMSE = 0.164, n = 100.
The maximum difference between Mw estimations from Ref. [5] and Mwg estimations in the present study for local magnitude conversion is found to be 0.3. This difference can be higher than 0.3 while considering a larger magnitude range. The GOR1 relationship for ML to Mwg has the highest accuracy in terms of the uncertainty of the regression coefficients when compared to SLR and GOR2 (Figure 12 and Table 2 and Table 3).

4.3. Duration Magnitudes into Mwg

Based on 376 global data earthquakes from the ISC database, the relationship between duration magnitude (MD) and Mwg is derived using GOR1 methodology with η = 1 and is given below. The plot of the regression is shown in Figure 13.
The relationship between MD and Mwg using η = 1 is stated as
Mwg,GCMT = 0.821(±0.002)MD − 0.832(±0.102).
Rxy = 0.81, RMSE = 0.140, n = 376.
For the conversion of MD to Mwg, it is also found that the GOR1 method has the highest accuracy compared to SLR and GOR2 (Table 2 and Table 3 and Figure 14).

4.4. Intensity Conversion Relation

Considering 29 earthquakes in India and adjacent areas from 1897 to 2016, a magnitude–intensity GOR1 relationship has been developed, with independent MMI (I0) in the range 5–12 and seismic moment magnitude (Mwg) determined from several sources, as follows:
Mwg,GCMT = 0.48 Imax + 3.07(0.1)
Figure 15 depicts the plot of the Intensity relationship.
Ref. [5] derived the GOR1 relationships between Mw and MMI for the same datasets and found a lower slope (0.389) than the present study (0.48). The maximum difference between Mw and Mwg estimations using MMI is found to be 0.5, which may lead to serious biased in the preparation of seismicity parameters. It is found that the GOR1 relationship derived between MMI and Mwg has the highest accuracy when compared to SLR and GOR2 (Table 2 and Table 3).
The GOR1 relationships derived above for the study region are compared to the similar GOR1 relationships derived by Ref. [5] for NE India and its bordering region. It is found that the regression relationships derived in this work for the NE India area are not the same as those defined by Ref. [5], as in this paper, the regression relationship is for the seismic moment magnitude (Mwg) while Ref. [5] is for the moment magnitude (Mw). For easy reference to the reader, we are reproducing the GOR1 relationships of Ref. [5] as below:
Mw,GCMT = 0.615 Ms,ISC + 2.32, 4.1 ≤ Ms,ISC ≤ 6.1
Mw,GCMT = 0.699 Ms,NEIC + 1.878, 4.2 ≤ Ms,NEIC ≤ 6.1
Mw,GCMT = 1.084 mb,ISC − 0.3106, 4.8 ≤ mb,ISC ≤ 6
Mw,GCMT = 1.104mb,NEIC − 0.495, 4.8 ≤ mb,NEIC ≤ 6.1
Mw,GCMT = 1.193ML − 0.943, 5.0 ≤ ML ≤ 6.6
Mw,GCMT = 0.742MD + 1.565, 4.2 ≤ MD ≤ 6.8
Appendix B contains abbreviations for various magnitude scales. A procedure for conversion of different magnitudes towards seismic moment scale Mwg has been suggested in Appendix C. Homogeneous earthquake catalog in terms of Mwg has been reported in Appendix A. A full data catalog can be obtained from the corresponding author.

5. Declustering of the Catalog

In general, the earthquake catalog is composed of foreshocks, mainshocks, and aftershocks Ref. [32]. Foreshocks and aftershocks should be excluded from the catalog for the evaluation of seismic hazards because they are dependent events. Several strategies for declustering a catalog have been proposed (e.g., [33,34,35,36]). We use a space and time window technique of Ref. [36] procedure for declustering the catalog (Figure 16). After declustering, there are 942 earthquake clusters, a total of 2231 (22.381%) occurrences eliminated from the homogenized catalog of 9968 events using the GOR1 method for the period 1737–2012. The seismic moment released by the clusters is about 3.7554% of the total seismic moment of the catalog.
After declustering of the homogenized catalog of Mw in Ref. [5], there are 1232 earthquake clusters, a total of 3178 (31.882%) occurrences eliminated from the total catalog of 9968 earthquake events for the period 1737–2012. The seismic moment released by the clusters is about 14.313% of the total seismic moment of the catalog. Thus, we can see that the number of clusters in Mw is more than in Mwg.

6. The Magnitude of Completeness (Mc)

Declustered earthquake catalogs for the periods 1737–1963, 1964–1990, 1964–2000, and 1964–2012 have been studied to distinguish temporal differences in earthquake happenings in NE India. Four graphs, corresponding to four different periods, have been plotted to show the relationships between various bin of magnitudes and the corresponding cumulative number of events having an earthquake magnitude greater than the corresponding magnitude of completeness (Figure 17).
The Magnitude of Completeness Mc has been calculated by the EMR method [37] for various catalog times using the ZMAP program. Mc values are seen to decline with the introduction of the newest information throughout time (Table 4). As the detection threshold for the sample area increased from 1964 onwards, the general trend of the ‘b’ value continuously decreased over time, resulting in the recording of a greater number of smaller magnitude earthquakes in proportion to large magnitude occurrences, as seen in Table 4.
In the Mw catalog, the EMR approach is used to determine Mc for various catalog time periods using the ZMAP application. Mc values are decreasing w.r.t increasing time periods. From 1964 onwards, the detection threshold for the sample region is constantly grown over time, resulting in the recording of a higher number of lower magnitude earthquakes to big magnitude events. The value of Mc for the Mw catalog is more than the Mwg catalog (Table 5). This shows that more earthquakes within a certain region are more complete on the Mwg scale than the Mw scale (Table 5).
Furthermore, the magnitudes of ‘Mc’, ‘b,’ and ‘a’ values for all the nine seismogenic zones have been determined. Only data on occurrences from 1964 to 2012 are used to calculate Mc values for distinct seismic zones, while the whole catalog for the period 1737–2012 is used for the estimate of ‘b’ and ‘a’ values. Table 6 displays the seismicity parameters for the various seismogenic zones.

7. Seismogenic Zones

The highest Mc value obtained is 4.1 for seismogenic zone IV, while the least Mc value found is 3 for seismogenic zone VI (Table 6). According to Table 6, the Mc value valid for the whole Northeast India area may be designated as 4.1 for the catalog term 1964–2012.

8. Discussion and Conclusions

A homogenous earthquake catalog is of critical importance to understanding the seismicity of a seismic region. A total of 9969 events during the period 1737–2012 for NE India have been considered in this study. As the Mw scale was mainly derived and validated for the Southern California region, therefore, a globally valid seismic moment magnitude scale Mwg is reported in recent literature. Improved GOR (GOR1) relationships for converting mb and MS to Mwg,GCMT has been developed for the study region. For converting surface wave magnitudes to Mwg magnitudes, conversion relations have been developed for the ranges.
4.1 ≤ Ms,ISC ≤ 6.1 and 4.2 ≤ Ms,NEIC ≤ 6.1, based on data 93 and 57 events, respectively. Similarly, Body wave magnitudes (mb) are converted into Mwg scales following the GOR1 methodology. The mb,ISC to Mwg,GCMT conversion relationship has been developed using 116 event data for magnitude range 4.8 ≤ mb,ISC ≤ 6.1, while the mb,NEIC to Mwg,GCMT conversion relationship is derived for magnitude range 4.8 ≤ mb,NEIC ≤ 6.1 based on 106 events.
A significant difference between Mw and Mwg estimations has been observed from various observed magnitudes such as mb, MS, ML and MD. Therefore, these differences in the seismic moment magnitude scales will lead to serious biased in the seismicity parameters and, consequently, in seismic hazard results (Table 5, [31]).
Only MMI intensity 5 and above data were utilized to construct an empirical relationship between intensity (Imax) and Seismic moment magnitude Mwg,GCMT. In total, 29 MMI intensities and Mwg,GCMT data pairs are included from the entirety of India. The MMI empirical relationship for intensities 5 and higher is mostly compatible with the Indian seismic zoning chart, which shows the range of values associated with the major seismic zones according to the seismic code (IS, 2002) released by the Bureau of Indian Standards (BIS) [38]. Developed intensity-seismic moment magnitude relationship should be used to convert historical earthquakes to DMS (Das Magnitude Scale) scale when magnitude information is unavailable.
The declustering of the homogenized catalog for the time being 1737–2012 is carried out using the Ref. [36] procedure, and there is a reduction of 2231 events (22.38%) during this process. The entire homogenized earthquake catalog has been classified into four catalog periods, namely 1737–1963, 1964–1990, 1964–2000, and 1964–2012. For the four catalog periods, the magnitude of the completeness values is obtained using the EMR process. The Mc value declined from 5.6 for 1737–1964 to 3.9 for 1964–2012, as projected by a rise in instrumentation for this area beginning in 1964. Improving Mc has a major impact on the computation of the ‘b’ value and the assessment of the seismic hazard for a given location.
Several earthquakes of Mwg ≥ 7.0 occurred in this area between 1737 and 2012. The maximum time elapsed between the occurrences of a major earthquake (Mwg ≥ 7.0) was 18 years (1970–1988) in the 1737–2012 catalog. The last such occurrence happened in the area in 1988. The period of low seismic activity appears to have begun after 1988 and continues to this day.
The study area has been subdivided into nine seismogenic zones considering Mwg based earthquake catalogs, focal mechanisms, and fault types. ‘Mc’, ‘b,’ and ‘a’ values have been determined for each of these areas. The catalog data for the period 1897–2012 show that there have been no earthquakes of magnitude ≥ 7.0 in zones I, VI, or IX. Because zone IX is located between IV and VII zones in which big earthquakes have occurred, the probability of a big earthquake occurring in this zone is low in the immediate future. Based on these findings, a period of quiescence has been identified in seismogeneous zones I and VI for large-scale earthquakes.
A complete and consistent unified seismic catalog has been developed in terms of Mwg following a robust statistical procedure that could help to understand the seismicity of the region in a better way. Preparing a homogenized earthquake catalog by changing the original magnitude scales into seismic moment magnitude scale Mwg, an obstacle has been removed for seismic hazard assessment of the study region. Our analysis suggests that the use of the Mwg scale improves seismicity parameters ‘b’ up to 30%, ‘a’ up to 17%, and ‘Mc’ up to 38% for the Northeast India region (Table 5). Hence, the variations in these parameters may have a significant impact on the seismic hazard results. Therefore, the use of the Mwg scale is recommended for all practical cases.

Author Contributions

Conceptualization, R.D. and P.; data curation. R.D. and P.; formal analysis, R.D., P. and C.M.; investigation, R.D. and P.; methodology, R.D. and P.; resources, R.D.; software, R.D.; validation, R.D. and P.; visualization, R.D. and P.; writing—original draft, R.D., P., C.M., S.J. and T.B.; writing—review and editing, R.D., P., C.M., S.J. and T.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by FONDECYT Grant 11200618.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Body and surface wave magnitudes of earthquake data have been collected from ISC (International Seismological Center) U.K. (http://www.isc.ac.uk/search/Bulletin) (last accessed August 2012) and moment magnitudes have been collected from GCMT (Global Centroid Moment Tensor database, http://www.globalcmt.org/CMTsearch.html) (last accessed October 2012). Earthquake data of 9969 events for the time period 1897–2012 were compiled from various databases agencies (e.g. ISC, NEIC, GCMT, IMD, NEIST). Historical seismicity data were considered during 1897 to 1962 from Ref. [26].

Acknowledgments

We thank the Reviewers for their constructive comments, which largely improve the quality of this contribution.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Let Mx and My be the true values, and mx and my be the derived values with δ and ε as magnitudes of errors in measuring for the independent and dependent variables. Then we can write:
mx = Mx + δ,
my = My + ε,
as well as the regression-like model:
my = α + βMx + ε,
with My = α + βMx, where α and β are the slopes and intercepts of a linear relationship between genuine values
Mx, ε, and δ are assumed to be distributed in the relationships on a regular and independent basis.
The observed value covariances σ m y   2 , σ m x m y and σ m x 2 are given by:
σ m y 2 = β 2 σ M x 2 + σ ε 2 ,
σ m x m y = β σ M x 2 ,
σ m x 2 = σ M x 2 + σ δ 2 ,
where the error variance ratio:
η = σ ε 2 σ δ 2
If, s m y 2 , s m x 2 and s m x m y are the sample covariances of my, mx and between my and mx, then:
s m y 2 = β ^ 2 σ ^ M x 2 + η σ ^ δ 2 ,
s m x m y   = β ^ σ ^ M x 2 ,
s m x 2 = σ ^ M x 2 + σ ^ δ 2
The estimators β ^ 2 , σ ^ M x 2 and σ ^ δ 2 maybe simply determined using the above simultaneous Equations (A8)–(A10). For example, the elimination of σ ^ M x 2 and σ ^ δ 2 , we get the quadratic equation:
β ^ 2 s m x m y β ^ s m y 2 η s m x 2 η s m x m y ,
Which yields
β ^ = s m y 2 η s m x 2 + s m y 2 η s m x 2 2 + 4 η s m x m y 2 2 s m x m y
The σ ^ M x 2 and σ ^ δ 2 are similarly derived as follows:
σ ^ M x 2 = s m y 2 η s m x 2 2 + 4 η s m x m y 2 s m y 2 η s m x 2 2 η ,
σ ^ δ   2 = s m y 2 + η s m x 2 s m y 2 η s m x 2 2 + 4 η s m x m y 2 2 η
The estimator for α can be obtained from the relation:
α ^ = m ¯ y β ^ m ¯ x
where m ¯ x and m ¯ y are the averages of the observed values.
Using Ref. [29] equations, the estimated variances of regression parameters β   ^ and α ^ , may be expressed as follows:
σ ^ β ^ 2 = σ ^ M x 2 n 1 η + β ^ 2 σ ^ δ 2 + σ ^ δ 2 2 n 1 η + β ^ 2 n 2 β ^ σ ^ δ 2 2 n 2 n 1 σ ^ M x 2 2
and
σ ^ α 2 = n 1 η + β ^ 2 σ ^ δ 2 n n 2 + m ¯ x 2 σ ^ β ^ 2 ,
where n is the sample size.

Appendix B. Descriptions of Various Notations Used in Our Study

NotationDetail
mb,ISC:Body Wave Magnitude from ISC
mb,NEIC:Body Wave Magnitude from NEIC
MS,ISC:Surface Wave Magnitude from ISC
MS,NEIC:Surface Wave Magnitude from NEIC
ML,IMD:Local Magnitude from Indian Meteorological Department
ML:Local Magnitude Scale
MD,NEIC:Duration Magnitude from NEIC
M0:Seismic Moment
MW:Moment Magnitude was given by Ref. [8]
Mwg:Seismic moment magnitude or Das magnitude Scale is given by Ref. [1]
Mwg,GCMT:Seismic moment magnitude determined by GCMT
MW,NEIC:Moment Magnitude determined by NEIC
MMI:Modified Mercalli Scale
GOR:Conventional General Orthogonal Regression
GOR1:General Orthogonal Regression gave by Ref. [19]
GOR2:Conventional General Orthogonal Regression
SLR:Standard Least square Regression
Mx:Theoretical True value corresponding to the observed independent variable
My:Theoretical True value corresponding to the observed dependent variable
η:Error Variance Ratio
MSE:Mean Square Error
MAE:Mean Average Error
Rxy:Correlation Coefficient
RMSE:Root Mean Square Error

Appendix C. A Scheme for Conversions of Different Magnitude into Mwg,GCMT

Code 1: Proxy Mwg,GCMT estimates from Mwg,GCMT and MS,NEIC
Code 2: Proxy Mwg,GCMT estimates from Mwg,GCMT and MS,ISC
Code 3: Proxy Mwg,GCMT estimates from Mwg,GCMT and ML,IMD
Code 4: Proxy Mwg,GCMT estimates from Mwg,GCMT and mb,NEIC
Code 5: Proxy Mwg,GCMT estimates from Mwg,GCMT and mb,ISC
Code 6: Proxy Mwg,GCMT estimates from Mwg,GCMT and MD,NEIC
Code 7: Proxy Mwg,GCMT estimates from Mwg,GCMT, and MS
Code 8: Proxy Mwg,GCMT estimates from Mwg,GCMT and MW,NEIC

References

  1. Das, R.; Sharma, M.; Choudhury, D.; Gonzalez, G. A Seismic Moment Magnitude Scale. Bull. Seismol. Soc. Am. 2019, 109, 1542–1555. [Google Scholar] [CrossRef]
  2. Beresnev, I.A. The reality of scaling law of earthquake-source spectra? J. Seismol. 2009, 13, 433–436. [Google Scholar] [CrossRef]
  3. Bormann, P.; Di Giacomo, D. The moment magnitude Mw and the energy magnitude Me: Common roots and differences. J. Seismol. 2011, 15, 411–427. [Google Scholar] [CrossRef]
  4. Bormann, P.; Saul, J. A Fast, Non-saturating Magnitude Estimator for Great Earthquakes. Seismol. Res. Lett. 2009, 80, 808–816. [Google Scholar] [CrossRef]
  5. Das, R.; Meneses, C. A unified moment magnitude earthquake catalog for Northeast India. Geomat. Nat. Hazards Risk 2021, 12, 167–180. [Google Scholar] [CrossRef]
  6. Choy, G.L.; Boatwright, J.L. Global patterns of radiated seismic energy and apparent stress. J. Geophys. Res. Solid Earth 1995, 100, 18205–18228. [Google Scholar] [CrossRef] [Green Version]
  7. Choy, G.L.; Kirby, S.H. Apparent stress, fault maturity and seismic hazard for normal-fault earthquakes at subduction zones. Geophys. J. Int. 2004, 159, 991–1012. [Google Scholar] [CrossRef] [Green Version]
  8. Hanks, T.C.; Kanamori, H. A moment magnitude scale. J. Geophys. Res. 1979, 84, 2348–2350. [Google Scholar] [CrossRef]
  9. Ekström, G.; Dziewoński, A.M.; Maternovskaya, N.N.; Nettles, M. Global seismicity of 2003: Centroid–moment-tensor solutions for 1087 earthquakes. Phys. Earth Planet. Inter. 2005, 148, 327–351. [Google Scholar] [CrossRef]
  10. Kanamori, H.; Anderson, D.L. Theoretical basis of some empirical relations in seismology. Bull. Seismol. Soc. Am. 1975, 65, 1073–1095. [Google Scholar]
  11. Thingbaijam, K.K.S.; Nath, S.K.; Yadav, A.; Raj, A.; Walling, M.Y.; Mohanty, W.K. Recent seismicity in Northeast India and its adjoining region. J. Seismol. 2008, 12, 107–123. [Google Scholar] [CrossRef]
  12. Ristau, J. Comparison of Magnitude Estimates for New Zealand Earthquakes: Moment Magnitude, Local Magnitude, and Teleseismic Body-Wave Magnitude. Bull. Seismol. Soc. Am. 2009, 99, 1841–1852. [Google Scholar] [CrossRef]
  13. Das, R.; Wason, H.R.; Sharma, M.L. Magnitude conversion to unified moment magnitude using orthogonal regression relation. J. Asian Earth Sci. 2012, 50, 44–51. [Google Scholar] [CrossRef]
  14. Das, R.; Wason, H.R.; Sharma, M.L. Homogenization of Earthquake Catalog for Northeast India and Adjoining Region. Pure Appl. Geophys. 2012, 169, 725–731. [Google Scholar] [CrossRef]
  15. Das, R.; Wason, H.R.; Sharma, M.L. General Orthogonal Regression Relations between Body-Wave and Moment Magnitudes. Seismol. Res. Lett. 2013, 84, 219–224. [Google Scholar] [CrossRef]
  16. Das, R.; Wason, H.R.; Sharma, M.L. Reply to ‘Comment on “Magnitude conversion problem using general orthogonal regression” by H. R. Wason, Ranjit Das and M. L. Sharma’ by Paolo Gasperini and Barbara Lolli. Geophys. J. Int. 2014, 196, 628–631. [Google Scholar] [CrossRef]
  17. Das, R.; Wason, H.R.; Sharma, M.L. Unbiased Estimation of Moment Magnitude from Body- and Surface-Wave Magnitudes. Bull. Seismol. Soc. Am. 2014, 104, 1802–1811. [Google Scholar] [CrossRef]
  18. Wason, H.R.; Das, R.; Sharma, M.L. Regression Relations for Magnitude Conversion for the Indian Region. In Advances in Indian Earthquake Engineering and Seismology; Sharma, M.L., Shrikhande, M., Wason, H.R., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 55–66. [Google Scholar] [CrossRef]
  19. Das, R.; Wason, H.R.; Gonzalez, G.; Sharma, M.L.; Choudhury, D.; Lindholm, C.; Roy, N.; Salazar, P. Earthquake Magnitude Conversion Problem. Bull. Seismol. Soc. Am. 2018, 108, 1995–2007. [Google Scholar] [CrossRef]
  20. Das, R.; Wason, H.R.; Sharma, M.L. Temporal and spatial variations in the magnitude of completeness for homogenized moment magnitude catalogue for northeast India. J. Earth Syst. Sci. 2012, 121, 19–28. [Google Scholar] [CrossRef] [Green Version]
  21. Yadav, R.B.S.; Bormann, P.; Rastogi, B.K.; Das, M.V.; Chopra, S. A Homogeneous and Complete Earthquake Catalog for Northeast India and the Adjoining Region. Seismol. Res. Lett. 2009, 80, 609–627. [Google Scholar] [CrossRef]
  22. Nath, S.K.; Mandal, S.; Das Adhikari, M.; Maiti, S.K. A unified earthquake catalogue for South Asia covering the period 1900–2014. Nat. Hazards 2017, 85, 1787–1810. [Google Scholar] [CrossRef]
  23. Pandey, A.K.; Chingtham, P.; Roy, P.N.S. Homogeneous earthquake catalogue for Northeast region of India using robust statistical approaches. Geomat. Nat. Hazards Risk 2017, 8, 1477–1491. [Google Scholar] [CrossRef] [Green Version]
  24. Anbazhagan, P.; Balakumar, A. Seismic magnitude conversion and its effect on seismic hazard analysis. J. Seismol. 2019, 23, 623–647. [Google Scholar] [CrossRef]
  25. Dutta, T.K. Seismicity of Assam-zones of tectonic activity. Seism. Assam Zones Tecton. Act. 1964, 2, 152–163. [Google Scholar]
  26. Gupta, H.K.; Rajendran, K.; Singh, H.N. Seismicity of Northeast India region: PART I: The database. J. Geol. Soc. India 1986, 28, 345–365. [Google Scholar]
  27. Madansky, A. The fitting of straight lines when both variables are subject to error. J. Am. Stat. Assoc. 1959, 54, 173–205. [Google Scholar] [CrossRef]
  28. Kendall, M.; Stuart, A. The Advanced Theory of Statistics; Griffin: London, UK, 1976; Volume 1, p. 102. [Google Scholar]
  29. Fuller, W.A. Measurement Error Models; John Wiley & Sons: New York, NY, USA, 2009. [Google Scholar]
  30. Carroll, R.J.; Ruppert, D. The use and misuse of orthogonal regression in linear errors-in-variables models. Am. Stat. 1996, 50, 1–6. [Google Scholar]
  31. Vanzi, I.; Marano, G.C.; Monti, G.; Nuti, C. A synthetic formulation for the Italian seismic hazard and code implications for the seismic risk. Soil Dyn. Earthq. Eng. 2015, 77, 111–122. [Google Scholar] [CrossRef]
  32. Contento, A.; Aloisio, A.; Xue, J.; Quaranta, G.; Briseghella, B.; Gardoni, P. Probabilistic axial capacity model for concrete-filled steel tubes accounting for load eccentricity and debonding. Eng. Struct. 2022, 268, 114730. [Google Scholar] [CrossRef]
  33. Utsu, T. A method for determining the value of “b” in a formula log n = a-bM showing the magnitude-frequency relation for earthquakes. Geophys. Bull. Hokkaido Univ. 1965, 13, 99–103. [Google Scholar]
  34. Gardner, J.K.; Knopoff, L. Is the sequence of earthquakes in Southern California, with aftershocks removed, Poissonian? Bull. Seismol. Soc. Am. 1974, 64, 1363–1367. [Google Scholar] [CrossRef]
  35. Reasenberg, P. Second-order moment of central California seismicity, 1969–1982. J. Geophys. Res. Solid Earth 1985, 90, 5479–5495. [Google Scholar] [CrossRef]
  36. Uhrhammer, R.A. Characteristics of northern and central California seismicity. Earthq. Notes 1986, 57, 21. [Google Scholar]
  37. Wiemer, S. A software package to analyze seismicity: ZMAP. Seismol. Res. Lett. 2001, 72, 373–382. [Google Scholar] [CrossRef]
  38. Standard, I. Criteria for earthquake resistant design of structures. Bur. Indian Stand. Part 1893, 1, 1–21. [Google Scholar]
Figure 1. Seismotectonic map exhibiting seismicity for Mwg ≥ 2 with epicenters and tectonic characteristics of the NE India region on a GIS platform.
Figure 1. Seismotectonic map exhibiting seismicity for Mwg ≥ 2 with epicenters and tectonic characteristics of the NE India region on a GIS platform.
Applsci 13 02812 g001
Figure 2. Seismogenic source zones (I-IX) are being considered for Northeast India.
Figure 2. Seismogenic source zones (I-IX) are being considered for Northeast India.
Applsci 13 02812 g002
Figure 3. The (MS,ISC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 3. The (MS,ISC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g003
Figure 4. Comparison of the uncertainty in regression parameters for (MS,ISC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 4. Comparison of the uncertainty in regression parameters for (MS,ISC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g004aApplsci 13 02812 g004b
Figure 5. The (MS,NEIC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 5. The (MS,NEIC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g005
Figure 6. Comparison of the uncertainty in regression parameters for (MS,NEIC – Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 6. Comparison of the uncertainty in regression parameters for (MS,NEIC – Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g006
Figure 7. The (mb,ISC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 7. The (mb,ISC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g007
Figure 8. Comparison of the uncertainty in regression parameters for (mb,ISC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering the function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 8. Comparison of the uncertainty in regression parameters for (mb,ISC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering the function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g008
Figure 9. The (mb,NEIC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 9. The (mb,NEIC − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g009
Figure 10. Comparison of the uncertainty in regression parameters for (mb,NEIC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 10. Comparison of the uncertainty in regression parameters for (mb,NEIC − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g010
Figure 11. The (ML − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 11. The (ML − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g011
Figure 12. Comparison of the uncertainty in regression parameters for (ML − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 12. Comparison of the uncertainty in regression parameters for (ML − Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g012
Figure 13. The (MD − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Figure 13. The (MD − Mwg) data, as well as the GOR1, GOR2, and SLR regression lines.
Applsci 13 02812 g013
Figure 14. Comparison of the uncertainty in regression parameters for (MD – Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Figure 14. Comparison of the uncertainty in regression parameters for (MD – Mwg) data pairs using GOR1, GOR2, and SLR approaches considering function of error variance ratio: (a) Correlation Coefficient determination (Rxy), (b) Root mean square error (RMSE), (c) uncertainty values of slope and (d) uncertainty of intercept.
Applsci 13 02812 g014
Figure 15. Plot showing the GOR connection for intensity scaling (MMI) to Mwg,GCMT.
Figure 15. Plot showing the GOR connection for intensity scaling (MMI) to Mwg,GCMT.
Applsci 13 02812 g015
Figure 16. Aftershocks are removed using magnitude-dependent spatial and temporal frames (dependent occurrences). The asterisks under the window lines are interpreted as dependent occurrences.
Figure 16. Aftershocks are removed using magnitude-dependent spatial and temporal frames (dependent occurrences). The asterisks under the window lines are interpreted as dependent occurrences.
Applsci 13 02812 g016
Figure 17. Four different plots showing the cumulative number of earthquake events vs. magnitudes greater than the corresponding magnitude for the following catalog time periods (a) 1737–1963, (b) 1964–1990, (c) 1964–2000, and (d) 1964–2012. In addition, these plots show Mc values determined by the EMR technique.
Figure 17. Four different plots showing the cumulative number of earthquake events vs. magnitudes greater than the corresponding magnitude for the following catalog time periods (a) 1737–1963, (b) 1964–1990, (c) 1964–2000, and (d) 1964–2012. In addition, these plots show Mc values determined by the EMR technique.
Applsci 13 02812 g017
Table 1. Different Seismogenic zones for Northeast India.
Table 1. Different Seismogenic zones for Northeast India.
Seismogenic ZoneMajor DivisionSubdivision
IIndo Burma Fault BeltNS Indo-Burma Fold Belt
IIIndo Burma Fault BeltNE-SW Indo Burma Fold Belt
IIIPlateau RegionSagging Fault Region
IVMishmi MassifNW-SE trending feature
VPlateau RegionTibetan Plateau
VIHimalayan Mountain BeltEastern MCT
VIIShillong MassifShillong Plateau
VIIIBengal BasinSylhet Fault
IXHimalayan Mountain BeltNE-SW trending Structure
Table 2. Comparisons of regression parameters of GOR1, GOR2 and SLR for Northeast India Region Datasets.
Table 2. Comparisons of regression parameters of GOR1, GOR2 and SLR for Northeast India Region Datasets.
Regression RelationMagnitude RangeSlope
(GOR1)
Intercept
(GOR1)
Slope
(GOR2)
Intercept
(GOR2)
Slope
SLR
Intercept
SLR
Rxy
GOR1
Rxy
GOR2
Rxy
SLR
RMSE
GOR1
RMSE
GOR2
RMSE
SLR
mb, ISC to Mwg
η = 0.2
4.8 ≤ mb, ISC ≤ 6.11.19−1.191.61−3.381.16−1.010.740.590.690.220.280.24
mb, NEIC to Mwg η = 0.24.8 ≤ mb, NEIC ≤ 6.11.21−1.391.68−3.891.18−1.270.720.550.680.20.260.22
Ms,ISC to Mwgη = 0.64.1 ≤ Ms, ISC ≤ 6.10.681.690.711.5250.641.890.970.880.90.080.190.18
Ms,NEIC to Mwgη = 0.64.2 ≤ Ms, NEIC ≤ 6.10.771.190.820.970.731.400.940.790.80.090.180.17
Intensity to Mwg
η = 1
5 to 120.483.070.493.000.443.360.980.610.620.130.690.68
Local Magnitude
η = 1
5.0 ≤ ML ≤ 6.61.31−1.891.47−2.761.24−1.490.890.740.770.160.250.24
Duration Magnitude η = 14.2 ≤ MD ≤ 6.80.820.831.007−0.1090.641.760.810.260.40.140.30.27
Table 3. The list of error metrics for regression parameters such as Error in Slope, Error in intercept, Mean Square Error, Mean Average Error, and Root Mean square Error corresponding to the derived regression relations for this study [16,31].
Table 3. The list of error metrics for regression parameters such as Error in Slope, Error in intercept, Mean Square Error, Mean Average Error, and Root Mean square Error corresponding to the derived regression relations for this study [16,31].
Regression RelationMagnitude RangeError in SlopeError in InterceptMSEMAERMSE
GOR1GOR2SLRGOR1GOR2SLRGOR1GOR2SLRGOR1GOR2SLRGOR1GOR2SLR
mb, ISC to Mwg
η=0.2
4.8 ≤ mb, ISC≤ 6.1±0.01±0.098±0.07± 0.35±0.49± 0.370.0490.070.060.170.220.190.220.280.24
mb, NEIC to Mwg η=0.24.8 ≤ mb, NEIC≤ 6.1±0.02±0.11±0.08± 0.37±0.60± 0.420.040.290.050.170.470.190.20.260.22
Ms,ISC to Mwg η=0.64.1 ≤ Ms, ISC ≤ 6.1±0.00±0.04±0.03±0.08±0.24± 0.180.010.040.040.070.160.150.080.190.18
Ms,NEIC to Mwg η=0.64.2 ≤ Ms, NEIC ≤ 6.1±0.00±0.05±0.05±0.13±0.326± 0.230.010.030.030.070.150.140.090.180.17
Intensity to Mwg
η=1
5 to 12±0.001±0.07±0.07±0.51±1.38± 0.500.020.440.430.10.530.520.130.690.68
Local Magnitude
η=1
5.0 ≤ ML ≤ 6.6±0.00±0.08±0.07±0.25±0.5± 0.370.020.060.560.020.110.020.160.250.24
Duration Magnitude η=14.2 ≤ MD ≤ 6.8±0.002±0.06±0.04±0.10213±0.35774± 0.2030.010.050.050.110.310.220.140.30.27
Table 4. Estimated Mc, ‘b,’ and ‘a’ values for various catalogs.
Table 4. Estimated Mc, ‘b,’ and ‘a’ values for various catalogs.
Catalog Time PeriodsMcba
1737–19635.6 ± 0.230.746.5
1964–19904.1 ± 0.170.987.56
1964–20004 ± 0.170.856.86
1964–20123.9 ± 0.230.816.83
Table 5. Difference in seismicity parameters between Mwg catalog and Mw catalog for different time periods.
Table 5. Difference in seismicity parameters between Mwg catalog and Mw catalog for different time periods.
MwgMw
Mc (1964–1990)4.1 ± 0.174.7 ± 0.13
b0.981.27
a7.568.98
Mc (1964–2000)4 ± 0.174.4 ± 0.17
b0.850.95
a6.867.57
Mc (1964–2012)3.4 ± 0.114.1 ± 0.13
b0.610.8
a5.736.96
No. of Clusters942 (22.381%)1232 (31.882%)
Table 6. Estimated ‘Mc’, ‘b,’ and ‘a’ values for various seismogenic zones.
Table 6. Estimated ‘Mc’, ‘b,’ and ‘a’ values for various seismogenic zones.
Seismogenic ZoneMc‘b’‘a’
I3.40.65.02
II3.50.665.51
III3.90.74.96
IV4.11.096.52
V3.10.53.98
VI30.453.84
VII3.10.694.87
VIII3.20.613.85
IX3.30.654.92
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pallavi; Das, R.; Joshi, S.; Meneses, C.; Biswas, T. Advanced Unified Earthquake Catalog for North East India. Appl. Sci. 2023, 13, 2812. https://doi.org/10.3390/app13052812

AMA Style

Pallavi, Das R, Joshi S, Meneses C, Biswas T. Advanced Unified Earthquake Catalog for North East India. Applied Sciences. 2023; 13(5):2812. https://doi.org/10.3390/app13052812

Chicago/Turabian Style

Pallavi, Ranjit Das, Sandeep Joshi, Claudio Meneses, and Tinku Biswas. 2023. "Advanced Unified Earthquake Catalog for North East India" Applied Sciences 13, no. 5: 2812. https://doi.org/10.3390/app13052812

APA Style

Pallavi, Das, R., Joshi, S., Meneses, C., & Biswas, T. (2023). Advanced Unified Earthquake Catalog for North East India. Applied Sciences, 13(5), 2812. https://doi.org/10.3390/app13052812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop