Size-resolved Compositional Analysis and Source Apportionment of Submicron Aerosol during Lockdown Period Using HR-ToF-AMS

The size-resolved compositional analysis of non-refractory submicron aerosol (NR-PM 1 ) was conducted using the Aerodyne High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS) instrument over Pune, India during the COVID-19 lockdown period. The aerosol composition data shows the predominant presence of organics (Org) in the mass fraction followed by sulfate, ammonium, nitrate, and chloride during the pre-lockdown and lockdown periods. The size-resolved analysis showed the unimodal size distribution of organic and inorganic constituents with peaks at 550 nm, implying the dominant presence of mixed and aged aerosol species. The stoichiometric neutralization analysis showed the almost neutralized nature of submicron aerosol with an average aerosol neutralization ratio (ANR) of 0.8. The back trajectories, cluster analysis, and potential source contribution function (PSCF) showed the industrial belt located in the western part of the study location to be the potential source regions of NR-PM 1 . Positive matrix factorization (PMF) analyses have been applied to investigate the source apportionments of organic aerosols (OA). Four distinct OA factors, i.e., hydrocarbon-like OA (HOA), biomass burning OA (BBOA), low-volatile oxygenated OA (LVOOA), and semi-volatile oxygenated OA (SVOOA) were identified during the study period. Among these factors, HOA contributes nearly a quarter to the OA mass, and OOA accounted for nearly 60% of the total OA mass. The high-resolution positive matrix factorization (HR-PMF) analysis and the elemental ratios of H/C, O/C, and OM/OC showed distinct characteristics during different periods. The density of organic aerosol has been estimated using the elemental ratios and found to be 1.14, 1.28, and 1.35 respectively during the different lockdown periods, similar to 1.30 g cm –3 as mentioned in the literature. This study provides new insights into the chemical composition and source apportionment of the organic fraction of submicron aerosols for the first time over Pune using HR-ToF-AMS and HR-PMF.


INTRODUCTION
The outbreak of novel coronavirus has led to the death of millions of individuals around the world, and countries implemented strict restrictions (lockdown) to prevent the spread of the virus and preserve the health of the citizens (Dhaka et al., 2020;Agarwal et al., 2020). The majority of the industrial and other pollution-emitting establishments and activities, including vehicular movement, shopping malls, hotels, and construction work were closed, which caused a nearly complete cessation of anthropogenic emission of pollutants. The implemented lockdown has created a unique experimental situation in India to investigate the role of anthropogenic pollution sources on air quality, which is not possible in normal conditions.
The near-complete cessation of anthropogenic activities during the lockdown periods has greatly influenced air quality over different parts of the world, and the drastic reduction in pollution levels acted as the background or baseline pollution level (Venter et al., 2020;Bao and Zhang, 2020;Collivignarelli et al., 2020;Dantas et al., 2020;Dutheil et al., 2020;Wang et al., 2020;Jephcote et al., 2021). Pollutants like particulate matter (PM), carbon monoxide, nitrogen dioxide, sulfur dioxide, and other trace gases showed a significant decrease in the lockdown periods (Dhaka et al., 2020;Sharma et al., 2020;Jain and Sharma, 2020;Mahato et al., 2020). However, some species like ozone increased after lockdown restrictions were imposed, and the occurrence of several haze events was observed in the North China Plain (NCP) despite the inactivity of most of the anthropogenic primary pollutant sources (Huang et al., 2020). These variabilities were attributed to the stagnant meteorological conditions and enhanced oxidative capacity in the atmosphere. But, these phenomena raise the fundamental question, what is the role of anthropogenic sources on air quality? To date, very few studies have focused on the variability in the ultrafine particle compositions over the Indian region during the lockdown period (Dave et al., 2021).
Aerosols consist of a wide range of chemical species, and organic compounds are found to be the dominating fraction of the submicron aerosols . Organic aerosols (OA) are generally classified into primary organic aerosols (POA) and secondary organic aerosols (SOA). POA fraction refers to the aerosols directly emitted, and SOA is formed in the atmosphere via homogeneous, heterogeneous chemistry and other complex processes like photo-oxidation, nucleation, aqueous-phase reactions, and gas-to-particle partitioning of precursor gases (Kanakidou et al., 2005;Ervens et al., 2011;Hallquist et al., 2009). The knowledge of particle composition, size distributions, temporal variations, and the contribution of different source factors are needed to understand the fine particle chemistry in different environments.
In this study, real-time measurements of the chemical constituents of non-refractory particulate matter with ≤ 1 µm aerodynamic diameter (NR-PM1) were conducted during the pre-lockdown and lockdown periods using High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS). The temporal variation and mass size distribution of the aerosol species could give insight into the source sectors and the pollution level assuming lockdown period concentration as the baseline. The stoichiometric neutralization analysis has also been done to investigate the acidic nature of the submicron aerosols. The unit mass resolution (UMR), and high-resolution (HR) positive matrix factorization (PMF) techniques were employed to resolve the measured organic spectra into characteristic source profiles and evaluate their contribution. Back trajectories, cluster analysis, and potential source contribution function (PSCF) were also performed to investigate the potential source regions of NR-PM1. The salient features of the size-resolved compositional analysis of submicron aerosol (NR-PM1) with HR-ToF-AMS during the lockdown periods are presented in this manuscript.

