Aerosols over the Foothills of the Eastern Himalayan Region during Post-monsoon and Winter Seasons

In this study a Microtops II sunphotometer is used for the first time over Birtamode, a rapidly growing city located in the eastern Himalayan foothills of Nepal, to measure aerosol optical depth (AOD). The average AOD for the observation period (October 2018–February 2019) was 0.68 ± 0.39 with the post-monsoon season having a higher value (0.74 ± 0.43) compared to the winter season (0.60 ± 0.32). The Angstrom exponent ( α ) for post-monsoon and winter are found to be 1.08 ± 0.10 and 1.11 ± 0.16 respectively. During the monitoring period, the majority of AOD values (47%) are above 0.60, indicating moderately polluted conditions in the region. Anthropogenic, biomass burning, and mixed aerosols are identified as the prevalent aerosol types in the study region. The observed aerosol classification is also explained in terms of CAMS near-real-time model datasets. The AOD values retrieved by MODIS, VIIRS, Himawari-8, and CAMS show a good correlation with the observed Microtops AOD with R 2 values ranging from 0.60 to 0.94. Moreover, different MODIS aerosol products (DB, DT, and combined DB-DT) are evaluated based on a comparison of Collection 6.1 AOD with the ground truth obtained from Microtops. The spatial distribution of AOD as observed by various satellites are compared and the vertical distribution of aerosol is also explained with the extinction coefficient provided by the CALIOP lidar onboard CALIPSO and


INTRODUCTION
Atmospheric aerosols are the solid particles and liquid droplets suspended in the atmosphere varying between 0.001 µm to 100 µm in size. These originate from both natural and anthropogenic activities. They are known to have both warming and cooling effects on the atmosphere as they perturb incoming solar radiation and outgoing terrestrial long wave radiation by scattering or absorption. Aerosols perturb radiative forcing indirectly through cloud properties (Penner et al., 2003). Aerosols such as mineral dust and black carbon are identified to cause warming by absorbing solar radiation (Kaufman et al., 2002;Lohmann and Feichter, 2005). The microphysical and chemical properties of atmospheric aerosols strongly influence the direct and indirect effects 2 METHODS

Study Area
Birtamode is the socio-economic hub of one of the most populated districts of Nepal and is situated in the foothills of the eastern Himalaya (see Fig. 1). It is a highly urbanized region located in the Jhapa district, where a significant increase in urban area has occurred (from 27.12 km 2 in 2001 to 70.92 km 2 in 2016). It has an annual average rainfall of 270 cm and the temperature ranges between 32 to 35°C in summer and between 8 to 15°C in winter (Rimal et al., 2017). Though the air quality in Birtamode seems to be deteriorating with the increased urbanization, no targeted studies have been undertaken to investigate the issue. There is no air quality monitoring system in Birtamode with the nearest air quality observatory 85 km away in the neighboring district. To understand the aerosol properties in the eastern Himalaya, we set up a base station in the Birtamode region (26.65°N, 88°E, 137 m a.s.l).

Instrumentation
Ground-based spectral AOD was measured from 13 October 2018 to 28 February 2019 at an interval of 30 minutes/ 1 hour from 6:30:00 to 16:30:00 local time (LT) (which is 5 hours 45 minutes ahead of GMT). We used a hand-held Microtops II sunphotometer (manufactured by Solar Light Company Inc., USA) containing five channels at 340, 440, 500, 675, and 870 nm wavelengths that provide information about the columnar load of aerosols. The 340 nm channel has a bandwidth of 2 nm, and the other channels have a bandwidth of 10 nm (Sharma, 2018;Tiwari et al., 2018). The full field of view of the instrument is 2.5°. The details of the working principle of Microtops II instrument have been explained elsewhere (Morys et al., 2001).
Here, we would like to provide a brief overview of the working principle of the instrument. The solar radiation in an atmosphere is greatly influenced by Mie scattering, induced by aerosols, Rayleigh scattering produced by air molecules, and absorption by water vapor and other gaseous molecules (Ichoku, 2002a;Vijayakumar and Devara, 2012). The instantaneous solar radiation intensity, I, reaching the ground and that at the top of the atmosphere, I0 is related by Bouger-Lambert-Beer law as: (1) where τ * (λ) is the total aerosol optical depth, and m is the optical air mass, which is a relative path that the light travels in a slant path and is given by secant of the zenith angle (Pawar et al., 2012;Vijayakumar and Devara, 2012). AOD at each wavelength is estimated from τ * (λ) by subtracting the contributions due to the Rayleigh scattering and absorption at each wavelength. Spectral AOD and α are the primary parameters for the calculation of columnar aerosol properties; particularly, α contains information about the fraction of the fine versus coarse-mode particles (Schuster et al., 2006) and can be retrieved with Angstrom's power law (Ångström, 1964): where τ(λ) is the estimated AOD at the wavelength λ, α is the Angstrom exponent and β is the turbidity coefficient. The value of α and β obtained in this study are calculated in the wavelength interval of 340-870 nm by applying the least square method to Eq. (3).

