Geospatial Modelling for Estimation of PM 2.5 Concentrations in Two Megacities in Peninsular India

Airborne particles finer than 2.5 microns (PM 2.5 ) constitute a major public health risk in India. Therefore, extensive scientific studies must be conducted to assess the PM 2.5 exposures of Indians and determine the “exposure-response function” specific to India. While Peninsular India includes three megacities with populations exceeding 10 million each, there are very few studies on air quality modelling in this region compared to North India. In this paper, the authors describe a Linear Mixed Effects (LME) model to estimate monthly-average PM 2.5 concentrations at a spatial resolution of 1 km 2 between 2016 and 2019 in the megacities of Bengaluru and Hyderabad with a total population of 23 million. This model is based on covariates such as aerosol optical depth (AOD), meteorological parameters, and Land-use-Land-cover (LULC) variables and is validated with extensive datasets from continuous and manual air quality monitoring stations through a 10-fold cross-validation process. The final LME model can explain more than 60 percent of the variation in the PM 2.5 concentrations in Bengaluru and Hyderabad. This model is then used to predict the monthly-average grid-wise PM 2.5 concentrations in more than 800 grids in each of these two cities to study the spatial and temporal patterns in PM 2.5 concentrations between 2016 and 2019. These spatiotemporal maps of PM 2.5 concentration are critical to overcoming the misclassification of exposure and will form a crucial input to much-needed PM exposure-response studies in these two megacities. This paper can serve as a useful framework for similar studies by showing the way to bridge the gaps in the current air quality monitoring network in Peninsular India.


INTRODUCTION
Fine particulate matter (PM 2.5 ) pollution is a primary global public health concern.It is the fourth global leading risk factor for premature mortality and accounts for 4.14 million deaths globally (Health Effects Institute, 2020).Several air pollution exposure studies in the last three decades suggest that prolonged exposure to PM 2.5 pollution is associated inter alia with respiratory and cardiovascular mortality and morbidity (Dockery et al., 1993;Krewski et al., 2005;Jerrett et al., 2005Jerrett et al., , 2009;;Pope et al., 2002Pope et al., , 2009;;Miller et al., 2007;Samet et al., 2000).However, the PM 2.5 attributable mortality burden varies sharply across countries based on their population, per capita income and levels of PM 2.5 exposure (Prabhakaran et al., 2020).According to the Global Burden of Diseases (GBD) estimates, almost 58% of the global PM 2.5 attributable deaths occurred in Asia's two most populous countries, India and China which experience the highest PM 2.5 attributable age-standardized death rates of 89-98 per lakh population (Health Effects Institute, 2020).10-fold out-of-sample cross-validated (CV) R 2 of 0.81 and 0.724 in the Mid-Atlantic states and Mexico City, respectively.Maheshwarkar and Sunder Raman (2021) have shown that the spatial variability in the PM 2.5 concentration in Madhya Pradesh is reflected more accurately in the LME model than in the CTM based model of the same area.Several studies have used machine learning algorithms such as random forest, support vector machine, extreme gradient boosting, elastic net, and neural networks to improve prediction accuracy (Stafoggia et al., 2017;Di et al., 2016;Mandal et al., 2020).
In this paper, the authors have used the LME model based on satellite-derived AOD, geographical and meteorological covariates to predict the monthly PM2.5 concentration between January 2016 and December 2019 at a spatial resolution of 1 km 2 for the metropolitan regions of two cities in peninsular India-Bengaluru and Hyderabad.These cities have experienced rapid growth in land area as well as population during the last two decades due to the boom in the IT services sector in the last two decades.This boom has led to rapid increases in population, urban density, and number of vehicles.However, there are very few studies on the spatio-temporal changes in PM 2.5 concentrations in these two megacities of Peninsular India since most of the studies conducted in India are based on the National Capital Region and the Indo-Gangetic Plain.
The study period of this research was set as 2016 since this was the first full year during which ambient air PM 2.5 concentrations were first measured in these two megacities.The end of the study period was fixed as December 2019 to avoid the impacts of the unprecedented total lockdowns imposed in Bengaluru and Hyderabad for several weeks commencing 25 March 2020.Similarly, AQMS data recorded during 2021 are not considered since complete city-wide lockdowns to control the impacts of the 2 nd wave of the COVID-19 pandemic were imposed in Bengaluru and Hyderabad starting in April and May 2021, respectively.

