Dispersion Modeling of the Transport and Dispersion of SO 2 Pollutants Emitted from a Power Plant in Tong Liang

Final Report for the Project Dispersion Modeling of the Transport and Dispersion of SO 2 Pollutants Emitted from a Power Plant in Tong Liang Contractors: Dr. Darko Koracin, Dr. John Watson, 1. Executive summary We have performed a dispersion modeling study of the transport and dispersion of SO 2 pollutants emitted from the Tong Liang power plant. Results from the dispersion modeling study will be used to support an ongoing epidemiological study in this area. The dispersion modeling study was conducted using the EPA regulatory dispersion model Industrial Source Complex 3 Short Term (ISC3ST). The simulation period covered two years, 2002 and 2003, for which required input meteorological data sets were available. The modeling results consist of annual average, monthly maximum, daily maximum, 3-hr maximum, and all monthly averages during the power plant operation of SO 2 concentrations for 2002 and 2003 at each of 170 specified receptors. Additionally, we defined 216 polar receptors covering six distances (100, 500, 1000, 2000, 4000, and 7000 m) and each of them having 36 azimuth locations. These 386 receptors provide details on the spatial distribution of concentrations that helps to understand overall structure and evolution of the pollutant plume in the area. These results represent a basis for statistical and graphical analysis. 2. Project motivation A comprehensive study was conducted by Perera et al. (2005) to assess the health effects of industrial-origin pollution on newborns in Tong Liang, China. The main objective of the study was to investigate the extent to which reductions in coal-burning emissions might result in the improvements to the health of children. Chow et al. (2005) provided further insight into the pollution impact by analyzing measurements of speciated PM2.5, carbon, and particle-bound polycyclic aromatic hydrocarbon (PAH) at three community exposure sites. They found high concentrations compared to U.S. cities and that the highest concentrations occurred in winter. These high levels result from a mixture of contribution from several sources, including the Tong Liang power station, other industries (e.g., strontium), distant power stations, domestic coal use, burning and cooking, and vehicle exhaust. 1

3. Background Tong Liang County has approximately 800,000 inhabitants and is located about 60 km north of Chongqing City in the upper Yangtze River Valley. Most of the population is located in the downtown area with a radius of only about 3 km. The local generator of electricity (46,000 kw) is a coal-fired power plant, which is located only one kilometer south of the city center. The power plant is operational from mid-november to mid-may when the main need for electricity is present. Figure 1 shows topography surrounding of the power plant and Figure 2 shows the plant in relation to the city and monitoring stations. Fig. 1. Topographical setup of the Tong Liang power plant. 2

Fig. 2. Geographical setup of the Tong Liang County with indicated power plant and the monitoring stations (A, B, and C). During its annual operation, the power plant consumes 5,000 tons of coal per month. Since the coal contains 3-5% of sulfur by weight, it is important to investigate the extent to which SO2 that results from its combustion can contribute to the total pollution impact in the area. Other directly emitted pollutant contributions to ambient levels can be estimated if their ratios to SO2 in the stack emissions are known. 3

4. Methodology 4.1 Dispersion model Industrial Source Complex Short Term 3 (ISCST3) model ISCST3 is a steady-state, Gaussian plume dispersion model designed for assessing pollutant concentrations from elevated point sources. According to the EPA Guideline on Air Quality Models (U.S. EPA, 1995), the ISCST3 model (classified as an A group preferred model) is recommended for industrial source complexes, rural or urban areas, flat or rolling terrain, transport distances less than 50 km, 1 hr to annual averaging time, and continuous pollutant emissions (Koracin et al., 2000). ISCST3 can be used to model primary pollutants from point, line, area, and volume sources. ISCST3 contains algorithms for effects such as building downwash and pollutant removal (wet scavenging and dry deposition), and is capable of modeling either gaseous species or particulate matter. 4.2 Basic concepts of atmospheric dispersion of pollutants After pollutants are emitted from a power plant stack, they form a plume and are transported by the wind. The central axis of the plume is aligned with the direction determined from the wind direction measured at a nearby weather station. The wind speed measured at the weather station is the model input that is used to determine the time of the plume arrival at each receptor. According to the Gausssian assumptions, the lateral and vertical spreads of the plume are selected according to atmospheric stability estimated from data provided by the weather station, upper-air stations, and the time of the day. High stability results in a less dispersed plume (usually during nighttime and or at lower temperatures). High stability causes high ground-level concentrations at receptors. On the other hand, during daytime and high temperatures with strong convective mixing of the atmosphere, the pollutants are efficiently mixed within a deep atmospheric boundary layer. This mixing results in lower ground-level concentrations. The concentrations are inversely proportional to the wind speed and consequently during low-wind conditions there is less atmospheric mixing and ground-level concentrations are high. With higher wind speeds, the pollutants are more mixed and the surface concentrations are low. Also, simulated concentrations are proportional to the emission rate the higher the emission rate the higher are the concentrations. Due to the continuous lateral and vertical dispersion, the pollutants are distributed in the space and the concentrations decrease with the distance from the power plant. The plume will affect receptors which are aligned with the wind direction measured at the weather station more than receptors that are down wind. 4.3 Model inputs We used an EPA approved dispersion model ISC3ST to simulate the transport and dispersion of pollutants emitted from the Tong Liang power plant (29.8350 N; 106.0516 E; 880 m MSL). Based on the given information, the emission parameters were set up as follows: 4

- Emission rate of SO 2 : 76.74 g s -1 - Release height above the ground: 85 m - Stack gas exit temperature: 423 K - Stack gas exit velocity: 2.97 m s -1 - Stack diameter: 1.59 m The input surface meteorology was obtained from the nearest weather station (29.8409 N; 106.0521 E; 1107 m AGL). Unfortunately, the station reported only every 6 hours and we had to interpolate the data to hourly intervals as required by the dispersion model for wind speed, wind direction, temperature, and cloudiness. Using the wind speed, solar radiation, cloudiness, and the time of day/night, we computed the stability classes (A-F), which are also a required input. Required upper-air data were obtained from the nearest radiosonde station in Chongqing (29.52 N; 106.48 E; elevation 351 m MSL). Within the global observational network, radiosonde stations send twice daily balloons with sensors to measure vertical profiles of the winds, temperature, humidity, and pressure from the surface to the top of the atmosphere. The main parameter obtained from the radiosonde measurements was the mixing depth (height of the atmospheric boundary layer in which the transport and dispersion of the emitted pollutants occur) as a required input to the dispersion model. 4.4 Required Model Setup Parameters According to specifications provided by the Columbia University (file DIST_BEAR.xls), coordinates and elevation of the source and 170 receptors used for the epidemiological study were specified. 5. Model results There are three major components, namely emissions, input meteorology, and the model structure, that influences the accuracy of the model results. We were given a constant emission rate; however, in reality the rate can change diurnally, monthly, and seasonally as well as due to unforeseen industrial operation changes. Since there is no documentation on these possible changes, the approximation of the constant emission rate has inherent uncertainty that influences magnitude of predicted concentrations in an unknown way. Regarding input meteorology, we have used available data. Hourly ground-level meteorological data from a weather station near the source as well as vertical sounding by radiosonde every 12 hours also near the source are preferred. Since the weather station data are recorded every six hours, we have interpolated the data into hourly intervals. This increases the uncertainty of estimated pollutant concentrations. Additionally, the data from the weather station should be quality assured, evaluated, and calibrated. The nearest radiosonde station is Chongqing, which is about 60 km southeast from the power plant. Considering complex terrain in the region, this separation can induce further uncertainty in the model results. Regarding the model, ISCST3 is an approved EPA model that can be applied to these conditions. Since ISCST3 is based on a Gaussian approximation of the plume spreading, it has uncertainties for complex terrain. 5

