Air mass physiochemical characteristics over New Delhi: impacts on aerosol hygroscopicity and cloud condensation nuclei (CCN) formation

Delhi is a megacity subject to high local anthropogenic emissions and long-range transport of pollutants. This work presents for the first time time-resolved estimates of hygroscopicity parameter (κ) and cloud condensation nuclei (CCN), spanning for more than a year, derived from chemical composition and size distribution data. As a part of the Delhi Aerosol Supersite (DAS) campaign, the characterization of aerosol composition and size distribution was conducted from January 2017 to March 2018. Air masses originating from the Arabian Sea (AS), Bay of Bengal (BB), and southern Asia (SA) exhibited distinct characteristics of timeresolved sub-micron non-refractory PM1 (NRPM1) species, size distributions, and CCN number concentrations. The SA air mass had the highest NRPM1 loading with high chloride and organics, followed by the BB air mass, which was more contaminated than AS, with a higher organic fraction and nitrate. The primary sources were identified as biomassburning, thermal power plant emissions, industrial emissions, and vehicular emissions. The average hygroscopicity parameter (κ), calculated by the mixing rule, was approximately 0.3 (varying between 0.13 and 0.77) for all the air masses (0.32± 0.06 for AS, 0.31± 0.06 for BB, and 0.32± 0.10 for SA). The diurnal variations in κ were impacted by the chemical properties and thus source activities. The total, Aitken, and accumulation mode number concentrations were higher for SA, followed by BB and AS. The mean values of estimated CCN number concentration (NCCN; 3669–28926 cm−3) and the activated fraction (af; 0.19–0.87), for supersaturations varying from 0.1 % to 0.8 %, also showed the same trend, implying that these were highest in SA, followed by those in BB and then those in AS. The size turned out to be more important than chemical composition directly, and the NCCN was governed by either the Aitken or accumulation modes, depending upon the supersaturation (SS) and critical diameter (Dc). af was governed mainly by the geometric mean diameter (GMD), and such a high af (0.71±0.14 for the most dominant sub-branch of the SA air mass – R1 – at 0.4 % SS) has not been seen anywhere in the world for a continental site. The high af was a consequence of very low Dc (25–130 nm, for SS ranging from 0.1 % to 0.8 %) observed for Delhi. Indirectly, the chemical properties also impacted CCN and af by impacting the diurnal patterns of Aitken and accumulation modes, κ and Dc. The high-hygroscopic nature of aerosols, high NCCN, and high af can severely impact the precipitation patterns of the Indian monsoon in Delhi, impact the radiation budget, and have indirect effects and need to be investigated to quantify this impact.

polluted clouds. Complete suppression of warm rain might also occur and get aggravated due to additional CCN activation above the cloud base (Braga et al., 2017). While rain suppression was observed in the case of polluted urban and industrial plumes (Rosenfeld, 2000), and smoke arising from forest fires (Rosenfeld, 1999), the precipitation tendency increases due to influx of giant CCN consisting of sea salt  and salt playas (Rudich et al., 2002) due to acceleration of the auto-conversion rate . To understand the impact of pollution on indirect radiative forcing and precipitation in highly polluted regions, the information on CCN number concentration is essential in global climate models (GCMs) and regional climate models (RCMs).
As per the fifth IPCC report (Boucher et al., 2013), the two most important factors governing CCN activation and number concentration are size and composition. The aerosol chemical composition impacts the aerosol hygroscopicity, which impacts the critical diameter and hence CCN activation. The hygroscopicity parameter (κ) is defined as the total water uptake ability of aerosols (Petters and Kredenweis, 2007). Further, the increase in relative humidity (RH) due to water uptake by aerosol can impact visibility (Lee et al., 2016;Liu et al., 2012); secondary particle formation (Ervens et al., 2011); and measurements of remote sensing (Wang and Martin, 2007;Brock et al., 2016), aerosol loading, and its chemical composition (Chen et al., 2018). Hence, it is essential to determine the hygroscopicity of aerosols, especially in the polluted regions of the world, where these impacts are expected to be highly significant.
Although the recent precipitation data during 1950-2011, averaged over July and August for Delhi, reveal a significant decreasing trend, there has been an increasing trend in the frequency of heavy rainfall events and a decrease in the frequency of wet and rainy days when it rains for a shorter period (Guhathakurta et al., 2015). These occurrences are most likely signatures of aerosols impacting the cloud-nucleating properties, which calls for detailed CCN data examination. The high uncertainties associated with radiative forcings, both direct and indirect, especially at the regional level, are a result of poor representation of the aerosol distributions in GCMs. This is critical for the Indian sub-continent, where the variability in aerosol microphysical properties is very high, at various spatial and temporal scales. These necessitate the measurement of long-term aerosol physiochemical properties, the hygroscopicity parameter, and CCN estimates. Detailed CCN and κ measurements have been carried out in different parts of the world (Rissler et al., 2004;Bougiatioti et al., 2011;Engelhart et al., 2012) and in India, in places like Kanpur Tripathi, 2014, 2015;Ram et al., 2014), Mahabaleshwar (Leena et al., 2016), and the eastern Himalayas (Roy et al., 2017). However, no CCN measurements or estimates have been developed so far for Delhi. There is only one study that has estimated aerosol hygroscop-icity based on PM 2.5 mass, RH, and visibility data (Wang and Chen, 2019).
In this work, for the first time for Delhi, time-resolved size distribution and chemical speciation measurements were carried out from 15 January 2017 to 31 March 2018 as a part of the Delhi Aerosol Supersite (DAS) campaign (Gani et al., 2019). The time-resolved hygroscopicity parameter and CCN estimates were derived using chemical speciation data from an aerosol chemical speciation monitor (ACSM), and number concentration data from a scanning mobility particle sizer (SMPS), measured from January 2017 to March 2018. Data were analysed to investigate the following hypotheses: (a) the precursors to secondary organic aerosol (SOA) formation critically impact the chemical composition over Delhi, (b) the emission sources significantly impact CCN formation by governing the size distributions and chemical composition, and thus hygroscopicity, and (c) physical properties impact CCN more compared to chemical properties directlyhowever, the physical properties are, in turn, shaped by the chemical properties.

