摘 要

高雄市空氣品質PSI > 100之比例在民國84年為10.3%，在91年為5.5 %，雖有改善趨勢，但懸浮微粒(PM10)事件日之比例由民國91年的1.0 %上升至民國93年的2.9 % ，卻有微幅劣化的現象。本研究以統計盒型圖對民國86-93年高雄市4個空品測站之PM10濃度作一趨勢探討，再以t-test及F-test進行事件日統計分析，並搭配模式(CAMx)模擬懸浮微粒(PM10)污染來源及成因。

民國86-93年長期趨勢分析顯示，本市四測站PM10月平均濃度在72.9-81.7 μg/m3之間，以小港站最高，楠梓站最低。PM10長期趨勢則處在微輻下降中，年下降輻度楠梓站為1.05 %，左營站為1.38 %，前金站為1.51 %，小港站為1.91 %。

在民國86-93年PM10事件日站日數之空間分佈上，則呈南向北遞減趨勢(小港>前金>左營>楠梓)。民國86-93年四測站事件日平均值(t-test)及變異數(F-test)在95% 信心水準下之統計檢定結果顯示：在空間分佈上，四測站逐時PM10濃度之變化約42 % 有顯著差異(或不同)，即空間分佈的差異性較明顯。

本研究分析民國86-93年期間高雄市空氣品質及氣象監測數據，並以空氣品質模式(CAMx)模擬本市四個環保署空氣品質監測站(楠梓、左營、前金、小港)之懸浮微粒(PM10)事件日時空分佈。空氣品質模擬減量結果顯示，本市固定源對粒狀物之減量效果為38.9 % (NOx：24.7 %；SO2：14.2 %)，明顯優於移動源8.7 % (NOx：7.3 %；SO2：1.4 %)及逸散源0.9 % (SO2：0.9 %) 之減量效果。港區船舶排放量對粒狀物之減量效果為3.5 %。

關鍵字：懸浮微粒、趨勢分析、統計檢定、CAMx模式

ABSTRACT

Although the fractions of station-days that the Pollutant Standard Index (PSI) exceed 100 (also referring to the episodes) in Kaohsiung City showed a decline trend from about 10.3% in 1995 to about 5.5% in 2002, the percentage of particulate PM10 events showed a increase trend: from 1.0 % in 2002 to 2.9 % in 2004. This study first statistically summarized the trends of PM10 concentrations using box plots for four air-quality monitoring stations in Kaohsiung during the period of 1997 to 2004, together with the t-test and F-test. The Comprehensive Air Quality Model with extensions (CAMx model) was then applied to analyze the source and the cause of the PM10 events.

The monthly averages of PM10 concentrations at four air-quality monitoring stations were 72.9 – 81.7 μg/m3 during the period of 1997 to 2004, highest at Hsiung-Kong and lowest at Nan-Chie. The long-term trend analyses show slightly decline results for yearly-averaged PM10 concentrations (1.05% at Nan-Chie, 1.38% at Tzuo-Yin, 1.51% at Chien-Chin, and 1.91% at Hsiung-Kong).

During 1997 to 2004, the PM10 episodes occurred most frequently, while the numbers of PM10 episodes decreased from south to north (i.e., Hsiung-Kong > Chien-Chin > Tzuo-Yin > Nan-Chie). The statistical tests using t-test for the mean and F-test for the variance with 95% confidence level show that the probability that the hourly PM10 concentrations differ insignificantly among the four stations is only about 42%. That is, the spatial difference of pollutant concentrations among four air-quality monitoring stations is rather significant in Kaohsiung.

The CAMx simulations show that contributions to ambient PM10 from stationary source is about 38.9% (NOx: 24.7%; SO2: 14.2%), 8.8% from mobile source (NOx: 7.4%; SO2: 1.4%), and 0.9% from fugitive emissions (SO2: 0.9%) in Kaohsiung. The contributions to ambient PM10 from the emissions in Kaohsiung Harbor are about 3.5%.


Keywords: Particulate matter, Trend analysis, Statistical test, CAMx model.

