Relationship between Mid-latitude Circulation Indices and Indian Northeast Monsoon Rainfall

—The study presents the results of the statistical relationship between seasonal northeast monsoon rainfall over Tamil Nadu state of India (TNR) and southeast India (SER) and mid-latitude circulation indices viz., zonal index (ZON) meridional index (MER) and the ratio of meridional to zonal index (M/Z) between the geographical area 35°N to 70°N at 500 hPa level over three sectors and hemisphere, based on 19 years (1971–1989) of data. The results indicate that northeast monsoon rainfall over India shows a strong antecedent relationship with the strength of ZON over all the sectors and hemisphere. The best association is observed during antecedent March over sector I (45°W–90°E) where direct and strong correlation coefficients of 0.69 and 0.64 are obtained with TNR and SER, respectively. Antecedent MAM (spring) season over sector I also shows a significant positive correlation with TNR/SER. Thus, the mid-latitude zonal circulation index may have possible use for the long-range forecasting of northeast monsoon rainfall over India.     

Numerical Simulation of the Sensitivity of Summer Monsoon Circulation and Rainfall over India to Land Surface Processes

—The influence of soil moisture and vegetation variation on simulation of monsoon circulation and rainfall is investigated. For this purpose a simple land surface parameterization scheme is incorporated in a three-dimensional regional high resolution nested grid atmospheric model. Based on the land surface parameterization scheme, latent heat and sensible heat fluxes are explicitly estimated over the entire domain of the model. Two sensitivity studies are conducted; one with bare dry soil conditions (no latent heat flux from land surface) and the other with realistic representation of the land surface parameters such as soil moisture, vegetation cover and landuse patterns in the numerical simulation. The sensitivity of main monsoon features such as Somali jet, monsoon trough and tropical easterly jet to land surface processes are discussed.¶Results suggest the necessity of including a detailed land surface parameterization in the realistic short-range weather numerical predictions. An enhanced short-range prediction of hydrological cycle including precipitation was produced by the model, with land surface processes parameterized. This parameterization appears to simulate all the main circulation features associated with the summer monsoon in a realistic manner.     

Variability of Convective Activity over the North Indian Ocean and its Associations with Monsoon Rainfall over India

The present study is an attempt to examine the variability of convective activity over the north Indian Ocean (Bay of Bengal and Arabian Sea) on interannual and longer time scale and its association with the rainfall activity over the four different homogeneous regions of India (viz., northeast India, northwest India, central India and south peninsular India) during the monsoon season from June to September (JJAS) for the 26 year period (1979 to 2004). The monthly mean Outgoing Long-wave Radiation (OLR) data obtained from National Oceanic and Atmospheric Administration (NOAA) polar orbiting spacecraft are used in this study and the 26-year period has been divided into two periods of 13 years each with period-i from 1979 to 1991 and period -ii from 1992 to 2004. It is ascertained that the convective activity increases over the Arabian Sea and the Bay of Bengal in the recent period (period -ii; 1992 to 2004) compared to that of the former period (period -i; 1979 to 1991) during JJAS and is associated with a significantly increasing trend (at 95% level) of convective activity over the north Bay of Bengal (NBAY). On a monthly scale, July and August also show increase in convective activity over the Arabian Sea and the Bay of Bengal during the recent period and this is associated with slight changes in the monsoon activity cycle over India. The increase in convective activity particularly over the Arabian Sea during the recent period of June is basically associated with about three days early onset of the monsoon over Delhi and relatively faster progress of the monsoon northward from the southern tip of India. Over the homogeneous regions of India the correlation coefficient (CC) of OLR anomalies over the south Arabian Sea (SARA) is highly significant with the rainfall over central India, south peninsular India and northwest India, and for the north Arabian Sea (NARA), it is significant with northwest India rainfall and south peninsular rainfall. Similarly, the OLR anomalies over the south Bay of Bengal (SBAY) have significant CC with northwest India and south peninsular rainfall, whereas the most active convective region of the NBAY is not significantly correlated with rainfall over India. It is also found that the region over northeastern parts of India and its surroundings has a negative correlation with the OLR anomalies over the NARA and is associated with an anomalous sinking (rising) motion over the northeastern parts of India during the years of increase (decrease) of convective activity over the NARA.     

Observed and Forecasted Intraseasonal Activity of Southwest Monsoon Rainfall over India During 2010, 2011 and 2012

