Earthquake responses in the high-latitude ionosphere

The data, obtained using the methods of partial reflections and ionosphere vertical sounding on the Kola Peninsula and in Scandinavia, at Tumannyi (69.0° N, 35.7° E) and Sodankyla (67.37°N, 26.63°E) observatories, have been analyzed in order to detect earthquake responses. The strong earthquakes have been considered: one earthquake with a magnitude of 7.7 occurred at 0819:25 UT on July 17, 2006, on the western coast of Indonesia (9.33° S, 107.26° E), and another earthquake with a magnitude of 6.2 occurred 2253:59 UT on May 26, 2006, on Yava (7.94° S, 110.32° E). These earthquakes, the epicenters of which were located in the same region and at identical depths (10 km), were observed under quiet conditions in the geomagnetic field (ΣK _p = 5.7 and 6.3) and during small solar flares. The response of the ionosphere to these flares was mainly observed in the parameters of the lower ionosphere in the D and E regions. It has been found out that the period of variations in the ordinary component of the partially reflected signal at altitudes of the E region increased before the earthquake that occurred on July 17, 2006. The f _min variations at Sodankyla observatory started 20 h before the earthquake. The periods of these variations were 3–6 h. The same periods were found in the variations in other ionospheric parameters (foEs and h’Es). The variations in the ordinary component of partially reflected signals with periods of 2–5 hours were observed on the day of another earthquake (May 26, 2006). Internal gravity waves with periods of several hours, which can be related to the earthquakes, were detected in the amplitude spectra of the ordinary component of partially reflected signals and in other parameters in the lower ionosphere.     

The influence of geomagnetic and solar variabilities on lower thermosphere density

Atmospheric density measurements near 200 km from the Satellite Electrostatic Triaxial Accelerometer (SETA) experiment are analyzed for geomagnetic and solar flux variability effects. Data from the SETA experiment, onboard two satellites, are available for the periods of May to November 1982, and July 1983 to March 1984. The data utilized the span ±79.5° latitude, and are available for both day (1030 LT) and night (2230 LT). Annual and semiannual density variations are removed and regression analyses are performed on the residuals using a series of lagged 3 h K_p indices to determine and remove geomagnetic fluctuations. Densities are found to increase by as much as 134% in response to an increase in the K_p index from 1 to 6. Monthly curves are generated for the K_p regression coefficients to delineate seasonal-latitudinal and day/night dependences, which reflect the effects of mean meridional advection of disturbances from high to low latitudes. Further analyses are performed comparing measured densities with MSISE-90 predictions. Results show that the model is able to capture many of the prominent features, but does not fully predict the level of variability for the individual disturbance periods analyzed. After the geomagnetic effects are removed, the residual densities are interpreted in terms of solar flux variability. The daily-averaged SETA density residuals are strongly correlated with long-term solar flux variability, and exhibit a much greater dependence on the 27-day solar rotation period than MSISE-90 predictions. Variations in residual density of the order of 10–20% occur in association with day-to-day and 27-day solar flux variations. The MSIS model does not accurately predict the magnitude of these short-term density variations in response to solar activity.     

Analogue model,relating Kp index to solar wind parameters

In a previous work the authors have developed a model, providing K_p as a function of the interplanetary magnetic field B_z component. They introduced a modified B_z function (denoted as B_zm), exhibiting a delayed reaction to B_z changes. The modified function B_zm was defined by using the analogy with a damping RC-circuit output voltage. The delaying reaction of B_zm to B_z was characterized by two time constants, one for rising and one for decreasing parts of B_z. The cross-correlation between K_p and B_zm has increased to 0.7, compared with −0.4 between K_p and B_z. In this paper, new dependences of K_p on solar wind velocity and dynamic pressure are included in the model to improve its accuracy. These solar wind parameters are found to correlate best with K_p. The hourly interpolated values are also added to the 3-h K_p values to increase the statistics. The new K_p data set is denoted as K_p1. The mean dependence of K_p on B_zm and dynamic pressure are approximated with parabolas, while the dependence on the velocity is linear. The constants in the model expression are obtained by using ACE data (1998–2000). The overall model error is estimated at 0.63 units K_p. The improvement over the previous simpler dependence in terms of the model error is about 30%.     

Tec climatology derived fromtopex/poseidonmeasurements