謝 誌……………………………………………….…..……...I
摘 要…………………………………………………..……...II
ABSTRACT………………….………………………..….…...III
目 錄………………………………………………….….......IV
表目錄………………………………………...………...….…VI
圖目錄………………………………………………………...VII
第一章 前言……………....................................................................1-1
1.1 研究動機……………………....………...…...………..…...…...1-1
1.2 研究目的………………………………...……….........…...…...1-1
1.3 研究架構及流程………...…………...………..…………...…...1-3
第二章 相關研究及文獻回顧……………………………………2-1
2.1 高雄市84-93年空氣品質趨勢及空氣污染物排放量概述……2-1
2.1.1 各測站PSI變化趨勢……….……….….…………….……2-1
2.1.2 各空氣污染物濃度變化趨勢…………...…………………2-2
2.1.3 高雄市空氣污染物排放量概估...…………………....……2-7
2.1.4 高雄地區歷年氣象概述…………….….…………..….…2-10
2.2 大氣懸浮微粒形成機制…………………………………....…2-15
2.3 大氣懸浮微粒之特性…………………………………………2-16
2.3.1 大氣懸浮微粒之定義………………..………...…………2-16
2.3.2 懸浮微粒分類與來源…………………....…….…………2-16
2.3.3 懸浮微粒化學組成…………………..…..…….…………2-18
2.4 大陸沙塵暴長程傳輸對空氣品質的影響……………………2-20
2.5 懸浮微粒受氣象條件之影響………...……..…..….…………2-21
2.6 懸浮微粒對人體及環境之影響…………...………….………2-23
2.6.1 懸浮微粒對人體健康之影響……….………...………….2-23
2.6.2 懸浮微粒對環境之影響………………….…...………….2-25
第三章 研究方法…………………………………………….3-1
3.1 大氣中懸浮微粒(PM2.5)採樣與分析…..…………...….….…...3-1
3.1.1 大氣採樣規劃………………....…….…...…………..…….3-1
3.1.2 採樣儀器與分析方法…………………....…………..…….3-2
3.1.3 品保與品管……………………………....…………..…….3-4
3.2 平均值(t-test)及變異數(F-test)之統計檢定方法………..….....3-8
3.3 空氣品質模式(CAMX)理論礎……....……………...…....…...3-11
3.3.1 模式之制御方程式…….…...………………..……..…….3-11
3.3.2 二次氣懸膠反應機制…………......……..…………....….3-12
3.3.3 模擬區域及網格分配…………………………………….3-14
3.3.4 模式輸入資料概述……………………………………….3-15
3.3.5 垂直擴散係數(Kv)之設定……………………………….3-18
3.3.6 排放資料之設定…………..……....…………..………….3-19
3.3.7 邊界及初始條件之設定…………..…......……………….3-21
第四章 結果與討論………………………………………….4-1
4.1 民國86-93年高雄市懸浮微粒時空變化趨勢分析…….…..….4-1
4.1.1 民國86-93年PM10逐時濃度趨勢分析………..……..…....4-2
4.1.2 民國86-93年PM10逐日濃度趨勢分析……….……..…....4-2
4.1.3 民國86-93年PM10逐季濃度趨勢分析…………………...4-2
4.1.4 懸浮微粒月長期趨勢分析………………………………...4-6
4.2 高雄市四空品測站懸浮微粒事件日統計結果.……………...4-13
4.2.1 民國86-93年懸浮微粒事件日統計……………..……….4-13
4.2.2 民國86-93年四空品測站懸浮微粒事件日之平均值(t-test)及變異數(F-test)統計檢定結果……….………………..….4-14
4.3 懸浮微粒事件日之模擬分析結果.………..……………..…...4-16
4.3.1 懸浮微粒採樣分析結果及模式換算因子之校正…….....4-16
4.3.2 懸浮微粒事件日模擬結果與污染源減量效益分析….....4-19
4.3.3 逆軌跡分析結果……………..…………………………...4-26
4.3.4 其他研究結果之比較…………………...….…..………...4-26
第五章 結論與建議………………………………………….5-1
5.1 結論…………………………………….………….…….…….5-1
5.2 建議………………………………..……………….………….5-2
參考文獻…………………………………………………….參-1
附錄A 高雄市小港測站PM2.5懸浮微粒之平均濃度(μg/m3).…附A-1
附錄B PM2.5 RAAS2.5-100 sampler流量校正記錄表……....…附B-1
附錄C 離子層析儀NO3-、SO42-、NH4+標準品製備檢量線圖…附C-1
附錄D 長期趨勢分析中之擾動值檢定圖…………………....…附D-1
附錄E 民國86-93年懸浮微粒事件日統計表……………..……附E-1
附錄F 民國86-93年(1997-2004)高雄市四空品測站於懸浮微粒事
件日平均值(t-test)及變異數(F-test)之統計檢定結果....附F-1
附錄G CAMX二次氣懸膠(NO3-、SO42-、NH4+)模擬逐時等濃度圖……………………………….…………………….….附G-1