Study Area
The geographical, climatological, demographic and topographical parameters related to the study areas in the megacities of Bengaluru and Hyderabad are shown in Table 1.The study areas include the respective municipal regions in Bengaluru (711 km 2 ) and Hyderabad (872 km 2 ).Based on the spatial orientation and shape of the land area within the respective municipal boundaries, the city of Bengaluru was divided into 801 grids (1 km × 1 km) matching with the spatial resolution of the Moderate Resolution Imaging Spectroradiometer (MODIS) -Multi-Angle Implementation of Atmospheric Correction (MAIAC) AOD, while Hyderabad was divided into 873 grids.The study areas in Bengaluru and Hyderabad are shown in Fig. 1.
Bengaluru and Hyderabad are located at different elevations (920 m and 545 m above Mean Sea Level, respectively).While Hyderabad remains warm throughout the year with maximum temperatures reaching 41 degrees centigrade, Bengaluru experiences moderate temperatures throughout the year with maximum temperatures limited to 34 degrees centigrade.In addition, Bengaluru experiences rain due to both southwest (SW) and northeast (NE) monsoons between June to September and October to November while Hyderabad receives most of the rainfall during the SW monsoon.
Bengaluru and Hyderabad have undergone rapid urbanization resulting in increased infrastructural activities and increased PM2.5 pollution in the last couple of decades.However, the AQMS in Bengaluru and Hyderabad are primarily clustered in specific parts of the city.Therefore, attributing the PM 2.5 exposure measured at these stations to people residing far from the AQMS would result  in exposure misclassification.This anomaly poses several difficulties in relating the personal exposure of the people to mortality and morbidity statistics related to chronic and acute illnesses.Therefore, the dose-response relationship between the health impact and personal exposure can be better established in developing countries like India with the help of a high-resolution spatiotemporal PM 2.5 exposure model.

PM 2.5 Data
The ambient PM 2.5 concentrations used in this study were collated from both manual and continuous AQMS established and maintained by the Government agencies such as the Central Pollution Control Board (CPCB), the Karnataka State Pollution Control Board (KSPCB) in Bengaluru and the Telangana State Pollution Control Board (TSPCB) in Hyderabad between 2016 and 2019.
In addition, we have collated daily PM 2.5 concentration measurements from continuous AQMS established by the Indian Institute of Tropical Meteorology (IITM) under the Modelling Air Pollution and Networking (MAPAN) project.The continuous and manual AQMS installed by the CPCB, KSPCB, and TSPCB measure PM 10 and PM 2.5 concentration based on beta ray attenuation and gravimetric method, respectively (CPCB, 2013).The MAPAN stations (one each in Bengaluru and Hyderabad) measure PM 2.5 using instruments calibrated based on U.S. EPA (Environmental Protection Agency) standards (Beig et al., 2021).The minimum detection limit of the continuous AQMS instrument used at these AQMS is 2 µg m -3 .To ensure the quality of the collated daily PM 2.5 data, suitable data filters are applied during this study.PM 2.5 measurements below 10 µg m -3 and not between µ * ± 3 σ * (*respective month-wise mean (µ) and standard deviation (σ)) were removed and considered missing (Mandal et al., 2020).In addition, PM 2.5 measurements greater than PM 10 measurements for co-located stations were removed and considered missing.
The monthly variability in the PM 2.5 concentrations recorded by continuous AQMS in Bengaluru and Hyderabad during the year 2019 is shown in Table 2.The PM 2.5 /PM 10 ratio provides information regarding the source of emission of particulates.While the annual PM 2.5 /PM 10 ratio in Bengaluru during 2019 varied between 0.36 and 0.52 in different AQM stations, this ratio varied in a narrow range (0.41-0.47) in Hyderabad.Therefore, the PM 2.5 and PM 10 measurements in the co-located AQMS were used to predict PM 2.5 from PM 10 measurements in Hyderabad, wherever PM 2.5  Hyderabad 49.5 40.2 41.3 33.3 21.1 15.5 10.0 10.7 15.4 57.2 76.6 73.5 36 measurements were not available.In this case, the PM 2.5 levels are derived from the measured PM 10 concentration using the LME model, with meteorological variables as covariates and month of the year as a random effect.The calibration model is cross-validated using the ten-fold crossvalidation method (Stone, 1974).In the ten-fold cross-validation method, the entire dataset is randomly divided into ten equal parts, where the nine parts are used to train the calibration model, and the left-out part is used to test the model.The ten-fold CV R 2 of the PM 2.5 calibration model is 85.1% and the root mean squared error (RMSE) is 9.8 µg m -3 (Fig. 2; Table 3).This PM 2.5 calibration model is then used to predict the 9480 daily PM 2.5 values from PM 10 measurements in AQMS across Hyderabad where PM 2.5 measurements are not available.In the next step, these