We have used the TOPEX data sets available from JPL to create a customizeddatabase of the TOPEX TEC measurements that contain data from 1992 through part of 1996.Thedata base includes time, geographic and geomagnetic coordinates of themeasurement,geomagnetic indices (Kp, previous K_p, HemisphericPower, andintegral of hemispheric power over the previous 36 h), solar index (F 10.7),andInternational Reference Ionosphere (IRI) results corresponding to the TOPEX measurements.Inthis paper we present global maps of TEC for low solar activity conditions (F 10.7 ⩽ 120) for quiet (integral of hemispheric power less than 800 GWh, roughly correspondingto K_p = 2), moderately disturbed (integral of hemispheric power greaterthan 800GWh but less than 1200 GWh, roughly corresponding to K_p = 3),anddisturbed (integral of hemispheric power greater than 1200 GWh, roughly corresponding to K_p = 4) geomagnetic conditions, derived by binning all appropriateTOPEX TECdata from 1992 to 1996. The analysis is performed in a Magnetic-Local-Time,Magnetic-Latitudecoordinate system. The most prominent feature of the global TEC maps is thefeaturecorresponding to the equatorial anomaly. The feature becomes wider in magnetic latitudeandmore pronounced in amplitude as the activity level increases. The equatorward shift of thecrests,with increased magnetic activity, can produce apparent decreases in TEC at their quiettimelocation for individual storms as evident in the conflicting conclusions of TECgeomagneticdependence studies of the 1960s. For the same activity level, TEC values in theequatorialanomaly are higher during equinox compared to solstice.     

Long term ionospheric electron content variations over Delhi

Ionospheric electron content (IEC) observed at Delhi (geographic co-ordinates: 28.63°N, 77.22°E; geomagnetic co-ordinates: 19.08°N, 148.91E; dip Latitude 24.8°N), India, for the period 1975/80 and 1986/89 belonging to an ascending phase of solar activity during first halves of solar cycles 21 and 22 respectively have been used to study the diurnal, seasonal, solar and magnetic activity variations. The diurnal variation of seasonal mean of IEC on quiet days shows a secondary peak comparable to the daytime peak in equinox and winter in high solar activity. IEC_max (daytime maximum value of IEC, one per day) shows winter anomaly only during high solar activity at Delhi. Further, IEC_max shows positive correlation with F_10.7 up to about 200 flux units at equinox and 240 units both in winter and summer; for greater F_10.7 values, IEC_max is substantially constant in all the seasons. IECmax and magnetic activity (A_p) are found to be positively correlated in summer in high solar activity. Winter IEC_max shows positive correlation with A_p in low solar activity and negative correlation in high solar activity in both the solar cycles. In equinox IEC_max is independent of A_p in both solar cycles in low solar activity. A study of day-to-day variations in IEC_max shows single day and alternate day abnormalities, semi-annual and annual variations controlled by the equatorial electrojet strength, and 27-day periodicity attributable to the solar rotation.     

Variability of the 557-nm atmospheric emission

The variability of diurnal, day-to-day, and monthly average rates of the 557.7-nm atmospheric emission (I _557.7) is considered. We use 1997–2010 airglow observation data obtained for the upper atmosphere over Eastern Siberia (52°N, 103°E). The variation coefficient K_V of corresponding quantities is taken as the variability index. For the 23rd solar cycle, we examine the resulting seasonal variation of K_V of monthly averaged I _557.7, the dependence of monthly averaged I _557.7 on solar activity for each month, the variation coefficient of diurnal values of I _557.7 for different seasons of the year, the variability of I _557.7 during some geophysical events, and the correlation of I _557.7 variations with global climatic indices.     

Determining dynamic parameters of different-scale ionospheric irregularities over northern Siberia

In 1995–1996, observations were carried out at Norilsk (geomagnetic latitude and longitude 64.2°N and 160.4°E) to determine dynamic parameters of irregularities in the high-latitude ionosphere. The short-baseline spaced-receiver method that has been implemented at the ionospheric facility of the Norilsk Integrated Magnetic–Ionospheric Station, provides a means of simultaneously measuring parameters of small-scale irregularities (spatial scale of 3–5 km) by the Similar-Fading Method (SFM), as well as of medium-scale irregularities (time scale of 10–30 min, spatial scale of hundreds of kilometres) by the Statistical Angle-of-arrival and Doppler Method (SADM). About 20 h of the observational data for the F2-layer under quiet geomagnetic conditions (K_p < 3), 20 h under disturbed conditions (K_p ≥ 3) and about 15 h for the sporadic E-layer (K_p ≈ 3) were processed. It has been found that the propagation directions and velocities of different-scale irregularities do not coincide. Small-scale irregularities of the F2-layer travel predominantly eastward or westward. The velocity of the F2-layer irregularities is about 100 m/s, and under disturbed conditions it is up to 200–250 m/s. Small-scale irregularities of the sporadic E-layer travel mostly in the northward direction. It is confirmed that the E_s-layer is characterised by high velocities of the irregularities (as high as 1000 m/s). Medium-scale irregularities with periods in the range of 10–30 min travel mostly in a southward direction with velocities of 20–40 m/s.     