Instrumentation
For a detailed assessment of aerosol physiochemical properties, an SMPS (TSI, Shoreview, MN), ACSM (Aerodyne Research, Billerica MA), and aethalometer (Magee Scientific Model AE33, Berkeley, CA) were operated at the Indian Institute of Technology (IIT) Delhi in Block 5, at a height of nearly 15 m, as a part of the DAS campaign. This sampling site in New Delhi was free from any source activity, except for a road, located 150 m away. The IIT campus is cleaner than the rest of the city. However, it lies in the heart of the city, and the outskirts of the campus experience fresh traffic influx. The IIT campus allows only limited access to vehicles and therefore has less traffic compared to the city in general. A temperature-controlled room was used to carry out the measurements. Two separate and thermally insulated sampling lines (3/4 in. outer-diameter stainless-steel tubes) with flows of 3 and 2 L min −1 equipped with PM 1 cyclone, in line with a water trap and a Nafion membrane diffusion dryer (Magee Scientific Sample Stream Dryer, Berkeley, CA), were used for (1) an SMPS and ACSM in conjunction with a flow controller and (2) an aethalometer, respectively. A brief description of the instruments is given below. A detailed description of the instruments is given in Gani et al. (2019) and Bhandari et al. (2020).
The SMPS was comprised of a differential mobility analyser (DMA; TSI 3081), an electrostatic classifier (TSI 3080), an X-ray aerosol neutralizer (TSI 3088), and a water-based condensation particle counter (CPC; TSI 3785). The ambient air was sampled in the size range 12-560 nm, with a time difference of 135 s between two scans. The sheath-to-aerosol flow ratio was 4 : 1, and the total flow drawn by the CPC was 1 L min −1 . The two dominant modes (Aitken and accumulation) are well captured within this size range. The ACSM sampled the inlet air at a flow rate of 0.1 L min −1 with a time resolution of ∼ 1 min. The calibration procedures and data processing are discussed in Gani et al. (2019). The ACSM collected time-resolved NRPM 1 (non-refractory PM 1 ) based on species that volatilize by 600 • C and included NO − 3 , Cl − , SO 2− 4 , NH + 4 , and organics.
2.2 Qualitative separation of organic aerosols as BBOA, HOA, and OOA The composition data presented in this work were collected in the DAS study. PMF (positive matrix factorization) analysis was conducted on the 15 months in the dataset. As a result, biomass-burning organic aerosol (BBOA) could be resolved as a separate factor only in spring 2018. This inability to resolve primary organic aerosol (POA) to separate factors, namely hydrocarbon-like organic aerosol (HOA) and BBOA, was attributed to the unit mass resolution of the instrument (Bhandari et al., 2020, and references therein). Owing to the lack of explicit BBOA and HOA separation in all seasons, the Ng et al. (2010) compilation of profiles was analysed, combined with the profiles identified in spring in Delhi. It was observed that spring 2018 profiles fell within the bounds of the uncertainty of the Ng et al. (2010) compilation. Thus, Ng et al. (2010) reference profiles were utilized for source attribution of each cluster. While factor profiles can differ across the world, taking regionally relevant profiles together with those usually employed as reference profiles for PMF analysis likely accounts for this variability. As a part of the analysis conducted here, the mean strength at the relevant m/z (s) (m/z 57 and 60) and the standard deviation (SD) of the profiles at these m/z (s) were utilized in the analysis. Organic aerosols were qualitatively segregated by comparing the m/z ratios of f 57 and f 60 with the reference profiles of BBOA (f 57 : 0.0337±0.00884, f 60 : 0.025±0.00521), HOA (f 57 : 0.0838 ± 0.00378, f 60 : 0.00227 ± 0.00214), and oxygenated organic aerosol (OOA) (f 57 : 0.00997 ± 0.00786, f 60 : 0.00571±0.00349), as reported by Ng et al. (2010). This was done by first calculating the cluster means of f 57 and f 60 for each cluster (Table 2). This was followed by an evaluation of residuals. Residuals represent the deviation of the cluster means from the reference profiles. The HOA, BBOA, and OOA residuals (R HOA , R BBOA , and R OOA ) were then calculated based on cluster means of f 57 and f 60 (CM f 57 and CM f 60 ), with respect to the corresponding means of reference profiles (RM f 57 and RM f 60 ), as given below for HOA in Eq. (1): The reference residuals for HOA, BBOA, and OOA (R Ref_HOA , R Ref_BBOA , R Ref_OOA ) were then calculated using standard deviations of reference profiles (SD f 57 _Ref_HOA , SD f 57 _Ref_BBOA , SD f 57 _Ref_OOA ), as given below for HOA in Eq. (2): The residuals of the cluster means were then compared with the reference residuals as per the six conditions described in detail in Sect. S1 in the Supplement and classified as HOA, BBOA, OOA, or mixed.

Estimation of κ and CCN
The ACSM data were used to calculate κ as per the following mixing rule in Eq.
(3) (Petters and Kredenweis, 2007): where i and κ i represent the volume fractions and individual hygroscopicity parameters of the various components. The inorganics were represented by (NH 4 ) 2 SO 4 , NH 4 Cl, and NH 4 NO 3 . The organic κ i was taken as 0.1 (Gunthe et al., 2009Dusek et al., 2010;Rose et al., 2011). κ i values were taken as 0.61 for (NH 4 ) 2 SO 4 , 1.02 for NH 4 Cl, and 0.67 for NH 4 NO 3 (Sullivan et al., 2009;Petters and Kredenweis, 2007). The density values to estimate the volume fraction of the inorganic constituents were taken as 1770 kg m −3 for (NH 4 ) 2 SO 4 , 1519 kg m −3 for NH 4 Cl, and 1720 kg m −3 for NH 4 NO 3 (Haynes et al., 2014). The density of organics was taken as 1500 kg m −3 (Bougiatioti et al., 2009). κ for BC was taken as zero, as reported in several studies (Hong et al., 2014;Leng et al., 2014;Wu et al., 2013). It should be noted that we assumed that the κ calculated from NRPM 1 data of ACSM represents the bulk hygroscopicity parameter in the absence of size-resolved measurements and is a limitation of this work. The difference due to the assumption cannot be accounted for and should be investigated in the future. However, it is reported that for κ > 0.1, CCN closures within 20 % can be achieved assuming bulk composition and internal mixing (Wang et al., 2010). Temperature, relative humidity (RH), and the calculated κ were then used to calculate the critical diameter (D c ) from the multicomponent κ-Köhler theory . The temperature and RH data are available from the RK Puram site (∼ 3-4 km aerial distance from the measurement site) maintained by the Central Pollution Control Board (CPCB), India. N CCN was then estimated by integrating the size distribution obtained from SMPS above D c . The CCN estimates were obtained for the supersaturation (SS) of 0.1 %, 0.15 %, 0.2 %, 0.35 %, 0.4 %, 0.5 %, 0.6 %, 0.7 %, 0.75 %, 0.8 %, 0.85 %, and 1 %. However, for the sake of detailed analysis, 0.1 %, 0.4 %, and 0.8 % were chosen: 0.1 % represents the condition when the effect of chemical composi-tion is expected to be the highest, 0.4 % represents the condition for convective clouds, and 0.8 % represents a highsupersaturation state, when almost all aerosols tend to get activated as CCN.

