Assessment of the Reverberations Caused by Predominant Air Pollutants on Urban Vegetation: A Multi-Site Study in Varanasi Located in Indo-Gangetic Plains

: Plant responses to air pollution have been extensively studied in urban environments. Nevertheless, detailed and holistic studies assessing their retaliation to air contaminants are still limited. The present study evaluates the effect of criteria pollutants (SO 2 , NO 2 , PM 10 and O 3 ) on the overall biochemistry and resource allocation strategy of plants in order to categorize the dominant roadside species ( Mangifera indica , Psidium guajava , Ficus religiosa , Azadirachta indica , Dalbergia sissoo , Cascabela thevetia and Bougainvillea spectabilis ) of the Indo-Gangetic Plains (IGP), with different morphologies and habits, into species that are tolerant and sensitive to the prevailing air pollutants. This study was performed at three different land-use sites (industrial, commercial and reference) in Varanasi for two seasons (summer and winter). It was inferred that NO 2 and PM 10 consistently violated the air quality standards at all the sites. The ﬁfteen assessed parameters reﬂected signiﬁcant variations depending upon the site, season and plant species whereupon the enzymatic antioxidants (superoxide dismutase and catalase) and resource utilization parameters (leaf area and leaf dry matter content) were remarkably affected. Based on the studied parameters, it was entrenched that deciduous tree species with compound leaves ( D. sissoo > A. indica ) were identiﬁed as the less sensitive, followed by a shrub ( C. thevetia > B. spectabilis ), while evergreen species with simple leaves were the most sensitive. It was also substantiated that the morphology of the foliage contributed more toward the differential response of the plants to air pollutants than its habit.