Biological accumulation of chlordane compounds in marine organisms from the northern North Pacific and Bering Sea

Sum of chlordane compounds (ΣCHL; cis-chlordane + trans-chlordane + cis-nonachlor + trans-nonachlor + oxychlordane) are concentrated gradually with trophic levels from zooplankton to Dall's porpoise (Phocoenoides dalli) through squid and fish. The order of bioconcentration factors (BCF: concentration in organism/concentration in seawater) in these organisms was ΣDDT (p,p′-DDE + p,p′-DDD + p,p′-DDT)>ΣCHL≦PCBs>ΣHCH (α-HCH + β-HCH + γ-HCH). Calculation for the concentration factor against food, namely biomagnification factor (BMF: concentration in organism/concentration in its food), was made for Dall's porpoise and thick-billed murre (Uria lomvia). The BMFs of these chemicals in thick-billed murre were lower than those of Dall's porpoise, suggesting the degradation and/or excretion of organochlorines through the uropygial gland with lipids. Moreover, the lowest BMF of ΣCHL in thick-billed murre among organochlorines may indicate that chlordane compounds (CHLs) are metabolized more rapidly by this seabird than Dall's porpoise.     

Stress-strain state of the lithosphere in the southern Baikal region and northern Mongolia from data on seismic moments of earthquakes

Investigation and understanding of the present-day geodynamic situation are of key importance for the elucidation of the laws and evolution of the seismic process in a seismically active region. In this work, seismic moments of nearly 26000 earthquakes with K _p ≥ 7 (M _LH ≥ 2) that occurred in the southern Baikal region and northern Mongolia (SBNM) (48°–54°N, 96°–108°E) from 1968 through 1994 are determined from amplitudes and periods of maximum displacements in transverse body waves. The resulting set of seismic moments is used for spatial-temporal analysis of the stress-strain state of the SBNM lithosphere. The stress fields of the Baikal rift and the India-Asia collision zone are supposed to interact in the region studied. Since the seismic moment of a tectonic earthquake depends on the type of motion in the source, seismic moments and focal mechanisms of earthquakes belonging to four long-term aftershock and swarm clusters of shocks in the Baikal region were used to “calibrate” average seismic moments in accordance with the source faulting type. The study showed that the stress-strain state of the SBNM lithosphere is spatially inhomogeneous and nonstationary. A space-time discrepancy is observed in the formation of faulting types in sources of weak (K _p = 7 and 8) and stronger (K _p ≥ 9) earthquakes. This discrepancy is interpreted in terms of rock fracture at various hierarchical levels of ruptures on differently oriented general, regional, and local faults. A gradual increase and an abrupt, nearly pulsed, decrease in the vertical component of the stress field S _v is a characteristic feature of time variations. The zones where the stress S _v prevails are localized at “singular points” of the lithosphere. Shocks of various energy classes in these zones are dominated by the normal-fault slip mechanism. For earthquakes with K _p = 9, the source faulting changes with depth from the strike-slip type to the normal-strike-slip and normal types, suggesting an increase in S _v . On the whole, the results of this study are well consistent with the synergism of open unstable dissipative systems and are usable for interpreting the main observable variations in the stress-strain state of the lithosphere in terms of spatiotemporal variations in the vertical component of the stress field S _v . This suggests the influence of rifting on the present-day geodynamic processes in the SBNM lithosphere.     

Distribution of ionospheric O<Superscript>+</Superscript> ion in synchronous altitude region