MODIS
The MODIS instruments onboard the Terra and Aqua satellites have been observing the earth since 2000 and 2002 respectively . The equator crossing times for Terra and Aqua are 10:30 and 13:30 local solar time (LST) respectively. The instrument records spectral radiance in 36 channels with the resolution varying between 250 m to 1 km (Levy et al., 2007).
The MODIS aerosol products are based on three operational algorithms. These are Deep Blue (DB) over land, Dark Target (DT) over land, and DT over ocean. The enhanced DB algorithm is a modification of the DB algorithm that retrieves AOD values over both dark and bright surfaces. The algorithm uses the surface reflectance values at 0.412, 0.47, and 0.65 µm channels and estimates the AOD values based on the pre-calculated reflectance values in a lookup table (Hsu et al., 2013). The DT algorithm is prepared separately for ocean and land. Over the land surface, the DT algorithm utilizes surface reflectance at 0.47, 0.66, 0.86, 1.24, and 2.11 µm channels to retrieve AOD. The 1.38 µm channel is used for cloud masking. AOD is retrieved at 10 km and 3 km resolutions which aggregate 20 × 20 and 6 × 6 half km pixels respectively. Based on the reflectance value at 0.66 µm, the retrieval algorithm removes the darkest 20% and brightest 50% of pixels . The collection 6.1 (C6.1) AOD products have been used in this study as they are better compared to the collection 6 (C6) products (Wei et al., 2019) due to the application of modified surface parameterization, which incorporates the urban surface properties (Gupta et al., 2016). Similar to the C6 AOD products, the C6.1 AOD products also contain the merged DB-DT AOD values, which provide information for all transition regions with variable brightness. These AOD values are based on the MODIS NDVI value and include DT and DB AOD of superior quality obtained by applying the associated quality flags (Mhawish et al., 2017).