Aerosol Optical Depth
MODIS is an instrument placed in the Terra and Aqua satellites launched by NASA that performs measurements in the visible to thermal infrared wavelengths.The local equatorial passing times of the Terra and Aqua satellites are 10.30AM and 1.30 PM (Indian Standard Time), respectively.For this analysis, the daily MODIS AOD at 550 nm (AOD 550 ) data is derived using the MAIAC algorithm (Lyapustin et al., 2018).The MAIAC algorithm is chosen since it has a relatively finer spatial resolution of 1 km 2 and better agrees with AERONET stations than other retrieval algorithms.However, the AOD 550 observations over Bengaluru and Hyderabad are not continuous due to intermittent cloud cover, particularly during the monsoon season.In the case of Bengaluru, the AOD 550 observations between July and October are almost completely absent.To impute missing entries in the MAIAC-AOD database, it is calibrated against the global atmospheric reanalysis based AOD from MERRA-2 (Modern-Era Retrospective analysis for Research and Applications version 2) using the Goddard Earth Observing System Model (GEOS) (Gelaro et al., 2017).MERRA2 reanalysis data is available at a spatial resolution of 0.5° × 0.65° and 1-hour temporal frequency (Randles et al., 2017).The gaps in the MODIS AOD are computed using the LME model with MERRA 2 AOD and geographical coordinates (latitude and longitude) of the centroids of the grids as covariates.The day of the year (DOY) is used as a random effect.The model is applied for every year between 2016 and 2019.
The year-wise R 2 and RMSE for the MAIAC AOD calibration models for Bengaluru and Hyderabad are given in Table 4.The year-wise linear models between the predicted and observed MAIAC AOD values are shown in Figs.3(a In this manner, the best-fit LME calibration models (one for each of the four years between 2016 and 2019) between the MAIAC AOD550 observations and MERRA2 reanalysis AOD 550 data with the month of the year as random effect were used to fill the gaps in the AOD data.The daily predicted MAIAC AOD 550 are then averaged on a monthly basis for incorporation into the final LME model.

Urban built-up
The impact of increasing built-up area on urban air pollution is well documented in the literature (Gaigne et al., 2010).Therefore, the percentage of urban built-up area in each city is one of the critical predictor variables of the city's particulate pollution (Sarrat et al., 2006).To extract the percentage of urban built-up within each 1 km × 1 km grid in the Bengaluru and Hyderabad study areas, two Landsat-8 (Collection 1) images in each year (one each from pre-and post-monsoon seasons) between 2014 and 2019 were downloaded from the U.S. Geological Survey (USGS) portal (https://earthexplorer.usgs.gov/).The bands 2-7 ((Blue, Green, Red, NIR, SWIR I, II) were stacked using the Q-GIS platform (QGIS Documentation, 2021).The stacked raster image was exported to the Google Earth Engine.The K-means clustering method was used to classify the image into 25 classes (Lloyd et al., 1982).Further, the classified raster image was exported to the QGIS platform and compared with the original satellite image to get the different class numbers with similar spectral signatures.
In the case of Bengaluru, the classes with similar signatures were merged to form a combined class, leading to the classification of land use into four basic categories: urban built-up area, vegetated land, water bodies, and the barren land.However, in the case of Hyderabad, the merging of the similar signatures resulted in five final categories: urban built-up, vegetation, waterbody, fallow Land, and barren Land.It is noteworthy that the classified raster images of Hyderabad contained mixed pixels.Therefore, the reclassification of the raster image was done by overlaying the road network extracted from the Open Street maps portal (OpenStreetMap, 2021).While the satellite images of Hyderabad for 2015 and 2020 were classified using supervised classification in QGIS, the rest were classified using the K-means clustering method in the Google Earth Engine (Gorelick et al., 2017).Hyderabad.Due to the 17% increase in the population of Bengaluru and a 12% increase in that of Hyderabad during the study period 2016-2019, Bengaluru and Hyderabad experienced increases in urban built-up areas of 12% and 11%, respectively (Macrotrends, 2022).The monthly increase in the built-up area in each grid is obtained by performing a cubic spline interpolation, assuming a gradual increase in the built-up area during the study period between 2015 and 2019.

Road density
The land use map for 2015 was obtained from the Bengaluru Developmental Authority (BDA), Karnataka, and the open street maps were used to extract the grid-wise road density in Bengaluru and Hyderabad respectively (BDA, 2015;OpenStreetMap, 2021).The road density is calculated as the sum of all the road lengths (primary, intermediary, tertiary) in the 1 km 2 grid divided by the grid area.The grid-wise road density was calculated using the QGIS platform (QGIS Documentation, 2021).

Meteorological Covariates
While the emission of pollutants is one of the key factors in the ambient concentrations of any pollutant, meteorological variables also interact with the pollutants via convection, advection, deposition, dispersion, and dilution.In this study, daily meteorological variables such as Temperature, Relative Humidity, Planetary Boundary Layer height, Surface Pressure, Wind Speed, and Wind direction are obtained from the Indian Monsoon Data Assimilation and Analysis (IMDAA) regional reanalysis data.This single-level IMDAA regional reanalysis data is maintained by the National Center for Medium-Range Weather Forecasting (NCMRWF) under the Ministry of Earth Sciences (MOES), Government of India (Rani et al., 2021).This reanalysis data has a spatial resolution of 12 km and a temporal resolution of one hour (Rani et al., 2021).Further, bilinear spatial interpolation was used on the daily IMDAA meteorological data files over Bengaluru and Hyderabad to obtain daily meteorological data at a spatial resolution of 1 km.The daily meteorological variables thus obtained were arithmetically averaged for each month.In the case of the planetary boundary layer (PBL) height, the monthly average of daily planetary boundary layer height between 6:00 AM and 6:00 PM IST is calculated.While the monthly wind speed was calculated using the scalar averages, the monthly wind direction was computed using vector averaging.

Model Development
As shown in Fig. 1, most of the AQMS in Bengaluru and Hyderabad are located near point or line emission sources.Since the PM 2.5 measurements at these AQMS exhibit a right-skewed distribution, the natural log transformation is applied to the response variable 'PM 2.5 ' to ensure the homoscedasticity and normality of the residuals (Kloog et al., 2011).Both geographical and meteorological variables are used in the model as covariates.Since the covariates that can enter, the model is screened using correlation and stepwise regression analysis, all the covariates considered for the model development are not included in the LME model.The correlations between PM2.5 concentrations and AOD are statistically significant at the 95% confidence level in Bengaluru and Hyderabad with Pearson correlation coefficients of 0.315 and 0.349, respectively.The predictor variables such as temperature, relative humidity, wind speed, planetary boundary layer height, surface pressure, and road density correlate better with PM 2.5 measurements than any other predictor variables in Bengaluru and Hyderabad.Since several of the Spatio-temporal predictor variables used in the analysis have strong autocorrelation amongst themselves, including all the predictor variables in the model would result in a singular solution.Therefore, the predictor variables with the least autocorrelation are selected based on the stepwise regression method.The stepwise regression method is a statistical method that includes variables into the model until the R 2 value reaches saturation, after which there is no incremental change in the R 2 with the inclusion of any more predictor variables (Kutner et al., 1983).Stepwise regression was performed between the log-transformed monthly average grid-wise PM 2.5 concentration (response variable) and the corresponding Spatio-temporal predictor variables between 2016 and 2019 in both Bengaluru and Hyderabad.
Several researchers have used each day of the study period as a random effect (Kloog et al., 2011(Kloog et al., , 2014;;Just et al., 2015).However, a major part of the AQMS data in Bengaluru and are the common covariates selected for both the Bengaluru and Hyderabad LME models.While surface pressure and road density are also having a highly significant impact (p-values of 0.000) on the monthly PM 2.5 concentrations in Bengaluru, wind speed has a highly significant impact in Hyderabad (p-value of 0.000).All covariates (except AOD) of the grid-level average PM 2.5 concentrations listed in Tables 5(a) and 5(b) are highly significant at the 99% confidence level, while AOD is statistically significant at 95% confidence level.
The ten-fold CV final LME model diagnostics for Bengaluru and Hyderabad are shown in Tables 6(a) and 6 The cross-validated results of the model suggest that the final LME model can explain at least 60% of the variability in the monthly-average PM 2.5 concentrations in Bengaluru and Hyderabad.The grid-wise annual average PM 2.5 concentrations between 2016 and 2019 are predicted using the predictor variables belonging to the individual grid cells of the study area in Bengaluru (Fig. 6(a)) and Hyderabad (Fig. 6(b)).As shown in Fig. 6(a), hotspots of PM 2.5 concentration in Bengaluru are seen over the Peenya industrial area, the city railway station, the K.R. market, the Central Silk Board area, Whitefield, Hebbal, and Kalyan Nagar in all four years (2016-2019) irrespective of the seasons.Except for Peenya, the other hotspots are in areas witnessing dense traffic due to commercial activities.The south-eastern corners of Bengaluru have lesser PM 2.5 pollution for all the years between 2016 and 2019 due to the higher vegetation cover coupled with lack of commercial activities (Figs. 4(a) and 6(a)).
As shown in Table 2(a), the PM 2.5 pollution levels recorded in the continuous AQMS in Bengaluru peaked during the months between December and February.In the case of Hyderabad, the PM 2.5 pollution peaked between November and February (Table 2(b)).Though the PM 2.5 concentrations in most of the AQMS in Bengaluru comply with India's annual average National Ambient air Quality Standard (40 µg m -3 ) for PM 2.5 concentration, the monthly average PM 2.5 levels during the winter months are much higher than the annual NAAQ standard (Table 2(a)).As expected, the PM 2.5 concentrations were low during the southwest monsoon period (June-September) in both Bengaluru and Hyderabad, and during the northeast monsoon period (October-November) in Bengaluru.Unlike Bengaluru, which experiences low PM 2.5 levels for six months between June and November, Hyderabad experiences low PM 2.5 levels only for four months (June-September) since the NE monsoon does not touch Hyderabad.
As shown in Fig. 6(b), the city of Hyderabad has hotspots of PM 2.5 pollution over the central parts   2(a) and 2(b)).In Bengaluru, the PM 2.5 concentration declined between 2016 and 2019.This can be attributed to the major flyover and metro construction works that were taking place between 2012 and 2016 (Chaturvedi, 2012;Sastry, 2012;Mukherjee, 2012;The New Indian Express, 2017).The decline in the contribution of construction dust to the overall PM 2.5 concentrations in Bengaluru was also reported in two studies conducted by CSTEP firstly in 2015 and then in 2019 (CSTEP, 2022b;Guttikunda et al., 2019).
The annual average PM 2.5 concentration in Bengaluru declined from 40.8 µg m -3 in 2018 to 33.2 µg m -3 in 2019.This sharp fall of 18.6% is also due to the steep fall of 36% in the annual average PM 2.5 concentration recorded in the Information Technology Park Ltd. (ITPL) in Whitefield after the closure of Graphite India Ltd from February 2019 pursuant to the order of the National Green Tribunal (The Hindu, 2019).Four other AQMS in Bengaluru also recorded a decline in annual average PM 2.5 concentrations ranging between 19% and 26% while some stations recorded a marginal increase.
In the case of Hyderabad, the PM 2.5 concentration were higher in 2016 and 2019 compared to 2017 and 2018.While other factors may also be at play, the impact of the steep reduction (75%) in the total precipitation in Hyderabad during 2017 (246 mm) compared to that recorded in 2016 (990 mm) is one of reasons for the average PM 2.5 concentration recorded in Hyderabad in 2017 being 10% higher than that in 2016 (IMD, 2010(IMD, -2020)).The increase in total precipitation from 246 mm in 2017 to 607 mm in 2018 and 682 mm in 2019 has played a major role in reducing the average PM 2.5 concentration in Hyderabad from 51.3 µg m -3 in 2017 to 48.8 µg m -3 and 42 µg m -3 , respectively (IMD, 2010(IMD, -2020)).All existing air pollution models in India use air quality data recorded by continuous AQMS only.Due to the paucity of continuous AQMS in Peninsular India even in the megacities, the ground-truthing carried out in earlier studies is inadequate to study the PM 2.5 pollution in the megacities of this region (Gupta et al., 2020).As shown in Figs.1(a) and 1(b), while the number of continuous AQMS in Bengaluru is more than in Hyderabad, PM 2.5 levels were not recorded in most of these stations in 2016 and 2017.However, continuous AQMS are expensive, costing more than Rs.30 million to procure and install per station (MoEFCC, 2020b).Therefore, the final LME models developed for Bengaluru and Hyderabad in the present study have used the PM 2.5 concentrations recorded in the manual AQMS as well as the continuous AQMS.The extensive datasets on PM pollution in the megacities of Bengaluru and Hyderabad collated from a wide variety of sources (CPCB, KSPCB, TSPCB, and IITM) to develop LME models to estimate monthly average PM 2.5 levels at a 1 km × 1 km grid level with a much higher level of performance compared to earlier studies is a major contribution of this study.

CONCLUSIONS
Several PM pollution modeling studies have been conducted in India's National Capital Region (NCR) and the Indo-Gangetic plain (IGP).However, only a handful of air pollution modeling studies are published for the megacities of Bengaluru and Hyderabad with a combined population of 23 million.The population of these two cities have doubled in the last 20 years with consequent changes in land use-land cover, transportation infrastructure, etc.
In this study, the monthly-average PM2.5 concentrations between 2016 and 2019 for 801 (1 km × 1 km) grids in Bengaluru and 873 grids in Hyderabad are derived using the LME model.The LME model shows a satisfactory performance with a ten-fold CV R 2 of 65.5% and 61.6% for Bengaluru, and Hyderabad, respectively.This paper is the first to use the LME model for the megacities of Bengaluru and Hyderabad and successfully captures the spatial and temporal variability in the PM 2.5 concentrations in both study areas.Therefore it is a valuable addition to the literature on air pollution research.
Due to the competing priorities for public finance in a developing country like India, a robust geospatial model relating ambient air pollution with remote sensing and meteorological parameters and validated with both continuous and manual AQMS data will be useful to monitor PM 2.5 pollution levels by filling up the gaps in AQMS data.Therefore, such models are particularly useful for monitoring and control of ambient air pollution in the megacities of peninsular India.
The LME models developed during this study indicate the importance of LULC changes and meteorological parameters in determining the ambient air PM concentrations in both megacities.Therefore, urban planners must take care to provide adequate green areas and public transportation (to reduce road density) as cities grow in population and size.
The grid-level PM 2.5 concentrations estimated using this research methodology can help researchers overcome the misclassification of exposure to PM 2.5 pollution in the megacities of Bengaluru and Hyderabad.As mortality and morbidity data for Indian cities are widely available every month, the PM 2.5 estimates developed with this model can be used to correlate with the monthly mortality and morbidity statistics in Bengaluru and Hyderabad.Therefore, this method is particularly suitable for exposure-response studies that must be carried out in India to assess the efficacy of India's Air Quality Standards.Ultimately, proliferation of such studies in other megacities is crucial to reduce mortality and morbidity due to air pollution-related diseases and the attainment of Sustainable Development Goals (esp.SDG 3.9).
) and 3(b) for Bengaluru and Hyderabad, respectively.As shown in Figs.3(a) and 3(b), the MAIAC AOD calibration LME models performed well for all the four years between 2016 and 2019, for both Bengaluru and Hyderabad.