Based on satellite observation data, using dynamics equation, the ionospheric O⁺ ion’s distribution in the synchronous altitude region for different geomagnetic activity indexK _p is studied by theoretical modeling and numerical analyzing, and semi-empirical models for the O⁺ ion’s density and flux versus longitude in the synchronous altitude region for differentK _p are given. The main results show that in the synchronous altitude region: (i) The O⁺ ion’s density and flux in day-side are larger than those in nightside. (ii) With longitude changing, the higher the geomagnetic activity indexK _p is, the higher the O⁺ ion’s density and flux, and their variation amplitude will be. The O⁺ ion’s density and flux whenK _{p ≥} 6 will be about ten times as great as that whenK _p = 0. (iii) WhenK _p = 0 orK _{p ≥} 6, the O⁺ ion’s density reaches maximum at longitudes 120° and 240° respectively, and minimum in the magnetotail. WhenK _p = 3−5, the O⁺ ion’s density gets to maximum at longitude 0°, and minimum in the magnetotail. However, the O⁺ ion’s flux reaches maximum at longitude 120° and 240° respectively, and minimum in the magnetotail for anyK _p value.     

Solar activity dependence of ionospheric electron content and slab thickness using different solar indices

Ionospheric electron content (IEC) and slab thickness () data for the period 1977 to 1980 from Lunping (23.03°N; 121.90°E subionospheric) have been examined for their solar activity dependence. Local noontime monthly means as well as values for the 5 QQ days in a month have been examined separately with different solar indices, namely: solar EUV flux (170–190 Å),S _{10.7 cm} flux and sun spot number (SSN) on a seasonal basis. Both IEC and parameters exhibit better correlation with solar EUV andS _{10.7 cm} fluxes than with SSN for all seasons. IEC increases linearly with both EUV andS _{10.7 cm} flux whereas with SSN it shows a distinct nonlinear relationship during all seasons in both monthly mean and 5 QQ days' values. This study indicates that for correlating and predicting the variations (especially the medium term) in the ionospheric parameters, both EUV andS _{10.7 cm} fluxes have an advantage over SSN.     

Statistical Analysis and Modeling of the Local Ionospheric Critical Frequency: A Mid-latitude Single-Station Model for Use in Forecasting

The hourly values of the F-layer critical frequency from the ionospheric sounder in Dourbes (50.1°N, 4.6°E) during the time interval from 1957 to 2010, comprising five solar cycles, were analyzed for the effects of the solar activity. The hourly time series were reduced to hourly monthly medians which in turn were used for fitting a single station foF₂ monthly median model. Two functional approaches have been investigated: a statistical approach and a spectral approach. The solar flux F_10.7 is used to model the dependence of foF₂ on the solar activity and is incorporated into both models by a polynomial expression. The statistical model employs polynomial functions to fit the F-layer critical frequency while the spectral model is based on spectral decomposition of the measured data and offers a better physical interpretation of the fitting parameters. The daytime and nighttime foF₂ values calculated by both approaches are compared during high and low solar activity. In general, the statistical model has a slightly lower uncertainty at the expense of the larger number of fitting parameters. However, the spectral approach is superior for modeling the periodic effects and performs better when comparing the results for high and low solar activity. Comparison with the International Reference Ionosphere (IRI 2012) shows that both local models are better at describing the local values of the F-layer critical frequency.     

Variability of the ionosphere

Hourly foF2 data from over 100 ionosonde stations during 1967–89 are examined to quantify F-region ionospheric variability, and to assess to what degree the observed variability may be attributed to various sources, i.e., solar ionizing flux, meteorological influences, and changing solar wind conditions. Our findings are as follows. Under quiet geomagnetic conditions (K_p<1), the 1-σ (σ is the standard deviation) variability of N_max about the mean is approx. ±25–35% at ‘high frequencies’ (periods of a few hours to 1–2 days) and approx. ±15–20% at ‘low frequencies’ (periods approx. 2–30 days), at all latitudes. These values provide a reasonable average estimate of ionospheric variability mainly due to “meteorological influences” at these frequencies. Changes in N_max due to variations in solar photon flux, are, on the average, small in comparison at these frequencies. Under quiet conditions for high-frequency oscillations, N_max is most variable at anomaly peak latitudes. This may reflect the sensitivity of anomaly peak densities to day-to-day variations in F-region winds and electric fields driven by the E-region wind dynamo. Ionospheric variability increases with magnetic activity at all latitudes and for both low and high frequency ranges, and the slopes of all curves increase with latitude. Thus, the responsiveness of the ionosphere to increased magnetic activity increases as one progresses from lower to higher latitudes. For the 25% most disturbed conditions (K_p>4), the average 1-σ variability of N_max about the mean ranges from approx. ±35% (equator) to approx. ±45% (anomaly peak) to approx. ±55% (high-latitudes) for high frequencies, and from approx. ±25% (equator) to approx. ±45% (high-latitudes) at low frequencies. Some estimates are also provided on N_max variability connected with annual, semiannual and 11-year solar cycle variations.     