The analysis of the modeled ground-level concentrations at specified receptors and dense polar grid for each year (2002 and 2003) includes: 1. Annual average SO 2 concentrations 2. Annual maximum monthly SO 2 concentrations 3. Average SO 2 concentrations for each month of the power plant operation 4. Annual maximum daily SO 2 concentrations 5. Annual maximum 3-hr SO 2 concentrations In order to understand and interpret the model results and concentration patterns, we analyzed input wind direction and wind speed measured at the weather station. Figure 3 shows frequency distributions of the wind direction for 2002 and 2003, respectively. Both figures show dominant winds from the north which implies that the concentrations of the transported and dispersed and pollutants will be at maximum mainly to the south of the power plant. It should be also noticed that the flows are mainly channeled in the north-south direction, so the highest concentrations will be along this line. This will be obvious for long-term averages such as for an annual cycle. For shorter times, high concentrations can be present for other directions which might occur with less frequency within the annual cycle, but still can make an impact over shorter intervals. Figure 3. Frequency distribution of the wind direction at the Tong Liang weather station for 2002 (upper panel) and 2003 (lower panel). 6

The analysis of the wind speed shows that during 2002 and 2003 the station recorded zero winds (calm conditions) for 53% and 60% of the hours, respectively, and that significantly influences the model results. This number of calm conditions is quite unusual. Future studies should provide more accurate hourly meteorological parameters that could improve model results. 5.1 Study-area receptors and additionally specified receptors in a polar system The use of an additional dense receptor grid in a polar system allows for a symmetric coverage of the entire area around the power plant and helps understanding of the main spatial characteristics of the simulated concentrations. Notice that this polar receptor grid is added to the already specified grid of 170 study-area receptors. These 170 study-area receptors would not generally cover the entire area around the source. Figures 4 and 5 show locations of the 170 study-area receptors and additional receptors in a polar system. The polar receptors were selected in the following way. At six distances from the source (100, 500, 1000, 2000, 4000, and 7000 m), receptors were chosen for each 10 degrees of azimuth. Therefore for each of the distances, there are 36 receptors. There are 216 additional polar receptors to cover area where there are no field study receptors. In summary, there are 170 receptors to represent the residency of study subjects and 216 polar receptors that we specified to examine the spatial distribution of plume impact. Consequently, there are total of 386 receptors that were used for the analysis. Figure 4. Locations of the 170 study-area receptors within the 5 km from the source. The power plant is positioned in the origin. 7

Figure 5. Locations of the additional 216 receptors in a polar system within 7 km from the source. These additional receptors facilitate spatial analysis. The power plant is positioned in the origin. 5.2 Statistics of extreme simulated SO 2 concentrations Table 1 summarizes highest SO 2 concentrations simulated for 2002 and 2003. Notice that the highest values occur at receptors downwind and in close proximity to the power plant. Almost all highest concentrations are at spatial location 35 with bearing 180 from the north and 377 m from the power plant. Notice that the bearing 180 is along of the prevailing wind direction, as noted above (see the following sections 5.3-5.6). Locations 35 and c1125 create favorable conditions for extremely high concentrations. A study subject c1125 at which one of the extremes is simulated is located 786 m from the power plant and also is aligned with the prevailing flow channel (bearing is 189 ). SO2 statistics of concentrations are listed in Appendix A for 2002, and Appendix B for 2003. Table 1. Extreme SO 2 concentrations simulated for 2002 and 2003. Concentrations 2002 2002 2003 2003 SO2(ug/m3) Location Location Max. annual mean 53.96 Study area 35 47.64 c1125 Max. monthly mean 106.04 Study area 35 113.77 Study area 35 Max. daily mean 282.4 Study area 35 319.49 Study area 35 Max. 3-hr mean 643.76 Study area 35 645.88 Study area 35 8

5.3 Annual average SO 2 concentrations Figure 6 shows spatial distribution of the annual average ground-level concentrations of SO 2 for 2002 and 2003, Figure 6. Simulated annual average SO 2 ground-level concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel). 9

The highest concentrations are mainly located to the south of the source and in the proximity of the source for several other directions as a consequence of the prevailing wind direction. Notice that the extreme concentrations are simulated in a narrow band due to the channeling of the winds in the area. Also, the receptors selected for the epidemiological health-impact study are almost entirely located to the north of the power plant where concentrations are generally lower. Due to the frequent calm conditions, high concentrations are also found near the power plant for other directions. During the calm conditions, the wind direction is undefined and the pollutants accumulate near the source. Appendices A and B show values of annual mean concentrations for 2002 and 2003, respectively for all receptors. 10

5.4 Annual maximum monthly SO 2 concentrations Annual maximum SO 2 monthly concentrations are shown in Fig. 7. Figure 7. Simulated maximum monthly average SO 2 ground-level concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel). The highest monthly average during the plant operation is included in Table 1. The maximum monthly concentration pattern follows the structure of the annual average concentrations (Fig. 6); however, the magnitude is higher and the band of the maximum concentrations is wider compared to the annual average pattern. This reflects 11

the differences in monthly wind frequencies as shown in Fig. 3. Since the concentrations can be smaller during other months with higher winds and less stable conditions, annual concentrations can be still fairly small. However, in most of the cases, receptors which have high monthly maximum concentration will also have high concentrations for other periods. Due to seasonal effects of the meteorology, most of the variability is present for individual monthly averages as shown in the appendices A and B and the supplemental data base. 5.5 Annual maximum daily SO 2 concentrations Figure 8 shows simulated maximum daily-averaged SO 2 ground-level concentrations. Figure 8. Simulated maximum daily average SO 2 ground-level concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel). 12

The figure indicates that on a daily interval there can be a significant contribution to the concentration field around the source. Notice also that the magnitude is generally an order of magnitude greater than the monthly average maximum. One can see that there are noticeable differences in maximum daily averages between the 2002 and 2003 results due to more variability in the daily meteorology. Also, the differences in the meteorology between 2002 and 2003 led to occurrence of some simulated concentrations to the north of the power plant which are not significant for 2002. 5.6 Annual maximum 3-hour SO 2 concentrations Figure 9 shows simulated maximum 3-hr averaged SO 2 ground-level concentrations for 2002 and 2003. Figure 9. Simulated maximum 3-hour average SO 2 ground-level concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel). 13

Three-hour average concentrations are more than twice the average daily concentrations for most directions within two kilometers from the source. During the calm no-wind conditions, the pollutants accumulate near the source and contribute to very high concentrations. For the dominant flow channeling in the southern region, high concentrations are estimated even beyond 6 km from the source. 14

6. Additional presentation of the simulation results in the proximity of the power plant overlaid with the Tong Liang city map In order to present detailed dispersion modeling results for the main city center around the power plant, we overlaid the same dispersion results (shown in Section 5) on the city map for a 2 km radius around the power plant. By this process, one can see the distribution of the concentrations with respect to the city infrastructure and locations of monitoring stations. Figure 10 shows the distribution of the mean annual SO 2 concentrations for 2002 and 2003 overlaid over the city map. Fig. 10. Annual average SO 2 concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel) overlaid on the Tong Liang city map. 15

The distribution of simulated concentrations is similar for both years to the south an west from the plant. However, in 2003 there was more plume transport over the northeast residential area. - Maximum monthly mean SO 2 concentrations Fig. 11. Maximum monthly average SO 2 concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel) overlaid over the Tong Liang city map. 16

For shorter time periods compared with the annual, there are more chances to transport and disperse pollutants in more sectors around the power plant. - Maximum daily mean SO 2 concentrations Fig. 12. Maximum daily average SO 2 concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel) overlaid over the Tong Liang city map. For daily averages, power plant contributions are evident at nearly every location. Notice that the concentrations are more uniform around the power plant (more in 2003 than in 2002) and they drop off quite rapidly with distance for all segments except in the southern section. 17

- Maximum 3-hr mean SO 2 concentrations Fig. 13. Maximum 3-hr average SO 2 concentrations (μg m -3 ) for 2002 (upper panel) and 2003 (lower panel) overlaid over the Tong Liang city map. For the shortest intervals, concentrations are uniform around the power plant. Concentrations are more uniform for 2003 than for 2002. The central area is mostly pollutant free since the plume from an elevated source needs to travel for some distance before dispersing to ground-level. 18