The monsoon seasons of 2010 and 2011, with almost identical seasonal total rainfall over India from June to September, are associated with slightly different patterns of intraseasonal rainfall fluctuations. Similarly, the year 2012, with relatively less rainfall compared to 2010 and 2011, also witnessed different intraseasonal rainfall fluctuations, leading to drought-like situations over some parts of the country. The present article discusses the forecasting aspect of monsoon activity over India during these 3 years on an extended range time scale (up to 3 weeks) by using the multimodel ensemble (MME), based on operational coupled model outputs from the ECMWF monthly forecasting system and the NCEP’s Climate Forecast System (CFS). The average correlation coefficient (CC) of weekly observed all-India rainfall (AIR) and the corresponding MME forecast AIR is found to be significant, above the 98 % level up to 2 weeks (up to 18 days) with a slight positive CC for the week 3 (days 19–25) forecast. However, like the variation of observed intraseasonal rainfall fluctuations during 2010, 2011 and 2012 monsoon seasons, the MME forecast skills of weekly AIR are also found to be different from one another, with the 2012 monsoon season indicating significant CC (above 99 % level) up to week 2 (12–18 days), and also a comparatively higher CC (0.45) during the week 3 forecast (days 19–25). The average CC between observed and forecasted weekly AIR rainfall over four homogeneous regions of India is found to be the lowest over the southern peninsula of India (SPI), and northeast India (NEI) is found to be significant only for the week 1 (days 5–11) forecast. However, the CC is found to be significant over northwest India (NWI) and central India (CEI), at least above the 90 % level up to 18 days, with NWI having slightly better skill compared to the CEI. For the individual monsoon seasons of 2010, 2011 and 2012, there is some variation in CC and other skill scores over the four homogeneous regions. Thus the slight variations in the characteristics of intraseasonal monsoon rainfall over India is associated with variations in predictive skill of the coupled models and the MME-based predictions of intraseasonal monsoon fluctuations for 2–3 weeks, providing encouraging results. The MME forecast in 2010 is also able to provide useful guidance, well in advance, about an active September associated with a delayed withdrawal of the monsoon and also the heavy rainfall over north Pakistan.     

Crustal structure of the Peninsular India

Summary Earthquake parameters for the forty aftershocks of the main Koyna earthquake of 10 December, 1967, have been determined. Depths of the foci of the earthquakes have been found to vary between 2 to 17 km. The velocities for the phasesP _g ,P ^*,P _n have been observed to be 5.78±0.00, 6.58±0.04, 8.19±0.02 km/sec, and forS _g ,S ^*,S _n to be 3.42±0.00, 3.92±0.01 and 4.62±0.01 km/sec respectively. A two-layered crustal model has been interpreted for the Peninsular shield with the average thickness of the granitic layer as 20 km and that for the basaltic layer as 18.7 km. A plot of the epicenters suggests a NNE to SSW orientation of the fault.     

Predictability of Indian Monsoon Circulation with High Resolution ECMWF Model in the Perspective of Tropical Forecast During the Tropical Convection Year 2008

To address some of the issues of project Year of Tropical Convection (YOTC) and the project ATHENA as ongoing international activities, an endeavor has been made for the first time to study the predictability of Indian summer monsoon in the backdrop of tropical predictability using 850 hPa atmospheric circulations with the high resolution (T1279) ECMWF model during the boreal summer of 2008 as one of the focus years of YOTC. The major findings obtained from the statistical forecast have been substantiated by the dynamical prediction in terms of the systematic error energy, its growth rate and the attribution of the dominant nonlinear dynamical processes to error growth. The systematic error energy of T1279 (16 km resolution) ECMWF model are generated in African landmass, India and its adjoining oceanic region, in near equatorial west Pacific and around the Madagascar region where the root mean square errors are observed and the zonal wind anomaly shows poor forecast skill. As far as the inadequate predictability of Indian summer monsoon by T1279 ECMWF model (revealed from the results of project ATHENA) is concerned, the systematic error energy and the error growth over Arabian Sea, in the eastern and western India due to the nonlinear convergence and divergence of error flux along with the erroneous Mascarene high may possibly be the determining factors for not showing any discernable improvement in Indian monsoon during the medium range forecast up to 240 h. This work suggests that the higher resolution of ECMWF model may not necessarily lead to the better forecast of Indian monsoon circulations during 2008 unless a methodology can be devised to isolate the errors due to the nonlinear processes that are inherent within the system.     

Influence of the Planetary Boundary Layer Physics on Medium-range Prediction of Monsoon over India