Himawari-8 AHI
The Advanced Himawari Imager (AHI) is installed on the geostationary earth-orbiting satellite-Himawari-8 which provides dynamic monitoring of the atmosphere for the Asia-Pacific region. The Japanese Aerospace Exploration Agency (JAXA) P-Tree data repository system provides access to Level 2 (L2) and Level 3 (L3) products that give the information of full disk at 500 nm with temporal resolution of 10 minutes and 30 minutes respectively. L2 AOD is based on five AHI bands (470, 510, 600, 800, and 1600 nm) (Fukuda et al., 2013). L3 AOD is a product obtained by processing L2 AOD and has two variables, AOD_Pure and AOD_Merged. AOD_Pure is obtained by removing cloudy pixels from L2 AOD, whereas, AOD_Merged is obtained by interpolating of AOD_Pure (Kikuchi et al., 2018). The data is accessible through ftp server after registering to JAXA Himawari Monitor P-Tree System (https://www.eorc.jaxa.jp/ptree/registration_top.html).

VIIRS
The Visible Infrared Imaging Radiometer Suite (VIIRS) is installed on the Suomi National Polar Orbiting Partnership (SNPP) satellite. SNPP is part of the environmental satellite system under operation by the Joint Polar Satellite System (JPSS) (Huang et al., 2016). The VIIRS aerosol products are derived at a spatial resolution of 0.75 km in the spectral region between 412 nm and 2250 nm (Jackson et al., 2013) aggregated to 6 km (i.e., 8 × 8 pixels). The VIIRS retrieval algorithm is based on the MODIS AOD retrieval algorithm architecture. The VIIRS aerosol product is available for download at the LAADS DAAC platform (https://ladsweb.modaps.eosdis.nasa.gov/).

CAMS
The Copernicus Atmosphere Monitoring Service (CAMS) is a component of the European Union's Earth observation program, Copernicus, which disseminates information about the Earth's environment based on modeling, satellite observations, and in-situ measurement data. CAMS uses the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). The IFS defines the sources of five tropospheric aerosols viz. dust, black carbon, organic, sea-salt, and sulfate aerosols, and incorporates five removal mechanisms viz. dry deposition, sedimentation, conversion of hydrophobic to hydrophilic BC and OM, transformation of SO2 to SO4, and wet deposition. A detailed description of the IFS mechanism is provided elsewhere (Morcrette et al., 2009). The near-real-time data available at 0 UTC and 12 UTC was used in this study. The data is available for download at the ECMWF data platform (https://apps.ecmwf.int/ datasets/data/cams-nrealtime/levtype=sfc/).

CALIPSO
The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) flies in a sunsynchronous orbit carrying the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, which crosses the equator at 2 pm and has a repeat cycle of 16 days. CALIOP receives backscatter signals at 532 nm and 1064 nm, and the profiles are available at a vertical resolution of 30 m and 60 m for altitudes below 8.2 km and between 8.2 km and 20.2 km respectively (Winker et al., 2013). The Level 2 aerosol vertical profile data (APro) was used to prepare the vertical profile of 532 nm Extinction Coefficient (bext). The data has a horizontal resolution of 5 km and was acquired from the data subsetter web app (https://subset.larc.nasa.gov/calipso/index.php).

Error Analysis
To avoid cloud contamination, Microtops measurements were undertaken on clear days or in instances when clouds were not occluding the solar disk. Once collated, any extremely high values suspected to have been affected by undetected cirrus clouds were removed. Also, AOD values beyond daily mean ±2 × standard deviation were discarded as outliers from further analysis. Inaccuracies in the fit of Angstrom's formula lead to curvature effect in the retrieved spectral AOD values. The comparison of α values obtained using different methods is used to minimize this effect. Hence, α values obtained from the linear fit method (LFM, use of five wavelengths) are correlated with those obtained from the Volz method (VM, use of two extreme wavelengths).
In this study, the CALIPSO data was averaged for the black box highlighted in Fig. 1. For the quality screening of the data, the cloud aerosol discrimination score (CAD_Score), extinction uncertainty (Extinction_Coefficient_Uncertainty_532), and extinction quality control flag of 532 nm (Extinction_QC_flag_532) were used. Features with a CAD_Score greater than -80 were removed and those layers which have a higher confidence that they are aerosols were retained. Only data (bext values) with Extinction_QC_Flag_532 values between 0 and 1 were used. Retrievals with an extinction coefficient of 99.9 km -1 were also removed following the method described by Winker et al. (2013).

Data and Collocation Strategy
This study being the first of its kind for the eastern Himalaya foothills region, MODIS aerosol products available with both the operational algorithms (DT and DB) were evaluated against ground-based observations. Moreover, the aerosol products provided by other polar orbiting and geostationary satellites and the CAMS model were also evaluated. The products selected and the quality filters used are tabulated in Table 1. As the timings of satellite overpass and the Microtops measurement don't match exactly, the collocation of these two types of data requires a spatiotemporal window (Jethva et al., 2019). The spatiotemporal collocation ensures that the air mass observed by the satellite sensors/model over the study region is also considered by the ground-based observations, which reduces the representation errors (More et al., 2013;Schutgens et al., 2017). Previous validation studies of MODIS aerosol products have used the spatiotemporal approach presented by Ichoku (2002b). In this study, we compare satellite AOD values averaged over a 0.3° × 0.3° grid box surrounding the ground station with the ground-based AOD averaged within ±30 min of MODIS and VIIRS overpass. Since Himawari-8 AHI has a higher temporal resolution of 10 min, the AHI AOD is averaged in the time window of ±30 min of ground-based observations. CAMS near-real-time analysis provides data at only two time-periods of 0 UTC and 12 UTC. So, daily averaged CAMS AOD is compared with the daily averaged ground observation. As the sunphotometer does not have the 550 nm channel, AOD values at 550 nm are calculated using Eq. (3).
To evaluate the uncertainty in the aerosol algorithms, the relative mean bias (RMB), the rootmean square error (RMSE), the mean absolute error (MAE), and the square of the Pearson correlation coefficient (R 2 ) are used. RMB values greater than 1 indicate overestimation whereas the RMB value less than 1 indicates the underestimation of AOD retrievals.
satellite model ground

RESULTS AND DISCUSSION
3.1 Synopsis of AOD, α, and β Fig. 2 shows the variation of Microtops AOD500 (unless otherwise stated, Microtops derived AOD hereafter refers to AOD500), α340-870 (unless otherwise stated, α here refers to α340-870), and β340-870 (unless otherwise stated, β here refers to β340-870) during the study period. For the study duration, AOD ranges between 0.13 and 1.93 with a mean value of 0.68 ± 0.39. Similarly, α values range between 0.34 and 1.49 with a mean value of 1.09 ± 0.13. β values range between 0.07 and 1.10 with a mean value of 0.30 ± 0.18. This variation in AOD, α, and β is due to differing atmospheric conditions driven by seasonal meteorology as well as changes in sources and types of aerosol emissions over the IGP (Tiwari et al., 2018(Tiwari et al., , 2016. Moreover, the seasonal average value of β shows that the aerosol loading during post-monsoon (0.33 ± 0.19) is higher compared to winter (0.27 ± 0.16), which is contrary to the results obtained from the central Himalayan foothills region, Pokhara (Regmi et al., 2020). The β values show significant correlation (R 2 = 0.964) with AOD and anticorrelation with α on majority of the days during study period.
The maximum recorded value of AOD during this study for October with the highest standard deviation is attributed to higher AOD values and variation during an observed haze episode. The percentage of total number of days with AOD and α values greater than the measured mean values are 42% and 52% respectively. Of the 1064 data points, 507 (47.65%) are above 0.6 indicating moderately turbid conditions around Birtamode compared to highly turbid areas like New Delhi (76%) (Tiwari et al., 2016) and Greater Noida (70%) . Nearly 18% of the days exhibited very high AOD values (i.e., AOD > mean +1σ) and on these days 88% of α values are greater than 1. This illustrates that the highly turbid conditions are due to the presence of finemode anthropogenic or biomass burning aerosols.
The AOD and α values from this work were compared with those over neighboring sites in the Indian subcontinent, and the results summarized in Table S1. The seasonal AOD values at Birtamode are comparable to those over Lumbini  and Kanpur (Kaskaoutis et al., 2012). The AOD values are lower compared to highly polluted urban areas like Delhi (Tiwari et al., 2016), Bhubaneswar (Mukherjee and Vinoj, 2020), and Varanasi (Tiwari et al., 2018). However, Birtamode was seen to have higher AOD values compared to other locations in the Indian sub-continent (Alam et al., 2012;Gogoi et al., 2009). Thus, Birtamode can be considered a moderately polluted semi-urban environment. Similarly, α values during the post-monsoon and winter periods are lower compared to several locations in the IGP and surrounding areas including Varanasi, Lumbini, Kanpur, Karachi, and Greater Noida. This suggests the presence of coarser aerosols over the Birtamode area. The monthly mean spectral variation of AOD (Fig. 3) shows significant differences between calendar months. Though the AODs at 870nm overall are quite similar (~0.4), October and February exhibit slightly higher AODs (~1) at the shorter wavelength, which are due to the presence of fine-mode aerosols (Kedia and Ramachandran, 2011;Tiwari et al., 2018). However, no different spectral variation of monthly mean AODs is found between December and January, and October and February. The wavelength dependence of seasonal AOD over Birtamode is comparable to other sites where data is available in the IGP Tiwari et al., 2018).
The frequency distributions of AOD, α, and β reveal a wide range for both seasons (Fig. S2) which indicates the variability in aerosol types and size distribution. The higher AODs recorded (> 0.6) with nearly 50% and 40% contributions in the post-monsoon and winter seasons, indicate moderately turbid conditions in Birtamode. Moreover, 67% and 31% of the observed β values greater than 0.2 and 0.4 respectively demonstrate that the post-monsoon period was more polluted compared to the winter period having 57% and 17% of β values greater than 0.2 and 0.4 respectively (Iqbal, 1983;Sapkota and Dhaubhadel, 2002). However, Varanasi, a nearby IGP location had 76% and 64% of higher AODs (> 0.7) in post-monsoon and winter respectively (Tiwari et al., 2018). The frequency distribution for AOD, α, and β values at Birtamode are unimodal during both seasons, with modal values close to 0.4, 1.2, and 0.2 respectively, which is attributed to the prevalence of certain aerosol types (Filonchyk et al., 2020).  Similarly, for Aqua MODIS (10 km product) the collocations were 62, 60, and 60 respectively for DB, DT, and combined DB-DT products. For 3 km products, the valid collocations were 75 and 73 respectively for Terra and Aqua. Moreover, the numbers of matchups were 90, 872, and 94 respectively for VIIRS, AHI, and CAMS. The relatively high number of collocations for AHI is due to its high temporal resolution of 10 minutes. MODIS C6.1 aerosol products (both 10 km and 3 km) have relative mean bias (RMB) values between 0.80 and 0.82 which suggests that they underestimate the AOD values by 18%-20%. Terra MODIS AOD retrieved from the DB algorithm (MOD04_L2_DB) has the highest R 2 of 0.945. Compared to Aqua MODIS, Terra MODIS shows better agreement with the ground observation which is indicated by the higher R 2 value. The evaluation of these products over other areas in the Asian region show results with lower R 2 values (Bilal et al., 2018;Huang et al., 2020;Sharma et al., 2021;Tian and Gao, 2019). This may be attributed to the lower number of data samples in the present study, and the difference in the spatial and temporal window size taken for collocation. The DT AOD and the merged DB-DT AOD show no difference for the Birtamode region. This is due to the merged product completely following the DT algorithm which indicates the region is moderately vegetated implied by the Normalized Difference Vegetation Index (NDVI) (Bilal et al., 2018). Similar observation of the DT AOD and DB-DT AOD have been reported with the C6 product for the Dhaka region (Mhawish et al., 2017).

Satellite and CAMS Retrieval Accuracy against Microtops Measurments
VIIRS AOD has an RMB value of 0.89 and an R 2 of 0.856 which is lower when compared to the MODIS aerosol products. A global evaluation study from 23 January 2013 to 31 December 2014 showed that the AOD retrieval generally underestimated AOD over India during agricultural burning seasons. This study had 24145 analysis samples with an accuracy of -0.0008, precision of 0.116, uncertainty of 0.116 and correlation coefficient (R) of 0.817 (Huang et al., 2016).
The correlation between AHI AOD and ground observations has the lowest R 2 value (0.6) and an RMSE value of 0.32. This overestimation may be attributed to the use of L2 products which is likely to be affected by cloud contamination. However, both L3 products, AODpure and AODmerged, don't show much improvement in the correlation against ground observations (Fig. S3). The global validation of L3 version 3 AOD shows that the R 2 values for different AERONET locations in nearby regions such as Bhola, Dhaka University, Pokhara, and Gandhi College are 0.66, 0.66, 0.77, and 0.38 respectively (Zhang et al., 2019). The underestimation of AOD by CAMS is comparable to that of MODIS. However, the R 2 value is lower (0.713) and the RMSE (0.21) and MAE (0.16) values are higher than MODIS products. A previous study shows that CAMS reanalysis performs better when the climatology is dominated by coarse mode particles relative to fine mode particles. The performance of CAMS near real-time data with the ground observation over the Birtamode region is better compared to other locations in the IGP such as Kanpur, Lahore, Gandhi College, and Jaipur where the reanalysis product was evaluated (Misra et al., 2020). Fig. 5 shows the spatial distribution of MODIS, VIIRS, AHI, and CAMS AOD over the eastern Himalaya foothills region for the post-monsoon and winter seasons. Although plume features are captured by various platforms, there are significant differences between them. The aerosol loads obtained from all platforms show higher AOD around the central IGP and the Bay of Bengal than the northeast areas including the Birtamode region. This is possibly due to anthropogenic activity including industrial activity in those regions. The evaluation results in the previous section show that the Terra MODIS (10 km product) has a higher accuracy than other products when compared to ground-based observation. So, the spatial distribution of VIIRS, AHI, and CAMS AOD were evaluated based on MODIS AOD retrievals. The percentage AOD difference in Fig. 6 is defined as "100 × (Other -MODIS)/MODIS", which describes the average difference between MODIS and other products. For consistency, the resolution of MODIS and the other products were kept constant by regriding to calculate the difference. During the post-monsoon, similar aerosol load distributions are found in MODIS and VIIRS plots over the IGP, although VIIRS tends to overestimate the value along the Brahmaputra and the Ganges rivers and underestimates the value along the eastern Himalaya. However, during winter, VIIRS overestimates the AOD value over eastern India Himalaya  and mostly underestimates values over the IGP. When AOD from AHI and MODIS are compared, the IGP region shows strong sensitivity to seasonal change; AHI AOD is significantly higher while the same region shows significantly lower MODIS AOD. AHI continues to overestimate along the foothills and the southern face of the eastern Himalayas and the Bay of Bengal. CAMS overestimates the AOD values over parts of the IGP and along the eastern Himalaya foothills while it underestimates the value over the Bay of Bengal during both seasons. The evaluation of AHI L3 data suggests that the AODpure and AODmerged products show improvement over both land and ocean which is demonstrated by lower overestimation as compared to AHI L2 data (Fig. S4).

Classification of Aerosols
Correlating AOD and α is one method for classifying aerosols (Kaskaoutis et al., 2009). To classify different aerosol types in Birtamode, scatter plots (Fig. 7) between AOD and α are used. This study used a similar threshold to those used in Varanasi (Tiwari et al., 2018) as AOD > 0.7 and α < 0.6 (Mostly dust; MD); 0.3 ≤ AOD ≤ 1 and α > 0.9 (Anthropogenic aerosol; AA); AOD > 1.0 and α > 1.0 (Biomass burning; BB) and those used in New Delhi  as AOD < 0.3 and α < 0.9 (CBC). The cases that don't belong to these categories are categorized as mixed type (MT) aerosols.
During the post-monsoon season, the contribution of AA type aerosol is the highest (~54%) which is followed by BB type aerosol (~22%). Similarly, AA and BB types during winter have contributions of 66% and 11% respectively. The high frequency of BB type aerosol during the post-monsoon may be attributed to extensive agricultural residue burning in the IGP . The CBC type has negligible contribution during both seasons. Similarly, a negligible contribution of MD (2.41%) is seen during winter which is attributed to the strong surface winds carrying dust from several construction sites in the vicinity. Fig. 8 shows the spatial distribution of CAMS aerosol types in the study domain. The organic matter AOD seems to have a major contribution (~50% each in post-monsoon and winter) of the total AOD, which is followed by sulfate AOD (38% and 40% in post-monsoon and winter respectively). The contribution of dust AOD and black carbon AOD is equal in post-monsoon (~5%) which is double i.e., ~3% by dust and ~6% by BC during winter. However, the contribution of sea-salt AOD is negligible during both post-monsoon and winter seasons. Though the temporal contribution of different aerosols varies, the hotspots of the black carbon, organic matter and sulfate remain the same as shown by the model output. A study over Kolkata, a neighboring megacity, with MERRA-2 reanalysis data showed that the seasonal average of sulfate AOD from 2001-2017 for both postmonsoon and winter is the highest followed by organic carbon, black carbon, and dust (Rawat et al., 2019). Given the higher emissions of SO2 from fossil fuel burning in the megacity versus rural  area, it is expected that sulfate would dominate the AOD over Kolkata versus our study domain that includes two megacities (Kolkata and Delhi) but is dominated by carbonaceous aerosols emitted from the surrounding rural environment.

Vertical Variation of Aerosol
The CALIPSO profiles for post-monsoon and winter seasons (Fig. 9) show the vertical distribution of aerosol within the domain where the plume from the IGP and the Bay of Bengal converge. Majority of aerosols are confined to altitude below 2 km during post-monsoon while there seems to be a slightly elevated layer with bext > 0.1 km -1 around 3 km. However, during winter, several elevated layers are present with bext > 0.2 km -1 above 2 km altitude. The CALIPSO profiles obtained for our study domain have been compared with other studies in the nearby locations. Unlike the present study, the vertical profile for the winter season in the central IGP doesn't show prominent features (Mahapatra et al., 2019). However, the study over Delhi NCR region shows that the daytime mean seasonal bext for the winter season has a secondary peak at an altitude of around 2.6 km which is comparable to this study with one of the secondary peaks lying at 3 km. During post-monsoon season, bext has a secondary peak at a relatively lower altitude (0.8 km) over Delhi NCR compared to the present study domain (Srivastava et al., 2014). Both Volume 22 | Issue 4 | 210152 studies do show a lowered elevated aerosol layer (implied via bext peak) during the post-monsoon season compared to the winter season. Climatological analysis of the planetary boundary layer (PBL) implies that PBL is lower over the IGP and the surroundings during the winter versus the post-monsoon season which should lead to higher elevated aerosol layers during the postmonsoon versus winter season (Putero et al., 2018;Raatikainen et al., 2014). The reason for the lowered elevated aerosol layer during the post-monsoon season versus winter needs further research. The Micropulse Lidar Network (MPLNET) -derived extinction coefficients over Kanpur (central IGP) during post-monsoon seasons show high values near the surface with no secondary peak value (Misra et al., 2012). For post-monsoon season, Dhaka, Agratala, Guwahati, Banmauk, Dhubri, Imphal, Aizawi, Shillong, Tawang, and Thimphu have high values of bext near the surface with secondary peaks at an altitude between 2.5-3 km similar to our findings. Whereas, Dibrugarh a city in the eastern Himalaya region has an elevated peak layer at 1 km altitude (Dahutia et al., 2019). The results indicate that even within the IGP there is a great seasonal variation, as well as the height of the elevated aerosol layers. The difference in the magnitude of bext values and the vertical profile for different locations may also be attributed to the different methodological approaches of using the CALIPSO dataset and the duration of studies. Fig. 10 shows the CAMS model derived vertical concentration profiles of organic matter, black carbon, and sulfate aerosols over the eastern Himalaya foothills region similar to that of CALIPSO analysis. Unlike CALIPSO observations, the model shows all these aerosols declining with altitude. Given the coarseness of the model vertical layers, it is expected that the model will not be able to resolve the layers obtained from the CALIPSO profiles. Model-based results do not show much difference in the vertical profiles of all three aerosol types between post-monsoon and winter season. Volume 22 | Issue 4 | 210152

CONCLUSIONS
The satellite-based observations of AOD show the hotspots over the eastern Himalayan foothills region. This study provides a comprehensive analysis of the ground-based AOD and various satellite-based and modeled aerosol products. A comparison of seasonal aerosol optical depth values over Birtamode, study site, with previous studies in the vicinity shows that the atmospheric condition of Birtamode is comparable to that of Lumbini and Kanpur. The AOD values are lower than the observed values at the highly polluted stations (Delhi and Varanasi) of the IGP region. The frequency distribution of AOD indicates that Birtamode is moderately polluted compared to other locations in IGP. However, the AOD and α values being higher compared with other cities in the Indian subcontinent, Birtamode may be considered as a polluted and aerosol-laden semi-urban area. The seasonal-average value of β signifies higher aerosol loading over the eastern Himalayan foothills region during the post-monsoon season, contrary to the observations recorded in the central Himalaya foothills region. Aerosols are classified based on the correlation of AOD and α, which shows that the contribution of anthropogenic aerosol type is the highest among all others. However, biomass burning contributions to the total aerosol loading are higher during the post-monsoon season. The evaluation of MODIS, VIIRS, AHI, and CAMS aerosol products was carried out for the first time in the eastern Himalaya region with Microtops ground observations. The MODIS AOD values from both Terra and Aqua satellites show a good correlation with groundbased observations compared to other platforms and provide the opportunity for further research based on MODIS values. Besides, evaluations of various MODIS AOD products such as DB, DT, and combined DB-DT have also been conducted. The DB product shows a better correlation for Birtamode despite the number of collocation points being slightly lower than DT and merged DB-DT products. The values obtained from the combined DB-DT algorithm didn't show differences from the values obtained from the DT algorithm, which is mainly due to the DT contribution in aerosol retrieval. AHI AOD (R 2 of 0.6) seems promising with a higher frequency of AOD points to monitor intra-day variation. The comparison of the spatial distribution of AOD shows that AHI and CAMS overestimate the value over the eastern Himalaya foothills, whereas VIIRS overestimates the value only during the winter season. The CAMS aerosol products were used to identify aerosol composition and vertical profiles over the eastern Himalaya foothills region. The model, as anticipated, was unable to capture the fine vertical aloft layers, unlike CALIPSO observations. CALIPSO vertical profiles show similar results as in north-eastern India, with the majority of the pollutants confined within 3 km of the boundary layer. However, the lowered elevated aerosol layer during the post-monsoon season versus the winter season requires further research.