3.1 Contrasting ISM spells : Relative Humidity and Horizontal Wind anomaly
The reason behind the investigation of VSC focussing on the active and break ISM spell is explained in this section. This section brings out the contrast in the background dynamical and thermodynamic conditions during ISM active and break days. Figure 2 shows the composite mean horizontal wind component above the shaded mean relative humidity (RH) for all the nine years of active (left panels) and break (right panels) days at four different pressure levels; 1000, 850, 600, and 300 hPa. Strong (weak) low-level moist (dry) westerly jet winds are predominant during active (break) ISM over peninsular India. The conspicuous strong cyclonic circulation over central India at 600 hPa level is the contrasting difference between active and break ISM. The high-level anti-cyclonic circulation-driven tropical easterly jet feature is omnipresent over peninsular India during both the ISM spells. RH in the two ISM spells mainly differs at mid-level and high-level. Likewise, horizontal wind speed and direction mainly show contrast at 850 hPa onwards. The detailed discussion of using these background conditions for the observed VSC much evident in section 3.5.
3.2 Contrast in VSC during the active and break ISM spells
In order to investigate the VSC pertinent to ISM active and break period, we composited the nadir profiles of CloudSat swath data for 2008–2016. Here we used a contoured frequency by altitude diagram (CFAD; Yuter and Houze, 1995) analysis to explain the statistical distribution of the Ze. CFAD is a 3-D contour plot of histogram height profile that summarises normalized histogram analysis of Ze profiles as a representative cloud vertical structure. It also overcomes the mismatch in space and time scales of radar data. Typical CFADs of Ze are derived for eight different regions (see Fig. 1) over India for the active and break ISM periods, and composite VST are shown in Fig. 3(a-h) and 4(a-h), respectively. The striking difference between active and break ISM days is seen over WI, NW, CI, and CE region CFADs. It shows continuous or un-truncated cloud vertical structure during active ISM days due to the ubiquities presence of Johnson tri-modal convection, thus the interaction among the three-level estates within the VSC (warm, mixed, and ice phase cloud regions). On the other hand, break in the VSC (gap exists in the mid-level mixed-phase) because the clouds are mainly present in either or both warm and ice cloud regions during break ISM days. Space-borne CPR does not detect signals from heavy rain due to the signal attenuation from bigger rain drops. Hence decreasing values of Ze with decreasing height can be observed below the bright band. Therefore, the ZCFAD, in general, is a curved bow or arc shape especially with higher frequency count of high Ze values centered in the mixed phase region between 4–11 km for all the regions except NI, EI, and SI (WB) region during the ISM active (break) period. Bow shaped ZCFAD indicates a favourable condition to have both ice (warm) microphysical growth process from the top of higher (lower) estates in the VSC (e.g., Fig. 3(a-d) and Fig. 4(f-g)). In fact, perfect bow shape druing active ISM over NW and CE (Fig. 3b and 3c) indicating deep convective updrafts responsible to have high frequency of Ze values of above 10 dBZe above the − 40°C level. Convective precipitation over NW can also be confirmed by the absence of a bright band. Where the VSC is continuous, bright band feature is observed between 4–5 km altitudes, except over NI, NW, and SI, due to the sharp gradient in the value of the hydrometeor dielectric constant at the melting layer resulting in enhanced radar reflectivity. Bow shape is typical to CloudSat's ZCFAD mainly due to the attenuation with CPR (94 GHz) encountering only with lower level bigger rain drops, below the meling level extending much down to ground.
Over the WI region, VSC has developed up to 15 km. While higher (lower) values of Ze, in the range of above 0 dBZe with a frequency count more than 50% close to the ground, for active (break) ISM period indicate the formation of precipitating (non-precipitating) clouds in this region. Both on active and break days, the CFAD of WI is similar to that of CI and CE, which is also indicated by the rainfall anomaly map in Fig. 1. It further indicates a similar mechanism of formation and growth of warm and mixed-phase clouds over these three regions. The exception is at the high level, which shows over WI, more than 60–80% cloud occurrence above 8 km, denoting the predominant ice cloud processes during active days. At the same time, only ~ 30 (50) % of cloud can exceed 10 km over CI (CE). During the active spell, 50% of the precipitating cloud occurrence is prominent up to 6 km (100% at 2–4 km) in Fig. 3 (a,c,d), with an increase of Ze values from 0 to 10 dBZe. Figure 3(a,c,d) also depicts that 40% of cloud occurrence is attributed to ice cloud processes as the cloud height can touch 12 km from the ground. Secondary non-precipitating cloud patch in the order of -10 to -30 dBZe existing at the mid-level (6–8 km) is also commonly associated with the three regions. However, this mid-level mixed-phase cloud is a predominant feature over EI during both ISM spells but peculiar to notice precipitation seldom reaching the ground. Similarly, besides the mid-levels weak or cloud void region in the VSC, above 20% cloud occurrence is mainly limited to up to 5 km at low level and 10–16 km at high level over WI, CI, and CE during break days. At a low and high level, the maximum cloud frequency concentrated around − 25 dBZe depicts the dominance of smaller cloud particles or weak layer-type clouds. Besides this maximum frequency, there exists a secondary maximum (50–60%) around 0–5 dBZe at low level, indicating local precipitation over CI during break ISM. This ground precipitation inferring secondary maximum at low-level clouds can be seen only over CI and WB during break ISM. Predominant high-level clouds during active ISM show extension of its VSC down to the ground during break ISM days, but non-precipitation clouds in nature (above 0 dBZe values show frequency count below 30% close to the ground). Over WB and EI regions, during both ISM spells, VSC is continuous and reaches up to 16 and 14 km (Fig. 3(f-g) and Fig. 4(f-g)), respectively. During the active spell, VSC of WB shows less than 30% of 0 dBZe value between 8–10 km, whereas relatively higher frequency of above 5 dBZe values between 6–10 km indicates intense convection during break ISM days. VSC over EI during break days (Fig. 4g) is extention of its active ISM featue of dominent mixed phase supported by higher level ice microphysical growth. CFAD over the SI region shows high cloud occurrence from 6–14 km due to the presence of thick ice clouds or descending cirrus clouds. The absence of low clouds portrays the absence of ISM rain over this region. Since the high-level clouds lag with respect to the convection center (here it’s CI) during active days (Jiang et al., 2011), it explains the presence of predominent high clouds and the complete absence of low or mid clouds over the SI region. Due to the same reason, even its continuous VSC during SI break days having weaker frequency count of ~ 30% for the absence of ISM rain over this region.
Over the NI region, in both the ISM spells, low-level cloud top exists ~ 8 km with a weaker 30% cloud occurrence above melting layer. Further, active days are dominated by non-precipitating clouds (below − 20 dBZe). It is observed that the most significant contrast between active and break ISM spell is observed over all regions except WB and EI regions (Fig. 3(f-g) and Fig. 4(f-g)).
3.3 Occurrence frequency of warm, mixed, and ice phase clouds during the ISM active and break days
This section offers the study of the occurrence frequency of clouds inside the three different cloud zones depending on the phase of the cloud particles, namely warm or low (below 0° isotherm), mixed or mid (0 to -40°C), and ice or high (<-40°C) cloud region. Histogram analysis of CPR measured Ze normalized with the maximum occurrence value serves the purpose as shown in Fig. 5. As revealed in the previous section, along with the structure of the cloud, further contrast can be seen in their number of occurrences during the two ISM spells. Over WI and NI, the distributions of cloud at the low and high levels are negatively (positively) skewed during active (break) days indicating the dominance of higher (lesser) Ze. Interestingly over WI, cloud presence at the mid-level is two times higher in the active spell than in the break spell. During break days NW region shows maximum cloud occurrence below 20% except at the higher level. A meager difference in cloud frequency is observed over CI and CE regions at mid and low levels from small to medium-sized cloud drops (-50 to -5 dBZe) in active and break spells. Big cloud drops of Ze > 0 dBZe are significantly dominant at the low and mid-level during active days. In contrast, maximum frequency at -35 dBZe during break days dictates the predominance of weak high-level clouds. Bimodal distributions at low-level are salient features for both WB and EI regions, revealing the presence of rainy and non-rainy cloud types. Almost similar frequency of cloud over WI and EI during active and break spells leads to a similar type of VSC in the two ISM spells found in the previous section. Unlike CI, a higher value of Ze (> 5 dBZe) is more dominant during the break spell. EI region consists of high Ze value (< 10 dBZe) at low levels and higher (> 60%) occurrence of cloud in both the spells at all the levels, which is in contrast to other regions. All regions show that high-level cloud occurrence is maximum compared to the other two levels irrespective of the monsoon spells. Due to the unfavorable large-scale condition, SI region is deprived of the low-level cloud occurrence during active days and the maximum Ze value is limited to 5 dBZe, unlike other regions. One must not ignore the lesser Ze value at low levels compared to mid-levels over most regions which restricts our study from comparing the two levels over the same region. It might occur due to the attenuation from the rainy cloud of the small wavelength of CPR observed from space.
3.4 Latitudinal variation of VSC during the ISM active and break spells
The contrast of cloud vertical structure with different latitudes during active and break days for 2008–2016 is shown in Figs. 6 and 7, respectively. Besides the same eight understudied regions, a 3º x 3º box around the equator (EQ) is also included in this section to understand the northward propagation of the cloud belt from the equator. During Active days due to the northward propagation of the cloud belt towards the CI region, CI, CE, WI, and WB regions show complete VSC due to the frequent occurrence of deep convective clouds. On the contrary, EQ SI regions show the dominance of either low cloud or high cloud with the isolated presence of deep cloud towards the north of the EQ. Few cloud occurrence over NI shows the dominance of warm clouds. Over the EI region, after 28° latitude, cloud height decreases as the dominance of mid-level clouds and high clouds increases. Due to the few cases of deep cloud over the NW region, some of the latitudes show void of cloud in the whole VSC. Maximum Ze of 20 dBZe reaches above 10 km over CI compared to other deep convective cloud regions. During break days, cloud occurrence decreases severely over CI and NW, showing the occurrence of only warm clouds or multi-layered clouds and cirrus clouds reduced to the presence of only cirrus over NW. A sustained shallow warm cloud with the high-level cirrus cloud over the WI region leads to the broken VSC. The maximum Ze over those regions is reduced to 10 dBZe or below than active days. In contrast, a complete VSC region can be found over EQ and EI due to the incessant presence of deep clouds along the latitudes, in contrast to the active day’s VSC. The deep cloud is persistent over WB, similar to active days. However, during break days, the maximum Ze of 20 dBZe is limited to certain latitude. CE and SI regions show the mixed occurrence of shallow and deep convective clouds. Mid-level congestus are predominant in the NI region.
3.5 Background condition behind the formation of VSC
To investigate the background condition of these contrasting VSC, horizontal wind is plotted above the shaded relative humidity (RH) at the left and right panel in Fig. 2 for active and break days, respectively. South westerly wind is more intense at the surface during active days than break days as shown in Fig. 2a and 2e especially from central to south India. The gradual reduction of RH from the surface to the upper level is prominent during break days in Fig. 2(e-h). RH also carries the contrast between the two spells since it is 10% lesser than CI and WI during break days. The cyclonic circulation generated over the head of the Bay of Bengal and above is prominent at 850, and 750 hPA (Fig. 2(b -c)) is completely missing on break days. Certainly, this cyclonic circulation is one of the important factors in active days to provide ample moisture, especially at the mid-level required for the generation and sustenance of deep clouds. The active (break) condition is generally associated with an increase (decrease) of cyclonic vorticity and decrease (increase) of surface pressure over the central Indian monsoon trough region, and strengthening (weakening) of the low-level jet. In Fig. 2g and 2h, RH at 750 and 550 hPA below 70% over most of the CI and WI region indicate the lack of moisture during break days, limiting the growth of clouds.
During the active period of ISM, the mid-tropospheric vorticity over the northern Arabian Sea (AS) and northern BoB experience anomalous intensification. Due to synoptic-scale cyclonic anomalies at midlevel, the West coastal zones of India favor medium systems. At the same time, BoB and Myanmar's west coast are favorable for larger systems to develop. Enhanced stratiform precipitation areas and new convection are triggered due to southwesterly low- to mid-level flow, which is advected directly across the mountainous coastlines (Medina et al., 2010). Grossman and Durran (1984) proposed that deep convection over the BoB could be initiated by forced ascent against the mountains due to the blocking of low-level flow by the mountain. Due to upstream blocking, low-level winds start de-accelerating causing gentle lifting offshore and large vertical or horizontal wind shear near the mountains. This mechanism inhibits deep convection. At the same time, the thermodynamic and dynamic structure in the mid-troposphere can trigger deep convection in a conditionally unstable situation.
During the ISM active period, the NW region receives rainfall from the abrupt, intense convective systems formed due to the formation of lows over Pakistan and the adjacent region. This low-pressure system draws surface moisture from the BoB branch of ISM by eastward-moving low-level jet (LLJ). During the break time, LLJ direction is changed to south-eastward. There is no LLJ through peninsular India (Joseph and Sijikumar, 2004), resulting in a deficit of moisture transport over this region. Also, due to the arid conditions of this region, moisture drawn from the offshore is dried out due to prevailing dry conditions. It results in the drying of the mid and lower troposphere.
It is observed that during the ISM active spell, CI and CE regions are under the influence of synoptic systems. These synoptic systems form over the BoB and move North-westward along the trough line (Rao, 1976; Mooley and Shukla, 1989; Goswami et al., 2003). These synoptic-scale disturbances intensify the low-level southwesterly jet and withdraw moisture from BoB. It impacts precipitation and the vertical development of clouds. Mainly, adverse impact is observed over the CI and CE region during the different spells of monsoon. During the break period, the VSC observed over CE and CI is the same as the active spell. However, the occurrence frequency of lower reflectivity indicates the presence of thin, non-precipitating clouds.
Small systems form over the western foothills of the Himalayas because of low-level moist air, which is also associated with instability. This low-level moist air is capped by dry westerly flow aloft and then is lifted over the foothills (Medina et al., 2010).
During the active spell of ISM, precipitating clouds are formed over WB due to the convergence of moist monsoonal air with Himalayan elevations driven by daytime heating over the Tibetan Plateau (Houze et al., 2007). During the break period, monsoonal wind patterns over this region are normal to the orography leading to enhanced deep convection over this region and the concave indentation of Himalayan orography.
Over EI, distinct rainfall patterns were observed due to low-level zonal wind fluctuations and their interaction with regional-scale topography (Fujinami et al., 2014). Over this region, broad stratiform cloud systems are formed as a part of Mesoscale convective systems in association with BoB depressions (Houze et al., 2007). Monsoon depressions provide a moist, maritime environment favorable for convection, allowing the mesoscale systems to develop larger stratiform clouds. SI shows similarity in the VSC for both the ISM spells but layered clouds in higher levels are seen during the break period (Fig. 2h, 3h). VSC properties over WB and EI are observed fairly consistently during the active and break spells. Due to the weakening of descending flow in the mid-troposphere from the Tibetan Plateau generated by mechanical and thermal effects during ISM months which helps deep convection to develop (Sato and Kimura, 2007). SI region comes under the rain shadow region during ISM. SI receives rainfall from the north-east or retreating monsoon from October to December (Rao, 1976).
3.6 Occurrence of different types of clouds
Figure 8 shows the normalized cloud type frequency over the selected Indian regions derived using CloudSat 2B-CLDCLASS product. St and Sc clouds are not well separated in the 2B-CLDCLASS product and hence combined for further analysis. It is also documented that over the Indian region occurrence of St and Ns Cloud is nominal (Subrahmanyam and Kumar, 2013). Figure 8h shows that during both spells of monsoon, mainly three clouds, mid-level (Ac and As) and high level (Ci) are showing a higher percentage of occurrences. It is observed that As clouds are dominant during the active period. Previous studies of cloud types over the Indian region suggested the presence of Ac cloud over Indian landmass during ISM months (Subrahmanyam and Kumar, 2013). Unlike many boundary layer clouds, mid-level clouds cannot be maintained by surface moisture fluxes. Hence, the presence of mid-level clouds over the Indian subcontinent indicates the presence of large-scale convection.
In Fig. 8a, it is observed that Ci clouds are present most prominently over the NW (WB, CE, and EI) region during the active (break) ISM period. During the ISM, monsoon depressions are formed in head BoB resulting in the formation of Ci clouds (from anvils of deep towers). Higher-level strong Tropical Easterly Jet (TEJ) spread these clouds over the ISM region (Sathiyamoorthy et al., 2004). It is to be noted that Ci cloud occurrence is underestimated by CloudSat since it is not sensitive to optically thin clouds.
Figure 8b and 8c show As and Ac cloud’s distribution over all selected Indian regions. It is seen that Ac and As clouds most frequently occur over CE, WB, and EI regions. As discussed above, VSC over these regions is under the influence of monsoon depressions and is favorable for the formation of mid-level clouds. It is also observed that during the active (break) period, As and Ac clouds show higher occurrence frequency over NW (WB) regions. The relative occurrence of Ci, As, and Ac follow similar distribution over WI, NW, CI, CE, NI, and WB during active days and over WB, EI, and SI during break days. It indicates the proportionate relationship between mid-level and high-level clouds.
Sc cloud is a drizzling cloud that shows a near similar appearance in both the monsoon spells but has a slight increase in appearance during the active period (Fig. 8d). Unlike other cloud types, Sc occurrence is more than 10 overall in the regions except for SI and NW, especially during active days. Sc clouds frequency is higher over EI during both the ISM spells (except NW). Cu cloud frequency is reduced for the WI region during the break period. This reduction of the frequency of Cu clouds and the presence of Sc clouds over WI is associated with descending motion of monsoonal convection and intense dry air mixing from Western Asia at ~ 650 − 500 hPa (Krishnamurti et al., 2010). This inhibits the vertical development of clouds over this region. Sc and Cu clouds are more favorable over WB (during break) and over EI (during both ISM spells).
Rain-bearing Ns clouds (Fig. 8f) are present over all the regions during the active spell of the monsoon. These clouds are present over WB and EI during the break period, with maximum frequency observed over the EI region. Over the SI region, Ns clouds are absent during both spells of the monsoon. Dc clouds (Fig. 8g) are present with a higher percentage over CI (WB and EI) during the active (break) period. Earlier studies (Zuidema, 2003; Subrahmanyam and Kumar, 2013) documented a high presence of Dc over the BoB region. Seesaw in the occurrence of Dc, Ns, and As clouds (Fig. 8g) is predominantly seen over southwest regions: WI, CI, and NW (Foothills of Himalaya regions: WB and EI) of monsoon trough line in the active (break) ISM spell. It indicates the movement of the monsoon trough, which leads to a change in cloud occurrence pattern and hence, resulting in rainfall pattern.
3.7 Vertical profiles of LWC and IWC
The vertical profiles of LWC and IWC over the eight regions are displayed in Figs. 9 and 10, respectively. Profiles with the bad flags in 2B-Geoprof and 2B-CLDCLASS datasets are excluded as these may impact the calculation of LWC. It is observed that during the active (break) ISM period, LWC profiles show a unimodal (bimodal) distribution with peaks at 1.5 km (above 5 km). WI region (Fig. 9a) shows a maximum value of LWC (185 mg m− 3) at 2 km altitude during the active period, while the minimum LWC is observed over the NW region (< 1 mg m− 3 Fig. 8b) during break spells. Additionally, NW region possesses the minimum LWC compared to other places in both spells. CI and CE regions (Fig. 9c and 9d) show the prominent contrasting feature with higher LWC during the active period and lower LWC during the break period. Fairly similar values of LWC during both the spells are observed over NI and WB regions (Fig. 9f and 9g). EI region shows increased LWC at 5 km for both spells of monsoon. It is to be noted that CloudSat measures profiles from the sea surface. Over Himalayan ranges, the mean surface height is ~ 3 km, hence LWC Profiles calculated over the Himalayan foothills regions show the effect of orography.
The IWC exhibits unimodal distribution with skewness at altitude ~ 9 km during active days. During the break period, skewness is observed towards higher altitudes (> 9 km). It indicates higher ice sedimentation during active spells than break spells. Core monsoon regions and NWI show one order high values of IWC during the active period than the break period. The structure of IWC profiles during both the ISM spells similarized with the ground-based cloud radar estimation over WGs by Sukanya and Kalapureddy (2019) but the difference in spatial resolution of the two radars causes the difference in IWC values. CI (Fig. 10c) and WB (Fig. 10f) regions show the highest values of IWC (160 (65) mg m− 3) during the active (break) period. Maximum IWC is observed at ~ 9 km over CI (Fig. 10c) and CE (Fig. 10d) during the active period and over WB (Fig. 10f) and EI (Fig. 10g) during the break period. Webster et al. (1986) showed that monsoon trough shifts towards the foothills of the Himalaya during the break period leading to stronger updrafts over these regions resulting in higher IWC.
3.8 Vertical profiles of IER and INC
The vertical profiles of IER and INC over the eight regions are illustrated in Fig. 10 by dotted and cross-marked lines, respectively. The highest value of IER is observed over CI (58 µm, Fig. 10c) during the active period, while over EI (48 µm, Fig. 10g) during break period. A peak at IER values is observed at 5–6 km during the active period and 12 km during the break period over all the regions except NI, EI, and SI. Unlike other regions, ten times and 1.5 times higher IER can be seen over NI and EI, respectively, during the break spell than in the active spell. The exact opposite structure of IER maxima can be seen over SI, having a peak of IER at 5 km during the break spell and 10–12 km during the active spell. The most prominent contrast of IER is evident over WI, CI, and CE, having ten times higher IER during active days, indicating dominating ice process than on break days. A higher percentage of cirrus cloud may lead to a secondary peak in IER at 12 km, especially for the regions experiencing a break period. During the active (break) period, it is observed that larger ice hydrometeors with IER greater than 24 µm are observed over CI, WI, CE, WB, and EI (EI, NI, and WB) at an altitude of 5–6 km.
INC profiles show that number concentrations increased with height for both spells and were observed to be highest in upper levels. A higher value of INC is observed over CI (90 L− 1, Fig. 10c) and WB (48 L− 1, Fig. 10g) during the active period and the break period, respectively. The most gradual increase is observed over the regions affected by the monsoon trough, i.e., CI and CE (WB and EI) during the active (break) ISM spell. Hence, the presence of synoptic disturbances over the region might lead to lifting these hydrometeors higher in altitude.
From Fig. 10, it can be concluded that large (small) sized hydrometeors present in lower (higher) heights with lower (highest) concentrations. Using the CloudSat satellite and ground observation data, Deng et al. (2014) concluded that IER has the best correlation with precipitation at 5–6 km level. This correlation decreases with an increase in height above 5 km. Aircraft observations show a higher number concentration of ice particles, the existence of small and large cloud droplets, graupel, columns and rimed needles within the temperature range of -8 to -3°C. This indicates that the increase in the ice concentration is due to the Hallett-Mossop (HM) process of secondary ice formation (Patade et al., 2016). INC reaching above 5 L− 1 is representative of efficient production of ice hydrometeors by ice processes.
Further, it is also concluded that the HM process is responsible for graupel formation. During the active (break) ISM period, IER over the regions CI, WI, CE, WB, and EI (EI, NI, and WB) cross the threshold (5 L− 1) for the HM process to start ice production. It suggests favorable conditions for HM in mixed-phase clouds. It also underlines the importance of the role of mixed and ice phase microphysics in ISM.
3.9 CPR observation complemented with ground-based cloud radar
IITM’s Ka-band radar (KaSPR) vertical looking observation deployed over Mandhardev, Western Ghats, India, is quality checked (Kalapureddy et al., 2018) to complement space-based CPR observation. In Fig. 11–13, HTI (a-b) and CFAD (c-d) of Ze profile from CPR and KaSPR are shown for three particular cases when CloudSat observation was available over the radar site. KaSPR has high resolution (25 m, 1 s) Ze profile, and CPR has 1.7 km as along-track resolution. Hence, the two radar profiles for each case are different. The presence of deep cloud, thick high-level cloud, and shallow warm cloud cases are shown here. In all the figures, the latitudinal variations of VSC are in well agreement with the time evolution of VSC from KaSPR both in the HTI (a-b) and CFAD (c-d) map. Figure 11 shows cloud top height above 14 km in both the radar observations. Both the CFAD analysis indicates a maximum frequency above 10 km. In Fig. 12, the descent of thick high clouds and merging with low-level clouds are prominent in both the observations. Both the instrument depicts the frequent generation and dissipation of shallow cumulus clouds below 4 km in Fig. 13. Due to high sensitivity towards smaller droplets, cloud boundaries and cloud parts having Ze <-30 dBZe are well distinguished in KaSPR compared to CPR. The same reason leads to different maximum frequencies of Ze in the two radars in the CFAD graph.