Fig. 4 .
Fig. 4. Land Use Land Cover classified images of (a) Bengaluru and (b) Hyderabad in March 2019.
(b), respectively.The RMSE, Prediction Error (MPE) and Relative Prediction Error (RPE) of the final LME model for Bengaluru were 8.39 µg m -3 , 6.4 µg m -3 , and 21.4%, respectively.The corresponding values in the case of Hyderabad are 11.3 µg m -3 , 8.2 µg m -3 , and 25%, respectively.The ten-fold CV R 2 values for Bengaluru and Hyderabad are 65.5% and 61.6%, respectively.The model validation parameters between the iterations of the 10-fold cross-validation procedure are consistent.The predicted versus measured values of PM 2.5 concentrations derived from the ten-fold CV LME models and the corresponding 95% confidence intervals are shown in Figs.5(a) and 5(b) for Bengaluru and Hyderabad, respectively.

Fig. 6
Fig. 6(b).Spatiotemporal distribution of model-derived annual average PM 2.5 concentrations in Hyderabad between 2016 and 2019.

Table 1 .
Key physical, geographical, and demographic information of Bengaluru and Hyderabad metropolitan areas.

Table 2 .
Monthly average PM 2.5 concentration for all the real-time AQMS in Bengaluru and Hyderabad.

Table 3 .
PM 2.5 calibration model validation parameters in Hyderabad..5 values are added to the PM 2.5 measurements database that is used to build the final LME model.However, this method is not followed in Bengaluru due to the wide range of PM 2.5 /PM 10 ratios between different stations.This variation in the PM 2.5 /PM 10 ratio in Bengaluru is mainly due to the variability in the sources of emission around air pollution monitoring stations located in different areas that has also been documented in a study conducted by CSTEP (2022a) in this city.Consequently, the datasets of monthly PM 2.5 data used to develop the LME model had only 420 values in the case of Bengaluru and 1145 values in the case of Hyderabad.