Measurement Site
The measurements were conducted at the Indian Institute of Tropical Meteorology (IITM, 18.54°N, 73.81°E, 560 m asl), Pune as shown in Fig. S1. Pune is situated on the lee-side of the Western Ghats and is about 100 km east of the Arabian Sea. The site is surrounded by hillocks of variable heights (up to 200 m), forming a valley-like appearance. In the last few years, Pune has been experiencing tremendous growth in terms of urbanization and vehicular activity. The HR-ToF-AMS (Aerodyne Research, Inc., USA) was installed on the second floor of the institute's building at about 10 meters above the ground. The inlet was installed on the rooftop at about 12 meters above the ground. The measurement was conducted from 10 March to 17 May 2020, and a comprehensive description of the observational site and its meteorological conditions has been mentioned elsewhere (Ali et al., 2012(Ali et al., , 2016Safai et al., 2014).

High-Resolution Time-of-Flight Aerosol Mass Spectrometer (HR-ToF-AMS)
There are various methodologies like ion chromatography (IC), particle into liquid sampler (PILS), Monitor for AeRosols and Gases in Ambient air (MARGA), and aerosol mass spectrometry (AMS) are available to characterize the chemical constituents of aerosols (ten Brink et al., 2007;Schaap et al., 2011;Acharja et al., 2020). AMS excels over other measurement methods due to its fast, sensitive, versatile analytical technique providing size-resolved compositional analysis of major organic and inorganic constituents in real-time time resolution with lesser positive and negative artifacts (Jayne et al., 2000;Jimenez et al., 2003;Allan et al., 2003). The HR-ToF-AMS data were analyzed in IGOR (Wavemetrics, Inc., Lake Oswego) using the tools named SeQUential Igor data RetRiEvaL (SQUIRREL) and Peak Integration by Key Analysis (PIKA) (Alfarra et al., 2004;Drewnick et al., 2005). A number of calibration methods were carried out at regular intervals during the measurement period. For example, m/z calibration was performed daily to inspect the mass spectra, and single ion signal (SIS) calibration was performed every alternate day to inspect the strength in single ion signal. The velocity calibration has also been performed with ammonium nitrate of different diameters to perform the size distribution calibrations. A detailed description of the HR-ToF-AMS instrument, its operating principle, and calibration method is discussed elsewhere, and a brief description is given in the Supplementary Section Aiken et al., 2008;Ulbrich et al., 2009;Salcedo et al., 2010;Lanz et al., 2010).

Source Apportionment Using Positive Matrix Factorization (PMF)
In order to perform the source apportionment of OA mass spectra and separate total OA into different aerosol types, the positive matrix factorization (PMF) tool is used in this study (Paatero and Tapper, 1994;Ulbrich et al., 2009). PMF is a bilinear unmixing model that does not require a priori information of the source profile or their temporal contribution and has been identified as an appropriate receptor modeling technique to quantify source contributions Canagaratna et al., 2007;Aiken et al., 2008;DeCarlo et al., 2006). It performs deconvolution of mass spectra (MS) into the summation of products of positively constrained mass spectral profiles and their corresponding time series under the assumption that the mass spectral profiles remain constant in time.
The measured mass concentration (Xij) can be defined as per the following equation  ( ) where i, j, and k represent time, variable, and particular factor, respectively. Xij, gik, fkj, and eij are elements of the measured organic matrix, time series of factor contributions, factor profiles, and residual matrix respectively. The model uses a least square algorithm (sum of squares of model error normalized to measurement error) to iteratively minimize the fit parameter Q, which is defined as the sum of the squares of the ratio of residual (eij) to measurement uncertainties 2 1 1 where m is the number of samples (time series), and n is the number of variables (m/z) in the matrix. PMF can result in mathematically accurate solutions through several linear combinations of the factor profile and time series, also known as rotations. Details of the PMF source apportionment technique and its application to the AMS dataset can be found in Ulbrich et al. (2009) andCanagratna et al. (2007).

