Spatial variation of seismic b-value is estimated in the Indo-Myanmar subduction zone of northeast (NE) India
using the homogeneous part of earthquake catalogue (1996–2015), recorded by International Seismological
Center (ISC), consisting of 895 events of magnitude MW?3.9. The study region is divided into 1°×1° square
grids and b-values are estimated at each grid by maximum likelihood method. In this study, the b-value varies
from 0.75 to 1.54 in the region. Significant variation of low b-value in the respective location may indicate high
stress accumulation in that region. Spatial variation reveals intermediate b-value anomalies around the epicenter
of the Mw=6.7 Manipur earthquake which occurred on 3rd January at 23:05 UTC (4 January 2016 at 04:35
IST). The variations of b-values are also estimated with respect to depth. The low b-value associated with the
depth range?15–55 km, which may imply crustal homogeneity and high stress accumulation in the crust. Since,
NE India lies in the seismic zone V of the country; this study can be helpful to understand seismotectonics in the
region. The study of seismicity patterns is crucial for evaluation of seismic
hazards and better understanding of different seismotectonic aspects
within a seismogenic zone (Bernard et al., 2006; Schorlemmer et al.,
2005; White et al., 2009; Sorbi et al., 2012). In the Gutenberg and
Richter (1944) frequency-magnitude distribution relation, the b-value is
an important seismic parameter that describes the characteristics of an
ensemble of earthquakes. The parameter ‘b’ is believed to be dependent
on the stress regime and tectonic character of a region (Allen et al.,
1965; Mogi, 1967; Scholz, 1968). Many studies reveal that the b-value
indicates stress distribution and level of heterogeneity in the Earth’s
crust (e.g. Mogi, 1962; Scholz, 1968; Wyss, 1973; Wiemer and
Katsumata, 1999; Schorlemmer et al., 2005; Mousavi, 2017). As such, it
can provide significant information on the seismotectonic and hazard
potential of a region. Laboratory studies show, b-value are governed by
several factors, such as, the degree of material heterogeneity (Mogi,
1962), the level of applied stress (Scholz, 1968), the degree of stress
concentration (Meredith and Atkinson, 1983), the chemical reactivity
of pore fluid (Meredith and Atkinson, 1983), and the pore fluid pressure
(Sammonds et al., 1992; Lockner and Byrlee, 1991; Power et al., 1995;
Wiemer and McNutt, 1997; Wyss et al., 1997; Wiemer et al., 1998; Jollyand McNutt, 1999; Murru et al., 1999; Sanchez et al., 2004; Bridges and
Gao, 2006; Farrell et al., 2009; Van Stiphout et al., 2009; Bachmann
et al., 2012; Mousavi, 2017). Increase in the material heterogeneity or
an increase in the thermal gradient results in high b-values (Mogi, 1962;
Warren and Latham, 1970), while b-value diminishes with rise in applied
shear stress or effective shear stress (Wyss, 1973; Scholz, 1968;
Kayal, 2008; Singh et al., 2008). The b-values systematically vary with
earthquake focal mechanism. For example, normal faulting events exhibit
the highest, thrusts events have the lowest and strike-slip events
show intermediate b-values, respectively (Schorlemmer et al., 2005).
The areas of crustal homogeneity and high stress correspond to smaller
b-values (b 1) imply crustal heterogeneity
and low stress (Bridges and Gao, 2006; Mousavi, 2017). In
volcanic region with earthquake swarms, b-values are rather high, and
can have values as high as 3.0 (McNutt, 2005). For tectonic events bvalues
are in the range of 0.8–1.2 (e.g., Schorlemmer et al., 2005;
Kamer and Hiemer, 2015). Spatial variations of b-value have been
studied by several researchers in a number of seismically active areas
(e.g., Wiemer and Katsumata, 1999; Singh et al., 2009). For NE India,
few researchers have performed the spatial variation of b-values accompanied
by depth (e.g., Bhattacharya et al., 2002; Bhattacharya and
Kayal, 2003; Khan et al., 2011). The Indo-Myanmar subduction zone of NE India lying between latitude
22–27°N and longitude 93–97°E, is tectonically one of the most
earthquake-prone zones in the world (Fig. 1). The region has experienced
several large earthquakes of magnitude M ; 6.0 in the past few
decades (Kundu and Gahalaut, 2012). Subduction of the Indian lithosphere
plate beneath the Burmese plate makes this region tectonically
very complex. In this paper, an attempt has been made to evaluate the
spatial variation, and depth distribution of b-value in the Indo-
Myanmar subduction zone before the occurrence of 3rd January 2016
Manipur earthquake of magnitude MW=6 (Gahalaut et al., 2016). The NE India (Fig. 1) is characterized by high seismic activity. The
region is flanked by the ongoing India-Asia collision in the north and
the Indo-Myanmar subduction in the east, the Brahmaputra river alluvium
in the Assam valley, the Shillong-Mikir plateau in between these
two arcs, and thick sediments of the Bengal basin in the south (Bilham
and England, 2001; Kayal et al., 2012). The Indo-Myanmar subduction zone of NE India is one of the most
seismically active fold and thrust belts in the world, which extends from
Mishmi hills to the Bay of Bengal in NE-SW, NNE-SSW and N-S directions,
where the Indian plate subducts under the Burmese micro-plate
and shows very complex deformation (Barman et al., 2014). Since
1897, the NE India has experienced 20 large (M?7.0) and two great
earthquakes (M?8.5); one on June 12, 1897 (Oldham, 1899) and the
other on August 15, 1950 (Tandon, 1954). The western Indo-Burma
fold and thrust belt is segmented into N-S blocks along the E-W transverse
zones and exhibits dextral slip between the Nagaland salient and
the Manipur recess and sinistral slip between Manipur recess and the
Tripura-Mizoram salient (Jade et al., 2007). Towards the east of the
Indo-Burma fold and thrust belt (IBFTB), the India Sunda relative plate
motion of about 36 mm/year is partitioned between the Churachandpur-
Mao fault through dextral strike-slip of 16 mm/year and the
remaining motion (?20 mm/year) is accommodated at Sagaing fault
through dextral strike-slip motion (Kundu and Gahalaut, 2013). The
Sagaing fault is one of the most prominent active strike-slip faults in the
Indo-Burmese arc. It extends over 1000 km from north to south. In this way, it touches the Andaman spreading center (ASC) towards its
southern termination (Hurukawa and Maung, 2011). The largest
earthquake that had occurred along, or near, the Myanmar portion of
this fault is the Burma earthquake of Mw ?8 (Kundu and Gahalaut,
2013). The earthquakes in the Indo-Myanmar arc are referred to as
subduction tectonics and normal, thrust and strike-slip faulting have
been reported for the region (Kayal, 1996; Kumar and Rao, 1995).
Kayal (2008) and Le Dain et al. (1984) reported that intermediate depth
earthquakes dip eastward beneath the Indo-Myanmar ranges and they
reach a depth of ?200 km. The tectonic stress is much complicated in
the region due to collision tectonics in the north and subduction tectonics
in the east. The detail of the regional stress was studied by several
researchers (Jade et al., 2007; Kundu and Gahalaut, 2012, 2013;
Gahalaut et al., 2016). In this study, an earthquake data set containing 2648 events has
been selected of different magnitude scales within Indo-Myanmar
subduction zone. The magnitude of these events ranges from 2.5 to 7.0.
The data-set was prepared from the International Seismological Center
(ISC) catalogue during the period 1964–2015. While selecting the
events, duplication of events was kept minimum. Among the magnitude
scales as prepared by ISC catalogue, MW, mb, ML, MD and MS are the
significant ones, and this entails different conversion relations. So far,
several conversion relations of magnitude have been suggested by different
authors for NE India. In this study, Das et al. (2012) and Bora
et al. (2013) have been followed to convert the seismicity catalogue
into a homogeneous one (moment magnitude, Mw). Still, there is an
apparent influence of Magnitude of completeness (MC). This is considered
to be crucial for all seismicity-based studies as it leads to reliable
result by optimizing number of events available (Wiemer and
Wyss, 2000). As per definition, Mc refers to minimum magnitude of
complete recording at which 100% of the earthquakes in a space-time
volume are detected (Wiemer and Wyss, 2000; Woessner and Weimer,
2005; Rydelek and Sacks, 1989). In general, there happens to be spatial
and temporal change in Mc which corresponds to the ability of detection
of networks augments with passage of time (Mousavi, 2017). Methods
like Entire-magnitude-range (EMR) or the best fit suggested by ‘Zmap’
(Wiemer, 2001) are very instrumental. Although, they result in high Mc,
but their accuracy is limited upto some extent. However, while conducting
a standard analysis, best combination (MC ?95, MC ?90
Maximum Curvature) method (Wiemer and Wyss, 2000) emerges to be
a good choice and it is widely accepted. For entire data, temporal
variation of MC has been estimated (Fig. 2a). A sample window of size
200 samples has been chosen based on the number of earthquakes and
their distribution with time in the catalog. From Fig. 2a, it can be seen
that MC decreases sharply from 1996. This can be seen from the cumulative
curve as well (Fig. 2b). This may be either due to network
change or the use of more stations in the recent years that resulted in
detection of more events. We have selected the part of the catalog from 1996, where cumulative curve shows constant and steep slope, and MC
is computed accordingly using maximum curvature method (Wiemer
and Wyss, 2000) and with bootstrapping (Schorlemmer et al., 2003).
The value is estimated to be 3.9 (Fig. 3). After considering MW?MC,
total 895 events have been considered as final database and are used to
calculate the b-values in the area. Epicenters of 895 earthquakes (Mw ?
3.9) are shown in Fig. 4. Epicenters of some strong earthquakes
(M?6.0) are also shown in Fig. 1. These events mostly show thrust/
reverse faulting; a few solutions with strike slip components (Kundu
and Gahalaut, 2013). The Gutenberg-Richter relation (Gutenberg and Richter, 1944) is
one of the well known empirical relations in seismology. It represents
the frequency of occurrence of earthquakes as a function of magnitude
as: log10N (M)=a?bM, where N (M) is the cumulative number of
earthquakes having magnitude equal to or larger than magnitude M. ‘a’
is a measure of seismic activity that depends on the size of the area and
constant ‘b’ describes the relative size distribution of earthquakes. Numerous
methods have been proposed to calculate out b-value and its’confidence limit (e.g. Aki, 1965; Utsu, 1965, 1992; Bender, 1983; Shi
and Bolt, 1982; Frohlich and Davis, 1993; Okal and Kirby, 1995; Scholz,
1997; Mousavi, 2017). In this study, the modified maximum likelihood
method of Aki (1965) and Utsu (1965) is utilized which is as follows:
Here, M denotes the average magnitude of a group of earthquakes
with M > MC, MC is the magnitude of completeness and Mbin is the
binning width. For spatial distribution, the study area is divided into
square grids having dimension of 1°×1° (geographical window) each
with an overlapping window of dimension 0.5°×0.5° (moving
window). To maintain the inherent grid to grid continuity of data
points, the moving window has been opted. Due to insufficient data, bvalue
calculations are not performed in the grids having Nmin (minimum
number of events) 120 km), the number of events are decreasing. This
may be due to low tectonic stress level at deeper depths, resulting the
increase in b-values. The b-values obtained in the present study are
compatible with the earlier reported b-value for the NE India (Khan
et al., 2011). In some segments of the study region reflect higher stress
accumulation compared to others, which may indicate the possibility of
future large rupture at these locations. Since, the Indo-Myanmar subduction
zone is an active region with the ongoing subduction of Indian
plate beneath the Burmese plate, the present study will provide useful
information to estimate seismic hazard assessment in the study region.
However, comprehensive catalogue will be essentially important in
future for more detailed analysis in the region. In this study, spatial and depth distribution of b-values are estimated
in the Indo-Myanmar subduction zone of NE India using earthquake
catalog from 1996 to 2015. Spatial variation map showed that the bvalue
varies from 0.75 to 1.54 within the region. Intermediate variation
of b-values were observed around the epicenter of the 3rd January,
MW=6.7, Manipur earthquake. These b-value anomalies can be interpreted
as a high stress accumulated zone before the Manipur earthquake.
Low b-value anomalies are observed along Sagaing fault and
Kabaw fault region. Low b-values (0.94–1) are found in the upper crust,
which indicate homogeneity and high stress accumulation in the upper
crust. Also, this low b-value may be due to high seismic moment release
in the uppermost crust. As b-value is directly related to the local and
regional factors including tectonic stresses, seismic moment release,
and earthquake focal mechanisms, our study on it could provide valuable
information of precursory signal of seismic hazard of the Indo-
Myanmar subduction zone. The authors JB and DB acknowledges financial support from
University Grant Commission, New Delhi, India, vide sanctioned No. F.
No.-43-522/2014 (SR). The figures were made using Generic Mapping
Tools (GMT) (Wessel and Smith, 1998). The Editor-in-Chief Prof. Mei-
Fu Zhou, Associate Editor Prof. Dapeng Zhao and three anonymous
reviewers are gratefully acknowledged for their constructive comments
and suggestions which helped us in upgrading manuscript significantly
Aki, K., 1965. Maximum likelihood estimate of b in the formula logN = a-bM and its
confidence limits. Bull. Earthq. Res. Inst. Tokyo. Univ. 43, 237–239.
Allen, C., Amand, P., Richter, C., Nordquist, J., 1965. Relation between seismicity and
geological structure in the southern California region. Bull. Seismol. Soc. Am. 55,
Bachmann, C.E., Wiemer, S., Goertz-Allmann, B.P., Woessner, J., 2012. Influence of
porepressure on the event-size distribution of induced earthquakes. Geophys. Res.
Lett. 39, L09302.
Bender, B., 1983. Maximum likelihood estimation of b-values for magnitude grouped
data. Bull. Seismol. Soc. Am. 73, 831–851.
Bernard, P., Lyon-Caen, H., Briole, P., Deschamps, A., Boudin, F., Makropoulos, K.,
Papadimitriou, P., Lemeille, F., Patau, G., Billiris, H., Paradissis, D., Papazissi, K.,
Castarede, H., Charade, O., Nercessian, A., Avallone, A., Pacchiani, F., Zahradnik, J.,
Sacks, S., Linde, A., 2006. Seismicity, deformation and seismic hazard in the western
rift of Corinth: new insights from the Corinth Rift Laboratory (CRL). Tectonophysics
Barman, P., Ray, J.D., Kumar, A., Chowdhury, J.D., Mahanta, K., 2014. Estimation of
present day inter-seismic deformation in Kopili fault zone of north-east India using
GPS measurements. Geomat. Nat. Hazards Risk 7, 586–599.
Bhattacharya, P.M., Majumder, R.K., Kayal, J.R., 2002. Fractal dimension and b-value
mapping in Northeast India. Curr. Sci. 82, 1486–1491.
Bhattacharya, P.M., Kayal, J.R., 2003. Mapping the b-values and its correlation with the
fractal dimension in the northeast region of India. J. Geol. Soc. India. 62, 680–695.
Bilham, R., England, P., 2001. Plateau Pop-up during the 1897 Assam earthquake. Nature
Bridges, D.L., Gao, S.S., 2006. Spatial variation of seismic b-values beneath Makushin
volcano Unalaska Island, Alaska. Earth Planet. Sci. Lett. 245, 408–415.
Bora, D.K., Baruah, S., Biswas, R., Gogoi, N.K., 2013. Estimation of source parameters of
local earthquakes originated in Shillong-Mikir plateau and its adjoining region of
Northeastern India. Bull. Seismol. Soc. Am. 103, 437–446.
Das, R., Wason, H.R., Sharma, 2012. Magnitude conversion to unified moment magnitude
using orthogonal regression relation. J. Asian Earth Sci. 50, 44–51.
Farrell, J., Husen, S., Smith, R.B., 2009. Earthquake swarm and b-value characterization
of the Yellowstone volcano-tectonic system. J. Volcanol. Geotherm. Res. 188,
Frohlich, C., Davis, S., 1993. Teleseismic b-values: or much ado about 1.0. J. Geophys.
Res. 98, 631–644.
Gahalaut, V.K., Stacey, S.M., Srinagesh, D., Kapil, S.L., Suresh, G., Saurav, S., Kumar, V.,
Dadhich, H., Patel, A., Prajapati, S.K., Shukla, H.P., Gautam, J.L., Baidya, P.R.,
Mandal, S., Jain, A., 2016. Seismological, geodetic, macroseismic and historical
context of the 2016 Mw 6.7 Tamenglong (Manipur) India earthquake. Tectonophysics
Gutenberg, B., Richter, C.F., 1944. Frequency of earthquake in California. Bull. Seismol.
Soc. Am. 34, 185–188.
Hurukawa, N., Tun, P.P., Bunichiro, S., 2012. Detailed geometry of the subducting Indian
Plate beneath the Burma Plate and subcrustal seismicity in the Burma Plate derived
from joint hypocenter relocation. Earth Planets Space 64, 333–343.
Hurukawa, N., Maung, P.M., 2011. Two seismic gaps on the Sagaing Fault, Myanmar,
derived from relocation of historical earthquakes since 1918. Geophys. Res. Lett. 38,
Jade, S., Mukul, M., Bhattacharya, A., Vijayan, M.S.M., Jaganathan, S., Kumar, A., Tiwari,
R.P., Kumar, A., Kalita, S., Sahu, S.C., 2007. Estimates of interseismic deformation in
northeast India from GPS measurements. Earth Plan. Sci. Lett. 263, 221–234.
Jolly, A.D., McNutt, S.R., 1999. Seismicity at the volcanoes of katmai national park,
Alaska, July 1995–December 1997. J. Volcanol. Geother. Res. 93, 173–190.
Kamer, Y., Hiemer, S., 2013. Comment on analysis of the b-values before and after the 23
October 2011 mw 7.2 Van-Ercis, Turkey, earthquake. Tectonophysics 608,
Kamer, Y., 2014. Comment on systematic survey of high-resolution b value imaging along
Californian faults: inference on asperities by T. Tormann et al. J. Geophys. Res. Solid
Earth 119, 2029–2054. http://dx.doi.org/10.1002/2014JB011147.
Kamer, Y., Hiemer, S., 2015. Data-driven spatial b-value estimation with applications to
California seismicity: to b or not to b. J. Geophys. Res. 120, 5191–5214.
Kayal, J.R., 2008. Microearthquake Seismology and Seismotectonics of South Asia.
Springer, pp. 273–275.
Kayal, J.R., Arefiev, S.S., Baruah, S., Hazarika, D., Gogoi, N., Gautam, J.L., Baruah, S.,
Dorbath, C., Tatevossian, R., 2012. Large and great earthquakes in the Shillong
plateau–Assam valley area of North east India region: pop-up and transverse tectonics.
Tectonophysics 532–535, 186–192.
Kayal, J.R., 1996. Earthquake source process in north east India: a review. J. Himala.
Geol. 17, 53–69.
Kayal, J.R., Arefiev, S.S., Barua, S., Hazarika, D., Gogoi, N., Kumar, A., Chowdhury, S.N.,
Kalita, S., 2006. Shillong plateau earthquakes in northeast India region: complex
tectonic model. Curr. Sci. 91, 109–114.
Khan, P.K., Ghosh, M., Pritam, P.C., Mukhargee, D., 2011. Seismic b-value and the assessment
of ambient stress in Northeast India. Pure Appl. Geophys. 168, 1693–1706.
Kumar, M.R., Rao, N.P., 1995. Significant trends related to the slab seismicity and tectonics
in Burmese arc region from Harvard CMT solutions. Phys. Earth Planet. Inter.
Kundu, B., Gahalaut, V.K., 2012. Earthquake occurrence processes in the Indo-Burmese
wedge and Sagaing fault region. Tectonophysics 524–525, 135–146.
Kundu, B., Gahalaut, V.K., 2013. Tectonic geodesy revealing geodynamic complexity of
the Indo-Burmese arc region, north east India. Curr. Sci. 104, 920–933.
Le Dain, A.Y., Tapponnier, P., Molnar, P., 1984. Active faulting and tectonics of Burma
and surrounding regions. J. Geophys. Res. 89, 453–472.
Lockner, D.A., Byrlee, J.D., 1991. Precursory AE patterns leading to rock fracture. In:
Hardy, H.R. (Ed.), Proc. 5th Conference on Acoustic Emisson/Microseismic Activity
in Geologic Structures and Materials, Trans-Tech.Publications, Clausthal-Zellerfeld,
Germany, pp. 1–14.
McNutt, S.R., 2005. Volcanic seismology. Ann. Rev. Earth Planet. Sci. 33, 461–491.
Meredith, P.G., Atkinson, B.K., 1983. Stress corrosion and acoustic emission during tensile
crack propagation in Whin Sill dolerite and other basic rocks. Geophys. J. R.
Astron. Soc. 75, 1–21.
Mogi, K., 1967. Regional variation in magnitude-frequency relation of earthquake. Bull.
Earthq. Res. Inst. 45, 313–325.
Mogi, K., 1962. Magnitude–frequency relation for elastic shocks accompanying fractures
of various materials and some related problems in earthquakes. Bull. Earthq. Res.
Inst. Univ. Tokyo 40, 831–853.
Monterroso, D.A., Kulhanek, O., 2003. Spatial variations of b-values in the subduction
zone of Central America. Geofís. Int. 42, 575–587.
Mousavi, S.M., 2017. Mapping seismic moment and b-value within the continental-collision
orogenic-belt region of the Iranian Plateau. J. Geodynamics 103, 26–41.
Murru, M., Montuori, C., Wyss, M., Privitera, E., 1999. The location of magma chambers
at Mt. Etna, Italy, mapped by b-values. Geophys. Res. Lett. 26, 2553–2556.
Nuannin, P., Kulhanek, O., Persson, L., 2005. Spatial and temporal b value anomalies
preceding the devastating off coast of NW Sumatra earthquake of December 26, 2004.
Geophys. Res. Lett. 32, L11307.
Okal, E.A., Kirby, S.H., 1995. Frequency-moment distribution of deep earthquake; implications
for the seismogenic zone at the bottom of slabs. Phys. Earth Planet Int. 92,
Oldham, R.D., 1899. Report on the great earthquake of the 12th June 1897. Geological
Survey of India Memoir 29, 1–379.
Power, J.A., Jolly, A.D., Page, R.A., McNutt, S.R., 1995. Seismicity and forecasting of the
1992 eruptions of crater peak vent, Mt. spurr, alaska: an overview, U. S. Geol. Surv.
Rastogi, B.K., 2016. A report on Mw 6.7 Manipur earthquake of January 4, 2016. J.
Earthquake Sci. Eng. 3, 29–35.
Ruff, L., Kanamori, H., 1980. Seismicity and the subduction process. Phys. Earth Planet
Rydelek, P.A., Sacks, I.S., 1989. Testing the completeness of earthquake catalogs and the
hypothesis of self-similarity. Nature 337, 251–253.
Sammonds, P.R., Meredith, P.G., Main, I.G., 1992. Role of pore fluids in the generation of
seismic precursors to shear fracture. Nature 359, 228–230.
Sanchez, J.J., McNutt, S.R., Power, J.A., Wyss, M., 2004. Spatial variations in the frequency-
magnitude distribution of earthquakes at Mount Pinatubo Volcano. Bull.
Seismol. Soc. Am. 94, 430–438.
Schorlemmer, D., Neri, G., Wiemer, S., Mostaccio, A., 2003. Stability and significance
tests for b-value anomalies: example from the Tyrrhenian Sea. Geophys. Res. Lett. 30,
Schorlemmer, D., Weimer, S., Wyss, M., 2004. Earthquake statistics at Parkfield:
1.Stationarity of b values. J. Geophys. Res. 109, B12307. http://dx.doi.org/10.1029/
Schorlemmer, D., Wiemer, S., Wyss, M., 2005. Variations in earthquake size distribution
across different stress regimes. Nature 437, 539–542.
Scholz, C.H., 1968. The frequency-magnitude relation of micro-fracturing in rock and its
relation to earthquakes. Bull. Seismol. Soc. Am. 58, 399–415.
Scholz, C.H., 1997. Size distributions for large and small earthquakes. Bull. Seismol. Soc.
Am. 87, 1074–1077.
Shi, Y., Bolt, B.A., 1982. The standard error of the magnitude frequency b value. Bull.
Seismol. Soc. Am. 72, 1677–1687.
Singh, C., Singh, A., Chadha, R.K., 2009. Fractal and b-value mapping in eastern
Himalaya and southern Tibet. Bull. Seismol. Soc. Am. 99, 3529–3533.
Singh, C., Bhattacharya, P.M., Chadha, R.K., 2008. Seismicity in Koyna-Warna Reservoir
Site in Western India: fractal and b-value mapping. Bull. Seismol. Soc. Am. 98,
Sorbi, M.R., Nilfouroushan, F., Zamani, A., 2012. Seismicity patterns associated with the
September 10th, 2008 Qeshm earthquake, South Iran. Int. J. Earth Sci. 101,
Tandon, A.N., 1954. A study of Assam earthquake of August, 1950 and its aftershocks.
Indian J. Metero. Geophys. 5, 95–137.
Utsu, T., 1965. A method for determining the value of b in a formula log n = a – bM
showing the magnitude frequency for earthquakes. Geophys. Bull. Hokkaido Univ.
Utsu, T., 1992. On seismicity. In: Report of the Joint Research Institute for Statistical
Mathematics. Institute for Statistical Mathematics, Tokyo, pp. 139–157.
Van Stiphout, T., Kissling, E., Wiemer, S., Ruppert, N., 2009. Magmatic processes in the
Alaska subduction zone by combined 3-D b value imaging and targeted seismic tomography.
J. Geophys. Res. 114, B11302.
Warren, N.W., Latham, G.V., 1970. An experiment study of thermal induced microfacturing
and its relation to volcanic seismicity. J. Geophys. Res. 75, 4455–4464.
Wessel, P., Smith, W.H.F., 1998. New, improved version of Generic Mapping Tools released.
EOS Trans. 79, 579.
White, B.J.P., Smith, R.B., Husen, S., Farrell, J., Wong, I., 2009. Seismicity and earthquake
hazard analysis of the Teton-Yellowstone region, Wyoming. J. Volcanol.
Geotherm. Res. 188, 277–296.
Wiemer, S., McNutt, S., 1997. Variations in frequency-magnitude distribution with depth
in two volcanic areas: Mount St. Helens, Washington, and Mt. Spurr, Alaska.
Geophys. Res. Lett. 24, 189–192.
Wiemer, S., McNutt, S.R., Wyss, M., 1998. Temporal and three-dimensional spatialanalyses
of the frequency-magnitude distribution near Long Valley Caldera, California.
Geophys. J. Int. 134, 409–421. http://dx.doi.org/10.1046/j.1365-246x.1998.
Wiemer, S., Katsumata, K., 1999. Spatial variability of seismicity parameters in aftershock
zones. J. Geophys. Res. 104, 13135–13151.
Wiemer, S., Wyss, M., 2000. Minimum magnitude of completeness in earthquake catalogs:
examples from Alaska, the western U.S. and Japan. Bull. Seismol. Soc. Am. 90,
Wiemer, S., 2001. A software package to analyze seismicity: ZMAP. Seismol. Res. Lett. 72,
Wiemer, S., Benoit, J., 1996. Mapping the b-value anomaly at 100 km depth in the Alaska
and New Zealand subduction zones. Geophys. Res. Lett. 23, 1557–1560.
Woessner, J., Wiemer, S., 2005. Assessing the quality of earthquake catalogues:estimating
the magnitude of completeness and its uncertainty. Bull. Seismol. Soc. Am. 95,
Wyss, M., 1973. Towards a physical understanding of the earthquake frequency distribution.
Geophys. J. R. Astron. Soc. 31, 341–359.
Wyss, M., Shimazaki, K., Wiemer, S., 1997. Mapping active magma chambers by b-value
beneath the off-Ito volcano. J. Geophys. Res. 102, 20,413–20,422.
Wyss, M., Hasegawa, J., Nakajima, J., 2001. Source and path of magma for volcanoes in
the subduction zone of northeastern Japan. Geophys. Res. Lett. 28, 1819–1822