—The present study emphasizes the importance of proper representation of boundary layer physics in a general circulation model. The Turbulent Kinetic Energy (TKE) closure scheme incorpo rates important processes of the Planetary Boundary Layer (PBL) compared to a simplistic first-order closure model. Hence the model which has the TKE closure scheme is capable of simulating important weather systems associated with summer monsoon, such as monsoon depressions and lows that form over the Indian subcontinent quite well compared to the first-order closure model. The present study indicates better performance of the global model with the TKE scheme in the prediction of the monsoon circulation, including the tracks of the depressions over the Indian subcontinent. Medium-range weather prediction has also improved with the use of the TKE closure. However further studies are necessary to improve the forecast, with emphasis on boundary layer processes.     

An Objective Approach for Prediction of Daily Summer Monsoon Rainfall over Orissa (India) due to Interaction of Mesoscale and Large-scale Synoptic Systems

Orissa State, a meteorological subdivision of India, lies on the east coast of India close to north Bay of Bengal and to the south of the normal position of the monsoon trough. The monsoon disturbances such as depressions and cyclonic storms mostly develop to the north of 15° N over the Bay of Bengal and move along the monsoon trough. As Orissa lies in the southwest sector of such disturbances, it experiences very heavy rainfall due to the interaction of these systems with mesoscale convection sometimes leading to flood. The orography due to the Eastern Ghat and other hill peaks in Orissa and environs play a significant role in this interaction. The objective of this study is to develop an objective statistical model to predict the occurrence and quantity of precipitation during the next 24 hours over specific locations of Orissa, due to monsoon disturbances over north Bay and adjoining west central Bay of Bengal based on observations to up 0300 UTC of the day. A probability of precipitation (PoP) model has been developed by applying forward stepwise regression with available surface and upper air meteorological parameters observed in and around Orissa in association with monsoon disturbances during the summer monsoon season (June-September). The PoP forecast has been converted into the deterministic occurrence/non-occurrence of precipitation forecast using the critical value of PoP. The parameters selected through stepwise regression have been considered to develop quantitative precipitation forecast (QPF) model using multiple discriminant analysis (MDA) for categorical prediction of precipitation in different ranges such as 0.1–10, 11–25, 26–50, 51–100 and >100 mm if the occurrence of precipitation is predicted by PoP model. All the above models have been developed based on data of summer monsoon seasons of 1980–1994, and data during 1995–1998 have been used for testing the skill of the models. Considering six representative stations for six homogeneous regions in Orissa, the PoP model performs very well with percentages of correct forecast for occurrence/non-occurrence of precipitation being about 96% and 88%, respectively for developmental and independent data. The skill of the QPF model, though relatively less, is reasonable for lower ranges of precipitation. The skill of the model is limited for higher ranges of precipitation. accepted September 2006     

Representing the Sensitivity of Convective Cloud Systems to Tropospheric Humidity in General Circulation Models

Convective cloud variability on many times scales can be viewed as having three major components: a suppressed phase of shallow and congestus clouds, a disturbed phase of deep convective clouds, and a mature phase of transition to stratiform upper-level clouds. Cumulus parameterization development has focused primarily on the second phase until recently. Consequently, many parameterizations are not sufficiently sensitive to variations in tropospheric humidity. This shortcoming may affect global climate model simulations of climate sensitivity to external forcings, the continental diurnal cycle of clouds and precipitation, and intraseasonal precipitation variability. The lack of sensitivity can be traced in part to underestimated entrainment of environmental air into rising convective clouds and insufficient evaporation of rain into the environment. As a result, the parameterizations produce deep convection too easily while stabilizing the environment too quickly to allow the effects of convective mesoscale organization to occur. Recent versions of some models have increased their sensitivity to tropospheric humidity and improved some aspects of their variability, but a parameterization of mesoscale organization is still absent from most models. Evidence about the effect of these uncertainties on climate change projections suggests that climate modelers should make improved simulation of high and convective clouds as high a priority as better representations of low clouds.     

Role of Ocean in the Variability of Indian Summer Monsoon Rainfall