RESULTS AND DISCUSSION
The fundamental approach to this work is to understand how the anthropogenic sources influence aerosol loading, especially during a unique experimental situation of near-complete cessation of anthropogenic activities. This study characterizes concentrations of the organic and inorganic constituents of submicron aerosol monitored using the HR-ToF-AMS during the period from 10 March to 17 May 2020. The whole period under consideration includes pre-lockdown and lockdown periods and is classified into three categories based on the efficacy of the pollution sources.
The first period spread over 10-24 March 2020 is categorized as pre-lockdown (P), characterized by the regular presence of natural and anthropogenic emission sources. The second period spread over 25 March-3 May 2020 and is categorized as lockdown-I (L1), characterized by the complete shutdown of anthropogenic emission sources like construction work, industrial activity, traffic, and other commercial activities except for essential services. However, it is mentioned that data from 25 March-29 April is not available as the instrument could not be operated in the laboratory due to strict restrictions on movement. Data during L1 was available from 30 April to 03 May 2020, and this dataset reflects the non-availability of anthropogenic pollution sources. The third period of this study spread over 4-17 May 2020 and is classified as lockdown-II (L2), characterized by partial availability of pollution sources like coal and mineral production, agricultural activities, and agro-based packaging industries. The characterization of aerosols during these periods reflects the impact of anthropogenic activities on air quality.