7. Additional statistics helpful for the health study Owing to the uncertainties in input data and model formulation, SO 2 concentration estimates are best used in a relative rather than in an absolute sense. Monthly and annual averages are also more accurate than shorter averaging periods. For example, one could look at the power station at-home exposure of a study during the first month or two of pregnancy and compare this with the exposure of the study subjects. These relations can be easily examined using the supplemental spread sheets. In Appendix C, we sorted the model results at all receptors by the magnitude of the simulated concentrations (highest to lowest) as a summary of the ranking for each type of average. In numbers, this means that the receptor with highest concentrations is assigned number 1, the next highest with 2, and so on. The receptor with the lowest concentration is assigned number 170. This represents the summary ranking of the receptors depending on the strength of the pollution impact taking into consideration each type of the average. For each category (annual mean, monthly maximum, daily maximum, and 3-hr maximum concentrations), we assign this ranking and sum up all ranks for each receptor. By that, we come up with the final ranking in which the lower the sum the higher overall impact can be expected at that particular receptor. We calculate the ranking for 2002 and 2003 and sum up the total ranking for both years. This information can be useful since the health data can be also ranked in the same manner and the rankings for the pollution impact can be correlated with the ranking for the health data. The details of the final ranking computation are: Start with the supplemental excel files with estimated concentrations for 2002 and 2003. Sort all receptors by the magnitude of annual average concentrations from the highest to the lowest values. Assign ranking numbers to all receptors: the receptor with the highest concentration is assigned 1, next highest concentration is assigned 2, and so on. The receptor with the lowest concentration is assigned with 170. Store these numbers for each receptor as annual ranking number. In the next step repeat the ranking procedure explained above, but for monthly maximum concentrations. Assign this new monthly maximum ranking number to each receptor. Add this monthly maximum ranking number to the annual average ranking number obtained in the previous step. In the next step repeat the ranking procedure explained above, but for daily maximum concentrations. Assign this new daily maximum ranking number to each receptor. Add this daily maximum ranking number to the annual average ranking number and monthly maximum ranking number obtained in the previous two steps. In the final step repeat the ranking procedure explained above, but for 3-hr maximum concentrations. Assign this new 3-hr maximum ranking number to each receptor. Add this daily maximum ranking number to the annual average ranking number, monthly maximum ranking number obtained in the previous three steps. The final ranking number for each receptor is the sum of the annual average ranking number, the monthly maximum average ranking number, the daily maximum ranking number, and the 3-hr maximum ranking number. 19

Consequently, the smaller the final ranking number at a receptor the higher the impact at that receptor. 8. Conclusions In order to support an ongoing epidemiological health-effect study, we have performed a dispersion modeling study of the transport and dispersion of SO 2 emitted from the Tong Liang power plant. The dispersion modeling study for the periods of power plant operation in 2002 and 2003 was conducted using the EPA regulatory dispersion model Industrial Source Complex 3 Short Term (ISC3ST). The inputs for the dispersion model were data from a nearby meteorological station and from the nearest upper-air radiosonde station in Chongqing located approximately 60 km from the power plant. The modeling results consist of annual average, monthly maximum, daily maximum, 3-hr maximum, and monthly averages of SO 2 concentrations for 2002 and 2003 at each of 170 study subject receptors. Additionally, we defined 216 polar receptors covering six distances (100, 500, 1000, 2000, 4000, and 7000 m) with 36 azimuth locations. These 386 receptors provide details on the spatial distribution of the concentrations that helps us to understand the overall structure and evolution of the plume for the entire area. Monthly averages are shown for the months when the power plant was operational (December through May). For the months when the power plant was not in operation there were no emissions and the simulated concentrations are equal to zero. Whatever concentrations are measured during the time when the power plant was not operational are due to other emission sources in the area and in the region. Due to the prevailing winds blowing from the north to the south, most of the plume is transported and dispersed in the area south of the plant. This is mainly evident for the longer term averages (annual). Since most of the study subject receptors were north of the plant, they do not show concentrations as high as those south of the plant. The other important meteorological condition is a large occurrence of calm and low wind speeds. This causes high concentrations in the vicinity of the plant where most of the emissions accumulate. Whether these meteorological conditions are accurately represented by the nearby weather station should be evaluated with additional colocated measurements and other stations. The analysis of the wind speed shows that during 2002 and 2003 the station recorded zero winds (calm conditions) for 53% and 60% of hours, respectively. The consequence of these calm conditions is that some of the highest concentrations are estimated in the vicinity of the power plant. Monthly averages show a wider spatial distribution of high concentrations than the annual average. This is due to the natural variability (for short time intervals days or hours) of the wind direction and atmospheric stability superimposed on the prevailing winds. This is even more evident when the winds are light and of variable directions. In these cases the dispersion is less pronounced and higher concentrations can occur in many directions around the plant. The simulated maxima in the monthly averages, daily averages, and 3-hr averages are over 100, 300, and 600 μg m -3, respectively. High concentrations are simulated with similar values for both years. As shown in Section 5, the meteorology (mainly through 20

the analysis of the wind speed and direction) appears to be similar in 2002 and 2003. This implies that the impact of the power plant can be expected to be similar from year to year. Sorting the model results for all receptors by the magnitude of the simulated concentrations (highest to lowest) can be used as a measure of exposure to the power station emissions. The ranking was done for each category (annual mean, monthly maximum, daily maximum, and 3-hr maximum concentrations), and then all ranks were summed up for each receptor and for both years. That way, we developed a ranking in which the lower the sum the higher overall impact can be expected at that particular receptor. We believe that this information can be useful since the health-impact data can be also ranked in the same manner and the rankings for the pollution impact can be correlated with the ranking for the health data. Similar rankings can be performed with additional information about the study subjects, e.g., plume exposure during the first month of pregnancy vs. exposure of other subjects. We believe that the results of this study can be further analyzed together with health data in possible future projects. This study results could be also used in the future to determine risk assessment and health exposure levels due to pollution from the power plant. References Chow, J. C., J. G. Watson, L.-W. A. Chen, D. Koracin, B. Zielinska, D. Tang, F. Perera, J. Cao, and S.C. Lee, 2006: Exposure to PM2.5 and PAHs from the Tong Liang, China - Epidemiological Study. J. Environ. Sci. Health, Part A (in print). Koracin, D., D. Podnar, V. Isakov, J. Chow, Y. Dong, A. Miller, and M. McGown, 2000: PM 10 dispersion modeling for Treasure Valley, Idaho. J. Air Waste Manage. Assoc., 50, 174-185. Perera, F., D. Tang, R. Whyatt, S. A. Lederman, W. Jedrychowski, 2005: DNA damage from polycyclic aromatic hydrocarbons measured by benzo[a]pyrene-dna adducts in mothers and newborns from Northern Manhattan, the World Trade Center Area, Poland, and China. Cancer Epidemiol. Biomarkers Prev., 14(3), 709-714. U.S. Environmental Protection Agency, 1995: Guideline on Air Quality Models. EPA- 450/2-78-027R; Research Triangle Park, NC. 21