表 目 錄

表2.1-1 民國89年高雄市各污染源排放總表……………..….….2-9
表2.1-2 高雄地區近十年各月份平均氣象資料統計表…………2-11
表2.6-1 SO2 與PM10 並存對人體健康之影響………………….2-25
表3.1-1 高雄市小港空氣品質監測網測站環境資料一覽表……..3-1
表3.1-2 離子層析儀之儀器偵測極限……………………………..3-7
表3.1-3 離子層析儀添加添加標準品回收率……………..……....3-8
表3.3-1 CAMx輸入資料一覽表……...…...…………………..…3-16
表3.3-2 CAMx土地使用分類…….………………………..….…3-17
表3.3-3 垂直擴散係數(Kv)之設定……………..………....…...…3-19
表3.3-4 懸浮微粒模擬之邊界及初始條件濃度值……........……3-22
表4.1-1 民國86-93年PM10趨勢分析係數表…………..………....4-8
表4.2-1 民國86-93年高雄市四空品測站懸浮微粒(PM10)事件日每年之天數………………………………………………...4-14
表4.2-2 民國86-93年高雄市四空品測站事件日懸浮微粒(PM10)逐時濃度統計檢定結果彙整……………………..…....….4-15
表4.3-1 小港地區PM2.5中之三種無機鹽類分析結果……...…...4-17
表4.3-2 污染源對PM10減量效益模擬結果………………...……4-25
表4.3-3 將NOx及SOx排放量各削減5%、10%及15%時PM10日平均值之下降輻度.………………………………..…….……4-25

圖 目 錄