Temporal and Diurnal Variation of Organic and Inorganic Constituents of NR-PM1
Mass concentrations of the organic (Org), and inorganic constituents of PM1 and PM2.5 during the pre-lockdown and lockdown periods are shown in Table 1, and a significant reduction in PM1 and PM2.5 constituents was observed during the lockdown periods. The mass fractions of NR-PM1 constituents are presented in Fig. 1, which shows that organic and SO4 2dominated the submicron particles, followed by NO3and NH4 + . Non-refractory Clwas present in a significantly low fraction over Pune. From the table, it can be seen that organic mass concentrations during the P, L1, and L2 periods were 10.2 ± 6.6, 4.1 ± 2.1, and 7.3 ± 3.3 µg m −3 , respectively. Organic mass fraction is ubiquitous in various atmospheric environments, and the presence of a high organic fraction of PM1 is consistent with those obtained in other parts of the world Jimenez et al., 2003;Cash et al., 2021;Petit et al., 2021;Via et al., 2022). The primary sources of organic aerosols (OA) include fossil fuel burning, domestic burning (cooking and heating), and uncontained burning of vegetation and agricultural waste (Kanakidou et al., 2005).
Organic mass concentrations decreased by about 60% and 29% during the L1 and L2 lockdown periods, respectively, compared to the pre-lockdown period. The average concentration of NH4 + , SO4 2-, and NO3also decreased by 18%, 22%, and 53% respectively, during the lockdown period L1, which could be due to an appreciable reduction in fossil fuel and biomass burning. The concentration of NH4 + , SO4 2-, and NO3during L2 was nearly equal to the values observed during the pre-lockdown period. The increase in L2 was due to the gradual resumption in services and sectors like coal and mineral production, agricultural activities, and agricultural-related packaging industries that are potential sources of fine aerosols .  The concentration of SO4 2always remained higher than NO3during the study periods. The major contribution from the thermal power plants of Mumbai (1580 MW) and Nashik (900 MW) located in the north to the northwest direction of the measurement site was the primary cause of higher sulfate than nitrate over Pune. The power sector remained operational throughout the pre-lockdown and lockdown periods and acted as the major anthropogenic pollution source sector. The nitrate concentration dropped significantly from 0.9 µg m −3 during pre-lockdown to 0.4 µg m −3 during L1. It moderately increased to 0.7 µg m −3 during L2, which is consistent with the temporal variation in NOx concentrations, as shown in Table 1. The NO3and NOx concentration decrease in lockdown is mainly due to the near-complete cessation of vehicular movement in L1. The concentration of NH4 + was 1.1 ± 0.5, 0.9 ± 0.6, and 1.2 ± 0.5 µg m −3 during the P, L1, and L2 periods respectively. The variation in the concentration of NH4 + is not drastic as organic and nitrate due to the unabated agro-activity of spraying fertilizer in the cultivated fields and increased soil ammonia emissions in the summertime. Non-refractory Clwas significantly low during the whole observational period with a minimum value of 0.1 ± 0.0 µg m -3 during L1 and a maximum value of 0.2 ± 0.1 µg m -3 during the pre-lockdown period.
Diurnal variations of the organic, SO4 2-, NO3 -, NH4 + , and Clduring pre-lockdown, L1, and L2 are shown in Fig. 1. Organic showed a peak in the morning hours during the pre-lockdown period, which is a primary feature of the urban environment and mainly associated with rush hour traffic emissions and secondary productions in the atmosphere. The relatively high photochemical production in the daytime may get compensated by the loss due to the dilution effect and increased boundary layer height. The night peak in organic concentration was reduced compared to the sharp peak in the morning. During L1, the organic concentration and the sharp morning peak significantly reduced and the nighttime peak at 20:00 hrs was the highest, that could be due to the nighttime chemistry and boundary layer dynamics. During L2, organics showed a morning peak that could be due to vehicular emission as some of the vehicles were on the roads due to essential services. The diurnal cycle of the secondary inorganic aerosols (sulfate, nitrate, ammonium, and chloride) has also been investigated in Fig. 2.
The concentration in SO4 2showed short diurnal variation but is not conspicuous. The diurnal variation of nitrate showed a higher concentration during office hours in the pre-lockdown period. It reduced significantly during the lockdown periods as the vehicular emissions were significantly reduced. The diurnal variation of ammonium and chloride showed no perturbation in concentrations. The monitored secondary aerosol species are in the 1 µm size range, and the temporal variations of these species are comparatively newly explored as several studies have explored the characterization and presence of secondary species mostly in the super-micron size range of 1-2.5 µm (Sun et al., 2020;Acharja et al., 2022).

Size Distributions of the Constituents of NR-PM1
The particle size distribution, chemical composition, and mass concentrations have a vital impact on the atmosphere and human health due to their size distributions (Seinfeld et al., 1998;Davidson et al., 2005). The HR-ToF-AMS measurements were used for the first time over Pune to assess the size-resolved compositions of aerosols in real-time. This yields information on the composition of aerosols with mass concentrations in the Aitken-mode (Da ~10-100 nm) and the accumulation mode (Da ~100-1000 nm) (Drewnick et al., 2005;DeCarlo et al., 2006).
The size-resolved measurements of NR-PM1 during the pre-lockdown and lockdown periods are shown in terms of aerosol aerodynamic diameter (dM/dlogDva) in Fig. 2. The mass size distribution peak of fine aerosols at the smaller aerosol size (Dva ≤ 200 nm) suggests the dominant presence of freshly formed aerosols over the region. A peak at a higher size range (500-600 nm) would refer to the aged and mixed regional organics particles . The monitored aerosols during the pre-lockdown and lockdown periods showed mass concentrations peak at the vacuum aerodynamic diameter (Dva) of nearly 550 nm, indicating the presence of more internally mixed, aged aerosols and less freshly emitted over Pune. The broader organic peak compared to the inorganic species illustrates the dominant nature of organic species in NR-PM1, which is consistent with what is mentioned in Section 3.1.