Appendix A Spatial grid (preceded by c), study subject (no prefix), monitoring location coordinates, annual average, maximum 3-hr average, maximum daily average, maximum monthly average, and monthly average SO 2 contribution (μg m -3 ) for 2002. June through November are not included because the power station was not operating during this period. 2002 Tongliang ISC3 CONC CONC CONC CONC MAX CONC CONC CONC CONC CONC CONC Rec.# STUDY_ID WAYP. Distance Bearing ANNUAL 3hr MAX 24hr MAX MONTH Jan-02 Feb-02 Mar-02 Apr-02 May-02 Dec-02 (m) (deg) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) 1 c1094 181 979 246 1.45 135.59 36.46 3.49 1.82 0.92 3.49 0.51 2.10 0.31 2 c1001 050 1221 339 0.81 269.84 98.77 3.54 0.00 0.42 3.54 0.00 1.09 0.01 3 c1002 052 832 346 0.73 210.50 75.15 2.85 0.00 0.54 2.85 0.00 1.18 0.00 4 c1003 54 383 337 1.84 263.08 82.91 5.67 0.02 2.14 3.79 0.00 5.67 0.01 5 c1004 51 1053 341 0.83 281.95 102.41 3.68 0.00 0.45 3.68 0.00 1.06 0.00 6 c1005 053 1486 28 0.16 74.21 10.12 0.61 0.61 0.00 0.00 0.18 0.02 0.17 7 c1006 085 1715 324 0.52 39.04 15.34 1.63 0.03 0.50 0.88 0.00 1.63 0.21 8 c1007 068 1311 342 0.58 199.77 73.35 2.67 0.00 0.27 2.67 0.00 0.67 0.00 9 c1008 071 209 241 1.60 252.09 80.75 5.03 1.39 0.00 5.03 0.45 3.10 0.00 10 c1009 070 840 40 1.97 322.99 90.72 5.48 5.48 0.00 0.00 2.73 0.12 4.01 11 c1010 084 1003 334 1.06 314.54 113.90 4.13 0.01 0.71 4.13 0.00 1.79 0.04 12 c1011 075 1084 340 0.84 289.44 105.30 3.77 0.00 0.45 3.77 0.00 1.06 0.00 13 c1012 069 724 344 1.02 297.68 104.80 3.91 0.00 0.79 3.91 0.00 1.71 0.00 14 c1013 123 1027 310 1.79 188.71 64.69 5.89 0.72 1.48 1.99 0.01 5.89 1.19 15 c1014 95 1227 37 0.92 153.97 42.55 2.69 2.69 0.00 0.00 0.98 0.03 2.03 16 c1015 80 1077 356 0.15 47.73 5.97 0.46 0.00 0.15 0.33 0.00 0.46 0.00 17 c1016 78 1651 358 0.06 20.39 2.69 0.20 0.00 0.05 0.11 0.00 0.20 0.00 18 c1017 82 558 328 2.11 281.63 105.56 6.71 0.20 2.35 3.89 0.00 6.71 0.22 19 c1018 90 581 346 1.03 243.96 83.44 3.36 0.00 0.98 3.36 0.00 2.13 0.00 20 c1019 79 987 19 0.17 144.05 19.64 0.79 0.79 0.00 0.00 0.11 0.16 0.01 21 c1020 76 1543 340 0.52 185.63 68.53 2.48 0.00 0.21 2.48 0.00 0.56 0.00 22 c1021 81 721 2 0.18 101.24 12.66 0.73 0.04 0.19 0.20 0.01 0.73 0.00 23 c1022 77 1509 10 0.05 25.90 3.53 0.14 0.14 0.01 0.01 0.01 0.12 0.00 24 c1023 104 1130 348 0.38 105.20 38.37 1.52 0.00 0.26 1.52 0.00 0.62 0.00 25 c1024 99 1502 11 0.05 32.60 4.45 0.18 0.18 0.01 0.01 0.01 0.12 0.00 26 c1025 100 1063 7 0.08 46.09 5.76 0.30 0.11 0.04 0.05 0.01 0.30 0.00 27 c1026 105 1204 338 0.80 274.40 100.37 3.59 0.00 0.40 3.59 0.00 0.99 0.00 28 c1027 94 686 89 4.06 360.20 183.02 10.02 0.92 1.71 10.02 5.64 6.58 0.74 29 c1028 86 1083 332 0.94 253.40 92.19 3.42 0.02 0.68 3.42 0.00 1.75 0.06 30 c1029 98 1618 13 0.05 42.70 5.82 0.23 0.23 0.00 0.01 0.01 0.09 0.00 31 c1030 130 798 0 0.18 85.85 10.73 0.66 0.02 0.19 0.24 0.00 0.66 0.00 32 c1031 132 890 56 2.43 241.08 73.90 7.00 3.68 0.00 0.38 7.00 2.29 1.98 33 c1032 88 585 316 2.84 347.58 136.92 9.75 1.11 2.98 3.06 0.01 9.75 1.01 34 c1033 87 846 307 2.27 237.64 66.07 6.81 1.41 1.70 2.77 0.04 6.81 1.53 35 c1034 146 1922 36 0.38 64.18 17.58 1.14 1.14 0.00 0.00 0.38 0.01 0.87 36 c1035 93 555 53 3.39 301.80 102.69 11.16 5.26 0.00 0.34 11.16 1.72 2.96 37 c1036 102 1521 342 0.46 162.88 60.11 2.20 0.00 0.19 2.20 0.00 0.51 0.00 38 c1037 91 649 31 0.75 98.04 24.36 2.00 2.00 0.00 0.00 1.50 0.10 1.13 39 c1038 92 559 41 2.53 384.73 103.55 6.15 6.15 0.00 0.01 5.36 0.18 4.21 40 c1039 89 624 12 0.19 100.09 12.51 0.59 0.42 0.06 0.05 0.09 0.59 0.00 41 c1040 115 1009 37 1.19 197.36 54.11 3.42 3.42 0.00 0.00 1.44 0.05 2.53 42 c1041 103 1232 340 0.71 249.03 91.19 3.27 0.00 0.34 3.27 0.00 0.84 0.00 43 c1042 161 1163 330 0.83 189.32 69.14 2.69 0.02 0.67 2.69 0.00 1.72 0.10 44 c1043 96 1350 37 0.80 134.42 37.22 2.36 2.36 0.00 0.00 0.82 0.02 1.79 45 c1044 97 1408 21 0.12 109.48 14.93 0.60 0.60 0.00 0.00 0.07 0.06 0.01 46 c1045 153 1915 345 0.24 77.36 28.70 1.12 0.00 0.09 1.12 0.00 0.30 0.00 47 c1046 129 904 353 0.28 62.97 16.11 0.79 0.00 0.28 0.79 0.00 0.71 0.00 48 c1047 117 1572 30 0.19 52.61 7.17 0.64 0.64 0.00 0.00 0.22 0.01 0.29 49 c1048 140 694 315 2.49 266.04 105.93 8.27 0.93 2.60 2.76 0.01 8.27 1.17 50 c1049 184 2315 342 0.25 86.64 32.23 1.22 0.00 0.07 1.22 0.00 0.26 0.00 51 c1050 136 1118 340 0.81 279.51 101.87 3.65 0.00 0.42 3.65 0.00 1.01 0.00 52 c1051 148 1322 8 0.05 28.73 3.59 0.18 0.11 0.02 0.03 0.01 0.18 0.00 53 c1052 133 342 32 0.63 158.58 40.37 2.68 0.78 0.00 0.00 2.68 0.12 0.45 54 c1053 186 1985 330 0.34 85.32 31.68 1.26 0.00 0.21 1.26 0.00 0.62 0.03 55 c1054 150 1741 358 0.05 18.23 2.51 0.18 0.00 0.04 0.10 0.00 0.18 0.00 56 c1055 119 1144 3 0.08 42.91 5.36 0.31 0.03 0.06 0.08 0.00 0.31 0.00 57 c1056 128 1054 324 1.06 106.04 41.93 3.20 0.11 1.20 1.80 0.00 3.20 0.41 58 c1057 141 355 339 1.59 206.85 65.57 4.98 0.01 1.94 3.13 0.00 4.98 0.00 59 c1058 178 197 206 2.68 206.43 49.20 7.04 1.06 1.14 7.04 2.49 4.62 0.48 60 c1059 184 2315 342 0.25 86.64 32.23 1.22 0.00 0.07 1.22 0.00 0.26 0.00 61 c1060 147 1387 12 0.06 44.98 6.13 0.24 0.24 0.01 0.01 0.01 0.13 0.00 62 c1061 154 1774 344 0.31 102.30 37.90 1.44 0.00 0.12 1.44 0.00 0.35 0.00 63 c1062 108 1521 3 0.04 23.96 3.00 0.18 0.02 0.03 0.05 0.00 0.18 0.00 64 c1063 122 1169 347 0.41 120.49 44.02 1.71 0.00 0.25 1.71 0.00 0.61 0.00 65 c1064 113 1216 28 0.22 94.73 12.92 0.83 0.83 0.00 0.00 0.28 0.04 0.26 66 c1065 118 1129 7 0.07 40.88 5.11 0.26 0.10 0.04 0.04 0.01 0.26 0.00 67 c1066 112 1040 23 0.19 147.82 20.16 0.86 0.86 0.00 0.00 0.18 0.10 0.06 68 c1067 135 1178 31 0.37 66.30 13.19 1.18 1.18 0.00 0.00 0.47 0.03 0.65 69 c1068 107 1523 4 0.04 23.55 2.94 0.17 0.02 0.03 0.05 0.00 0.17 0.00 70 c1069 160 1626 341 0.45 161.43 59.68 2.18 0.00 0.17 2.18 0.00 0.48 0.00 22

