Details of volcanic eruptions in the tropics during 1998–2005 having significant stratospheric impact
Particulates in the Upper Troposphere and Lower Stratosphere (UTLS) gained considerable interest due to their role in the dehydration of tropospheric air entering the stratosphere [1,2] as well as their potential to influence the radiation budget of the Earth-atmosphere system . The upper troposphere region, which is also conducive for the formation of cirrus, plays a major role in the transport of water vapour and other chemical constituents into the stratosphere. The physical processes responsible for maintaining the observed aerosol distribution in the tropical UTLS, the process with which it interacts with cirrus clouds and the effect of these particulates on the radiation budget are not fully understood . Studies have shown that the microphysical (such as particle shape, size, and size distribution) as well as the chemical properties of particles in the UTLS region [5-7] are mainly governed by the strength of tropospheric convection and the prevailing dynamics of the underlying troposphere. The formation and persistence of cirrus clouds in the upper troposphere is mainly governed by the concentration of available condensation nuclei in this region and their physical and chemical properties. These clouds are believed to be a significant contributor to atmospheric greenhouse effect [8,9] as well as hypothesized to play a major role in the dehydration of the lower stratosphere [1,2] and thus becomes an important factor governing global climate, through their positive feedback.
The last four decades of the 20th century have been marked by relatively intense volcanic activity  and hence long-term measurements of aerosols during this period mostly characterize the volcanically perturbed aerosols system rather than ‘background’ conditions [11-13]. During this period the increase in aerosol loading in the stratosphere could accelerate the heterogeneous chemistry of sulfate aerosols leading to a decrease in ozone amount [14-17], altering the NO2 concentration [18-20] and hence modifying the earth’s radiation budget . However, long-term studies on stratospheric aerosols show that on a global scale the stratospheric aerosol loading has returned to the pre-eruption levels (prevailed in the late 1970s) after the eruption of Mt. Pinatubo in 1991 [13,22] for the first time in 1998 and continued to remain almost at the same level for a couple of years. Only very few studies are carried out on the characteristics of these background stratospheric aerosols using observational data [7,13,23]. These studies, however, have shown that the global distribution of stratospheric aerosols will be significantly influenced by the atmospheric dynamics which includes periodic variations, such as the Quasi-Biennial Oscillation (QBO), seasonal cycles and long-term secular changes in addition to small perturbations due to the feeble volcanic eruptions and also due to degassing from the Earth’s crust. Identifying these secular trends in the background stratospheric aerosol system is crucial to predict future aerosol levels . While
The quasi-biennial oscillation in stratospheric zonal wind (QBOU) is found to influence significantly the distribution of volcanic stratospheric aerosols mainly over the tropics [28-30]. Even though, recent observational studies [26,31] revealed significant seasonal and inter-annual variations in the stratospheric aerosol load during the volcanically quiescent period, studies on the influence of these types of periodic oscillations on the stratospheric aerosol distribution over the tropics during volcanically quiescent periods are very rare . This study involves an attempt to study the features of particulates in the UTLS region during the relatively quiescent volcanic period (1998-2005) using global aerosol data from Stratospheric Aerosol and Gas Experiment (SAGE-II) archive and lidar data from Gadanki [13.5ºN, 79.2ºE]. Formation of semitransparent cirrus (STC) is very common in the upper troposphere. The characteristics of these STCs and their contribution to particulate scattering in UTLS region are also investigated.
2. Extinction/Backscatter data from LIDAR and SAGE-II
The biaxial, monostatic dual polarization Lidar at the National Atmospheric Research Laboratory (NARL), Gadanki, is used to study the scattering properties of atmospheric particulates in the UTLS region. This Lidar  is equipped with Nd:YAG laser (Model: PL8020, Continuum, USA) emitting linearly polarized pulses with 7 ns width and 20 Hz repetition rate at its second harmonic wavelength of 532 nm with a pulse energy of 550 mJ. The basic beam emerging from the laser source with a divergence of 0.45 m rad is expanded using a 10X beam expander to reduce the divergence to <0.1 m rad before transmitting vertically into the atmosphere. The time series of backscattered photons from different altitudes corresponding to each transmitter pulse are received using a 350 mm diameter Schmidt-Cassegrain telescope having a field of view of ~1 m rad. Both the transmitted beam and vertically looking receiving telescope are configured with a fixed horizontal separation of ~3 m. For this lidar configuration, as the lowest altitude at which the full overlap of the transmitter beam with the receiver field-of-view (beam-filled condition) is encountered around 7 km, the data from altitudes above 8 km only are used for retrieving the aerosol properties. A polarized beam splitter in the receiver beam path splits the beam into co-polarized and cross-polarized components which are detected independently using two identical photomultiplier tubes operated in photon counting mode and acquired with a bin width of 2 μs corresponding to an altitude resolution of 300 m. These photon-number profiles corresponding to each transmitted pulse are summed over 250 s to achieve a good signal to noise ratio up to altitude above ~40 km.
The SAGE-II onboard the Earth Radiation Budget Satellite (ERBS) employs solar occultation technique to measure the attenuation of solar radiation at the Earth's limb between the satellite and the Sun due to scattering and absorption by different atmospheric species . These measurements provide the altitude profile of the volume extinction coefficients of atmospheric particulates which includes particles of thin sub-visual cirrus clouds and aerosols at four different wavelengths in the visible and near-IR range (1020, 525, 453 and 385 nm) with a horizontal resolution of about 200 km and a vertical resolution of 0.5 km . This sensor takes 30 occultation observations on a single day, which are equally spaced in longitude round the globe but vary in latitude by a few degrees giving a near global coverage over a period of 2540 days. Details regarding the SAGE instrumentation and algorithms are discussed in earlier publications [32,34]. The upper limit of the particulate extinction measurable by SAGE sensor at 1020 nm is ~2 x 10-2 km-1, which is much larger than that in the UTLS region (~2 x 10-4 km-1) under volcanically quiescent period. Cirrus cloud with extinction greater than this value, generally referred to as ‘opaque clouds’ [35,36], are not measurable by the SAGE-II sensor. Presence of such clouds, limits the SAGE measurements below tropopause.
3. Estimation of particulate extinction from lidar data
The lidar data (backscattered signal) on different nights from Gadanki are used to derive the altitude profiles of particulate backscatter coefficient (βp) and volume depolarization ratio (δ) [37-39]. In the lidar system, received backscattered signal (at 532 nm wavelength) is separated into co-polarized and cross-polarized components (⊥ and ‖ channels, respectively) and recorded separately in two channels. The data in these two channels are analyzed separately employing the Fernald’s algorithm  to estimate the total backscattering coefficients (β⊥ and β║, respectively) taking 30 km as the reference altitude where the aerosol contribution is assumed to be negligible. For this inversion a value of 40 sr-1 is assigned for the lidar ratio (SP) and its variation with altitude depending on δ is also accounted appropriately. With this correction the value of SP reduces to ~26 sr-1 within the STC  in the upper troposphere. Further incorporating the correction for multiple scattering the value of SP within the STC reduces to 20 sr-1 (which is used to study the properties of STCs). The molecular backscatter coefficients for the two polarized components are estimated from the mean molecular number density profile taking a molecular depolarization factor (δm) of 0.028 . The molecular backscatter coefficient of the co-polarized component (βm⊥) is related to that of the cross-polarized component (βm║) as βm⊥ =δmβm║. Subtracting β m⊥ and βm║ from the altitude profiles of β⊥ and β║, respectively obtained from lidar data employing the Fernald’s algorithm, the altitude profiles of particulate backscatter coefficient, βp⊥ and βp║, are estimated. The respective backscatter ratios for the co-polarized (R⊥) and cross-polarized (R║) components are estimated  as R⊥=β⊥/βm⊥ and R║=β║/βm║. As far the net atmospheric backscattering is concerned, the ‘‘unbiased’’ or effective backscatter ratio (R) is to be defined, to quantify the gross property of the medium, which on mathematical simplification can be written as R(h)=[R⊥(h)+δmR║(h)]/(1+δm). The volume depolarization ratio is obtained from the ratios of R⊥ and R║ as δ(h)= [δm R║(h)]/ R⊥(h). This ratio is a good indicator for distinguishing the cirrus based on ‘particle habit’. While for small spherical particles, the values of δ will be relatively small, its value increases significantly as they become large and non- spherical. Using this property of cloud particles, structure and altitude extent of cirrus can be estimated from each lidar profile, which will be used to study the temporal variation of cirrus properties during the entire period of lidar observation. Based on a detailed scrutiny of a number of profiles at different cloud conditions a threshold value of δ ≥0.04 is assigned for discriminating the STC . If the value of δ exceeds this threshold value it is classified as STC. The effective particulate backscatter coefficient, βP, is the sum of βp⊥ and βp║. The altitude profile of particulate extinction coefficient, αP, is estimated by multiplying the altitude profile of βp with the corresponding profile of SP.
4. Altitude structure of backscatter/extinction, scattering ratio, and depolarization from lidar and comparison with SAGE measurements
Figure 1 shows the altitude profiles of mean particulate backscatter coefficient (βp), effective backscatter ratio (R) and volume depolarization ratio (δ) for a few nights at Gadanki obtained from Lidar data during the year 1999 as typical samples. In general, βp and R show a general decrease with increase in altitude (eg. 05 and 12 April 1999) in the troposphere and stratosphere. But on a few nights a significant enhancement is observable over a small region between 9 and 17 km. This sharp increase in βp is due to strong scattering from ice particles of thin STC layers, formed at these altitudes either by
A detailed error analysis  showed that the estimated values of βp is less sensitive to the variability in SP. For a given uncertainty of 25% in Sp, the maximum uncertainty in βp is 10% in the absence of clouds, ~15% for thin cirrus and ~30% for thick cirrus. For the same uncertainty in SP, the maximum uncertainty in the retrieved backscatter coefficient and effective backscatter ratio are around 0.6%, 2% and 10%, respectively, for clear atmosphere, atmosphere with thin cirrus and atmosphere with thick cirrus. Including the possible errors in the lidar signal inversion associated with the uncertainty in the molecular backscatter coefficient, the resultant error in the derived optical depth would be ~20%. As the signal-to-noise ratio is >2 up to ~45 km, for altitudes <30 km the system induced errors will be significantly small (<1%) compared to that from other sources. The error due to the influence of background and system noise is negligible (< 0.001%) compared to that due to the uncertainty in Sp.
Altitude profiles of αp at 525 and 1020 nm from Global SAGE-II aerosol data archive (version 6.2) during the period 1998-2005 are obtained through the NASA website
5. Semitransparent cirrus in the upper troposphere
For studying the role of deep convection on the genesis of STC, the monthly mean spatial distribution of the frequency of occurrence of deep clouds (FD) with TIR brightness temperature < 235 K derived from KALPANA-1 during January-December 2005 are presented in Figure 4. A comparison of Figures 3 and 4 suggests that that the longitudinally extended band of high FSTC (~40-60%) in the region 10°S-20°S over the western tropical Indian Ocean and around the equator over the eastern Indian Ocean in January (Figure 3) is closely associated with the deep convection linked to the Inter tropical convergence Zone. However, the values of FSTC far exceed FD over all these regions. During this month, the highest values of FSTC (~60%) are observed at the north of Madagascar in the western Indian Ocean and over Sumatra/Indonesia in the eastern Indian Ocean, where the frequency of occurrence of very large deep convection is quite large. A region of less cloudiness (with FSTC<15%) is observed in the central Arabian Sea and Indian Peninsula centered around 5ºN-15ºN, which runs parallel to the equatorial band of high STC occurrence and is well separated from the deep convective regions. A similar STC-free zone is also observed in the southeast Indian Ocean centered on 20ºS. Clearly, these regions with low occurrence of STC are caused by the large subsidence in the upper troposphere associated with the descending limb of the Hadley circulation cell.
6. General features of STC over the Indian region
Lidar studies from Gadanki indicate that the occurrence of STC in this geographical region is the largest from May-October, associated with the formation of intense convection and the subsequent onset of Asian summer monsoon (ASM). Figure 5 shows the mean feature of STC occurrence and its altitude extend averaged for the period 1998‐2003. While the STCs occurring during the period May- October are generally thick and optically denser , those occurring during the rest of the year (dry months) are relatively thin (both geometrically and optically). As the Gadanki region is almost free from deep convection (Figure 4) during winter, the STCs forming during this period could be of
On examining the frequency of occurrence of STC during the period 1998‐2003 (Figure 6a) it can be seen that the frequency of occurrence of SVC is much larger than TC and DC. In general, STCs are observed in the altitude region 10 to 18 km with a preferred altitude (frequent occurrence) region between 14 and 16 km (Figure 6b) Thin clouds occur more frequently than thick clouds (Figure 6c). Though the vertical extent of STC generally vary from 0.4 to ~ 4.0 km in majority of the cases it is less than 1.7 km. Though the volume depolarization,δ, in these clouds varies in the range 0.03 to 0.6 its distribution peaks (Figure 6d) in the lowest value. The value of δ for SVC and TC are generally very small compared to DC. The particulate depolarization (δp) of STC (Figure 6e) generally varies from zero to unity. The distribution of δp peaks around a value of 0.15. The properties of STC vary significantly with cloud temperature (or altitude). Figure 7 shows the variation of thickness, depolarization and optical depth of STC with cloud temperature. As can be seen from Figure 7a, the thickness of STC is a maximum for temperatures in the range 55° to 75°C. Above and below this temperature range the cloud thickness decreases. Similarly the depolarization is maximum around 75°C and decreases steadily with increase in temperature. The value of δ also shows a sharp decrease when the cloud temperature decreases below 75°C. The decrease in δ with increase in temperature can be attributed to the melting and evaporation of ice crystals and subsequent blunting of their edges . Decrease in temperature leads formation of particles with sharp edges. But when the temperature decreases below a threshold value, the particle size  becomes small (needle type). This can lead to a decrease in δ. Figure 7c shows that on an average τc increases with increase in temperature. The cloud becomes more opaque at higher temperatures. This has important implication in the radiative effects of STCs  in the context of global change.
7. Mean annual variation of particulate scattering in the UTLS region
Figure 8a shows a contour plot of the logarithm of mean βp as a function of month and altitude for the period 1998–2003. The values of βp are relatively large during the May to September period and small during the winter months. This is due to the influence of STC. Relatively high values of βp in the UT region during the monsoon period are due to the presence of relatively dense cirrus and low values during winter are due to the presence of SVC. In contrast to the UT region, βp in the LS is generally large (as high as 6 × 10-8 m-1 sr-1) during the winter (November to January) and pre-monsoon (April–May) months and low (as low as 10-9 m-1 sr-1) during the summer (July and August) months. Prominent peaks are observed during May–June and November–January periods with low values in July–August and February. The mean annual pattern of αp at different altitudes from the lidar derived βp (Figure8a) is presented in Figure 8b. Similar grading scheme is used in both these plots to make an easy direct visual comparison of the pattern. The major features of the annual variation of βp and αp in different altitudes are very similar in these two plots, even though the lidar ratio is assumed to be variable with altitude depending on aerosol properties.
Microphysical properties like the size and shape of particles in the UTLS region can be delineated from the depolarization of backscattered radiation [57,58]. Figure 8c shows a contour plot of monthly mean δ in the altitude region 8–28 km. The value of δ varies in the range 0.03 to 0.6 in UT region and from 0.03– 0.04 in the LS region. High values of δ are generally confined to a narrow altitude region (14–16 km) during winter while these extend to a wider altitude region (12– 18 km) during summer monsoon period. High values of δ (>0.2) are observed in the altitude region 14– 16 km during April–October period. Values of δ exceeding 0.04 observed in the UT (above 10 km) are mainly due to presence of highly non-spherical ice particles associated with STC . The overall low values of δ suggests that particles in the LS region are very small and tend to become more spherical in nature.
Figure 9 shows a contour plot of αp similar to Figure 8b but generated from SAGE-II derived mean particulate extinction at 525 nm and 1020 nm over a small geographical grid size of 10-16°N and 73-86°E centered around Gadanki. As thin STCs in the UT region significantly attenuates the SAGE-II wavelengths, (especially that at 525 nm) there will be a large data gap at the lower altitudes. The region bound between X-axis and the rectangular vertical bars (shown white) are the data gaps. As the attenuation for 1020 nm is less than that for 525 nm this wavelength can penetrate to lower altitudes to yield useful data. The data gap is relatively less for 1020 nm. However, in generating the contours, the data gap is appropriately interpolated. General similarity of the pattern in Figure 9a and 9b suggests that the interpolation did not influence the major features of Figure 9b. Except for an overall decrease in the values of αp derived from SAGE-II (at 525nm) compared to those derived from lidar data, the major spatio-temporal features in Figure 9 also matches well with those of Figure 8b. Thus, the inferences derived from lidar data is reconfirmed by SAGE-II observations during the same period.
To make the features more concise, the month to month variation of altitude weighted δ (altitude integral of δ normalized to the slab thickness) for different altitude regions with a slab thickness of 2 km are examined in the UT and LS region. Though a sharp definition of UT and LS region is rather difficult, for the present analysis we use the term UT for the altitude region from 10 to 16 km and LS from 18 to ~30 km. Though the cold point tropopause shows a small variation with time of the day and day of the year, it lies always in the range 16–18 km (a transition region), such that the region defined as UT is always below cold point and region defined as LS is above. Figure 10a shows the mean annual variation of the cold point tropopause altitude and tropopause temperature.
The annual variation of δ for different altitudes in UT is shown in Figure 10b and that in LS in Figure 10c. Though, in general, on an average, the value of δ in LS region is < 0.04, it shows a pronounced oscillation with prominent peaks in May and September–January period. Above 20 km, the value of mean δ (~0.03) is very close to the molecular depolarization, it is relatively large (in the range 0.03 – 0.04) in the altitude region 18–20 km just above the cold point tropical tropopause. Thus, in the LS region, the particle size and non-sphericity decreases with increase in altitude. Note that, the value of δ encountered in this region is less than the threshold value (0.04) used for identifying the cirrus particles. The value of δ is largest in the altitude region 14–16 km, just below the mean level of tropical tropopause. The general similarity in the annual variation of δ in the altitude regions 14–16 km and 10-12 km indicates that the particle microphysics in these altitudes are strongly coupled. Even though the cold point tropopause altitude shows a small variation from month-to-month, the mean level lies around 16.5 km. Figure 10d shows the mean annual variation of the altitude weighted δ around the cold point tropopause for a slab thickness of 2 km. The observed general similarity of the annual variation of δ in the entire region from 10 to 20 km indicates that the annual variation of the particle habit (size and shape) in the UT and LS regions are strongly coupled. The observed general decrease in δ from UT to LS suggests that the particles tend to become small and more regular in shape with increase in altitude.
8. Effect of tropospheric convection on UTLS particulates
During the ASM period, the upper troposphere is significantly influenced by the tropospheric convection. This in turn could influence the microphysical properties of UTLS particulates. The outgoing long wave radiation (OLR), which is directly influenced by cloud cover and in turn by convection, could be effectively used as a proxy to the strength of tropospheric convection  and the seasonal variation of deep convection. The value of mean OLR is generally high for clear sky condition (~460 W m-2) and decreases with increase in cloud cover and cloud vertical extent. A threshold value  of 200 Wm-2 for OLR which when declines below could be treated as in index for deep convection. The percentage of occurrence of daily mean OLR with its value < 200 Wm-2, within a grid size of 2.5° × 2.5° over a geographical region 12.5° –15° N in latitude and 77.5° –80° E in longitude obtained from NCEP/NCAR reanalysis provided by the Climate Diagnostic Center through their Web site at http://www.cdc.noaa.gov/ averaged in different months for the period 1998–2003, which is almost negligible during the December to April period, starts increasing from May and shows a broad peak during the period June to October (Figure 11a).
The strength of tropospheric convection also can be assessed from the convective available potential energy (CAPE), defined as the altitude integrated buoyancy of lifted air parcel from the level of free convection to the level of neutral buoyancy , is estimated from the altitude profile of potential temperature as CAPE = ∫ [g(qv qT)/qT]dz, where qv is the virtual potential temperature of the air parcel, qT is the environmental potential temperature and g is acceleration due to gravity. The integral is taken from the level of free convection to the level of neutral buoyancy. The mean values of CAPE, which is a measure of stability in the atmosphere as far as vertical displacements are concerned , in different months for Chennai (13°N, 80.2°E), a station located very close to Gadanki, obtained through India Meteorology Department (Web site http://www.weather.uwyo.edu/upperair/ sounding.html.) are presented in Figure 11b. The CAPE also shows a prominent peak during the period April to May along with a small secondary peak in September. Intense convection is closely related to thunderstorm activity. Climatologically averaged thunderstorm activity in different latitude belts over the Indian subcontinent was investigated in detail in earlier study . The annual variation of the mean number of thunderstorm days for the latitude belt 10–15°N (Figure 11c) shows two prominent peaks in May and September–October period. Deep convection and high thunderstorm activity leads to the formation of thick convective clouds. The main reason for low lidar observation statistics during the monsoon period is the presence of these convective clouds which impede the observations. Though these clouds will have a large spatial extent and persist for several days, some gap region which is devoid of thick clouds (favouring lidar observation) starts developing after a few days The outflow from adjacent convective anvils spreading over to this gap region leads to the formation of STCs.
Being originated from convective outflow [43,44], especially during the summer monsoon period, particles of STCs in the UT region during this period will be relatively large and highly non-spherical [8,45]. Presence of these large non-spherical particles leads to an increase of δ in this region, as is observed. Subsequent uplift of some of these particles along with tropospheric air across the tropopause leads to an increase in δ along with the integrated backscatter (ΙβP) in the region just above the cold point tropopause. As the particles and precursor gases are not directly injected into the stratosphere but diffuse slowly across the tropopause they confine to a small region just above the cold point tropopause. Note that during the period July to August when CAPE as well as thunderstorm activity shows a small decrease, δ (Figure 12) and ΙβP (Figure 13) in the 18–20 km also shows a decrease. However, the STCs during this period shows large spatio-temporal variations  introducing large day-to-day variability in ΙβP and δ in the UT and LS region as indicated by the large error bars during the May-October period.
Though during the December–February period CAPE is relatively small and thunderstorm activity also is less, a broad peak in ΙβP at 18–20 km is observable during this period (Figure 13).This feature though appears to be contradicting to the above hypothesis proposed for the summer monsoon period, the real mechanism responsible for the transport of tropospheric air into the stratosphere during winter is different from that during the monsoon period. During winter the intrusion of tropospheric air into the lower stratosphere occurs through vertical ascent of air driven by the strength of the ascending branch of the Hadley cell circulation in the troposphere and Brewer- Dobson circulation (BD) in the UTLS region. The BD circulation is primarily driven by mechanical forcing (westward directed wave drag) arising from the dissipation of planetary-scale waves (breaking of Rossby waves) in the extra tropical stratosphere . Because of filtering by the large-scale stratospheric winds, vertical propagation of planetary waves into the extra-tropical stratosphere occurs primarily during winter and this seasonality in wave forcing accounts for the winter maximum in the BD circulation [65,66]. In association with these processes, the cold point tropopause at tropics is pushed to a higher altitude during winter and they play a major role in exchange of mass from troposphere to stratosphere.
As can be seen from Figure 10a the tropopause is cooler and higher during the January–March period and warmer and lower during August–October period which is in accordance with what is reported for tropical locations [67-70]. These Seasonal changes in the cold point tropopause altitude also contribute to the mass influx to stratosphere [71-73]. Transport of constituents like water vapor, aerosols, trace gases etc across the tropopause which acts as a permeable membrane through which continuous mixing between tropospheric and stratospheric air takes place is well demonstrated by various investigators [74-77]. Over and above the penetration of deep convection, this annual adjustment of tropopause altitude is an important mechanism for the stratosphere troposphere exchange.
9. Contribution of STC to particulate scattering in the UTLS region
Fine ice crystals of STCs originating either through the outflow from deep convective anvils or through freeze drying of moist air lifted up to the tropopause by normal convection contribute significantly for scattering in the UTLS region. In addition, other aerosol particles originating from the surface (mainly through bulk to particle conversion and those of vegetative origin) as well as those formed in the upper troposphere through gas- to-particle conversion (mainly of various industrial gases) also contribute for the particulate loading in the UTLS region. It would be worth examining the relative contributions of the two components in the UTLS region to make a quantitative assessment of the contribution of STCs to the scattering properties in the UTLS region. Figure 8a shows the contour plots of month-to-month variation of mean βp derived from lidar data during the period 1998–2003 (including both STC-contaminated and STC-free). Note that, as this study period was mostly devoid of major volcanic eruptions  the stratospheric aerosol loading can be considered to be in its background level.
While βp in the UT is a maximum during the May –September period and minimum during October–November, in the lower stratosphere it is minimum during summer (July and August) and maximum during winter. However, the winter high in the LS region was attributed  to the transport of tropospheric air (containing aerosols and precursor gases) in conjunction with the tropical upwelling and B D circulation while the observed high during May–June is due to the upward influx of particles (including ice crystals of STCs) from the UT region. As the uncertainty associated with lidar-derived βp is small compared to αp , for the lidar based study the altitude profile of βp and ΙβP for a desired altitude region (layer integrated) is used for studying the annual pattern of particulate scattering in the UTLS region. Figure 14 shows the contour plots depicting the annual variation of βp at different altitudes for the STC-contaminated case and the STC‐free case separately.
The cloud cover being quite large during the summer monsoon period no lidar data were available without STC during the July and August months, which lead to a data-gap for the STC-free case presented in Figure 14b. As can be seen, the variations of βp in Figure 14a is very similar to those in Figure 8a generated by considering both the STC-contaminated and STC-free profiles even though the absolute magnitude in the UT region is slightly small in the latter case especially during the October–November period. High values of βp are observed in the UT region for the STC-contaminated case during summer monsoon period. Figure 14b shows that, βp in the UT region is very small for the STC-free case compared to that of STC-contaminated case for the same period. The values of βp in the UT region are high during winter and spring and low during summer and autumn for the STC-free case. This annual pattern of βp is significantly different from that of STC-contaminated case. From this, it is quite reasonable to infer that the prominent peak of βp in the UT region observed during the summer monsoon period (in Figures 8a and 14a) is due to the influence of STC. However, the enhancement in βp in the UT region due to STC-contamination is relatively small during the winter months
As seen from Figure 5 the occurrence of STC is mostly confined to the uppermost part of the troposphere, above ~10 km. Because of their presence these STCs directly contribute to the particulate scattering in this region. Over and above, these prevailing STCs in the upper troposphere can also modify the scattering property of particulates (aerosols) above the cloud-top as well as below the cloud-base. These effects in the UT and LS regions can be inferred by examining the ΙβP at four different altitude regions, 8–10 km (UT1), 12–16 km (UT2), 18–21 km (LS1) and 21–25 km (LS2). Among these, the LS1 and LS2 regions are above the STC-top while UT1 is below the STC-base. The month-to-month variation of mean ΙβP for the three cases; viz, (i) including all profiles (ii) considering only STC‐contaminated profiles and (iii) Considering only STC-free profiles, are presented in Figure 15. In the absence of STC the ΙβP in the UT2 region is significantly small. The STC contribution to ΙβP in this region works out to be around 93±5%. The ΙβP shows an annual variation with relatively high values during the winter/dry months and low during the summer monsoon period for the STC-free case. When STCs prevail, the ΙβP in the UT2 region increases significantly especially during the summer monsoon period. Even though the enhancement in ΙβP due to the presence of STC is seen during the winter/dry months also, its magnitude is relatively less. It would be worth in this context to note that, during the monsoon period, the UT2 region is dominantly influenced by dense STCs originating from the outflow of convective anvils . The particle associated with these STCs will be relatively large and highly non-spherical [8,45] and hence their contribution to βp and δ will be significantly large .
It is quite interesting to note that the presence of STC in the UT2 region significantly enhances the ΙβP in the lower stratosphere also especially during the May–October period. In the LS1 and LS2 regions ΙβP shows a double peak structure with prominent maxima during May–June and October–January for the STC-contaminated cases (in the UT region), while it is very small from May to October for the STC-free case. A significant amount from the abundant ice crystals (of STC) present in the UT region during the summer monsoon period will be lifted up along with air-mass across the tropopause. Consequently, ΙβP in the lower stratosphere also increases. However, unlike the case during winter, this increase in βp (or ΙβP) is rather confined to a small region just above the tropopause (the LS1 region). All the three curves in Figure 15 for the LS1 and LS2 regions almost overlap each other during the December to April period. This indicates that the STCs in the UT2 region during this period did not contribute significantly to ΙβP in the stratosphere. In general, the influence of STC to the mean backscattering coefficient (ΙβP /Δh, where Δh is the slab thickness) in LS1 region is relatively more than that in LS2, particularly during the summer monsoon period. Thus, presence of dense STC in the UT2 region (during the summer monsoon period) enhances the particulate scattering in the lower stratosphere, even though this region is practically free from STC.
As lidar data yields reliable βp values only from the region above the altitude of ‘full beam overlap’ the effect of STCs below the cloud-base can be examined by studying the scattering property of the medium over a narrow altitude region 8–10 km (UT1), which is mostly free from STCs (Figure 6b). The mean annual variation of ΙβP in the UT1 region for the three cases, (i) when the UT2 region is STC-free, (ii) the UT2 region is STC-contaminated and (iii) with both these cases combined, is presented in Figure 15d. Except for the STC-free case, though the ΙβP in this region also shows an annual variation similar to that for the UT2 region, the amplitude of this variation is significantly small. Notwithstanding the fact that the values of ΙβP during the November to April period for the three cases are comparable, for the STC-free condition it shows a decrease during the May–October period. This shows that the presence of dense STCs in the UT2 region enhances the particulate scattering in the region below the STC-base and the magnitude of this enhancement is much smaller than that in the UT2 region (where the ΙβP increases by a factor >20). But when these STCs are thin (optically as well as geometrically), especially during the November to April period, this contribution is almost negligible. Thus, the prevailing STCs in the UT2 region during the summer monsoon period significantly enhance the scattering from the region above the cloud-top as well as below the cloud-base. The enhancement in particulate scattering in the LS region above the cloud-top could be mainly due to the lofting of STC-particles across the tropopause, which joining with the prevailing LS aerosols modify the volume scattering properties in this region. Though the Brewer–Dobson circulation is weak during summer the troposphere–stratosphere exchange during this period could be aided by an increase in wave activity depending on the prevailing atmospheric condition. Strong convection prevailing over the Indian landmass during the summer–monsoon period, which can transport abundant moisture to the upper troposphere inducing cirrus formation, is also a major source for gravity waves .
A study on the influence of STC in the particulate extinction in UT and LS region using SAGE-II data over the Indian longitude sector (70-90ºE) from 30ºS to 30ºN also revealed similar results. However,, as the occurrence of STC is less frequent in the southern hemispheric off-equatorial region, the increase in τp in the UT and LS regions due to the influence of STC is relatively small compared to that in the equatorial region .
10. Mean latitude variation of the altitude structure of τP over the Indian longitude sector
The latitudinal structure of the annual pattern of τp in the tropical UTLS region over the Indian longitude Sector is examined using the altitude profiles of particulate extinction at 525 nm obtained from SAGE-II data archive in the latitude region 30°S to 30°N for the period 1998-2005. The profiles are grouped in different latitude bands each having a width of 5°. Contours presenting the annual variation of τp in different latitudes bands in the UTLS region are presented in Figure16. The mean annual variation of the altitude of the lapse rate tropopause for each of these latitude bands is superposed on the respective contours.
For the latitudinal region between 0-15°N, the tropopause is cooler (higher) during the December-May period and warmer (lower) during July-October period, in accordance with that reported for tropical locations by prior investigators [67,68]. For the latitudinal sector 15-20°N, the tropopause is higher during April-June and lower during the July-August period. At latitudes > 20°N, a pronounced maximum in tropopause altitude can be observed during the boreal summer and minimum during boreal winter. Though the tropopause altitude varies with latitude, time of the day and day of the year, on an average, this altitude mostly lies in the range 16-18 km.
In general, αp is relatively large (>10-7 m-1) in the UT region for the region north of 20°S. In both the hemispheres αp shows two peaks. The summer peak (in the respective hemispheres) is more prominent, compared to the winter peak. South of 20°S, αp in the UT region is relatively small and does not show any pronounced annual variation. The summer-winter contrast in the UT region is almost insignificant beyond 25°S. The winter peak in τp becomes relatively weak in the region north of 15°N and becomes almost insignificant beyond 25°N. The summer-winter contrast in αp is well pronounced beyond 20°N. In contrast to the UT region, the values of αp in the LS region are relatively small (<10-7 m-1), at least by a factor of two. The annual variation of extinction shows a weak summer winter contrast in the LS region for all the latitudinal sectors between 0-30°N. Due to the influence of STCs, which occur at random in different altitudes, the standard error of αp (expressed as percentages of mean αp) is generally very large in the UT region. The mean standard error in the attitude region 10-15 km ranges from 20% to 60%. This error decreases progressively with increase in altitude. In the LS region, the error is very small and mostly confined to values in the range 5-15%. In the transition region 15-20 km, the mean error is of the order of 30% which is less than that in the UT region.
Figure 16 shows that the annual variation of αp in the UT region over the southern hemisphere is distinctly different from that over the northern hemisphere. In this hemisphere high values of αp (in the UT region) remains fairly confined to latitudes north of 20°S. This difference can mainly be attributed to the corresponding difference in the pattern of deep convection in these two hemispheres. While the southward migration of ITCZ (during the boreal winter) is fairly confined to latitudes north of 20°S, northward migration (during the boreal summer) extends beyond ~20°N. More over, the convection over the northern hemispheric land mass is much stronger than that over the southern hemispheric oceanic regions. Because of these features, the increase in αp in the UT region is also confined to the geographical region north of 20°S. The value of αp in the LS region is relatively small and show similar variations in both the hemispheres. This indicates that the features of aerosol transport in the equatorial (±15°latitude) LS region are fairly symmetric with respect to equator.
11. Tropospheric convection over the Indian region
The strength of tropospheric convection which can be indexed based on the CAPE could in turn also be related to the thunderstorm activity. The climatology of monthly mean thunderstorm days (TSD) for the Indian region was studied  using the data for the period 1970-1980. The monthly mean values of CAPE at different latitude sectors for the period 1998-2003 over this region is estimated from the respective daily values in each month obtained from the website of IMD (www.weather.uwyo.edu /upperair/ sounding.html). The annual variation of these two parameters for different latitude bands are shown in Figure 17. On an average, the annual pattern of TSD matches well with that of CAPE in all these latitude belts except for the fact that the secondary peak in October is less prominent in CAPE for latitudes north of 15°N. Both these parameters show a significant positive correlation with coefficient exceeding 0.7. Two prominent peaks observed during April-May and October in the near equatorial region merges to become a broad peak in the latitude region 20-25°N (off-equatorial) and subsequently becomes a well defined sharp peak around June-July in the latitude region 25- 30°N. Examining the annual variation of αp in the light of the annual variation of CAPE/TDS, it can be seen that the occurrence of the high values of αp (during May and October in Figure 16 h,i, j and those during June, July and August in Figure 16k,l), coincides with the peak in thunderstorm activity (or CAPE) in the respective latitudinal belt. During high convective activity while high values of αp are observed very close to the tropopause in the latitude belt 0-15°N, significant high values of extinction are observed just above the tropopause in the latitude region between 15-30°N. Intense thunderstorm activity along with deep convection during boreal summer is highly favorable for the formation of dense STCs very close to tropopause leading to an increase in αp in the UT region. Penetration of particles from these STCs to higher altitudes increases αp in the LS1 region, as is observed in Figure 16.
12. Influence of dynamics in the latitude variation of particulates in the UTLS region
To delineate the latitude variation of αp in the UTLS region over the Indian longitude sector (70-90°E), the annual variation of monthly mean τP for a latitude band of 5°width in the latitude region 30ºS to 30ºN, for the three altitude regions 12-16 km (UT), 18-21 km (LS1) and 21-30 km (LS2) are examined separately. Figure 18 shows the contour plots of the mean τP in the UT, LS1 and LS2 regions with month along x-axis and latitude along y-axis. In the UT region, τP shows a general decrease with increase in latitude from equator, with its gradient showing a pronounced variation from month-to-month. The summer-winter contrast (with relatively low values during winter and high values during summer) is well discernable in the UT region (Figure 18a). Relatively high values of τP observed in the UT region between 15°S-15°N during May to February period are mostly due to presence of dense STCs resulting from the outflow of convective anvils. These are the periods when convective activity in the troposphere is very strong (associated with the southwest and northeast monsoons) in this region. Above 15°N, relatively high values of τP are more-or-less confined to the June-August period when the monsoon trough usually reaches its extreme north over the continent. However, beyond 15°S, the values of τP decreases significantly with increase in latitude, compared to that observed in the northern hemisphere.
The latitudinal variation of the annual pattern of τP in the UT region shows fairly good correspondence with CAPE and mean thunderstorm activity (Figure 17) over the northern hemisphere. The annual variation of TSD in the latitude region 8-15°N shows two peaks; one during the April-May period and then another during October. Further north, these two peaks get closer and merge to form a single peak in the latitude belt 20-30°N during the June-September period. This corresponds well with the latitudinal variation of the annual pattern of τP in the UT region (Figure 18a). Increase in convective activity (as well as convective outflow) over the Indian region associated with the northward migration of ITCZ along with the development of deep convection over the Bay of Bengal (particularly over its northern parts) during this period aid the formation of abundant dense STCs in the UT region. The annual variation of high altitude cloud amount over the Indian region also shows the presence of large deep convective cloud systems reaching very high altitudes and even penetrating the tropopause has been reported . As the strength of the tropospheric convection decreases significantly during winter and the ITCZ shifts to the southern hemisphere sector, the probability of occurrence of dense STC over the southern hemisphere increases significantly. The deep convective clouds over this region get confined to a small geographical region between equator and 10°. Most of the STCs occurring over the Indian region during this period will be of
The latitude variation of the annual pattern of τP in the LS1 and LS2 regions are presented in Figure 18b and Figure 18c respectively. In general, the mean τP in the LS1 (18-21 km) region is in the range of 0.0008-0.003, which is one order in magnitude less than that in UT. In this altitude region, τP clearly shows relatively high values in the off-equatorial (north of 15°) regions and low values in the equatorial regions. Studies  on stratospheric aerosol optical depth during the decay phase of the volcanic aerosol at mid-latitudes using SAGE-II data have shown the presence of sinusoidal variation in the aerosol optical depth superimposed on an exponential decay with maximum and minimum occurring during the local winter and summer respectively. Over tropics the amplitude of these oscillations are significantly small and hence gets submerged in the disturbances caused by subsequent minor volcanic eruptions. As the period selected for the present analysis is volcanically quiescent, these oscillations are well discernable over the tropics. In the region 10-15°N, the annual variation of τP shows a relatively high value during winter and low value during summer. Note that, this variation is very similar to the annual variation of particulate backscatter observed in the lidar data from Gadanki . Above 15°N, high values of τP are observed in the LS1 region during the May-August period centered around 20-25°N. The particulate extinction just above the tropopause also shows a significant enhancement in the latitude region between 20-25°N (Figure 18b). This could be due to the penetration of the top of the high altitude semitransparent cirrus clouds above the cold point (tropopause). The two peaks observed in CAPE (and TSD) near equator (Figure 17) during April and October becomes more prominent and gets closer with increase in latitude (towards north) and merges to become a strong broad peak centered around July. This shows that the convective activity in the 20-30°N is very strong and the outflow occurs very close to the tropopause. The frequency of occurrence of cloud top altitude (observation from CALIPSO) shows a maximum value of 17±0.5 km in this latitude region during the June to September period which is at least 1 km larger than the maximum value observed at other latitudes (between 30°S to 30°N) at any period . Thus convection over the Indian land mass during the summer monsoon period in the 20-30°N latitude band is the strongest one in the entire latitude region 30°S to 30°N at any time during the year. This is a characteristic feature for the Indian longitude region. Strong convection plays a major role in transporting particulates from the upper tropospheric cirrus cluster to lower stratosphere causing a pronounced increase in τP in the LS1 region. Note that, such a feature in convection and hence in τP (in the LS1 region) is not observed in the southern hemisphere. In LS2, τP shows a distinct latitude variation with relatively high values near the equator up to 15° in both the hemispheres and low values over the off-equatorial regions (>15°). Relatively high values of τP are observed in the LS2 region during the January-June period in the latitude region 15°S-10°N and low values during the rest of the period. This pattern slowly reverses with increase in latitude. Beyond 10°N, τP shows a pronounced winter peak with low values during the March-May period. The annual variation of τP in the LS2 region over the southern hemisphere is quite similar to that in the northern hemisphere. High values of τP observed between 15°S and 15°N in the LS2 region confirms the presence of a Tropical Stratospheric aerosol Reservoir (TSR) during the study period. Earlier studies carried out by several investigators [29,83] revealed the presence of this band structure (with high aerosol loading) in the 21- 30 km altitude region over the equatorial region during volcanically perturbed period. Examining the aerosol climatology in the LS region using the SAGE-II data both during the volcanically perturbed period as well as during the near background conditions,
13. Influence of moderate volcanic eruptions on αp in the LS region
The background stratospheric aerosol layer usually referred to as the Junge layer  consists of liquid droplets composed of a mixture of sulfuric acid and water. This layer will be quite prominent subsequent to major volcanic eruptions, such as El Chichón (Mexico, 1982) and Mount Pinatubo (Philippines, 1991), which are powerful enough to inject large amount of SO2 into the stratosphere . After oxidation, sulfate aerosols are formed at these altitudes. These particles are subsequently distributed globally depending on the latitude of the eruption. Removal of these aerosols is rather difficult. It takes several months to years (depending on size) to scavenge these volcanic aerosols. The last major eruption (Mount Pinatubo) took place in 1991 and the stratospheric aerosol layer returned to its “background” level around 1997. There is no major increase in stratospheric aerosol loading after 1997 [13,25]. However, increase in the anthropogenic SO2 emission has been proposed as a plausible mechanism responsible for the observed small increasing trend in stratospheric background level in the recent past . An overview of the current understanding on stratospheric aerosol science can be found in
Though, in general, most of the particles generated through the gas-to-particle conversion process in the stratosphere will be small and nearly spherical in nature, there could be a few larger size particles in the lower stratosphere associated with moderate and intense volcanic eruptions, leading to a pronounced enhancement in δ and αp in this region. These volcanic perturbations are clearly distinguishable from STCs in the UT region based on the amount of enhancement, its temporal structure as well as the duration of enhancements. Moreover, while δ (and αp) of STCs vary significantly at shorter time scales [57,87], the stratospheric cloud formed through volcanic emissions will be stable for a longer period. In addition, the values of δ associated with STC will be significantly larger (as they are mostly composed of non-spherical ice crystals) than those of volcanic clouds. The vertical structure of δ for the volcanic cloud will remain fairly stable at short time scales typically over a night. Thus a long lasting enhancement of δ in the LS region observed by the lidar is an indicator for assessing the volcanic impact on stratospheric aerosols. In order to illustrate this in detail and depict the difference in the nature of perturbations due to STC and volcanic cloud, a contour map of δ for a typical night (25 November 2002) during the eruption period of Reventador (started in November 2002 and lasted up to January 2003) is presented in Figure 19a. In order to accommodate the large variations, the contouring of δ in this figure is performed in two bands; one from 0.04 to 0.2 in steps of 0.02 and the other above 0.2 at 0.2 interval. Two enhanced layers (of δ) one between 15 and 16 km and the other around 19km, are distinctly seen in this figure. The value of δ and its temporal variations are very large for the lower layer while these are very small for the upper layer. The lower layer disappears after mid-night while the upper layer continues to persist up to the end of lidar observations. This shows that the lower layer is an STC while the upper one is the volcanic cloud. This volcanic cloud will persist on subsequent night also while STC may or may not be present. To illustrate this the sequences of lidar profiles observed on a few nights during the eruption period of Reventador are presented in Figure 19b.The features of the layer located around 19 km (in the LS region) is rather steady in all these profiles while those of the layer below ~17 km (due to STC), is highly variable. The stratospheric layer is strongest in December and started decaying in January 2003 and became almost insignificant by March 2003, while the STC layers appear at random.
The perturbation in the lower stratosphere over the equatorial and off-equatorial regions during the period 1998-2005, when the volcanic activity was relatively small, are examined in detail  using the altitude profiles of αp over the tropics from SAGE-II data archive (Figure 20a) and the altitude profiles of αp and δ obtained from lidar data (Figure 20b and 20c) at Gadanki. Zonal-averaged monthly mean αp at different altitudes in the lower
stratosphere for four different latitude belts (over tropics) are used mainly to have sufficient number of profiles in each belt as well as owing to the fact that in the lower stratosphere the spatial variability along longitude could be minimal (because of efficient mixing in the zonal direction and strong horizontal transport prevailing in the region). As expected the aerosol extinction decreases with increase in altitude with significant loading confined to the altitude region 18–27 km. Short-lived enhancement in αp are distinctly seen in the lower stratosphere (Figure 20a) in different years. These signatures could be attributed to the influence of various minor volcanic eruptions during the study period. Those eruptions which could make discernable impact in the stratosphere are marked at the start of the respective eruption along the X-axis of each panel and further details of these eruptions are summarized in Table 1. Only those eruptions which occurred between 30°S and 30°N are included in Table 1. The Volcanic Explosivity Index (VEI) for different eruptions  is obtained from the Web site of Smithsonian Institution–Global Volcanism Program (GVP). The eruptions of Ulawun (eruption d), Ruang (eruption e), Reventador (eruption f), and Manam (eruption j) are relatively strong with VEI ~4 (Table 1), and the signature of these eruptions are well discernable in the lower stratospheric (18–20 km) aerosol extinction (Figure 20a). These perturbations can be distinctly seen both over the equatorial as well as off-equatorial regions.
Figure 20b shows the temporal variation of the altitude structure of αp obtained from lidar during the period 1998– 2005 in the form of a contour plot. This figure also depicts the signatures of various minor volcanic eruptions similar to that depicted in the zonal mean values obtained from SAGE-II data, in the latitude sector 0–15°N (Figure 20a). Small-scale features are more pronounced in lidar data (Figure 20b) mainly because of the fact that it corresponds to a point observation while that in Figure20a is the zonal average. In volcanically quiescent periods the sulfur bearing gases, SO2 and OCS emitted from the earth’s surface are transported across the tropopause , photolyzed, and oxidized to sulfuric acid before condensation to form sulfuric acid and water droplets in the stratosphere [24,91], which are, by their liquid nature, spherical. The depolarization caused by these particles will be very small and hence the resultant δ in the stratosphere will be very close to that of the molecules. In general, it is in the range 0.03 to 0.04. But, during major volcanic eruption abundant amount of precursor gases will be injected into the lower stratosphere along with a few fine particulates. Because of this influx of particles and gases there will be an increase in the number density of particles as well as an increase in the size of these particles. The size spectrum of the stratospheric particles also shifts toward the larger size regime following the volcanic eruption . As some of those particles that are directly injected into the lower stratosphere during volcanic eruption could be non-spherical in nature, an increase in δ in the stratospheric aerosol would be expected.
Figure 20c shows a contour plot of δ at different altitude in the lower stratosphere over Gadanki during the period 1998–2005. This plot generated adopting the same procedure as that used for generating Figure 20b from αp profiles, clearly shows a few short-lived δ enhancements in the lower stratosphere. The sporadic increase in δ is associated with the eruption of a few moderately intense volcanic eruptions. The duration of the increase in αp also matches well that of δ. The disturbances caused by eruption of volcanoes Ruang (eruption e) and Reventador (eruption f) are relatively stronger (δ ranging from 0.05 to 0.2 in the altitude region 18–21 km during November 2002 to February 2003). In a few cases the enhancements in αp and δ does not match exactly. This could be due to the fact that the volcanic locations are at different distance from Gadanki as well as the prevailing transport process could be different at various occasions.
14. Long term variations of τp in the LS region
Zonal mean values of τp in the altitude region 18–28 km averaged for the four latitude belts, each of width 15°, in each month are used to examine the temporal variations in both the hemispheres. Time series plots of these for the equatorial (0–15°S and 0–15°N) and off-equatorial regions (15–30°N and 15–30°S) are shown in Figures 21a and 21b, respectively, along with a similar plot of τp obtained from lidar at Gadanki in Figure 21c. In the equatorial region the temporal variations in τp are very similar in both the hemispheres. The mean level of τp shows an abrupt increase after 2002  though the oscillations around this mean level is fairly similar to those before 2002. Over the off-equatorial regions the mean level of τp shows a rather steady (gradual) increase from 1998 to 2005. A sharp increase in τp is observed towards the trailing edge of the data in the year 2005. During the period 1998–2002, the value of τp in the equatorial region is a minimum and is close to ~0.0025. The mean τp in the altitude region 18–28 km obtained from lidar at Gadanki also shows similar feature as that observed from SAGE-II over the equatorial region. In general the value of τp in the LS region obtained from lidar data is larger than that obtained from SAGE-II. This could be due to the fact that the lidar observation is a single point measurement while the SAGE-II data used in this analysis is zonally and meridionally averaged for the equatorial and off-equatorial regions.
All the plots in Figure 21 show a general increasing trend in stratospheric particulate (aerosol) optical depth during the period 1998–2005 in addition to the periodic variations. On the basis of lidar observations during the period 2000–2009, a similar increase in the integrated stratospheric backscatter coefficient (in the altitude region 20–25 km) at the rate
of 4.8% and 6.3% per year (with respect to its value in 2002) was reported  at Hawaii (19.5°N, 155.6°W) (Mauna Loa Laboratory) as well as at Boulder (40°N) (Colorado). The increase in τp at Hawaii and Boulder was attributed to the increase in global coal consumption since 2002, mainly from China, and subsequent increase in emission of SO2. Anthropogenic aerosols produced through gas-to-particle conversion of precursor gases like sulfates and ammonia transported to the upper troposphere  through intense convection  in the tropics and subsequently across the tropopause also could possibly be a contributing factor for this increase. An increase in tropical upwelling (Brewer- Dobson circulation) because of global warming also was suggested to be a plausible mechanism for the observed increasing trend [93,94] in stratospheric βp after 2002.
Even though on an average the stratospheric particulate loading is in its background level during the period 1998–2005, it was influenced particularly by a few moderate volcanic eruptions mainly after September 2002. While the period before September 2002 was absolutely quiet (with low particulate loading), the later period was mildly disturbed. The variation in stratospheric particulate loading need not solely be represented by a corresponding variation in tephra emissions Their could be some other causative mechanisms, such as increase in anthropogenic emissions as well as the increase in tropical upwelling, which could influence the stratospheric particulate loading.
15. Periodic variation of τp in the LS region
In addition to the general increasing trend, the values of τp in Figure 21 shows a seasonal cycle with winter maximum and summer minimum modulated by a long-period oscillation. These oscillations could primarily be due to the influence of large-scale atmospheric waves. The time series data of τp is spectrum analyzed to bring out the characteristics of the prevailing periodic variations. Before subjecting the data to spectral analysis, the linear trend is removed from the original data. The residual part is Fourier analyzed and the resulting amplitude spectra for different latitude bands (0–15°N, 0–15°S, 15–30°N and 15–30°S) are presented in Figures 22a and 22b. These amplitude spectra reveal the presence of a strong annual component (~12 months) along with a quasi-biennial component (~30 month) both in the equatorial and off-equatorial regions . The spectral amplitude of QBO is as strong (significant) as that of annual oscillation (AO). Figure 22c shows the amplitude spectrum obtained from the lidar derived values of τp (at Gadanki). Even though Figure 22c shows more significant peaks in the short-period regime, the spectral amplitude is more pronounced for semi-annual (SAO) and annual (AO) components. This spectrum also shows a secondary peak around 46 months followed by troughs at 23 and 92 months. The spectral amplitude for 30 month periodicity is larger than those at the troughs on either side of this secondary peak. Though the period for the peak amplitude (46 months) is much larger than that expected for the stratospheric QBO, on the basis of the inference derived from Figures 22a and 22b, as well as owing to the fact that the spectral amplitude at 30 months is not a minimum, the characteristics of the 30 month periodicity is examined in the later part to delineate its altitude structure in the lidar data. The SAGE-II data did not show the signature of SAO, which could be due to the inherent smoothing out of this component while taking the zonal mean. The higher spectral amplitude of AO compared to that of QBO is quite expected (as it is true for the wind field also).
16. Quasi-biennial oscillations in τp and zonal wind in the LS region
For this study the high resolution Radiosonde data of zonal wind in the lower stratosphere from an equatorial station (where the quasi-biennial oscillation, QBO, signature is expected to be maximum), Singapore (1°22″N, 103°55″E), obtained from Web site http://www.geo.fu-berlin.de/en/met/ag/strat/produkte/qbo/ are used. The time series data of monthly mean zonal wind at Singapore at 20 hPa and 30 hPa levels during the study period are presented in Figure 23. This figure shows that the quasi-biennial oscillation in zonal wind (QBOU) is in the easterly phase during the 1998, 2000–2001, 2003 and 2005. The westerly phase during 1999 and the easterly phase during 2000–2001 are relatively broad. Different phases of QBOU are indentified from this time series to study the influence of QBOU in lower stratospheric aerosols.
16.1. Latitude variation of τp in LS region over the Indian longitude sector in two different phases of QBOU
The latitude variation of τp in the altitude range 21–28 km in the band 0–30°N (averaged for every 5°) over the Indian longitude sector (70–90°E) is examined during the consecutive easterly and westerly phases of QBOU in 1998 and 1999 separately. Figure 24 shows the latitudinal variation of τp during these two phases of QBOU. This mean τp is obtained by averaging the particulate optical depth in individual months when the QBOU phase has reversed completely (the wind speed has reached its highest value). While the value of τp is relatively high during the westerly phase of QBOU in the equatorial region up to 15°N, relatively high values are observed during the easterly phase of QBOU in latitudes beyond 15°N. The lidar derived τp is relatively low during the easterly phase of QBOU compared to that during its westerly phase. This is in good agreement with that from the latitude variation of τp for these 2 years derived from SAGE-II data.
The value of τp shows a pronounced decrease with latitude during westerly phase of QBOU compared to that during its easterly phase. Interestingly, during the easterly phase of QBOU, while the latitude variation of τp is relatively small in the equatorial region (0– 15°N), beyond 20°N it decreases sharply with increase in latitude. By examining the altitude profile of backscatter ratio obtained from lidar data at Mauna Loa at Hawaii (19.5°N, 155.6°W), Barnes and Hofmann  reported an enhancement of backscatter ratio (enhanced aerosol loading) in the altitude region 21–30 km during the easterly phase of QBOU. This feature agrees well with the above observation derived from the latitude variation of τp at Hawaii.
16.2. Latitude variation of zonal mean τp in the LS region over the tropics in different phases of QBOU
On the basis of the measured particulate loading in the stratosphere, out of the 8 year period considered for the present analysis, the period 1998–2002 was absolutely quiescent while the period 2003–2005 was mildly disturbed. The variation of zonal mean τp in the latitude region 30°S to 30°N (averaged for every 5°) in the alternative easterly and westerly phases of QBOU during the volcanically quiescent and mildly disturbed periods are presented in Figures 25a and 25b respectively. As seen from Figure 23 during the study period (8 years) the QBOU completes approximately 3.5 cycles, with two easterly (E1 and E2) and two westerly (W1 and W2) phases during the period 1998–2002 and two easterly (E3 and E4) and one westerly (W3) phase during the mildly disturbed period of 2003–2005.
During the former half (very quiet period), the mean τp in the equatorial region shows a general enhancement during westerly phase of QBOU (W1 and W2) while in the off-equatorial region of northern hemisphere it shows an enhancement during the easterly phase of QBOU . When the stratospheric QBOU was in its westerly phase the mean values of τp in both the hemispheres decreases steeply with increase in latitude on either side of equator. During the easterly phase of QBOU, though the value of τp decreases with increase in latitude beyond 15°N, it is fairly uniform in the equatorial region (15°S to 15°N) with a bite-out over the equator. During the latter half (mildly disturbed period), the mean τp decreases rapidly on either side of equator with increase in latitude in both the phases of QBOU. Figure 25b shows a general enhancement in τp during W3 period in the equatorial and off-equatorial regions of both the hemispheres compared to that during the E3 period. However, during E4, τp in general is relatively large in the equatorial and northern hemispheric off-equatorial region compared to that during E3 and W3. Note that, in the above analysis two consecutive easterly and westerly cycles of QBOU are considered and the inferences arrived based on the latitudinal structure of τp in these two phases. The inference may apparently be contradicting if one considers the other two pairs, westerly of the first pair and easterly of the next (eg W3 and E4), particularly during the mildly disturbed phase of 2003–2005 (Figure 25b). This is mainly due to the fact that the background τp has a strong increasing trend during this half in addition to the periodic variations. This can override the periodic variations associated with QBOU. However, this contradiction is totally absent in the former half (1998–2002) when the background τp was fairly steady. During this period τp in the equatorial region is high in both the westerly phases (W1 and W2) compared to that during the easterly phases E1 and E2, irrespective of sequence of pairs considered. This shows that the influence of QBOU in τp can be clearly delineated only during the very quiescent volcanic periods, when the background stratospheric aerosol loading remains in its steady background level.
16.3. Altitude structure of αp in the lower stratosphere in different phases of QBOU
Figure 26a shows the mean profile of αp in the altitude region 18–28 km obtained from lidar data at Gadanki for the two periods, January to September 1998 and November 1998 to September 1999, when the QBOU was in its easterly and westerly phases, respectively. Figure 26a clearly shows that the value of αp is consistently larger during the westerly phase of QBOU than the corresponding values during its easterly phase. This shift in the altitude profile toward a higher value side during the westerly phase is not due to the influence of any trend because the mean level of the stratospheric particulate optical depth in the 0°–15°N region remains fairly the same during the period 1998–2002 and an increase in this level occurred only after 2002 (Figure 21).
For a more detailed study on the effect of QBOU on stratospheric aerosols, the zonal mean altitude profile of αp in different phases of QBOU during the study period is examined for the equatorial and off -equatorial regions separately. The mean profiles thus obtained for the latitude region 0–10°N and 20–30°N (averaged zonally, as well as along the respective latitude bands) are shown in Figures 26b and 26c, respectively. Over the equatorial regions the altitude profile of aerosol extinction coefficient in the altitude region 22–27 km during the easterly phase of QBOU are consistently lower than that during the subsequent westerly phase. In the off-equatorial region, the aerosol extinction coefficient in the altitude region 25–32 km is found to be relatively low during the westerly phase. It may be noted that the earlier studies also have shown that aerosol extinction in the lower stratosphere can be influenced by the phase of the quasi biennial oscillation [11,28,29,95] in this region even though the difference in the nature of the latitudinal dependence over the equatorial and off-equatorial regions were not addressed in detail.
Association between the altitude structure of QBO in particulate extinction (QBOa) and QBOU for the entire study period is examined by subjecting the time series of zonal mean monthly average αP (obtained from SAGE-II) and zonal mean zonal wind, U (obtained from NCEP reanalysis), at each altitude bin from 20 to 30 km to Fourier analysis after removing the linear trend. To illustrate the QBO features, the spectral components corresponding to the periods other than 22-45 month are removed from the respective trend removed time series data using the corresponding amplitudes and phases. The altitude structure of the residual amplitudes thus obtained for αp and U over the equatorial and off-equatorial regions are presented in Figure 27. The left side panels show the residual amplitude for αp and the right side panels the same for zonal wind. Both these panels show clearly the presence of a QBO. While the phase of the biennial oscillation in wind shows a clear downward propagation with time as expected, the phase of QBOa changes many times with increase in altitude. This change in the phase of QBOa with altitude is analogous to that reported for the stratospheric ozone  and was attributed to the secondary meridional circulation (SMC).
17. QBO in αp and secondary meridional circulation
To study the altitude structure of the QBOa, the time series of zonal-averaged monthly mean aerosol extinction coefficient at different altitudes (at 0.5 km interval) are subjected to Fourier analysis after removing the linear trend. This analysis of αp at different levels showed a prominent peak in its amplitude corresponding to the periodicity of QBO (30 month) and AO (12 month). Figure 28 shows the altitude profile of the amplitude and phase of QBOa from 18 to 32 km in the equatorial (0–15°N and 0–15°S) and off-equatorial (15–30°N and 15–30°S) regions. In general, the amplitude of QBOa decreases with increase in altitude (as that observed in the raw aerosol extinction data). Over the equatorial region, the altitude variation of QBOa amplitude (Figures 28b and 28c) shows three prominent peaks in the region 18–22 km, 23– 27 km and 28–32 km with peak amplitude centered around 20 km, 25 km and 30 km, respectively. This feature is remarkably symmetric about the. equator, in both the hemispheres. Another interesting feature to be noted is that the phase remains fairly constant around these peaks. In the equatorial region, the phase of QBOa around 25 km is ~15 months and that at around 30 km and 20 km are ~28 and 0 months, respectively. Thus the phase difference of QBOa in the upper and lower regime with respect to that at 25 km is around 13 and 15 months, which corresponds to 156° and 180°, respectively (for the 30 month periodicity one month in phase corresponds to 12°). This shows that the QBOa around 25 km is almost out of phase with that in the upper (28–32 km) and lower (18–22 km) regime .
Over the off-equatorial region the altitude structure of the amplitude and phase of QBOa are found to be different from that over the equatorial region. Figures 28a and 28d show that the amplitude of QBOa is significant in the lower and upper altitude regions with a minimum value around 24.5 km in the northern hemisphere and around 21 km in the southern hemisphere. While the amplitude of QBOa shows a broad maximum in the altitude region 24–28 km with a peak around 26 km in the southern part, correspondingly it shows a maximum around 29 km in the northern part.
Over the equatorial region, the observed features of QBOa in the three altitude regions could be attributed to the influence of the secondary meridional circulation (SMC) induced by the vertical shear of QBO in stratospheric zonal wind. Through a detailed analysis of the QBO modulation of the meridional wind in the stratosphere, Ribera et al.  demonstrated the existence of discrete zones of meridional wind convergence and divergence over the equator. These convergence (divergence) zones during the westerly (easterly) phase of QBOU are located at the lower and upper limits of the maximum zonal wind shear and maximum temperature anomaly layers. The existence of convergence and divergence zones, which are directly related to the rising and sinking motions, forms two circulation cells (SMC) in the consecutive vertical levels, quasi-symmetric about the equator.
To examine the cycle-to-cycle variation of the QBO in aerosol extinction coefficient and its phase structure at the three altitude regions, the monthly mean values of αp at the central altitudes (of the three regions) 20 km, 25 km and 30 km for the equatorial region 0–15°N and 0–15°S are estimated during the period 1998 to 2005 and this data is subjected to wavelet analysis (after de-trending). Contour plots in Figure 29 shows the time history of the amplitudes for different periodicities (ranging from 2 to 45 months) in αp for the study period. The spectral characteristics at the respective altitudes are similar in both the hemispheres. The annual component repeats coherently in all the three altitudes indicating that they follow the same pattern in the entire altitude range 18–32 km. But the temporal pattern of QBO differs. When the QBO in αp is in its positive phase (maximum amplitude) in the lower (20 km) and upper (30 km) regimes, it is in the negative phase (minimum amplitude) in the middle regime (25 km). The QBOa at 25 km is out of phase with that in the upper and lower regimes during the entire study period. However, the amplitudes of AO and QBO in αp at 30 km are relatively small in the first 40 months compared to rest of the period. This analysis confirms the temporal repeatability of QBOa phase structure in these three altitude regions inferred from the mean αp at different altitudes averaged for the entire period of analysis shown in Figure 28. The phase of QBOU (Figure 23) matches with that of QBOa in the middle regime (25 km) while it is in opposite phase with those in the lower (20 km) and upper (30 km) regime.
Studies by Fadnavis and Beig  on the influence of SMC induced by QBO in the spatiotemporal variations of ozone over the tropical and subtropical regions however showed maximum amplitude of the QBO signal in ozone concentration at the two pressure levels, 30 hPa and 9 hPa, where the QBO manifests in opposite phase. This feature is in good agreement with the amplitude and phase of the QBO signal observed in aerosol extinction coefficient at 25 km and 30 km. From Figures 28b and 28c it can be seen that while the maximum amplitude of QBO signal in aerosol extinction coefficient centered around 20 and 30 km coincides with SMC cells centered at ~80 hPa and ~9 hPa as illustrated by Plumb and Bell  and Punge et al. , the maximum QBO signal in aerosol extinction coefficient in the altitude region 23–27 km centered around 25 km coincides with the level of maximum zonal wind shear (centered ~24 hPa) where the divergence and convergence zones (circulation cells) intersect. Similarly, Dunkerton  studied the QBO anomalies in ozone, methane and water vapor and observed a double peak structure in their amplitudes. Numerical simulations of QBO anomalies also showed influence of SMC in modulating the distribution of tracers leading to the formation of two peaks in their number density near 24 and 32 km [101-103]. However, the observed maximum amplitude of QBO in αp centered around 20 km is not reproduced in these model simulations.
The altitude structure of aerosol extinction in the tropical UTLS region and its variability in different time scales is examined with particular stress for the Indian longitude sector, using the data from a dual polarization lidar (at 532 nm wavelength) located at Gadanki along with data obtained from SAGE-II onboard ERBS. Organized convection associated with the Asian summer monsoon (ASM) and highly dynamic ITCZ makes this region unique in the global scenario. Altitude profiles of particulate extinction for a period of eight years from 1998-2005 (volcanically quiescent period) obtained from SAGE-II in the latitude region 30°S to 30°N and the altitude profiles of particulate backscatter, particulate extinction and volume depolarization for the same period derived from lidar data at Gadanki are used for this purpose. The lidar data are used mainly to study the influence of local features such as STCs, local convection etc in the properties of particulates in the UTLS region while the satellite data are used to study the spatial features in the tropical UTLS from 30°S to 30°N. A comparison of aerosol extinction from SAGE-II with that from lidar on a profile basis as well as on monthly mean basis over a small geographical grid size (6° in latitude and 13° in longitude) centered at Gadanki showed a very good agreement between the two in the major features such as mean altitude structure as well as the mean annual pattern.
Lidar studies show a significant increase is particulate scattering in the upper troposphere with a high value of δ over the Indian subcontinent during the ASM period when the tropospheric convective activity is the highest. This feature is clearly associated with particle formation due to homogenous/heterogeneous nucleation and condensation and subsequent freezing to form non-spherical ice crystals (associated with STC). The particulate backscatter coefficient in the LS region is maximum during winter and minimum during the ASM period. The winter peak is closely associated with the increase in vertical mass flux in conjunction with the tropical upwelling and Brewer Dobson circulation.
Occurrence of STCs in the UT region is very common in this geographical region. The frequency of occurrence of STC is the largest during the ASM period, when most of the observed STCs are optically dense and geometrically thick. The most favoured altitude for STC is 14-16 km. Most of the observed STCs are optically and geometrically thin. Thin STCs are generally observed at higher altitudes, very close to the cold point. Most of the dense STCs are found to be associated with organized convection while the thin STCs could be of
The convective activity prevailing in the troposphere significantly influences the microphysical properties of particulates in the UTLS region. The particulate scattering and optical depth in the UT region shows a general decrease with increase in latitude on either side of the equator with a well pronounced summer-winter contrast. However, organized convection during the ASM period enhances the particulate loading in the UT in the northern latitudes even beyond 25°N. While the particulate optical depth in the 18–21 km region (lowest part of the stratosphere) is relatively low in the equatorial region, it shows an increase in the off-equatorial region mainly due to this enhancement in particulate concentration above the cold point, particularly over the Indo-Gangetic Plain, during this period. At a higher altitude (21–30 km) it shows a different pattern, with high values near the equator and low values in the off-equatorial region. This confirms the existence of a stratospheric aerosol reservoir. This spatial distribution could be attributed to horizontal advection in the lower regime (rapid transport from near equatorial region to higher latitudes) as well as lofting to higher altitudes over the equatorial region (B-D circulation).
Spectral analysis of zonal mean particulate optical depth in the stratosphere (18-32km) revealed the existence of a strong QBO both in the equatorial and off-equatorial regions. The phase of the QBO signal in particulate extinction (QBOa) around 25 km is found to be in opposite phase with that in the upper (28-32 km) and lower regime (18-22 km), illustrating the existence of a secondary meridional circulation (SMC) produced due to vertical shear of QBO phase in zonal wind (QBOU). While the particulate optical depth in the lower stratosphere is relatively large during the westerly phase of QBOU in the equatorial region, relatively high values are observed during the easterly phase of QBOU in the off-equatorial region. During the westerly phase of stratospheric QBOU, the mean particulate optical depth rapidly decreases with increase in latitude on either side of equator in both the hemispheres. During the easterly phase, this remains fairly steady between ±15° latitude, with a small bite-out around the equator, and decreases steadily for latitudes beyond 15°.
We thank the NASA Langley Research Center (NASA-LaRC) and the NASA Langley Aerosol Research Branch for providing the SAGE-II data through the web site ftp://ftp-rab.larc.nasa.gov/pub/sage2/v6.20. The authors are thankful to the technical and scientific staff of National Atmospheric Research Laboratory (NARL), Gadanki for their dedicated efforts in conducting the Lidar observations. One of the authors, Dr. K. Parameswaran would like to acknowledge CSIR for providing grant through the Emeritus Scientist scheme.