Aerosol Neutralization Ratio (ANR)
The composition of atmospheric particles plays a significant role since it affects the nature of the acidity and hygroscopicity of aerosols. Direct measurements of the neutralization nature of aerosols are challenging and rare. The neutralization nature of aerosols is derived using the stoichiometric neutralization analysis. The aerosol neutralization ratio (ANR) of aerosols gives insight into the composition and characteristics of the submicron particles Acharja et al., 2022). The ANR is defined as follows ( ) ( ) where NH4 + (m), SO4 2-(m), NO3 -(m), and Cl -(m) are the measured mass concentrations of ammonium, sulfate, nitrate, and chloride, respectively. The average ANR values during the prelockdown and lockdown periods were estimated using the aerosol composition data. During P, L1, and L2 the ANR was 0.75 ± 0.17, 0.87 ± 0.17, and 0.82 ± 0.16, respectively, showing the submicron aerosols to be mostly neutralized by ammonium during the measurement period. The presence of ammonium nitrate, ammonium sulfate, and ammonium chloride can be seen from this neutralized nature of ammonium. The ANR values showed that NH4 + is insufficient to neutralize all the major inorganic anionic species. The remaining fractions of the anionic species could be present in the form of organic compounds like organosulfur and organonitrate (Farmer et al., 2010). However, the chemical nature, concentrations, and formation mechanisms of these compounds during the lockdown period need to be investigated in future studies.

Source Apportionment of OA Using the Unit Mass Resolution Positive Matrix Factorization (UMR-PMF)
The total measured organic species were deconvolved into the characteristic of different organic aerosol components, following the methodologies described in Ulbrich et al. (2009) and Zhang et al. (2005). Using the unit mass resolution positive matrix factorization (UMR-PMF) technique, the number of factors varied from 1 to 5 based on the diagnostic results of Q/Qexp with FPeak = 0 and seed = 0. An optimal solution of four factors was chosen in this study.
The identified factors are broadly categorized as hydrocarbon-like organic aerosol (HOA) and oxygenated organic aerosol (OOA). The mass spectral distribution of HOA and OOA are physically and chemically different, and the different species in the mass spectra act as a tracer of different factors. For example, m/z 57 (mostly C4H9 + ) is a tracer of HOA, m/z 44 (mostly CO2 + ) of LV-OOA, m/z 43 (mostly C2H3O + ) of SV-OOA, and m/z 60 (mostly C2H4O2 + ), and 73 (mostly C3H5O2 + ) of BBOA respectively. The mass spectra of the four-factor solution during the pre-lockdown and lockdown periods are shown in Figs. 3(a-c).
The second factor shows significant peaks at m/z 60 (C2H4O2 + ) and 73 (C3H5O2 + ). These species correspond to the fragmentation of levoglucosan and other sugars considered as a good tracer of biomass burning emissions (Cubison et al., 2011;Mohr et al., 2009;Simoneit et al., 1999;Hu et al., 2013). The peaks at m/z 60 and 73 were insignificant in other factor profiles, implying that the factors were accurately estimated and not splitted. Biomass burning in terms of agricultural waste is not a very common phenomenon in Pune, however, cooking practices using biomasses are common in nearby residential areas. On average, this BBOA factor accounts for 15% of measured organic aerosol mass during the study period that increased from March to May as biomass burning increased in summer in the agricultural lands surrounding Pune.
The third and fourth factors represent the oxygenated organic aerosol (OOA) that generally dominates the OA fraction worldwide (Alfarra et al., 2004;Lanz et al., 2010). These OOA species contain aged organics that are assumed to be formed in the atmosphere from gaseous anthropogenic and biogenic precursors by photochemical and aqueous oxidation. In summer, when the temperature, solar radiation, and ozone concentration are high, affecting the volatilities variability of OOA, these OOA factors can be segregated into volatility-dependent fractions Crippa et al., 2013).
One fraction shows a more prominent peak at m/z 44 and is aged, highly oxidized, and low volatile OA, whereas the other OOA is freshly oxidized, locally produced, and shows a prominent peak at m/z 43 . The ratio of m/z 44 to m/z 43 shows the different degrees of oxidation, and the high concentrations at m/z 44 compared to m/z 43 reflect the aging process and increased degree of oxidation in the atmosphere. The ion signal at m/z 44 generally comes from the thermal decomposition of carboxylic acids and other highly oxygenated organic compounds. These species are strongly processed material and low volatile in nature, so termed as LVOOA. These species contributed the most to the OOA aerosols, contributing nearly 35% to the OA mass during the study period, similar to the inorganic SO4 2aerosols that are primarily generated from combustion sources, marine biological processes, and are long-range transported to the study location.
In contrast to this aged, highly oxidized factor, other factor shows the highest contribution at m/z 29 (CHO + ), 43 (C2H3O + ), and 55 (C3H3O + ). These oxygenated species are comparatively freshly emitted with low oxidation level suggesting a fresh, more locally produced aerosol that was quickly formed from regional precursors due to increased photochemistry during this season. These species could have originated from the nearby urban or industrial areas and started getting oxidized during transportation. The SVOOA contributes nearly 25% to the measured OA mass concentrations during the study period. Since HR-ToF-AMS ionizes molecules with 70 eV electron impact, the generated mass spectra for the aerosol species ensemble reflect the bulk chemical composition. On average, the primary factor contributed nearly 35-40% compared to the secondary factor contributing nearly 60-65% to the OA concentrations. Previous studies have shown that BBOA contributes less to the aerosol mass fraction compared to LVOOA and SVOOA aerosols in summer .
Freshly emitted HOA and SVOOA species contribute comparatively less, inferring the less active nature of primary emitting sources in the lockdown period. The dominating nature of secondary organic aerosol species also reveals the size distribution nature of aerosols, where more oxidized (high f44/43 ratio) aerosol generally remains in the accumulation mode.