Appendix A - Simulated SO 2 concentrations (μg m -3 ) for 2002 continued (II). 70 c1069 160 1626 341 0.45 161.43 59.68 2.18 0.00 0.17 2.18 0.00 0.48 0.00 71 c1070 110 1118 21 0.16 139.04 18.96 0.77 0.77 0.00 0.00 0.11 0.10 0.02 72 c1071 180 1647 227 8.04 273.81 145.13 18.03 3.41 3.17 18.03 13.36 8.98 3.75 73 c1072 106 1535 342 0.46 160.77 59.35 2.18 0.00 0.19 2.18 0.00 0.50 0.00 74 c1073 114 1128 25 0.19 128.92 17.58 0.82 0.82 0.00 0.00 0.21 0.07 0.11 75 c1074 130 798 0 0.18 85.85 10.73 0.66 0.02 0.19 0.24 0.00 0.66 0.00 76 c1075 106 1535 342 0.46 160.77 59.35 2.18 0.00 0.19 2.18 0.00 0.50 0.00 77 c1076 121 1170 6 0.07 38.99 4.87 0.26 0.07 0.04 0.05 0.01 0.26 0.00 78 c1077 131 576 0 0.29 142.09 17.76 1.10 0.02 0.34 0.37 0.00 1.10 0.00 79 c1078 116 1224 36 0.80 132.90 36.36 2.35 2.35 0.00 0.00 0.87 0.02 1.75 80 c1079 109 1348 15 0.09 74.68 10.18 0.40 0.40 0.00 0.01 0.03 0.11 0.00 81 c1080 109 1348 15 0.09 74.61 10.17 0.40 0.40 0.00 0.01 0.03 0.11 0.00 82 c1081 142 332 330 2.18 373.89 135.81 8.61 0.16 2.60 2.38 0.00 8.61 0.03 83 c1082 137 1154 341 0.73 253.81 92.65 3.34 0.00 0.37 3.34 0.00 0.89 0.00 84 c1083 155 2166 349 0.11 26.14 9.51 0.45 0.00 0.05 0.45 0.00 0.20 0.00 85 c1084 125 1238 321 1.06 93.27 35.38 3.63 0.12 1.10 1.36 0.00 3.63 0.52 86 c1085 174 1254 39 1.14 191.45 54.19 3.34 3.34 0.00 0.00 1.19 0.05 2.55 87 c1086 152 1560 353 0.11 27.79 5.97 0.34 0.00 0.08 0.34 0.00 0.28 0.00 88 c1087 213 1796 18 0.07 69.73 9.51 0.38 0.38 0.00 0.00 0.02 0.04 0.00 89 c1088 126 1066 325 1.00 100.01 39.28 2.89 0.09 1.10 1.90 0.00 2.89 0.34 90 c1089 134 300 3 0.27 173.13 21.64 1.20 0.01 0.27 0.25 0.01 1.20 0.00 91 c1090 111 1089 36 0.94 155.77 42.43 2.74 2.74 0.00 0.00 1.10 0.03 2.01 92 c1091 127 1203 347 0.39 115.58 42.29 1.64 0.00 0.24 1.64 0.00 0.58 0.00 93 c1092 108 1521 3 0.04 23.96 3.00 0.18 0.02 0.03 0.05 0.00 0.18 0.00 94 c1093 124 1064 309 1.70 193.29 62.15 5.52 0.74 1.35 1.91 0.01 5.52 1.19 95 c1095 151 1710 352 0.11 27.14 6.68 0.36 0.00 0.07 0.36 0.00 0.25 0.00 96 c1096 1600 355 0.08 22.76 3.79 0.24 0.00 0.07 0.21 0.00 0.24 0.00 97 c1097 149 1563 1 0.05 23.02 2.88 0.19 0.01 0.04 0.07 0.00 0.19 0.00 98 c1098 180 1647 227 8.04 273.81 145.13 18.03 3.41 3.17 18.03 13.36 8.98 3.75 99 c1099 159 1331 2 0.06 31.91 3.99 0.24 0.01 0.05 0.08 0.00 0.24 0.00 100 c1100 206 363 336 1.85 248.47 87.65 6.10 0.03 2.22 3.38 0.00 6.10 0.01 101 c1101 143 572 29 0.59 106.85 21.79 1.50 1.50 0.00 0.00 1.38 0.15 0.68 102 c1102 144 562 44 2.99 430.03 122.23 7.06 7.06 0.00 0.02 6.66 0.36 4.72 103 c1103 145 1119 13 0.10 70.33 9.59 0.38 0.38 0.01 0.01 0.03 0.19 0.00 104 c1104 158 1312 355 0.12 31.53 4.77 0.34 0.00 0.11 0.30 0.00 0.34 0.00 105 c1105 162 1046 334 1.00 300.97 109.27 3.96 0.01 0.65 3.96 0.00 1.65 0.03 106 c1106 156 1337 354 0.13 30.45 6.04 0.36 0.00 0.11 0.36 0.00 0.35 0.00 107 c1107 138 1116 320 1.28 114.38 42.58 4.36 0.18 1.36 1.55 0.00 4.36 0.65 108 c1108 139 980 296 2.45 268.30 74.29 3.80 2.54 3.05 3.80 0.20 3.67 2.26 109 c1109 195 925 49 2.52 319.65 112.33 6.28 6.28 0.00 0.05 4.44 0.79 4.23 110 c1110 169 1144 3 0.08 42.91 5.36 0.31 0.03 0.06 0.08 0.00 0.31 0.00 111 c1111 170 1009 12 0.11 66.61 9.08 0.36 0.36 0.02 0.02 0.03 0.25 0.00 112 c1112 163 1162 321 1.14 99.97 38.34 3.83 0.14 1.22 1.47 0.00 3.83 0.56 113 c1113 190 1629 6 0.03 19.58 2.45 0.13 0.04 0.02 0.03 0.00 0.13 0.00 114 c1114 176 1432 37 0.73 123.13 34.12 2.16 2.16 0.00 0.00 0.73 0.02 1.65 115 c1115 172 1052 20 0.17 142.89 19.48 0.78 0.78 0.00 0.00 0.11 0.13 0.01 116 c1116 177 948 36 1.12 183.84 49.69 3.20 3.20 0.00 0.00 1.44 0.04 2.33 117 c1117 183 1846 331 0.40 109.66 40.66 1.56 0.00 0.23 1.56 0.00 0.66 0.03 118 c1118 171 1111 20 0.16 136.35 18.59 0.75 0.75 0.00 0.00 0.10 0.11 0.01 119 c1119 191 1574 360 0.05 22.73 2.84 0.20 0.00 0.04 0.08 0.00 0.20 0.00 120 c1120 179 658 343 1.20 340.27 118.42 4.44 0.00 0.97 4.44 0.00 2.14 0.00 121 c1121 187 1894 344 0.28 92.10 34.17 1.31 0.00 0.10 1.31 0.00 0.32 0.00 122 c1122 167 949 325 1.18 126.24 49.33 3.46 0.12 1.33 2.20 0.00 3.46 0.38 123 c1123 192 1526 11 0.05 31.77 4.33 0.17 0.17 0.01 0.01 0.01 0.11 0.00 124 c1124 182 1415 321 0.92 85.16 37.51 3.24 0.08 0.88 1.15 0.00 3.24 0.44 125 c1125 199 786 189 52.11 323.85 181.17 85.92 47.23 44.49 85.92 61.04 49.92 41.63 126 c1126 209 821 42 2.33 372.25 107.71 6.39 6.39 0.00 0.00 3.37 0.20 4.61 127 c1127 189 1800 3 0.03 16.80 2.10 0.13 0.01 0.02 0.04 0.00 0.13 0.00 128 c1128 173 1192 30 0.31 75.88 10.35 1.02 1.02 0.00 0.00 0.39 0.03 0.49 129 c1129 166 1034 332 1.01 266.96 96.86 3.60 0.02 0.74 3.60 0.00 1.91 0.07 130 c1130 188 1629 358 0.06 20.97 2.74 0.20 0.00 0.05 0.11 0.00 0.20 0.00 131 c1131 175 1212 37 0.93 156.59 43.26 2.74 2.74 0.00 0.00 1.01 0.03 2.06 132 c1132 165 1016 306 1.78 209.54 54.72 5.23 1.10 1.32 2.15 0.03 5.23 1.39 133 c1133 194 1116 39 1.33 221.81 62.57 3.85 3.85 0.00 0.00 1.48 0.06 2.90 134 c1134 168 1037 342 0.79 266.11 96.57 3.49 0.00 0.44 3.49 0.00 1.02 0.00 135 c1135 200 706 318 2.23 254.67 101.19 7.43 0.64 2.54 2.50 0.00 7.43 0.97 136 c1136 185 2077 343 0.27 91.10 33.84 1.29 0.00 0.09 1.29 0.00 0.29 0.00 137 c1137 164 1178 319 1.31 117.04 50.66 4.59 0.18 1.30 1.50 0.00 4.59 0.67 138 c1138 193 1130 17 0.13 112.40 15.33 0.61 0.61 0.00 0.00 0.06 0.14 0.00 139 c1139 196 932 50 2.48 295.26 108.11 5.98 5.98 0.00 0.06 4.65 0.93 3.95 140 c1140 198 344 305 4.39 468.63 174.65 14.15 4.09 2.23 6.41 0.20 14.15 0.43 23