Asian summer monsoon sets in over India after the Intertropical Convergence Zone moves across the equator to the northern hemisphere over the Indian Ocean. Sea surface temperature (SST) anomalies on either side of the equator in Indian and Pacific oceans are found related to the date of monsoon onset over Kerala (India). Droughts in the June to September monsoon rainfall of India are followed by warm SST anomalies over tropical Indian Ocean and cold SST anomalies over west Pacific Ocean. These anomalies persist till the following monsoon which gives normal or excess rainfall (tropospheric biennial oscillation). Thus, we do not get in India many successive drought years as in sub-Saharan Africa, thanks to the ocean. Monsoon rainfall of India has a decadal variability in the form of 30-year epochs of frequent (infrequent) drought monsoons occurring alternately. Decadal oscillations of monsoon rainfall and the well-known decadal oscillation in SST of the Atlantic Ocean (also of the Pacific Ocean) are found to run parallel with about the same period close to 60 years and the same phase. In the active–break cycle of the Asian summer monsoon, the ocean and the atmosphere are found to interact on the time scale of 30–60 days. Net heat flux at the ocean surface, monsoon low-level jetstream (LLJ) and the seasonally persisting shallow mixed layer of the ocean north of the LLJ axis play important roles in this interaction. In an El Niño year, the LLJ extends eastwards up to the date line creating an area of shallow ocean mixed layer there, which is hypothesised to lengthen the active–break (AB) cycle typically from 1 month in a La Niña to 2 months in an El Niño year. Indian monsoon droughts are known to be associated with El Niños, and long break monsoon spells are found to be a major cause of monsoon droughts. In the global warming scenario, the observed rapid warming of the equatorial Indian ocean SST has caused the weakening of both the monsoon Hadley circulation and the monsoon LLJ which has been related to the observed rapid decreasing trend in the seasonal number of monsoon depressions.     

A Note on the Maintenance of Asian Summer Monsoon Circulation

In this paper, the maintenance of Asian summer monsoon circulation is compared in the National Centres for Environmental Prediction/National Centre for Atmospheric Research (NCEP/NCAR) reanalysis and National Centre for Medium Range Weather Forecasting (NCMRWF), India operational analysis. The time mean summer monsoon circulation is bifurcated into stable mean and transient eddy components. The mean component of the summer monsoon circulation is elucidated through the study of large-scale balances of kinetic energy, heat and moisture. Although the balance requirements are fairly satisfied by both NCEP/NCAR and NCMRWF fields, a major difference is noticed. Strong monsoon circulation is evinced by NCEP/NCAR over the Oceanic regions and NCMRWF over the landmass. The various mean budgets under consideration indicated this feature invariably.     

The Assimilation of GPS Radio Occultation Data and its Impact on Rainfall Prediction along the West Coast of India during Monsoon 2002

In this study, the Weather Research and Forecasting (WRF-2.0.3.1) model with three-dimensional variational data assimilation (3DVAR) was utilized to study a heavy rainfall event along the west coast of India with and without the assimilation of GPS occultation refractivity soundings in the monsoon period of 2002. The WRF model is a next-generation mesoscale numerical weather prediction system designed to serve both operational forecasting and atmospheric research communities. The Global Positioning System (GPS) radio occultation (RO) refractivity data, processed by UCAR, were obtained from the CHAMP and SAC-C missions. This study investigates the impact of thirteen GPS occultation refractivity soundings only, as assimilated into the WRF model with 3DVAR, on the rainfall prediction over the western coastal mountain of India. The model simulation, with the finest resolution of 10 km, was in good agreement with rainfall observations, up to 72-h forecast. There are some subtle but important differences in predicted rainfalls between the control run CN (without the assimilation of refractivity soundings) and G13 (with the assimilation of thirteen GPS RO soundings). In general, the assimilation run G13 gives a better prediction in terms of both rainfall locations and amounts at later times. The moisture increments were analyzed at the initial and forecast times to assess the impact of GPS RO data assimilation. The results indicate that remote soundings in the forcing region could have significant impacts on distant downstream regions. It is anticipated, based on this study, that considerably occultation soundings available from the six-satellite constellation of FORMOSAT-3/COSMIC would have even more significant impacts on weather prediction in this region.     

Observation of Sea Breeze Front and its Induced Convection over Chennai in Southern Peninsular India Using Doppler Weather Radar