High-Resolution Positive Matrix Factorization (HR-PMF) of Organic Aerosol (OA)
In the above section, the organic species were deconvolved into the characteristic of different organic aerosol components using the unit mass resolution (UMR) spectra. Here, we performed PMF analysis on the high-resolution (HR) AMS organic aerosol (OA) mass spectra. The factors are shown in Figs. 4(a-c), and their physical interpretability has been discussed here.
The HR-PMF analysis divides the total OA into four factors similar to UMR-PMF analyses. The elemental ratios (H/C and O/C) of organic aerosols are used, which provides crucial information for understanding the phase state and atmospheric aging of aerosols (Aiken et al., 2008;Zhang et al., 2011). These elemental ratios exhibited large variation during prevailing conditions, e.g., H/C varied from 1.36 to 1.93 with the highest value in the pre-lockdown period and O/C varied from 0.17 to 0.98 with the highest value in the lockdown period. The increased O/C ratio suggested that the bulk OA composition had relatively more oxygenated compounds during lockdown. The organic mass-to-carbon (OM/OC) ratio during different lockdown periods varied from 1.28 to 2.46. The slight decrease in O/C and increase in H/C during pre-lockdown to lockdown periods is likely due to emissions of primary hydrocarbon-like organic aerosols from fuel combustions which were almost non-existent during lockdown periods. The variability in O/C during the study period from other periods could be due to the impact of the aqueous phase processing of aerosols in clouds that increases the O/C of organics.
Using the elemental ratios, the density of organics (ρorg) can be estimated following the methodologies of Kuwata et al. (2012): 12 1 1.6 7 5 4.15 In which ρorg is expressed in g cm -3 , and H/C and O/C are the average values during different measurement periods of this study. We find the predicted ρorg during pre-lockdown, lockdown-I, and lockdown-II periods to be 1.14, 1.28, and 1.35 respectively, to be similar to the reference value i.e., 1.30 g cm -3 Canagaratna et al., 2015).
The different factors are characterized by their different elemental ratios, like hydrocarbonloke OA (HOA) showed high H/C = 1.93 and low O/C = 0.17 value, whereas oxygenated OA (LVOOA and SVOOA) showed high O/C and low H/C ratio. A number of previous studies have also found similar properties in the elemental ratios of OA factors. The different volatility profiles of these factors follow the pattern of the worldwide datasets reported in Jimenez et al. (2003) and . Primary organic aerosol species like HOA shows higher H/C and lower O/C values, and a higher N/C ratio is observed when higher interaction of organics with NH3 occurs. The variability in H/C, O/C, and OM/OC ratios for HOA, LV-OOA, SV-OOA, and BBOA during the study period are shown in Fig. 4, and all these elemental ratios showed distinct characteristics during different periods.