Appendix A - Simulated SO 2 concentrations (μg m -3 ) for 2002 continued (III). 141 c1141 207 112 142 0.01 1.30 0.24 0.02 0.00 0.00 0.01 0.02 0.01 0.00 142 c1142 208 547 52 3.40 304.08 108.28 10.79 5.65 0.00 0.28 10.79 1.48 3.27 143 c1143 211 1224 36 0.80 132.90 36.36 2.35 2.35 0.00 0.00 0.87 0.02 1.75 144 c1144 214 1502 343 0.43 149.54 55.17 2.04 0.00 0.19 2.04 0.00 0.49 0.00 145 c1145 204 670 340 1.42 414.86 144.70 5.31 0.00 1.10 5.31 0.00 2.53 0.00 146 c1146 201 1064 336 0.97 315.74 114.75 4.11 0.01 0.57 4.11 0.00 1.42 0.02 147 c1147 202 1221 321 1.08 94.74 35.07 3.67 0.12 1.13 1.38 0.00 3.67 0.53 148 c1148 203 772 309 2.49 233.56 82.41 7.78 1.38 2.03 2.96 0.03 7.78 1.48 149 c1149 212 1371 31 0.29 54.52 10.08 0.92 0.92 0.00 0.00 0.34 0.02 0.50 150 c1150 205 263 315 2.76 472.36 176.63 11.52 1.39 2.02 2.39 0.01 11.52 0.04 151UN HOTEL 15 987 350 0.36 83.97 30.39 1.28 0.00 0.30 1.28 0.00 0.70 0.00 152STATION A 16 889 349 0.48 118.03 42.38 1.72 0.00 0.39 1.72 0.00 0.88 0.00 153STATION C 17 439 16 0.22 120.06 15.01 0.68 0.40 0.03 0.03 0.27 0.68 0.00 154STATION B 24 1515 339 0.56 199.43 73.58 2.65 0.00 0.23 2.65 0.00 0.61 0.00 155TINQING (S 18 8589 198 0.48 46.36 14.04 1.23 0.03 0.06 0.58 1.23 1.11 0.02 156 ANFONG ( 19 4290 178 43.53 263.94 184.09 56.83 52.57 38.39 56.83 41.89 35.00 51.01 157ARLONG (F 20 14673 77 0.03 12.31 1.61 0.10 0.00 0.00 0.01 0.06 0.10 0.00 158 NGLONG ( 22 14426 77 0.03 12.59 1.64 0.10 0.00 0.00 0.01 0.06 0.10 0.00 159 ANJIU (SI 23 17529 353 0.00 0.78 0.10 0.01 0.00 0.00 0.00 0.00 0.01 0.00 160 ATHER ST 25 656 4 0.19 115.21 14.40 0.79 0.07 0.18 0.17 0.01 0.79 0.00 161UDY AREA 26 1005 229 6.70 354.60 142.18 17.62 3.52 1.33 17.62 9.24 8.26 2.13 162UDY AREA 27 3242 334 0.19 62.75 23.37 0.90 0.00 0.06 0.90 0.00 0.22 0.00 163UDY AREA 28 2093 344 0.24 78.18 29.05 1.12 0.00 0.08 1.12 0.00 0.27 0.00 164UDY AREA 29 1372 14 0.08 63.67 8.68 0.34 0.34 0.01 0.01 0.02 0.12 0.00 165UDY AREA 30 687 93 3.95 347.20 176.21 9.89 0.86 1.57 9.89 4.47 7.43 0.65 166UDY AREA 31 1644 29 0.15 57.00 7.77 0.54 0.54 0.00 0.00 0.17 0.01 0.19 167UDY AREA 32 930 58 2.34 214.17 76.65 7.32 2.78 0.00 0.62 7.32 2.73 1.35 168UDY AREA 33 622 38 1.94 307.13 80.88 4.98 4.98 0.00 0.00 3.59 0.10 3.49 169UDY AREA 34 760 101 2.63 234.07 100.80 7.25 0.36 0.55 6.05 2.01 7.25 0.24 170UDY AREA 35 377 181 53.96 643.76 282.40 106.04 26.57 44.83 106.04 85.54 52.49 27.06 24

Appendix B - Spatial grid (preceded by c), study subject (no prefix), monitoring location coordinates, annual average, maximum 3-hr average, maximum daily average, maximum monthly average, and monthly average SO 2 contribution (μg m -3 ) for 2003. June through November are not included because the power station was not operating during this period. 2003 Tongliang ISC3 CONC CONC CONC CONC CONC CONC CONC CONC CONC CONC Rec.# STUDY_ID WAYP. Distance Bearing ANNUAL 3hr MAX 24hr MAX MXMONTH Jan-02 Feb-02 Mar-02 Apr-02 May-02 Dec-02 (m) (deg) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) (ug/m3) 1 c1094 181 979 246 2.44 336.83 74.42 6.6 0.11 1.27 6.08 6.6 4.27 0.05 2 c1001 050 1221 339 1.37 173.11 48.67 8.44 0.68 0.02 0.22 0.09 8.44 0.69 3 c1002 052 832 346 1.01 168.53 62.43 5.11 0.67 0.21 0.68 0.07 5.11 0.7 4 c1003 54 383 337 2.81 546.32 206.04 18.18 0.32 0.01 1.1 0.78 18.18 0.33 5 c1004 51 1053 341 1.22 129.58 44.36 7.12 0.73 0.04 0.32 0.1 7.12 0.73 6 c1005 053 1486 28 0.52 120.64 30.9 1.47 0.12 0.1 1.47 1.3 0.82 0.07 7 c1006 085 1715 324 1.04 178.16 38.63 3.99 0.64 0.23 0.08 1.77 3.99 1.02 8 c1007 068 1311 342 1.02 127.08 33.3 6.17 0.55 0.04 0.21 0.05 6.17 0.55 9 c1008 071 209 241 1.4 240.06 55.57 3.76 0.04 2.69 2.54 1.63 3.76 0 10 c1009 070 840 40 3.61 327.92 118.51 8.43 1.75 2.77 8.17 8.43 4.58 1.5 11 c1010 084 1003 334 1.59 251.65 67.02 9.47 0.76 0.04 0.3 0.36 9.47 0.8 12 c1011 075 1084 340 1.31 133.95 42.4 7.77 0.74 0.03 0.29 0.11 7.77 0.75 13 c1012 069 724 344 1.3 228.96 84.93 7.09 0.76 0.14 0.73 0.14 7.09 0.77 14 c1013 123 1027 310 2.36 250.19 111.69 6.95 1.06 0.73 1.03 6.69 6.95 1.26 15 c1014 95 1227 37 1.98 159.31 68.62 4.79 1.26 1.43 4.11 4.79 2.54 0.79 16 c1015 80 1077 356 0.84 127.51 42.47 2.09 0.71 0.75 1.14 0.01 1.57 2.09 17 c1016 78 1651 358 0.78 226.8 88.17 3.32 0.48 0.52 0.66 0 0.82 3.32 18 c1017 82 558 328 2.95 431.91 123.63 17.66 0.44 0.25 0.72 2.3 17.66 0.47 19 c1018 90 581 346 1.5 309.91 115.34 8.35 0.64 0.21 1.09 0.14 8.35 0.66 20 c1019 79 987 19 0.53 139.93 55.83 2.47 0.12 0.05 2.47 0.75 0.53 0.04 21 c1020 76 1543 340 1.22 199.89 50.38 7.83 0.5 0.01 0.14 0.04 7.83 0.5 22 c1021 81 721 2 1.01 169.66 43.81 2.56 0.91 0.98 2.56 0.03 1.71 1.37 23 c1022 77 1509 10 0.33 73.9 17.15 0.83 0.22 0.23 0.83 0.04 0.51 0.59 24 c1023 104 1130 348 0.63 85.33 31.56 2.82 0.48 0.24 0.49 0.02 2.82 0.62 25 c1024 99 1502 11 0.3 74.79 18.38 0.84 0.18 0.19 0.84 0.05 0.51 0.42 26 c1025 100 1063 7 0.62 97.89 29.76 1.61 0.53 0.57 1.61 0.05 0.59 1.26 27 c1026 105 1204 338 1.45 179.25 53.4 8.92 0.7 0.01 0.22 0.12 8.92 0.71 28 c1027 94 686 89 3.87 433.16 181.3 8.61 8.61 0 7.54 6.38 0.71 5.55 29 c1028 86 1083 332 1.47 275.4 68.97 8.79 0.63 0.07 0.25 0.45 8.79 0.71 30 c1029 98 1618 13 0.23 71.81 17.49 0.74 0.1 0.11 0.74 0.07 0.49 0.18 31 c1030 130 798 0 1 173.21 30.17 2.11 0.93 1.01 2.11 0.02 1.76 1.69 32 c1031 132 890 56 2.03 160.24 47.44 5.02 1.26 1.17 4.01 5.02 3.15 0.68 33 c1032 88 585 316 3.56 372.92 137.07 16.88 0.17 0.92 0.85 7.55 16.88 0.26 34 c1033 87 846 307 2.3 271.72 95.52 7.85 0.33 0.61 1.76 6.26 7.85 0.41 35 c1034 146 1922 36 0.97 107.65 32.53 2.3 0.78 0.66 1.73 2.3 1.47 0.35 36 c1035 93 555 53 3.82 331.08 111.77 8.74 1.53 2.19 8.74 8.67 6.54 1 37 c1036 102 1521 342 0.96 143.14 36.18 6.01 0.45 0.03 0.15 0.03 6.01 0.46 38 c1037 91 649 31 2.4 289.11 116.59 7.76 0.67 0.79 7.76 5.66 2.67 0.4 39 c1038 92 559 41 4.79 424.36 172.1 12.24 1.57 3.13 12.24 10.73 6.87 1.43 40 c1039 89 624 12 0.88 267.23 102.6 4.31 0.37 0.32 4.31 0.45 0.78 0.32 41 c1040 115 1009 37 2.45 199.55 90.81 5.91 1.25 1.72 5.54 5.91 3.06 0.97 42 c1041 103 1232 340 1.28 159.64 43.24 7.8 0.65 0.02 0.22 0.08 7.8 0.66 43 c1042 161 1163 330 1.32 286.29 66.94 7.72 0.5 0.11 0.21 0.6 7.72 0.65 44 c1043 96 1350 37 1.78 140.87 59.31 4.28 1.25 1.29 3.53 4.28 2.32 0.7 45 c1044 97 1408 21 0.33 74.11 28.53 1.27 0.05 0.01 1.27 0.53 0.55 0.01 46 c1045 153 1915 345 0.48 63 16.73 2.87 0.25 0.04 0.12 0.01 2.87 0.27 47 c1046 129 904 353 0.85 111.48 38.58 2.59 0.7 0.65 1.16 0.01 2.59 1.2 48 c1047 117 1572 30 0.58 125.44 33.45 1.48 0.16 0.18 1.46 1.48 0.93 0.12 49 c1048 140 694 315 3.22 350.46 129.66 13.61 0.42 1.02 0.85 7.62 13.61 0.62 50 c1049 184 2315 342 0.75 143.51 35.88 4.96 0.24 0.01 0.06 0.01 4.96 0.24 51 c1050 136 1118 340 1.3 136.38 41.95 7.77 0.72 0.03 0.27 0.1 7.77 0.72 52 c1051 148 1322 8 0.47 67.11 22.9 1.1 0.36 0.4 1.08 0.03 0.51 1.1 53 c1052 133 342 32 3.3 464 173.54 12.44 0.8 0.43 12.44 6.68 4.15 0.11 54 c1053 186 1985 330 0.77 139.42 39.03 4.73 0.25 0.04 0.07 0.19 4.73 0.41 55 c1054 150 1741 358 0.76 232.03 90.21 3.36 0.44 0.49 0.6 0 0.78 3.36 56 c1055 119 1144 3 0.8 139.11 54.07 2.43 0.69 0.75 1.37 0.01 0.72 2.43 57 c1056 128 1054 324 1.74 329.86 78.05 8.16 0.59 0.45 0.24 2.64 8.16 0.89 58 c1057 141 355 339 2.63 557.91 211.28 17.22 0.23 0.01 1.12 0.58 17.22 0.24 59 c1058 178 197 206 2.13 241.86 45.56 4.68 0.62 4.68 3.1 4.17 3.55 0.18 60 c1059 184 2315 342 0.75 143.51 35.88 4.96 0.24 0.01 0.06 0.01 4.96 0.24 61 c1060 147 1387 12 0.3 68.43 22.92 1.01 0.17 0.18 1.01 0.08 0.5 0.31 62 c1061 154 1774 344 0.63 89.55 22.39 3.84 0.31 0.04 0.12 0.01 3.84 0.32 63 c1062 108 1521 3 0.72 194.99 75.8 2.95 0.5 0.55 0.85 0.01 0.5 2.95 64 c1063 122 1169 347 0.64 82.81 30.62 3.08 0.47 0.18 0.41 0.03 3.08 0.56 65 c1064 113 1216 28 0.73 156.86 42.71 2.23 0.18 0.16 2.23 1.84 0.96 0.1 66 c1065 118 1129 7 0.59 92.17 26.23 1.45 0.5 0.54 1.45 0.04 0.56 1.31 67 c1066 112 1040 23 0.63 130.21 51.48 2.53 0.13 0.04 2.53 1.27 0.7 0.03 68 c1067 135 1178 31 1.05 187.76 56.7 2.82 0.33 0.41 2.82 2.68 1.36 0.26 69 c1068 107 1523 4 0.66 169.36 65.84 2.59 0.47 0.51 0.85 0.01 0.47 2.59 70 c1069 160 1626 341 1.08 177.93 44.48 6.88 0.44 0.01 0.13 0.03 6.88 0.44 25

Appendix B - Simulated SO 2 concentrations (μg m -3 ) for 2003 continued (II). 71 c1070 110 1118 21 0.48 114.08 44.83 2.05 0.09 0.03 2.05 0.83 0.58 0.02 72 c1071 180 1647 227 4.62 273.27 88.2 9.73 2.95 4.76 5.29 9.73 8.36 3.79 73 c1072 106 1535 342 0.96 143.85 36.29 5.99 0.45 0.02 0.15 0.03 5.99 0.45 74 c1073 114 1128 25 0.64 130.49 43.26 2.3 0.13 0.07 2.3 1.45 0.78 0.04 75 c1074 130 798 0 1 173.21 30.17 2.11 0.93 1.01 2.11 0.02 1.76 1.69 76 c1075 106 1535 342 0.96 143.85 36.29 5.99 0.45 0.02 0.15 0.03 5.99 0.45 77 c1076 121 1170 6 0.63 99.25 33.99 1.63 0.53 0.58 1.35 0.03 0.57 1.63 78 c1077 131 576 0 1.16 156.64 53.3 3.01 0.84 0.9 2.92 0.03 3.01 0.98 79 c1078 116 1224 36 1.77 169.05 67.21 4.31 1.01 1.2 3.81 4.31 2.3 0.68 80 c1079 109 1348 15 0.28 69.13 27.6 1.15 0.1 0.09 1.15 0.18 0.47 0.11 81 c1080 109 1348 15 0.28 69.03 27.55 1.15 0.1 0.09 1.15 0.18 0.47 0.11 82 c1081 142 332 330 3.36 520.1 197.86 21.8 0.11 0.04 0.96 1.83 21.8 0.12 83 c1082 137 1154 341 1.19 126.72 36.91 7.06 0.66 0.03 0.27 0.08 7.06 0.67 84 c1083 155 2166 349 0.25 40.57 11.52 1.07 0.16 0.1 0.15 0 1.07 0.35 85 c1084 125 1238 321 2.01 232.03 75.11 6.93 1.25 0.55 0.2 4.44 6.93 1.63 86 c1085 174 1254 39 2.4 192.92 68.39 5.74 1.88 1.88 4.6 5.74 2.97 0.98 87 c1086 152 1560 353 0.53 87.15 33.88 1.47 0.37 0.36 0.48 0 1.23 1.47 88 c1087 213 1796 18 0.18 44.74 17.08 0.69 0.03 0.02 0.69 0.17 0.44 0.02 89 c1088 126 1066 325 1.61 335.03 76.82 7.95 0.52 0.37 0.24 2.12 7.95 0.79 90 c1089 134 300 3 0.95 270.2 96.29 3.57 0.09 0.09 3.57 0.09 3.11 0.09 91 c1090 111 1089 36 2.04 193.93 79.74 4.97 1.04 1.37 4.61 4.97 2.58 0.78 92 c1091 127 1203 347 0.62 78.11 28.84 2.98 0.45 0.17 0.39 0.02 2.98 0.54 93 c1092 108 1521 3 0.72 194.99 75.8 2.95 0.5 0.55 0.85 0.01 0.5 2.95 94 c1093 124 1064 309 2.15 226.51 100.44 6.24 0.93 0.65 1.11 6.1 6.24 1.1 95 c1095 151 1710 352 0.43 70.47 24.94 1.17 0.3 0.27 0.36 0 1.17 1.1 96 c1096 1600 355 0.64 144.51 56.18 2.23 0.42 0.44 0.56 0 1.08 2.23 97 c1097 149 1563 1 0.8 233.69 90.85 3.45 0.52 0.57 0.79 0 0.6 3.45 98 c1098 180 1647 227 4.62 273.27 88.2 9.73 2.95 4.76 5.29 9.73 8.36 3.79 99 c1099 159 1331 2 0.8 190.34 73.99 3 0.61 0.67 1.05 0.01 0.65 3 100 c1100 206 363 336 2.91 557.37 210.8 19 0.25 0.01 1.09 0.89 19 0.26 101 c1101 143 572 29 2.36 331.14 128.53 8.55 0.61 0.5 8.55 5.25 2.48 0.24 102 c1102 144 562 44 5.01 459.76 166.29 12.28 1.71 3.51 12.28 11.09 7.5 1.61 103 c1103 145 1119 13 0.39 93.83 37.22 1.6 0.2 0.19 1.6 0.18 0.48 0.26 104 c1104 158 1312 355 0.71 124.88 48.55 2.11 0.54 0.57 0.78 0 1.32 2.11 105 c1105 162 1046 334 1.56 240.84 66.26 9.34 0.74 0.04 0.28 0.33 9.34 0.78 106 c1106 156 1337 354 0.65 102.83 39.97 1.79 0.5 0.5 0.69 0 1.38 1.79 107 c1107 138 1116 320 2.31 242.29 90.73 8.05 1.3 0.68 0.27 5.28 8.05 1.68 108 c1108 139 980 296 1.56 303.84 102.71 4.91 0.08 0.04 3.32 3.3 4.91 0.01 109 c1109 195 925 49 3.42 323.78 87.68 7.76 2.13 2.84 6.92 7.76 4.57 1.53 110 c1110 169 1144 3 0.8 139.11 54.07 2.43 0.69 0.75 1.37 0.01 0.72 2.43 111 c1111 170 1009 12 0.47 111.66 43.76 1.91 0.27 0.27 1.91 0.18 0.5 0.36 112 c1112 163 1162 321 2.08 250.85 77 7.51 1.16 0.59 0.23 4.5 7.51 1.53 113 c1113 190 1629 6 0.5 115.65 44.96 1.8 0.35 0.39 0.73 0.01 0.43 1.8 114 c1114 176 1432 37 1.66 130.14 54.07 3.99 1.23 1.21 3.2 3.99 2.19 0.65 115 c1115 172 1052 20 0.5 126.43 50.03 2.25 0.1 0.03 2.25 0.79 0.55 0.03 116 c1116 177 948 36 2.39 223.69 96.33 5.8 1.06 1.56 5.72 5.8 2.96 0.89 117 c1117 183 1846 331 0.93 146.19 45.58 5.83 0.3 0.03 0.08 0.16 5.83 0.43 118 c1118 171 1111 20 0.46 114.4 45.06 2.01 0.09 0.03 2.01 0.71 0.54 0.03 119 c1119 191 1574 360 0.81 239.62 93.15 3.53 0.52 0.57 0.77 0 0.67 3.53 120 c1120 179 658 343 1.51 274.95 102.11 8.48 0.78 0.11 0.79 0.19 8.48 0.8 121 c1121 187 1894 344 0.59 88.52 22.13 3.67 0.27 0.03 0.11 0.01 3.67 0.28 122 c1122 167 949 325 1.83 375.11 88.48 9.39 0.49 0.42 0.3 2.33 9.39 0.72 123 c1123 192 1526 11 0.29 75.61 17.77 0.82 0.18 0.19 0.82 0.05 0.51 0.42 124 c1124 182 1415 321 1.86 194.57 70.32 5.87 1.4 0.47 0.15 4.26 5.87 1.78 125 c1125 199 786 189 47.64 353.07 232.46 94.21 63.55 43.47 62.28 94.21 56.86 38.44 126 c1126 209 821 42 4 378.15 120.64 9.22 2.01 3.2 8.83 9.22 5.16 1.71 127 c1127 189 1800 3 0.66 210.55 81.86 3.04 0.4 0.44 0.63 0 0.41 3.04 128 c1128 173 1192 30 0.92 179.28 51.98 2.59 0.26 0.31 2.59 2.36 1.2 0.19 129 c1129 166 1034 332 1.52 289.25 70.85 9.04 0.66 0.07 0.28 0.49 9.04 0.73 130 c1130 188 1629 358 0.78 225.24 87.56 3.31 0.48 0.53 0.67 0 0.83 3.31 131 c1131 175 1212 37 2.01 161.76 69.9 4.86 1.26 1.45 4.19 4.86 2.57 0.8 132 c1132 165 1016 306 1.79 205.66 74.59 5.61 0.39 0.46 1.62 4.88 5.61 0.45 133 c1133 194 1116 39 2.66 223.88 81.19 6.37 1.8 2.07 5.42 6.37 3.3 1.12 134 c1134 168 1037 342 1.13 124.26 44.99 6.44 0.7 0.05 0.35 0.09 6.44 0.71 135 c1135 200 706 318 3.1 326.69 112.41 14.06 0.43 0.96 0.63 6.48 14.06 0.63 136 c1136 185 2077 343 0.67 115.56 28.89 4.31 0.26 0.02 0.08 0.01 4.31 0.26 137 c1137 164 1178 319 2.45 248.91 102.13 7.68 1.64 0.69 0.27 6.02 7.68 2.03 138 c1138 193 1130 17 0.39 104.24 41.54 1.76 0.1 0.07 1.76 0.4 0.47 0.07 139 c1139 196 932 50 3.19 299.02 79.31 7.22 1.97 2.61 6.46 7.22 4.33 1.42 140 c1140 198 344 305 3.96 510.08 175.22 17.84 0.07 0.11 0.98 10.57 17.84 0 26