how the atmosphere over eastern himalaya, india is polluted ...darjeeling) and ends at darjeeling....

18
Aerosol and Air Quality Research, 17: 2206–2223, 2017 Copyright © Taiwan Association for Aerosol Research ISSN: 1680-8584 print / 2071-1409 online doi: 10.4209/aaqr.2017.01.0048 How the Atmosphere over Eastern Himalaya, India is Polluted with Carbonyl Compounds? Temporal Variability and Identification of Sources Chirantan Sarkar 1 , Abhijit Chatterjee 1,2* , Dipanjali Majumdar 3 , Arindam Roy 1 , Anjali Srivastava 5 , Sanjay K. Ghosh 2,4 , Sibaji Raha 2,4 1 Environmental Science Section, Bose Institute, Kolkata 700054, India 2 National Facility on Astroparticle Physics and Space Science, Bose Institute, Darjeeling 734101, India 3 National Environmental Engineering Research Institute, Kolkata 700107, India 4 Center for Astroparticle Physics and Space Science, Kolkata 700091, India 5 National Environmental Engineering Research Institute, Nagpur 440020, India ABSTRACT A study was conducted on atmospheric carbonyl compounds for the first time over a Himalayan atmosphere in India. Samples were collected from a high altitude hill station, Darjeeling (27.01°N, 88.15°E, 2200 masl) during June 2012 to May 2013. Temporal variation, meteorological influence, source apportionment and ozone formation potential etc were studied for acetaldehyde, formaldehyde, acetone, butanaldehyde, propanaldehyde, benzaldehyde, crotonaldehyde, valeraldehyde, isovaleraldehyde, hexanaldehyde, p-tolualdehyde and o-tolualdehyde. High concentration of Acetone (81.6 ± 63.5 μg m –3 ) was observed which could be due to the higher photochemical production from its precursor volatile organic compounds emitted from tea plants and tea processing units. The concentration of acetaldehyde (20.7 ± 47.6 μg m –3 ) and formaldehyde (11.6 ± 16.3 μg m –3 ) were found to be comparable with metro cities of India and world. The average annual concentration of total carbonyl compounds was 174.2 ± 184.6 μg m –3 with maximum during postmonsoon (456.9 ± 199.5 μg m –3 ) and minimum during winter season (72.2 ± 42.9 μg m –3 ). Meteorological parameters like temperature and surface reaching solar radiative flux played the major roles for the seasonal variation of the carbonyl concentration over the hill station. The average ratio of formaldehyde to acetaldehyde over Darjeeling was found to be 1.64 ± 1.43 well representing a typical urban atmosphere at this part of Himalaya. Positive matrix factorization model showed that the biogenic emissions from tea plants and vehicular emissions were the major sources of carbonyl compounds over the hill station. Keywords: Himalaya; Carbonyl compound; Principal component analysis; Darjeeling. INTRODUCTION Atmospheric carbonyl compounds (CCs) are oxygenated volatile organic compounds, receiving worldwide attention (Andreini et al., 2000; Bakeas et al., 2003; Cavalcante et al., 2006; Pal et al., 2008; Uchiyama et al., 2012; Ho et al., 2013) since the last few decades. They are ubiquitous and are responsible for the degradation of poor air quality through the formation of tropospheric ozone and peroxy acetyle nitrate which in turn forms photochemical smog. They also adversely affect human health including carcinogenesis upon acute or chronic exposure (Jacob, 2000). Exposure to atmospheric aldehydes and ketones cause * Corresponding author. Tel.: +91-9051585800; Fax: +91-33-2355 8386 E-mail address: [email protected] serious health effects. The World Health Organisation (WHO) classified some CCs as toxic air pollutants (TAPs), such as formaldehyde and acetaldehyde for their adverse health effects and carcinogenic consequences on human beings (WHO, 1987; CEPA, 1993). CCs have drastic environmental effect too as they are responsible for the formation of atmospheric free radicals through photo oxidation during day time (Jacob, 2000) contributing a major part of the photochemical smog (Possanzini et al., 1996; Atkinson, 2000; Solberg et al., 2001; Possanzini et al., 2003). These compounds initiate the peroxy-acetyl nitrate or PAN formation which is responsible for photochemical smog. The presence of CCs in the atmosphere is thus of great concern in respect of their adverse effects on human health and environment. Many studies have been performed to identify their sources and it was found that anthropogenic emissions (Kean et al., 2001; Kim et al., 2008) and biogenic emissions (Seco et al., 2007; Guo et al., 2009) are the major sources

Upload: others

Post on 07-May-2021

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Aerosol and Air Quality Research, 17: 2206–2223, 2017 Copyright © Taiwan Association for Aerosol Research ISSN: 1680-8584 print / 2071-1409 online doi: 10.4209/aaqr.2017.01.0048 

How the Atmosphere over Eastern Himalaya, India is Polluted with Carbonyl Compounds? Temporal Variability and Identification of Sources Chirantan Sarkar1, Abhijit Chatterjee1,2*, Dipanjali Majumdar3, Arindam Roy1, Anjali Srivastava5, Sanjay K. Ghosh2,4, Sibaji Raha2,4 1 Environmental Science Section, Bose Institute, Kolkata 700054, India 2 National Facility on Astroparticle Physics and Space Science, Bose Institute, Darjeeling 734101, India 3 National Environmental Engineering Research Institute, Kolkata 700107, India 4 Center for Astroparticle Physics and Space Science, Kolkata 700091, India 5 National Environmental Engineering Research Institute, Nagpur 440020, India ABSTRACT

A study was conducted on atmospheric carbonyl compounds for the first time over a Himalayan atmosphere in India. Samples were collected from a high altitude hill station, Darjeeling (27.01°N, 88.15°E, 2200 masl) during June 2012 to May 2013. Temporal variation, meteorological influence, source apportionment and ozone formation potential etc were studied for acetaldehyde, formaldehyde, acetone, butanaldehyde, propanaldehyde, benzaldehyde, crotonaldehyde, valeraldehyde, isovaleraldehyde, hexanaldehyde, p-tolualdehyde and o-tolualdehyde. High concentration of Acetone (81.6 ± 63.5 µg m–3) was observed which could be due to the higher photochemical production from its precursor volatile organic compounds emitted from tea plants and tea processing units. The concentration of acetaldehyde (20.7 ± 47.6 µg m–3) and formaldehyde (11.6 ± 16.3 µg m–3) were found to be comparable with metro cities of India and world. The average annual concentration of total carbonyl compounds was 174.2 ± 184.6 µg m–3 with maximum during postmonsoon (456.9 ± 199.5 µg m–3) and minimum during winter season (72.2 ± 42.9 µg m–3). Meteorological parameters like temperature and surface reaching solar radiative flux played the major roles for the seasonal variation of the carbonyl concentration over the hill station. The average ratio of formaldehyde to acetaldehyde over Darjeeling was found to be 1.64 ± 1.43 well representing a typical urban atmosphere at this part of Himalaya. Positive matrix factorization model showed that the biogenic emissions from tea plants and vehicular emissions were the major sources of carbonyl compounds over the hill station.

Keywords: Himalaya; Carbonyl compound; Principal component analysis; Darjeeling. INTRODUCTION

Atmospheric carbonyl compounds (CCs) are oxygenated volatile organic compounds, receiving worldwide attention (Andreini et al., 2000; Bakeas et al., 2003; Cavalcante et al., 2006; Pal et al., 2008; Uchiyama et al., 2012; Ho et al., 2013) since the last few decades. They are ubiquitous and are responsible for the degradation of poor air quality through the formation of tropospheric ozone and peroxy acetyle nitrate which in turn forms photochemical smog. They also adversely affect human health including carcinogenesis upon acute or chronic exposure (Jacob, 2000).

Exposure to atmospheric aldehydes and ketones cause * Corresponding author.

Tel.: +91-9051585800; Fax: +91-33-2355 8386 E-mail address: [email protected]

serious health effects. The World Health Organisation (WHO) classified some CCs as toxic air pollutants (TAPs), such as formaldehyde and acetaldehyde for their adverse health effects and carcinogenic consequences on human beings (WHO, 1987; CEPA, 1993). CCs have drastic environmental effect too as they are responsible for the formation of atmospheric free radicals through photo oxidation during day time (Jacob, 2000) contributing a major part of the photochemical smog (Possanzini et al., 1996; Atkinson, 2000; Solberg et al., 2001; Possanzini et al., 2003). These compounds initiate the peroxy-acetyl nitrate or PAN formation which is responsible for photochemical smog. The presence of CCs in the atmosphere is thus of great concern in respect of their adverse effects on human health and environment.

Many studies have been performed to identify their sources and it was found that anthropogenic emissions (Kean et al., 2001; Kim et al., 2008) and biogenic emissions (Seco et al., 2007; Guo et al., 2009) are the major sources

Page 2: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2207

of carbonyl compounds. Anthropogenic sources like vehicular emissions, industrial emission, solvent evaporations (Friedrich and Obermeier, 1999), biomass burning (Crutzen et al., 1985; Tanner et al., 1996) etc. produce significant amount of CCs in the atmosphere. CCs may also form from the directly emitted volatile organic compounds by secondary photochemical transformation (Possanzini et al., 1996; Lary and Shallcross, 2000; Christensen et al., 2000; Bakeas et al., 2003; Obermeyer et al., 2009).

Most of the CCs have atmospheric lifetime of a few hours with respect to photolysis and OH photo-oxidation reactions (Seinfeld and Pandis, 1998; Macdonald et al., 2001) during day time except acetone which has much higher life time of several days (Atkinson 2000). Night time reactions with ozone and nitrate radical act as the major sink of CCs (Slemr and Junkermann, 1996; Finlayson-Pitts and Pitts, 2000). CCs are also lost from the atmosphere through dry and wet deposition (Finlayson-Pitts and Pitts, 2000).

Several studies were performed on CCs across the world (Moussa et al., 2006; Ohura et al., 2006; Kim et al., 2008; Duan et al., 2012), but very few studies have been made so far in India. Studies were made over metropolitan cities like Kolkata (Dutta et al., 2009; Majumdar et al., 2014), Mumbai (Majumdar and Srivastava, 2012) and Agra (Khare, 1997) in India to determine the ambient concentration of few CCs and their emissions from waste dumping sites. As far as Himalaya is concerned, a case study was performed on atmospheric CCs in Bhutan (Wangchuk et al., 2015). Earlier we have reported high ambient level of black carbon aerosols (Sarkar et al., 2015), volatile organic compounds (Sarkar et al., 2014), submicron aerosols (Adak et al., 2014), anthropogenic aerosols (Chatterjee et al., 2010, 2012; Roy et al., 2016) over Darjeeling at eastern Himalaya. This prompted us to make a year-long study on other major part of carbonaceous pollutants like carbonyl compounds over this part of Himalaya where the loading of atmospheric pollutants are of major concern due to increased anthropogenic activities related to tourists and rapid pace of urbanization and socio-economic development.