Sea breeze, the onshore wind over a coastal belt during daytime, is a welcoming weather phenomenon as it modulates the weather condition by moderating the scorching temperature and acts as a favourable mechanism to trigger convection and induce precipitation over coastal and interior locations. Sea breeze aids dispersal of pollutants as well. Observational studies about its onset, depth of circulation and induced precipitation have been carried out in this paper for the period April to September, 2004–2005 using a S-band Doppler Weather Radar functioning at Cyclone Detection Radar Station, India Meteorological Department, Chennai, India. The onset of sea breeze has been observed to be between 0900 and 1000 UTC with the earliest onset at 0508 UTC and late onset at 1138 UTC. The frequency is greater during the southwest monsoon season, viz., June – September and the frequency of initial onset is greater in north Chennai. The modal length of sea breeze is between 20 and 50 km with extreme length as high as 100 km also having been observed. Though the inland penetration is on average 10 to 20 km, penetration reaching 100 km was also observed on a number of cases. The induced convection could be seen in the range 50–100 km in more than 53% of the cases. The mean depth of sea breeze circulation is 300–600 m but may go well beyond 1000 m on conducive atmospheric conditions.     

The Role of Radiative Transfer in Maintaining the Indian Summer Monsoon Circulation

—The radiative-convective feedback and land-sea thermal forcing play significant roles in maintenance of the summer monsoon circulation over the Indian sub-continent. In this study, the role of radiative transfer in maintaining the monsoon circulation is examined with numerical sensitivity experiments. For this purpose, a sixteen layer primitive equation limited area model is used to perform numerical simulations with and without atmospheric radiative transfer processes parameterized in the model. The initial values and boundary conditions for the numerical integrations of the model are derived from operational analyses of the ECMWF, UK. The results show that the radiative transfer is essential in maintaining the intensity of the low level Somali Jet as well as the upper level Tropical Easterly Jet (TEJ) over the Indian sub-continent and adjoining seas. The meridional circulation over the region is also well simulated. As a result, enough moisture transports from the warm equatorial region to simulate more realistic orographic precipitation in the windward side of the mountains along the West coast of India. Without radiative transfer processes in the model atmosphere the simulated monsoon circulation weakens, moisture transport decreases and the precipitation lessens.     

Application of the time-predictable model in Peninsular India for future seismic hazard assessment

It has been the belief among Earth scientists that the Peninsular Shield is aseismic, as the region attained stability long ago. However, the earthquake at Koyna (10 December 1967), Bhadrachalam (13 April 1969), Broach (23 March 1970), Hyderabad (30 June 1983), Khillari (30 September 1993), Jabalpur (22 May 1997), Gujarat (26 January 2001), and additional ones of smaller magnitudes, altered this concept. This area has experienced many widely distributed shallow earthquakes, some of them having large magnitudes. It is now widely accepted that seismic activity still continues with moderate events. Therefore, a need has arisen to take into consideration recent seismological data to assess the future seismic status of Peninsular India. Earthquake generation model has been studied to develop the statistical relations with surface wave magnitude (M _S ≥ 4.5). Five seismogenic sources showing clustering of earthquakes and including at least three main shocks of magnitude 4.5 ≤ M _S ≤ 6.5 giving two repeat times, have been identified. It is mainly based on the so-called “regional time-predictable model”. For the considered region it is observed that the time interval between two consecutive main shocks depends on the preceding main shock magnitude (M _p ) not on the following main shocks magnitude M _f suggesting the validity of time predictable model in the region.     

Estimation of Maximum Earthquakes in Northeast India

We attempt to estimate possible maximum earthquakes in the northeast Indian region for four seismic source zones, namely EHZ, MBZ, EBZ, and SHZ, which encapsulates the various seismogenic structures of the region and also for combined source zones taken as a single seismic source regime. The latter case exhibits a high maximum earthquake estimate of M_W 9.4 (±0.85) through Bayesian interpretation of frequency magnitude distribution with Gamma function implicating a moderate deviation from the standard Gutenberg Richter model at the higher magnitudes. However, tapering Gutenberg Richter models with corner magnitudes at M_W 8.01, 8.7 and 9.1, respectively indicated maximum values corresponding to M_W 8.4, 9.0, and 9.3. The former approach was applied to each of the source zones wherein the data are presented in parts according to the data completeness, thereof. EHZ, MBZ, EBZ and SHZ are seen with maximum earthquakes of M_W 8.35 (±0.59), 8.79 (±0.31), 8.20 (±0.50), and 8.73 (±0.70), respectively. The maximum possible earthquakes estimated for each individual zone are seen to be lower than that estimated for the single regime. However, the pertaining return periods estimated for the combined zone are far less than those estimated for the demarcated ones.     

Micro pulse lidar observations of mineral dust layer in the lower troposphere over the southwest coast of Peninsular India during the Asian summer monsoon season