Air mass characterization
To characterize the air masses, the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used (Draxler and Rolph, 2003) to determine the major pathways of aerosols reaching Delhi. The 5 d back-trajectory analysis was done at the receptor site at a height of 500 m. The cluster analysis was then performed seasonally to identify the cluster mean trajectories per season. These mean trajectories were then again re-clustered to identify three main clusters based on the directions of the mean cluster trajectories: the Arabian Sea (AS) branch (16.5 % of total trajectories), the Bay of Bengal (BB) branch (13 % of total trajectories), and the southern Asian (SA) branch (70.5 % of total trajectories). The BB branch was further classified as B (54 %) and B.reg (45 %), where B represents the air masses reaching the sea, while B.reg represents the air masses that aligned towards reaching the Bay of Bengal but did not hit the sea. The SA branch was partitioned into L (17.5 %), R1 (54 %), R2 (18 %), and R3 (11 %). L represents the local trajectories originating within India, mainly from Delhi, Punjab, and Haryana. R1 represents trajectories coming from Pakistan and Afghanistan. R2 represents trajectories originating from Iran. R3 is representative of all trajectories beyond these, including a portion of South Africa, the Mediterranean Sea, and Turkey. The seasonal clusters for the winter, spring, summer, and monsoon of the year 2017 and winter and spring of the year 2018 are shown in Fig. S1 in the Supplement. The re-clustering is shown in Fig. 1. All the chemical speciation data from ACSM and size distribution data from SMPS were then categorized as per the classification discussed above and used in the following discussion.

Aerosol ageing estimation
The NO x , toluene, and benzene data inventory for the entire campaign was taken from the CPCB for RK Puram, whenever available. The data were not available for branch B. To determine the photochemical ageing of aerosols, toluene and benzene concentrations were used to calculate the life (in hours), as per Nault et al. (2018) in Eq. (4): where [OH] = 1.5×10 6 molecules cm −3 (Nault et al., 2018), k toluene = 2.3 × 10 −12 exp(−190/T ), and k benzene = 1.8 × 10 −12 exp(340/T ) (Atkinson et al., 2006) are the rate constants for each aromatic compound, toluene i (o) = 1.85 and benzene i (o) = 2.31. [OH] is not constant and varies considerably temporally and spatially, but due to the unavailability of data of its variation for Delhi, it was assumed constant for ageing calculation.

Results and discussion
The HYSPLIT analysis revealed that the north-western direction is the most dominant direction, which is representative of SA air masses, and within it, R1 is the most dominant, indicating that overall, the emissions from Pakistan and Afghanistan and the sources en route govern Delhi's aerosol characteristics. However, the chemical signatures were potentially different for the various clusters, which explains the variation in aerosol properties with time. Due to the different nature of sources and pathways, aerosol properties vary, resulting in different hygroscopic properties and CCN forming potential. These aspects are discussed in the following sections.

Introduction to characteristics and sources of air masses
Out of the three main branches, the SA branch was the most anthropogenically contaminated, followed by BB and AS branches, as indicated by the mean NRPM 1 mass concentrations: 125.2 ± 91.6, 45.9 ± 23.3, and 32.5 ± 20.6 µg m −3 , respectively (Fig. 2a). The total NRPM 1 loading for the SA branch followed the sequence: NRPM 1(L) < NRPM 1(R3) < NRPM 1(R1) < NRPM 1(R2) . Amongst the SA branches, L was associated with the lowest organic (52.8 ± 40.6 µg m −3 ) and inorganic (42.1 ± 33.1 µg m −3 ) content, while R2 had the maximum organic (85.4 ± 59.8 µg m −3 ) and inorganic (57.9±47.5 µg m −3 ) content. A summary of the overall characteristics is given in Table 1. The prominent sources for the SA air mass include metal processing industries (Haryana and Delhi NCR), coke and petroleum refining (Punjab), thermal power plants (Pakistan, Punjab, and NCR Delhi), agricultural-residue burning (Punjab and Haryana), soil dust (Pakistan, Punjab) (Jaiprakash et al., 2017), and coal mines in Pakistan, where non-ideal burning of (NH 4 ) 2 SO 4 occurs (Chakraborty et al., 2015). When comparing the BB branches, total NRPM 1 for B (41.9 ± 20.8 µg m −3 ) was slightly less than that for B.reg (47.3±25.4 µg m −3 ), and this can be attributed to the fact that the B.reg air mass does not travel over water (originates adjacent to the coastline) but is subject to its influence, while the B air mass travels over water and is therefore cleaner. B.reg had a slightly higher inorganic and organic content than that of B.

