An Investigation of High An Annual Progress Report Resolution Spatio–Temporal (2016) Variation of Morphological, Submitted by Microphysical and Rainfall the Leading Institution: Properties of Precipitating Science Systems and its Social Impact: College, An Integrated Multi Sensor and Jotsoma, Multi Institutional Approach , 797002, INDIA

A study under the Research, Innovation & Quality Improvement (RI & QI) component of RUSA An investigation of high resolution spatio –temporal variation of morphological, microphysical and rainfall properties of precipitating systems and its social impact: An integrated multi sensor and multi institutional approach

Project Progress Report 1st Year (2016) Submitted to Ministry of Human Resource Development (Higher Education) Rashtriya Uchchatar Shiksha Abhiyan (RUSA)

By , Jotsoma, Kohima, Nagaland, 797002 (The Lead Institution)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 1

Content

Preface 04 General Information 06 Financial status for 1st Year 07 Outcome of the project for 1st year 08 Scientific Objectives of the Approved Project 09 Executive Summary of the Scientific Report 10 Scientific Progress Report-2016 14 Chapter 1: 15 1.1 Introduction 16 1.2 Data and Methodology 18 1.3 Results 19 1.3.1 Spatial and seasonal variation of rainfall contribution 19 1.3.2. Intra-seasonal variation of rainfall contribution 23 1.3.3 Spatial and seasonal variation of extreme rain intensity 27 1.3.4 Regional variability of the characteristics of convective 29 systems associated with extreme rain intensity 1.3.5 Spatial & seasonal variation of ice microphysical properties. 32 1.4 Summary and conclusion 35 Chapter 2: 37 2.1 Introduction 38 2.2 Result 38 2.2.1 Spatiotemporal variation of occurrence of MCSs during the 38 premonsoon and monsoon season (Diurnal variation) 2.2.2 Diurnal variation of occurrence of MCSs over two different 41 climatic regions, A and B during the premonsoon and monsoon seasons

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 2

2.2.3 Diurnal variation of the convective activity as observed by 43 INSAT 3D infrared thermal channel 2.3 Summary and Conclusions 45 Chapter 3 46 3.1. Introduction 47 3.2 Study Area 49 3.3 Data products 51 3.4 Methodology 52 3.5 Results 53 3.5.1 Spatial variation of hail events 53 3.5.2 Vertical profiles of radar reflectivity during hail events 56 3.5.3 Monthly and diurnal variation of hail events 58 3.5.4 Interannual variation of occurrence of detected and 60 reported hail events 3.5.5 Case study of the hail events 61 3.6 Summary and Conclusion 63 Future Plan 65 References 66 Annexures (I-IV) 71

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 3

Preface

The project entitled “An investigation of high resolution spatio–temporal variation of morphological, microphysical and rainfall properties of precipitating systems and its social impact: An integrated multi sensor and multi institutional approach” was approved in the 9th PAB meeting on 1st December 2015 with a total budget of Rs. 227.34 Lakhs. Kohima Science College an autonomous government post graduate college is a lead institute for the approved project. The main scientific objectives of the approved projects are (i) To study the spatial and temporal (seasonal and diurnal) variability of the morphological, microphysical and rainfall properties of the precipitating systems (ii) To study the extreme weather systems in terms of these properties and atmospheric energetic.(iii) To develop a comprehensive “Precipitation Features & Atmospheric Data Base” with the help of satellite and ground based observations and (iv). To study the social impact of weather related disasters and its socio –psychological response. A total amount of Rs. 113.67 Lakhs, including the state share has been received by the college. At the outset, an eight members expert committee was formed, with Prof. Animesh Maitra, Department of Electronics & Radio Physics, Calcutta University, Kolkata as a chairman. It was followed by appointment of three research scholars and one project assistant and nomination of nine associate members for the project. The formation of expert committee, appointment of project staff and nomination of associate members have government approval. The nominated associate members are faculty members of different government colleges and institutes within and outside Nagaland. The associate members are eligible for access to the data repository developed under the project. They are also eligible for sponsorship to participate in project related academic and research activities. During the last one year, a data repository consisting of satellite data over Indian subcontinent and adjoining ocean and ground reporting of thunderstorms, hail and lightning over north eastern region during 1998- 2014 is developed. The data collection work is in progress. The data set is analyzed and a detail scientific progress report is presented along with the annual progress

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 4 report. Significant spatiotemporal variation in the properties of convective systems is observed over the study region. Few scientific results are published in a reputed journal (International Journal of Climatology, a Journal of Royal Meteorological Society UK, published by Wiley International with a impact factor of 3.61). Work is in progress to prepare more manuscripts for submission to scientific journals. Five research papers are also presented in a regional and in a international conference. Project staff including the principal investigator have participated in various training programs. Five associate members of the project are sponsored to participate in a short term course entitled “ Landslides and debris flow systems: Prediction, control and reclamation” under the Global Initiatives of Academic Network (GIAN) at Kohima during March 7th -11th, 2017. Two main stations, one at Kohima Science college, Kohima, Nagaland and other at Indian Statistical Institute, Giridih (Jharkhand) and twelve substations at various colleges in Nagaland are identified for placing the instruments. The work on installation of instruments at various locations is in progress. A mutual agreement is signed with Prof. Alexander G Keul, Environmental Psychology, Psychology Department, University of Salzburg, Austria to study the social impact of weather related disasters and its socio –psychological response. A set of questionnaires is developed for the analysis of socio –psychological response. A survey is in progress to collect the data in the state of Nagaland and west Bengal During the last one year, a total of Rs. 79.87 Lakhs have been utilized and balance amount is under commitment. On behalf of all the members of the project I would like to take the opportunity to present a detail progress report of the project. Looking forward for an enriching experience in coming years

Dr. Sanjay Sharma Principal Investigator RUSA (Research, Innovation & Quality Improvement) Project Kohima Science College, Jotsoma, Kohima, Nagaland

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 5

General Information

1 Title of the Project An investigation of high resolution spatio – temporal variation of morphological, microphysical and rainfall properties of precipitating systems and its social impact: An integrated multi sensor and multi institutional approach 2 Name and address of the Kohima Science College, Jotsoma, Kohima, Institution Nagaland, 797002 (Lead Institution)

3 Principle investigator with Dr. Sanjay Sharma address Assistant Prof. ( Sr.), Department of Physics Kohima Science College, Jotsoma, Nagaland, 797002 4 Staff under the project 1. Partha Roy (Research Scholar) 2. Rupraj Bishwasharma (Research Scholar) 3. Imolemba (Research Scholar) 4. Vitsiavi Nyuthe (Project Assistant) 5 Nominated Associate 09 ( Please refer Annexure I) members

6 Nominated members of 8 members (Please refer Annexure II) expert committee

7 Duration of the project Three Years

8 Starting date of the project 01-01-2016

9 Total outlay of the Project Rs. 227.34 Lakhs

10 Amount Received Rs. 113.67 Lakhs

11 Amount Utilized Rs. 79.87 Lakhs

12 No of Publications. 01 (Annexure III A)

13. Paper presented in the 05 (Annexure III B) conferences.

14. Participation to short term 03 (Annexure III C) courses/training programs

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 6

Financial Status for 1st Year (2016)