A large aerosol plume with optical depth exceeding 0.7 engulfs most parts of the Arabian Sea during the Asian summer monsoon season. Based on Micro Pulse Lidar observations during the June–September period of 2008 and 2009, the present study depicts, for the first time, the existence of an elevated dust layer occurring very frequently in the altitude band of 1–3.5 km over the west coast of peninsular India with relatively large values of linear depolarization ratio (δ_L). Large values of δ_L indicate the dominance of significantly non-spherical aerosols. The aerosol optical depth of this layer (0.2) is comparable to that of the entire atmospheric column during dust-free days. Back-trajectory analysis clearly shows the advection of airmass from the arid regions of Arabia and the west Arabian Sea, through the altitude region centered around 3 km. This is in contrast to the airmass below 1 km originating from the pristine Indian Ocean region which contains relatively spherical aerosols of marine origin with δ_L generally <0.05.     

On cumulus-scale transport of horizontal momentum in monsoon depression over India

A diagnostie method of cumulus parameterization is suggested in which vertical transport of horizontal momentum by cumulus-scale is derived by making use of large-scale vorticity as well as divergence budget equations. Data for composite monsoon depression over India available from our earlier studies used to test the method. As a first approximation, the results are obtained using only the vorticity budget equation.The results show that in the southwest sector of the monsoon depression, which is characterized by maximum cloudiness and precipitation, there is an excess of cyclonic vorticity in the lower troposphere and anticyclonic vorticity in the upper troposphere associated with the large-scale motion. The distribution of eddy vertical transport of horizontal momentum is such that anticyclonic vorticity is generated in the lower troposphere and cyclonic vorticity aloft. Cumulus-scale eddies thus work against the large-scale system and tend to off-set the large-scale imbalance in vorticity.     

Multiscale analysis of Asian Monsoon over the past 640 ka

The empirical mode decomposition(EMD) method is used to re-analyse the high-resolution and precisely-dated stalagmite record from Chinese caves over the past 640 ka. Results show that(1) the variation in the Asian Monsoon can be completely decomposed into ten quasiperiod oscillations, among which the precession and semiprecession band oscillations are the most prominent periodicities, with contribution rates of 31.1% and 30.7%, respectively;(2) the cross-spectrum analysis of the semiprecession component and bi-hemisphere insolation(BHI) are strongly correlated, indicating an amplified response of precipitation and temperature variability to the interhemispheric insolation in the low-latitude regions, thus further affecting the intensity of the Asian Monsoon;(3) on millennial timescales, obvious oscillations at the 5 ka and 1–2 ka bands roughly correspond to the classical Bond and Dansgaard-Oeschger(DO) cycles. Additionally, a strong correlation is found between the detrended stalagmite δ18 O records and Ca/Sr sequence from the North Atlantic(especially at the 5 ka band). This result means that the 5 ka cycle is characteristic of the glacial-interglacial cycle since the middle and late Pleistocene and may imply that climate change on the millennial timescale is the result of an interaction between global ice volume and insolation.     

Seismic source characterization for the Shillong Plateau in Northeast India

The present study aims at understanding the seismotectonic province of the Shillong Plateau (SP) by identifying the potential seismic source zones within a radius of 500 km from the centre of the SP. From existing literature and earthquake (EQ) data, the seismotectonic region is found to vary in terms of seismicity, tectonic features, geology, thickness of overburden, rupture characteristics and rate of movement. Thus, entire 500-km-radius seismotectonic region is divided into four seismic source zones: namely (1) the Shillong Plateau–Assam Valley Zone (SP-AVZ), (2) the Indo-Burma Ranges Zone (IBRZ), (3) the Bengal Basin Zone (BBZ) and (4) the Eastern Himalaya Zone (EHZ). EQ catalogues for each source zone is analysed for completeness of magnitude and time. Seismic parameter b estimated using a maximum likelihood method is found to be 0.91 ± 0.03, 0.94 ± 0.02, 0.80 ± 0.03 and 0.89 ± 0.03 for the SP-AVZ, IBRZ, BBZ and EHZ, respectively. In addition, the maximum likelihood method is used to estimate the mean annual activity rate, maximum possible magnitude (m _max), return period and probability of exceedance for the four zones. The b values estimated suggest that the BBZ is seismically more active; however, the rate of occurrence of EQs is highest in the IBRZ. Findings from this study are an indication of the relative contribution from each of the four seismic source zones towards a seismic hazard of the SP.