圖1.3-1 本研究架構及流程圖....………………….……..…..…..1-3
圖2.1-1 高屏地區空氣品質現況(84-93年)………………………2-1
圖2.1-2 高雄市一般測站84年至93年PSI＞100指標污染物逐月變化圖..……………………...………………………..….2-2
圖2.1-3 高屏地區84-93年逐月PM10濃度變化…….…..……......2-2
圖2.1-4 高屏地區84-93年逐月SO2濃度變化…………...…..…..2-3
圖2.1-5 高屏地區84-93年逐月NO2濃度變化……………..……2-3
圖2.1-6 高屏地區84-93年逐月O3濃度變化………………….…2-4
圖2.1-7 高屏地區84-93年逐月CO濃度變化……………...……2-4
圖2.1-8 高雄市各測站PM10逐月濃度變化趨勢 (1/2)….....……2-5
圖2.1-9 高雄市各測站PM10逐月濃度變化趨勢 (2/2)…...…..…2-6
圖2.1-10 民國93年春季高雄測站風瑰圖…………..………...…2-13
圖2.1-11 民國93年夏季高雄測站風瑰圖…….…………...….…2-13
圖2.1-12 民國93年秋季高雄測站風瑰圖…………..……...……2-14
圖2.1-13 民國93年冬季高雄測站風瑰圖………………..…...…2-14
圖3.1-1 PM2.5採樣設備示意圖……………………………...……3-2
圖3.2-1 高雄市四空品測站分佈………………………….….…3-10
圖3.3-1 高雄市模擬區域範圍示意圖…………………….….…3-15
圖4.1-1 盒型圖所代表的意義………………….……….….….…4-1
圖4.1-2 民國86-93年高雄市空品測站PM10濃度逐時變化趨勢.4-3
圖4.1-3 民國86-93年高雄市空品測站PM10濃度逐日變化趨勢4-4
圖4.1-4 民國86-93年高雄市空品測站PM10濃度逐季變化趨勢.4-5
圖4.1-5 楠梓測站民國86-93年懸浮微粒月趨勢分析圖…..……4-9
圖4.1-6 楠梓測站民國86-93年懸浮微粒月趨勢分析實測值與模擬值比較圖………………………………………….…4-9
圖4.1-7 左營測站民國86-93年懸浮微粒月趨勢分析圖………4-10
圖4.1-8 左營測站民國86-93年懸浮微粒月趨勢分析實測值與模擬值比較圖…………………………………………...4-10
圖4.1-9 前金測站民國86-93年懸浮微粒月趨勢分析圖.……...4-11
圖4.1-10 前金測站民國86-93年懸浮微粒月趨勢分析實測值與模擬值比較圖……………………………………….......4-11
圖4.1-11 小港測站民國86-93年懸浮微粒月趨勢分析圖……...4-12
圖4.1-12 小港測站民國86-93年懸浮微粒月趨勢分析實測值與模擬值比較圖…………………………………………...4-12
圖4.3-1 民國93年小港地區PM2.5中之三種無機鹽實測濃度..4-18
圖4.3-2 民國93年小港站PM10推估值與實測值之比較…..…4-18
圖4.3-3 小港站PM10推估值與實測值之相關性比較……….…4-18
圖4.3-4 民國89年3月25日10時PM10事件日 NO3-等濃度線(μg/m3)分佈圖…….……………………………..……4-19
圖4.3-5 民國89年3月25日10時PM10事件日SO42-等濃度線(μg/m3)分佈圖……….……………………………..…4-20
圖4.3-6 民國89年3月25日10時PM10事件日NH4+等濃度線(μg/m3)分佈圖……….……………………………..…4-20
圖4.3-7 民國89年3月24-26日楠梓站懸浮微粒逐時濃度圖..4-21
圖4.3-8 民國89年3月24-26日左營站懸浮微粒逐時濃度圖..4-21
圖4.3-9 民國89年3月24-26日前金站懸浮微粒逐時濃度圖..4-22
圖4.3-10 民國89年3月24-26日小港站懸浮微粒逐時濃度圖..4-22
圖4.3-11 民國89年3月24-26日PM10模擬值與實測值比較….4-23
圖4.3-12 民國89年3月24-26日小港站SO2逐時濃度比較…...4-23
圖4.3-13 民國89年3月24-26日小港站NOx逐時濃度比較…...4-24
圖4.3-14 高雄市四空品測站89年3月24-26日350 m逆軌跡圖.4-27
圖4.3-15 高雄市四空品測站89年3月24-26日600 m高空逆軌跡圖……………………………………...………………4-28