Back Trajectory and Cluster Analysis
The influence of air mass origin on the species concentrations was investigated by tracing the path of the air masses reaching the study location. During the study period, 24-hour backward trajectories were obtained using the HYSPLIT model (Stohl, 1996;Draxler and Rolph, 2003) at an altitude of 500 m above ground level (AGL) every hour for each day. This height was chosen which minimizes the effect of surface friction and represents winds in the lower boundary layer (Begum et al., 2005;Gao et al., 1993). The estimated back trajectories were then grouped into clusters using the trajectory endpoints (lat, long), and each cluster represents different geographical origins (Dumka et al., 2013).
These clusters determine and demarcate the relative importance of regional and distant source regions. Various receptor models have been developed, combining meteorology with the measured chemical compositions to determine the preferred pathways that give rise to the observed high particle concentrations at the sampling site (Hopke, 2003). The primary criterion of trajectory clustering is to minimize the variability among the trajectories within a cluster and maximize the variability between the clusters (Vinoj et al., 2010). The optimal number of clusters was assessed from the variation of total spatial variance (TSV), and when TSV increased rapidly, that cluster number has been used in the analysis (Draxler, 1992).
In our analysis, 3 clusters were identified as the final simulated cluster trajectories. The back trajectory and cluster analysis of the measurement period are shown in Fig. 5. Based on the trajectory analysis, the most common wind directions contributing to PM concentrations were from the west and northwest directions, carrying air masses primarily from land and ocean. Cluster analysis revealed three major clusters representing the mean trajectory pathways. In the pre-lockdown period, the majority of air parcels were coming from two regions, the north, and northwest directions. While during the lockdown period, air masses were predominantly from the northwest. These air mass trajectories are consistent with the previous studies reported by Ali et al. (2012) over the same study location. Cluster 1 was the major cluster representing 57% of air masses that might have traveled through Mumbai to reach Pune, which carried more pollutants to Pune. Cluster 2 represented 27.8% of air masses from the northwest direction, primarily originating from the coastal regions, the Arabian Sea, and to some extent from Gujarat and the Thar Desert. Cluster 3 represented 14.4% of regional air masses near Pune, depicting major source regions near neighboring areas. The mass concentrations of the different NR-PM1 components and the air mass direction combinations might give an insight into the geographical origins of chemical species. Fig. 6 shows the PSCF results in the form of maps of the measurement location on which the weighted PSCF (WPSCF) values ranging from 0 to 1 are displayed in terms of a color scale. For PSCF analysis, the domain covered by the trajectories was divided into grids with 0.2° × 0.2° resolutions. These cells are indicative of areas of potential contributions for PM pollutants, and grid cells with higher PSCF values are more likely to contribute to the constituent pollutants. The spatial distribution of WPSCF values shows that WPSCF > 0.7 was over the north and northwest parts of Pune, indicating these locations to be the major source areas of PM. From the distribution of WPSCF values, it can be seen that the sources most likely to affect Pune are located in the (1) industrialized and populated city of Mumbai, indicating shorter transport pathways (2) the coastal and ocean region of the Arabian Sea (3) industrial regions of Rajasthan and Gujarat region.

Weighted Potential Source Contribution Function (WPSCF)
The PSCF plots highlight the dominant oceanic origin of the constituents of NR-PM1, and some chemical species are transported or formed above the ocean rather than directly emitted by sea spray (Ovadnevaite et al., 2012). Under the influence of air masses from the Arabian Sea, emissions of dimethyl sulfide (DMS) and organosulfur compounds by microorganisms could contribute to the observed Org and SO4 concentrations (Luria et al., 1996). Previous studies have also shown the contribution of sulfur and Org matter from the marine air in PM1 chemical compositions (Fomba et al., 2014). These species may also be released by maritime traffic across coastal regions, shipping emissions, and anthropogenic activities in distant continental parts, whose emissions may be carried toward the ocean and brought back to Pune by western winds. During the transport of Org over the ocean surface, SOA formation also enhances, which under specific atmospheric conditions can be initiated by chlorine atoms (Hallquist et al., 2009;Hossaini et al., 2016).

Fig. 5.
The 24 h backward trajectories of air masses at 500 m above ground level arriving at the measurement location. All the near trajectories were merged to a mean trajectory to represent the entire groups by cluster analysis.
Sources located in Mumbai may influence PM by direct emissions through industrial activities, biomass and waste burning, dense populations, and high traffic emissions under prevailing meteorological conditions. Recirculating narrow and short trajectories from these adjoining regions result in high source contributions at the receptor site. PSCF maps showed a high frequency of potential source location (WPSCF > 0.7) over industrialized areas of Maharashtra. Air masses from the southern part of the measurement location impact less to the pollutant concentrations.
The northern source regions are responsible for bringing significant dust in the summer season from regions like Rajasthan. Sources in the nearby regions also appear to be important sources for PM especially due to regional sources like biomass burning, waste burning, brick kilns, construction activities, and thermal power plants. Thus, from the assessment of clusters and spatial variation of WPSCF, the potential source sectors can be estimated, helping in policymaking to control emissions.