The present study is focused on the investigation of the ambient concentrations and temporal variations of carbonyl compounds, identification of the sources of the carbonyl compounds and the impact of carbonyl compounds on the tropospheric ozone formation over Darjeeling, a high altitude hill station at eastern part of Himalaya in India.

DESCRIPTION OF THE STUDY AREA AND SYNOPTIC METEOROLOGY

The study was carried out at Darjeeling (27°01′N, 88°15′E, 2200 m asl), a well known tourist destination and a hill station over eastern Himalaya in India with the population of ~100,000 people. The geographical position of the sampling site and adjacent regions are shown in a map in Fig. 1. The densely populated Indo-Gangetic Plain (IGP) is located ~50 km south of Darjeeling. The main town centre is located ~200 m below the sampling site and ~250 m below the ridge top. The area in and around the main town centre is populated mainly by several hotels and

resorts, offices, schools and colleges, bus and car stands, shopping malls, and hence is the main source of anthropogenic CCs. The road with high vehicular density is the National Highway (NH55) which starts from a plain land city called Siliguri (~2000 m below and south of Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs by both coal and diesel along the national highway. The Darjeeling railway station and the railway workshop are situated in the main town centre. The measurement was carried out at a height of ~10 m above ground in the campus of Bose Institute. The road adjacent to the campus starts from the main town centre and ends at a tourist spot which is ~300 m west of our campus with minimum traffic density. There are several tea gardens and tea processing units/plants in and around the sampling site and a total of ~90 tea gardens covering 17500 hectares of land present in Darjeeling. About 9 millions kilogram of tea is produced per year by these tea gardens engaging more than 50% of the people of Darjeeling.

Table 1 summarizes the seasonal average of surface meteorological parameters during the entire sampling days. The study period was segregated into four seasons; winter (December–February), premonsoon (March–May), monsoon (June–September) and postmonsoon (October–November). The annual average temperature was 11.3 ± 4.4°C with the maximum during monsoon (15.9°C) followed by premonsoon (12.6°C), postmonsoon (11.2°C) and winter (5.3°C). The average temperature during monsoon was ~3 times higher than winter. Darjeeling had a typical humid atmosphere with an annual average humidity of 86.1 ± 5.0%. Wind speed was highest in premonsoon (1.3 ± 0.1 m s–1) while the average wind speed during postmonsoon, monsoon and winter was similar (~0.7 m s–1). The surface reaching solar radiation flux was also highest during premonsoon season (280 ± 104.9 w m–2) and clear sky conditions were observed mostly during this season (10 out of 12 sampling Days). The minimum solar radiation flux was recorded during winter season (177.1 ± 91.6 w m–2). The total rainfall for the entire sampling days was recorded as 493.7 mm. However, ~67% rain occurred during monsoon (329.8 mm) season only.

METHODOLOGY Sampling Protocol and Analytical Procedure

The USEPA TO-11A (USEPA, 1996, 1999) method was followed for sampling and analysis of carbonyls due to its simplicity and good reproducibility. Analysis was done in the laboratory of CSIR (Council of Scientific and Industrial Research, Govt of India)-National Environmental Engineering Research Institute, Kolkata Zonal Laboratory.

Air samples were collected from June 2012 to May 2013 on day (7:00 am–7:00 pm) and night (7:00 pm–7:00 am) basis at every 8 day interval. One weekly sampling was not done during February and March 2013 respectively due to some technical and logistic issue. The sampling and analysis of CCs were done following the method used by Sirju and Shepson (1995), Lee et al. (2001) and Pang et al. (2001).

Page 3: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2208

Fig. 1. The geographical location of the study area over Darjeeling.

Table 1. Seasonal variation of the meteorological parameters during the study period.

Monsoon Postmonsoon Winter Premonsoon Annual Temp (°C) 15.9 ± 1.1 11.2 ± 2.6 5.3 ± 1.8 12.6 ± 1.9 11.3 ± 4.4 RH (%) 93.3 ± 6.2 84.5 ± 9.1 82.0 ± 12.8 84.7 ± 12.9 86.1 ± 5.0 SR (W m-2) 201.0 ± 95.0 224.0 ± 97.2 177.1 ± 91.6 280.0 ± 104.9 220.6 ± 44.0 WS (m s–1) 0.8 ± 0.1 0.7 ± 0.1 0.8 ± 0.1 1.3 ± 0.1 0.9 ± 0.2 Total Rain (mm) 329.8 0.1 37.7 126.1 493.7

Page 4: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2209

Briefly, air samples were collected by drawing air with a sampling pump (Envirotech, India) through DNPH-coated cartridges with build in ozone scrubber (BPE-DNPH cartridge sigma- Aldrich 54728-U SUPELCO). The sample collection time was 12 h; the flow rate was ~100 mL min–1 (sample volume ~7.2 L). The cartridges were sealed and well wrapped in aluminium foil and stored in refrigerator at 4°C immediately after sampling untill further analysis. A total of eighty six samples were collected during the entire sampling event. One laboratory blank and one field blank were included during sampling and analysis of carbonyls compounds.

2,4-Dinitrophenylhydarazone derivatives were eluted using formaldehyde free HPLC grade aetronitrile. 2–3 mL of acetonirile was used to elude the derivatives and stored in the vial in a refrigerator for analysis. The 2,4-dinitrophenylhydrazone derivatives were analyzed using a Dionex HPLC. A five point calibration curve was used to calibrate the HPLC using a standard calibration mixture of 15 carbonyl species procured from Supelco, USA. 10 µL sample was injected into a reversed-phase C18 (RP-HPLC) column (Supelco make, C18 columns and guard column) and the derivatives were separated using gradient elution programme. A UV detector operating at 360 nm was used to quantify the derivatives. A mixture of acetonitrile and tetrahydrofuran (10% THF in water) were used as a solvent for the gradient elution programme with the flow rate of 0.3 mL min–1. Quality Assurance and Quality Control

To ensure the quality of air sampling method we have performed the breakthrough experiment and the cartridge collection efficiency were also determined. For this experiment two DNPH cartridges were connected in series and air sampling was done with the same flow rate and same duration mentioned earlier. Formaldehyde was chosen as the reference carbonyl for this experiment. Breakthrough was assumed to be occurred when the concentration of formaldehyde in the second cartridge of the series exceeded 10% of the first cartridge concentration. The breakthrough volume was 10 liter of air and ~99% of the CCs were found

to be collected in the first cartridge when the sampling volume was under the breakthrough volume. The US-EPA method of detection limit development (https://www3.epa. gov/ttn/amtic/files/2009conference/swift.pdf) was followed for finding the method detection limit (MDL). The MDL was calculated as 3 times of the standard deviation of seven replicate measurements of known standard solution of each carbonyl species. The MDL values of each CC are given in Table 2. The precision of determination was calculated similarly by repeated analysis of a known standard solution and the precision was found to vary within ± 5% for each CC concentration.

RESULTS AND DISCUSSION General Characteristics of the Carbonyl Compounds and Comparison with Other Studies

The annual average concentrations, standard deviations, maximum and minimum of individual and total CCs are given in Table 2. The average annual concentration of total CCs was 174.2 ± 184.6 µg m–3. Acetone had the maximum abundance with the annual average concentration of 81.6 ± 63.5 µg m–3 ranging between as low as 12.7 µg m–3 and as high as 391.1 µg m–3 contributing 47% to the total annual CCs concentration. Acetaldehyde had the second highest annual average concentration of 20.7 ± 47.6 µg m–3 ranging from 1.1 to 315.1 µg m–3 followed by crotonaldehyde, butanaldehyde, formaldehyde and others. The percentage contributions of all the CCs are shown in Fig. 2. Acetaldehyde contributed 12% whereas each of formaldehyde, butanaldehyde and crotonaldehyde contributed 7% to the total annual CCs concentrations.

The most important observation in this study is the very high concentration of acetone which was found to be 4 and 6 times higher than formaldehyde and acetaldehyde respectively. This high concentration of acetone could be explained by the contribution of variety of sources like vehicular emission (Dong et al., 2014), biomass burning (Lewis et al., 2013), solvent evaporation (Feng et al., 2005; Pang and Lee, 2010), biogenic emissions (Singh et al., 1994; Janson and Serves, 2001) as well as secondary photochemical

Table 2. Minimum, maximum and average concentrations of measured carbonyl compounds. Concentrations are given in µg m–3.

Annual average SD Max Min MDL Formaldehyde 11.6 16.3 91.0 2.1 0.221 Acetaldehyde 20.7 47.6 315.1 1.1 0.139 Acetone 81.6 63.5 391.1 12.7 0.080 Propanaldehyde 2.9 5.4 24.3 BDL 0.133 Butanaldehyde 12.2 27.1 123.9 BDL 0.083 Crotonaldehyde 13.1 29.3 151.2 0.8 0.126 Benzaldehyde 9.3 26.6 152.3 BDL 0.134 Isovaleraldehyde 5.7 12.4 63.4 BDL 0.187 Valeraldehyde 8.0 15.0 85.0 BDL 0.219 o-Tolualdehyde 1.0 3.6 25.9 BDL 0.074 p-Tolualdehyde 2.7 5.3 23.4 BDL 0.247 Hexanaldehyde 5.5 11.6 62.6 BDL 0.214 Total 174.2 184.6 795.5 22.4

Page 5: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2210

Fig. 2. Contributions of individual carbonyls to the total carbonyl compounds concentration.

productions (Singh and Hanst, 1981; Chatfield et al., 1987; Singh et al., 1994, 1995; Reissell et al., 1999; Orlando et al., 2000) in the atmosphere of Darjeeling. Moreover acetone has longer lifetime (~15 Days) compared to other carbonyl compounds (few hours). Thus it could also be transported from long distant source regions like densely populated Indo-Gangetic Plane (IGP), sustain in the atmosphere for a longer period of time and enhance the ambient level of acetone concentration over the hill station. However an important factor behind the high level of acetone over Darjeeling could be the biogenic emission of its precursors followed by their atmospheric degradation.

Darjeeling is world famous for the production of flavoured tea. Phytofluene is a carotenoid pigment mainly found in the tea leaves (Premachandra et al., 1977; Ho et al., 2015) undergoes epoxidation and photo-oxidation and produce a volatile organic compound called geranyl acetone ((5E)-6,10-dimethylundeca-5,9-dien-2-one) which is an important aroma of black, green and oolong tea (Ho et al., 2015). Geranyl acetone having very short life time of less than 1 hour (Fruekilde et al., 1998) could further get degraded and ozonolysed to form acetone. The reaction pathway of the production of acetone from phytofluene is given in Fig. 3.

The population of tea leaves depends on the pruning and skiffing of tea plants. Pruning and skiffing helps the tea plants to provide stimulus for vegetative growth, to correct past defects in bush architecture, to reduce the incidence of pests and diseases and to regulate the crop. Winter is the pruning time of tea plants over Darjeeling resulting the decrease in the population of tea leaves. New leaves began to grow during the premonsoon season whereas heavy rainfall during monsoon helps the tea plants to grow quickly and their growth continues until the postmonsoon season before the next pruning time. Acetone concentration over Darjeeling follows this seasonal growth cycle of the

tea plants showing maximum concentration in postmonsoon followed by monsoon, premonsoon and winter.

The ambient level of formaldehyde, acetaldehyde and acetone over Darjeeling have been compared with other studies and shown in Table 3. These three carbonyl compounds are the common amongst the carbonyls measured by most of the studies made earlier. The highest level of acetone over Darjeeling than other stations across the globe could be easily explained in terms of the biogenic emissions as discussed above. It is observed from Table 3 that the concentrations of formaldehyde and acetaldehyde were much higher than the study made at Bhutan at eastern Himalaya. Formaldehyde shows ~5 times and acetaldehyde shows ~20 times higher loading in the atmosphere of Darjeeling compared to Bhutan. The concentration of formaldehyde over Darjeeling was 2–5 times higher than some of the urban cities like Shimizu (Japan), Beirut (Lebanon), Niteroi (Brazil) and ~10 times higher than that in Agra (India). But, formaldehyde over Darjeeling shows much lower concentration than that in Kolkata, a metro city located at Indo-Gangetic Plain in eastern India. The ambient level of acetaldehyde over Darjeeling was found to be higher than most of the cities in India and abroad. Earlier, we observed higher concentrations of carbonaceous aerosols and VOCs (BTEX: Benzene, Toluene, Ethylbenzene, Xylene) over Darjeeling from various anthropogenic activities and their concentrations were comparable and higher than some of the metro-cities in India (Sarkar et al., 2014, 2015). The sources of carbonaceous aerosols and BTEX are similar to that of acetaldehyde over Darjeeling. Seasonal Variation of the Carbonyl Compounds

The seasonal variation of total CCs is shown in Fig. 4 and the seasonal average concentrations of individual CCs are given in Table 4. Maximum concentration for total CCs

Page 6: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2211

Fig. 3. Formation pathway of acetone from phytofluene via geranylacetone.

was observed during the postmonsoon season (456.9 ± 199.5 µg m–3) followed by monsoon (167.5 ± 142.2 µg m–3), winter (72.2 ± 42.9 µg m–3) and premonsoon (73.5 ± 38.2 µg m–3). The CCs concentrations during premonsoon and winter were found to be almost equal. The average CCs concentration during postmonsoon was ~2.5 times higher than monsoon and ~6 times higher than winter and premonsoon. The seasonal variations of most of the CCs were similar to that of total CCs except crotonaldehyde. Crotonaldehyde had the maximum concentration during monsoon season (23.6 ± 42.6 µg m–3).

Due to the massive tourist influx, Darjeeling experiences huge emissions of carbonaceous pollutants from various anthropogenic activities like vehicular emission and fossil fuel burning during premonsoon and postmonsoon which are the peak tourist seasons in this high altitude station.

Black carbon aerosols during premonsoon and VOCs during postmonsoon were also observed to be higher in our previous studies (Sarkar et al., 2014, 2015). Anthropogenic activities like vehicular emission and biomass burning have been reported as the significant sources of CCs in various studies all over the world (Possanzini et al., 1996; Kean et al., 2001; Ho et al., 2002a, b; Wang et al., 2010; Chen et al., 2014). In addition to that significant emission flux of many CCs has been observed from the waste dumping sites (Majumdar et al., 2014). The human, animal and other domestic wastes are treated unscientifically over this hill station. Thus one would expect high concentration of CCs during both these seasons. Instead, we observed very low concentration of CCs during premonsoon compared

to postmonsoon season (Table 4). The phenomena can be explained by the effect of

meteorological influence on the concentration of the carbonyl species. To better understand the effect of meteorological parameters we have made correlation analysis between CCs and micro-meteorological parameters over different seasons and shown in Table 5. It is observed that ambient temperature and surface reaching solar radiation could play different roles on the CCs concentration during different seasons. The incoming solar radiation flux during premonsoon season (280.0 ± 104.9 W m–2) was much higher than any other seasons. Average temperature during premonsoon was also higher than postmonsoon and winter season (Table 1). The higher temperature and intense solar radiation favour the photolysis of the carbonyl species (Atkinson, 1990; Anderson et al., 1996; Christensen et al., 2000). On the other hand ambient temperature and solar radiation during postmonsoon were comparatively lower and thus could help in building up the CCs over the hill station.

The second highest concentration of total CCs was found during monsoon (Table 4). The ambient temperature might play the major role in increasing the concentration of CCs in this season. The ambient temperature during monsoon was significantly higher than other seasons (Table 1). Thus fuel and solvent evaporation could be responsible for enhancing the CCs concentration over the study area as reported by the other studies too (Lee et al., 2001; Chen et al., 2014). Crotonaldehyde had the second highest concentration (23.6 ± 42.6 µg m–3) after acetone during monsoon. Crotonaldehyde is used in industrial sectors like industrial organic chemicals

Page 7: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2212

(Salgado, 2008; Grosjean, 1990), solvent in production of dyes, pesticides, pharmaceuticals, rubber antioxidants etc. Hence higher temperature during monsoon season could help to evaporate crotonaldehyde from those small scale industries situated at the foothill regions and could be transported to the study area as it has relatively higher atmospheric lifetime (half life ~11 hr).

The minimum CCs concentration was observed during the winter season. This observation is quite different from other studies performed over the metropolitan cities or rural areas in plane land regions in different parts of the world. In those studies the seasonal variation was mainly explained by mixing layer/boundary layer dynamics. During the hotter months, the increase in the height of the boundary layer decreases the surface level concentrations. But during the colder months the boundary layer remains much closer to the surface, trapping the pollutants near the surface and the concentration of the pollutants gets increased. But the scenario of high altitude stations like Darjeeling is quite opposite from the plane land regions. During the winter season, the boundary layer remains well below Darjeeling (2.2 km) and the atmosphere of Darjeeling remains somewhat pristine as the pollutants from the low land region are not able to reach to the hill top. During this season anthropogenic activities related to tourists also get decreased significantly and hence the local generation of CCs. Due to very low temperature, the evaporation of solvents also become insignificant. Day/Night Ratios of Carbonyl Compounds

The day/night concentration ratio for total and most of the CCs was > 1 in all the seasons as shown in Fig. 5. The higher anthropogenic emissions and the photolysis of VOCs could significantly enhance the day time concentrations of CCs. The NO3 radicals could also act as a sink and decrease the concentration of CCs during night time (Atkinson, 2000). But acetone and butanaldehyde were found to have slightly higher night time concentrations during postmonsoon and winter. One of the major sources of acetone and butanaldehyde is biomass/wood burning (Zhang and Smith, 1999). The ambient temperature during night remains much lower during postmonsoon and winter compared to other seasons when biomass burning activities get increased by the local people (tourists, local habitants as well as villagers). This night time massive biomass/wood burning could act as additional source to increase the night time concentration of these CCs during postmonsoon and winter. Effect of Meteorology on Temporal Variation of the Carbonyl Compounds

We have used an automated weather station to continuously monitor different meteorological parameters throughout the whole sampling period. The weather station is situated at a roof top at 15 m above the ground level and 5 m away from sampling area. The weather station was equipped with the different sensors like temperature and relative humidity sensor (Model No. HMP45C-L, Campbell Scientific, USA) wind speed and win direction measurement sensor (Anemometer: Model No. 05103, Young), total

Tab

le 3

. Com

pari

son

of f

orm

alde

hyde

, ace

tald

ehyd

e an

d ac

eton

e co

ncen

trat

ions

ove

r D

arje

elin

g w

ith

othe

r st

atio

ns.

Stu

dy a

rea

Dar

jeel

ing,

In

dia

Kun

glun

g,

Bhu

tan

Agr

a,

Indi

a K

olka

ta,

Indi

a S

him

izu,

Ja

pan

Bei

rut,

Leb

anon

Ans

an,

Sou

th

Kor

ea

Nan

-Chi

e an

d K

aohs

iung

, T

aiw

an

Bei

jing

, C

hina

S

hang

hai,

Chi

na

Hyy

tiälä

, F

inla

nd

Rom

e, I

taly

Sao

pa

ulo,

B

razi

l

Nit

erói

, B

razi

l

Sam

plin

g pe

riod

June

–201

2 to

May

–20

13

May

–Nov

, 20

13

(exc

ept

Oct

)

Mon

soon

, 19

95–9

6

Mar

ch–

June

, 20

06

Aug

, 20

00

Aug

–S

ep,

2004

Aug

–Sep

, 20

05

May

–Dec

, 20

06

Aug

, 20

06

Jan–

Oct

, 20

07

Mar

ch–

Apr

il,

2003

June

–Jul

y,

1994

, Jan

–M

arch

19

95

Win

ter,

19

99

Jan,

20

10

Nat

ure

of th

e st

udy

area

H

ill s

tatio

n,

Him

alay

an

Him

alay

an

Rur

al

Urb

an

Indu

stri

al

harb

orU

rban

In

dust

rial

co

mpl

ex

Urb

an

Urb

an

Urb

an

For

est

Urb

an

Urb

an

Urb

an

For

mal

dehy

de 1

1.6

2.4

1.7

19.8

2.

84.

624

.118

.335

.719

.4

0.5

21

21.1

3.

2 A

ceta

ldeh

yde

20.7

1.

2 -

13.7

2.

32.

136

.614

.915

.715

.9

0.4

17

34

6.5

Ace

tone

81

.6

7.7

-5

38.4

1.2

8.8

11.9

1.

3 16

.1-

7.3

Ref

eren

ce

Thi

s st

udy

Wan

gchu

k et

al.,

201

5 K

hare

et

al.,

1997

D

utta

et

al.,

2009

Ohu

ra e

t al

., 20

06

Mou

ssa

et a

l.,

2006

Kim

et

al.,

2008

Wan

g et

al.,

2010

Dua

n et

al

., 20

12

Hua

ng e

t al

., 20

08

Hel

l’en

et

al.,

20

04

Pos

sanz

ini

et a

l., 1

996

Mon

tero

et

al.,

20

01

Och

s et

al

., 20

11

Page 8: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2213

Fig. 4. Seasonal variation of total carbonyl compounds shown in box-whisker plot. (Lower, middle and upper boundary of the box correspond to 25, median and 75 percentile; lower and upper whiskers correspond to minimum and maximum and square in the box correspond to mean).

Table 4. Average concentrations of individual carbonyl compounds during different seasons.

Monsoon ± SD Postmonsoon ± SD Winter ± SD Premonsoon ± SD Formaldehyde 9.1 9.8 32.5 26.8 3.5 2.8 7.1 4.4 Acetaldehyde 13.3 26.4 79.2 84.3 2.2 2.1 5.0 4.8 Acetone 90.8 64.4 138.4 78.1 57.6 41.5 51.9 33.7 Propionaldehyde 0.9 0.6 10.2 9.2 0.7 0.7 2.0 2.2 Butanaldehyde 4.9 4.7 56.1 40.2 0.3 0.3 0.7 0.8 Crotonaldehyde 23.6 42.6 22.0 33.2 4.6 5.4 2.2 2.0 Benzaldehyde 10.2 18.3 31.8 51.7 0.2 0.2 0.5 0.8 iso-Valeraldehyde 2.4 4.2 24.3 19.8 0.9 0.8 1.0 0.7 Valeraldehyde 6.4 6.0 28.8 25.1 1.4 1.2 1.0 0.7 o-Tolaldehyde 0.3 0.9 4.3 7.6 0.1 0.1 0.2 0.4 p-Tolualdehyde 3.4 6.2 6.4 7.8 0.4 0.3 1.3 1.5 Hexanaldehyde 2.9 4.5 22.8 18.3 0.3 0.5 1.1 0.9

incoming solar radiation measurement sensor (Pyranometer: Model No. SP Lite2, Kipp and Zonen) and rain gauge (TE525 tipping bucket rain gauge). The data was recorded at an interval of 1 min. We have used the meteorological data of those days when the carbonyl sampling was done.

In order to evaluate the effect of meteorology on the atmospheric loading of CCs over Darjeeling we have made correlation analysis between each CCs concentration with meteorological parameters like temperature, wind speed, wind direction, total rainfall and solar radiation. But except temperature and solar radiation we did not find any significant correlation between carbonyls and other meteorological parameters. The correlation coefficients (R) for temperature and solar radiation are given in Table 5 where values higher than 0.5 are given in bold numbers.

Surface reaching solar radiation is the most important meteorological parameter that governs the photochemistry of oxygenated and non-oxygenated volatile organic compounds in the lower troposphere. UV radiation (λ < 320 nm) photo

chemically decomposes ozone to produce OH radical. O3 + hν (≤ 320 nm) → O(1D) + O2 (1) O(1D) + M → O(3P) + M (where M is O2 or N2) (2) O(1D) + H2O → 2OH· (radical) (3)

The OH radical then reacts with different VOCs to produce carbonyl compounds by the following equations (Koppmann, 2008).

CH4 + OH· 2O , NO HCHO (Formaldehyde) + H2O + NO2 (4)

C2H6 + OH· 2O , NO CH3CHO (Acetaldehyde) + H2O + NO2 (5)

Page 9: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2214

Hence surface reaching solar radiation is a major factor for secondary production of CCs in the atmosphere. Again, photolysis of CCs by surface reaching solar radiation is also the main sink of the atmospheric CCs (Koppmann, 2008). The following equations show the photolysis pathway of formaldehyde, acetaldehyde and acetone in the atmosphere.

HCHO + hν → H• + HCO• → H2 + CO (6) CH3CHO + hν → CH3• + HCO• → CH4 + CO → CH3CO• + H (7) CH3COCH3 + hν → CH3• + CH3CO• → 2CH3• + CO (8)

We observed good/moderate negative correlations between surface reaching solar radiation flux and concentration of all CCs during premonsoon season suggesting photochemical degradation of CCs under higher surface reaching radiation. But we did not observe such negative correlations during other seasons. Thus it could be said that the CCs during postmonsoon, winter and monsoon are not photochemically degraded over Darjeeling. However we observed good positive correlations between surface reaching solar radiation flux and concentration of formaldehyde and acetaldehyde. We observed higher methane concentration during postmonsoon over the same sampling site and in the same year (Ganesan et al., 2013) and thus it could suggest the photochemical formation of formaldehyde from methane through the reactions 1–4 stated above. The measurement of ethane was beyond our scope and thus the photochemical formation mechanism of acetaldehyde remain unexplained. However in detail studies are needed to establish the photochemical formation/degradation of carbonyls.

Good positive correlations between ambient temperature and CCs were observed during postmonsoon. It was reported in few studies that increasing temperature increases the emission of formaldehyde and acetaldehyde from wood and building materials (Parthsarathy et al., 2011; Suzuki et al., 2014). Thus temperature also played positive role to building up the CC concentration over the hill station during postmonsoon season.

However, we did not observe any significant correlations between CCs and meteorological parameters during monsoon season. Formaldehyde to Acetaldehyde Ratio (C1/C2)

The concentration ratio of formaldehyde to acetaldehyde (C1/C2) can be used to define the nature of the study area. Over the urban areas, the ratio varies between 1–2 (Feng et al., 2005; Cerón et al., 2007; Huang et al., 2008; Dutta et al., 2009; de M. Ochs et al., 2011), whereas over the rural areas the ratio becomes much higher, ~10 (Shepson et al., 1991; Possanzini et al., 1996). The probable reason for different C1/C2 ratios is the different sources of formaldehyde and acetaldehyde over different areas. Formaldehyde and acetaldehyde both have similar anthropogenic sources over urban areas like combustion of fossil fuel. But over rural areas, formaldehyde has significant biogenic source (Possanzini et al., 1996) whereas acetaldehyde does not.

Tab

le 5

. Cor

rela

tion

coe

ffic

ient

s (R

) of

car

bony

l com

poun

ds w

ith

tem

pera

ture

and

sol

ar r

adia

tion

dur

ing

diff

eren

t sea

sons

.

F

orm

al

Ace

tal

Ace

tone

Pro

pana

l B

utan

al

Cro

tona

l B

enza

lIs

oval

V

aler

alo-

Tol

p-T

olH

exan

alT

otal

M

onso

on

Day

T

emp

–0.0

8 0.

28

–0.2

0 0.

20

–0.2

4 –0

.09

–0.3

6 –0

.23

–0.4

1 0.

20

0.20

–0

.14

–0.1

2 S

R

0.16

0.

29

–0.0

5 0.

10

0.08

0.

06–0

.30

0.14

–0.3

9 –0

.03

0.46

0.

07

0.05

N

ight

T

emp

–0.4

9 –0

.15

–0.1

2 –0

.07

–0.0

8 –0

.16

–0.4

7 –0

.03

0.13

0.00

–0

.04

0.07

–0

.25

Pos

tmon

soon

D

ay

Tem

p 0.

83

0.70

0.

44

–0.0

1 0.

38

–0.2

0–0

.35

0.23

0.65

0.

53

0.45

0.

39

0.64

S

R

0.75

0.

73

–0.1

4 –0

.21

–0.3

2 –0

.12

–0.0

8 0.

160.

37–0

.35

0.44

0.

38

0.38

N

ight

T

emp

0.77

0.

45

0.29

–0

.21

0.04

0.

01

–0.3

9 0.

75

–0.0

9 0.

39

–0.5

6 –0

.13

0.29

W

inte

r D

ay

Tem

p 0.

18

0.29

0.

07

–0.5

0 0.

07

–0.5

6 –0

.40

–0.4

8–0

.16

–0.3

9–0

.05

–0.2

80.

00

SR

0.

11

0.51

0.

44

–0.0

3 0.

01

–0.3

60.

24

0.30

0.24

0.27

0.

09

–0.2

80.

44

Nig

ht

Tem

p –0

.60

–0.1

5 0.

04

–0.3

0 0.

16

–0.6

0 0.

12

–0.0

3–0

.23

0.82

–0

.21

–0.0

5–0

.05

Pre

mon

soon

D

ay

Tem

p 0.

15

0.21

–0

.28

0.10

–0

.47

–0.0

40.

04

–0.1

10.

25–0

.50

0.50

0.

69

0.06

S

R

–0.5

9 –0

.62

–0.3

0 –0

.04

–0.3

1 –0

.52

–0.3

4 –0

.34

–0.4

0 0.

38

–0.4

5 –0

.34

–0.4

1 N

ight

T

emp

0.20

0.

29

0.56

0.

41

–0.4

5 –0

.27

–0.1

3 0.

51

0.08

0.14

0.

58

0.65

0.

60

Page 10: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2215

Fig. 5. Day/Night ratio of carbonyl compounds during different seasons.

Hence over the urban areas the ratio varies between 1–2 when formaldehyde and acetaldehyde both have similar anthropogenic sources whereas the ratio increases significantly when rural biogenic sources play major role in production of formaldehyde and not acetaldehyde.

In the present study the annual average C1/C2 ratio value over Darjeeling was 1.64 ± 1.43 and comparable with that over an Indian mega city Kolkata (Day time: 1.32–1.59, night time: 1.44–1.89, Dutta et al., 2009). Thus Darjeeling represents a typical urban atmosphere over this part of Himalaya. The ratio values for different seasons are given in Table 6. The ratio was found to be maximum during premonsoon followed by winter, monsoon and postmonsoon. The higher C1/C2 ratio could be explained by the higher photochemical activities during the premonsoon season. The higher reaction rate constant of acetaldehyde (1.5 × 1011 cm3 molecule–1 s–1) than formaldehyde (0.9 × 1011 cm3 molecule–1 s–1) with respect to OH radical reaction (Atkinson and Arey, 2003) could be responsible for more photochemical loss of acetaldehyde during premonsoon and hence the ratio increased in this season. On the other hand, acetaldehyde generated from vehicular emissions could not get photochemically degraded during postmonsoon as the solar radiation flux was much lower in this season and hence the ratio value decreased to 0.7.

In addition to vehicular emission and secondary photochemical production biomass burning is also a significant source of formaldehyde and acetaldehyde as reported in Akagi et al. (2011), Lee et al. (1997), Holzinger et al. (1999) and several other studies. During colder seasons the anthropogenic activities like vehicular emission decreased due to the least tourist activity in this season. But biomass

burning is a common practice by the local people in winter producing more formaldehyde than acetaldehyde (Holzinger et al., 1999; Akagi et al., 2011) in the atmosphere which actually elevated C1/C2 ratio during the winter season. Thus, in winter, direct emission of formaldehyde from biomass burning dominates over the meteorological control like solar radiation fluxes. Acetaldehyde to Propanaldehyde Ratio (C2/C3)

The ratios of acetaldehyde to propanaldehyde (C2/C3) on annual and seasonal scale are given in Table 6. C2/C3 ratio is often used to determine whether the carbonyls have any anthropogenic sources (Ho et al., 2002b) as propanaldehyde is mainly considered to be anthropogenic origin (Anderson et al., 1996). This ratio is found to be higher in rural atmosphere than the urban atmospheres (Feng et al., 2005). Maximum value for C2/C3 ratio was observed during postmonsoon followed by monsoon, winter and premonsoon. The annual average C2/C3 ratio in the present study was 8.8 ± 11.8 which was higher than the study made over Kolkata (5.05, Dutta et al., 2009). However the minimum ratio value in premonsoon over Darjeeling (even lower than annual average ratio over Kolkata) indicates the major anthropogenic activities during this season. Source Apportionment Using Positive Matrix Factorization Model (PMF)

PMF (version 5.0) was used to apportion the sources of carbonyl compounds over Darjeeling. Positive Matrix Factorization has been extensively used to identify the sources of pollutants in many studies in recent years (e.g., Jorquera and Rappengluck, 2004; Latella et al., 2005; Xie

Page 11: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2216

Table 6. The average formaldehyde/acetaldehyde (C1/C2) and acetaldehyde/propanaldehyde (C2/C3) ratios during different season.

 C1/C2 C2/C3

Mean ± SD Mean ± SD Monsoon 1.10 1.16 10.00 10.51 Postmonsoon 0.70 0.77 15.71 20.83 Winter 1.84 1.02 6.52 9.17 Premonsoon 2.52 1.75 3.77 2.38 Annual 1.64 1.43 8.76 11.78

and Berkowitz, 2006; Brown et al., 2007; Song et al., 2007; Yuan et al., 2009). This source apportionment model does not require any prior knowledge on the exact emission profiles of the individual pollutants but source contributions can be determined based on observations at the receptor site. In brief, a parameter Q is calculated by the model which is the measure of the goodness of the fit parameters. Q (true) is calculated including all data points where Q (robust) is calculated excluding data points not fit by the model. The difference between Q (true) and Q (robust) is the measure of the impact of data points with high scaled residuals whereas high uncertainties result in similar Q values. Q (robust) is used to choose the optimal run from the multiple runs. If the number of sources is estimated properly, Q (robust) will have little variability between the runs and thus we get stable Q. The variability of Q (robust) is the indication of the variability of initial base run results. More details about the PMF model can be found in several other studies (Paatero and Tapper, 1994; Paatero, 1997; Reff et al., 2007).

Six factors were resolved from the analysis. The percentage contribution of each factor and the percentage contribution of each CC in all the factors are shown in Fig. 6.

Factor 1: The first factor was loaded mainly with acetone (92%) and the contributions of other CCs were almost negligible in this factor. We have already discussed in the previous section about the formation of acetone from its precursor VOC emitted from the tea plants. Thus this factor could be attributed to the biogenic emission from tea gardens. The contribution of this factor was maximum (43%) amongst all the other sources.

Factor 2: The contribution of factor 2 was 16%. Maximum loading of butanaldehyde and hexanaldehyde were present in this factor. A significant amount of acetone was also present in this factor. Several studies have identified the presence of these CCs in diesel exhaust emissions (Grosjean et al., 2001; Kean et al., 2001; Jacober et al., 2008; Dong et al., 2014). High number of heavy duty trucks and 4 wheelers driven by diesel engines are the major transports over Darjeeling. Thus diesel driven vehicular emission could be a significant source of CCs over the study area. Factor 2 could be attributed to the diesel exhaust emissions.

Factor 3: This factor had the third highest (15%) contribution among the sources identified by PMF. Crotonaldehyde and acetone were the major CCs of this factor. Both crotonaldehyde and acetone are generated from industrial sources (Graedel et al., 1986; Grosjean 1990; Salgado, 2008). Many industries like paint manufacturing

and petrochemical refinery uses acetone as the chemical solvent (Battig and Grandjean, 1964; Feng et al., 2005). Though these types of industries are not present over the hill station but these CCs could be transported (due to their higher atmospheric lifetime) to the study region from the polluted foot hill regions and other IGP regions which is the main industrial region in India. Thus factor 3 could be attributed to the industrial source.

Factor 4: This factor has similar contribution (15%) like factor 2 containing formaldehyde, acetaldehyde and hexanaldehyde mostly. Formaldehyde, acetaldehyde and hexanaldehyde have the common primary source of gasoline driven vehicular emission (Corrêa et al., 2003; Jacobson, 2007; Tunsaringkarn et al., 2012). In addition to the vehicular emission, formaldehyde and acetaldehyde could also be produced from secondary photochemical production from other primary VOCs. Thus it is not possible to assign factor 4 to a single particular source rather we have assigned factor 4 as gasoline exhaust/secondary production.

Factor 5: Contribution of Factor 5 was much less (8%) compared to the other factors. Iso-valeraldehyde and valeraldehyde are the major CCs present in this factor. These compounds are emitted during the processing/ fermentation (oxidation) of tea leaves, the intermediate step for the preparation of black tea and oolong tea (Runeckles and Tso, 2013; Panda, 2016). There are many tea processing industries present over Darjeeling. Thus factor 5 could be attributed to the tea processing industries.

Factor 6: Propanaldehyde is the main contributing CC in this factor along with hexanaldehyde and acetone. This factor contributes only 3% to the total CC concentration. Propanaldehyde, hexanaldehyde and acetone together could be emitted from the waste dumping areas as reported by many other studies (US DHHS, 1993; Eitzer, 1995; Jay and Stieglitz, 1995, Majumdar et al., 2014). It has been discussed in the previous sections that the sewage system is very poor and the waste disposal system is very unscientific over Darjeeling. Thus factor 6 could be attributed to the human waste disposal. Photochemical Reactivity and Ozone Formation Potential of Carbonyls

CCs react with OH radical to produce tropospheric ozone depending on the availability of NOx and non-methane VOC (NMVOC) concentrations (Carter, 1994). The propene-equivalent concentration (prop-equiv) (Chameides et al., 1992; So and Wang, 2004) and the maximum incremental reactivity (MIR) (Carter, 1994) scale were introduced to

Page 12: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2217

Fig. 6. (a) Percentage contribution of different sources for CCs over Darjeeling (b) Source profiles of several factors estimated from PMF model.

estimate the reactivity of each CC towards OH radical reaction and to calculate the contribution of each CC for the photochemical production of tropospheric O3. The propane equivalent concentration of a carbonyl compound is given by

OH

OH 3 6

k jprop-equiv j conc j

k C H (9)

where prop-equiv (j) is the propane equivalent concentration of jth carbonyl species expressed in µg m–3, conc (j) is the concentration of jth carbonyl species expressed in µg m–3, kOH(j) is the rate constant of jth carbonyl species with respect to OH radical reaction. The individual rate constant values are taken from Atkinson and Arey (2003). kOH(C3H6) is the OH rate constant of propane.

The reaction with the OH radical is the primary step for the oxidation of CCs, thus prop-equiv concentration is an

indicator of the relative contribution made by each CC to the total flux of VOCs being oxidized in the area. The prop-equiv concentration of each carbonyl species are summarized in Table 7.

The significance of prop-equiv (j) of jth carbonyl species lies in the measure of the concentration of CC(j) on an OH-reactivity based scale normalized to the reactivity of propene. The total prop-equiv concentration for all the CCs measured varied over a wide range from 368 µg m–3 during winter to 6252 µg m–3 during postmonsoon season. The mean annual average prop-equiv concentration was 1945 µg m–3. Crotonaldehyde had the maximum annual pro-equiv concentration of 578.5 µg m–3 accounting for ~30% of total concentration followed by acetaldehyde (15%) and butanaldehyde (14%). It could be noted that though acetone had the maximum concentration in all the seasons but due to very low reactivity toward OH radical, acetone had the least prop-equiv concentration of 12.7 µg m–3. This accounted for only < 1% to the total concentration.

Page 13: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2218

Table 7. The propane equivalent concentration of total and individual CCs during different seasons.

1012 × KOH Annual Monsoon Postmonsoon Winter Premonsoon Formaldehyde 9.37 89.6 78.2 279.2 30.1 61.4 Acetaldehyde 15.00 285.0 183.6 1090.2 29.8 68.4 Acetone 0.17 31.3 14.2 122.6 9.0 8.1 Propanaldehyde 20.00 52.7 17.1 187.1 12.9 36.2 Butanaldehyde 24.00 267.9 108.4 1236.1 5.9 16.1 Crotonaldehyde 48.00 578.5 1039.0 970.8 204.6 95.0 Benzaldehyde 12.00 102.5 112.2 350.1 1.9 5.9 Isovaleraldehyde 27.00 141.4 59.0 602.8 23.0 23.8 Valeraldehyde 28.00 204.4 165.2 740.8 35.7 26.1 o-Tolualdehyde 18.00 16.1 5.4 70.7 1.2 4.0 p-Tolualdehyde 13.00 32.5 40.7 75.8 4.7 15.9 Hexanaldehyde 30.00 151.3 79.4 627.3 9.4 30.2 Total 1953.4 1902.4 6353.4 368.4 391.1

O3 formation potential (OFP) of CCs can be calculated as the product of the MIR coefficient of a CC and the concentration of that CC. The maximum incremental reactivity (MIR) of a carbonyl is operationally defined as the change in maximum ozone production by a small increase in emissions of that carbonyl species. MIR coefficient is dimensionless and can be defined as gram of O3 produced per gram of carbonyl species. MIR coefficients for CCs were taken from Carter (2007). OFP was thus calculated as follows:

O3 formation potential (j) = Concentration (j) × MIR coefficient (j) (10)

Table 8 represents the OFP values of each CC in µg m–3.

The total annual average OFP was found to be 558 µg m–3

and this high value of OFP indicates that carbonyl compounds could play significant role in tropospheric ozone formation over Darjeeling. The OFP varied from 127 µg m–3 during premonsoon to 1770 µg m–3 during postmonsoon season suggesting more secondary tropospheric ozone formation through the oxidation of volatile organic compounds during postmonsoon. Acetaldehyde had the maximum OFP followed by crotonaldehyde and formaldehyde. Acetaldehyde contributed 23.5% of total ozone formation followed by crotonaldehyde (21.5%) and formaldehyde (19.3%). Rest of the carbonyls contributed ~35% to the total tropospheric ozone formation. Hence, we observed that acetaldehyde and crotonaldehyde had higher propane equivalent concentration as well as higher OFP indicating that they are the most potential carbonyl compounds to form tropospheric ozone over Darjeeling. CONCLUSIONS

The first ever year-long study on atmospheric carbonyl compounds over Indian Himalaya was conducted where samples were collected from a high altitude station, Darjeeling (27.01°N, 88.15°E, 2200 masl). The major findings of the study are as follows: ● Major carbonyl compounds like formaldehyde and

acetaldehyde over Darjeeling were found to be comparable

and higher than some of the metro cities in India and other parts of the world. Carbonyl compound concentrations over Darjeeling were also higher than Bhutan, a high altitude state at eastern Himalaya.

● Acetone over Darjeeling had the maximum abundance and its concentration was much higher compared to other carbonyls. It was even higher than that over other parts of the world. In addition to various anthropogenic sources, biogenic emissions of volatile organic compounds from tea plants followed by their photochemical transformation also played significant role in building up such high acetone in the atmosphere of Darjeeling.

● Carbonyl compounds showed distinct seasonal variation with the maximum concentrations during postmonsoon followed by monsoon, winetr~premonsoon. Surface reaching solar radiation played major role in photochemical production/degradation of carbonyl compounds over Darjeeling. In spite of higher vehicular activities, photochemical degradation due to higher solar radiation flux reduced the atmospheric loading of carbonyl compounds during premonsoon. On the other hand, low solar radiation along with higher vehicular activities helped in accumulating carbonyl compounds during postmonsoon.

● Most of the carbonyl compounds showed higher concentrations during day time except acetone and butanaldehyde during winter and postmonsoon. Biomass/wood burning by the local habitants of Darjeeling during colder nights could enhance the night time concentrations of acetone and butanaldehyde.

● Formaldehyde to acetaldehyde (C1/C2) concentration ratio over Darjeeling represents Darjeeling as a typical urban station at eastern part of Himalaya. Acetaldehyde to propanaldehyde (C2/C3) ratio showed higher anthropogenic activities during premonsoon.

● Higher propane equivalent concentration and ozone formation potential suggests that more tropospheric ozone could be produced through the oxidation of volatile organic compounds during postmonsoon and acetaldehyde and crotonaldehyde were the potential carbonyl compounds to form tropospheric ozone.

● Six major sources of carbonyl compounds were identified by PMF model analysis. Biogenic emissions from tea

Page 14: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2219

Table 8. Maximum incremental reactivity and ozone formation potential of individual carbonyl species.

Ozone formation (µg m–3) MIR Annual Monsoon Postmonsoon Winter Premonsoon

Formaldehyde 9.24 107.5 84.1 300.1 32.3 66.0 Acetaldehyde 6.34 131.3 84.6 502.3 13.7 31.5 Acetone 0.35 28.6 31.8 48.5 20.2 18.2 Propanaldehyde 6.83 19.6 6.4 69.6 4.8 13.5 Butanaldehyde 5.75 70.0 28.3 322.8 1.5 4.2 Crotonaldehyde 9.14 120.1 215.6 201.5 42.5 19.7 Benzaldehyde –0.67 –6.2 –6.8 –21.3 –0.1 –0.4 Isovaleraldehyde 4.79 27.3 11.4 116.6 4.4 4.6 Valeraldehyde 4.89 38.9 31.4 141.0 6.8 5.0 o-Tolualdehyde –0.54 –0.5 –0.2 –2.3 0.0 –0.1 p-Tolualdehyde –0.54 –1.5 –1.8 –3.4 –0.2 –0.7 Hexanaldehyde 4.18 23.0 12.1 95.3 1.4 4.6 Total 558.0 496.8 1770.5 127.4 166.1

plantation and vehicular emission (gasoline and diesel exhaust) were the major sources of carbonyl compounds over Darjeeling. The high atmospheric loading of carbonyl compounds

over Darjeeling is of a serious concern from the point of view of human health and sensitive ecosystem over this part of Indian Himalaya. Thus it is imperative to control vehicular pollution, solvent usage, adopt proper and scientific waste management system etc which are the major causes of carbonyls and other gaseous and particulate carbonaceous pollutants as well over Darjeeling. ACKNOWLEDGEMENTS

Authors would like to thank Science and Engineering Council, Department of Science and Technology, Government of India for supporting the study under IRHPA (Intensification of Research in High Priority Areas) scheme. Authors also like to thank University Grant Commission, Govt of India for providing fellowship to C. Sarkar (first author). The study was undertaken by Bose Institute in collaboration with CSIR-NEERI. Field study conducted by Bose Institute and sample analysis performed by CSIR-NEERI, Kolkata Zonal Laboratory. Thanks are due to Sabyasachi Majee, Bose Institute for his consistent support in sampling and Debasish Sengupta and Pamela Chowdhury, NEERI for their assistance during sample analysis. Authors would also like to thank D. K. Roy, Bose Institute for his overall logistic support. REFERENCES Adak, A., Chatterjee, A., Singh, A.K., Sarkar, C., Ghosh,

S. and Raha, S. (2014). Atmospheric fine mode particulates at eastern Himalaya, India: Role of meteorology, long-range transport and local anthropogenic sources. Aerosol Air Qual. Res. 14: 440–450.

Akagi, S.K., Yokelson, R.J., Wiedinmyer, C., Alvarado, M.J., Reid, J.S., Karl, T., Crounse, J.D. and Wennberg, P.O. (2011). Emission factors for open and domestic biomass burning for use in atmospheric models. Atmos.

Chem. Phys. 11: 4039–4072. Anderson, L.G., Lanning, J.A., Barrell, R., Miyagishima,

J., Jones, R.H. and Wolfe, P. (1996). Sources and sinks of formaldehyde and acetaldehyde: An analysis of Denver's ambient concentration data. Atmos. Environ. 30: 2113–2123.

Andreini, B.P., Baroni, R., Galimberti, E. and Sesana, G. (2000). Aldehydes in the atmospheric environment: evaluation of human exposure in the north-west area of Milan. Microchem. J. 67: 11–19.

Atkinson, R. (1986). Kinetics and mechanisms of the gas-phase reactions of the hydroxyl radical with organic compounds under atmospheric conditions. Chem. Rev. 86: 69–201.

Atkinson, R. (1990). Gas-phase tropospheric chemistry of organic compounds: A review. Atmos. Environ. 24: 1–41.

Atkinson, R. (2000). Atmospheric chemistry of VOCs and NOx. Atmos. Environ. 34: 2063–2101.

Atkinson, R. and Arey, J. (2003). Atmospheric degradation of volatile organic compounds. Chem. Rev. 103: 4605–4638.

Bakeas, E.B., Argyris, D.I. and Siskos, P.A. (2003). Carbonyl compounds in the urban environment of Athens, Greece. Chemosphere 52: 805–813.

Battig, K. and Grandjean, E. (1964). Industrial solvents and avoidance conditioning in rats. Arch. Environ. Health 9: 745–749.

Brown, S.G., Frankel, A. and Hafner, H R. (2007). Source apportionment of VOCs in the Los Angeles area using Positive Matrix Factorization. Atmos. Environ. 41: 227–37.

Carter, W.P. (1994). Development of ozone reactivity scales for volatile organic compounds. J. Air Waste Manage. Assoc. 44: 881–899.

Carter, W.P. (2007). Development of the SAPRC-07 chemical mechanism and updated ozone reactivity scales. California Air Resources Board, Research Division.

Cavalcante, R.M., Campelo, C.S., Barbosa, M.J., Silveira, E.R., Carvalho, T.V. and Nascimento, R.F. (2006). Determination of carbonyl compounds in air and cancer risk assessment in an academic institute in Fortaleza,

Page 15: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2220

Brazil. Atmos. Environ. 40: 5701–5711. CEPA (California Environmental Protection Agency)

(1993). Acetaldehyde as a toxic air contaminant. Part A: Exposure; Part B: Health assessment. Air Resources Board, Stationary Source Division, Sacramento, CA, USA.

Cerón, R.M., Cerón, J.G. and Muriel, M. (2007). Diurnal and seasonal trends in carbonyl levels in a semi-urban coastal site in the Gulf of Campeche, Mexico. Atmos. Environ. 41: 63–71.

Chameides, W.L., Fehsenfeld, F., Rodgers, M.O., Cardelino, C., Martinez, J., Parrish, D., Lonneman, W., Lawson, D.R., Rasmussen, R.A., Zimmerman, P., Greenberg, J., Middleton, P. and Wang, T. (1992). Ozone precursor relationships in the ambient atmosphere. J. Geophys. Res. 97: 6037–6055.

Chatfield, R.B., Gardner, E.P. and Calvert, J.G. (1987). Sources and sinks of acetone in the troposphere: Behavior of reactive hydrocarbons and a stable product. J. Geophys. Res. 92: 4208–4216.

Chatterjee, A., Adak, A., Singh, A.K., Srivastava, M.K., Ghosh, S.K., Tiwari, S., Devara, P.C. and Raha, S. (2010). Aerosol chemistry over a high altitude station at north eastern Himalayas, India. PLoS One 5: 11122.

Chatterjee, A., Ghosh, S.K., Adak, A., Singh, A.K., Devara, P.C. and Raha, S. (2012). Effect of dust and anthropogenic aerosols on columnar aerosol optical properties over Darjeeling (2200 m asl), eastern Himalayas, India. PLoS One 7: e40286.

Chen, W.T., Shao, M., Lu, S.H., Wang, M., Zeng, L.M., Yuan, B. and Liu, Y. (2014). Understanding primary and secondary sources of ambient carbonyl compounds in Beijing using the PMF model. Atmos. Chem. Phys. 14: 3047–3062.

Christensen, C.S., Skov, H., Nielsen, T. and Lohse, C. (2000). Temporal variation of carbonyl compound concentrations at a semi-rural site in Denmark. Atmos. Environ. 34: 287–296.

Corrêa, S.M., Martins, E.M. and Arbilla, G. (2003). Formaldehyde and acetaldehyde in a high traffic street of Rio de Janeiro, Brazil. Atmos. Environ. 37: 23–29.

Crutzen, P.J., Delany, A.C., Greenberg, J., Haagenson, P., Heidt, L., Lueb, R., Pollock, W., Seiler, W., Wartburg, A. and Zimmerman, P. (1985). Tropospheric chemical composition measurements in Brazil during the dry season. J. Atmos. Chem. 2: 233–256.

de M. Ochs, S., Albuquerque, F.C., Pontes Massa, M.C.G. and Pereira Netto, A.D. (2011). Evaluation of C1–C13 carbonyl compounds by RRLC-UV in the atmosphere of Niterói City, Brazil. Atmos. Environ. 45: 5183–5190.

Dong, D., Shao, M., Li, Y., Lu, S., Wang, Y., Ji, Z. and Tang, D. (2014). Carbonyl emissions from heavy-duty diesel vehicle exhaust in China and the contribution to ozone formation potential. J. Environ. Sci. 26: 122–128.

Duan, J., Guo, S., Tan, J., Wang, S. and Chai, F. (2012). Characteristics of atmospheric carbonyls during haze days in Beijing, China. Atmos. Res. 114: 17–27.

Dutta, C., Som, D., Chatterjee, A., Mukherjee, A.K., Jana, T.K. and Sen, S. (2009). Mixing ratios of carbonyls and BTEX in ambient air of Kolkata, India and their associated

health risk. Environ. Monit. Assess. 148: 97–107. Eitzer, B.D. (1995). Emissions of volatile organic chemicals

from municipal solid waste composting facilities. Environ. Sci. Technol. 29: 896–902.

Feng, Y., Wen, S., Chen, Y., Wang, X., Lü, H., Bi, X., Sheng, G. and Fu, J. (2005). Ambient Levels of carbonyl Compounds and Their Sources in Guangzhou, China. Atmos. Environ. 39: 1789–1800.

Finlayson-Pitts, B.J. and Pitts, J.N.J. (2000). Chemistry of the Upper and Lower Atmosphere. Academic Press, San Diego.

Friedrich, R. and Obermeier, A. (1999). Anthropogenic emissions of volatile organic compounds. In Reactive Hydrocarbons in the Atmosphere, Hewitt, N.C. (Ed.), Academic Press, San Diego, USA, pp. 1–39.

Fruekilde, P., Hjorth, J., Jensen, N.R., Kotzias, D. and Larsen, B. (1998). Ozonolysis at vegetation surfaces: A source of acetone, 4-oxopentanal, 6-methyl-5-hepten-2-one, and geranyl acetone in the troposphere. Atmos. Environ. 32: 1893–1902.

Ganesan, A.L., Chatterjee, A., Prinn, R.G., Harth, C.M., Salameh, P.K., Manning, A.J., Hall, B.D., Mühle, J., Meredith, L.K., Weiss, R.F. and O'Doherty, S. (2013). The variability of methane, nitrous oxide and sulfur hexafluoride in Northeast India. Atmos. Chem. Phys. 13: 10633–10644.

Graedel, T.E., Hawkins, D.T. and Claxton, L.D. (1986). Atmospheric Chemical Compounds. Sources, Occurrence, and Bioassay, Academic Press, San Diego, CA, pp. 295–314.

Grosjean, D. (1990). Atmospheric chemistry of toxic contaminants. 3. Unsaturated aliphatics: Acrolein, acrylonitrile, maleic anhydride. J. Air Waste Manage. Assoc. 40: 1664–1668.

Grosjean, D., Grosjean, E. and Gertler, A.W. (2001). On-road emissions of carbonyls from light-duty and heavy-duty vehicles. Environ. Sci. Technol. 35: 45–53.

Guo, S., Wen, S., Wang, X., Sheng, G., Fu, J., Hu, P. and Yu, Y. (2009). Carbon isotope analysis for source identification of atmospheric formaldehyde and acetaldehyde in Dinghushan Biosphere Reserve in South China. Atmos. Environ. 43: 3489–3495.

Hellén, H., Hakola, H., Reissell, A. and Ruuskanen, T.M. (2004). Carbonyl compounds in boreal coniferous forest air in Hyytiälä, Southern Finland. Atmos. Chem. Phys. 4: 1771–1780.

Ho, C.T., Zheng, X. and Li, S. (2015). Tea aroma formation. Food Sci. Hum. Wellness 4: 9–27.

Ho, K.F., Lee, S.C. and Chiu, G.M.Y. (2002a). Characterization of selected volatile organic compounds, polycyclic aromatic hydrocarbons and carbonyl compounds at a roadside monitoring station. Atmos. Environ. 36: 57–65.

Ho, K.F., Lee, S.C., Louie, P.K.K. and Zou, S.C. (2002b). Seasonal variation of carbonyl compound concentrations in urban area of Hong Kong. Atmos. Environ. 36: 1259–1265.

Ho, S.S.H., Ip, H.S.S., Ho, K.F., Dai, W.T., Cao, J. and Ng, L.P.T. (2013). Technical note: Concerns on the use

Page 16: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2221

of ozone scrubbers for gaseous carbonyl measurement by DNPH-coated silica gel cartridge. Aerosol Air Qual. Res. 13: 1151–1160.

Holzinger, R., Warneke, C., Hansel, A., Jordan, A., Lindinger, W., Scharffe, D.H., Schade, G. and Crutzen, P.J. (1999). Biomass burning as a source of formaldehyde, acetaldehyde, methanol, acetone, acetonitrile, and hydrogen cyanide. Geophys. Res. Lett. 26: 1161–1164.

Huang, J., Feng, Y., Li, J., Xiong, B., Feng, J., Wen, S., Sheng, G., Fu, J. and Wu, M. (2008). Characteristics of carbonyl compounds in ambient air of Shanghai, China. J. Atmos. Chem. 61: 1–20.

Jacob, D.J. (2000). Heterogeneous chemistry and tropospheric ozone. Atmos. Environ. 34: 2131–2160.

Jacobson, M.Z. (2007). Effects of ethanol (E85) versus gasoline vehicles on cancer and mortality in the United States. Environ. Sci. Technol. 41: 4150–4157.

Jakober, C.A., Robert, M.A., Riddle, S.G., Destaillats, H., Charles, M.J., Green, P.G. and Kleeman, M.J. (2008). Carbonyl emissions from gasoline and diesel motor vehicles. Environ. Sci. Technol. 42: 4697–4703.

Janson, R. and De Serves, C. (2001). Acetone and mono terpene emissions from the boreal forest in northern Europe. Atmos. Environ. 35: 4629–4637.

Jay, K. and Stieglitz, L. (1995). Identification and quantification of volatile organic components in emissions of waste incineration plants. Chemosphere 30: 1249–1260.

Jorquera, H. and Rappengluck, B. (2004). Receptor modeling of ambient VOC at Santiago, Chile. Atmos. Environ. 38: 4243–63.

Kean, A.J., Grosjean, E., Grosjean, D. and Harley, R.A. (2001). On-road measurement of carbonyls in California light-duty vehicle emissions. Environ. Sci. Technol. 35: 4198–4204.

Khare, P., Satsangi, G.S., Kumar, N., Kumari, K.M. and Srivastava, S.S. (1997). Surface measurements of formaldehyde and formic and acetic acids at a subtropical semiarid site in India. J. Geophys. Res. 102: 18997–19005.

Kim, K., Hong, Y., Pal, R., Jeon, E., Koo, Y. and Sunwoo, Y. (2008). Investigation of carbonyl compounds in air from various industrial emission sources. Chemosphere 70: 807–820.

Koppmann, R. (2008). Volatile organic compounds in the atmosphere. John Wiley & Sons.

Lary, D.J. and Shallcross, D.E. (2000). Central role of carbonyl compounds in atmospheric chemistry. J. Geophys. Res. 105: 19771–19778.

Latella, A., Stani, G., Cobelli, L., Duane, M., Junninen, H., Astorga, C. and Larsen, B.R. (2005). Semicontinuous GC analysis and receptor modelling for source apportionment of ozone precursor hydrocarbons in Bresso, Milan, 2003. J. Chromatogr. A 1071: 29–39.

Lee, M., Heikes, B.G., Jacob, D.J., Sachse, G. and Anderson, B. (1997). Hydrogen peroxide, organic hydroperoxide, and formaldehyde as primary pollutants from biomass burning. J. Geophys. Res. 102: 1301–1309.

Lee, S.C., Ho, K.F., Chan, L.Y., Zielinska, B. and Chow, J.C. (2001). Polycyclic aromatic hydrocarbons (PAHs)

and carbonyl compounds in urban atmosphere of Hong Kong. Atmos. Environ. 35: 5949–5960.

Lewis, A.C., Evans, M.J., Hopkins, J.R., Punjabi, S., Read, K.A., Purvis, R.M., Andrews, S.J., Moller, S.J., Carpenter, L.J., Lee, J.D., Rickard, A.R., Palmer, P.I. and Parrington, M. (2013). The influence of biomass burning on the global distribution of selected non-methane organic compounds. Atmos. Chem. Phys. 13: 851–867.

Macdonald, A.M., Makar, P.A., Anlauf, K.G., Hayden, K.L. and Bottenheim, J.W. (2001). Summertime formaldehyde at a highelevation site in Quebec. J. Geophys. Res. 106: 32361–32374.

Majumdar, D. and Srivastava, A. (2012). Volatile organic compound emissions from municipal solid waste disposal sites: A case study of Mumbai, India. J. Air Waste Manage. Assoc. 62: 398–407.

Majumdar, D., Ray, S., Chakraborty, S., Rao, P.S., Akolkar, A.B., Chowdhury, M. and Srivastava, A. (2014). Emission, speciation, and evaluation of impacts of non-methane volatile organic compounds from open dump site. J. Air Waste Manage. Assoc. 64: 834–845.

Montero, L., Vasconcellos, P.C., Souza, S.R., Pires, M.A.F., Sanchez-Ccoyllo, O.R., Andrade, M.F. and Carvalho, L.R.F. (2001). Measurements of atmospheric carboxylic acids and carbonyl compounds in São Paulo City, Brazil. Environ. Sci. Technol. 35: 3071–3081.

Moussa, S.G., El-Fadel, M. and Saliba, N.A. (2006). Seasonal, diurnal and nocturnal behaviors of lower carbonyl compounds in the urban environment of Beirut, Lebanon. Atmos. Environ. 40: 2459–2468.

Obermeyer, G., Aschmann, S.M., Atkinson, R. and Arey, J. (2009). Carbonyl atmospheric reaction products of aromatic hydrocarbons in ambient air. Atmos. Environ. 43: 3736–3744.

Ohura, T., Amagai, T. and Fusaya, M. (2006). Regional Assessment of Ambient Volatile Organic Compounds in an Industrial Harbor Area, Shizuoka, Japan. Atmos. Environ. 40: 238–248.

Orlando, J.J., Nozière, B., Tyndall, G.S., Orzechowska, G.E., Paulson, S.E. and Rudich, Y. (2000). Product studies of the OH- and ozone-initiated oxidation of some monoterpenes. J. Geophys. Res. 105: 11561–11572.

Paatero, P. and Tapper, U. (1994). Positive matrix factorization: A non-negative factor model with optimal utilization of error estimates of data values. Environmentrics 5: 111–126.

Paatero, P. (1997). Least squares formulation of robust non-negative factor analysis. Chemom. Intell. Lab. Syst. 37: 23–35.

Pal, R., Kim, K., Hong, Y. and Jeon, E. (2008). The pollution status of atmospheric carbonyls in a highly industrialized area. J. Hazard. Mater. 153: 1122–1135.

Panda, H. (2016). The complete book on cultivation and manufacture of tea. Asia Pacific Business Press Inc. pp. 497–498.

Pang, X. and Lee, X. (2010). Temporal variations of atmospheric carbonyls in urban ambient air and street canyons of a Mountainous city in Southwest China. Atmos. Environ. 44: 2098–2106.

Page 17: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017 2222

Parthasarathy, S., Maddalena, R.L., Russell, M.L. and Apte, M.G. (2011). Effect of temperature and humidity on formaldehyde emissions in temporary housing units. J. Air Waste Manage. Assoc. 61: 689–695.

Possanzini, M., Di Palo, V., Petricca, M., Fratarcageli, R. and Brocco, D. (1996). Measurements of lower carbonyls in Rome ambient air. Atmos. Environ. 30: 3757–3764.

Premachandra, B.R., Venkatakrishna, S. and Cama, H.R. (1977). Distribution of carotenoid pigments in tea leaves. Tea Q 47: 28–31.

Reff, A., Eberly, S.I. and Bhave, P.V. (2007). Receptor modeling of ambient particulate matter data using positive matrix factorization: Review of existing methods. J. Air Waste Manage. Assoc. 57: 146–154.

Reissell, A., Harry, C., Ascchmann, S.M., Atkinson, R. and Arey, J. (1999). Formation of acetone from the OH radical- and O3- initiated reactions of a series of monoterpenes. J. Geophys. Res. 109: 13869–13879.

Roy, A., Chatterjee, A., Tiwari, S., Sarkar, C., Das, S.K., Ghosh, S.K. and Raha, S. (2016). Precipitation chemistry over urban, rural and high altitude Himalayan stations in eastern India. Atmos. Res. 181: 44–53.

Runeckles, V.C. and Tso, T.C. (2013). Structural and Functional Aspects of Phytochemistry: Recent Advances in Phytochemistry (Vol. 5). Elsevier, pp. 296–298.

Salgado, M.S., Monedero, E., Villanueva, F., Martín, P., Tapia, A and Cabañas, B. (2008). Night-time atmospheric fate of acrolein and crotonaldehyde. Environ. Sci. Technol. 42: 2394–2400.

Sarkar, C., Chatterjee, A., Majumdar, D., Ghosh, S.K., Srivastava, A. and Raha, S. (2014). Volatile organic compounds over Eastern Himalaya, India: Temporal variation and source characterization using Positive Matrix Factorization. Atmos. Chem. Phys. Discus. 14: 32133–32175.

Sarkar, C., Chatterjee, A., Singh, A. K., Ghosh, S. and Raha, S. (2015). Characterization of black carbon aerosols over Darjeeling - A high altitude Himalayan station in eastern India. Aerosol Air Qual. Res. 15: 465–478.

Seco, R., Peñuelas, J. and Filellaet, I. (2007). Short-chain oxygenated VOCs: Emission and uptake by plants and atmospheric sources, sinks, and concentrations. Atmos. Environ. 41: 2477–2499.

Seinfeld, J.H. and Pandis, S.N. (1998). Atmospheric chemistry and physics: From air pollution to climate change. 2nd Edition, John Wiley & Sons, Inc, New Jersey.

Shepson, P.B., Hastie, D.R., Schiff, H.I., Polizzi, M., Bottenheim, J.W., Anlauf, K., Mackay, G.I. and Karecki, D.R. (1991). Atmospheric concentrations and temporal variations of C1-C3 carbonyl compounds at two rural sites in central Ontario. Atmos. Environ. 25: 2001–2015.

Singh, H.B., and Hanst, P.L. (1981). Peroxyacetylnitrate (PAN) in the unpolluted atmosphere: An important reservoir for nitrogen oxides. Geophys. Res. Lett. 8: 941–944.

Singh, H.B., O'hara, D., Herlth, D., Sachse, W., Blake, D.R., Bradshaw, J.D. and Crutzen, P.J. (1994). Acetone in the atmosphere: Distribution, sources, and sinks. J.

Geophys. Res. 99: 1805–1819. Singh, H.B., Kanakidou, M., Crutzen, P.J. and Jacob, D.J.

(1995). High concentrations and photochemical fate of oxygenated hydrocarbons in the global troposphere. Nature 378: 50–54.

Sirju, A.P. and Shepson, P.B. (1995). Laboratory and field investigation of the DNPH cartridge technique for the measurement of atmospheric carbonyl compounds. Environ. Sci. Technol. 29: 384–392.

Slemr, J. and Junkermann, W. (1996). Temporal variations in formaldehyde, acetaldehyde, and acetone and budget of formaldehyde at a rural site in Southern Germany. Atmos. Environ. 21: 3667–3676.

So, K.L. and Wang, T. (2004). C3–C12 non-methane hydrocarbons in subtropical Hong Kong: Spatial–temporal variations, source–receptor relationships and photochemical reactivity. Sci. Total. Environ. 328: 161–174.

Solberg, S., Dye, C., Walker, S.E., Simpson, D. (2001). Long-term measurements and model calculations of formaldehyde at rural European monitoring sites. Atmos. Environ. 35: 195–207.

Song, Y., Shao, M., Liu, Y., Lu, S. H., Kuster, W., Goldan, P. and Xie, S.D. (2007). Source apportionment of ambient volatile organic compounds in Beijing. Environ. Sci. Technol. 41: 4348–4353.

Suzuki, M., Akitsu, H., Miyamoto, K., Tohmura, S.I. and Inoue, A. (2014). Effects of time, temperature, and humidity on acetaldehyde emission from wood-based materials. J. Wood Sci. 60: 207–214.

Tanner, R.L., Zielinska, B., Uberna, E., Harshfield, G. and McNichol, A.P. (1996). Concentrations of carbonyl compounds and the carbon isotopy of Formaldehyde at a coastal site in Nova Scotia during the NARE summer intensive. J. Geophys. Res. 101: 28961–28970.

Tunsaringkarn, T., Prueksasit, T., Kitwattanavong, M., Siriwong, W., Sematong, S., Zapuang, K. and Rungsiyothin, A. (2012). Cancer risk analysis of benzene, formaldehyde and acetaldehyde on gasoline station workers. J. Environ. Eng. Ecol. Sci. 1: 1.

Uchiyama, S., Inaba, Y. and Kunugita, N. (2012). Ozone Removal in the Collection of Carbonyl Compounds in Air. J. Chromatogr. A 1229: 293–297.

US DHHS (1993). US Department of Health and Human Services. Hazardous substances data bank (HSDB, online database). National Toxicology Information Program, National Library of Medicine, Bethesda, MD.

USEPA (1996). Determination of Formaldehyde in Ambient Air Using Adsorbent Cartridge Followed by High Performance Liquid Chromatography (HPLC), EPA Method: TO-11A.

USEPA (1999). Compendium of Methods for the Determination of Toxic Organic Compounds in Ambient Air, Method TO-11A: Center for Environmental Research Information, Office of Research and Development, Cincinnati.

Wang, H.K., Huang, C.H., Chen, K.S., Peng, Y.P. and Lai, C.H. (2010). Measurement and source characteristics of carbonyl compounds in the atmosphere in Kaohsiung

Page 18: How the Atmosphere over Eastern Himalaya, India is Polluted ...Darjeeling) and ends at Darjeeling. Darjeeling Himalayan Railway (a world heritage site) is a narrow gauge railway runs

Sarkar et al., Aerosol and Air Quality Research, 17: 2206–2223, 2017  2223

city, Taiwan. J. Hazard. Mater. 179: 1115–1121. Wangchuk, T., He, C., Dudzinska, M.R. and Morawska, L.

(2015). Seasonal variations of outdoor air pollution and factors driving them in the school environment in rural Bhutan. Atmos. Environ. 113: 151–158.

World Health Organisation (WHO). Air quality guidelines for Europe, WHO European Series No. 23, Copenhagen.

Xie, Y.L. and Berkowitz, C.M. (2006). The use of positive matrix factorization with conditional probability functions in air quality studies: An application to hydrocarbon emissions in Houston, Texas. Atmos. Environ. 40: 3070–3091.

Yuan, B., Shao, M., Lu, S. and Wang, B. (2010). Source profiles of volatile organic compounds associated with solvent use in Beijing, China. Atmos. Environ. 44: 1919–1926.

Zhang, J. and Smith, K.R. (1999). Emissions of carbonyl compounds from various cookstoves in China. Environ. Sci. Technol. 33: 2311–2320.

Received for review, January 22, 2017 Revised, July 11, 2017

Accepted, July 12, 2017