PM 1 chemical composition of different air masses
Mass closure between SMPS size distribution data and the sum of ACSM species together with BC was achieved (R 2 =  0.83) as detailed in our parallel paper (Gani et al., 2019). The NRPM 1 species (NH + 4 , Cl − , NO − 3 , SO 2− 4 , POA, and OOA) and BC varied significantly for the different air masses, both in terms of the mass of species (Fig. 2c) and the diurnal patterns (Fig. 3), leading to different aerosol chemistry and chemical reactions. A summary of the average mass of each species for all air masses is detailed in Table 1. In brief, both POA and OOA, followed by NO − 3 , SO 2− 4 , and Cl − , dominated the PM composition for the SA air mass, while OOA, followed by SO 2− 4 and OOA, was dominant for BB and AS air masses. High chloride was a special feature of the SA air mass which was not apparent in the other two branches.
NH + 4 was assumed to be the dominant cation based on high aerosol neutralization ratio (ANR) values (mean values ranging from 0.95 to 0.85). The ANR is defined as the normalized ratio of the measured NH + 4 concentration to the NH + 4 concentration needed for full neutralization of the anions and calculated as per Eq. (5) (Zhang et al., 2007): Detailed ANR values are given in Table S1 in the Supplement. ANR values revealed that while AS, B, and L branches were completely neutralized, B.reg, R1, R2, and R3 were only partly neutralized, indicating that minor components of sulfate, chloride, and nitrate may be bound to non-volatile salts such as NaNO 3 , NaCl, or Na 2 SO 4 or are associated with organics as organosulfates, organochlorides, or organonitrates, evidence for which is shown in a previous DAS study   (Fig. S6). In all cases, an increase in r 2 for combined NO − 3 +Cl − as compared to individual ions indicates that both HNO 3 and HCl were synchronously neutralized by NH 3 . If Cl − and NO − 3 are present in the fine mode, they are expected to be in the form of their respective ammonium salts  (Harrison and Pio, 1983). Thus, the dominating salts were (NH 4 ) 2 SO 4 for AS, (NH 4 ) 2 SO 4 and NH 4 NO 3 for the BB air mass, and NH 4 Cl for SA and its sub-branches.
The organic speciation revealed that AS organics were BBOA; BB (both B and B.reg) organics were a mix of the three; and SA organics were BBOA, wherein L, R1, and R2 organics were BBOA while R3 organics were both HOA and BBOA.
The NO x emissions (in µg m −3 ) for the SA air mass (96.88 ± 127.22) were the highest, followed by those in BB (38.30 ± 64.79) and then those in the AS (36.71 ± 68.13) air mass. BB is representative of B.reg only, as NO x data for the B air mass were not available. NO x data for AS were also scarce. The SA aerosols exhibited less ageing (4.38 ± 4.49 h) compared to B.reg. (11.58 ± 3.45 h), but both were representative of aged aerosols. Ageing was not calculated for AS due to very little data availability.
To determine the presence of biomass burning and traffic emissions, BCwb (wood-burning component) and BCff (traffic component) for all air masses were determined, based on aethalometer data, as per Sandradewi et al. (2008). The contributions of BCwb and BCff are summarized in Table S4. Since fossil fuel sources are active year-round, there was a strong presence of BCff, ranging from 70 % to 86 %. However, biomass burning is only active during certain specific times for short durations and is very prominent in the northwestern direction for the SA air masses. It was observed that the more distant air masses exhibited a higher BCwb contribution compared to those originating within proximity. Hence, while L was associated with 13.9 % BCwb, R3 exhibited 29.2 % BCwb. The BCwb contribution for A and BB air masses was 21 %. It can thus be concluded that both biomass burning and traffic emissions are important sources contributing to the chemical composition of the various air masses.

The southern Asian air mass
This air mass ranked highest in NH + 4 concentration compared to other branches. Locally (i.e. for L), the source for NH + 4 may be attributed to NH 3 gas from the nearby agricultural fields of the Indian Agricultural Research Institute (IARI) (Sharma et al., 2014). For R1, R2, and R3, sharp spikes in early morning hours seen in the diurnal patterns of NH + 4 indicate its formation from ammonia as a result of industrial exhaust of untreated ammonia. This is because its diurnal variation is very similar to the diurnal of NH 3 emissions of an industrial origin (Wang et al., 2015). A very prominent feature of the SA that made it distinct from the other two air masses was the presence of high chloride ). High Cl − in the SA branch can be attributed to several factors: (a) Khewra salt mines in Pakistan that might contribute to high Cl − in other branches compared to L; (b) locally, plastic burning, refuse burning, and soil dispersion; (c) biomass burning, which is a very prominent feature of the SA branch, as indicated by f 57 and f 60 values and also by a large number of fire counts from MODIS fire-count data , dominantly in Punjab, Haryana, and a few places in Pakistan; (e) coal-based thermal power plants in Delhi, Punjab, Haryana, and Pakistan; and (f) small-and medium-scale metal processing industries in Delhi, Punjab, and Haryana, where HCl is used in a pickling process of hot and cold rolling of steel sheets and acid recovery from fume generation is not practised (Jaiprakash et al., 2017).
As far as the increase in chloride with the increasing length of trajectories is concerned, the most plausible explanation is biomass burning. It is pointed out in Sect. 3.2 that BCwb contribution increases as the air mass trajectories become distant, a feature similar to chloride emissions. There was a marked similarity in the diurnal patterns of NH + 4 and Cl − ions such that the sequence of [Cl − ] for SA sub-branches was also valid for [NH + 4 ], indicating the formation of NH 4 Cl. NH 4 Cl may also be emitted directly from cement plants (Cheney et al., 1983) in Punjab. The equilibrium constant for Cl − is more sensitive to ambient temperature than NO − 3 , as a result of which, during the daytime, a large amount of NH 4 Cl dissociates to form NH 3 and HCl if the temperature exceeds 10 • C (Kaneyasu et al., 1999). The diurnal patterns of both [NH + 4 ] and [Cl − ] exhibited a sharp decrease after 08:00 LT (all times listed in the paper are in local time, UTC+5:30) in the morning, which is obvious, since NH + 4 for the SA air mass is mostly associated with Cl − . At the same time  showed an increase between 09:00 and 10:00 and then started decreasing, but the rate of decrease was lower than NH 4 Cl. This is expected, as NH 4 NO 3 is more stable than NH 4 Cl (Kaneyasu et al., 1999). During winter, the ambient temperature drops slightly below 10 • C in the morning hours, increases sharply after 08:00 to reach a maximum at 14:00, and then again starts decreasing and reaches around 10 • C at midnight (Gani et al., 2019). Since during winters, the air mass comes mostly from the north-western direction of the SA air mass, it is evident that the formation and dissociation of NH 4 Cl were governed by the ambient temperature at the receptor site.
The reduction in chloride concentrations at midday can also be attributed to sulfate substitution mechanism, when sulfate formation enhances, and was also marked by the ratio R[SO 2− 4 ] ratio greater than 2 for L, R1, R2, and R3. This is valid especially for L, wherein [SO 2− 4 ] increased significantly. However, the diurnal patterns of [SO 2− 4 ] and [NH + 4 ] did not resemble each other, indicating that (NH 4 ) 2 SO 4 may be present in small amounts, but primarily SO 2− 4 is associated elsewhere. Hence, SO 2− 4 in combined form can be expressed in two ways: (a) small amounts of (NH 4 ) 2 SO 4 and (b) mostly in combination with K + . Thus, the sharp jump in [SO 2− 4 ] in locally originated air masses in the late morning and afternoon hours may be attributed to SO 2 emissions. SO 2 emissions in India are primarily attributed to power generation plants that make use of coal combustion as the chief source (Reddy and Venkataraman, 2002), followed by transportation. Such coal-based power plants are located in the IGP, with a high concentration in Haryana. SO 2 dissolves readily in water and can form sulfite ion, which in the presence of ozone can form sulfate ion (Erickson et al., 1977). H 2 SO 4 formed from the reaction of SO 2 and ozone can react with NH 3 to form NH 4 HSO 4 , which combines with NH 3 again to form (NH 4 ) 2 SO 4 (Stelson and Seinfeld, 1982;Seinfeld, 1986). Since the ozone spiked during the daytime (10:00-16:00), more sulfate formation was seen when ozone was maximum. The diurnal variation for ozone is explained in Gaur et al. (2014) for Kanpur, where a spike in ozone levels was seen during 10:00-16:00. A peak in sulfate concentration was also previously observed for the foggy period in Kanpur at 10:00 due to the resumption of photochemical activity after fog dissipation (Chakraborty et al., 2015). SO 2− 4 in SA branches may also combine with K + , as K + is produced in biomass burning. Evidence for the presence of K + along with SO 2− 4 in the accumulation mode is reported in Fuzzi et al. (2007). SO 2− 4 and NO − 3 emissions may also be associated with secondary formation for R1, R2, and R3 due to industrial emissions from metal product manufacturing industries in Punjab and Haryana; large-scale manufacturing of porcelain insulators; switchgear in Islamabad (Jaiprakash et al., 2017); and steel rolling mills in Iran, Iraq, and Turkey and Punjab.
The NO − 3 levels were very high for SA. The high NO − 3 in SA can be explained by the non-ideal burning of NH 4 NO 3 and NO x emissions due to mining equipment in the coal mines in Pakistan, leading to high NO − 3 formation (Chakraborty et al., 2015). It is mentioned in Sect. 3.2 that NO x emissions in the SA branch are very high. NH + 4 neutralizes NO − 3 simultaneously with Cl − ; however, the correlation of [NO − 3 ] with [ns-NH + 4 ] is moderate for SA. Therefore, it is expected that NO − 3 might be associated with K + and Na + , since biomass burning results in K + and Na + emissions (Fuzzi et al., 2007). K + and Na + exhibit a high affinity for nitrate during neutralization reactions, thus aiding in particulate nitrate formation (Bi et al., 2011). This is in addition to other nitrate sources that are discussed above along with SO 2− 4 sources for SA. The BC concentrations were highest for SA, followed by those in AS and then those in BB. The biomass burning in the SA air mass could be a major source of BC besides power plants, cement plants, local traffic, and industries. The POA emissions for SA followed the order BC (L) < BC (R3) < BC (R1) < BC (R2) . The spikes during the early morning hours and nighttime of the POA diurnal profile may be attributed to lower boundary layer heights during the two periods. BC and POA were well correlated (r 2 = 0.77) for R3, indicating that they come from primary emissions. The diurnal profiles for all branches were similar, which shows a decline as the day proceeds, followed by an increase as the night proceeds. OOA was present significantly in all the three branches but is maximum for SA. Its diurnal variation resembled that of NO − 3 (r 2 = 0.78), indicative of its semi-volatile nature. OOA and NO − 3 correlations were strongest for L (r 2 = 0.91), followed by R2 (r 2 = 0.85), R3 (r 2 = 0.81) and R1 (r 2 = 0.75).