SUMMARY AND CONCLUSIONS
In this study, we present the comprehensive temporal variation and physicochemical characterization of Org and inorganic constituents of NR-PM1 measured with an HR-ToF-AMS instrument during the pre-lockdown and lockdown periods in India. This kind of real-time high temporal resolution measurement of submicron particles is particularly relevant in India, where air pollution affects the population to a great extent. This assessment would characterize the submicron aerosol (NR-PM1) possibly impacted by different source regions during the COVID-19 lockdown restrictions. The following conclusions can be made from this study: The composition data of NR-PM1 over Pune showed that Org dominated the mass fractions of submicron aerosols with nearly 60% contribution. The mass concentration was followed by sulfate, ammonium, and nitrate. Non-refractory chloride contributed very little to the PM1 mass fractions.
During the lockdown period, the Org and nitrate concentrations decreased by 60% and 53% from the pre-lockdown period. The significant reduction in the primary emission from traffic, coal combustion, and regional biomass burning might have impacted the reduction in concentrations.
The size-resolved submicron aerosol measurements showed the peaks of Org and inorganic constituents at an aerodynamic diameter of 550 nm, referring to the higher presence of aged, oxidized constituents and less freshly emitted aerosols in Pune. These peaks could also be due to the regional secondary aerosol formation from the gaseous precursors and strong oxidation in the urban environment.
Unit mass resolution positive matrix factorization (UMR-PMF) analysis was applied to the OA mass spectra and four OA factors were identified, hydrogenated organic aerosol (HOA), biomass burning organic aerosol (BBOA), low volatile oxygenated organic aerosol (LVOOA), and semivolatile oxygenated organic aerosol (SVOOA).
The oxygenated organic aerosol (OOA) contributed more (60%) to the OA mass compared to hydrogenated organic aerosol (HOA) (40%). Among the OA factors, LVOOA contributed the most, similar to the inorganic SO4 2species. The HOA factor corresponded primarily to vehicular emission, whereas the BBOA factor corresponds to biomass burning. The two OOA species LVOOA and SVOOA correspond to the secondary formation mainly from the photochemical and aqueous phase oxidation of the precursors.
The diurnal variations of the constituents showed the dominant role of vehicular emission and boundary layer dynamics affecting the ambient concentrations. High emissions from vehicles mostly contributed to the nitrate concentration and led to the two concentration peaks in nitrate during office hours in the morning and evening. Whereas during the lockdown period, the concentrations showed almost constant values at all hours with no peaks or perturbation in the diurnal signal.
Back trajectories, cluster analysis, and PSCF analysis identified Mumbai and the Arabian Sea as potential source regions of air masses and responsible for the regional aerosol loading during the observation period.
High-Resolution Positive Matrix factorization (HR-PMF) analysis divides the total OA into four factors, similar to UMR analyses. The elemental ratios (H/C and O/C) of organic aerosols are used, which provides crucial information to understand the phase state and ageing of aerosols. These elemental ratios exhibited large variation during prevailing conditions, e.g., H/C varied between 1.36 to 1.93 with the highest value in the pre-lockdown period and O/C varied between 0.17 to 0.98 with the highest value in the lockdown period. The increased O/C ratio suggested that the bulk OA composition had relatively more oxygenated compounds during lockdown.
Using these ratios, we estimated the ρorg during pre-lockdown, lockdown-I, and lockdown-II periods to be 1.14, 1.28, and 1.35 g cm -3 respectively. The stoichiometric neutralization analysis showed the aerosol neutralization ratio (ANR) to be nearly 0.8, implying the higher presence of ammonium salts like ammonium sulfate and ammonium nitrate over Pune.

Data Availability
All the data used to prepare the manuscript are archived on the Prithvi (IITM) super-computer and can be provided upon request to the corresponding authors.