Meridional distribution of226Ra in the eastern Pacific along GEOSECS cruise tracks

We present the distribution of²²⁶Ra in eight vertical profiles from the eastern Pacific. The profiles are located along a meridional trend near 125°W, from 43°S to 29°N. Surface²²⁶Ra concentrations are about 7 dpm/100 kg, except for the two stations south of 30°S where the higher values are due to the Antarctic influence. Deep waters show a distinctive south-to-north increase in the²²⁶Ra content, from about 26 to 41 dpm/100 kg near the bottom. Unlike in the Atlantic and Antarctic Oceans, the effect of²²⁶Ra injection from bottom sediments is clearly discernible in the area. The presence of this primary²²⁶Ra can be traced up to at least 1–1.5 km above the ocean floor, making this part of the sea bed among the strongest source regions for the oceanic²²⁶Ra. Numerical solutions of a two-dimensional vertical advection-diffusion model applied to the deep (1.2–4 km)²²⁶Ra data give the following set of best fits: upwelling velocity(V_z) = 3.5m/yr, vertical eddy diffusivity(K_z) = 0.6cm²/s, horizontal (north-south) eddy diffusivity(K_y) = 1 × 10⁷cm²/s, and water-column regeneration flux of²²⁶Ra(J) = 3.3 × 10^?5dpmkg^?1yr^?1 as an upper limit. These parametric values are in general agreement with one-dimensional (vertical) model fits for the Ra-Ba system. However, consideration of²²⁶Ra balance leads us to suspect the appropriateness of describing the vertical exchange processes in the eastern Pacific with constantV_z and K_z. If future modeling is attempted, it may be preferable to treat the area as a diffusion-dominant mixing regime with depth-dependent diffusivities.     

Global electron content during solar cycle 23

Global electron content (GEC) as a new ionospheric parameter was first proposed by Afraimovich et al. [2006]. GEC is equal to the total number of electrons in the near-Earth space. GEC better than local parameters reflects the global response to a change in solar activity. It has been indicated that, during solar cycle 23, the GEC dynamics followed similar variations in the solar UV irradiance and F _10.7 index, including the 11-year cycle and 27-day variations. The dynamics of the regional electron content (REC) has been considered for three belts: the equatorial belt and two midlatitude belts in the Northern and Southern hemispheres (±30° and 30°–65° geomagnetic latitudes, respectively). In contrast to GEC, the annual REC component is clearly defined for the northern and southern midlatitude belts; the REC amplitude is comparable with the amplitude of the seasonal variations in the Northern Hemisphere and exceeds this amplitude in the Southern Hemisphere by a factor of ～1.7. The dayside to nightside REC ratio, R(t), at the equator is a factor of 1.5 as low as such a GEC ratio, which indicates that the degree of nighttime ionization is higher, especially during the solar activity maximum. The pronounced annual cycle with the maximal R(t) value near 8.0 for the winter Southern Hemisphere and summer Northern Hemisphere is typical of midlatitudes.     

Flow directions in ash-flow tuffs: a comparison of geological and magnetic susceptibility measurements,Tshirege member (upper Bandelier Tuff), Valles caldera,New Mexico,USA

This study concludes that the elongation axis (K ₁) of the ellipsoid of anisotropic magnetic susceptibility (AMS) is a suitable proxy for flow axis in ashflow tuffs. 153 oriented samples (176 specimens) were studied from 18 sites in the 1.1 Ma Tshirege member of the Bandelier Tuff. These sites are distributed around the Valles caldera at distances of 5–25 km outside of the rim.K ₁ axes correlate well with postulated radial flow axes at 13 sites.K ₁ also agrees with measured geological flow indicators, mainly imbricated larger clasts, at 7 sites. At 2 of the 5 sites where significant disagreement is seen between theoretical radial flow directions and measuredK ₁ axes, theK ₁ axes correspond well with geological flow indicators, indicating that the divergence of flow from the predicted radial flow pattern is real. Two major topographic buttresses are suggested as the cause of flow divergence for the Tshirege ash flows: the San Pedro buttress northwest of the caldera, and the San Miguel buttress in the southeast. In situK ₁ axes plunge about 7° toward the source at two-thirds of the sites; therefore the plunge ofK ₁ is a plausible in situ indicator for thedirection of flow. Multiple flow zones in sections of several meters thickness indicate changes of flow direction that are both rapid and large during ash-flow emplacement. These observations raisre the question of how best to represent ‘mean’ flow directions in ash-flow sheets: by eigenvector methods, by vector-sum methods, or by modes. A method for measuring imbrication of larger clasts using apparent dips in vertical joints is outlined. Imbrication, determined in this way at one-third of the sites, dips toward the source, i.e., up-flow. The minimum (K ₃) axis of the AMS ellipsoid correlates with the flow foliation rather than with the larger clast imbrication. The flow axes of ash flows correspond with theK ₁ axes, not with the declination ofK ₃ axes as suggested by some authors. Initial dip of the sampled ash flows is not large and does not affect the paleomagnetic remanence direction, which is reversed with a mean ofD=173.5°,I=-38.4°, α₉₅=3.4°N=18. This mean is not different at the 95% confidence level from that of earlier workers. The mean pole, at 098.0°E, 74.8°N,A ₉₅=3.3°,N=18, is about 15° far-sided relative to the expected time-averaged geomagnetic pole, suggesting a history of emplacement too short to adequately average secular variation.     

