<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JEP</journal-id><journal-title-group><journal-title>Journal of Environmental Protection</journal-title></journal-title-group><issn pub-type="epub">2152-2197</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jep.2012.37079</article-id><article-id pub-id-type="publisher-id">JEP-21141</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Temporal Trends in Ambient SO&lt;sub&gt;2&lt;/sub&gt; at a High Altitude Site in Semi-Arid Western India: Observations versus Chemical Transport Modeling
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>immy</surname><given-names>Francis</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Physical Research Laboratory, Ahmedabad, India</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>timmyf@prl.res.in</email></corresp></author-notes><pub-date pub-type="epub"><day>18</day><month>07</month><year>2012</year></pub-date><volume>03</volume><issue>07</issue><fpage>657</fpage><lpage>680</lpage><history><date date-type="received"><day>March</day>	<month>18th,</month>	<year>2012</year></date><date date-type="rev-recd"><day>April</day>	<month>19th,</month>	<year>2012</year>	</date><date date-type="accepted"><day>May</day>	<month>22nd,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Ambient sulphur dioxide (SO
  <sub>2</sub>) measurements have been performed at a high altitude site in the semi arid region of western India, Gurushikhar, Mt. Abu (24.6&#176;N, 72.7&#176;E, 1680 m ASL), during different sampling periods span over Sep-Dec 2009 and Feb-Mar 2010. A global three dimensional chemical transport Model, GEOS-Chem, (v8-03-01) is employed to generate the SO
  <sub>2</sub> profile for the entire region for the different sampling months which in turn is used to explain the major features in the measured SO
  <sub>2</sub> spectra via correlating with HYSPLIT generated wind back trajectories. The mean SO
  <sub>2</sub> concentrations recorded at the sampling site varied for the different sampling periods (4.3 ppbv in Sep-Oct 2009, 3.4 ppbv in Nov 2009, 3.5 ppbv in Dec 2009, 7.7 ppbv in Feb 2010 and 9.2 ppbv in Mar 2010) which were found to be strongly influenced by long range transport from a source region surrounding 30&#176;N, 75&#176;E—the one projected with the highest SO
  <sub>2</sub> concentration in the GEOS-Chem generated profiles for the region—lying only a few co-ordinates away. A diurnal cycle of SO
  <sub>2</sub> concentration exists throughout the sampling periods, with the greatest day-night changes observed during Feb and Mar 2010, barely detectable during Sep-Oct 2009, and intermediate values for Nov and Dec 2009 which are systematically studied using the time series PBL height and OH radical values from the GEOS-Chem model. During the sampling period in Nov 2009, a plume transport to the sampling site also was detected when a major fire erupted at an oil depot in Jaipur (26.92&#176;N, 75.82&#176;E), located few co-ordinates away. Separate runs of the model, performed to study the long range transport effects, show a drop in the SO
  <sub>2</sub> levels over the sampling region in the absence of transport, throughout the year with Jan to Apr seen to be influenced the lowest by long range transport while Jul and Dec influenced the highest.
 
</p></abstract><kwd-group><kwd>SO&lt;sub&gt;2&lt;/sub&gt;; GEOS-Chem; HYSPLIT; IOCL Oil Fire; PBL Height; OH Radical</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The present day increased fossil fuel and coal usage in Asia, to meet the energy needs of a rapidly emerging economy, has raised concerns about the potential threats of sulphur dioxide (SO<sub>2</sub>) emission from this region [1-5] to the global climate. Further, there exist many recent reports [6-8]—based on direct observations as well as chemical transport modeling studies—of SO<sub>2</sub> from Asia out-flowing to the western Pacific and influencing the climatology of this region.</p><p>In urban areas it is the fossil fuel combustion which contributes to the majority of the anthropogenically produced SO<sub>2</sub>, releasing over 90% of the sulphur content in these fuels to the atmosphere as sulphur dioxide [9,10]. With the known prominent roles of SO<sub>2</sub> in the chemistry of the global troposphere [11-18] and in the sulphur cycle including its potential to interfere with the oxidizing power of the atmosphere, to cause acid rain (e.g., [19, 20]), and to modify the atmospheric radiative forcing pattern via gas-phase chemistry and particle formation, there is an impending need to monitor systematically the emission and transport pattern SO<sub>2</sub> over different region in Asia, with differing environmental conditions.</p><p>As a special case, the monitoring of ambient SO<sub>2</sub> levels coupled with wind trajectory analysis at a location devoid of local emissions and situated at free tropospheric altitudes (seasonal entry into the free troposphere when the planetary boundary layer (PBL) height is low) can help improve the present understanding of the free tropospheric transport patterns of SO<sub>2</sub> from source to remote locations. Keeping this as one of the major goals, ambient SO<sub>2</sub> concentrations were measured at a high altitude site in the semi arid region of western India, Gurushikhar, Mt. Abu (24.6˚N 72.7˚E, 1680 m Above Sea Level (ASL)), during Sep-Oct’09, Nov’09, Dec’09, Feb’10 and Mar’10. In this study a global three dimensional chemical transport model, GEOS-Chem, (v8-03- 01) is also employed, to generate the SO<sub>2</sub> profiles for the entire region for the different sampling months to identify the major SO<sub>2</sub> source regions.</p></sec><sec id="s2"><title>2. Materials and Methods</title><sec id="s2_1"><title>2.1. Site Description: Gurushikhar, Mount Abu (24.6˚N, 72.7˚E, 1680 m ASL)</title><p>The sampling site for this study, Gurushikhar, Mt. Abu (24.6˚N, 72.7˚E, 1680 m ASL) is a high altitude site in the semi-arid region of western India, having a more or less clean atmosphere with minimal local emissions. Gurushikhar is the highest mountain peak in the southern end of Aravali range of mountains in western India with an annual average rainfall of about 600 - 700 mm precipitating only during SW-monsoon. The Gurushikhar site, because of its high elevation, can enter the free tropospheric zone during the winter and post-monsoon months when the Planetary Boundary Layer (PBL) height is low over the region for low incoming solar radiation (insolation) rates, while during summer it gets accommodated within the PBL. This special characteristic makes the Mt. Abu site suitable for the kind of study mentioned above viz., assessing the free tropospheric transport pattern of SO<sub>2</sub> from source to remote locations.</p></sec><sec id="s2_2"><title>2.2. Experimental Setup</title><p>A primary UV fluorescence SO<sub>2</sub> Monitor (Thermo—43i Trace Level Enhanced (TLE)) was employed for the monitoring of ambient SO<sub>2</sub> levels during the different sampling periods. The performances of similar systems are reported elsewhere [21-23]. Calibrations were performed routinely with a standard SO<sub>2</sub> gas (2 ppmv with N<sub>2</sub> balance gas, Spectra, USA) and zero air using a dynamic gas calibrator (Thermo—146i). The integration time for the SO<sub>2</sub> analyser was set at five minutes and the data was periodically downloaded to a computer for analysis.</p></sec><sec id="s2_3"><title>2.3. The GEOS-Chem Model</title><p>The GEOS-Chem global 3-D chemical transport model (v8-03-01; http://acmg.seas.harvard.edu/geos/) employed for this study uses GEOS-5 assimilated meteorological observations—including winds, convective mass fluxes, mixed layer depths, temperature, clouds, precipitation, and surface properties—from NASA Global Modeling and Assimilation Office (GMAO). The GEOS-5 data have a temporal resolution of 6-h (3-hour resolution for surface fields and mixing depths) and horizontal resolution of 0.5˚ latitude &#215; 0.667˚ longitude, with 72 levels in the vertical extending from surface to approximately 0.01 hPa. The GEOS-Chem model (v8-03-01) simulates the Ozone-NOx-VOC-aerosol chemistry [24,25] and elaborative simulation evaluations are provided in [25-28].</p><p>In the present work, simulations were carried out for Jan 2009-Apr 2010 at 4˚ &#215; 5˚ resolutions with the runs initialized on 1<sup>st</sup> Jan 2009 using GEOS-Chem fields generated by a one year spin-up simulation at 4˚ &#215; 5˚ resolution. The emission inventories EMEP [<xref ref-type="bibr" rid="scirp.21141-ref29">29</xref>], BRAVO [<xref ref-type="bibr" rid="scirp.21141-ref30">30</xref>], EDGAR [<xref ref-type="bibr" rid="scirp.21141-ref31">31</xref>], Streets inventory [<xref ref-type="bibr" rid="scirp.21141-ref32">32</xref>], CAC (http://www.ec.gc.ca/pdb/cac/cac_home_e.cfm), and EPA/ NEI05 were included in the model runs by setting all of them to “True” in the standard input.geos file. Unless mentioned otherwise, the parameters set in the input.geos file for the model runs are identical to the ones in the standard input.geos file, distributed with the GEOSChem codes (v8-03-01), with the simulation start/end dates and the ND49 diagnostics—to generate the time series data—modified.</p></sec><sec id="s2_4"><title>2.4. GAMAP Visualization Toolkit</title><p>The GEOS-Chem Model outputs during different runs are plotted using the Global Atmospheric Model Analysis Package (GAMAP) (Version 2.15; http://acmg.seas.harvard.edu/gamap/) which is a selfcontained, consistent, and user-friendly software package written in IDL (Interactive Data Language). The package can read outputs from chemical tracer models (CTM’s) and visualize it via line plots, 2D plots, 2D animations, or 3D iso-contour surface plots.</p></sec><sec id="s2_5"><title>2.5. The HYSPLIT Trajectory Model</title><p>The wind back trajectory analysis for this work employs the HYSPLIT (HYbrid Single-Particle Lagrangian Integrated Trajectory) model, developed by NOAA ARL [33, 34] which uses a modelled vertical velocity scheme and can depict the vertical motion of the relevant air parcel. The 7-day back trajectories plotted for altitudes 50, 500 and 3000 m above the ground level (AGL) are used to identify the origin and the transect of the air parcels reaching the sampling site.</p></sec></sec><sec id="s3"><title>3. Results and Discussions</title><p>Ambient SO<sub>2</sub> measurements were performed during Sep-Oct’09, Nov’09, Dec’09, Feb’10 and Mar’10. The GEOS-Chem chemical transport model (v8-03-01) is used to generate the SO<sub>2 </sub>profile for the entire region surrounding the sampling site for the different sampling months to identify the major source regions. Now by analysing the back trajectories generated using HYSPLIT for the sampling periods—performed at different altitudes to investigate the source and the transport pathways of the air parcel, prior to its arrival to the sampling site— and correlating them with the GEOS-Chem generated SO<sub>2</sub> plots showing the source regions, the major features in the measured spectra are explained.</p><sec id="s3_1"><title>3.1. SO<sub>2</sub> Plots for Different Sampling Months from GEOS-Chem</title><p>The GEOS-Chem model (v8-03-01) runs were performed for Sep 2009 to Apr 2010, with the inventories discussed previously, to find out the major SO<sub>2 </sub>source regions in the vicinity of the sampling site. The output averaged for the first half of each sampling month was then plotted using GAMAP and is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. The plots showed a region surrounding 30˚N, 75˚E with the highest SO<sub>2</sub> concentration in the vicinity. With the sampling site, Gurushikhar located only a few co-ordinates away, a long range transport to the sampling site from this high SO<sub>2</sub> region (surrounding 30˚N, 75˚E) is feasible under favourable wind conditions, including free tropospheric transport, for the high elevation of the sampling site.</p></sec><sec id="s3_2"><title>3.2. Time Series SO<sub>2</sub> Profile from GEOS-Chem</title><p>The time series SO<sub>2</sub> profiles for the nearest co—ordinates to the Mt. Abu site possible with the model (26˚N, 75˚E), for the ground level (0 - 0.3 km AGL), also was generated during the GEOS-Chem runs by turning ON the ND49 diagnostics in the input.geos file. This model generated SO<sub>2</sub> time series were then plotted for the sampling periods to make comparison with the experimental results. The data is also used to study the diurnal variabilities and its dependence on different atmospheric parameters such as Planetary Boundary Layer (PBL) height and OH radical concentrations.</p></sec><sec id="s3_3"><title>3.3. SO<sub>2</sub> Measurements at Mt. Abu</title><p>The measured SO<sub>2 </sub>time series profiles had distinctly different features during different sampling periods span over Sep-Oct’09, Nov’09, Dec’09, Feb’10, and Mar’10, and were found to be dictated by the transport patterns and climatology prevailed in the respective months. The major features in the SO<sub>2</sub> spectra for different sampling months and its correlation with various atmospheric parameters are discussed below. A comparison of the experimental SO<sub>2</sub> time series data with GEOS-Chem model generated one’s also is made for the different sampling periods.</p><sec id="s3_3_1"><title>3.3.1. Major Features: September-October 2009</title><p>During Sep-Oct 2009, when the winds were mostly from the South-West direction (as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>), with the 7 day back trajectories originating mostly from Arabian sea where SO<sub>2</sub> concentration is a minimum in the GEOS-Chem generated plots (<xref ref-type="fig" rid="fig1">Figure 1</xref>), the SO<sub>2</sub> time series showed a near constant continuum with very few spikes (<xref ref-type="fig" rid="fig3">Figure 3</xref>). The mean SO<sub>2</sub> concentration for this sampling period is 4.3 ppbv. The spikes in the spectrum seen on the 2<sup>nd</sup> and 3<sup>rd</sup> of Oct 2009 correlate very well with a change in wind trajectory pattern (<xref ref-type="fig" rid="fig2">Figure 2</xref>) when an air parcel from the high SO<sub>2</sub> source region surrounding 30˚N, 75˚E reached the sampling site.</p><p>The SO<sub>2</sub> profile had minimal diurnal variability (<xref ref-type="fig" rid="fig4">Figure 4</xref>) which is explained as due to the low PBL heights as well as high atmospheric OH radical concentration during these months over the sampling region (as will be discussed in section 3.4). The PBL height for this month over the sampling region is amongst the lowest of the year and the sampling site Gurushikhar is more or less within the free troposphere, nullifying the PBL height variation effects and hence the observed minimum diurnal variation for the SO<sub>2</sub> levels. Also the high OH radical concentration due to a higher water vapour content and incoming solar radiation (insolation) rate—the two parameters necessary for the production of OH radical in the atmosphere via photo chemistry—possibly makes sure sufficient oxidation of SO<sub>2</sub> even when the concentration is high to minimise the diurnal variation in the SO<sub>2</sub> oxidation efficiency. The very negligible contribution from long range transport during this sampling period—as almost all the back trajectories originated in the Arabian sea where the SO<sub>2</sub> concentration is a minimum in the GEOS-Chem generated plots, resulting in non dependence of the measured SO<sub>2</sub> levels on possible diurnal variability of wind pattern—can be another possible reason for the minimal diurnal variation in the SO<sub>2</sub> levels. The GEOS-Chem generated plots for the sampling period (shown in <xref ref-type="fig" rid="fig3">Figure 3</xref>) however showed some diurnal variation in SO<sub>2</sub> levels<sub> </sub>which may be attributed to the difference in the altitudes of the sampling site and the co-ordinates used in the model calculations (due to the constraint that at 4˚ &#215; 5˚ resolution, the nearest possible grid point available with the model was 26˚N, 75˚E, where the terrain height may be different).</p></sec><sec id="s3_3_2"><title>3.3.2. Major Features: November 2009</title><p>During Nov 2009 when the wind back trajectories (<xref ref-type="fig" rid="fig5">Figure 5</xref>) extended mostly to North-East, the SO<sub>2</sub> profile showed short duration spikes superimposed on the back ground values. The mean of the SO<sub>2</sub> concentrations recorded (<xref ref-type="fig" rid="fig6">Figure 6</xref>) for this sampling period is 3.4 ppbv. For this month, the wind back trajectories showed occasional long range transport from the high SO<sub>2</sub> region (surrounding 30˚N, 75˚E) along with those from the low SO<sub>2</sub> regions (e.g., regions lying in the North-East of the sampling site) in the GEOS-Chem generated plots.</p><p>Because of the wind pattern specific for this month, only short duration random air parcels arrive at the sampling site at Mt. Abu from the high SO<sub>2</sub> region surrounding 30˚N, 75˚E—as shown by the back trajectory analysis (<xref ref-type="fig" rid="fig5">Figure 5</xref>)—explaining the observed short duration spikes in the SO<sub>2</sub> profile.</p><p>A diurnal pattern is observed (<xref ref-type="fig" rid="fig7">Figure 7</xref>) in the measured SO<sub>2</sub> spectra, which may be attributed to the combined effects put in by the insolation dependent variation in OH radical levels and PBL height as well as the diurnal variation in the wind pattern carrying polluted air parcel to the sampling site from long distances.</p><p>The spikes as well as the concentration levels decreased from 12<sup>th</sup> Nov 2009, as is the case with the GEOS—Chem generated SO<sub>2</sub> time series shown also in <xref ref-type="fig" rid="fig6">Figure 6</xref>—which has association with the wind trajectory change from high SO<sub>2</sub> regions surrounding 30˚N, 75˚E to relatively low SO<sub>2</sub> regions lying North-East of the sampling site.</p><p>Detection of SO<sub>2</sub> Plume during a Major Oil Fire Incident A major fire had erupted at an oil depot of Indian Oil Corporation Limited (IOCL) at Jaipur, Rajasthan (26.92˚N, 75.82˚E) on 29th October 2009 at around 7:30 PM (IST). Approximately 60,000 kilolitres of oil were burned in this fire incident which lasted for almost a week till 6<sup>th</sup> Nov 2009, resulting in a significant rise in the air pollution levels over the region.</p><p>Studying the possible long range transport of SO<sub>2</sub> containing plumes to the high altitude sampling site at Mt. Abu from the disaster site, located only a few co-ordinates away was another objective of the sampling for this month. Three short duration high intensity spikes were recorded (<xref ref-type="fig" rid="fig8">Figure 8</xref>) each on 2<sup>nd</sup> (around 0125 hrs) 4<sup>th</sup><sup></sup></p><p>(around 2230 hrs) and 5<sup>th</sup> (around 2120 hrs) of Nov 2009 which was found associated with this major oil fire from back trajectory analysis—which showed air parcels reaching the sampling site from the disaster site (26.92˚N, 75.82˚E), during the above days (<xref ref-type="fig" rid="fig5">Figure 5</xref>).</p></sec><sec id="s3_3_3"><title>3.3.3. Major Features: December 2009</title><p>During Dec 2009, the winds—as seen in the back trajectory plots shown in <xref ref-type="fig" rid="fig9">Figure 9</xref>—were mostly from North and North-East of the sampling site and occasionally carried high SO<sub>2</sub> from the source region surrounding</p><p>30˚N, 75˚E to produce many spikes in the measured spectra (<xref ref-type="fig" rid="fig1">Figure 1</xref>0). The mean SO<sub>2</sub> concentration for this sampling period is 3.5 ppbv. The observed shift in concentration from higher during the initial sampling days to lower towards the end of the sampling period is explained by the observed shift in wind back trajectories from high SO<sub>2</sub> source regions (surrounding 30˚N, 75˚E) to low SO<sub>2</sub> regions (lying in the eastern side of the sampling site).</p><p>A diurnal pattern is observed during this sampling period (<xref ref-type="fig" rid="fig1">Figure 1</xref>1) with the minimum around 1000 hrs IST, which is explained as due to the combined effects put in by the diurnal variability in the PBL height and OH radical levels. A possible periodic variation in the wind pattern carrying polluted air to the sampling site from long distances also could contribute to the observed diurnal variability.</p><p>The GEOS-Chem generated SO<sub>2</sub> time series (seen in <xref ref-type="fig" rid="fig1">Figure 1</xref>0) also reproduced the trend observed for this month including the hump in the concentration values between 10<sup>th</sup> and 15<sup>th</sup><sub> </sub>of this month.<sub> </sub>&#160;</p><p>The GEOS-Chem based calculations showed (as will be discussed in section 3.5) that it is in the month of December every year, the contribution of long range transported</p><p>SO<sub>2</sub> over Mt. Abu attains its peak.</p></sec><sec id="s3_3_4"><title>3.3.4. Major Features: February 2010</title><p>The sampling period in Feb 2010 was characterized by high SO<sub>2</sub> concentrations (<xref ref-type="fig" rid="fig1">Figure 1</xref>2) and clear diurnal variability. The trend observed is high SO<sub>2</sub> concentrations during the initial days of the sampling followed by a dip and the concentration again rose towards the end of the sampling period. This is well explained by comparing backward wind trajectories (<xref ref-type="fig" rid="fig1">Figure 1</xref>3) with GEOSChem generated SO<sub>2</sub> profiles (<xref ref-type="fig" rid="fig1">Figure 1</xref>). During the initial days of sampling all three components of the wind trajectories (plotted for altitudes 50 m, 500 m and 3000 m AGL) extended to the high SO<sub>2</sub> regions (surrounding 30˚N, 75˚E) in the GEOS-Chem generated plots. But during the time when dip in SO<sub>2</sub> values were observed, the only trajectory at 50 m AGL carried polluted air parcel from the high SO<sub>2</sub> region (30˚N, 75˚E) whereas the higher altitude trajectories (500 m and 3000 m AGL) came from the West Asian region where no major hot spots of SO<sub>2</sub> prevail in the GEOS-Chem generated plots. The fact that the Mt. Abu site enters the free tropospheric</p><p>zone during the winter time explains the non-influence of pollutant carrying winds in the lower layers of the atmosphere on the ambient SO<sub>2</sub> levels for this month.</p><p>The mean SO<sub>2</sub> concentration for this sampling period is 7.7 ppbv. The strong diurnal pattern observed for this sampling period (<xref ref-type="fig" rid="fig1">Figure 1</xref>4) is attributed to the large diurnal PBL height variation as seen in the GEOS-Chem generated PBL values for this month (discussed in section 3.4).</p><p>For this sampling period, the GEOS-Chem generated SO<sub>2</sub> profile (seen in <xref ref-type="fig" rid="fig1">Figure 1</xref>2) deviated in magnitude from experimental observations, suggesting the prominent roles of local sources as well as parameters not incorporated in the model in controlling the ambient SO<sub>2</sub> levels during this month.</p></sec><sec id="s3_3_5"><title>3.3.5. Major Features: March 2010</title><p>During the sampling period in Mar 2010 also, relatively high SO<sub>2</sub> concentrations were recorded (<xref ref-type="fig" rid="fig1">Figure 1</xref>5). Occasional long range transport of SO<sub>2 </sub>to the sampling site is indicated by the back trajectory analysis which showed (<xref ref-type="fig" rid="fig1">Figure 1</xref>6) arrival of air parcel from both the high SO<sub>2</sub> region surrounding 30˚N, 75˚E and low SO<sub>2</sub> regions such as the Arabian sea and West Asia. The mean SO<sub>2</sub> concentration for this sampling period is 9.2 ppbv. A well defined diurnal variability (<xref ref-type="fig" rid="fig1">Figure 1</xref>7) is a characteristic feature for this period, which is attributed to the very large diurnal variability in the PBL heights as seen in the GEOS-Chem generated time series PBL values (discussed in section 3.4). The diurnal pattern resembled the one observed in Feb 2010 with the whole graph shifted by ~3 hrs to the left along the time axis. This shift is explained as due to the increase in the insolation rate from February to March, to cause the PBL rise to take place in the morning hours itself compared to February. Possible diurnal variation in the wind pattern carrying polluted air to the sampling site also could contribute to the observed diurnal variability.</p></sec></sec><sec id="s3_4"><title>3.4. Explaining the Diurnal Variations in SO<sub>2</sub> Conc. Using GEOS-Chem Model</title><p>The diurnal variation patterns observed in the measured SO<sub>2</sub> profiles over Mt. Abu as well as those in the GEOS-Chem generated SO<sub>2</sub> time series data for the co-ordinates 26˚N, 75˚E (shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>8)—the nearest grid point available with the model at 4˚ &#215; 5˚ resolution runs—for different months of the year 2009 can be well explained based on the time series PBL height and OH radical values generated using the GEOSChem model. The diurnal pattern seen in the time series PBL height and OH radical values generated using GEOS-Chem model at 26˚N, 75˚E are shown in Figures 19 and 20 respectively.</p><p>The GEOS-Chem generated SO<sub>2</sub> profiles (<xref ref-type="fig" rid="fig1">Figure 1</xref>8) over the co-ordinates (26˚N, 75˚E) showed that during Jan, Feb, Mar, Apr and May there is a sudden dip in the SO<sub>2</sub> concentration between 0600 hrs and 0800 hrs IST. This may be explained on the basis of the diurnal variation pattern seen in the time series PBL height (<xref ref-type="fig" rid="fig1">Figure 1</xref>9) and the OH radical profiles (<xref ref-type="fig" rid="fig2">Figure 2</xref>0). During the above months (Jan, Feb, Mar, Apr and May), the PBL height remain below 100 m during night time till 0600 hrs IST followed by a 10 - 20 fold increase in PBL height by 0800 hrs, which continue to rise till 1400 hrs, remain more or less constant till 1600 hrs, followed by a decrease in its value to reach the minimum by 2000 hrs IST</p><p>and continue to remain at this minima till 0600 hrs IST. The spread in the PBL values increases from Jan to Jun with maximum spread for the month of Jun, followed by a decrease. The minimum diurnal variability in the PBL height was observed during Jul to Sep.</p><p>Unlike the diurnal variation pattern seen in the PBL heights, where during Jan-May there is a transient increase in the values between 0600 hrs and 0800 hrs, the OH variation during the months Jan-Mar are at a relatively slow pace. The sudden manifold increase in OH levels between 0600 hrs and 0800 hrs IST takes place during May, Jun and Jul when the relative humidity (RH) in the atmosphere is sufficiently high over this region. The GEOS-Chem generated OH radical profiles showed highest concentrations during the months, Jun, Jul, Aug, and Sep followed by a decrease till Dec and then the accenting mode for OH radical levels in the atmosphere start from Jan. The minimum OH concentrations were recorded during Dec when both the insolation as well the atmospheric water content—the two major parameters</p><p>controlling the OH formation in the atmosphere—over this region are on the lower side.</p><p>As can be seen from the above discussion, the PBL height and the OH radical maxima doesn’t overlap but falls on different months to give the relative influence of these parameters to the SO<sub>2 </sub>oxidation efficiency a strong seasonality and is discussed in detail below.</p><sec id="s3_4_1"><title>Relative Influence of PBL Height and OH Radical on the SO<sub>2</sub> Levels</title><p>It is evident from the above discussion that interplay between seasonal PBL height and the OH concentration variability can modulate the diurnal pattern of the SO<sub>2</sub> concentrations. The fact that the rates at which these two parameters vary—both seasonal as well as diurnal are different, leads to a seasonal dependence for the diurnal variability of SO<sub>2</sub> levels.</p><p>For example, during the months Feb, Mar and Apr the effect put in by the sudden change in PBL height in the morning hours leads to the sudden decrease in the SO<sub>2</sub> values between 0600 hrs and 0800 hrs IST, where the OH radical variation is not sufficient to cause much of a change in the SO<sub>2</sub> variability pattern. Whereas, the sudden changes in the OH concentrations occurring during Jun-Sep is the major player controlling the SO<sub>2</sub> diurnal pattern during these months.</p><p>To visualize the above mentioned relative contributions clearly, 3-dimensional graphs were plotted with the GEOS-Chem generated OH concentration along X-axis, PBL height along Y-axis and the SO<sub>2</sub> concentration along the Z-axis for different months of the year 2009 (<xref ref-type="fig" rid="fig2">Figure 2</xref>1). During all the months, the SO<sub>2</sub> concentration decreased with PBL height increase as well as the OH radical concentration, but with varying magnitude. During Jan 2009, a PBL height variation from the lowest of the month to ~600 m takes the SO<sub>2</sub> value to a constant for all the OH values lying above ~2e<sup>+6</sup>. Similarly even with the highest OH concentration available with the month (~7e<sup>+6</sup>), the SO<sub>2</sub> concentration is not pulled down for PBL values below ~400 m. The random peaks observed in the SO<sub>2</sub> values for Jan has probably got to do with the influence of other parameters randomly controlling the SO<sub>2</sub> concentration over the region, and the transient long range transport is one such possible parameter.</p><p>During Feb 2009, for all values of OH concentration above ~4e<sup>+6</sup>, an increase in PBL height above ~700 m is sufficient to pull the SO<sub>2</sub> values to the minimum for that month. Below an OH level of ~4e<sup>+6</sup>, a combined effect of OH radical as well as PBL height variation is found to influence the SO<sub>2</sub> concentration. Similarly for OH concentrations below ~2e<sup>+6</sup>, even with the highest PBL height recorded for the month, the SO<sub>2</sub> concentration remained high. Compared to Jan 2009, the graph plotted for Feb 2009 has fewer uncorrelated humps for the SO<sub>2</sub> values, indicating the possible shift in the transport pattern.</p><p>The month, Mar 2009 saw a clear interplay between the PBL and OH values in deciding the diurnal pattern of ambient SO<sub>2 </sub>over this region. This is one such month where both PBL height and OH levels are significant. As a result, for OH values between ~1e<sup>+6</sup> and 4e<sup>+6</sup>, the increase in PBL height from the minimum of the month to ~1500 m, reduces the SO<sub>2</sub> concentration at a steady rate to reach the minimum of the month. For all OH values above ~4e<sup>+6</sup>, the SO<sub>2</sub> concentration drops to the minimum, for PBL heights above ~700 m. Interference from</p><p>other players such as transient long range transport is a minimum in this month as evidenced by the absence of any random humps in the SO<sub>2</sub> values.</p><p>The behaviour observed in Apr, May and Jun were somewhat similar to the one in Mar. The scenario started changing from Jul and became clearly visible by Aug, when the atmospheric OH concentration over the sampling region reached its peak values of the year. In Aug, the OH concentration was sufficient to pull down the SO<sub>2</sub> values to the minimum of the month for all PBL values above ~150 m, showing the major role the OH radical play in the diurnal pattern of SO<sub>2</sub> for this month. The non observation of any humps in the SO<sub>2</sub> values correlate very much with the minimum transient transport expected</p><p>for these months (as will be discussed in section 3.5), underlining the argument that the third major player controlling the ambient SO<sub>2</sub> levels over this region is the long range transport.</p><p>During Sep 2009 and Oct 2009, when the boundary layer heights are on the lower side, the OH radical and the transport are the only major players, resulting in many humps in the SO<sub>2</sub> values in addition to the contribution from heterogeneous phase oxidation on water droplets. During Nov and Dec, an overlap of low PBL heights with low OH concentration along with higher transport events (as discussed below) keep the SO<sub>2</sub> values on the higher side during most of the time.</p></sec></sec><sec id="s3_5"><title>3.5. Effect of Long Range Transport on the Ambient SO<sub>2</sub> Levels over Mt. Abu</title><p>In order to understand the possible effects of long range</p><p>transport on the measured SO<sub>2</sub> values over Mt. Abu, GEOS-Chem model runs were performed with all the parameters keeping identical to the previous runs but now with the Transport mechanism turned OFF. Now the percentage difference in the SO<sub>2</sub> concentration at 26˚N, 75˚E—the nearest grid point co-ordinates available with the GEOS-Chem model at 4˚ &#215; 5˚ resolutions in the absence of transport is calculated as follows:</p><p>Percent Difference in the SO<sub>2</sub> Concentration in the Absence of Transport = (WithOutTransport – FullInput) &#215; 100/FullInput</p><p>where:</p><p>WithOutTransport → SO<sub>2</sub> concentration recorded when the Transport turned OFF</p><p>FullInput → SO<sub>2</sub> concentration recorded when Transport also was kept ON</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>2 shows the percent difference in the SO<sub>2</sub> concentration in the absence of the transport mechanism over the co-ordinates (26˚N, 75˚E). Once the transport is turned OFF, sources located within the 4˚ &#215; 5˚ grid containing</p></sec></sec></body><back><ref-list><title>References</title><ref id="scirp.21141-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">J. Crawford, J. Olson, D. Davis, et al., “Clouds and Trace Gas Distributions during TRACE-P,” Journal of Geophysical Research, Vol. 108, No. D21, 2003, pp. 8818-8830. doi:10.1029/2002JD003177</mixed-citation></ref><ref id="scirp.21141-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">T. Larssen, E. Lydersen, D. Tang, et al., “Acid Rain in China,” Environmental Science &amp; Technology, Vol. 40, No. D2, 2006, pp. 418-425.</mixed-citation></ref><ref id="scirp.21141-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">T. Ohara, H. Akimoto, J. Kurokawa, et al., “An Asian Emission Inventory of Anthropogenic Emission Sources for the Period 1980-2020,” Atmospheric Chemistry and Physics Discuss, Vol. 7, 2007, pp. 4419-4444.  
doi:10.5194/acpd-7-6843-2007</mixed-citation></ref><ref id="scirp.21141-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">D. G. Streets, T. C. Bond, G. R. Carmichael, et al., “An Inventory of Gaseous and Primary Aerosol Emissions in Asia in the Year 2000,” Journal of Geophysical Research, Vol. 108, No. D21, 2003, pp. 8809-8831.  
doi:10.1029/2002JD003093</mixed-citation></ref><ref id="scirp.21141-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">D. G. Streets and S. T. Waldhoff, “Present and Future Emissions of Air Pollutants in China: SO2, NOX, and CO,” Atmospheric Environment, Vol. 34, No. 3, 2000, pp. 363-374. doi:10.1016/S1352-2310(99)00167-3</mixed-citation></ref><ref id="scirp.21141-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">T. Francis, “Effect of Asian Dust Storms on the Ambient SO2 Concentration over North-East India: A Case Study,” Journal of Environmental Protection, Vol. 2, No. 6, 2011, pp. 778-795. doi:10.4236/jep.2011.26090</mixed-citation></ref><ref id="scirp.21141-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">D. C. Thornton, A. R. Bandy, B. W. Blomquist, et al., “Sulfur Dioxide as a Source of Condensation Nuclei in the Upper Troposphere of the Pacific Ocean,” Journal of Geophysical Research, Vol. 101, No. D1, 1996, pp. 1883 -1890. doi:10.1029/95JD02273</mixed-citation></ref><ref id="scirp.21141-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">D. C. Thornton, A. R. Bandy, B. W. Blomquist, et al., “Transport of Sulfur Dioxide from the Asian Pacific Rim to the North Pacific Troposphere,” Journal of Geophysical Research, Vol. 102, No. D23, 1997, pp. 28489-28499.  
doi:10.1029/97JD01818</mixed-citation></ref><ref id="scirp.21141-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">C. M. Benkovitz, M. T. Scholtz, J. Pacyna, et al., “Global gridded Inventories of Anthropogenic Emissions of Sulfur and Nitrogen,” Journal of Geophysical Research, Vol. 101, No. D22, 1996, pp. 29239-29253.  
doi:10.1029/96JD00126</mixed-citation></ref><ref id="scirp.21141-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">P. A. Spiro, D. J. Jacob and J. A. Logan, “Global Inventory of Sulfur Emissions with 1? × 1? Resolution,” Journal of Geophysical Research, Vol. 97, No. D5, 1992, pp. 6023-6036. doi:10.1029/91JD03139</mixed-citation></ref><ref id="scirp.21141-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">M. O. Andreae, “Ocean-Atmosphere Interactions in the Global Biogeochemical Sulfur Cycle,” Marine Chemistry, Vol. 30, 1990, pp. 1-29.  
doi:10.1016/0304-4203(90)90059-L</mixed-citation></ref><ref id="scirp.21141-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">M. O. Andreae and W. A. Jaeschke, “Exchange of Sulfur between Biosphere and Atmosphere over Temperate and Tropical Regions,” In: R. W. Howarth and J. W. B. Stewart, Eds., Sulfur Cycling in Terrestrial Systems and Wetlands, John Wiley, New York, 1992.</mixed-citation></ref><ref id="scirp.21141-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">T. S. Bates, J. D. Cline, R. H. Gammon, et al., “Regional and Seasonal Variations in the Flux of Oceanic Dimethyl Sulfide to the Atmosphere,” Journal of Geophysical Research, Vol. 92, No. C3, 1987, pp. 2930-2938.  
doi:10.1029/JC092iC03p02930</mixed-citation></ref><ref id="scirp.21141-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">C. F. Cullis and M. M. Hirschler, “Atmospheric Sulphur: Natural and Man-Made Sources,” Atmospheric Environment, Vol. 14, No. 11, 1980, pp. 1263-1278.</mixed-citation></ref><ref id="scirp.21141-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">J. R. Freney, M. V. Ivanov and H. Rhodhe, “The Sulphur Cycle, in the Major Biogeochemical Cycles and Their Interactions, SCOPE 24,” John Wiley, New York, 1983.</mixed-citation></ref><ref id="scirp.21141-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">H. Rodhe and I. Isaksen, “Global Distribution of Sulfur Compounds in the Troposphere Estimated in a Height/Latitude Transport Model,” Journal of Geophysical Research, Vol. 85, No. C12, 1980, pp. 7401-7409.  
doi:10.1029/JC085iC12p07401</mixed-citation></ref><ref id="scirp.21141-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">E. S. Saltzman and D. J. Cooper, “Shipboard Measurements of Atmospheric Dimethyl Sulfide and Hydrogen Sulfide in the Caribbean and Gulf of Mexico,” Journal of Atmospheric Chemistry, Vol. 7, No. 2, 1988, pp. 191-209.  
doi:10.1007/BF00048046</mixed-citation></ref><ref id="scirp.21141-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">O. B. Toon, et al., “The Sulfur Cycle in the Marine Atmosphere,” Journal of Geophysical Research, Vol. 92, No. D1, 1987, pp. 943-963.  
doi:10.1029/JD092iD01p00943</mixed-citation></ref><ref id="scirp.21141-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">L. Gerhardsson, A. Oskarsson and S. Skerfving, “Acid Precipitation—Effects on Trace Elements and Human Health,” The Science of the Total Environment, Vol. 153, No. 3, 1994, pp. 237-245.  
doi:10.1016/0048-9697(94)90203-8</mixed-citation></ref><ref id="scirp.21141-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">D. W. Schindler, “Effects of Acid Rain on Freshwater Ecosystems,” Science, Vol. 239, No. 4836, 1988, pp. 149-157.</mixed-citation></ref><ref id="scirp.21141-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Y. Igarashi, Y. Sawa, K. Yoshioka, et al., “Monitoring the SO2 Concentration at the Summit of Mt. Fuji and a Comparison with Other Trace Gases during Winter,” Journal of Geophysical Research, Vol. 109, No. D17, 2004, Article ID: D17304. doi:10.1029/2003JD004428</mixed-citation></ref><ref id="scirp.21141-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">W. T. Luke, “Evaluation of a Commercial Pulsed Fluorescence Detector for the Measurement of Low-Level SO2 Concentrations during the Gas-Phase Sulfur Intercomparison Experiment,” Journal of Geophysical Research, Vol. 102, No. D13, 1997, pp. 16255-16265.  
doi:10.1029/96JD03347</mixed-citation></ref><ref id="scirp.21141-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">M. Luria, J. F. Boatman, J. Harris, et al., “Atmospheric Sulfur Dioxide at Mauna Loa, Hawaii,” Journal of Geophysical Research, Vol. 97, No. D5, 1992, pp. 6011-6022.  
doi:10.1029/91JD03126</mixed-citation></ref><ref id="scirp.21141-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">I. Bey, D. J. Jacob, R. M. Yantosca, et al., “Global Modeling of Tropospheric Chemistry with Assimilated Meteorology: Model Description and Evaluation,” Journal of Geophysical Research, Vol. 106, No. D19, 2001, pp. 23073-23095. doi:10.1029/2001JD000807</mixed-citation></ref><ref id="scirp.21141-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">R. J. Park, D. J. Jacob, B. D. Field, et al., “Natural and Transboundary Pollution Influences on Sulfate-Nitrate-Ammonium Aerosols in the United States: Implications for Policy,” Journal of Geophysical Research, Vol. 109, No. D15, 2004, Article ID: D15204.  
doi:10.1029/2003JD004473</mixed-citation></ref><ref id="scirp.21141-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">J. Wang, D. J. Jacob and S. T. Martin, “Sensitivity of Sulfate Direct Climate Forcing to the Hysteresis of Particle Phase Transitions,” Journal of Geophysical Research, Vol. 113, No. D11, 2008, Article ID: D11207.  
doi:10.1029/2007JD009368</mixed-citation></ref><ref id="scirp.21141-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">T. Duncan Fairlie, D. J. Jacob and R. J. Park, “The Impact of Transpacific Transport of Mineral Dust in the United States,” Atmospheric Environment, Vol. 41, No. 6, 2007, pp. 1251-1266. doi:10.1016/j.atmosenv.2006.09.048</mixed-citation></ref><ref id="scirp.21141-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">D. Chen, Y. Wang, M. B. McElroy, et al., “Regional CO Pollution and Export in China Simulated by the High-Resolution Nested-Grid GEOS-Chem Model,” Atmospheric Chemistry and Physics, Vol. 9, No. 11, 2009, pp. 3825-3839. doi:10.5194/acp-9-3825-2009</mixed-citation></ref><ref id="scirp.21141-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">V. Vestreng and H. Klein, “Emission Data Reported to UNECE/EMEP. Quality Assurance 769 and Trend Analysis and Presentation of WebDab,” MSC-W Status Report, Norwegian Meteorological Institute, 2002.</mixed-citation></ref><ref id="scirp.21141-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">H. Kuhns, M. Green and V. Etyemezian, “Big Bend Regional Aerosol and Visibility Observational (BRAVO) Study Emissions Inventory,” Desert Research Institute, 2003.</mixed-citation></ref><ref id="scirp.21141-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">J. G. J. Olivier and J. J. M. Berdowski, “Global Emissions Sources and Sinks, in the Climate System,” The Netherlands A.A. Balkema Publishers/Swets &amp; Zeitlinger Publishers, 2001.</mixed-citation></ref><ref id="scirp.21141-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">D. G. Streets, Q. Zhang, L. Wang, et al., “Revisiting China’s CO Emissions after the Transport and Chemical Evolution over the Pacific (TRACE-P) Mission: Synthesis of Inventories, Atmospheric Modeling, and Observations,” Journal of Geophysical Research, Vol. 111, No. D14, 2006, Article ID: D14306.  
doi:10.1029/2006JD007118</mixed-citation></ref><ref id="scirp.21141-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">R. R. Draxler and G. D. Rolph, “HYSPLIT—HYbrid Single-Particle Lagrangian Integrated Trajectory Model,” 2012. http://www.arl.noaa.gov/ready/hysplit4.html</mixed-citation></ref><ref id="scirp.21141-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">G. D. Rolph, “Real-Time Environmental Applications and Display System (READY),” 2003.</mixed-citation></ref></ref-list></back></article>