Lightning and precipitation have received a lot of attention in the scientific literature, not only because it is an impressive atmospheric phenomenon but also because it is associated with severe storms that may cause important damages to agriculture, electric power networks, property and life      . Stability indices are a measure of the atmospheric static stability. Their values are used to quickly assess the potential of the atmosphere to produce convection and consequently severe weather. LI indicates the likelihood of severe thunderstorms. Lifted index (LI) is a measure of upper level instability and takes into account elevated convection which tends to occur when upper level disturbances move across unstable equilibrium environment aloft. Negative LI indicates the possibility of convection. Values less than −9˚C reflect extreme instability. A value of LI between 0˚C and −2˚C indicates a small chance of having a severe thunderstorm. Air mass thunderstorms can occur even when the LI is slightly positive    . K Index has been proved useful in indicating the probability of air mass thunderstorms. As the K Index increases, so does the probability of having an air mass thunderstorm. K Index takes into account three components which are important for convective instability: temperature lapse rate, lower tropospheric humidity and the vertical extent of the moist layer. Values of K Index lower than 15˚C  indicate very small likelihood of thunderstorms, while K Index values higher than 35˚C indicate high probability for the development of numerous and/ or severe thunderstorms.
Large body of literature is available for thunderstorm and lightning research  -  . Also, over the Indian region, attempts have been made to understand relation of thunderstorms to wet bulb temperature     , lightning- rainfall relationship  and also, thunderstorm climatology  . The analysis of satellite data in this study focuses on the association of rainfall and stability index with lightning parameter over the hot spot region of Indo-Gangetic Plains.
2. Data and Source Description
The Lightning Imaging Sensor (LIS) is a satellite-borne instrument used to detect the lightning flashes occurring in the tropics   . It is a National Aeronautics and Space Administration’s (NASA) Earth Observing System (EOS) instrument on the Tropical Rainfall Measuring Mission (TRMM) Satellite. The LIS is designed with a higher sensitivity and spatial accuracy than the Optical Transient Detector (OTD). It detects lightning flash with storm-scale (spatial) resolution (4 to 7 km) over a large region (600 × 600 km) of the Earth’s surface. LIS was launched on 28 November 1997 into low Earth orbit to circle the Earth at an altitude of 350 km. The TRMM satellite was boosted from an average altitude of about 350 km before August 2001 to about 400 km after August 2001. The higher altitude (at an orbital inclination of 35˚) gives a correspondingly larger field of view for each of the sensor pixels and swath and longer sampling duration that increases the total flash counts  . Lightning flash illuminates a cloud top area significantly larger than the field of view of an LIS pixel. Also, the LIS detection threshold settings were not changed after boost. Thus, there is no change in detection efficiency. The inclination of the orbit is 35˚, which allows the LIS to observe the lightning activity in the tropical regions of the globe. It detects the lightning activity with detection efficiency of 90% with negligible regional bias  for flash with detection efficiency 93% ± 4% and 73% ± 11% during night and day respectively. LIS records both cloud-to-ground and intra-cloud lightning in terms of total flash count.
The present study involves LIS flash rate density data for a period of 10 years (1998-2007) over the IGP domain (21˚ - 32˚N and 73˚ - 89˚E) region as shown in Figure 1. The LIS flash count data were retrieved  . The monthly flash rate density is obtained from total lightning flash count and gridded with a spatial resolution of 0.5˚ × 0.5˚ for a 10-year period from 1998-2007. The monthly flash rate density (Fl·Km−2·Yr−1) over a year at each grid (0.5˚ × 0.5˚) is calculated by dividing flash count with view time and grid and by multiplying with total seconds of a year. India Meteorological Department categorizes four seasons namely; pre-mon- soon (March-May), monsoon (June-September), post-monsoon (October-November) and winter (December-February). Along with monthly flash rate density a 10- year period (1998-2007), rainfall data for the same period and region are extracted from Climatic Data Centre of National Oceanic and Atmospheric Administration (NOAA). The data for LI and K-index for Patna in the IGP region for the same period were obtained from University of Wyoming (http://weather.uwyo.edu/upperair/).
Figure 1. Map of Indo-Gangetic plains. Lifted index and K index for Patna station.
3. Results and Discussions
3.1. Spatial Variations of Flash Rate Density over the Indo-Gangetic Plains
Figure 2 shows the spatial variation of flash rate density over the IGP region during the period of study. The spatial variation shows that the maximum flash rate density is recorded over the northern region of IGP (40 fl·km−2·yr−1) to that of eastern IGP region. Thunderstorms followed by lightning are very frequent during the pre-monsoon season (March, April, and May) over central, east and northwest India and adjoining Bangladesh. During pre-monsoon season a westerly jet stream lies south of the Tibetan Plateau at upper levels while there is ascent over the Indian land mass and descent over the surrounding seas  . Anstrong solar heating at the lower levels, resulting in a large-scale circulation with inflow of warm moist air from the Bay of Bengal (BoB) takes place across the eastern coast, depth of which is about 500 m as detected at an inland station, at about 80 km from BoB   . Hence the air in the south eastern part is humid and moist while it is relatively dry and cold in the northwestern part, forming the dry line which has a significant role in the formation of severe thunderstorm activity.
During the monsoon, Madden Julian Oscillation (MJO) is largest element of the intra-seasonal (30 - 90 days) variability in the tropical atmosphere. The monsoon trough of low pressure at surface runs from the head of the BoB in the east to Ganganagar in the northwest India and is a semi-permanent feature of the monsoon circulation. The monsoon trough favors the formation of thun-
Figure 2. Spatial variation of average flash rate density over the Indo-Gangetic plains during the period 1998-2007.
derstorms and lightning in the north. When the oscillating MJO comes over the northern most regions that are over the location of the monsoon trough, it produces maximum thunderstorm activity. The landmass distribution and northward movement of the MJO contribute to increasing flash rate density with the latitude. The maximum activity observed over the north India as seen in Figure 2 may be due to formation of thunderstorms due to retreating ITCZ and synoptic scale systems  . When ITCZ approaches the equator from its northernmost convection associated with the large-scale circulation dominates over the landmass and orography distribution which produces a bimodal type of variation as seen in Figure 3.
3.2. Annual Variation of Flash Rate Density and Rainfall Rate over Indo-Gangetic Plains
Figure 3 and Figure 4 shows annual variation of average FRD and precipitation rate (mm/day). From these figures, it is seen that during the annual course, the curve for lightning activity show bimodality and rainfall activity exhibit a peak. The FRD shows a consistent increase from February that attains a maximum in May and thereafter it starts decreasing consistently till the end of August and again increases to a second lower peak in September. Pre-monsoon monthly (March-May) convective instability is shown in Figure 5(a) and correlation of air temperature (at 700 mb) with FRD is shown in Figure 5(b). FRD peaked in pre-monsoon before the rainfall peak (monsoon) as seen in Figure 3. This is attributed to environmental lapse rates (>7.0˚C/km) during pre-monsoon and al-
Figure 3. Annual variation of average flash rate density over the Indo-Gangetic plains during the period (1998-2007).
Figure 4. Annual variation of average rainfall rate over the Indo-Gangetic plains during the period (1998-2007).
so, evident from the altitudinal air temperature variation analysis (Figure 5(a)) and a high correlation (R = 0.70, p ≤ 0.0001) of air temperature (at 700 mb) with FRD during pre-monsoon (Figure 5(b)).
The bimodality of lightning activity is related to the variation in air temperature in the tropics. Variations in the incoming shortwave radiation from the Sun are attributed  to be ultimately responsible for variations in global meteorological parameters on the seasonal timescale. The rainfall activity also shows a similar increase where steep increase is in April and reaches its maximum in July, and then it shows decrease and reaches to its minimum in November. This observed peak of rainfall activity is the significant climatological feature pattern
of the monsoon rainfall and it is the difference between the ascending and descending phase (onset and withdrawal phase) of the monsoon  . The onset of the monsoon in early June is followed by the progressive increase in rainfall over the entire country. This phase attains its maximum in the month of July  . A careful examination of this figure suggests that there is a lag of nearly 2 months in the occurrence of peaks in lightning and rainfall activity. The plausible explanation for this lag is given below:
Prior to the onset phase of the monsoon, the Indian land mass region experiences maximum convective instability caused by insolation   . Hence the overlying boundary layer becomes more strongly buoyant in relation to its surroundings. The larger cloud buoyancy leads to the stronger continental up-
Figure 5. (a) Variation of air temperature with height in km; (b) Scatter plot of flash rate density and air temperature at 700 mb level for 10 year period (1998-2007).
drafts, which favors charge separation and lightning     .  made a study of the electrical properties of large number of deep convective thunderstorms that occurred over the Pune region through the measurement of point discharge current and electric field. Their studies revealed that the high values of electric field and point discharge current in the pre-monsoon season is an indication of the electrically active state of the thundercloud, under such conditions the cloud region electric field is usually intense to cause dielectric breakdown of the atmosphere and thus trigger lightning. Electrified thunderclouds are usually occurs at the end of April and throughout the May. The frequency of occurrence thunderclouds gradually reduces in the first week of June indicating the advancement of the monsoonal rainfall (Figure 3).
It is seen from Figure 6 that the negative LI is seen over March to September. The air is unstable and suggests the possibility of convection. LI are indicators of occurrence of thunderstorms. This is because air parcel rising is much warmer than its surroundings which can accelerate rapidly and create severe thunderstorms. There was sufficient moisture in the atmosphere which leads to an increase in the number of thunderstorm clouds     . The value of K Index lower than 15˚C  indicate very small likelihood of thunderstorms, while K Index values higher than 35˚C indicate high probability for the development of numerous and/or severe thunderstorms (as seen in Figure 7). The intense conditional instability over the region develops mixed phase clouds (water + ice) which subsequently produce the lightning activity.
Figure 6. Annual variation of lifted index over Patna station for 10 year period.
Figure 7. Annual variation of K-index over Patna station for 10 year period.
The paper presents the association of rainfall and stability index with lightning activity over Indo-Gangatic plain (IGP) region. The results indicate the following concluding remarks.
1) The spatial variation of flash rate density (FRD) and AOD is found to be higher by 40 fl·km−2·yr−1 over the northern region of IGP to that of eastern IGP region.
2) The annual variation of FRD exhibits bimodal distributions and the precipitation rate shows only one peak in its distributions. The time-lag of 2 months in the peak occurrence of FRD and rainfall rate may be attributed to the environmental lapse rates of more than 7.0˚C/km (thunderstorm condition) during pre- monsoon as is evident from the temperature profiles. Also, high correlation (R = 0.70, p =< 0.0001) of air temperature (at 700 mb) and FRD during pre-monsoon support this hypothesis.
3) The annual variation in LI distribution shows negative value from March- September and also, K-index during pre-monsoon and withdrawal phase of monsoon is due to intense insolation, availability of moisture and on account of large CAPE. The intense conditional instability over the region develops mixed phase clouds (water + ice) which subsequently produce the lightning activity.
4) The thermal heating diminishes on account of monsoonal rainfall and revival of monsoon rain even after break period, the LI is not strong enough to lift the moisture above freezing level. These warm clouds are maritime clouds with intermediate height, moderate updraft and less cloud electrification and hence minimum lightning activity in the monsoon season.
The Indian Institute of Tropical Meteorology (IITM), Pune is supported by the Ministry of Earth Sciences (MoES), Government of India, New Delhi. The authors are thankful to Professor Ravi S. Nanjundiah, Director of IITM, Pune, for his kind support and valuable guidance to carry out this research work.
 Soriano, L.R., Pablo, F. and Diez, E.G. (2001) Relationship of Precipitation and Cloud-to-Ground Lightning in the Iberian Peninsula. Monthly Weather Review, 129, 2998-3003.
 Ezcurra, A., Areitio, J. and Herrero, I. (2002) Relationships between Cloud to Ground Lightning and Surface Rainfall during 1992-1996 in the Spanish Basque Country Area. Atmospheric Research, 61, 239-250.
 Katsanos, D., Lagouvardos, K., Kotroni, V. and Argiriou, A. (2007) Relationship of Lightning Activity with Microwave Brightness Temperatures and Spaceborne Radar Reflectivity Profiles in the Central and Eastern Mediterranean. Journal of Applied Meteorology, 46, 1901-1912.
 Chaudhari, H.S., Sawaisarje, G.S., Ranalkar, M.R. and Sen, P.N. (2010) Thunderstorms over a Tropical Indian Station, Minicoy: Role of Vertical Wind Shear. Journal of Earth System Science, 119, 603-615.
 Litta, A.J., Idicula, S.M., Mohanty, U.C. and Prasad, S.K. (2012) Comparison of Thunderstorms from WRF-NMM and WRF-ARW Models over East Indian Region. Scientific World Journal, 2012, Article ID: 951870.
 Chakraborty, R., Saha, U., Singh, A.K. and Maitra, A. (2017) Association of Atmospheric Pollution and Instability Indices: A Detailed Investigation over an Indian Urban Metropolis. Atmospheric Research, 196, 83-96.
 Price, C. and Rind, D. (1994) Modeling and Global Lightning Distributions in a General Circulation Model. Monthly Weather Review, 122, 1930-1939.
 Williams, E.R., Chan, T. and Boccippio, D.J. (2004) Islands as Miniature Continents: Another Look at the Land-Ocean Lightning Contrast. Journal of Geophysical Research, 109, D16206.
 Manohar, G.K., Kandalgaonkar, S.S. and Tinmaker, M.I.R. (1999) Thunderstorm Activity over India and the Indian Southwest Monsoon. Journal of Geophysical Research, 104, 4169-4188.
 Kandalgaonkar, S.S., Tinmaker, M.I.R., Kulkarni, J.R., Nath, A. and Kulkarni, M.K. (2005) Spatio-Temporal Variability of Lightning Activity over the Indian Region. Journal of Geophysical Research, 110, D11108.
 Tinmaker, M.I.R., Aslam, M.Y., Ghude, S.D. and Chate, D.M. (2016) Lightning Activity with Rainfall during El Nino and La Nina Events over India. Theoretical and Applied Climatology, 1-10.
 Chate, D.M., Tinmaker, M.I.R., Aslam, M.Y. and Ghude, S.D. (2017) Climate Indicators for Lightning over Sea, Sea-Land Mixed and Land-Only Surfaces in India. International Journal of Climatology, 37, 1672-1679.
 Christian, H.J., et al. (1999) Global Frequency and Distribution of Lightning as Observed by the Optical Transient Detector (OTD). 11th International Conference on Atmospheric Electricity, Guntersville, 7-11 June 1999, 726-729.
 Bond, D.W., Steiger, S., Zhang, R., Tie, X. and Orville, R.E. (2002) The Importance of NOx Production by Lightning in the Tropics. Atmospheric Environment, 36, 1509-1519.
 Cecil, D.J., Buechler, D.E. and Blakeslee, R.J. (2014) Gridded Lightning Climatology from TRMM-LIS and OTD: Dataset Description. Atmospheric Research, 135-136, 404-414.
 Boccippio, D.J., Goodman, S.J. and Heckman, S. (2000) Regional Differences in Tropical Lightning Distribution. Journal of Applied Meteorology, 39, 2231-2248.
 Albrecht, R.I., Goodman S.J., Petersen W.A., Buechler, D.E., Bruning E.C., Blakeslee R.J. and Christian, H.J. (2011) The 13 Years of TRMM Lightning Imaging Sensor: From Individual Flash Characteristics to Decadal Tendencies. 14th International Conference on Atmospheric Electricity, Rio de Janeiro, 8-12 August 2011, 1-4.
 Weston, K.J. (1972) The Dry Line of Northern Indian and Its Role in Cumulonimbus Convection. Quarterly Journal of the Royal Meteorological Society, 98, 519-532.
 Sikka, D.R. and Gadgil, S. (1980) On the Maximum Cloud Zone and the ITCZ over Indian Longitude during Southwest Monsoon. Monthly Weather Review, 108, 1840-1853.
 Williams, E.R. (1994) Global Circuit Response to Seasonal Variations in Global Surface Air Temperature. Monthly Weather Review, 172, 1917-1929.
 Williams, E.R., et al. (2002) Contrasting Convective Regimes over the Amazon: Implications for Cloud Electrification. Journal of Geophysical Research, 107, 8082.
 Rutledge, S.A., Williams, E.R. and Keenan, T.D. (1992) The Down under Doppler and Electricity Experiment (DUNDEE): Overview and Preliminary Results. Bulletin of the American Meteorological Society, 73, 3-16.
 Halverson, J., Rickenbach, T., Roy, B., Pierce, H. and Williams, E.R. (2002) Environmental Characteristics of Convective Systems during TRMM-LBA. Monthly Weather Review, 130, 1493-1509.
 Cifelli, R., Petersen, W.A., Carey, L.D. and Rutledge, S.A. (2002) Radar Observations of the Kinematic Microphysical and Precipitation Characteristics of Two MCSs in TRMM LBA. Journal of Geophysical Research, 107, 1-44.
 Schulz, P. (1989) Relationships of Several Stability Indices to Convective Weather Events in Northeast Colorado. Weather Forecasting, 4, 73-80.