Variations in radon activity in the crustal fault zones: Spatial characteristics

The data of the profile gas emanation survey conducted on three spatial scales in separate regions of the Mongolia-Baikal seismic belt are generalized to establish the regularities of the spatially heterogeneous distribution of soil radon activity above the active faults in the Earth’s crust. It is shown that the shapes, sizes, and contrast of the near-fault radon anomalies are complicated by erosion and weathering; however, the critical role in their formation is played by the structural-geological controls, which determine the internal structure and recent activity of the fault zones. As a consequence, the cross-fault shape of the studied radon anomalies is vitally controlled by four structural situations, which correspond to the combinations of the structural type of the fault (localized/distributed) and the presence/absence of the fine filler material in the zone controlled by the fault. The cross-fault dimension of the emanation anomaly is commensurate or slightly larger than the width of the fault zone comprising all the fractures and joints associated with the formation of the main fault, which, due to the low permeability of the tectonites, is in most cases marked by the lowest concentration of soil radon. The contrast of the emanation anomalies, which we suggest to estimate in terms of a relative parameter K _Q , gravitates to certain levels of this parameter. This provides the basis for distinguishing five groups of the fault zones with low (K _Q ≤ 2), moderate (2 < K _Q ≤ 3), increased (3 < K _Q ≤ 5), high (5 < K _Q ≤ 10), and ultrahigh (K _Q > 10) radon activity. The previous studies show that for increasing the efficiency of the emanation survey in the fault zones, it is advisable to set up long profiles, reduce the measurement step in the vicinities of the main faults, specify the threshold of identifying the anomalies at the arithmetic mean level over the profile, and use the relative parameter K _Q for comparing and estimating the faults in terms of the intensity of their radon activity.     

Latitudinal dependence of the seasonal variation of particulate extinction in the UTLS over the Indian longitude sector during volcanically quiescent period based on lidar and SAGE-II observations

The altitude profiles of particulate extinction in the upper troposphere and lower stratosphere (UTLS) obtained from SAGE-II in the latitude region 0–30°N over the Indian longitude sector (70–90°E) are used to study the latitudinal variation of its annual pattern in this region during the volcanically quiescent period of 1998–2003. The SAGE-II data is compared with the lidar measurements from Gadanki (13.5°N, 79.2°E) when the satellite had an overhead occultation pass over a small geographical grid centered at this location. The particulate optical depth (τ_p) in the UT region shows a general decrease with increase in latitude and a pronounced summer–winter contrast with relatively low values during winter and high values during summer. In general, these variations are in accordance with the latitudinal variation of convective available potential energy (CAPE) and thunderstorm activity, which are good representative indices of tropospheric convection. While the particulate extinction (and τ_p) in the 18–21 km (LS₁) region is relatively low in the equatorial region up to 15°N, it shows an increase in the off-equatorial region, beyond 15°N. While the annual variation of τ_p in the LS₁ region is almost insignificant near the equator, it is rather well pronounced in latitude region between 10 and 15°N with relatively high values during winter and low values during summer. Beyond 20°N, this shows a prominent peak during summer. At a higher altitude, the 21–30 km (LS₂) region, the latitude variation of τ_p shows a different pattern with high values near the equator and low values in the off-equatorial region confirming the existence of a stratospheric aerosol reservoir. Low values of τ_p at lower regime (LS₁) near the equator could be due to rapid transport of particulates from the near equatorial region to higher latitudes, while the equatorial high at upper regime (LS₂) could be due to lofting and subsequent accumulation.