Fund Expenditure Remark Received (Rs.) Head Item Amount (lakhs) (Rs.) (lakhs) `102.30 Equipment  Parsivel Disdrometer (02 19.00 Total 31.04 and units) balance of 12.04 (Vendor: Sutron Hydromet Systems lakhs to be utilized 11.37 Pvt. Limited, New Delhi) (Order placed on Sep th  Automatic rain gauges (10 26 , 2016 ) units)

(Vendor: Sutron Hydromet Systems Pvt. Limited, New Delhi

One unit of Disdrometer and 09 units of rain gauges are delivered at Kohima. The one unit each of Disdrometer and rain gauge is yet to be delivered at Indian Statistical Institute Giridih, Jharkhand. Lightning detector Sensors (04 8.60 Total 14.10 and Unit) balance of 5.5 lakhs (Vendor : Dhruva Technologies to be utilized Pvt. Ltd, New Delhi) (Order placed on Oct 20th, 2016) Micro Rain Radar 22.00 Utilized (Order Vendor: Electrotek International placed on Oct 20th, Inc, Chennai) 2016) Computers ( work station, Laptops) 9.00 Utilized , Printers, UPSs, Short Projectors

Salary Three Research scholars and one 9.28 Utilized project Assistant Contingency Travel, data download, Data 5.45 Utilized collection, Stationary Institutional To the Host Institute 6.54 Utilized Overhead 113.67 79.87 Balance amount to be utilized after installation of the instruments : Rs. 17.54 lakhs

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 7

Outcome of Project for the 1st Year

Sl No Progress Action Taken

1 Physical Progress  A government nominated expert committee for the review of the project is formed.  Three Research Scholars and one project Assistant are appointed.

 9 Associated members of the RUSA (RI &QI) project (faculty members from various Colleges/institute) are nominated.  After due process, order for the Instruments have been placed. Custom Duty Exemption (CDE) certificates for

the imported instruments are provided to the vendors.  Multi institutional site selection for the installation of the instruments is completed.

2. Scientific Progress  Satellite data are down loaded TRMM satellite : 1998-2015 INSAT 3D satellite : 2014-2016 GPM satellite : 2014-2016  Meteorological Data Hail and Lightning data from IMD* : 1998-2016 *IMD: India meteorological Department  Analysis of the satellite and other collected data is in progress ( 70 % of analysis is completed).  MOU is signed with the Prof. Alexander Keul Department of Environmental Psychology, Salzburg University, Salzburg, Austria. It resulted in the development of the questionnaires for the survey to asses the psychological response of the people to the severe weather ( Please refer Annexure IV for the questionnaires) Survey has already started in Nagaland and West Bengal. In coming days survey will also start in the state of Assam ( 40% of the sample collection is completed).  One paper is accepted for publication in “International Journal of Climatology ( IF:3.61) ( A Journal of Royal Meteorological Society, Great Britain, published by Wiley international).  Total 05 papers are presented in the conferences.  Five associate members and one Project Assistant are sponsored to participate in a workshop in Nagaland University  Detail scientific progress report is submitted.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 8

The Scientific Objectives of the Approved Project

(i) To study the spatial and temporal (seasonal and diurnal) variability of the morphological, microphysical and rainfall properties of the precipitating systems. (ii) To study the extreme weather systems in terms of these properties and atmospheric energetic. (iii) To develop a comprehensive “Precipitation Features & Atmospheric Data Base” with the help of satellite and ground based observations. (iv). To study the social impact of weather related disasters and its socio –psychological response.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 9

Executive Summary of the Scientific Results

The basic objective of the presented report is to study the spatial and temporal (seasonal and diurnal) variation of the convective systems in association with the rainfall and hail during the premonsoon and monsoon seasons. The scientific report consists of three parts namely (i) Seasonal and intra-seasonal variation of rainfall contribution by different types of convective systems over the south Asian region (ii). Diurnal variation of the occurrence of convective systems over north eastern part of India and adjoining region (iii) Detection of Hail Features (HFs) by satellite onboard microwave sensors over the north eastern part of India. The study is carried out over the south Asian region in general with special reference to the north eastern part of India and its adjoining area. Broadly study area consists of Indian subcontinent and surrounding oceanic region namely Bay of Bengal, Arabian sea and Indian ocean. The study is carried out by using the space borne sensors on Tropical Rainfall Measuring Mission (TRMM) a polar orbit satellite. The observations from TRMM-Precipitation Radar (PR) and TRMM Microwave Imager (37 GHz channels) are utilized for the present study. The observations by the TRMM sensor are supplemented by the observations from INSAT -3D onboard infrared thermal sensors. To study seasonal change, the premonsoon (March-May) and monsoon (June-September) months are considered. To study intra- seasonal change, the active and break days during the monsoon season are considered. The hail features as detected by the satellite onboard microwave sensors are analyzed during the premonsoon season. Overall results are summarized as follows. Chapter 1: Seasonal and intra-seasonal variation of rainfall contribution by different types of convective systems over the south Asian region

 Convective systems are classified into three categories namely, Large Convective Systems (LCSs), Deep Convective Systems (DCSs) and intense Convective Systems (ICSs). Rain characteristics from these three systems are studied in terms of rainfall contribution and rain intensity. The analysis is carried out during the premonsoon and monsoon seasons as well as during the active and break periods of the monsoonseason.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 10

 (i). Rainfall contribution by LCSs, in both the seasons, is predominantly over the oceanic region, with maximum over the northern Bay of Bengal and central Arabian Sea. Over the land during the premonsoon its maximum contribution is along the Himalayan foothills and during the monsoon its maximum contribution is over the central India region along with the Himalayan foothills. (ii). Rainfall contribution by DCSs, during the premonsoon, is predominantly over the oceanic region with maximum contribution over the central Bay of Bengal and central Arabian sea, whereas during the monsoon, it is predominantly over the land with maximum contribution over the Sindh region of Pakistan. (iii). Rainfall contribution by ICSs, in both the seasons, is predominantly a land phenomena. During the premonsoon the maximum contribution is over the eastern India and during the monsoon the maximum contribution is over the western Himalaya Indentation region.  (i). Rainfall contribution by LCSs over the land during the active monsoon days is predominantly over the core monsoon region, whereas during the break days rainfall contribution is shifted towards the Bangladesh plain, north east India and central & eastern Himalaya foothills region. Rainfall contribution by LCSs over ocean region is maximum over Bay of Bengal during the active days, and during the break days , it shifts to equatorial Indian oceanic region (ii). The nature of spatial variation of rainfall contribution by DCSs during active and break days are similar to LCSs. (iii). As far as contribution of rainfall from ICSs during active and break days is concerned, there is no appreciable shift in the region of the maximum rainfall during these days ( which is over the northern part of the East Coast of India and Western Himalaya Indentation). The pattern remain same in both the periods albeit with varying rainfall contribution.  The extreme rain intensity (irrespective of type of convective systems namely LCSs, DCSs and ICSs) is generally higher over the land compared to the ocean. During both the seasons, over the land, the most extreme rain intensity are found over the Himalaya-foothills (particularly eastern Himalaya foothills), northern part of east-coast of India, Gangetic West Bengal, Chota-Nagpur-plateau, Bangladesh-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 11

plain and Meghalaya-plateau albeit with relatively higher values during the premonsoon seasons. Over the ocean the most extreme rain intensity is observed over the Bay of Bengal and Arabian sea.  The characteristics of ICSs have strong regional variability. The ICSs are stronger during the premonsoon compared to the monsoon seasons and also ICSs over the land are stronger compared to the ocean. Over land, during the premonsoon and monsoon seasons, the maximum radar reflectivity at mixed-phase region (at 7 Km height ) are observed over the northern part of East Coast of India ( 4 dB ) and estern-Himalaya Indentation ( 43dBZ) respectively , Over the ocean, during the premonsson as well as the monsoon season, the maximum radar reflectivity at mixed-phase region are observed over the Bay of Bengal ( 44 dBZ and 40dBZ).  The preferred locations of ICSs are associated with a relatively higher value of the IWC in the upper part of the mixed-phase region. Over the ocean, with insignificant occurrence of ICSs, the parameter has relatively low values in the mixed-phase region.

Chapter 2: Diurnal variation of the occurrence of convective systems over north eastern part of India and adjoining region

 The observation from TRMM as well as INSAT 3D suggests that there is a significant regional variability in the diurnal characteristics of the convective systems over the study region. The foothills region of Himalya including the Assam valley have maximum occurrence during night to morning hours, whereas the plain region of Bangladesh, Gangetic west Bengal, northern part of east coast have maximum occurrence during afternoon to night hours.

Chapter 3: Detection of Hail Features (HFs) by satellite onboard microwave sensors over the north eastern part of India.

 On the basis of threshold value of the Polarization corrected temperature of 37

GHz channels (PCT37), hail features (HFs) are classified into three categories namely, T-1 (with hail detection probability of 24%) , T-2 (with hail detection

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 12

probability of 45% ) and T-3 (with hail detection probability of 70%).  T-1, T-2, and T-3 HFs have about 35, 45, and 49 dBZ radar reflectivity at 9 Km height (within the mixed-phase region). The nature of vertical structures of reflectivity for T-1, T-2, and T-3 HFs. correlates well with hail reporting at ground. The strongest vertical structure of T-3 HFs indicates the strong updraft and large ice particles presence in the mixed - phase region. It also represents that the strong mixed-phase microphysical processes (i.e., freezing of raindrops and riming) are involved for production of hail/graupel.  The two stations in the plain region, namely Agartala and Dhubri, detected the maximum occurrence of most severe HFs (type T-3). The stations in the valley regions namely Tezpur, Mohanbari, North Lakhimpur, Pasighat, and Imphal have not detected the most severe HFs. The spatial variability in the HFs is amply supported by the vertical profiles of reflectivity and its value at mixed phase region.  During the premonsoon, the maximum occurrence of HFs is found in April. The occurrence of HFs is minimum in March.  HFs show strong diurnal variation. with maximum occurrence during the afternoon hours  Occurrence of HFs show noticeable year to year variation. There is a decreasing trend during the period 1998 - 2013. The trend values are -0.70, -0.44, -0.14 and - 0.10 for T-1, T-2 and T-3 HFs and ground reporting respectively.

Overall significant spatiotemporal variability is observed, with respect to rainfall contribution by different types of convective systems . The study provides a better insight into the spatiotemporal characteristics of convective systems over the region. The present results will be helpful for the better understanding of severe weather and also to improve rainfall estimation from satellite on-board sensors over the region. The study will also be useful for hydro-meteorological applications, particularly for flood forecasting and landslide hazards.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 13

SCIENTIFIC PROGRESS REPORT

(2016)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 14

Chapter 1

Seasonal and intra-seasonal variation of rainfall contribution by different types of convective systems over the south Asian region

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 15

Seasonal and intra-seasonal variation of rainfall contribution by different types of convective systems over the south Asian region

1.1 Introduction

Convective systems play a crucial role in the hydrological cycle and energy budget of the planet and influence its climate variability. Spatio-temporal variation of their population, along with their vertical extent and microphysical properties influence the variation in the earth’s water cycle and also pose a challenge to parameterize the rainfall processes. Intense convective systems are often associated with strong lightning and hails in the storm updraft region (Dye et al., 1989, Saunders et al., 1991; Ushio et al., 2001; Petersen et al., 2005; Luo et al., 2011; Wiens et al., 2005).

The significant part of the South Asia is affected by hazardous convective systems during the premonsoon season (Zipser et al., 2006; Romatschke et al., 2010; Choudhury et al., 2015). The South-West monsoon system, an ensemble of convective systems, is one of the largest meteorological system (Gadgil, 2003), which affects the socioeconomic condition of a sizeable population of the South Asia. The launch of Tropical Rainfall Measuring Mission (TRMM; Simpson et al., 1988) in 1997 provided a great opportunity to study the properties of convective systems. Significant studies on the properties of these systems have been carried out by using the TRMM observations at a global scale (Nesbitt et al., 2000; Alcala and Dessler, 2002; Nesbitt and Zipser, 2003; Boccippio et al., 2005; Liu and Zipser, 2005; Cecil et al., 2005; Nesbitt et al., 2006; Zipser et al., 2006; Nesbitt and Anders 2009; Liu, 2011; Liu et al., 2011; Liu et al., 2012; Xu and Zipser 2012; Houze et al., 2015; Hamada et al., 2014; Hamada et al., 2015). By using the TRMM observations, several studies of this nature have also been carried out over the South Asia. By using the TRMM-PR near surface reflectivity observations, Bhatt and Nakamura (2005) showed that, during the monsoon season large portion of the rainfall is concentrated over the south facing slopes of the Himalaya. Similarly, during the monsoon season, by using the vertical profile of reflectivity from the TRMM-PR over the Himalaya region, Houze et al. (2007) reported that 40 dBZ echo top > 10 km and 40 dBZ echo area > 1000 km2 occurs preferentially in the north western concave Indentation of the barrier. Islam and

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 16

Uyeda (2008) studied TRMM-PR derived vertical profile of rain intensity in and around Bangladesh during various rainy seasons. They reported that, the premonsoon rainfall is characterized by higher rain intensity and echo top height compared to the monsoon and post monsoon seasons. With the help of high resolution numerical simulation model and observation from TRMM-PR, Medina et al. (2010) studied the terrain and land cover effects on the summer monsoon convection in the Himalaya region. They reported that, intense convective echoes occur in the western Himalaya region, whereas broad stratiform echo occurs in the eastern-Himalaya. Through numerical analysis they showed that these variations are a result of region specific orographically modified flows and land surface flux feedbacks. By using the TRMM-PR data Romatschke et al. (2010) reported that the preferred location of deep as well as wide convective cores changes distinctly from India’s East-coast in the premonsoon to Western-Himalaya-foothills in the monsoon. By using the same data set and time period, contribution of different size of the convective systems to the precipitation were studied during the premonsoon (Romatschke and Houze, 2011 a) and monsoon (Romatschke and Houze, 2011 b). During the premonsoon, it was reported that over the land, most of the rain falls from the medium sized systems. Over the Bay-of-Bengal, the dominant systems are larger. It was shown that, during the monsoon along the Western-Himalaya, precipitation falls mainly from the small but from the highly convective nature. Medium systems are favoured over the east Arabian- Sea and large systems are favoured over the Bay-of-Bengal. By using the TRMM sensors, Qie et al. (2014) reported that DCSs (20 dB echo top ≥ 14 km) and IDCSs (40 dB echo top ≥ 10 km) are most frequent over the southern Himalayan front (SHF) especially in the western most SHF. The DCSs over the Tibetan-plateau are relatively weak in convective intensity and small in size but occur frequently. The Oceanic DCSs possess the tallest cloud top and the largest in size but their convective intensity is significantly weaker. By using the TRMM observations, Choudhury et al. (2015) reported that the most intense MCSs occurs during the premonsoon over the Chota-Nagpur- Plateau and adjoining region with markedly reduced intensity in monsoon. They further pointed out that the Eastern-Himalaya and adjoining area consist of relatively weak MCSs and the region of most intense MCSs is associated with the large ice particles and ice

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 17 water content compared to the region of weak MCSs. Bhat and Kumar (2015) studied cumulonimbus towers and intense convective cells by using the vertical reflectivity profiles from the TRMM-PR. They reported that the frequency of the occurrence of cumulonimbus towers (CbTs; 20 dB radar echo ≥ 12 km height and at least 9 km thick) is highest over the foothills of Himalaya, plain of northern India and Bangladesh and minimum over the Arabian-Sea and equatorial-Indian-ocean. They reported marginal land, ocean contrast for CbTs.

The basic objective of the study is to investigate the f rainfall contribution by various type of convective systems such as large, deep and intense over the South Asian region

1.2 Data and Methodology:

To study the spatio-temporal variability of rainfall contribution of Convective Systems (CSs) over the South Asian region (Figure-1.1),

35°N R-1 30°N R-3 26°N R-2 22°N R-5 18°N R-7 R-4 14°N

10°N Latitude 6°N 2°N

2°S R-6

6°S 10°S 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude

Figure 1.1 : Physiographic map of South Asian region.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 18

For the present work 17 years (1998- 2014) of TRMM precipitation radar (PR) observation during the premonsoon (March - May) and the monsoon (Jun - September) seasons are utilized. Moreover, the intraseasonal change in rainfall characteristics of CSs is examined from the monsoon active days to break days. Total numbers of active days (108) and break days (148) from 1998 to 2014 are collected from as given by Pai et al (2015). PFs are defined by grouping the contiguous pixels with nonzero near-surface rainfall rate from TRMM PR 2A25 v7 data product (Liu et al., 2008).

CSs are categorized as a large, deep and intense. They are defined as follows

(i) Large Convective Systems (LCSs): When near-surface rainfall area of CSs are greater than 10000 Km2 (Hirose and Nakamura, 2005),

(ii). Deep Convective Systems (DCSs) : When maximum height of 20 dBZ echoes within the CSs are greater than 12 Km ( Xu., 2013)

(iii) Intense Convective Systems (ICSs) : When maximum height of 40 dBZ e echoes within the CSs are greater than 7 Km (Xu., 2013) respectively.

Rainfall contribution by each type of CSs are presented in each 2o x 2o grid. Rainfalls contributions are calculated by dividing rain volume of each CSs to total rain volume in 2o x 2o region.

To study the characteristics of extreme rain intensity from CSs which are defined as CSs in which maximum near-surface rainfall intensity (mm hr-1) are higher than the corresponding 99.9th percentile on a 2o x 2o grid.

1.3 Results:

1.3.1 Spatial and seasonal variation of rainfall contribution:

Table- 1.1 shows the regional and seasonal variability of rainfall contribution by each type of CSs. Regionally, LCSs and DCSs both show that a larger proportion of

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 19 rainfalls come from them over the ocean than over the land. Seasonally, rainfall contribution by LCSs shows a significant seasonal variability over land compared to ocean. Over land, LCSs (DCSs) contribute about 13% (~ 12%) rainfalls during the premonsoon and, it changed to about 18% (~ 16%) during the monsoon. Over ocean, LCSs (DCSs) contributes about 34% (~ 32%) rainfalls during the premonsoon and, it changed to about 34% (~ 22%) during the monsoon. Rainfall contribution by ICSs shows significant regional variability, it is generally higher over land compared to ocean during each seasons. Over the land During the premonsoon and monsoon , about 11% and 9% of rainfall contribution is from ICSs respectively.

Table-1.1 : Seasonal variation in rainfall contribution by each type CSs over land and ocean.

Land Ocean ______P-M M P-M M ______All CSs 796053 1948973 1523190 2220862 LCSs 13.15% 17.68% 34.12% 34.04% DCSs 11.55% 16.37% 32.13% 21.96% ICSs 10.64% 8.67% 4.31% 3.41% ______P-M: Premonsoon, M: Monsoon

During the premonsoon and the monsoon seasons the spatio-temporal variation of rainfall contribution by LCSs is shown in Figure 1.2 (a, b) respectively. Over land, premonsoon rainfall is mostly dominated by isolated thunderstorm over central-India, but in monsoon significant rainfall comes from organized LCSs. Rainfall contribution by LCSs is increased (decreased) in monsoon over the Meghalaya-plateau, equatorial-Indian- ocean and south-Myanmar-coast (west-coast of India). Over all, it is observed that, LCSs contribute large portion of rainfall over the north Bay-of-Bengal than those in high rainfall regions, such as equatorial-Indian-ocean, Meghalaya-plateau, Eastern-Himalaya-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 20 foothills, and west-coast of India during each seasons. It is interesting to see that, over north Arabian-sea (near Gujarat and Karachi coast), rainfall mostly comes from small CSs in the premonsoon and in the monsoon it is from LCSs.

(a) Premonsoon (b) Monsoon

Large PFs rainfall contribution (%) Large PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 20 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.2: shows the spatial variability of rainfall contribution by LCSs (a) during the Premonsoon and (b) during the Monsoon.

During the premonsoon and the monsoon seasons the spatio-temporal variation of rainfall contribution by DCSs is shown in Figure 1..3 (a, b) respectively. Over land, local rainfall received by DCSs shows the significant spatial variability from premonsoon to monsoon seasons. During premonsoon, DCSs contribute most of local rainfall along east- coast of India, Gangetic West Bengal, Chota-Nagpur-plateau and Bangladesh-plain and during the monsoon, there is less rainfall over Thar-desert but most of the rainfall over this region is from few DCSs. Over the ocean, DCSs contribute maximum percentage of local rainfall over north Bay-of-Bengal and central Arabian-sea during premonsoon (Figure 1.3.a) and during monsoon, it is observed over north Arabian-sea (near Gujarat and Karachi coast) (Figure 1.3.b). During premonsoon, over west-coast of India and southeast Bay-of-Bengal (near south-Myanmar-coast) significant proportion of local rainfalls receive

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 21 from DCSs, but in monsoon it is from shallow CSs.

(a) Premonsoon (b) Monsoon

Deep PFs rainfall contribution (%) Deep PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.3: shows the spatial variability of rainfall contribution by DCSs (a) during the Premonsoon and (b) during the Monsoon.

During the premonsoon and the monsoon seasons the spatio-temporal variation of rainfall contribution by ICSs is shown in Figure 1.4 (a, b) respectively. The ICSs contribute a larger proportion of rainfall over land than ocean. Over land, in the premonsoon, maximum rainfall contribution is over north of east-coast of India, Chota- Nagpur-plateau, and adjoining Bangladesh-plain (Figure 1.4.a) and in monsoon, which is shifted in western-Himalaya-Indentation region (Figure 1.4.b). It is interesting to see that CSs which are intense but not deep, contribute higher portion of rainfall than CSs those are deep but not intense over western-Himalaya-Indentation region during both seasons.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 22

(a) Premonsoon (b) Monsoon

Intense PFs rainfall contribution (%) Intense PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.4: shows the spatial variability of rainfall contribution by ICSs (a) during the Premonsoon and (b) during the Monsoon.

Overall it is observed that the rainfall contribution by LCSs, in both the seasons, is predominantly over the oceanic region, with maximum over the northern Bay of Bengal and central Arabian sea. Rainfall contribution by DCSs, during premonsoon, is predominantly over the ocean with maximum contribution over the central Bay of Bengal and central Arabian sea, whereas rainfall contribution by DCSs, during the monsoon, is both over land and ocean albeit with maximum contribution over the land, the Sindh region of Pakistan. Rainfall contribution by ICSs, in both the season, is predominantly over the land with maximum over the eastern India during the premonsoon and maximum over the western Himalaya Indentation region monsoon.

1.3.2. Intra-seasonal variation of rainfall contribution:

The previous section shows clear spatial contrasts in rainfalls contribution from each type of CSs over land and ocean. This section compares rainfall contribution by each type of CSs, during the monsoon active and break days. The Table-1.2 shows the

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 23 significant intra-seasonal change of rainfall contribution of CSs types over land and ocean. The intra-seasonal change suggests that over the land, larger percentage of rainfall contribution from each type of CSs is observed during the monsoon active days compared to monsoon break days. On the contrary, ocean region experiences an increase in rainfall contribution by LCSs and DCSs during monsoon break days, whereas rainfall contributions by ICSs is relatively larger during the active monsoon days.

Table-1.2: Intra-seasonal variation of rainfall contribution of each type of CSs over land and ocean.

Land Ocean ______Active Break Active Break ______All PFs 120088 151688 109260 164018

LCSs 23.65% 14.72% 27.72% 34.48% DCSs 20.80% 14.00% 16.80% 21.20% ICSs 20.90% 19.00% 6.00% 3.25% ______

Figure 1.5 (a, b) shows the spatial distribution of rainfall contribution by LCSs during active and break days, respectively. Over land, the central-India and north western part of India experiences the greater decrease in local rainfall contribution by LCSs in break monsoon days. Over Ocean, over south Bay-of-Bengal, south equatorial-Indian- ocean and southeast Arabian-sea (near west-coast), a significant increase in rainfall contribution by LCSs is observed in monsoon break days.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 24

(a) Monsoon (Active days) (b) Monsoon (Break days)

Large PFs rainfall contribution (%) Large PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.5: shows the spatial variability of rainfall contribution by LCSs (a) during Active days and (b) during Break days.

Rainfall contribution by DCSs during active and break days, are shown in figure 1.6 (a, b) respectively. Over land, the central Himalaya-foothills region experiences much smaller intra-seasonal changes in DCSs rainfall contribution compared to eastern Himalaya-foothills regions. Though, the eastern Himalaya-foothills regions receive maximum rainfall during break monsoon days, rainfall contribution by DCSs is significant in monsoon active days. Even though, during monsoon active days, monsoon trough locates over central-India and north western part of India, DCSs contribute relatively larger portion of local rainfall over the north western part of India.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 25

(a) Monsoon (Active days) (b) Monsoon (Break days)

Deep PFs rainfall contribution (%) Deep PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.6: shows the spatial variability of rainfall contribution by DCSs during (a) d Active days and (b) Break days.

Rainfall contribution by ICSs during active and break days, are shown in figure 1.7 (a, b) respectively. Spatial distribution on rainfall contribution by ICSs does not show any significant change. Majority of rainfall contribution by ICSs persist over western- Himalaya-Indentation region. Over northern part of east-coast and southeast peninsula experience a reasonably increase in local rainfall contribution by ICSs in break monsoon days.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 26

(a) Monsoon (Active days) (b) Monsoon (Break days)

Intense PFs rainfall contribution (%) Intense PFs rainfall contribution (%) 35°N 90 35°N 90

30°N 80 30°N 80 26°N 26°N 70 70 22°N 22°N 18°N 60 18°N 60 14°N 14°N 50 50

10°N 10°N

Latitude Latitude 6°N 40 6°N 40

2°N 30 2°N 30 2°S 2°S 20 20 6°S 6°S

10°S 10 10°S 10 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.7: shows the spatial variability of rainfall contribution by ICSs (a) during Active days and (b) during Break days.

Overall, the rainfall contribution by LCSs over the land during active monsoon days is predominantly over the core monsoon region, whereas during the break days rainfall contribution is shifted towards the Bangladesh plain, north east India and central & eastern Himalaya foothills region. Rainfall contribution by LCSs over ocean region with maximum at Bay of Bengal during the active days shift to maximum at equatorial Indian ocean during the break days. The nature of spatial variation of rainfall contribution by DCSs during active and break days are similar to LCSs. As far as contribution of rainfall from ICSs during active and break days is concerned, there is no appreciable shift in the region of maximum rainfall during these days.

1.3.3 Spatial and seasonal variation of extreme rain intensity :

Spatial variation of extreme rain intensity (irrespective of the type of convective systems namely LCSs, DCSs and ICSs) during the premonsoon and monsoon seasons are shown in figure 1.8 (a, b) respectively. The extreme rain intensity are generally higher

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 27 over land compared to ocean. Most extreme rain intensity are found over the land along Himalaya-foothills, northern east-coast, Gangetic West Bengal, Chota-Nagpur-plateau, Bangladesh-plain and Meghalaya-plateau during both seasons. Significant rain intensity is also observed the coastal regions, including south west-coast of India and Myanmar-coast during both the seasons. Extreme rain intensity is higher during the premonsoon compared to the monsoon. Over the ocean, in both the seasons, extreme rainfall is observed mainly over the Bay of Bengal and Arabian sea .

(a) Premonsoon (b) Monsoon

99.9-ile value of R ( mm/h ) 99.9-ile value of R ( mm/h ) max max 35°N 240 35°N 240 220 30°N 220 30°N 200 26°N 200 26°N 180 22°N 180 22°N 18°N 160 18°N 160 14°N 140 14°N 140 10°N 120 10°N 120 Latitude Latitude 6°N 100 6°N 100 2°N 80 2°N 80 2°S 60 2°S 60 6°S 40 6°S 40 10°S 20 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 10°S 20 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

Figure 1.8: shows the spatial variation of extreme rain intensity during the (a) Premonsoon and (b) Monsoon.

The extreme rain intensity is generally higher over the land compared to the ocean. Most extreme rain intensity (irrespective of type of convective systems namely LCSs, DCSs and ICSs) are found over the land along Himalaya-foothills, northern east-coast, Gangetic West Bengal, Chota-Nagpur-plateau, Bangladesh-plain and Meghalaya-plateau during both the seasons albeit with relatively higher values during the premonsoon seasons.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 28

1.3.4 Regional variability of the characteristics of convective systems associated with extreme rain intensity

In this subsection, the characteristics of the convective systems (e.g., size and vertical extent ) associated with the extreme rain intensity, are examined in four sub-regions over land and in three sub-regions over ocean (Table-1.3). Over the land the selected sub- regions are (i) R-1 (280N – 350N & 690E – 750E) (ii) R-2 (200N – 260N & 830E – 900E), (iii) R-3 (240N – 300N & 910E – 980E) and (iv) R-4 (90N – 180N & 740E – 790E). Over the ocean the selected sub-regions are (i) R-5 (140N – 220N & 870E – 940E) (ii) R-6 (100S – 00N & 750E – 820E), (iii) R-7 (120N – 180N & 650E – 720E).

Over land, during the premonsoon, maximum fraction of extreme rain intensity from the mesoscale convective systems (MCSs; near-surface rainfall area greater than 2 2000 Km ) are observed over the R-1 (~ 88%), followed by the R-2 (~ 73%), the R-3 (~

54%) and the R-4 (~ 40%) and during monsoon, it is observed over the R-3 (~ 68%), followed by the R-4 (~ 61%), the R-2 (~ 56%) and the R-1(~ 48%). Over ocean, during the premonsoon, maximum fraction of extreme rain intensity from MCSs are observed over the R-7 (~ 100%), followed by the R-6 (~ 97%) and the R-5 (~ 91%) and during monsoon, it is observed over the R-5 (~ 98%), followed by the R-7 (~ 95%) and the R-6 (~ 94%).

Over land, during the premonsoon, maximum fraction of extreme rain intensity by storms having strong mixed-phase process (40dBZ echo top > 7 Km, i.e., freezing of raindrops and riming) (Xu and Zipser, 2012) are observed over the R-2 (~ 80%), followed by the R-1 (~ 46%), the R-4 (~ 20%) and the R-3 (~ 12%) and during monsoon, maximum fraction is observed over the R-1 (~ 47%), followed by the R-2 (~ 23%), the R-

3 (~ 6%) and the R-4 (~ 2%). Over ocean, during the premonsoon, maximum fraction of extreme CSs by storms involving strong mixed-phase process are observed over the R-7

(~ 64%), followed by the R-5 (~ 59%), and the R-6 (~ 1%) and during monsoon, maximum fraction is observed over the R-5 (~ 40%), followed by the R-7 (~ 16%) and the R-6 (~ 0%).

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 29

Table-1.3: Regional variability of extreme PFs properties over the seven sub-regions.

______R-1 R-2 R-3 R-4 R-5 R-6 R-7 ______P-M, M ______Occurrence 52, 60 15, 75 65, 143 15, 61 22, 114 109, 126 11, 90

Area> 2000 km2 (%) 88.5, 48.3 3.3, 56.0 53.8, 68.5 40.0, 60.6 91.0, 98.2 97.2, 94.4 100.0, 95.5

40dBZ >7 Km (%) 46.1, 46.7 80.0, 22.7 12.3, 5.6 20.0, 1.6 59.0, 40.3 1.0, 0.0 63.6, 15.5 ______P-M: Premonsoon, M: Monsoon

The vertical profiles of maximum radar-reflectivity factor for the extreme rain intensity in the sub-regions over land and ocean during the premonsoon and monsoon seasons are shown in figure 1.9 (a, b). Over land, during the premonsoon, maximum radar reflectivity at mixed-phase region ( 7 Km) are observed over the R-2 (~ 48dBZ), followed by the R-1 (~ 40dBZ), the R-4 (~ 38dBZ) and the R-3 (~ 32dBZ) and during monsoon, it is observed over the R-1 (~ 43dBZ), followed by the R-2 (~ 38dBZ), the R-3

(~ 34dBZ) and the R-4 (~ 29dBZ). Over ocean, during the premonsoon, maximum radar reflectivity at mixed-phase region are observed over the R-5 (~ 44dBZ), followed by the

R-7 (~ 39dBZ), and the R-6 (~ 33dBZ) and during monsoon, maximum radar reflectivity at mixed-phase region are observed over the R-5 (~ 40dBZ), followed by the R-7 (~

32dBZ), and the R-6 (~ 32dBZ)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 30

(a) Premonsoon (b) Monsoon

Vertical profile of extreme PFs Vertical profile of extreme PFs

16 R-1 16 15 R-1 R-2 15 14 R-2 R-3 14 13 13 R-3 R-4 R-4 12 12 R-5 R-5 11 R-6 11 10 R-6 R-7 10 R-7 9 9 8 8

Height(Km) 7 Height(Km) 7 6 6 5 5 4 4 3 3 2 2 20 25 30 35 40 45 50 55 60 20 25 30 35 40 45 50 55 60 Max Radar Reflectivity (dBZ) Max Radar Reflectivity (dBZ) Figure 1.9: shows the vertical profiles of maximum radar-reflectivity factor for the extreme rain intensity during the (a) Premonsoon and (b) Monsoon.

Further, the spatial variation of the most extreme value of ETH40dBZ (max) (at 99 percentile; top 1%) during the premonsoon and monsoon seasons are provided in Figure 1.10 (a, b) respectively. During the premonsoon, over the land, its value varies in the range of 4-16 km with the maximum value over the Chota-Nagpur-plateau, Bangladesh- plain, central and northern part of the East-coast of India and minimum over the Tibetan- plateau. Over the ocean its value varies in the range of 8-14 km with maximum value over the northern Bay-of-Bengal near the coastal area and minimum over the Indian- ocean. During the monsoon season, over the land its value varies in the range of 8-16 km with its maximum value shifted to the Western-Himalaya-Indentation and minimum value

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 31 over the southern Myanmar and adjoining region. During this season, the region of Chota-Nagpur-Plateau, Bangladesh-plain and northern East-coast are no longer the region of maximum value. Over the ocean, its value varies in the range of 7-11 km, with a maximum value over the Bay-of-Bengal very near to the eastern coast and the minimum value over the Indian-ocean. It is observed that over the Bay-of-Bengal, its value decreases considerably compared to the premonsoon season. Overall significant land, ocean contrast is observed in their values,

(a) 99th percentile values of ETH ( km ) (b) 99th percentile values of ETH ( km ) 40dBZ(max) 40dBZ(max) 35°N 35°N 4 8 4 6 12 12 9 10 11 10 12 8 8 8 10 11 14 4 8 16 14 11 10 10 30°N 12 10 10 30°N 16

14

9 14 16 13 13 13 12 14 12

9 12 12 13 11 10 26°N 26°N 15 14 12 13 14 16 12 11 13

10 14 16 16 9 22°N 22°N 10 12 12 14 10 12 13

16 11 14 8 4 14 18°N 9 11 8 18°N 8 14 8 12 10 9 9 9 9 16 13 16 10 14 10 10 12 14°N 11 12 14 9 8 10 11 14°N 8 12 10 11 8 13 12 12 9 9 8 12 10°N 10°N 10 14 10 Latitude 7

8 8 8 Latitude 8 11 9 14 10 7 10 7 6°N 8 6°N 7 7 7 8 2°N 2°N 7 7 7 2°S 7 7 2°S 8 7

6 6°S 7 6°S 6 10°S 10°S 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

(a) (b)

th Figure 1.10: Spatial variation of 99 percentile value of ETH40dBZ(max) during the (a) pre monsoon and (b) monsoon.

1.3.5 Spatial & seasonal variation of ice microphysical properties.

To have a better insight of the regional and seasonal variability of the, ice microphysical properties, the vertical profiles of the ice water content ( IWC) are considered for further analysis. During the premonsoon and monsoon seasons, the contour plots of median value of the IWC for height vs longitude cross section over the o o o o o o o o o o o R1 (22 N–35 N, 65 E–100 E), R2 (6 N–22 N, 65 E–100 E) and R3 (10 S – 6 N, 65 E–

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 32

100o E) are shown in Figure 1.11 (a, b) and 1.11 (c, d) and 1.11 (e, f) respectively. To study the microphysical process in the mixed phase region, which is conducive for development of severe weather (lightning and hail features in the convective systems), a reference height of 10 km is selected. Over the R1, during the premonsoon, significant value of the IWC is observed at around 10 km (the mixed-phase region) in between 85o- 95o E, which coincide with the longitude range of the Meghalaya-plateau, Chota-Nagpur- Plateau and Bangladesh-plain. This is the region where the higher occurrence of lightning and hail features are significantly prevalent during the premonsoon (Roy et al. 2017). During the monsoon significant value of the IWC is observed at around 10 km in between the 65o-70o E, which coincide with the longitude range of the Western-Himalaya- Indentation and adjoining area, the region where ICSs are most prevalent along with a higher occurrence of lightning and hail features (Roy et al. 2017) . During premonsoon to monsoon, shifting of the higher amount of the IWC in the mixed-phase region from the Meghalaya-plateau, Chota-Nagpur-plateau and Bangladesh-plain to the Western- Himalaya-Indentation is in agreement with the seasonal shifting of the preferred location of lightning and hail features (Roy et al. 2017).Over the R2, during the premonsoon, significant value of the IWC is observed at around 10 km (in the mixed-phase region) in between the 75o-90o E, which coincide with the longitude range of the southern part of West-coast of India, Sri Lanka and the northern part of East-coast of India. It is the region where occurrence of ICSs is maximum during the premonsoon, particularly towards the northern part of East-coast. During the monsoon, lower value of the IWC in the mixed- phase region over the same longitude range is consistent with the reduced occurrence of

ICSs with lightning and hail features over this region (Roy et al. 2017) . Over the R3, during the premonsoon and monsoon, there is relatively low value of the IWC at 10 km (in the mixed phase region) compared to the other two regions. The lower value of the IWC at the mixed-phase heights may be attributed to the weak updrafts over the oceanic region (Zipser and Lutz, 1994). Overall, the results are in agreement with the spatio- temporal variation of IDCSs with lightning and hail features (Roy et al. 2017). The preferred locations of ICSs are associated with a relatively higher value of the IWC in the upper part of the mixed-phase region. Over the ocean, with insignificant occurrence of

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 33

ICSs, the parameter has relatively low values in the mixed-phase region.

(b) R ( average over 220 N - 350 N ) (a) R ( average over 220 N - 350 N ) 1 1 20 20 18 18 50 50 16 100 16 100 400 100 200 150 100 14 250

14 500 50 150 100 150 250 300 250 12 12 150 350 350 350 250 250 400 300 10 350 450 400 300 250 10 350 400 500 450 8 250

8 500 400 150 500 450 250 Height) km ( Height) km ( 250 50 6 200 6 50 4 4 2 2 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude b) Longitude (a)

(d) R ( average over 60 N - 220 N ) 2 (c) R ( average over 60 N - 220 N ) 2 20 20 18 50 18 50 16 100 50 16 100 250 14 250 100 14 200 350 150 150 12 200 250 250 300 12 250 400 10 300 300 400 400 250 350 300 250 300 300 10 350 300 400 8 250

400 Height) km ( 300 100 150 8 200

Height) km ( 6 50 100 150 6 50 4 4 2 2 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E Longitude (c ) Longitude (d) (e) R ( average over 100 S - 60 N ) (f) R ( average over 100 S - 60 N ) 3 3 20 20 18 18 16 50 100 100 100 16 50 14 50 14 150 100 12 250 200 500 300 350 12 150 10 450 250 400 350 400 300 400 10 8 350 350 350

Height) km ( 400 350 150 200 8 350

6 50 Height) km ( 250 100 4 6 50 2 4 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E 2 Longitude 65°E 70°E 75°E 80°E 85°E 90°E 95°E 100°E (e) (f) Longitude -3 Figure 6: The height longitudinal cross section of IWC (mg m ) (a) over the R1 during the premonsoon (b) over the R1 during the monsoon (c) over the R2 during the premonsoon (d) over the R2 during the monsoon (e) over the R3 during the premonsoon and (f) over the R3 during the monsoon.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 34

1.4 Summary and conclusion

 Convective systems are classified into three categories namely, Large Convective Systems (LCSs), Deep Convective Systems (DCSs) and intense Convective Systems (ICSs). Rain characteristics from these three systems are studied in terms of rainfall contribution and rain intensity. The analysis is carried out during the premonsoon and monsoon seasons as well as during the active and break periods of the monsoonseason.  (i). Rainfall contribution by LCSs, in both the seasons, is predominantly over the oceanic region, with maximum over the northern Bay of Bengal and central Arabian Sea. Over the land during the premonsoon its maximum contribution is along the Himalayan foothills and during the monsoon its maximum contribution is over the central India region along with the Himalayan foothills. (ii). Rainfall contribution by DCSs, during the premonsoon, is predominantly over the oceanic region with maximum contribution over the central Bay of Bengal and central Arabian sea, whereas during the monsoon, it is predominantly over the land with maximum contribution over the Sindh region of Pakistan. (iii). Rainfall contribution by ICSs, in both the seasons, is predominantly a land phenomena. During the premonsoon the maximum contribution is over the eastern India and during the monsoon the maximum contribution is over the western Himalaya Indentation region.  (i). Rainfall contribution by LCSs over the land during the active monsoon days is predominantly over the core monsoon region, whereas during the break days rainfall contribution is shifted towards the Bangladesh plain, north east India and central & eastern Himalaya foothills region. Rainfall contribution by LCSs over ocean region is maximum over Bay of Bengal during the active days, and during the break days , it shifts to equatorial Indian oceanic region (ii). The nature of spatial variation of rainfall contribution by DCSs during active and break days are similar to LCSs. (iii). As far as contribution of rainfall from ICSs during active

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 35

and break days is concerned, there is no appreciable shift in the region of the maximum rainfall during these days ( which is over the northern part of the East Coast of India and Western Himalaya Indentation). The pattern remain same in both the periods albeit with varying rainfall contribution.  The extreme rain intensity (irrespective of type of convective systems namely LCSs, DCSs and ICSs) is generally higher over the land compared to the ocean. During both the seasons, over the land, the most extreme rain intensity are found over the Himalaya-foothills (particularly eastern Himalaya foothills), northern part of east-coast of India, Gangetic West Bengal, Chota-Nagpur-plateau, Bangladesh- plain and Meghalaya-plateau albeit with relatively higher values during the premonsoon seasons. Over the ocean the most extreme rain intensity is observed over the Bay of Bengal and Arabian sea.  The characteristics of ICSs have strong regional variability. The ICSs are stronger during the premonsoon compared to the monsoon seasons and also ICSs over the land are stronger compared to the ocean. Over land, during the premonsoon and monsoon seasons, the maximum radar reflectivity at mixed-phase region (at m height ) are observed over the northern part of East Coast of India ( 4 dB ) and Western-Himalaya Indentation ( 43dB ) respectively , Over the ocean, during the premonsson as well as the monsoon season, the maximum radar reflectivity at mixed-phase region are observed over the Bay of Bengal ( 44 dB and 40dBZ).  The preferred locations of ICSs are associated with a relatively higher value of the IWC in the upper part of the mixed-phase region. Over the ocean, with insignificant occurrence of ICSs, the parameter has relatively low values in the mixed-phase region.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 36

Chapter 2

Diurnal variation of the occurrence of

convective systems over north eastern

part of India and adjoining region

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 37

Diurnal variation of the occurrence of convective systems over north eastern part of India and adjoining region

2.1 Introduction

The mechanism of the rainfall peak associated with overland afternoon instability are well understood and observed, however the diurnal cycle of rainfall associated with MCSs remain poorly understood and modelled (Nesbitt and Zipser 2003). Most of the previous work comparing the tropical rainfall diurnal cycle over land and ocean surface agree that the amplitude of the diurnal cycle over continent is larger than that over the open ocean. (Gray and Jacobson 1977). However studies conducted over different region of the tropics, have found significant differences in the characteristics of the diurnal cycle. Over land many studies using surface rain accumulation ( (Oki and Musiake 1994; Dai et al. 1999) and surface weather report of precipitation frequency (Dai 2000) link the timing of the diurnal precipitation frequency maximum to afternoon boundary layer destabilization caused by day time insolation. However many studies noted that there are land area with midnight to early morning maxima of the precipitation , which may be linked to local effects such as complex terrain and see breeze circulations (Oki and Musiake 1994; Yang and Slingo 2000) or the long nocturnal life cycle of MCSs ( Dai et al. 1999). There is no comprehensive study of the diurnal variation of the convective activity over the study region. In the present section an attempt is being made to study the diurnal variation of the occurrence of MCSs and DCSs over the study region.

2.2 Result

2.2.1 Spatiotemporal variation of occurrence of MCSs during the premonsoon and monsoon season (Diurnal variation)

The spatio-temporal variation of the occurrence of MCSs during during 0-6, 6- 12, 12- 18 and 18-00 mean standard time (MST) in the pre-monsoon seasons are shown in Figure 2.1 (a-d) respectively. It is observed that during the premonsoon season the

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 38 occurrence MCSs is predominantly over the eastern Himalaya foothills and Meghalaya Plateau, with the exception of 18-00 MST, when the MCSs occurrence is also observed over the plain region Bangladesh and Gangetic west Bengal.

MCSs population (Preonsoon) ( 00 - 06 MST ) 30°N MCSs population (Preonsoon) ( 06 - 12 MST ) 30°N

2

4

6 28°N 4 6 2 6 28°N 2 4

2 4 4 2 4 2 2 2 26°N 2 2 2 2

2 2 2 26°N

2 2 4

Latitude 24°N Latitude 24°N 2

22°N 22°N

20°N 80°E 85°E 90°E 95°E 100°E 20°N Longitude 80°E 85°E 90°E 95°E 100°E Longitude

(b) (a)

MCSs population (Preonsoon) ( 12 - 18 MST ) MCSs population (Preonsoon) ( 18 - 00 MST ) 30°N 30°N 2

6 2 4 2 6

4 2 6 2 2 4 28°N 4 28°N 2

2 4 2 2 2 2 2 2 4 2

2 2 2

2 2 2 26°N 26°N 2 2

2

2

2

Latitude Latitude 24°N 24°N 2 2

2 2 2

22°N 22°N

2 2

20°N 20°N 80°E 85°E 90°E 95°E 100°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

(c) (d)

Figure 2.1: Spatio temporal variation of the MCSs population over the study region during the pre-monsoon season at (a) from 00-06 MST (b) from 06-12 MST (c) from 12-18 MST and (d) from 18-00 MST.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 39

Further, the spatio-temporal variation of the occurrence of MCSs during 0-6, 6-12, 12- 18 and 18-00 mean standard time (MST) in the monsoon seasons are shown in figure 2.2 (a-d) respectively. Overall it is observed that over the Himalayan foot hill region , Patkai Bum and Meghalaya plateau there is significant diurnal variation whereas over the plain region of Bangladesh and Gangetic west Bengal diurnal variation in the occurrence of MCSs is insignificant. During 00-06 MST, the significant occurrence of MCSs is over the Central Himalaya Foothill region ( around 14 number of MCSs) , Patkai Bum region (around 20 number of MCSs ) and Meghalaya Plateau region (around 10 number of MCSs). Over the plain region of Bangladesh and Gangetic west Bengal the occurrence is significantly less varying in the range of 2-4 MCSs. During 06-12 MST, overall there is relative decrease in the occurrence of MCSs but still significant occurrence over the Patkai Bum region ( around 14 number of MCSs). During 12-18 there is further decrease in the occurrence of MCSs over the Patkai Bum region ( around 8 number of MCSs). Thereafter during 18-00 MST there is again increase in the occurrence of MCSs over Himalaya foothills, Meghalaya Plateau and Patkai Bum region.

MCSs population (Monsoon) ( 00 - 06 MST ) MCSs population (Monsoon) ( 06 - 12 MST ) 30°N 30°N 6 2

4

4 4 2 4 2 2 6 6 2 2 2 6 2 2 4

4 4 4

28°N 8 14 14 2

2 28°N 4 10 8 8 6 4 8 6 2 6 26 2 6 4 2 2 14 6 4 6 4 6 6 4 2 8 2 6 4 8 2 4 6 6 6 6 2 10 10 2 6 6 4 6 14 10 2 8 2 10 2 10 2 2 2 8 2 2 6 4 4 4 6 2 4 8 26°N 4 6 6 6 26°N 4 2 10 6 4 4 4 4 2 8 6 10 2 4 2 2 8 6 6 4 2 10 2 2 4 4 8 10 2 4 2 2 10 6 2 6 2 4 4 4 2 Latitude 4 4 2 2 2 24°N Latitude 4 4 2 2 2 4 4 24°N 2 2 2 6 2 4 6 2 2 4 4 2 4 4 2 6 2 2 2 2 2 2 4 2

4 6 8 22°N 2 2 2 6

4 4 2 22°N 4

6 2 4 2 4 2 2 8 2 2 8 6 4 6 4 2 4 2 2 10 4 4 2 6 20°N 20°N 80°E 85°E 90°E 95°E 100°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude (b) (a)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 40

MCSs population (Monsoon) ( 12 - 18 MST ) MCSs population (Monsoon) ( 18 - 00 MST ) 30°N 30°N 2 2 2 8 8 4 2 4 2 8 2 4 10 6 4 2 2 46 2 6 4 4 6 48 2

2 6 2 2 2 2 6 2 2 4 2 4 2 8 4 2 4 4 2 4 6 28°N 6 28°N 2 14

10 6 4 2 14 2 6 8 6 8 4 4 2 4 8 2 8 2 10 2 8 2 4 2 8 2 2 2 2 4 8 4 8 2 4 8

2 6 4 4 4 4 8 6 6 2 2 4 8 6 4 2 4 6 4 6 2 2 2 4 4 2 2 26°N 6 8 26°N 4 2 2 4 4 2 2 2 2 4 2 2 2 6 4 4 2 4 2 8 6 6 2 4 4 4 2 4

2 4 2 4 4 2 4 2

6 2 2

2 6 2 2 Latitude Latitude 4 2 2 4 4 4 24°N 8 24°N

4 2

2 8 4 6 4 4 4 6 2 2 4

2

2 2 2 4 2 4 6 4 6 2 6 2

4 2 4 2 4 2

4 6 4 4 2 4 2 6 2 2 4 4 2 4 2 4 6 4 2 2 4 2 4 22°N 2 22°N 6

4 4 2 4 6 4 2 2 4 4 2 2 2 2 2 4 4 2 2 4 2 6 2 6 4 4 6 2 2 2 4 2 6 20°N 20°N 80°E 85°E 90°E 95°E 100°E 80°E 85°E 90°E 95°E 100°E Longitude Longitude

(c) (d)

Figure 2.2 : Spatio temporal variation of the MCSs population over the study region during the monsoon season at (a) from 00-06 MST (b) from 06-12 MST (c) from 12-18 MST and (d) from 18-00 MST.

Overall it is observed that during the premonsoon season the occurrence MCSs is predominantly over the eastern Himalaya foothills and Meghalaya Plateau, with the exception of 18-00 MST, when the MCSs occurrence is also observed over the plain region Bangladesh and Gangetic west Bengal. As far as diurnal variability of the MCSs over the study region during monsoon season is concerned it is observed that over the foothill region there is strong diurnal variability with early morning maxima compared to the plain region. With afternoon maxima.

2.2.2 Diurnal variation of occurrence of MCSs over two different climatic regions, A and B during the premonsoon and monsoon seasons

Over the study region two distinct climatic zones are considered in terms of intensity and population characteristics of MCSs, namely the region A (220N- 260 N; 840E- 900E; the Gangetic West Bengal, plain of Bangladesh, Chota Nagpur plateau) and the region B (240N- 280 N; 940E-1000E; the North of Patkai- upper Brahmaputra basin, Pegu

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 41

Yoma- Shan Plateau) as shown in Figure 2.3. The region A is characterized as a region with the most intense MCSs with moderate MCS population. The region B is characterized as a region with moderately intense MCSs with the maximum MCS population ( Choudhury et al. 2015)

30 1 28 3 2 7 4 26 8 9 B

24 A 6 10

Latitude(degree) 5 22 11 20 80 82 84 86 88 90 92 94 96 98 100 Longitude (degree) Figure 2.3: Physical map of the study region (200 N-300 N; 800 E-1000 E)

The temporal or diurnal variation of the frequency occurrence of MCS during the premonsoon and monsoon seasons over region A and B are shown in Figure 2.4 (a, b) respectively. The total number of MCSs in each season is provided in the respective figure panels. There is distinct difference in the diurnal variation over these two regions. Over region A the diurnal variability is stronger compared to region B and also there is different characteristics of the diurnal pattern over these two regions. Over the region A there is maximum occurrence during the late evening to mid night and minimum during morning to mid day. Over the region A there is time lag of maximum occurrence of MCSs between the pre monsoon and monsoon seasons. During the pre monsoon the maximum occurrence is at around 18-23 hours whereas during the monsoon the maximum

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 42 occurrence is at around 16 hrs. Over the region B the diurnal variability is weak and also the maximum occurrence is during 03-06 hrs in contrast to region A.

MCSs MCSs

10 10 9 Premonsoon (112) 9 Premonsoon (379) Monsoon (801) Monsoon (1498) 8 8 7 7 6 6 5 5 4 4 3 Frequency[%] 3 Frequency[%] 2 2

1 1 0 0 0 3 6 9 12 15 18 21 23 0 3 6 9 12 15 18 21 23 Time [MST] Time [MST]

(b) (a)

Figure 2.4 : Diurnal variation of the occurrence of MCSs population during the pre- monsoon and monsoon seasons over the (a) region A and (b) region B region.

Overall it is observed that the diurnal variation of the occurrence of MCSs have a regional variability over the study region. The continental region A has the continental characteristics where the occurrence of MCSs have afternoon maxima whereas continental region B has oceanic character where the occurrence of MCSs have early morning maxima. This is also supported by the fact that the amplitude of the diurnal variation of the occurrence of MCSs over region A is larger compared to the region B.

2.2.3 Diurnal variation of the convective activity as observed by INSAT 3D infrared thermal channel

Further an analysis of thermal imagery data of geostationary satellite INSAT 3D is also analyzed. For this purpose half an hour images of infrared thermal band (48 images/day) are utilized. At present INSAT-3D data are analyzed during the

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 43 premonsoon months (March-May) for three years (2014-2016). It is to be mentioned that , lower the value of the parameter more the convection. The median value of the brightness o o temperature at 12 µm band Tb12 is estimated at 2 x2 grid scale. In the present report spatial variation of the the median value of the Tb12 at different time of the day is presented. The spatial variation of 6 hourly averaged value of Tb12 at 00, 06, 12 and 18 hrs (IST) are shown in figure 2.5 (a-d) respectively. Significant spatiotemporal variation of the Tb12 value is observed. It is observed that during the night (00 hrs) the convective systems are predominantly over the oceanic region (Bay of Bengal) and also over the foothills of central and eastern Himalaya, whereas during the afternoon (18 hrs) convective systems are predominantly over the land region. Overall, deepest convective systems are observed over the Bay of Bengal during the night time. It is interesting to mention that over the land region, the night time maximum occurrence of MCSs over the foothills of Himalaya and after noon maximum occurrence over the Bangadesh plain as observed by the TRMM (Figure 2.1, 2.2 and 2.3 ) is in agreement with the high temporoal resolution observations from the INSAT 3D (Figure 2.5).

(a) (b)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 44

(c) (d)

Figure 2.5: Spatio temporal variation of the median value of Tb12 during the premonsoon season at (a) from 00-06 MST (b) from 06-12 MST (c) from 12-18 MST and (d) from 18-00 MST.

2.3 Summary and Conclusions

A climatological study of diurnal variability of the properties of MCSs is carried out during 1998 to 2012 over the eastern and north eastern region of India along with the adjoining area (200-300N; 800-1000E), a region of complex topography. The MCSs properties are studied with the help of TRMM sensors. Significant spatiotemporal variability of the properties of MCSs is observed. The salient features of the present study are as follows –  The observation from TRMM as well as INSAT 3D suggests that there is a significant regional variability in the diurnal characteristics of the convective systems over the study region. The foothills region of Himalya including the Assam valley have maximum occurrence during night to morning hours, whereas the plain region of Bangladesh, Gangetic west Bengal, northern part of east coast have maximum occurrence during afternoon to night hours.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 45

Chapter 3

Detection of Hail Features (HFs) by satellite onboard microwave sensors over

the north eastern part of India

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 46

Detection of Hail Features (HFs) by satellite onboard microwave sensors over the north eastern part of India

3.1. Introduction

Thunderstorm with hail is a meteorological hazard, which causes significant property damage. Hail is labelled severe, when it is greater than 20 mm in diameter. Diameter greater than 8 mm is used to differentiate hail from graupel. Hail storms are manifestation of severe and deep layer moist convection in the atmosphere. The essential ingredients for the occurrence of severe storms are the presence of an environment of deep layer instability in a warm and humid air mass with high value of Convective Available Potential Energy (CAPE). The general atmospheric conditions under which hail storms occurs are as follows –

(a) High instability with Cb clouds growing to very high levels.

(b) Presence of large vertical currents inside the clouds. Large vertical current is associated with large thermal instability. The hail size is proportional to the updraft velocity in the thunder clouds.

(c) High moisture content in the atmosphere. The greater the moisture content, larger is the size of the hail stones.

(d) Lower freezing level is conducive to hail storms. One of the factors favourable for hail storms in north and central India in winter and early spring is the lower freezing level. Large hails are generally reported with wet bulb zero heights in the range of 2-3 km.

(e) Earlier analysis of upper tropospheric conditions on the day of hail storms show that the occurrence of hails over central and southern peninsula was generally associated with deep upper tropospheric westerly trough and the associated Southern Tropical Jet (STJ) core coming down to lower latitude.

Many researchers have studied the hail phenomena and associated thunderstorms over different parts of India and Bangladesh (Eliot, 1899; Misra and Prasad, 1980; Chaudhaury and Mazumdar, 1983; Ramanamurthy, 1983; Chowdhury and Banerjee, 1983;

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 47

Pandharinath and Bhavanarayana, 1990; Kumar, 1992; Nizamuddin, 1993; Chakravorty and Bhowmik, 1993). Eliot (1899) studied the 597 hail events over India and reported that 153 hail events were with hail stone ≥ 3cm diameter. Chowdhury and Banerjee, (19 3) studied the eight year’s hail data over the north-eastern state and reported 30 hail storms/year Misra and Prasad (1980) pointed out that over the west coast and central inland India (150-200N), the frequency of hail event is relatively low. Chaudhaury and Mazumdar, (1983) pointed out that expectancy % of hail (number of days with hail/no. of days with thunderstorms) decreases from 5 to less than 2 % from March to May in north eastern states and Bangladesh. Pandharinath and Bhavanarayana (1990) carried out a case study of hail storms from 11-13 March 1981 in Andhra Pradesh (16 0N) which inflicted heavy damage to life and property in the state. Nizamuddin (1993) reported around 228 hail days over India during 1982-1989, primarily concentrated over the central high land and Himalayan region. Ramanamurthy (1983) studied the hail frequency for 100 years period and they reported one hail event/year in north-eastern India and 5 to 10 hail events /year in northern India and in the Himalayan region.

Conventionally the above mentioned hail storms studies are based on the ground reporting of the events by the staffers and volunteers. The midlatitude regions boast the existence of dense network of hail pad, whereas tropical region lacks the good network of hail paid or other hail detecting device. Therefore satellite based detection of hail features is an alternative option with promising potential. Satellite-borne passive microwave radiometers record brightness temperature depressions due to the scattering of upwelling radiation by large ice hydrometeors (graupel, hail). Spencer et al. (1983, 1987) and Spencer and Santek (1985) examined satellite measurements of 37-GHz brightness temperature as an indicator of intense convection and severe thunderstorms

In particular, Spencer et al. (1987) matched storm brightness temperatures to U.S. severe weather reports and Spencer and Santek (1985) mapped low–brightness temperature storm events globally. Cecil et al. (2005) and Zipser et al. (2006) have used Tropical Rainfall Measuring Mission (TRMM) measurements of several parameters as proxies to examine intense thunderstorms. These parameters include minimum 85- and 37-GHz brightness temperatures, maximum radar reflectivity at certain altitudes and maximum

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 48 heights attained by certain reflectivity, lightning flash rates, and other measures. More recently Ceil (2009, 2011; Cecil and Blankenship 2012) has utilized the passive microwave brightness temperature as a proxies for the detection of hail storms.

Conventionally the hail storms studies are based on the ground reporting of the events by the staffers and volunteers. The tropical region lacks the good network of hail paid or other hail detecting device. Therefore, satellite based detection of hail features is an alternative option with promising potential. Satellite-borne passive microwave radiometers record brightness temperature depressions due to the scattering of upwelling radiation by large ice hydrometeors (graupel, hail). Spencer et al. (1987) examined satellite measurements of 37-GHz brightness temperature as an indicator of intense convection and severe thunderstorms. Recently Ceil (2009) has utilized the minimum 37-GHz brightness temperatures, as proxies for the detection of hail storms. Motivated by the encouraging results of Cecil (2009), a study is carried out over north-east India to study the hail features in CSs.

The main focus of the present chapter is to study the spatial and monthly variability of the thunderstorms with hail features over the north east part of India by using the 37 (V, H) GHz channels of TRMM Microwave Imager (TMI) and TRMM-PR along with the ground reporting in association with the atmospheric convective parameters. the U.S. severe storm database has been used to quantitatively link TRMM Microwave Imager (TMI) measurements to the occurrence of large hail. Motivated by the encouraging results of Cecil (2009; 2011), a study is carried out over a north-east India to study the hail features in CSs.

3.2 Study Area

The present study is carried out over the north-east part of India within a domain of 23 – 29 0N and 89-100 0E as shown in Figure 3.1. The selected region has a complex topography, where height of the terrain is changing from mean Sea Level (MSL) to mountains of average height of 5000 meters. The ground hail reporting data are considered over the eleven IMD stations over the region. The details of the stations, namely their name, geo-location and MSL height of each station is provided in Table 6.1.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 49

Topographically region consists of Eastern Himalya range, Patkai hill range and there is a valley region, Assam valley along west - east direction.

29°N

28°N PSG MBR

27°N TZP NLP

DHB GHT 26°N SHL Latitude CPJ 25°N SCL IPM

24°N AGT

23°N 89°E 91°E 93°E 95°E 97°E Longitude

Figure 3.1: Physiographic map of the study region Table- 3.1: Geospatial information about the reporting station and reported thunderstorms

Station Name Lat (deg) Long (deg) Height from Total TS MSL (meter) Agartala (AGT) 23.30 91.26 16 189 Cherapanji (CPJ) 25.30 91.70 1313 277 Dhubri (DHB) 26.15 90.13 35 417 Guwahati (GHT) 26.11 91.44 54 1047 Imphal (IMP) 24.82 93.95 781 901 Mohonbari (MBR) 27.47 94.91 111 706

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 50

North Lakhimpur (NLP) 27.24 94.10 102 490 Pasighat (PSG) 28.07 95.33 157 253 Shillong (SHL) 25.56 91.88 1598 411 Silchar (SCL) 24.82 92.80 20 1040 Tezpur (TZP) 26.63 92.80 90 1298

3.3 Data products

In this present study, we collected hail reports from journal records available at the Regional Meteorological Center at Guwahati from 1998 to 2014. These reports are based on 11 stations (Figure 6.1) over North-East India. Surface-based hail reports depend on storm spotter reports. These reports generally have some uncertainty associated with the precise times and locations. Therefore, we identified multiple reports within an hour at a station and it considers a single report, because the hails are coming from same storm.

The University of Utah TRMM level-2 cloud and precipitation feature datasets is utilized (Liu et al., 2008) in this study. This dataset was developed through the collocation of TRMM-Precipitation Radar (TRMM-PR), TRMM microwave Imager (TMI), visible and infrared scanner (VIRS), and lightning imaging sensor (LIS). For the present study, radar detected Precipitation Feature (RPFs) database is utilized. The PFs are defined by contiguous TRMM-PR 2A25 v6 (Iguchi et al., 2000) near-surface raining pixels. The PFs with at least four contiguous pixels are only taken into account to eliminate noise. Six parameters are used to study the characteristics of storms, namely, maximum echo top heights of 20 and 40 dBZ (ETH20dBZ(max); ETH40dBZ(max)), flash counts (FC), minimum polarization corrected brightness temperatures at 85 and 37 GHz (PCT85(min); PCT37(min)) (Spencer et al., 1989) and area of the system. In addition to this parameters, vertical profiles of maximum reflectivity area of the system is defined as number of PR raining pixels multiplied by resolution of TRMM-PR pixel. ETH20dBZ (max) is an indicator of precipitation height, which represent how high the updraft can lift the precipitation size particles, whereas, ETH40dBZ(max) is an indicator of the convective intensity or the updraft speed of the convective cells (Xu and Zipser 2012). PCT85 and PCT37 are scattering signatures of precipitation size ice particles (Cecil, 2009; Nesbitt et al., 2000). Moreover,

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 51 the vertical profile of maximum radar reflectivity is a direct indicator of the vertical structure of the convective core (Xu and Zipser, 2012).

National Centers for Enviromental Prediction (NCEP) reanalysis on a 2.50 x 2.50 grid resolution (Kalnay et al. 1996) are used for this study. We use 4-times daily data to create composite large-scale patterns.

3.4 Methodology

Cecil (2009) provided an approach to quantify the fraction of precipitation features (PFs) with minimum PCT threshold that has corresponding large hail (about 2 cm) in the comparison with the United States storms reports. PFs with minimum PCT37GHz had most effective for identifying hail storms. About 24%, 44%, and 70% of PFs with minimum

PCT37GHz below 255, 220 and 180K have the large hail reports. Minimum 37-GHz PCT depends on vertical profile of hydrometeor contents, types and sizes in the PFs. There are many combinations of them which can be resulted the same brightness temperature over different meteorological regimes (Cecil 2011). So, probability of hail reports of PFs with minimum PCT37GHz over the United States are not directly use to other regimes. Cecil (2011) scaled the minimum 37-GHz PCT values to the equivalent subtropical land (United States) values using a linear best fits coefficients between minimum 37-GHz PCT and mixed phase region ice water content. Here, we scaled the tropical India minimum 37- GHz PCT toward an equivalent subtropical land value using the scale coefficients from Cecil (2011). Using the fit coefficients 257, 210, and 156 K of minimum 37-GHz PCT from the tropical India would be equivalent to 255, 220, and 180 K values from subtropical land.

For the present study, three types of Hail Features (HFs) are defined using its scaled minimum PCT37. They are type one (T-1), type two (T-2), and type three (T-3).

Criteria of minimum PCT37GHz are given in Table 3.2. Three types of HFs are observed at 11 stations within 100 Km from stations. Here, we mainly focus on the characteristics of these storms during the premonsoon season (March-April-May) from 1998 to 2014.

Table- 3.2: Criteria for three different types of HFs

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 52

______

Criteria Probability of hail Type of Total No. No. of reporting for large HFs of Detected Detected HFs hail ( >2cm) HFs with MCSs ______

210- ≤ PCT37GHz(min)< 257 K 24 % T-1 1131 194 156- ≤ PCT37GHz(min)< 210 K 45 % T-2 166 31 PCT37GHz (min) < 156 K 70 % T-3 31 13 ______

In the present study the HFs are further associated with TRMM-PR observations in terms of Echo Top Height (ETH), The detail data base is described by Liu et al. (2009). For the present study data are downloaded from the University of Utah website.

3.5 Results

3.5.1 Spatial variation of hail events

The Table 3.1 also shows the number of reported thunderstorms over each selected station. The significant spatial variability is observed over the study region. The maximum thunderstorms are reported over the Tezpur station (1298) and followed by Guwahati (1047) and Silchar (1040) in Assam. The minimum reported thunderstorms are over Agartala (189) Tripura in and closely followed by Pasighat (253) in Arunachal Pradesh

At the outset, with the help of TMI observation, the occurrence of HFs are identified over each station at different thresholds criteria of PCT37 namely as provided in Table 3.2. The number of detected HPFs for different thresholds over each station are shown in Figure 3.2 (a-c) respectively. The HPFs are considered within the 100 km diameter, with the location of the station at the centre. It is observed that there is change in the characteristics of the spatial variation of the occurrence of HFs for different threshold criteria. The most significant occurrence of type T-1 HFs is found over Guwahati Station (163), followed by Shillong (159) and Cherapungi (143) stations. Apart from them, over Dhubri (129) and Silchar (115) also occurrence is significant. It is also observed that towards the upper part of the Assam valley stations (Tezpur, North Lakhimpur, Mohanbari and Pasighat), the occurrence of type T-1 HFs is relatively less, varying in the range of 56-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 53

91, compared to other stations. For the type -2 HFs, overall occurrence is reduced and the characteristic of spatial variation is changed. Now the maximum occurrence is over the Agartala (33) and Dhubri (33). Over the upper part of the Assam valley, the pattern is similar to type -1 HFs, albeit with significantly reduced occurrence. Further for type -3 HFs, overall occurrence is further reduced and yet the characteristic is same i.e. maximum occurrence is over the Agartala (11) followed by Dhubri (6). Over the upper part of Assam valley, the type-3 HFs are absent. Further the spatial variation of the ground reporting of the hail events at different stations are shown in the figure 3.2 (d). Except over Tezpur station, the pattern of the spatial variation of the ground reporting of hail is similar to the type-3 HF (Figure 3.2 c). In both the cases the maximum occurrence is over Agartala (Type-3 HFs11; Ground reporting: 22).

T-1 T-2 29°N 29°N 91 28°N 28°N 1 56 2 59 1 27°N 82 27°N 2 129 163 33 19 26°N 26°N

159 19 Latitude Latitude 143 28 25°N 115 48 25°N 22 6

24°N 24°N 86 33

23°N 23°N 89°E 91°E 93°E 95°E 97°E 89°E 91°E 93°E 95°E 97°E Longitude Longitude

(a) (b)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 54

T-3 Hail Reports 29°N 29°N

28°N 0 28°N 2 0 6 0 6 27°N 27°N 0 14 6 10 12 26°N 4 26°N 14

2 Latitude Latitude 4 6 25°N 25°N 8 7 3 1 24°N 24°N 11 22

23°N 23°N 89°E 91°E 93°E 95°E 97°E 89°E 91°E 93°E 95°E 97°E Longitude Longitude (c) (d)

Figure 3.2: Spatial variation of occurrence of thunder storm with (a) T-1 HFs (b) T-2 HFs (c) T- 3 HFs. and (d) ground reporting.

Overall it is observed that HFs with low probability of ground reporting (T-1, probability- 24%) occur most frequently in Guwahati with 186 features. In addition, HFs also show more frequent occurrence in Shillong with 180 features, Cherrapunji with 175 features, Dhubri with 168 features and Agartala with 130 features. HFs with high probability of ground reporting (T-3, probability-70%) occur mostly in Agartala with 11 features. Moreover, the occurance is more frequent in Dhubri with 6 features, Guwahati with 4 features, Cherrapunji with 4 features and Shillong with 2 features. These HFs are relatively less at Pasighat, Tezpur, Mohanbari, North Lakhimpur, Imphal. These stations are mainly dominated by HFs with less probability of ground reporting (T-1, probability – 24%). It is important to mention that, the TMI detected HFs are for larger hail size ( > 2 cm). there is possibility that the ground reporting of the hail event with smaller hail size may not be detected by the TMI. Further due to poor temporal resolution, there is large probability to miss the detection of hail event. Imphal has not detected the most severe HFs as well as associated MCSs. This result is in reasonably good agreement with the ground reporting of the hail.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 55

3.5.2 Verical profiles of radar reflectivity during hail events The height profiles of the median value of radar reflectivity factor (dBZ) for T-1, T-2 and T-3 HPFs are shown in Figure 3.3. Convective vertical profiles show that T-3 HFs have strongest vertical profiles and T-1 HFs have weakest vertical profiles throughout all the heights. T-1, T-2, and T-3 HFs have about 35, 45, and 49 dBZ radar reflectivity at 9 Km height (within the mixed-phase region). These differences in vertical structures correlates with hail reporting probability at ground for T-1, T-2, and T-3 HFs, as T-3 HFs have a highest hail detection probability (around 70%) and T-1 HFs have lowest probability of detection, (around 24%). The strongest vertical structure of T-3 HFs indicates strong updraft and large ice particles presence in the mixed - phase region. It also represents that strong mixed-phase microphysical processes (i.e., freezing of raindrops and riming) involved for production of hail/graupel. Overall, for T-1 the gradient is sharper than the T-2 and T-3. The higher value of dBZ at mixed phase region support the defined category of the hail storms, which suggests that T-2 and T-3 are more conducive for hail stone. Overall, the correlation coefficient for PCT37 vs ETH40dBZ are found to be - 0.69.

18 T-1 17 T-2 16 T-3 15 14 13 12 11 10 9

Height(Km) 8 7 6 5 4 3 2 20 25 30 35 40 45 50 55 60 Max reflectivity factor (dBZ)

Figure 3.3 : Vertical profiles of radar reflectivity factor for T-1, T-2 and T- 3 HFs

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 56

The top 50 % (50th percentile) and 1% (99th percentile ) values of Echo top height of 40dBZ each type of PCT37 are provided in the Table 3.3. It is observed that the maximum value of the parameters in both the cases is for T-3 HFs. This is in agreement th th with the top 50 % (50 percentile) and 1% (99 percentile ) values of PCT 37 (Table 3.3) , which have the minimum value for T-3 HFs.

Table 3.3: 50th and 90th percentile values of various parameters for three different types of HPFs

HPFs Parameters 50th percentile 90th percentiles T-1 T-2 T-3 T-1 T-2 T-3 Max-ETH-40 dBZ 7.75 11.50 14.75 10.50 14.50 16.50 (Km)

Min-PCT37GHz (K) 245.73 192.90 149.67 222.87 171.21 126.55 Area (Km2) 295.84 314.33 1164.87 4178.70 5473.00 28104.80

The mean vertical profiles of the radar reflectivity factor over 8 stations are shown in Figure 3.4. A strong spatial variation in the characteristics of vertical profile is observed over these stations. It is observed that over Agartala and and Dhubri the 40 dBZ echoes are observed upto 9.50 and 9.25 km respectively, suggesting a strong mixed phase process (rimming ) over these two stations. The strong mixed phase process is an indicator of greater probability of hail storms over these stations, which is in good agreement with the PCT37 observations. On the other hand over Mohanbari and Pasighat stations the 40 dBZ echoes are observed upto 7.50 and 5.00 km respectively, suggesting a weak mixed phase process over these two stations, which is not a favourable situation for the large hail storms over these two stations. The result over these two stations is also in agreement with the PCT37 observations.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 57

16 Agartala 15 Cherapanji 14 Dhubri 13 Guwahati 12 Mohonbari 11 Pasighat Shillong 10 Silchar 9

8 29°N Height(Km) 7 28°N PSG MBR 6 27°N TZP NLP 5 DHB GHT 26°N SHL 4 Latitude CPJ 25°N SCL 3 IPM

24°N AGT 2 20 25 30 35 40 45 50 55 23°N 89°E 91°E 93°E 95°E 97°E Max reflectivity factor (dBZ) Longitude

Figure 3.4: Mean vertical reflectivity profile over eight stations

3.5.3 Monthly and diurnal variation of hail events During the pre-monsonn, the monthly variation of the occurrence of satellite detected T-1, T -2 and T-3 type of HFs along with the reported hail events are shown in Figure 3.5 (a, b) respectively. It is observed that for all types of storms the maximum occurrence is during the month of April albeit with different occurrence for each case. The similar pattern is observed for reported hail events also. These results are consistent with the previously reported work on hail storm occurrence. Chaudhary and Mazumdar (1983) shows that the expectancy percentage of the hail derease from 5% to less than 2% from

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 58

March to May for north eastern India and Bangladesh. Chowdhuiry and Banerjee (1983) examined eighth year years data for north eastern India. They found that 30 hail storms occurs per annum. The seasonal variation shows the influence of westerly trough in the pre-monsoon seasons with maximum in April. They further find a jet stream maximum, upper lecvel divergence accompanying the trough. They also showed the late evening hail storms maximum. They further pointed out that shear in the 850-300 mbar level is four time greater than average for the pre-monsoon months. Ramamurthy (1983) studied the frequency of hail for a 100 years period. They showed there is one event per year in north - east India. Further during the pre-monsoon season the diurnal variation (at the increment of 03 hrs) of the occurrence of T-1, T-2 and T-3 type of HFs along with the reported hail storms are shown in Figure 3.6 (a, b) respectively. It is observed that the T-2 and T-3 HFs occurs predominantly during 2100-2400 hrs with minimum occurrence during 0600-0900 hrs. The reported hail storms are during 1800 to 2100 with minimum occurrence at 2400 to 0300 hrs.

600 Hail reports 500 483 T-1 390 400 T-2 400 T-3 300 258

Storms 200

Hailreports 200 97 100 65 57 29 21 21 4 1 9 0 0 March April May March April May

Figure 3.5: Monthly variation of (a) reported hail storm events (b) Satellite detected thunderstorm for Type1, Type 2 and Type 3 HPFs

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 59

50 45 50 45 40 40 35 35 30 30 25 25

20 20 % HPFs % 15 15 % Hail % reports 10 10 5 5 0 0 0 2 4 6 8 10 12 14 16 18 2021 0 3 6 9 12 15 18 21 IST (hrs) Local time (IST)

Figure 3.6: Diurnal variation of (a) reported hail storm events (b) Satellite detected thunderstorm for Type1, Type 2 and Type 3 HPFs

3.5.4 Interannual variation of occurrence of detected and reported hail events The inter-annual variation of the occurrence of T-1, T-2, T-3 HFs and reported hail storms are shown in Figure 3.7 (a-d) respectively. The peaks ( 1999, 2004, 2010, and low (2007,2011, 2013 ) of the satellite detected and ground reporting are matching

30 120 coefficient (-0.70) coefficient (-0.44) 110 25 100 20 90 15 80

T-2 Storms T-2 10 T-1 Stroms T-1 70 60 5 50 0 1998 2000 2002 2004 2006 2008 2010 2012 2014 1998 2000 2002 2004 2006 2008 2010 2012 2014 Years Years (a) (b)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 60

15 15 14 Coefficient ( - 0.14 ) 14 Coefficient (-0.10) 12 12 10 10 8 8

6 6 T-3 Storms T-3 4 Hailreports 4 2 2 0 0 1998 2000 2002 2004 2006 2008 2010 2012 2014 1998 2000 2002 2004 2006 2008 2010 2012 2014 Years Years (c) (d)

Figure 3.7 : Trend of annual variation of Satellite detected thunderstorm for (a) T-1, (b) T-2 (c) T-3 HPFs and (d ) reported hail storm events.

3.5.5 Case study of the hail events

Three case studies of the reported hail events over the Agartala stations on May 28th 1998 , April 7th , 2003 and April 10th , 2010 are presented. On these three days near simultaneous collocated observation by TRMM sensors are available. The simultaneous observations of PCT37, near surface radar reflectivity and longitude/latitude vs height cross section of reflectivity during 1998, 2003 and 2010 are provided in figure 3.8 (a-c), 3.10 (d-f), 3.10 (g-i) respectively. The hail storms in 1998 and 2010 were multi cellular bow echoes events and a hail storms in 2003 was a single cell event. It is observed that during all these three events, the value of PCT37 and surface reflectivity are in the range of 120 -140 K and 50-55 dBZ respectively. It is also observed that in 1998, 2003 and 2010 the 40 dBZ echoes were observed upto 16 km, 13 km and 10 km respectively, indicating the presence of strong mixed phase process during the hail events.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 61

Agartala: May 28th , 1998 Agartala: April 7th, 2003 Agartala: April 10th, 2010

Figure 3.8 (a-c) Figure 3.8 (d-f) Figure 3.8 (g-i)

2870 1998 - 05 - 28 ( 20:99 UTC) 30906 2003 - 04 - 17 (21:45 UTC) K K 82036 2012 - 04 - 10 ( 8:37 UTC) K 300 300 26°N 26°N 26°N 300 PCT 37 GHz (K) PCT 37 GHz (K) PCT 37 GHz (K) 280 280 Min 37 PCT(K) - 147.63 K Min 37 PCT(K) - 151.39 K 280 Max 40 dBZ - 13.50 Km Max 40 dBZ - 15.25 Km 260 260 260 25°N Min 37 PCT(K) - 111.40 K 25°N 25°N Max 40 dBZ - 16.75 Km 240 240 240

220 220 220 24°N 24°N 24°N 200 200

Latitude 200

Latitude Latitude 180 180 180 23°N 23°N 160 160 23°N 160 140 140 140 120 22°N 22°N 120 22°N 120 88°E 90°E 92°E 94°E 88°E 90°E 92°E 94°E 88°E 90°E 92°E 94°E Longitude Longitude Longitude (a) (d) (g)

30906 2003 - 04 - 17 (21:45 UTC) dBZ 82036 2012 - 04 - 10 ( 8:37 UTC) dBZ 2870 1998 - 05 - 28 ( 20:99 UTC) dBZ 26°N 55 26°N 55 26°N 55 NSZ (dBZ) NSZ (dBZ) NSZ (dBZ) 50 50 50 25°N 25°N 45 45 25°N 45 40 40 B B 40 24°N 24°N

B Latitude 24°N 35 Latitude 35

Latitude 35 A A A 30 30 23°N 23°N 30 23°N 25 25 25 22°N 20 22°N 20 88°E 90°E 92°E 94°E 88°E 90°E 92°E 94°E 22°N 20 Longitude Longitude 88°E 90°E 92°E 94°E (b) Longitude (e) (h)

2870 1998-05-28 20.99 UTC 82036 2012-04-10 08.37 UTC dBZ 30906 2003-04-17 21.45 UTC dBZ 20 55 dBZ 20 55 20 55 18 18 50 18 50 50 16 16 16 45 45 45 14 14 14 12 40 12 40 12 40 10 35 10 35 10 35 Height(Km) 8 Height(Km) 30 8

Height(Km) 30 8 6 30 6 4 25 6 4 25 25 2 20 4 23.35 23.4 23.45 23.5 23.55 23.6 23.65 2 20 A Latitude (degree) B 91.5 91.7 91.9 92.1 92.3 92.5 2 20 A Longitude (degree) B 91.6 91.8 92 92.2 92.4 92.6 92.8 A B (c) Longitude (degree) (f) (i)

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 62

3.6 Summary and Conclusion The present work is carried out to study the spatial and temporal variability of the occurrence of HFs over the north east India. The analysis is carried out over the 11 stations. The results are as follows

 On the basis of threshold value of the Polarization corrected temperature of 37

GHz channels (PCT37), hail features (HFs) are classified into three categories namely, T-1 (with hail detection probability of 24%) , T-2 (with hail detection probability of 45% ) and T-3 (with hail detection probability of 70%).  T-1, T-2, and T-3 HFs have about 35, 45, and 49 dBZ radar reflectivity at 9 Km height (within the mixed-phase region). The nature of vertical structures of reflectivity for T-1, T-2, and T-3 HFs. correlates well with hail reporting at ground. The strongest vertical structure of T-3 HFs indicates the strong updraft and large ice particles presence in the mixed - phase region. It also represents that the strong mixed-phase microphysical processes (i.e., freezing of raindrops and riming) are involved for production of hail/graupel.  The two stations in the plain region, namely Agartala and Dhubri, detected the maximum occurrence of most severe HFs (type T-3). The stations in the valley regions namely Tezpur, Mohanbari, North Lakhimpur, Pasighat, and Imphal have not detected the most severe HFs. The spatial variability in the HFs is amply supported by the vertical profiles of reflectivity and its value at mixed phase region.  During the premonsoon, the maximum occurrence of HFs is found in April. The occurrence of HFs is minimum in March.  HFs show strong diurnal variation. with maximum occurrence during the afternoon hours  Occurrence of HFs show noticeable year to year variation. There is a decreasing trend during the period 1998 - 2013. The trend values are -0.70, -0.44, -0.14 and - 0.10 for T-1, T-2 and T-3 HFs and ground reporting respectively.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 63

Climatologically, there is a spatial and temporal variability in the occurrence of HPFs over the study region. With preferential occurrence over the plain regions. The patterns of seasonal and diurnal occurrence of MCS with hail features as detected by satellite derived parameters are consistent with the ground report over the stations.

Compared to other parameters, the ETH 40dBZ is better correlated with PCT37 with a CC values of -0.69 respectively. It is suggested that these values can also be considered as a proxies for the detection of hail features. The present work is still in a developing stage. Regional or regime dependent variability of hail detection can lead to similar brightness temperatures coming from different hydrometeor profiles. This methodology provides an approach for objective climatologies that do not rely on surface reports or spotter networks, which vary greatly from region to region”.

-000-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 64

Future Plan

(i) Installation of Instruments at different locations Kohima Science College: Micro Rain Radar Parsivel Disdrometer Lightning Detector ( LD -100) Electric Field Mill (EFM-350) Automatic Rain gauge. Indian Statistical Institute, Giridih: Parsivel Disdrometer (Jharkhand) Lightning Detector ( LD -100) Electric Field Mill (EFM-350) Automatic Rain gauge Fazal Ali College, Mokokchung (Nagaland): Automatic Rain gauge Phek College, Phek (Nagaland) Automatic Rain gauge Sao Chang College, Tuensang (Nagaland) Automatic Rain gauge , Dimapur (Nagaland) Automatic Rain gauge Zunheboto College, Zunheboto (Nagaland) Automatic Rain gauge Mounttiyi College, Wokha (Nagaland) Automatic Rain gauge Mon College, Mon, (Nagaland) Automatic Rain gauge Zisaji College, Kiphre (Nagaland) Automatic Rain gauge (ii). Analysis of Socio-Psycho data Data over west Bengal and Nagaland and Assam (iii) Downloading and Analysis of satellite data INSAT 3D, GPM and CloudSat

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 65

Reference

Alcala CM, Dessler AE. 2002. Observation of deep convection in the tropics using the tropical rainfall Measuring Mission precipitation Radar. J. Geophys. Res. 107: D24, 4792, doi:10.1029/2002JD002457. Bhatt BC, Nakamura K. 2005. Characteristics of monsoon rainfall around the Himalayas revealed by TRMM precipitation radar. Mon. Wea. Rev. 133:149-165. Bhat GS, Kumar S. 2015. Vertical structure of cumulonimbus towers and intense convective clouds over the south Asian region during the summer monsoon seasons. J. Geophys. Res. Atmos. 120: Doi10.1002/2014JD022552. Boccippio DJ, Peterson WA, Cecil DJ. 2005. The tropical convective spectrum. Part I: Archetypal vertical structures. J. Clim. 18:2744-2769. Cecil DJ, Goodman SJ, Boccippio DJ, Zipser EJ, Nesbitt SW.2005. Three years of TRMM precipitation features. Part I: Radar, radiometric and lightening characteristics. Mon. Wea. Rev. 133: 543-566. Cecil DJ. 2009. Passive microwave Brightness temperature as proxies for hail storms. J. Appl. Meteorol. Climatol.48:1281-1286. Cecil DJ, Blackenship CB. 2012. Toward a global climateology of severe hail storms as estimated by satellite passive microwave imager. J. Climate. 25:687-703. Chakabarty KK, Bhowmik SK.1993.Study of unusual hailstorm over Bombay, Mausam, 44: 292-295. Chaudhary A , Mazumdar AB. 19 3. Nor’westers and the synoptic climatology of hot weather season in northeast India . Vayu Mandal. 13: 48-55. Chowdhury A, Banerjee AK. 1983. A study of hailstorms over northeast India. Vayu Mandal. 13: 91-95. Choudhury H, Roy P, Kalita S, Sharma S .2015. Spatio-temporal variability of the properties of mesoscale convective systems over a complex terrain as observed by TRMM sensors. Int. J. Climatol., DOI: 10.1002/joc.4516. Dai A, Giorgi F, Trenberth KE, 1999. Observed and model simulated diurnal cycle of precipitation over the contiguous United states. J. Geophys. Res. 104: 6377-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 66

6402. Dai A, 2001.Global precipitation and thunderstorm frequencies. Part –II: Diurnal variation. J. Clim. 14: 1112-1128. Dye JE, Winn WP, Jones JJ, Breed DW.1989. The electrification of new Mexico thunderstorms 1. Relation between precipitation development and the onset of electrification. J. Geophys Res. 94:8643-8656. Eliot J. 1899.Hailstorms in India during the period 1883- 1897 with a discussion of their distribution. India Meteor, Memoirs 6(4): 237-315. Gadgil S. 2003. The Indian monsoon and its variability, Annu. Rev. Earth Planet. Sci. 31: 429-467. Gray WM, Jacobson RW, 1977. Diurnal variation of deep cumulus convection. Mon. Wea. Rev.105: 1171-1188. Hamada A, Murayama Y, Takayabu YN. 2014. Regional characteristics of extreme rainfall extracted from TRMM PR measurements. J. Clim. 27: 8151-8169. Hamada A, Takayabu YN, Liu C, Zipser EJ. 20115. Weak linkage between the heaviest rainfall and tallest storms. Nature Communications doi 10.1038/ncomms7213. Hirose, M., and K. Nakamura, 2005, Spatial and diurnal variation of precipitation systems over Asia observed by the TRMM Precipitation Radar, J. Geophys. Res. 110, D05106. Houze RA, Wilton DC, Smull BF. 2007. Monsoon convection in the Himalayan region as seen by the TRMM Precipitation Radar, Q. J. Roy. Meteorol. Soc. 133: 389- 411. Houze RA, Rasmussen KL, Zuluaga MD, Brodzik SR.2015. The variable nature of convection in the tropics and subtropics: a legacy of 16 years of the tropical rainfall measuring mission satellite. Rev. Geophys. 53: doi: 10.1002/2015RG000488. Islam MN, Uyeda Hiroshi. 2008.Vertical variation of rain intensity in different rainy periods in and around Bangladesh derived from TRMM observations. Int. J. Climatol. 28 : 273-279. Kumar A, 1992. A climatological study of thunderstorms at Lucknow airport. Mausam, 43:

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 67

441-444. Liu C, Zipser EJ. 2005. Global distribution of convection penetrating the tropical tropopause. J. Geophys. Res. 110: D23104,doi:10.1029/2005JD006063. Liu C, Zipser EJ, Cecil DJ, Nesbitt SW, Sherwood S, 2008. A Cloud and Precipitation Feature Database from Nine Years of TRMM Observations. J. Appl. Meteor. Climatol., 47: 2712–2728. Liu C. 2011. Rainfall contributions from precipitation systems with different sizes, convective intensities and durations over the tropics and subtropics. J. Hydrometeorol. 12: 394-412. Liu C, Cecil D, Zipser EJ. 2011. Relation between lightening flash rate and passive microwave brightness temperatures at 85 and 37 GHz over the tropics and subtropics. J. Geophys. Res. 116: doi 10.1029/2011JD016463. Liu C, Cecil D, Zipser EJ, Kronfeld K, Robertson R. 2012. Relation between lightning flash rate and radar reflectivity vertical structure in thunderstorms over the tropics and subtropics. J. Geophys. Res. 117: doi 10.1029/2011JD017123. Luo Y, Zhang R, Qian W, Luo Z, Hu X. 2011. Inter-comparison of deep convection over the Tibetan plateau-Asian monsoon region and subtropical North America in Boreal summer using CloudSat/CALIPSO data, J. Clim. 24: 2164-2177. Medina S, Houze RA, Kumar A, Niyogi D. 2010. Summer monsoon convection in the Himalaya region: terrain and land cover effects. Q. J. R. Meteorol. Soc. 136: 593- 616. Misra PK, Prasad SK. 1980. Forecasting hailstorms over India. Mausam. 31: 385-396. Nesbitt SW, Zipser EJ, Cecil DJ. 2000. A census of precipitation features in the tropics using TRMM: radar, ice scattering, and lightning observations. J. Clim, 13: 4087– 4106. Nesbitt, SW, Zipser EJ. 2003. The diurnal cycle of rainfall and convective intensity according to three years of TRMM measurements. J. Clim. 16: 1456-1475. Nesbitt SW, Cifelli R, Rutledge SA. 2006. Storm morphology and rainfall characteristics of TRMM precipitation features. Mon. Wea. Rev. 134: 2702–2721. Nesbitt SW, Anders AM. 2009. Very high resolution precipitation climatologies from

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 68

the tropical rainfall measuring precipitation radar. Geophy. Res. Lett. 36:L15815 doi: 10.1029/2009GL038026. Nizamuddin S.1993. Hail occurrences in India. Weather. 48: 90-92. Oki T, Musiake K, 1994. Seasonal change of the diurnal cycle of precipitation over Japan and Malaysia. J. Appl. Meteor. 33:1445-1463. Pai, D. S., L. Sridhar, and R. Kumar, 2015, Active and break events of Indian summer monsoon during 1901–2014, Clim. Dyn., doi:10.1007/s00382-015-2813-9. Pandharinath N, Bhavanarayana V.1990. Hailstorm over Telangana. Mausam, 41: 433-4 Petersen WA, Christian HJ, Rutledge SA. 2005. TRMM observations of the global relation between ice water content and lightening. Geophys. Res. Lett. 32: L14819, doi:101029/2005GL023216. Qie X, Wu X, Yuan T, Bian J, Lu D. 2014. Comprehensive pattern of deep convective systems over the Tibetan-plateau –south Asia Monsoon region based on TRMM data. J. Clim. 27: 6612-6626. Ramanamurthy BV.1983. Some cloud physical aspects of local severe storms. Vayu Mandal.13: 3-11. Romatschke U, Medina S, Houze JR. 2010. Regional, seasonal and diurnal variation of extreme convection in the south Asian region. J. Clim. 23: 419-439. Romatschke U, Houze Jr R. 2011a. Characteristics of precipitating convective systems in the premonsoon season of south Asia. J. Hydrometeorol. 12: 157-180. Romatschke U, Houze Jr R. 2011b. Characteristics of precipitating convective systems in the south Asian monsoon. J. Hydrometeorol. 12: 3-26. Roy P, Biswasharma R, Deshamukhya A, Sharma S, Gairola R M. 2017. A study of the spatio-temporal variability of the properties of intense precipitation features over South Asian region: An integrated multi sensor approach. Intl. J. Climatol. doi: 10.1002/joc.5027 (In Press). Saunders CPR, Keith WD, Mitzeva RP.1991. The effect of liquid water on thunderstorm charging. Journal of Geophysical Research 96:11007-11017. Simpson J, Adler RF, North GR. 1988. A proposed tropical rainfall measuring mission (TRMM) satellite. Bull. Amer. Meteor. Soc. 69: 278-295.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 69

Spencer RW, Santek DA.1985. Measuring the global distribution of intense convection over land with passive microwave radiometry. J. Climate Appl. Meteor. 24: 860– 864. Spencer RW, Olson WS, Rongzhang W, Martin DW, Weinman JA, Santek DA. 1983. Heavy thunderstorms observed over land by the Nimbus 7 scanning multichannel microwave radiometer. J. Appl. Meteor. 22: 1041–1046. Spencer RW, Howland MR, Santek DA. 1987. Severe storm identification with satellite microwave radiometry: An initial investigation with Nimbus-7 SMMR data. J. Climate Appl. Meteor. 26: 749–754. Ushio T, Heckman SJ, Boccippio DJ, Christian HJ. 2001. A survey of thundersotorm flash rates compared to cloud height using TRMM satellite data. J. Geophys. Res. 106: 24089-24095. Wiens KC, Rutledge SA, Tessendorf SA. 2005. The 29 June 2000 supercell observed during STEP. Part II: Lightening and charge structure . J. Atmos. Sci. 62: 4151- 4177. Xu W, Zipser EJ. 2012. Properties of deep convection in tropical continental, monsoon and oceanic rainfall regimes. Geophy. Res. Lett. 39: L07802, doi: 10.1029/2012GL.051242. Xu, W., 2012, Precipitation and convective characteristics of summer deep convection over East Asia observed by TRMM, Mon. Wea Rev. doi:10.1175/MWR-D-12- 00177.1. Yang GY, Slingo J, 2001, The diurnal cycle in the Tropics. Mon. Wea. Rev. 129: 784-801. Zipser EJ, Lutz KR. 1994. The vertical profile of radar reflectivity of convective cells: A strong indicator of storm intensity and lightening probability. Mon. Wea. Rev. 122:1751-1759. Zipser EJ, Cecil DJ, Liu C, Nesbitt SW, Yorty, DP. 2006. Where are the most intense thunderstorms on earth? Bull. Am. Meteorol. Soc. 87: 1057–1071.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 70

Annexure I: Associate members of RUSA project (RI &QI), Nagaland

Sl. Name Designation Address No.

1. . Prof. Srimanta Pal Associate Electronics and Member Communication Science Unit (ECSU), Indian Statistical Institute, Kolkata 2. Kekhriele Nakhro Associate Assistant Professor Member Department of Zunheboto Government College Zunheboto, Nagaland

3 R. Bendangtemjen Associate Assistant Professor Member Department of Geography Mokokchung, Nagaland

4 Nungsungtula Associate Assistant professor Member Department of Physics, Saochang Government College, Tuensang, Nagaland

5 Forchiba Kichu Associate Assistant professor Member Department of Physics, Phek Government College, Phek, Nagaland 6 Utpal Misra Associate Head, Member Department of physics, Patkai Christian College, Dimapur

7. Meripeni Ezung Associate Assistant Professor Member Department of Physics Kohima Science College, Jotsoma Kohima, Nagaland

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 71

8. Selie Puro Associate Assistant Professor Member Department of Geography Kohima Science College, Jotsoma Kohima, Nagaland 9. Meniele K.Nuh Associate Assistant Professor Member Department of Geology, Kohima Science College, Jotsoma Kohima, Nagaland

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 72

Annexure II: Expert Committee for RUSA project (RI &QI), Nagaland

Sl. Name Designation Address No.

1. . Prof. Animesh Maitra Professor, Institute of Radio Physics and Electronics, University of Calcutta 92 Acharya Prafulla Chandra Road Kolkata 700009 2. Shri Hargobinda Pathak Member Deputy Director General of Meteorology (Retd) Regional Meteorological Centre, Guwahati Indian Meteorological Department, Guwahati-781017 3. Prof. Glen Thong Member Department of Geology Member Nagaland University, Kohima Kohima 797004 4. Prof. M.S. Rawat Member Head, Department of Member Geography Nagaland University, Lumami Lumami-798601

5. Dr. Anungla Aier Member Principal Member Kohima Science College (Autonomous), Jotsoma Kohima Nagaland-797002 6. Dr. Sanjay Sharma Member Assistant Prof. Member Secretary Secretary ( PI –RUSA R & I Project) Department of Physics, Kohima Science College , Jotsoma, Kohima Nagaland- 797002 7. State Project Director RUSA, Special Department of Higher Education Nagaland Invitee Kohima, Nagaland-797004

8. Nominee Nagaland State Special Department of Home, Disaster Management Authority Invitee Nagaland. (NSDMA) Kohima, Nagaland-797004

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 73

Annexure III A: Publication

Sl. Authors Title Reference Impact factor No (Thomson Reuter)

1. Roy P, Biswasharma A study of the International 3.61 A Journal of R, Deshamukhya A, spatio-temporal Journal of Royal Sharma S, Gairola R variability of the Climatology, Meteorological M properties of intense doi: Society, UK precipitation features 10.1002/joc. (Published by Wiley over South Asian 5027. International ) region: An (2017) integrated multi sensor approach. In Press

Annexure III B: Paper presented in the conferences

Sl. Author Title Place and date No

1 Sharma S Study of Hail Features over India Conference on the North East India by using 37 Radar Meteorology, GHz TRMM Microwave Imager IIT Kharagpur, Channels and TRMM January 8 - 11, 2017 Precipitation Radar

2 Roy P Seasonal and Intraseasonal India Conference on variability of rainfall Radar Meteorology, characteristics of convective IIT Kharagpur, systems over the South Asian January 8 - 11, 2017 region using TRMM-PR

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 74

3. Biswasharma R Properties of Deep convective India Conference on clouds over the northeastern part Radar Meteorology, of India and adjoining region as IIT Kharagpur, observed by CloudSat onboard January 8 - 11, 2017 Cloud Profiling Radar

4. Roy P A study of the spatio-temporal North East Space variability of the properties of Application Centre intense precipitation (NESAC) -Academia features over the South Asian and Students region: An integrated multi sensor Interaction meet, approach, Umiam, Shillong,, June 24th , 2016.

5. Biswasharma R A study of convective systems North East Space over over north eastern India by Application Centre using INSAT – 3D observations (NESAC) -Academia and Students Interaction meet, Umiam, Shillong, June 24th , 2016.

Annexure III C: Participation to short term courses/training programs

Sl. Name Name of the Resource Persons Place and No Training date

program

1. Sharma S Forecasting Dr. C A Doswell ESSL Science severe and Training convection I USA Wiener Neustadt, Austria June 20th -24th , 2016.

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 75

2. (i). Roy P Short term (i). Prof. V IIT, Khragpur,

course on Radar Chandrasekar Colarado th (ii) Bishawasharma R Meteorology State University, USA January 8 , 2017 (iii) Imolemba (ii), Dr, N. Bharadwaj

(iv) Sharma S North Western University, USA

(iii) Prof. U Ushio, Osaka University, Japan,

(iv)

3. Project Assistant Short term (i). Prof. Martin Haigh, Nagaland course on Oxford Brook University, (i) Nyuthe V “Landslide and University, UK. March 7th - Associate Members debris flow th systems: (ii) Prof. G. T. Thong 11 , 2017 (ii). Meniele K.Nuh Prediction, Nagaland University, control and (iii).Selie Puro reclamation” (iii) Prof M S Rawat

(iv) Kekhriele Nakhro Nagaland University (v) R. Bendangtemjen [Under Global (vi) Forchiba Kichu Initiatives of Academic Network (GIAN) initiative, MHRD]

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 76

Annexure IV: Severe Weather Survey North Eastern India 2016 © A.G.Keul & S.Sharma final version 29-9-16 The following questions deal with your daily weather information and with severe weather. Thank you very much for your interest to take part in this anonymous survey ! 1. Are you weather-exposed in your daily life?

 yes, always  yes, often  sometimes  no, not at all 2. Do you follow the daily weather report(s) in the media?

 yes, always  yes, often  sometimes  no, not at all 3. Where do you get your daily weather information? [please mark the most important sources]  newspaper  internet  cellphone  radio  television  family  other:………  other:………  other:……… 4. Is the quality of the weather report in radio/TV OK for you?  yes  not always  no 5. Do you easily get your local weather from media weather reports?  yes  not always  no 6. Should media weather reports give precautions?  yes  only in severe danger  no 7. How dangerous do you think are any of the following severe weather phenomena for you? [Give each of them a number between 10=very dangerous and 0=not dangerous at all] ….. cyclones ….. heat/drought ….. landslides ….. hail …..tornadoes ….. floods ….. lightning ….. heavy rainfalls …..wildfires …..thunderstorms 8. Are you personally afraid about severe weather?  yes, very  at times  no, never 9. Do you feel well-informed about severe weather?  yes, very  at times  no, never 10. Did you learn about severe weather protection in your school education?  yes  no 11. Areas called „high“/“low“ on a weather map:

 are hotter/colder  are more windy/calm  have different air pressure  I don’t know

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 77

12. Do you know different cloud types?

 yes  no

Give an example: ………………………………….…… 13. Can you find out by your own observation whether severe weather (excessive rain, hail, lightning, storm) is upcoming?  yes, certainly  not always  no  I don’t know I4. When three seconds pass between a lightning flash and its clap of thunder, this lightning flash had a distance of ………. kilometer(s) / ……….. miles from my own position.  I don’t know I5. Who is hit by lightning, is killed instantly  yes  no  don’t know I6. You are safer from lightning inside a building.  yes  no  I don’t know I7. Out in the open, lie flat on the ground in a thunderstorm.  yes  no  I don’t know I8. Re-animation/resuscitation can help people hit by lightning.  yes  no  I don’t know 19. In which seasons storms are most likely in Northeastern India? Mark the correct season(s):  December-February  March-May  June-August  September- November  I don’t know 20 In which season lightning is most likely in Northeastern India? Mark the correct months:  January  February  March  April  May  June  July  August  September  I don’t know 21 In which season usually monsoon rainfalls occur in Northeastern India? Mark the correct months:  April  May  June  July  August  September  October  November  December  I don’t know 22 In which season(s) hail is most likely in Northeastern India? Mark the correct season(s):  December-February  March-May  June-August  September- November  I don’t know 23 In Northeastern India, falling hailstones: [Mark correct answers]

 do no damage at all  can damage crops  can damage cars/houses

 can injure people/animals  can kill people/animals  I don’t know 24 In Northeastern India, hailstones can reach a size of ….. cm [/……inches?]. 25 In India, landslides are mostly caused by:  lightning  heavy rainfall  earthquakes 26 When a big landslide occurs in India, usually:

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 78

 you have enough time to escape  you have no time to escape  I don’t know 27. Inside a house, you are safe from a landslide:  yes  no  I don’t know 28. Do you see a change in the occurrence of severe weather in the last 15-20 years?

 yes, increasing  yes, decreasing  no change  I don’t know 29. Was your house/property ever hit by lightning?

 yes, severe damage  yes, moderate damage  no 30. Was your house/property/crops ever hit by a flood?  yes, severe damage  yes, moderate damage  no 31. Was your house/property/crops ever hit by a storm?

 yes, severe damage  yes, moderate damage  no 32. Was your house/property/crops ever hit by hail?

 yes, severe damage  yes, moderate damage  no 33. Was your house/property/crops ever hit by heavy rain?

 yes, severe damage  yes, moderate damage  no 34. Was your house/property/crops ever hit by a landslide?

 yes, severe damage  yes, moderate damage  no 35. Do you feel personally prepared for possible risks of severe weather?  yes  partly  no 36. Do you hold an insurance for risks of severe weather?  yes, for most  for some  no 37. Did you receive government support after a severe weather calamity?

 yes, full  yes, partly  none 38. Should severe weather protection be a part of the school education curriculum?

 yes  no  I don’t know 39. Do you know the homepage of the National Disaster Management Authority of India?

[http://www.ndma.gov.in]  yes  no 40. Did you read safety tips on disaster protection by the National Disaster Management Authority of India?  yes  no

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 79

41. Should there be a voluntary rescue training available for all people?  yes  no

I am:  female /  male My age: …... years My profession: ………………..……… In my household, I live together with ……… adults and ………. children/youngsters. My home town/city: ……………………....…….. My federal province: ………………………..…… I live in a:  single house/hut  double or row house  multistorey house  high-rise My highest education level: [only mark the highest]

 elementary/primary school  apprenticeship  middle/secondary school

 higher secondary/board exam  college/university/vocational (Bachelor, Master, PhD) In my home town/city:  I do professional rescue work  I do voluntary rescue work  I was emergency-trained by a helping organization Many thanks for your valuable help! Please return the filled-in questionnaire directly to the person who gave it to you. If you are interested in our results, you can leave us your e-mail address.

-000-

Kohima Science College, Jotsoma, Kohima, Nagaland, INDIA 80