Andersen Instruments, Inc.,1999. RAASTM OPERATOR’S MANUAL.
A., Thuillier, R. H. and Magliano, K., 1996, Descriptive analysis of PM2.5 and PM10 at Regionally Representative Locations during SJVAQS/AUSPEX , Atmos. Environ. 30, 2079-2112.
Batterman, S. A., 1988. Development of crustal profiles for receptor modeling. Atmos. Environ. 22, 1821-1828.
Brook, J. R., Dann, T. F., and Burnett, R. T., 1997. The relationship among TSP, PM10 , PM2.5 , and inorganic consitituents of atmospheric particulate matter at multiple Canadian locations. J. Air ＆ Waste Manage. Assoc. 47, 2-19.
CAMx-4.0, 2004. User’s Guide for Comprehensive Air Quality Model with Extensions (CAMx), Version 4.0, Environmental International Corp., Novato, CA.
Chaloulakou, A., Kassomenos, P., Spyrellis, N., Demokritou, P., Koutrakis, P., 2003. Measurements of PM10 and PM2.5 particle concentrations in Athens, Greece. Atmos. Environ. 37, 649-660.
Chang, M. E., Cardelino, C., 2000. Application of the Urban Airshed Model to Forecasting Next-Day Peak Ozone Concemtrations in Atlanta, Georgia. J. Air & Waste Manage. Assoc. 50, 2010-2024.
Charlson, R. J., Wigley, T. M. L.,1994. Sulfate Aerosol and Climatic Change. Scientific American, No. 2, 48-57.
Carmichael, G. R., 1996. Seasonal Variation of Aerosol Composition at Cheju, Korea. Atmos. Environ. 30, 2407-2416.
Chen, K. S., Lin, C. F., Chou, Y. M., 2001. Determination of Source Contributions to Ambient PM2.5 in Kaohsiung City Using a Receptor Model, J. Air & Waste Manage. Assoc. 51, 489-498.
Chu, S. H.,1997. Meteorological Conditions Conducive to Regional High Particulate Matter Episodes. Air and Waste Management Association’s 90th Annual Meeting and Exhibition, June 8-13, Toronto, Ontario, Canada 97-MP112.06.
Cope, M., Hess, D., Lee, S., Azzi, M., Carras, J., Wong, N. and Young, M.,1999. Development of the Australian Air Quality Forecasting System : Current Status. International Conference on Urban Climatology, 6.
Cryer, J. D., 1986. Time Series Analysis. Duxbury Press, Boston.
Devore, J. L., 1991. Probability and Statistics for Engineers and the Sciences. Prentice Hall, Englewood Cliffs, New Jersey, USA.
EPA, 1990. User’s Guide for the Urban Airshed Model-Volume I; User’s Manual for UAM (CB-IV). EPA-450/4-90-007A, U.S. Environment Protection Agency, Research Triangle Park, NC.
Gehring, R. and Buchmann, B., 2003. Characterising Seasonal Variations and Spatial Distribution of Ambient PM10 and PM2.5 Concentrations Based on Long-Term Swiss Monitoring Data, Atmos. Environ. 37, 1135-1147.
Giovannoni, J. M., 1993. Modeling of SO2, Pb and Cd Atmospheric Deposition over a One-year Period. Atmos. Environ. 27, 1793-1808.
Hoek, G. and Brunekreff, B., 1995. Effect of photochemical air pollution on acute respiratory symptoms in children. Am. Rev. Respir. Dis. 151, 27-32.
Holland, D. M., Principe, P. P., Sickles, J. E., 1999. Trends in Atmospheric Sulfur and Nitrogen Species in the Eastern United States for 1989-1995. Atmos. Environ. 33, 37-47.
Holton, J. R., 1992. An Introduction to Dynamic Meteorology. 3rd ed., Academic Press, New York, N. Y., USA.
Hoppel, W. A., Fitzgerald, J. W., Frick, G. M., and Larson, R. E., 1990. Aerosol Size Distributions and Optical Properties Found in the Marine Boundary Layer Over the Atlantic Ocean. Journal of Geophysical Research, 95, No.D4, 3659-3686.
Hurley, P. J., Blockley, A. and Rayner, K., 2001. Verification of a Prognostic Meteorological and Air Pollution Model for Year-Long Predictions in the Kwinana Industrial Region of Western Australia, Atmos. Environ. 35, 1871-1880.
IMSL, 1999. IMSL Math/Library, Vol. 1, Visual Numerics, Taipei, Taiwan.