The Bay of Bengal air mass
[Cl − ] was lower in this air mass compared to that in SA. For both of the BB branches, fossil fuel combustion was the most likely source of Cl − , as fossil fuel emissions dominate the IGP. Cl − for BB was not correlated with ns-NH + 4 and may be present in the form of methyl chloride, methylene chloride, carbon tetrachloride, and tetrachloroethene (Ho et al., 2004).
Fossil fuel combustion from coal plants along the IGP can be explained as a common source for both SO 2− 4 and NH + 4 ions, leading to (NH 4 ) 2 SO 4 formation (r 2 This was also seen in the diurnal profiles of [SO 2− 4 ] and [NH + 4 ], both of which exhibited a sharp spike in the early morning hours between 10:00 and 16:00. SO 2 emissions, as explained for SA via photochemical oxidation by O 3 in combination with NH 3 , can lead to (NH 4 ) 2 SO 4 formation. For the B.reg branch, the diurnal profiles of NH + 4 and SO 2− 4 exhibited double spikes (M pattern, which is a typical feature of NO x profile for traffic emissions) during heavy traffic hours (06:00-08:00 and around 16:00-19:00), indicating (NH 4 ) 2 SO 4 formation. NO from automobile exhaust can also form NH 3 in three ways catalytic convertors (Gandhi and Shelf, 1991), which in combination with SO 2 , formed due to pyrolysis of sulfide fuels and subsequent oxidation, can lead to (NH 4 ) 2 SO 4 formation. The correlation of [NO − 3 ] with [ns-NH + 4 ] was appreciably high for both B and B.reg. For both these branches, the fossil fuel combustion resulting in NO − 3 emissions in combination with NH + 4 (Rajput et al., 2015;Pan et al., 2016) can lead to NH 4 NO 3 formation. This is also evident from the diurnal profile of NO − 3 that shows a very similar pattern to NH + 4 and is expected to be in the form of NH 4 NO 3 . The diurnal profiles of NO − 3 show a decline as the temperature increases during the day and NO − 3 converts back to HNO 3 due to its semivolatile nature.
The BC concentration in BB air masses was considerably lower than in the SA air masses. The missing points in the diurnal variability in BC for the B region are on account of the unavailability of aethalometer data. The BC in the IGP can be emitted from industries (as for B), traffic (as for B.reg), and natural sources (Derwent et al., 2001). For B and B.reg, the B.reg branch was subjected to a longer duration of anthropogenic influence compared to B, which also spent considerable time on the water; hence after the early morning hours, when the various fresh emissions start increasing, the magnitude of POA for B.reg exceeds B. However, POA for BB was very low compared to SA. For B.reg, the spike in OOA during daytime hours was very similar to that of odd oxygen (O 3 + NO 2 ) for Delhi. The O x profile for Delhi is given in Tiwari et al. (2015), indicating its production by local photochemistry despite the increase in boundary layer height in the afternoon. The similarity in OOA and O x diurnal profiles was also noted for Kanpur (Chakraborty et al., 2016).

The Arabian Sea air mass
Chloride amounts were very low for AS compared to that in SA. Biomass burning as indicated by f 57 and f 60 measurements seems to be the main Cl − contributor to AS and might be associated with K + , which is also emitted along with it.
Similar to the cases of the L branch in the case of the SA air mass and the BB air mass, the power stations in Gujarat and Rajasthan lead to SO 2 emissions. Since the power plants in this region, over which the AS air mass traverses, are relatively small in number, the SO 2 concentration is much lower compared to that in the BB air masses. SO 2 emissions subsequently lead to (NH 4 ) 2 SO 4 formation, which was the main salt present in this branch and was also evident from the high correlation between the two ions. (NH 4 ) 2 SO 4 may be formed due to emissions from both power plants and traffic (similar to B.reg). Traffic emissions can be understood from the M pattern in diurnal profiles of NH + 4 and SO 2− 4 , though the variation was not very pronounced and might be suppressed due to power plant emissions. The traffic signal was more clearly implied by the diurnal profile of (NH 4 ) 2 SO 4 for AS, as seen in Fig. S2. The correlation of [NO − 3 ] with [ns-NH + 4 ] was very poor for AS, indicating that NO − 3 might be associated with K + and Na + , similar to B.reg.
Both BC and POA for the AS air mass were less than in the SA air mass. However, compared to BB, BC was slightly higher and POA was comparable. BC was likely of an industrial origin. The POA diurnal profile was similar to the other air masses. Similar to B.reg, the OOA diurnal pattern resembled that of odd oxygen, where the odd-oxygen profile is reported in Tiwari et al. (2015).
Thus, the direct emission sources and the precursors (SO x , NO x , NH 3 , O 3 , and O x ) that lead to particulate matter formation strongly impacted the chemical properties of aerosols. The chemical properties of aerosol also impacted the hygroscopicity of aerosols, as is discussed in the following section.

Impact of chemical composition on the hygroscopicity of air masses
This study provides the first long-term estimation of aerosol hygroscopicity in the PM 1 regime. The mean κ was approximately the same for all the air masses, which is ∼ 0.3 (0.32 ± 0.06 for AS, 0.31 ± 0.06 for BB, and 0.32 ± 0.10 for SA) and in line with the global average value of 0.27 ± 0.21 for continental aerosols Petters and Kreidenweis, 2007;Pöschl et al., 2009;Pringle et al., 2010). Including BC in κ calculations leads to a difference of 10 % in κ on average, shifting the mean κ of 0.32 to 0.29. The BC mass fraction and volume fractions were 10 % and 9 %, respectively. Thus, the change in κ due to the introduction of BC was not significant. κ varied from 0.13 to 0.77, and there was a difference in the diurnal variation in the hygroscopicity parameter for the various air masses (Fig. 4). A similar finding was observed in China, with a mean κ of 0.3, varying in the range 0.1-0.5 . Recently, κ of 0.42 ± 0.07 was also reported for PM 2.5 for Delhi based on beta attenuation monitor (BAM) measurements of PM 2.5 (Wang and Chen, 2019), using an indirect method in the absence of direct measurements. Thus, the dependence of κ on size cannot be underestimated for Delhi and should be dealt with in the future.
Variation in κ with size has been seen in several places in the world; e.g. at 97 % RH, mass growth factors of 6.95 and 9.78 were reported for the size ranges 0.53-1.6 and 1.6-5.1 µm, respectively, on the Slovenian coast (Turšič et al., 2006). Aitken mode κ was 0.25, while the accumulation mode κ was 0.45 for Beijing . An increase in κ for the higher size regime points to the fact that the organic fraction is higher in the smaller size range, while the inorganic fraction increases substantially with size.
In the PM 1 size range, while κ of 0.1 indicates secondary organic aerosol, κ varied from 0.01 to 0.8 for biomassburning aerosols in lab studies (Petters et al., 2009). κ varied from 0.15 to 0.25, with lower values (around 0.16) being observed during the night, when biomass-burning particles prevailed during wintertime in Athens, Greece (Psichoudaki et al., 2018). Thus, κ values for Delhi can represent both secondary formation and biomass burning. This is true for Delhi, which had both POA and OOA in all the air masses, while BBOA was present in AS and SA air masses, as detailed in the preceding sections on chemical properties.
An important observation for all branches is that when the inorganic volume fraction (of dominant salt) increased (Fig. S2), during the times when κ was high (Fig. 4), or when the organic volume fraction decreased (Fig. S2), a dip in D c (Fig. 6) was seen, implying that a larger size regime was available for activation. The diurnal variation in κ (Fig. 4) more strongly followed the diurnal pattern of the dominant inorganic salts for a cluster (Fig. S2), since the hygroscopicity parameters for inorganic salts are considerably higher than those of organics. Pearson correlation coefficient (r) values between κ and the salt volume fractions revealed that the diurnal patterns of κ were governed dominantly by volume fractions of (NH 4 ) 2 SO 4 (r: 0.85) for AS, moderately by (NH 4 ) 2 SO 4 (r: 0.55) and NH 4 NO 3 (r: 0.49) for BB, and dominantly by NH 4 Cl for SA air masses. For the two BB branches, κ of the B branch was governed dominantly by (NH 4 ) 2 SO 4 (r: 0.78) and moderately by (NH 4 ) 2 SO 4 (r: 0.57) and NH 4 NO 3 (r: 0.54) for B.reg. For SA air masses, κ of R1, R2, and R3 was governed dominantly by NH 4 Cl (r values of 0.71, 0.89, and 0.95) and jointly by NH 4 Cl (r: 0.65) and NH 4 NO 3 (r: 0.73) for L.
High volume fractions of (NH 4 ) 2 SO 4 and NH 4 NO 3 may be attributed to SO x , NO x , and NH 3 emissions due to power plant emissions and traffic. In the SA sub-branches (Fig. 4), the spike in κ during the early morning (07:00-08:00) exhibited the sequence κ (R3) > κ (R2) > κ (R1) >= κ (L) , and the lower spike in late evening (18:00-22:00) exhibited the sequence R3 > R2 < R1 < L and is attributed to NH 4 Cl formation. This implies that during the morning, R3 aerosols were most hygroscopic, while L aerosols were least hygroscopic, while after 09:00, L aerosols were most hygroscopic and R3 aerosols were the least hygroscopic. The flatter curve of κ can be attributed to two factors: (a) the chloride contribution of distant trajectories decreased very steeply with time compared to the local emissions, and (b) κ of L was also supplemented substantially by NH 4 NO 3 . Thus, the source activities by which the chemical properties of aerosols are shaped impacted the hygroscopicity parameter tremendously. This consequently impacted the size regime of aerosols available for activation and is discussed in the following section.    Table S2.
High values of N CCN for Delhi are consistent with other polluted regions in the world. The relevant statistics for two highly polluted sites, namely Beijing and Kanpur, are presented by Gunthe et al. (2011) and, respectively. N CCN was 7660 ± 3460 and 900-27 000 (in cm −3 ) at 0.46 % SS and in the range 0.18 %-0.6 % SS, respectively. The high N CCN is consistent with the high N CN . Correspondingly, N CN was 16 800 ± 9100 and ∼ 50 000 (cm −3 ). Even though high number concentrations of CCN and condensation nuclei (CN) have been reported, the a f was not seen to be so high. The a f in Beijing was 0.54 ± 0.23 and 0.66 ± 0.23 at 0.46 % SS and 0.86 % SS, respectively. For Kanpur, a f was reported as ∼ 0.018-0.54 for 0.18 %-0.60 % SS. However for Delhi, the a f ranged from  (2019), which states that for Delhi, activation of a 0.1 µm particle requires SS ∼ 0.18 ± 0.015 %, compared with ∼ 0.3 % for Beijing; 0.28 %-0.31 % for Asia, Africa, and South America; and ∼ 0.22 % for Europe and North America. The high activated fractions of aerosol can impact the precipitation patterns in Delhi and may be responsible for the short, intense precipitation events and decrease in overall rainfall. However, no study to date has validated this growing trend with CCN measurements or estimates, and this needs to be investigated in the future. The a f and N CCN for all air masses increased as expected with an increase in supersaturation. The variation in CCN and the activated fraction with SS are shown in Figs. 5 and S3. The figures clearly show that even though N CCN for SA was far greater compared to BB and AS, the activated fractions were fairly close. B.reg had higher N CCN but a close a f compared to B. Similarly, L had the lowest N CCN among all SA branches but the highest a f . These features elucidate that many factors are at play and impact N CCN and a f differently. To determine the governing parameters impacting both N CCN and a f , the diurnal patterns of N CCN , a f , and D c at 0.1 %, 0.4 %, and 0.8 % SS; N CN , N Aitken , and N Accumulation ; κ; and geometric mean diameter (GMD) are shown in Figs. 6, S4, and S5.
The N CCN curve for SA showed a sharp diurnal feature which was not as prominent for the other two. However, at noon, at 0.4 % and 0.8 % SS, while N CCN dipped for SA, it rose for AS and BB (Fig. 6). Furthermore, with the increase in supersaturation, the dip in CCN of SA increased such that (a) at 0.1 % SS, SA CCN > AS CCN , BB CCN ; (b) at The explanation for this observation is that at 0.1 % SS, the CCN was governed by the accumulation mode, but at 0.4 % and 0.8 % SS, it was governed more by the Aitken mode. It is the Aitken mode that dominantly governs the total number concentration, and hence it can be said that CCN is governed by CN at a higher SS and by the accumulation mode at low SS. This is because as supersaturation increases, D c increases. At low SS, D c is high (D c > 100 nm), almost always at a 0.1 % SS, for all branches, as shown in Fig. 6d; hence the size distribution that is integrated to get the CCN involves the accumulation mode only. At 0.4 % SS, the D c was around 40-47 nm; therefore a considerable fraction of the Aitken mode and accumulation mode was available for activation. At 0.8 % SS, the D c was ∼ 25-30 nm; therefore, the contribution of the Aitken mode further increases and the accumulation mode is also available as usual. These findings are also true for BB branches (B and B.reg) and SA branches (L, R1, R2, and R3), as shown in Fig. S4. It is the low value of D c relative to other places that is responsible for high CCN. As explained in Sect. 3.4, D c is largely associated with κ. The D c at other places such as Kanpur varied from 50 to 200 nm, for SS ranging from 0.18 to 0.60 , compared to which D c for Delhi is lower (17-142 nm, for SS ranging from 0.1 % to 0.8 %), implying that a larger regime is available for activation. It is pertinent to mention here that the dip in CCN for SA at midday and the peak at the same time for AS and BB can be attributed to the following. (a) The dip in N Accumulation for SA at midday was much more prominent compared to AS and BB. (b) At midday, the N Aitken for SA also decreased, while it increased for AS and BB. Thus, the dip in N CCN was strengthened by the simultaneous dip of both the Aitken and accumulation modes, while the peak in N CCN for AS and BB was a manifestation of the dominant peak in Aitken modes. Similar features were also exhibited for sub-branches of BB and SA. A deeper insight reveals that the dip in number concentration at midday for SA was most aptly seen in the diurnal pattern of POA (which is the most dominant NRPM 1 species) and, to quite a good extent, in other NRPM 1 species barring the SO 2− 4 ion. Similarly, the peak in the Aitken mode for AS and BB can be attributed to NH + 4 , SO 2− 4 , and OOA concentrations (the dominating species in the respective branches). Thus, the source activities and trajectory pathways impact CCN concentration at the receptor site.
The diurnal pattern for a f showed a dip in the midday hours for all the air masses even though N CCN at midday for AS and BB peaked during midday hours. The time of dip in a f (more prominent at 0.4 % and 0.8 % SS) occurred earlier for AS compared to that for BB and SA (Fig. 6). This was governed by the time of dip in the GMD of the three branches. The GMD diurnal variation was very similar to a f , and the r 2 values between GMD and a f listed in Table S3 also point to the same result. Thus, even though a dip in D c should correspond with an increase in a f and vice versa, this does not usually happen, as the change in D c is less compared to the shift in size distribution such that not only the highest number concentration values but also the diameter at which this occurs change, thereby changing the number available for activation. For example, in Fig. 7 two different size distributions for the AS air mass at 02:00 and 11:00 are compared. At 02:00 and 11:00, the following characteristics were noted: (a) GMD 02:00 : 83.78 ± 16.58 nm (GSD -geometric standard deviation), which is considerably higher than GMD 11:00 : 47.18 ± 13.01 nm; (b) GSD 02:00 : 1.69 ± 0.13 nm, which is nearly the same as GSD 11:00 : 1.62 ± 0.14 nm; (c) D c02:00 : 44.45 ± 3.83 nm, Figure 7. Comparison of two size distribution profiles at different times of day for AS branch at 02:00 and 11:00. which is slightly higher than D c11:00 : 43.48 ± 2.28 nm (at 0.4 % SS); (d) CN 02:00 : 15 849 ± 9269 cm −3 , which is lower than CN 11:00 : 25 873 ± 9840 cm −3 ; (e) CN_Aitken 02:00 : 9819 ± 5945 cm −3 , which is lower than CN_Aitken 11:00 : 22 376 ± 9693 cm −3 ; (f) CN_Accumulation 02:00 : 5729 ± 3684 cm −3 , which is higher than CN_Accumulation 11:00 : 3311 ± 1569 cm −3 ; (g) a f02:00 : 0.75 ± 0.12, which is considerably higher than a f11:00 : 0.38 ± 0.17 at 0.4 % SS; and (h) CCN 02:00 : 11 595 ± 6710 cm −3 , which is higher than CCN 11:00 : 8946 ± 3899 cm −3 at 0.4 % SS. The a f at 02:00 is higher than that at 11:00. A D c decrease should correspond to a CCN increase, but the magnitude of D c = 0.97 nm is very small. The decrease in GMD ( GMD = 36.6 nm), on the other hand, is very high, with negligible changes in GSD. At 02:00, since D c was considerably less than GSD, most of the particles are counted for activation. At 11:00, since D c and GMD were very close, nearly 50 % of particles are not available for activation. The very high Aitken mode concentration at 11:00 (higher than that at 02:00) was not available for activation at both the times. Thus, N CCN here is governed by the accumulation mode, which was higher at 02:00, thus making N CCN higher. In this scenario, N CCN and a f go hand in hand. However, there also exists a second possibility; e.g. for the SA branch, N CCN , at 0.4 % SS, at 08:00 is higher than that at 05:00, while the a f at 08:00 is lower than that at 05:00 (Fig. 6). N CCN was governed at these times by the Aitken mode, which was higher at 08:00, while a f is governed by the GMD, which was lower at 08:00. Therefore, it is established that CCN is governed by CN (dominantly by either the Aitken or accumulation mode as the case may be), while the a f is governed by GMD.
At this junction, it is also pertinent to mention how chemical dispersion and parameters governing CCN are interconnected. The standard deviation of κ (σ (κ)) around κ is often used as an estimate of the degree of heterogeneity (chemical dispersion) of particles (Psichoudaki et al., 2018;Lance et al., 2013). The chemical dispersion for all air masses and their sub-branches is shown in Figs. 6, S4, and S5. The chemical dispersion for the SA air mass during the early hours (06:00-08:00) coincided with chloride emissions, and late at night, after 20:00, it coincided with POA and OOA emissions. During the time of high chloride emissions, κ also peaked, since inorganics are associated with high hygroscopicity, while during the late hours, κ dropped due to an increase in organics associated with low hygroscopicity. The diurnal patterns of the activated fraction, GMD, and chemical dispersion were also similar. This implies that higher heterogeneity shifts GMD towards a high value, thereby increasing the available regime for activation, and vice versa. There was no discernible pattern noted for the other air masses.
It is hereby established that aerosol physical and chemical properties, and their time evolution, are tightly linked with each other. The indirect impact of chemical composition on CCN and a f is manifested in two ways: (a) NRPM 1 species impact the diurnal patterns of Aitken and accumulation modes, which in turn impact CCN, and (b) NRPM 1 species impact κ, hence subsequently CCN, by impacting the size regime available for activation.