Introduction
Rapid industrialization and urbanization are the leading issues behind the alteration and deterioration of air quality and other environmental facets. Developing countries are quickly urbanizing, with a population increase from 751 million (30% of the world's population) in 1950 to 4.2 billion (55%) in 2018, an expansion of 4.6 times. The population is expected to reach 6.4 billion (68%) in 2050 [1]. Urbanization is essential for economic upturns, poverty eradication and overall human development. Furthermore, growing urbanization has posed a substantial obstacle to the United Nations Sustainable Development Goals [2]. One of the impediments is urban air pollution, which is a major public health concern worldwide. Every year, ambient air pollution leads to 4.2 million fatalities worldwide [3].
Attributed to their perennial nature and use in roadside plantations, trees and shrubs experience consistent exposure to the air pollutants [4]. Differential plant responses to air pollutants prove efficient proxies for monitoring air pollutants and providing abatement measures [5]. Tree plantation is a cost-effective measure to lower the pollutant burden, both gaseous and particulate, as trees provide a surface for absorption, deposition and Gases 2023, 3 58 dispersion, thereby reducing their human exposure [6]. The most widely assessed parameter for tracing the responses of perennial vegetation to air pollutants is the air pollution tolerance index (APTI), which compiles the changes occurring in the pH of the leaf extract, chlorophyll, ascorbic acid (AsA) and relative water content (RWC) of a leaf [7,8]. However, it is well-established that various abiotic adversities, including air pollution, generate oxidative stress in the exposed plants, leading to the increased generation of ROS, which ultimately harms plant functioning and causes cell death [9]. Plants involve different defense mechanisms to combat oxidative stress which are responsible for their tolerance/ susceptibility toward pollutants. The impacts of air pollution on plants manifest as a decrease in photosynthetic pigments, changes in metabolism, enzyme activities and impairments in physiological and metabolic processes [10]. As a consequence of air pollution, several biochemical anomalies occur in plants which encompass the generation of oxidative stress due to the increased production of ROS and a decline in photosynthetic pigments, which further adversely affects the photosynthetic processes [11]. Numerous enzymatic and non-enzymatic antioxidant defense mechanisms can lessen the negative effects of oxidative stress [10].
The present study was conducted with the aim of analyzing the impact of criteria pollutants, namely, sulfur dioxide (SO 2 ), nitrogen dioxide (NO 2 ), particulate matter (PM 10 ) and ozone (O 3 ), on seven perennial species (five trees and two shrubs), based on their dominance in roadside plantation at different geographical locations of Varanasi city (Uttar Pradesh) with variable land-use patterns.

Experimental Site Details
The experiment was conducted from December 2021 to May 2022, and the sites were selected on the basis of land-use patterns, namely, the Ramnagar Industrial area (R.I.A) (25 •  Varanasi experiences a prolonged hot and humid climate with a short low-temperature regime. There are three different seasons (summer, monsoon and winter). In this study, two prominent seasons were selected based on their temperature variations and pollutant loads, namely, summer (March-May) and winter (December-February). During the study period, the maximum temperature averaged between 24.24 • C and 40.70 • C, while the minimum temperature was between 9.39 • C and 26.52 • C. The total rainfall during the experimental period was recorded at 48.3 mm. The meteorological data for the experimental site were collected from the Indian Meteorology Division (IMD) observatory at Banaras Hindu University, Varanasi, India.
The monsoon season experiences a dip in pollutant load due to frequent rainfall (a total rainfall of 662.6 mm from June 2022 to August 2022); hence, it is not considered in the study.

Plant Details
Seven perennial species were selected from the experimental sites based on their dominance in roadside plantation. Among these species, five were trees and two were shrubs. The trees included the simple-leaved evergreen species Mangifera indica L. and Psidium guajava L., simple-leaved deciduous species Ficus religiosa L., compound-leaved semi-evergreen species Azadirachta indica A. Juss. and compound-leaved deciduous species Dalbergia sissoo Roxb. The shrub species included needle-leaved Cascabela thevetia L. and broad-leaved Bougainvillea spectabilis Willd. (Table 1). Prior to sample collection, each plant species was geotagged, and samples were collected from the same plant during each season at all the sites.

Plant Details
Seven perennial species were selected from the experimental sites based on their dominance in roadside plantation. Among these species, five were trees and two were shrubs. The trees included the simple-leaved evergreen species Mangifera indica L. and Psidium guajava L., simple-leaved deciduous species Ficus religiosa L., compound-leaved semi-evergreen species Azadirachta indica A. Juss. and compound-leaved deciduous species Dalbergia sissoo Roxb. The shrub species included needle-leaved Cascabela thevetia L. and broad-leaved Bougainvillea spectabilis Willd. (Table 1). Prior to sample collection, each plant species was geotagged, and samples were collected from the same plant during each season at all the sites.

Pollutant Monitoring
Throughout the entire experimental period, PM 10 monitoring was conducted for alternate weeks for an 8 h period at two sub-sites of each experimental location, using a high-volume sampler (Envirotech APM 460BL, Envirotech Instrument Pvt. Ltd., New Delhi, India). The sampler flow rate was fixed at 1.1 m 3 min −1 with a glass fiber filter paper (TISCH Scientific, USA) that had a pore size of 2.7 µm and dimensions of 20.3 × 25.4 cm, following the prescribed guidelines for National Ambient Air Quality Standards (NAAQS) by the Central Pollution Control Board (CPCB) [12]. The PM concentration was calculated using the methodology of Mukherjee and Agrawal [13].
The concentration of gaseous pollutants was estimated passively using wet chemical methods, following the methodology mentioned in the CPCB guidelines for the NAAQS [12]. SO 2 was estimated using the improved West and Gaeke method [14], NO 2 with the modified Jocob and Hochheiser method [15] and O 3 was quantified using the method of Hangartner et al. [16].

Plant Sampling
In each area, samples were collected randomly from three biological replicates of each species to minimize the local microclimatic effects. Fully expanded leaves with no sign of pest infection and damage (n = 10) were collected between 8.00 and 10.00 h every month during each season. The leaves were collected from the outer canopy and from all sides at a height of 1.5-2 m above the ground surface. The samples were carefully sealed in plastic bags and brought to the laboratory under cold conditions. The samples were washed properly and stored at −20 • C till further processing.

Photosynthetic Pigments
The protocols of Maclachlan and Zalik [17] and Duxbury and Yentsch [18] were used to estimate the photosynthetic pigments. Leaf extracts were extracted in 80% acetone, and their absorbance at 480, 510, 645 and 663 nm displayed carotenoids (carotene and xanthophyll), chlorophyll b (Chl b) and chlorophyll a (Chl a), respectively. The calculations used are as follows: Total chlorophyll (mg g −1 FW) = Chlorophyll a + Chlorophyll b where V is the volume of the extract (mL), W is the fresh weight of the sample (g) and d is the path length of light (cm).

Enzymatic and Non-Enzymatic Antioxidants
Enzyme extraction from the leaf samples was carried out as per the protocol of Takshak and Agrawal [19], with slight modifications. A leaf sample (0.1 g) was homogenized in 10 mL (0.1 M) potassium phosphate buffer (pH 7.8) consisting of 1% polyvinylpyrrolidone (PVP), 0.5% Triton-X 100 and 0.5 mM ethylenediaminetetraacetic acid (EDTA) and centrifuged twice at 13,000× g at 4 • C for 30 min. The supernatant was used further for the analysis of enzymatic antioxidants and proteins.
Catalase (CAT) activity was determined by following the method of Aebi [20], using potassium phosphate buffer (0.1 M) at pH 7.0 with 200 mM H 2 O 2 and 100 µL extracted enzyme. The decreasing absorbance was recorded at 240 nm. To assess the CAT activity, 0.036 mM −1 cm −1 was used as extinction coefficient in the relevant calculations. Superoxide dismutase (SOD) activity was estimated by a 50% reduction of nitro blue tetrazolium (NBT), following the methodology of Pandey et al. [21] with appropriate amendments. The reducing mixture contained 0.2 mM methionine, 2.25 mM NBT, 60 µM riboflavin, 3 mM EDTA, 0.1 M potassium phosphate buffer (pH 7.8), 1.5 M sodium carbonate (Na 2 CO 3 ) and enzyme extract. The absorbance was recorded at 560 nm.
To estimate the ascorbic acid (AsA) content in the leaf extract, a modified Keller and Schwager [22] method was utilized. A 10 mL solution, prepared by adding 0.75% oxalic acid and 0.05% EDTA, was used to extract AsA. The AsA concentration was estimated using its standard graph, following the spectrophotometric estimation of reduction of a 2,6-Dichlorophenolindophenol (DCPIP) reagent by the leaf extract at 520 nm. Phenol was estimated using the modified method of Bray and Thorpe (1954) [23], in which leaf samples were extracted in 80% acetone to estimate the concentration of total polyphenolics in the leaf extract. The content of total polyphenolics was determined using a Gallic acid standard graph, following the spectrophotometric determination of a blue-color complex generated by interaction of leaf extract with Folin-Ciocalteu reagent (1 N) in the presence of sodium carbonate (5%) at 650 nm. The leaf sample was homogenized in methanol (10 mL) with a drop of conc. HCl and 1% CaCO 3 to quantify anthocyanins [24]. The absorbance of the supernatant was measured at 535 and 650 nm after centrifugation. Flavonoid content was estimated following the protocol of Flint et al. [25].

Physio-Chemical Parameters
For the determination of the leaf pH, 1 g leaf sample was homogenized in 10 mL of double-distilled water and centrifuged. The supernatant was used to measure the pH using a pH meter (CyberScan-510, Thermo Scientific, Eutech Instruments, Waltham, MA, USA), following the protocol of Mukherjee and Agrawal [13]. The membrane stability index (MSI) was interpreted using the method of Gupta et al. [26], measuring the difference in the electrical conductivity (EC) of the leaf discs. The leaves were punched uniformly, avoiding the mid rib and in equal number. The leaf discs were incubated at room temperature for 2 h in 10 mL double-distilled water. The conductivity (EC1) of the incubated solution was recorded by a portable conductivity meter (Model 306, Systronics Limited, Ahmedabad, India). The solution was autoclaved at 120 • C and a pressure of 15 psi for 15 min, and the electrical conductivity (EC2) was then recorded again after cooling to room temperature. The following formula was used to calculate the MSI: Leaf discs were punched from fresh leaves and weighed to estimate fresh weight. After soaking the discs in distilled water for 16 h, the saturated weight of the discs was quantified. The discs were oven-dried and weighed again. The relative water content (RWC) was calculated using the formula of Smart and Bingham [27], as follows: The net water content (NWC) was calculated from the difference between the leaf fresh weight and the leaf dry weight.

Resource Utilization
The leaf dry matter content (LDMC) was estimated by dividing the dry weight of the leaves by their fresh weight. The leaf area of fresh leaf samples was calculated using Image J software (National Institutes of Health, Bethesda, MD, USA), and leaf mass per unit area (LMA) was calculated as the ratio of the dry weight of the leaves by the leaf area.

Statistical Analysis
The data from the leaf response parameters to air pollutants were analyzed using a oneway univariate analysis of variance (ANOVA) to investigate the effect of site on different parameters. Tukey's post hoc test was used after the one-way ANOVA for distinct measurements. A three-way ANOVA was used to assess the effect of the interaction between size and season on all of the evaluated plant attributes. To gain a better understanding of the relationship between air contaminants and plant performance, the Varimax rotation criterion with Kaiser normalization was used to perform a principal component analysis (PCA). In order to clusterize the plants based on the extent of the correlation and the similarity between them, a hierarchal correlation cluster analysis (HCA) was performed. A Pearson correlation measured the similarity between the plants. All the statistical analyses were performed using IBM SPSS Statistics21 and Origin 2023.

Pollutant Concentrations
The concentrations of the criteria pollutants for December 2021 to May 2022 are provided in Table 2

Photosynthetic Pigments
The total chlorophyll content is reflective of a plant's physiological health as it performs a key role in the process of photosynthesis [30]. Significant effects of site and season were observed on the chlorophyll a/b and carotenoid contents of the studied plants except A. indica and D. sissoo in the winter season (Table 3). The total chlorophyll content was found to be negatively correlated with an increasing pollutant concentration (R.I.A < B.R.S < B.H.U) in all the plants (Figure 2). The decline was most severe in P. guajava for R.I.A (68.59%) and B.R.S (53.79%), while it was the least severe in A. indica for B.R.S (1.87%) and R.I.A (6.93%) in the winter season. In the summer season, the reduction was pronounced in M. indica for R.I.A (65.94%), and it was at a minimum in B. spectabilis (19.27%). All the plants at R.I.A showed a higher variation for chlorophyll a/b, and the effect was more prominent during the winter season ( Figure 2). The order of reduction during the winter season was F. religiosa > P. guajava > C. thevetia > A. indica > D. sissoo > M. indica > B. spectabilis. For the summer season, P. guajava was found to be the most affected plant, followed by F. religiosa and C. thevetia. Conversely, the carotenoid content was found to be positively correlated with the pollutant load in all the species except A. indica and D. sissoo, which demonstrated reduced carotenoid contents (15.76-25.48%) for both seasons (Figure 2). In the summer season, B. spectabilis showed a maximum increase in the carotenoid content at B.R.S (59.05%) and R.I.A (44.29%), while in the winter season, M. indica showed a maximum increase at R.I.A (37%).

Enzymatic and Non-Enzymatic Antioxidants
To scavenge ROS, a plant's antioxidant defense system involves both enzymatic and non-enzymatic components. The enzymatic antioxidants revealed significant impacts of location and season. (Table 3) The SOD activity reflected a similar trend to catalase, i.e., it showed a significant increase with an increasing concentration of pollutants ( Figure 3).

Figure 2.
Seasonal and spatial variations in the photosynthetic pigment-total chlorophyll, chlorophyll a/b and carotenoid contents in different plant species during the study period (Mean ± SE). Bars with different letters show significant differences (p < 0.05) between the sites at specific season.

Enzymatic and Non-Enzymatic Antioxidants
To scavenge ROS, a plant's antioxidant defense system involves both enzymatic and non-enzymatic components. The enzymatic antioxidants revealed significant impacts of location and season. (Table 3)  The SOD activity reflected a similar trend to catalase, i.e., it showed a significant increase with an increasing concentration of pollutants ( Figure 3).

Physio-Chemical Parameters
The physio-chemical properties are often altered to a significant level upon exposure to air pollutants. All analyzed parameters established significant effects of site and season on the plants under consideration except for pH (Table 3). Here, the pH of the leaf extract showed a reduction during the winter season from R.I.A (11.70-24.

Physio-Chemical Parameters
The physio-chemical properties are often altered to a significant level upon exposure to air pollutants. All analyzed parameters established significant effects of site and season on the plants under consideration except for pH (Table 3). Here, the pH of the leaf extract showed a reduction during the winter season from R.I.A (11.70-24.80%), followed by The MSI was reduced during the winter season in all the plants ( Figure 5). The decline was the highest in the plants from R.I.A and maximum in the evergreen simple leaf species, followed by the semi-evergreen simple leaf, while it was minimum in the compound leaf and shrub species (P. guajava > M. indica > F. religiosa > A. indica > D. sissoo > C. thevetia > B. spectabilis). In B.R.S, both the shrubs showed an increase in MSI, while others showed decline, with highest decline found in D. sissoo (5.13%).
The RWC also declined in all the species at both the sites during winter season ( Figure 5). The reduction ranged between 1.66 and 6.54% except for C. thevetia and M. indica, which demonstrated an increase in RWC. In B.R.S, the decline ranged between 13.03 and 47.21% except in M. indica, P. guajava and F. religiosa. During the summer season, a uniform increment was observed at both the sites (R.I.A: 2.72-29.83%; B.R.S: 1.01-27.50%) and in all the plants, being maximum in M. indica and minimum in D. sissoo at both sites. The NWC showed low significance with respect to sites and seasons (Table 3). cline was the highest in the plants from R.I.A and maximum in the evergreen simp species, followed by the semi-evergreen simple leaf, while it was minimum in the pound leaf and shrub species (P. guajava > M. indica > F. religiosa > A. indica > D. sisso thevetia > B. spectabilis). In B.R.S, both the shrubs showed an increase in MSI, while o showed decline, with highest decline found in D. sissoo (5.13%).  (Table 3).

Resource Utilization Parameters
Plants' resource utilization strategies are reflected by their leaf functional among which LMA and LDMC are of prime importance [31]. In the present study, s icant effects of site and season on all the plants were observed for resource utiliz (Table 3).

Resource Utilization Parameters
Plants' resource utilization strategies are reflected by their leaf functional traits, among which LMA and LDMC are of prime importance [31]. In the present study, significant effects of site and season on all the plants were observed for resource utilization ( Table 3).
The leaf area was reduced at both the sites, the reduction being more at R.I.A  Figure 6). Conversely, during the summer season, LDMC decreas at both sites with an exception of P. guajava and A. indica. Figure 6. Seasonal and spatial variations in the resource utilization-leaf area, leaf mass per a and leaf dry matter content in different plant species during the study period (Mean ± SE). Bars w different letters show significant differences (p < 0.05) between the sites at specific season.

Discussion
Trees have been extensively used in urban air pollution abatement and are cons ently exposed to the increasing air pollution with an aid to provide a surface for the sorption and adsorption of the contaminants present in the surroundings. Henceforth is essential to choose suitable and efficient plant species for designing green belt or urb plantation based on their overall response to the prevailing pollutants. The present stu provides a detailed assessment of seven commonly used perennial plant species at th different land-use areas in two dominant seasons of the Indo-Gangetic plains and w conclude their susceptibility to air pollutants based on the HCA data (Figure 7).

Discussion
Trees have been extensively used in urban air pollution abatement and are consistently exposed to the increasing air pollution with an aid to provide a surface for the absorption and adsorption of the contaminants present in the surroundings. Henceforth, it is essential to choose suitable and efficient plant species for designing green belt or urban plantation based on their overall response to the prevailing pollutants. The present study provides a detailed assessment of seven commonly used perennial plant species at three different land-use areas in two dominant seasons of the Indo-Gangetic plains and will conclude their susceptibility to air pollutants based on the HCA data (Figure 7). Gases 2023, 3 69  This study concerns the relative performance of the dominant urban plant species in different land uses from the most polluted (R.I.A) location to a comparatively less polluted (B.H.U) location. Local weather conditions play a major role in determining the pollutant concentration over a region. The dust load was found to be at a maximum in the winter season at all the sites, attributing to the foggy conditions, wetness of the foliar surface and slow-moving wind. A consistently high concentration of PM in the commercial and industrial sites might be due to the heavy tailpipe emissions of the transportation sector and re-suspended soil dust in both seasons, with an additional increase in the burning of biomass during the winter season [32]. Conversely, the regional wind movement is high in summer season, which may be responsible for the re-suspension of dust from the foliar surface and hence less accumulation of the dust load [33]. The concentration of gaseous pollutants was also invariably higher during the winter season at all sites, which was found to be the highest at R.I.A and the lowest at B.H.U. This rise in the concentration of pollutants can be correlated with their restricted dispersion at the higher atmospheric levels which, in turn, enhances their ground level concentration [34]. Due to a sudden drop in temperature, a stronger and prolonged inversion layer develops, which worsens the air quality in the lower troposphere.

Photosynthetic Pigments
The total chlorophyll content was observed to decrease significantly with an increasing pollutant load, found to be minimum at R.I.A and maximum at B.H.U. This lends credence to the assertion by Wei et al. [35], which states that the chloroplast is the primary site of strike by air pollutants, mainly PM, SO 2 and NO 2 . Their entry into leaf tissues is regulated by stomatal apertures where they denature chlorophyll upon entering and reduce its concentration. It is also established that the phytotoxic behavior of SO 2 leads to a reduction in chlorophyll concentration via the interruption of the chloroplast layer, which will ultimately cause leaching of the pigment [36]. As they are acidic in nature, SO 2 and NO 2 will cause the acidification of the chlorophyll which, in turn, form pheophytin through the dispersion of magnesium by protons obtained from the acidic pollutants [37] and will further contribute to the reduction of the pigment concentration [38]. The increased PM 10 concentration at the industrial and commercial site might also have enhanced the degradation of chlorophyll content as a result of adhered heavy metals and the polycyclic aromatic hydrocarbons present on them, which have solubilized in the cell sap [39]. In the summer season, the concentration of chlorophyll was higher than in the winter due to the shorter photoperiod and higher PM load imparting a shading effect. It has been proposed that a plant's sensitivity is linked to the synthesis and breakdown of chlorophyll. Thus, plants with a high chlorophyll content are considered to be resistant to air pollution. In our study, the decline in the chlorophyll content was observed to be the most in P. guajava and M. indica, while it was the least in A. indica and B. spectabilis, thereby establishing that simple leaf and evergreen species are more susceptible to prevalent air pollutants in the area.
By enhancing the concentration of carotenoids, plants appear to be attempting to shield chlorophyll from damage. When cells are exposed to air pollution, the cell membrane is harmed, and the amount of thylakoid and chlorophyll in the chloroplast is reduced. In this study, carotenoid contents were significantly increased with the pollutant load in both the shrubs, while trees with compound leaves (A. indica and D. sissoo) showed a decline in carotenoids with an increasing pollutant load in both seasons. The decrease in carotenoids in the polluted areas might reflect the sensitivity of the carotenoids to SO 2 , as reported by Chauhan [4]. The trees with simple leaves (M. indica and F. religiosa) showed an insignificant increase in carotenoid content in summer, while the concentration increased with the pollutant load in the winter season. The increased concentration of carotenoids, the accessory photosynthetic pigments in both seasons at polluted sites, may be attributed to a reduction in the photooxidative stress on the primary photosynthetic pigments, which is consistent with the findings of Pellegrini et al. [40]. Further, the concentration of carotenoids can be correlated with the morphology of the leaves such that the concentration increased in the simple leaves and decreased in the compound leaves. The species with simple leaves showed more sensitivity to the pollutants due to their higher surface area for pollutant impingement compared to the compound leaves.

Enzymatic and Non-Enzymatic Antioxidants
To counteract oxidative damage, plants have devised enzymatic and non-enzymatic antioxidant defense machinery to scavenge hazardous ROS [41]. Antioxidative enzymes such as SOD and CAT may scavenge destructive ROS molecules or resist plants by triggering a non-enzymatic antioxidant mechanism [42]. SOD combats oxidative stress by converting superoxide radicals in mitochondria and chloroplasts into O 2 and H 2 O 2 [43]. Further, the CAT is involved in the protection of plant peroxidation by converting damaging intracellular H 2 O 2 to H 2 O and O 2 without using cell energy [44]. Plants increase their antioxidants in order to enhance their defensive property against the applied stress and also in response to the oxidative stress generated in them [45]. In this study, it was observed that the concentration of both the enzymes catalase and superoxide dismutase increased with an increase in the pollutant concentration. The enhancement was maximum for the simple leaf species (P. guajava and M. indica) during both seasons, reflecting a higher oxidative stress in them and a corresponding increase in the level of antioxidative enzymes, which is in line with the study of Tripathi and Gautam [46]. AsA acts as a protectant against the SO 2 -induced reactive oxygen species (H 2 O 2 , O 2− and OH − ), thereby safeguarding carbon fixation enzymes in addition to the inactivation of chlorophyll [47]. It in addition to leaf pH, is important in determining the sensitivity of a plant to SO 2 [48]. Its reducing power increases with an increasing pH and decreases with a reduction in the pH value. Thus, AsA may protect chloroplasts and chlorophyll activities from pollutants due to its pH-dependent reducing activity. In our study, the increase in AsA was maximum for M. indica and F. religiosa for both seasons, indicating higher oxidative stress in them.
Biological and non-biological stresses tend to accumulate the phenolic metabolites in plant tissues [49]. They are vital antioxidants that defend against the development and progression of the oxidation chain and reactive oxygen species by lowering or suppressing lipid auto-oxidation or degrading peroxides [50]. Despite the fact that phenols and flavonoids are present in plants as byproducts of related pathways, their concentrations may be elevated in response to oxidative stress. In our study, both the flavonoid and phenolic contents increased maximally in the shrub species (C. thevetia and B. spectabilis), thereby establishing a higher resistance in response to exposed stress. This is in accordance with the observations of Massad et al. [51] that the phenol production is boosted under stress at the expense of other primary metabolites, and this protects plants from stress caused by air pollution.

Physio-Chemical Parameters
Our study revealed a significant decline in the leaf pH with an increase in the pollutant concentration. This might be due to the impregnation of different air pollutants into the leaf cells, whereupon the plants with maximum decrease may be highly sensitive to air pollutants. The decline in pH can also be ascribed to the increased solubilization of SO 2 in the cell sap, which consequently lowers the pH [8]. Pollutant absorption results in acidification of intracellular pH, which is frequently seen in sensitive species. The reduction was minimum in both evergreen simple leaf species (M. indica and P. guajava) and maximum in shrubs, making them sensitive. Our findings are consistent with findings by Joshi et. al. [52], which revealed a direct correlation between a low leaf pH and the susceptibility of different plants to air pollution. Plants in a polluted environment (season and location) decrease their leaf pH in an attempt to convert carbohydrates (hexose) to ascorbic acid, hence inducing a mechanism to combat the generated oxidative stress. As a consequence of the pollutant load, membrane lipids frequently undergo phase transitions, affecting the permeability of the membrane and causing solute leakage, which is measured by MSI [26]. The MSI declined with an increasing pollutant load, and it was minimum in simple leaf evergreen species (P. guajava). A significant decline was observed in the RWC of the plants subjected to air pollution stress, reflecting their sensitivity to the subjected stress. This occurs in order to maintain the osmotic potential of the cells. The simple leaf species (M. indica, P. guajava and F. religiosa) maintained a relatively higher RWC among all species during the study period. The plants attempt to raise their RWC to maintain their physiological homeostasis in order to withstand abiotic stress, such as air pollution. The RWC was reported to analog with the protoplasmic permeability, which was frequently observed in the sensitive species [7]. The early senescence of leaves is a result of water and dissolved nutrient loss caused by pollutant-induced enhanced permeability in cells [53]. As a result, it is anticipated that plants with a higher RWC in contaminated sites will be resistant to pollutants [54].

Resource Utilization Strategy
The resource utilization strategy of the plants under study reflected that the tolerant/less sensitive species invest more resources in build up biomass at the expense of reducing their leaf area by increasing the LMA [13]. This is in line with our observation that the leaf area decreased with an increasing pollutant load while the LMA increased except in the simple leaf species (P. guajava, M. indica and B. spectabilis), making them sensitive to the exposed pollution concentration. LDMC is a significant determinant of a plant species' resource utilization strategy, i.e., its position in a fundamental trade-off between rapid assimilation and growth at one extreme and the efficient conservation of resources inside well-protected tissues at the other [55]. The present study showed a seasonal difference in the LDMC of the plants under study. However, in both seasons, the simple leaf species (M. indica and B. spectabilis) were noted to be more sensitive when compared to the compound leaf species (A. indica).
On the basis of the Pearson correlation coefficient, the hierarchical correlation cluster analysis (HCA) established the differential responses of plants under study (Figure 7). The species with simple leaves (M. indica, P. guajava and F. religiosa) were classified as more sensitive in both the polluted sites (R.I.A and B.R.S) irrespective of their habit (evergreen or semi-evergreen), followed by the compound leaf species (A. indica and D. sissoo). Among both the sites, the similarity was the highest between P. guajava and F. religiosa, while it was the lowest between M. indica and B. spectabilis at R.I.A. It was maximum between M. indica and P. guajava, whereas minimum for M. indica and B. spectabilis for B.R.S. For B.H.U, the similarity was the highest between C. thevetia and B. spectabilis and the lowest between F. religiosa and D. sissoo.

Principal Component Analysis
The principal component analysis shown in Figure 8 was performed individually for each of the three sites: the components were extracted with a 56.64%, 57.32% and 57.57% explained variance for R.I.A, B.R.S and B.H.U, respectively. The PCA analysis was performed for the winter season only as it was the period with the highest pollutant concentration. For R.I.A, the PC1 (eigenvalue 5.35) with a 33.43% variance was loaded with TC, Caro, Cat, SOD, RWC, LA and LDMC ( Figure 8A). The PC2 with an eigenvalue of 3.71 and a 23.21% variance showed higher loading for Chla/b, AsA, Flav, Antho, Phe, pH, MSI, NWC and LMA. Most of the antioxidants (enzymatic and non-enzymatic) and leaf area showed a higher correlation with the simple leaf species, i.e., M. indica, P. guajava and F. religiosa, while the compound leaf species (D. sissoo and A. indica) were more correlated with the photosynthetic pigments, LMA and LDMC. Both the shrub species were negatively correlated with all the parameters except RWC and AsA. For B.R.S, PC1 and PC2 showed eigenvalues of 4.78 and 4.39, respectively, with a 29.89% variance for PC1 and a 27.43% variance for PC2 ( Figure 8B). PC1 was loaded with Caro, Cat, SOD, AsA, Phe and MSI, whereas PC2 reflected loading of TC, Chla/b, Flav, Antho, pH, RWC, NWC, LA, LMA and LDMC. F. religiosa, P. guajava, and B. spectabilis showed a higher positive correlation with the resource utilization parameters and a negative correlation with the other analyzed parameters. In B.H.U, PC1, with an eigenvalue of 5.52 and a 34.50% variance, was loaded with Cat, AsA, Flav, Antho, Phe, MSI, NWC, LA and LDMC, while PC2, with an eigenvalue of 3.69 and a 23.07% variance, demonstrated a higher loading of TC, Caro, Chla/b, SOD, pH, RWC and LMA ( Figure 8C). D. sissoo, B. spectabilis and C. thevetia showed a higher correlation with the enzymatic antioxidants, physio-chemical and resource utilization parameters, while P. guajava and M. indica showed less correlation with all the studied parameters.  Figure 8C). D. sissoo, B. spectabilis and C. thevetia showed a higher correlation with the enzymatic antioxidants, physio-chemical and resource utilization parameters, while P. guajava and M. indica showed less correlation with all the studied parameters.  The differential response of plants appears to depend more on the leaf morphology than their habit. Leaf morphology is responsible for the per area load to exposed pollutants. The habit played a secondary role in response to the air pollutants in which the deciduous and semi-deciduous trees appeared to be less sensitive; this can be attributed to their leafless periods, during which their exposed surface to pollutants is drastically reduced and hence their adsorption and absorption are also reduced.

Conclusions
Seven dominant perennial species of Indo-Gangetic plains were assessed for their differential responses to air pollutants based on their morphological and biochemical parameters. The concentrations of NO 2 and PM 10 were consistently above the air quality standards in Varanasi. Pollution load and seasonal fluctuations had a major impact on differences in leaf functional characteristics. The analyzed photosynthetic pigments, enzymatic and non-enzymatic antioxidants, pH, MSI, RWC and resource utilization strategy of the plant revealed that the prevailing pollutant concentration had a more negative influence on evergreen trees with simple leaves than deciduous trees with compound leaves. The habit of the trees may have contributed to the observed responses indirectly, as the trees with simple leaves and evergreen habit will have a longer leaf lifespan which, in turn, will be expose them to the prevailing pollutant load for a longer period; therefore, more pollutant will be adsorbed/absorbed on their foliage. The trees with compound leaves and deciduous habit will have dual ease for the pollutant load, as they have less exposed surface and a shorter leaf lifespan, reducing the exposure time period and hence minimizing the impact of pollutants. Based on the assessed responses, it is inferred that in spite of different habits, plants with simple leaves tend to appear more sensitive to air pollutants than compound leaves. As the habit evolved from evergreen to deciduous with semi-evergreen in between, a resistance was induced against prevailing air pollution. Among both the perennial species, shrubs were more tolerant than trees when exposed to pollutant stress. Depending on the assessed parameters, C. thevetia among the shrubs and D. sissoo among the trees, followed by A. indica, were found to be less sensitive to the prevailing air pollutant concentration. It is highly recommended to promote mixed plantation with more compound leaf species in order to efficiently ameliorate the air pollutants from different land-use areas.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to privacy issues.