Katada, T., Carmichael, G. R., Peters, L. K., 1984, Numerical Simulation of the Transport of Chemically Reactive Species under Land and Sea-Breeze Circulations. J. Climate Appl. Meteor. 23, 1153-1172.
Kim, J. Y., Ghim, Y. S., Kim, Y. P. and Dabdub, D.,2000. Determination of domainfor diagnostic wind field estimation in Korea. Atmos. Environ. 34, 595-601.
Kitada, T. and Kitagawa, E.,1990. Numerical Analysis of the Role of Sea Breeze Fronts on Air Quality in Coastal and Inland Polluted Areas. Atmos. Environ. 24A, 1545-1559.
Kuang, Y.C., 2001. Simulation and Sensitivity Analysis of The Secondary Aerosols over Taiwan Island, Journal of the Chinese Institute of Environmental Engineering 11, 97-104.
Lam, G. C. K., Leung, D. Y. C., Niewiadomski, M., Pang, S. W., Lee, A. W. F. and Louie, P. K. K., 1999. Street-level concentrations of nitrogen dioxide and suspended particulate matter in Hong Kong. Atmos. Environ. 33, 1-11.
Lin, J.J. and Tai, H.S., 2001. Concentration and Distributions of Carbonaceous Species in Ambient Particle in Kaohsiung City, Taiwan. Atmos. Environ. 35, 2627-2636.
Lin, J. J., 2002. Characterization of the Major Chemical Species in PM2.5 in the Kaohsiung City, Taiwan, Atmos. Environ. 36, 1911-1920.
Lin, T. S., 2001. Long-Range Transport of Yellow Sand to Taiwan in Spring 2000 : Observed Evidence and Simulation. Atmos. Environ. 35, 5873-5882.
Lu, H. C. and Fang, G. C., 2002. Estimating the Frequency Distributions of PM10 and PM2.5 by the Statistics of Wind Speed at Sha-Lu, Taiwan. The Sci. of the Total Environ. 298, 119-130.
Luhar, A. and Hurley, P.J., 2003. Evaluation of TAPM, a Prognostic Meteorological and Air Pollution Model, Using Urban and Rural Point-Source Data, Atmos. Environ. 37, 2795-2810.
Mamane, Y., 1986. Estimate of Municipal Refise Incinerator Contribution to Philadelphia Aerosol Using Single Particle Analysis. Atmos. Environ. 24B, 127-135.
Marcazzan, G. M., Valli G. and Vecchi, R., 2002. Factors Influencing Mass Concentration and Chemical Composition of Fine Aerosols during a PM High Pollutant Episode, The Sci. of the Total Environ. 298, 65-79.
Melas, D., Ziomas, I., Klemm, O. and Zerefos, C. S.,1998. Flow Dynamics in Athens Area under Moderate Large-scale Winds. Atmos. Environ. 32, 2202-2222.
Milfrod, J. B., Gao, D., Sillman, S., Blosse, P., Russell, A. G., 1994. Total Reactivity Nitrogen (NOy) as an Indicator for the Sensitivity of Ozone to NOx and Hydrocarbons. J. Geophys. Res. 99, 3533-3542.
Nester, K., 1995. Influence of Sea Breeze Flows on Air Pollution over theAttica Peninsula. Atmos. Environ. 29, 3655-3670.
Noll, K. E. and Fang, K. Y. P., 1989. Development of a Dry Deposition Model for Atmospheric Coarse Particles. Atmos. Environ. 23, 585-594.
Ohta, S. and Okita, T., 1990. A Chemical Characterization of Atmospheric Aerosol in Sapporo. Atmos. Environ. 24A, No. 4, 815-822.
Park, I.S., Lee, S.J, Kim, C.H., Yoo, C. and Lee, Y.H., 2004. Simulation urban-scale air pollutions and their predicting capabilities over the Seoul metropolitan area. J. Air & Waste Manage. Assoc. 54, 659-710.
Pisano, J. T., Mckendry, I., Steyn, D. G., and Hastie, D. R., 1997. Vertical nitrogen dioxide and ozone concentrations measured from a tethered balloon in the lower Fraser valley. Atmos. Environ. 31, 2071-2078.
Schwartz, J., Dockery, D. W. and Neas, L. M.,1996. Is Daily Mortality Associated Specifically With Fine Particles?. J. Air & Waste Manage. Assoc. 46, 927-939.
Shrestha, A. B., Wake, C. P. and Dibb, J. E., 1997. Chemical Composition of Aerosol and Snow in the High Himalaya During the Summer Monsoon Season. Atmos. Environ. 31, 2815-2826.
Slinn, W. G. N., 1982. Predictions for Particle Deposition to Vegetative Canopies. Atmos. Environ. 16, 1785-1794.
Slinn, S. A., Slinn, W. G. N., 1980. Predictions for Particle Deposition on Natural Waters. Atmos. Environ. 14, 1013-1016.
Stout, J. E., 2001. Dust and Environment in the Southern High Plains of North America. J. of Arid Environ. 47, 425-441.
Suppan, P., Fabian, P., Vyras, L. and Gryning, S. E.,1998. The Bahaviour of Ozone and Peroxyacetyl Nitrate Concentrations for Wind Regimes during the Medeaphot- trace Campaign in the Greater Area of Athens, Gress. Atmos. Environ. 32, 2098-2102.
Tang, I. N., Munkelwitz, H. R. and Davis, J. G., 1977. Aerosol Growth Studies Ⅱ：Preparation and Growth Measurements of Monodisperse Salt Aerosols. J. of Aerosol Sci., 8, 149-159.
Triantafyllou, A. G., 2001. PM10 pollution episodes as a function of synoptic climatology in a mountains industrial area. Environ. Pollution, 112, 491-500.
Vega, E., Reyea, E., Ruiz, H., Gabriela, J., Sanchez, G., Gonzaiez U., 2004. Analysis of PM2.5 and PM10 in the Atmosphere of Mexico City during 2000-2002. Air ＆ Waste Manage. Assoc. 54, 786-798.
Vijayakumar, G., Parameswaran, K. and Rajan, R., 1998. Aerosols in the Atmospheric Boundary Layer and its Association with Surface Wind Speed at a Coastal Site. J. of Atmos. and Solar-terrestrial Physics, 60, 1531-1542.
Waggoner, A. P., Weiss, R. E., Ahlquist, N. C., Covert, D. S., Will, S. and Charlson, R. J.,1981. Optical Characteristics of Atmospheric Aerosols. Atmos. Environ. 15, 2617-2624.
Watson, J. G., Chow, J. C., Lu, Z., Fujita, E. M., Lowenthal, D. H., and Lawson, D. R., 1994. Chemical Mass Balance Source Apportionment of PM10 during the Southern California Air Quality Study. Aerosol Sci. Technol. 21, 1-36.
Wedepohl, K. H., 1987. The Composition of the Upper Earth’s Crust and the Natural Cycles of Selected Metals. Metals and their Compounds in the Environment.
Whitby, K. T., Husar, R. B., and Liu, B. Y. H., 1972. The Aerosol Size Distribution of Los Angeles Smog. Aerosols and Atmospheric Chemistry, G. M. Hidy, Ed., Academic Press, New York, 237-264.
Willmott, C. J., Ackleson, S. G., Davis, R. E., Feddema, J. J., Klink, K. M., Legates, D. R., O’Donnell, J. Rowe, C. M., 1985. Statistics for the Evaluation and Comparisons of Model. J. Geophys. Res. 90, 8995-9005.
Yamada, T., Mellor, G. L., 1975. A Simulation of the Wangara Atmospheric Boundary Layer Data. J. Atmos. Sci. 32, 2309-2329.
王秋森，1994，「石化工廠產的粒狀空氣污染物受體模式之建立」，國科會研究報告。
TEDS5.0, 2002，Taiwan Emission Data System (TEDS) 5.0，中鼎工程公司。
何宜達，2004，「高屏地區臭氧事件日光化學模式解析及氣象條件之探討」，國立中山大學環工所博士論文。
李俊璋，1992，「台北市空氣中懸浮微粒物理化學分析及學童肺功能之研究」，國立台灣大學環境工程研究所碩士論文。
李清勝，1998，「導致臺灣地區懸浮微粒高污染之氣象分析與預報」，行政院環境保護署研究計畫報告，EPA-87-E3L1-03-01。
余嘉裕、卓盈旻、涂建翊，2002，「東亞沙塵暴的時空特徵分析」，環境保護，25，178-198.
吳啟文，1996，「台灣中部都會區氣懸微粒粒徑分佈之污染物特性分析」，國立中央大學環境工程研究所碩士論文。
吳家興，1995，「台灣北部地區國中學生氣喘的研究-空氣污染的影響」，國立台灣大學公共衛生學研究所碩士論文。
林鉅富，1999，「高雄地區大氣中懸浮微粒PM2.5特性及來源之探討」，國立中山大學環境工程研究所碩士論文。
林能暉、黃景祥、彭啟明、胡翠華、方淑慧、李貞瑩，2002，「沙塵暴事件之空氣品質統計數值模式之建立及應用」，環境保護，25，117-133。
林清和、施百鴻、林啟燦、丁杉龍，1999，「南高屏地區混合層內平均風場之求取」，第十六屆空氣污染控制技術研討會論文集，665-670。
林能暉，2000，「北中南空品區O3 與PM10 污染之氣象條件研究」，國科會/環保署科技合作研究計劃報告，NSC 89-EPA-Z-008-007。
林鉅富，1999，「高雄地區大氣中懸浮微粒PM2.5特性及來源之探討」，國立中山大學環境工程研究所碩士論文。
柳中明、楊之遠，2002，「影響台灣空氣品質的大陸沙塵天氣：2002年個案研究」，環境保護，25，153-177。
柳中明、蘇維中，1997，「區域氣象環境與高臭氧之相關分析」，大氣科學25，27-49。
柳中明、蘇維中，1997，「區域氣象環境與高臭氧之相關分析」，大氣科學，25，27-48。
袁中新、張瑞正、袁菁、楊宏宇、林文印、李崇垓、李崇德，1999，「能見度與懸浮微粒物化特徵之相關性探討」，氣膠研討會論文集，67-75。
張順欽，2002，「亞洲沙塵暴對空氣品質影響與測報」，環境保護，25，134-152。
許文昌、李崇德，1992，「台灣都會區氣懸微粒污染來源之推估」，第九屆空氣污染控制技術研討會論文集，189-202。
陳順宇，2000，「迴歸分析」三版，華泰書局。
陳康興、周明顯、王文正、彭彥彬、管琪芬、黃耀田、楊峰杰、廖琇怡，2005，「九三年度空氣品質監測站污染源來源分析計畫」，高雄市政府環境保護局。
陳凱爾，2001，EXCEL 在統計學上的應用，五南圖書出版有限公司。
張順欽、吳權芳，2001，「大陸沙塵暴對台彎地區空氣品質之影響」，大陸沙塵暴對台灣地區空氣品質之影響與預測研討會。
黃宗正、李正綱、曾錦富等，1992，「台中發電廠南方空氣中懸浮微粒特性研究」，第九屆空氣污染控制技術研討會論文集。
黃明雄、李崇德、林立偉，1998，「墾丁地區氣膠污染來源推估及在不同天氣類型氣流軌跡的傳輸」，第十五屆空氣污染控制技術研討會論文集，595-602。
望熙榮、張智泳，1999，「中南部地區空氣品質監測站代表性評估」，第十六屆空氣污染控制技術研討會論文集，609-614。
莊秉潔，1993，「地表層擴散係數、風速、雨量、及混合層對台北都會區早上八點懸浮為例濃度之影響」，中國環境工程學刊，第三卷，第二期，95-102。
莊秉潔、黃沛霖、陳建隆、林忠賢、林宏乾、蔡徵霖，2000，「軌跡模式加入地形參數之影響研究」，第十七屆空氣污染控制技術研討會論文集，473-477。
程萬里，2000，「中部空品區空氣污染事件之天氣型態暨氣象資料庫整合建立」，環保署研究計劃報告。
楊忠盛，1998，「台北都會區懸浮微粒特性及來源之探討」，國立台灣大學環境工程研究所碩士論文。
楊定國，1999，「台中地區都會區、近郊區、郊區日夜間懸浮微粒、金屬離子及酸性陰離子之研究」，東海大學環境科學研究所碩士論文。
蔣本基等，1996，「台灣地區懸浮微粒空氣污染問題及防制之研究」，環保署研究報告。
蔡瀛逸、鄭曼婷，1998，「高污染狀態下大氣二次氣膠組成之探討」，第十五屆空氣污染控制技術研討會，711-718。
蔡瀛逸，1992，「大氣懸浮微粒中鉛元素之研究」，第九屆空氣污染控制技術研討會論文集，245-256。
蕭慧娟，1996，「空氣污染概論」，空氣污染防治專責人員訓練教材。
環保暑監資處，1998，「台灣地區細懸浮微粒(PM2.5)監測初步分析報告」。