Conclusion
Long-term measurements of NRPM 1 species and size distribution data were carried out in New Delhi. The air masses originated from SA (L, R1, R2, and R3), BB (B and B.reg), and AS. κ was estimated using the mixing rule, and the bulk κ was assumed for the entire size distribution. Using κ and size distribution data, CCN estimates were obtained. The SA air mass was the most contaminated air mass, followed by BB and then AS. This resulted in higher NRPM 1 , N CN (both Aitken and accumulation modes), N CCN , and a f for SA, followed by those in BB and then those in AS. The most dominant salts turned out to be (NH 4 ) 2 SO 4 for AS, (NH 4 ) 2 SO 4 and NH 4 NO 3 for BB, and NH 4 Cl for SA. The AS, B, and L branches were completely neutralized, while B.reg, R1, R2, and R3 were partially neutralized. The diurnal variations in NRPM 1 species were governed by source activities' aerosol precursors, like SO x , NO x , NH 3 , O 3 , and O x . The high PM 1 concentration for Delhi, which exceeds the National Ambient Air Quality Standards, can be mitigated only by controlling both the primary emissions and precursors. To address the situation justly, the following is lacking: (a) data listing measurements of PM 1 emissions from various industries in India and Asia, (b) a description of the chemical constituents of aerosol that are emitted, both qualitatively and quantitatively, and (c) a definition of emission limits and compliance with them.
The mean κ ∼ 0.3 was the same for all air masses, with the diurnal variation in κ governed by chemical species and, thus, source activities. The κ diurnal trends impacted the D c diurnal trend, which in turn affected the available regime for activation. The N CCN diurnal patterns were driven by the accumulation mode at a lower SS and Aitken mode with an increase in SS, depending upon D c , which decreases with an increase in SS. The D c obtained for Delhi was lower than that seen at other places in the IGP, for example, Kanpur. The activated fraction for Delhi was very high (0.71 ± 0.15 at 0.4 % SS for R1), with the means of activated fractions varying between 0.19 and 0.87, for SS varying from 0.1 % to 0.8 %, whereby their diurnal patterns were governed by GMD. A CCN measurement study with a CCN counter in the future can help verify the estimates, and a closure ratio may be determined. However, in the absence of long-term cloud condensation nuclei counter (CCNC) measurements, the importance of these findings cannot be neglected. These results can serve as valuable inputs to GCMs to better quantify precipitation. The high NRPM 1 loading and activated fractions are bound to significantly impact precipitation over Delhi, impact the aerosol radiation budget, and have indirect effects and need to be investigated thoroughly in the future.
These investigations may answer the short intense precipitation events occurring over Delhi and the decrease in the overall rainfall over the past half-century.
Author contributions. LHR, JSP, GH, and ZA designed the study. ZA, SG, and SB carried out the data collection. ZA carried out data processing and analysis. ZA and GH carried out the interpretation of the results. ZA wrote the paper and was assisted by SB, LHR, and GH in reviewing the paper.