大氣污染物干沉降速度和通量的計算方法比較一以南京仙林地區(qū)為例趙雄飛;王體健;黃滿堂;楊帆;楊丹丹;金龍山【摘要】目的基于不同方法對大氣污染物干沉降速度和通量的估計存在差異，開展比較研究.方法2016年9月至2017年9月,在南京大學(xué)仙林校區(qū)，基于75m觀測塔，對大氣中常見的六種污染物二氧化硫(SO2)、一氧化氮(CO)、二氧化氮(NO2)、臭氧(O3)、一氧化碳(CO)、細(xì)顆粒物(PM2.5)的濃度和氣象要素進行連續(xù)觀測.利用三層阻力模型計算大氣污染物的干沉降速度，利用濃度法和梯度法計算干沉降通量并對兩種方法進行比較結(jié)果SO2、NO、NO2、O3、CO、PM2.5的平均干沉降速度分別是0.270、0.019、0.089、0.449、0.038、0.147cm/s.干沉降速度具有明顯的日變化特征，一般情況下，白天大于夜間，在午后出現(xiàn)最大值.整個觀測期間，采用濃度法計算得到的SO2、NO、NO2、O3、CO、PM2.5干沉降通量分別為0.034、0.008、0.037、0.263、0.354、0.049pg/(m2-s),采用梯度法得到的干沉降通量分別為0.04、0.00193、0.035、0.278、0.192、0.063pg/(m2-s)結(jié)論對于NO、O3、PM2.5,濃度法和梯度法計算的干沉降通量具有較好的一致性.梯度法估計干沉降通量時很大程度上依賴于大氣污染物濃度梯度測量的準(zhǔn)確性，濃度法估計干沉降通量則更多依賴于干沉降速度計算的準(zhǔn)確性.【期刊名稱】《裝備環(huán)境工程》【年(卷)，期】2019(016)006【總頁數(shù)】9頁(P129-137)【關(guān)鍵詞】干沉降通量;干沉降速度;濃度法;梯度法【作者】趙雄飛;王體健;黃滿堂;楊帆;楊丹丹;金龍山【作者單位】南京大學(xué)大氣科學(xué)學(xué)院，南京210023;南京大學(xué)大氣科學(xué)學(xué)院，南京210023;南京大學(xué)大氣科學(xué)學(xué)院,南京210023;南京大學(xué)大氣科學(xué)學(xué)院,南京210023;南京大學(xué)大氣科學(xué)學(xué)院,南京210023;南京大學(xué)大氣科學(xué)學(xué)院,南京210023【正文語種】中文【中圖分類】X511物質(zhì)在大氣和下墊面之間的交換過程中，大氣沉降過程扮演著非常重要的角色。大氣沉降過程主要分為干沉降與濕沉降過程。干沉降是大氣中污染物清除的主要過程之一，直接影響空氣中污染物的濃度和空氣質(zhì)量狀況的時空分布［1］。大氣污染物的干沉降通量已經(jīng)在世界范圍內(nèi)得到了廣泛的研究，污染物沉降到地表的過程不僅與酸雨污染有著密不可分的聯(lián)系，同時也對下墊面的生態(tài)環(huán)境產(chǎn)生極大的影響。目前，發(fā)達(dá)國家的大氣污染物沉降觀測繼續(xù)朝著網(wǎng)絡(luò)化、系統(tǒng)化方向不斷發(fā)展，并開發(fā)和應(yīng)用了適用于不同空間尺度和時間精度需要的污染物排放、傳輸、沉降模型［2-7］，利用數(shù)值模型模擬不同生態(tài)系統(tǒng)的污染物沉降總量和沉降負(fù)荷，進而為相應(yīng)的污染物減量策略的制定提供科學(xué)依據(jù)。影響干沉降過程的因素主要分為氣象要素和污染物本身在大氣中的濃度水平。其中影響干沉降過程的主要氣象因素有空氣動力學(xué)特征量（摩擦速度、下墊面粗糙度長等）、流場分布、大氣穩(wěn)定度、風(fēng)、溫度、濕度的平均場和湍流場結(jié)構(gòu)參數(shù)、逆溫層高度等［22］。影響干沉降過程的因素還有污染物自身性質(zhì)。對于氣體而言，主要是化學(xué)活潑程度、溶解度、與沉降表面處于平衡狀態(tài)時的分壓強等主要因素；對于顆粒物污染物而言，則主要是粒子直徑、密度、形狀、狀態(tài)、表面性質(zhì)、化學(xué)成分等主要因素。總體而言，無論是對于氣體還是顆粒物，污染物本身在大氣中的濃度是決定其干沉降通量的最主要因素。對于不同的下墊面情況，也存在有不同的干沉降過程，例如對于下墊面中含有植被的干沉降過程，植物結(jié)構(gòu)（如葉子空間密度、葉面面積與形狀、葉子多孔性、葉孔阻力、葉子表面形狀等）、植被生長狀況、植被表皮葉質(zhì)滲出等情況、植被的軟硬以及其所含污染物的本底情況也是影響干沉降通量的主要因素［8-15］。大氣污染物干沉降通量測量方法主要有渦動相關(guān)法［16-17］、梯度法［18-20］以及濃度法［21］。大氣干沉降過程的效率受氣候背景條件、邊界層湍流狀況、下墊面植被類型等多種因素的影響［21］。文中以江蘇省南京市仙林地區(qū)為例，結(jié)合鐵塔觀測數(shù)據(jù)與數(shù)值模擬來研究大氣中污染物的干沉降過程。主要研究內(nèi)容包括采用三層阻力模型來計算不同種類污染物的干沉降速度，進而計算其干沉降通量，同時利用氣象資料計算大氣垂直湍流擴散系數(shù)，結(jié)合污染物濃度梯度進而計算其干沉降通量。將兩種方法計算所得的干沉降通量結(jié)果進行比較，并分析不同方法的優(yōu)劣。1研究區(qū)域與方法1.1觀測點情況觀測地點位于江蘇省南京市東北方向的南京大學(xué)仙林校區(qū)地球系統(tǒng)區(qū)域過程綜合觀測試驗基地75m氣象塔處（東經(jīng)118°，北緯32°）。氣象塔位于海拔60m的小山上，下墊面具有城郊混合特征，其正下方為平坦、裸露的標(biāo)準(zhǔn)氣象觀測場，山坡種植有多種樹木，山頂種有草本植物。圖1觀測點位置以及觀測裝置1.2氣象參數(shù)觀測方法利用中國華云集團生產(chǎn)的CAWS-600國家級氣象自動觀測站同時進行風(fēng)速及地表溫度數(shù)據(jù)的測量。CAWS-600系統(tǒng)采用的集采核芯具有很好的穩(wěn)定性，數(shù)據(jù)具有高精度、高可靠性和高準(zhǔn)確性。連續(xù)收集風(fēng)速以及地表溫度的同時，在氣象觀測塔上的25、50、75m三層分別安裝傳感器連續(xù)收集實時的溫濕度等數(shù)據(jù)。1.3污染物濃度觀測方法在氣象塔25、50、75m三層高度各放置了一臺MAS-100小型空氣站，對大氣中六種污染物SO2、NO、NO2、O3、CO、PM2.5的濃度進行連續(xù)監(jiān)測，并將收集到的濃度數(shù)據(jù)同時上傳到遠(yuǎn)程數(shù)據(jù)平臺。MAS-100型小型空氣站由深圳創(chuàng)龍清研科技有限公司制造，該型儀器可實時收集多種大氣污染物和環(huán)境參數(shù)，實現(xiàn)高時空分辨率的數(shù)據(jù)采集，并通過4G信號遠(yuǎn)程傳輸，實現(xiàn)云端數(shù)據(jù)庫儲存或與智能手機連接。采用高靈敏度電化學(xué)傳感器和顆粒物傳感器，利用電化學(xué)反應(yīng)分辨氣體成分，結(jié)合嵌入式技術(shù)和網(wǎng)絡(luò)通訊技術(shù)，可以連續(xù)監(jiān)測大氣中空氣污染物濃度和溫濕度等參數(shù)。1.4干沉降通量計算1.4.1濃度法利用污染物濃度觀測和數(shù)學(xué)模型相結(jié)合的辦法進行干沉降通量的計算，即采用三層阻力模型計算氣體與氣溶膠粒子的干沉降速率，與觀測到污染物的濃度分別相乘，即可得到相對應(yīng)的干沉降通量。對2016年9月1日到2017年9月1日之間采集到的數(shù)據(jù)進行計算，得到一年時間范圍內(nèi)的干沉降速度，從而計算干沉降通量。對于氣體和顆粒物而言，其干沉降通量的計算方法均為：式中：F為氣體或顆粒物的干沉降通量外g/(m2?s);C為50m高度處的污染物濃度，口g/m3;vd代表其干沉降速度，cm/s。使用濃度法進行干沉降通量計算時，其重點是干沉降速度的計算。干沉降速度的不同主要與當(dāng)?shù)叵聣|面環(huán)境以及氣象條件有關(guān)。對于氣體污染物和顆粒物污染物而言,干沉降速度的計算是不同的，計算方法分別為：式中：Ra、Rb、Rc分別為空氣動力學(xué)阻力、片流層阻力、近地面阻力。Ra是因為大氣湍流運動所引起的污染物從大氣中輸送至地表附近過程中所受到的阻力，根據(jù)近地層中污染物的質(zhì)量輸送以及類似熱量輸送的假設(shè)，Ra可以表示為：式中：z為計算干沉降速度時所選取的參考高度，m；u*是摩擦系數(shù)，m/s；z0為下墊面的粗糙度長，m；k為常數(shù)，通常根據(jù)不同地形所對應(yīng)的地形特征經(jīng)驗獲得，取值為0.4;^c為與質(zhì)量通量有關(guān)的無量綱穩(wěn)定度修正參數(shù)。Wc通常表示為：式中：，L為莫寧-奧布霍夫長度。1.4.2梯度法按照湍流擴散的梯度輸送理論來處理大氣污染物的垂直散布問題，是以垂直湍流擴散系數(shù)及其時空變化來表征擴散能力的。即表示垂直湍流所引起的物質(zhì)輸送的速率，取決于該物質(zhì)分布的不均勻程度（即污染物梯度的大?。┮约爱?dāng)前氣象環(huán)境流場下所具備的擴散能力。在常通量層假設(shè)條件下，根據(jù)湍流輸送理論，污染物輸送通量［12］可以表示為：式中：K為垂直湍流擴散系數(shù)，m2/s。本次試驗中，cd分別選取75m與25m的濃度數(shù)據(jù)進行計算，dz為兩層采樣器之間的高度差50m。在中性大氣中，K可以通過風(fēng)速廓線直接求得，但在實際大氣中，該條件則很難被滿足。因此常使用無量綱的穩(wěn)定度函數(shù)Wm來進行矯正［22］：式中：K為馮卡門常數(shù)，取值范圍一般為0.35~0.41，在文中取值為0.4;u*為摩擦速度，m/s;z為測量濃度污染物所對應(yīng)的高度，m;L為莫寧-奧布霍夫長度，可以利用10m處的湍流資料直接獲得。對于湍流擴散系數(shù)，在大氣最低的50~100m層內(nèi)，可用布辛格等人［23］提出的公式：2結(jié)果與分析2.1污染物濃度變化2016年9月到2017年9月南京仙林地區(qū)觀測點鐵塔50m處小型空氣站監(jiān)測到的氣體濃度的時間序列如圖2所示。計算表明，觀測期間，南京仙林地區(qū)大氣中SO2、NO、NO2和03的平均濃度分別為(12.9±5.2)、(39.7±20.3)、(46.4±17.7)、(67.2±34.3)|Jg/m3,CO和PM2.5的平均濃度分別為(0.93±0.3)mg/m3和(33.61士23.5)pg/m3。在全年污染物濃度變化中，S02在3—6月份出現(xiàn)濃度高值；NO濃度波動不大，與其他城市變化趨勢相似［24］;NO2濃度在冬春季節(jié)較高，其余時間較低；O3濃度則呈現(xiàn)明顯的夏秋高于冬春的特點；CO濃度在夏秋季節(jié)比較高，其余時間較低；PM2.5濃度則在2016年12月和2017年1、2、3、4、7月出現(xiàn)了較高的峰值。圖2仙林地區(qū)大氣污染物濃度變化2.2污染物的干沉降速率采用三層阻力模型方法，計算得到六種大氣污染物SO2、NO、NO2、O3、CO、PM2.5的干沉降速率分別是0.27、0.019、0.089、0.449、0.038、0.147cm/s。由圖3可知，仙林地區(qū)所處的類城郊下墊面區(qū)域，除NO之外，其他五個物種的干沉降速度都有明顯的日變化，即白天大于夜間?？赡茉蛟谟诖髿夥€(wěn)定度和表面阻力的日變化，夜間由于沒有太陽輻射作用，在地面輻射冷卻作用下，近地面氣溫逐漸降低，出現(xiàn)逆溫層，大氣呈現(xiàn)出穩(wěn)定狀態(tài)，導(dǎo)致大氣垂直擴散過程較難進行，使得干沉降速度在夜間呈現(xiàn)較低的態(tài)勢。相反，白天強烈的太陽輻射會使得近地面迅速升溫，導(dǎo)致近地面空氣的溫度迅速增加，形成了氣溫隨高度增加而增加的狀況。在太陽輻射的作用下，大氣由夜間的穩(wěn)定狀態(tài)逐漸轉(zhuǎn)向為不穩(wěn)定狀態(tài)，大氣湍流活動逐漸增強，湍流層阻力較小，使得污染物的干沉降速度出現(xiàn)明顯的增加［25］。由圖3可見，在上午8—9時，以O(shè)3、SO2、PM2.5為代表的污染物，干沉降速度開始增長，在中午到下午的時間段達(dá)到全天中最大值。圖3南京仙林地區(qū)大氣污染物干沉降速度的日變化對于氣體和顆粒物而言，影響其干沉降速度的因素是不同的。對于氣體而言，其干沉降速度不僅與下墊面特征以及氣象要素有關(guān)，同時太陽輻射也會改變植被氣孔開閉狀況以及植被表面溫度和相對濕度，所以粘附層阻力和表面阻力在很大程度上會決定污染物的干沉降速度日變化特征。對于顆粒物來說，粘附層阻力和顆粒物自身的重力沉降在其沉降過程中占有主導(dǎo)地位［26］。不同季節(jié)六種大氣污染物的干沉降速度見表1。從污染物干沉降速度的季節(jié)變化上來看，SO2干沉降速度具有明顯的冬春季較大、夏秋季較小的變化特征。這是由于大氣越穩(wěn)定，湍流運動越弱，湍流通量也就越小，而平均風(fēng)速在近地面常通量層中的垂直切變卻比不穩(wěn)定時強得多，所以SO2干沉降速度與一年四季中的大氣穩(wěn)定度有著對應(yīng)的關(guān)系。在大氣穩(wěn)定時，SO2的干沉降速度小，在大氣處于不穩(wěn)定狀態(tài)時，干沉降速度大。NO由于其亨利系數(shù)非常小，其表面阻力非常大，干沉降速度接近于0，幾乎沒有季節(jié)變化。對于NO2而言，其表面阻力非常大，在污染物干沉降過程中占主導(dǎo)地位，決定著沉降速度的大小和量級。表面阻力的季節(jié)變化也影響著NO的干沉降速度的季節(jié)變化，呈現(xiàn)出春冬較低，夏季較高的變化趨勢。O3和SO2一樣，表面阻力與湍流層阻力和粘附層阻力相當(dāng)，O3干沉降速度的季節(jié)變化是氣象要素和近地面綜合作用的結(jié)果。O3的干沉降速度呈現(xiàn)出夏季較大、冬季較小的變化趨勢。CO的表面阻力非常大，干沉降速度比較低，季節(jié)變化不明顯。PM2.5的干沉降速度在春秋天呈現(xiàn)出一個雙峰的變化趨勢，這是由于顆粒物干沉降速度主要受氣象要素的變化而產(chǎn)生的結(jié)果。表1南京仙林地區(qū)下墊面上干沉降速度月均值cm/sSO2NONO203COPM2.51月0.4910.0190.0770.4270.0380.1582月0.5040.0190.0780.4370.0380.1753月0.1600.0190.0780.4370.0380.1714月0.1670.0200.0800.4700.0380.1455月0.1700.0200.0810.4820.0390.1436月0.1670.0200.0800.4580.0390.1297月0.1720.0200.0810.4960.0390.1398月0.1680.0200.0800.4600.0390.1219月0.1680.0200.0800.4700.0390.13110月0.1660.0200.0800.4640.0380.15011月0.1580.0190.0770.4250.0380.15212月0.4830.0190.0770.4200.0380.1532.3干沉降通量計算結(jié)果結(jié)合以上2016年9月到2017年9月時間段內(nèi)的NO、NO2、SO2、O3、CO、PM2.5的濃度，以及分別對應(yīng)的小時干沉降速度，計算得到濃度法干沉降通量。六種污染物干沉降通量平均值分別為SO20.034pg/(m2-s).NO0.008|jg/(m2?s)、NO20.037|jg/(m2?s)、O30.263|jg/(m2?s)、CO0.354|jg/(m2?s)、PM2.50.049pg/(m2-s)o為保證干沉降結(jié)果的有效性，對觀測到的數(shù)據(jù)進行了篩選和質(zhì)量控制，以剔除地面排放大于沉降的情況。其主要方法為對污染物濃度在75m以及25m的濃度作差，所有75m小于25m濃度的數(shù)據(jù)均被剔除。利用地面的氣象觀測數(shù)據(jù)計算垂直湍流擴散系數(shù)，結(jié)合濃度梯度，從而求得干沉降通量。對于本次鐵塔觀測的六種污染物來說，其垂直湍流擴散系數(shù)與污染物種類無關(guān)，均由氣象條件所決定。六種污染物梯度法干沉降通量分別為：SO20.04|jg/(m2?s)、NO0.00193pg/(m2-s).NO20.035pg/(m2-s).O30.278pg/(m2?s)、CO0.192pg/(m2?s)、PM2.50.063pg/(m2-s)o將文中結(jié)果與文獻報道結(jié)果相比，三峽庫區(qū)2015年采用濃度法計算得到的NO2干沉降通量為0.010Mg/(m2-s),實驗中計算所得的NO2干沉降通量為0.035pg/(m2-s)o由于本實驗位于東部地區(qū)，相對于三峽庫區(qū)NO2濃度為(11.7±3.9)pg/m3，南京地區(qū)的NO2年均濃度為(46.4±17.7川g/m3，所以相對應(yīng)的NO2的沉降通量也較大［27］。環(huán)渤海地區(qū)2012年采用濃度法計算的SO2干沉降通量為0.05~0.25pg/(m2-s)o與文中的SO2干沉降通量0.034pg/(m2-s)相比，文中SO2干沉降通量較小，其原因可能在于近年來中國地區(qū)的SO2減排措施導(dǎo)致其濃度下降，從而導(dǎo)致干沉降通量的下降［28］。對于梯度法計算得到的干沉降通量，南京江北地區(qū)農(nóng)田下墊面在農(nóng)作物生長旺盛期間的O3干沉降通量均值為0.35Mg/(m2-s)，與文中O3全年干沉降通量0.278pg/(m2-s)具有較好的一致性［29］。2.4兩種方法的比較通過濃度法和梯度法計算的干沉降通量平均后的月均變化趨勢如圖4所示，利用兩種方法計算得到的六種污染物SO2、NO、NO2、O3、CO、PM2.5的平均干沉降通量分別為0.037,0.005,0.036,0.271,0.273,0.056pg/(m2-s)o通過濃度法和梯度法計算的六種大氣污染物的干沉降通量結(jié)果見表2,可以看出，NO和CO的干沉降通量計算結(jié)果的偏差相對較大。圖4仙林地區(qū)大氣污染物月均干沉降通量分布表2南京仙林地區(qū)大氣污染物干沉降通量pg/(m2-s)干沉降通量SO2NONO2O3COPM2.5濃度法0.0340.0080.0370.2630.3540.049梯度法0.0410.0020.0350.2780.1920.063絕對偏差-0.0070.0060.002-0.0150.162-0.014相對偏差/%-20.5875.005.41-5.7045.76-27.59進一步選取兩種方法計算結(jié)果均有效的時次進行對比分析，以六種大氣污染物的濃度法干沉降通量為橫軸，梯度法干沉降通量為縱軸，對NO、NO2、SO2、O3、CO、PM2.5六種大氣污染物進行兩種方法計算結(jié)果相關(guān)性分析，結(jié)果如圖5所示。計算得到六種污染物的濃度法干沉降通量-梯度法干沉降通量的相關(guān)系數(shù)分別為0.438、0.020、0.015、0.225、0.101、0.392。文中分別采用濃度法與梯度法對大氣中六種污染物SO2、NO、NO2、03、CO、PM2.5的干沉降通量進行了計算，對于兩種方法的結(jié)果進行對比，其偏差可能來源于以下幾個方面。1）由于對六種污染物的濃度是采用高靈敏傳感器測量得到的，相對傳統(tǒng)的點式測量儀器，傳感器本身存在一定的測量誤差。2）利用濃度法計算干沉降通量時，由于鐵塔觀測點的下墊面具有城郊混合的特征，在計算干沉降速度時，近地面阻力的參數(shù)化存在一定的誤差。3）利用梯度法計算干沉降通量時，一方面由于傳感器測量污染物濃度本身存在一定的誤差，另一方面由于鐵塔上只有三層觀測高度，兩層之間高度差最大只有50m，相對較小的濃度梯度使得計算的干沉降通量存在一定的誤差。3結(jié)論筆者在南京仙林地區(qū)進行了為期1年的大氣污染物的干沉降觀測，利用三層阻力模型計算干沉降速度，結(jié)合濃度法和梯度法估算大氣中常見的六種污染物SO2、NO、NO2、O3、CO、PM2.5的干沉降通量，得到如下主要結(jié)論。1）在南京仙林地區(qū)城郊結(jié)合下墊面，大氣污染物干沉降速度分別為：CO0.038cm/s、SO20.27cm/s、NO0.019cm/s、NO20.089cm/s、O30.449cm/s、PM2.50.147cm/s。干沉降速度具有明顯的日變化特征，一般情況下白天大于夜間，在午后出現(xiàn)最大值。圖5濃度法、梯度法沉降通量對比2）利用濃度法和梯度法計算得到六種污染物SO2、NO、NO2、O3、CO、PM2.5的平均干沉降通量分別為0.037、0.005、0.036、0.271、0.273、0.056Mg/（m2-s）o兩種方法對NO和CO干沉降通量計算結(jié)果的偏差相對較大。3）綜合濃度法和梯度法對干沉降通量的計算結(jié)果，發(fā)現(xiàn)對于NO、O3和PM2.5,兩者的相關(guān)性比較好，其相關(guān)系數(shù)分別為0.438、0.225、0.392;對于NO2、SO2、CO,兩者的相關(guān)性比較差，相關(guān)系數(shù)分別為0.020、0.015、0.101。在進行干沉降通量計算的時候，應(yīng)該根據(jù)實際儀器設(shè)備的情況對不同方法進行合理的選擇。參考文獻：【相關(guān)文獻】蘇航，王自發(fā)，朱彬，等.北京奧運會前后靜穩(wěn)天氣條件下SO2和NO2干沉降模擬[J].氣候與環(huán)境研究，2010，15(5):636-642.APSIMONHM，BARKERBM，KAYINS.ModellingStudiesoftheAtmosphericReleaseandTransportofAmmoniainAnticyclonicEpisodes[J].AtmosphericEnvironment1994，28(4):665-678.SAKURAIT，F(xiàn)UJITASI，HAYAMIH，etal.ACaseStudyofHighAmmoniaConcentrationintheNighttimebyMeansofModelingAnalysisintheKantoRegionofJapan[J].AtmosphericEnvironment，2003，37(31):4461-4465.DOREAJ，THEOBALDMR，KRYZAM，etal.ModellingtheDepositionofReducedNitrogenatDifferentScalesintheUnitedKingdom[C]//Proceedingofthe29thNATO/CCMSInternationalTechnicalMeetingonAirPollutionModelingandItsApplication.Aveiro，Portugal，2007.HELLSTENS，DRAGOSITSU，PLACECJ，etal.ModellingtheSpatialDistributionofAmmoniaEmissionsintheUK[J].EnvironmentalPollution，2008，154(3):370-379.WALKERJ，SPENCEP，KIMBROUGHS，etal.InferentialModelEstimatesofAmmoniaDryDepositionintheVicinityofaSwineProductionFacility[J].AtmosphericEnvironment，2008，42(14):3407-3418.KRYZAM,DOREAJ,BA,AA>M,etal.ModellingDepositionandAirConcentrationofReducedNitrogeninPolandandSensitivitytoVariabilityinAnnualMeteorology[J].JournalofEnvironmentalManagement,2011,92(4):1225-1236.WESELYML,COOKDR,HARTRL,etal.SpeerMeasurementsandParametrizationofArticleSulfurDepositionoverGrass[J]JGeophysRes,1985,90:2131-2143GIORGIF.AParticleDry-depositionParameterizationSchemeforUseinTracerTransportModels[J].JournalofGeophysicalResearchAtmospheres,1986,91(D9):9794-9806.PADROJ，DENHARTOGG，NEUMANNHH.AnInvestigationoftheADOMDryDepositionModuleUsingSummertimeO3MeasurementsAboveaDeciduousForest[J].AtmosphericEnvironmentPartAGeneralTopics，1991，25(8):1689-1704.ERISMANJW,VANPULA,WYERSP.ParametrizationofSurfaceResistancefortheQuantificationofAtmosphericDepositionofAcidifyingPollutantsandOzone[J].AtmosphericEnvironment,1994,28(16):2595-2607.PADROJ.SummaryofOzoneDryDepositionVelocityMeasurementsandModelEstimatesoverVineyard,Cotton,GrassandDeciduousForestinSummer[J].AtmosphericEnvironment,1996,30(13):2360-2369.WESELYML,DOSKEYPV,SHANNONJD.DepositionParameterizationsfortheIndustrialSourceComplex(ISC3)Model[R].Illinois:EnvironmentalResearchDivision,ArgonneNationalLaboratory9700SouthCassAvenueArgonne,2002.ZHANGL,BROOKJR,VETR.ARevisedParameterizationforGaseousDryDepositioninAir-qualityModels[J].AtmosphericChemistryandPhysics,2003,3(6):2067-2082.KORP,KHARRATR.ModelingofAsphalteneParticleDepositionfromTurbulentoilFlowinTubing:ModelValidationandaParametricStudy[J].Petroleum,2016,2⑷：393-398.PILEGAARDK,HUMMELSH0JP,JENSENNO.FluxesofOzoneandNitrogenDioxideMeasuredbyEddtCorrelationoveraHarvestedWheatField[J].AtmosphericEnvironment,1998,32(7):1167-1177.GEROSAG,VITALEM,FINCOA,etal.OzoneUptakebyanEvergreenMediterraneanForest(QuercusIlex)inItaly.PartI:MicrometeorologicalFluxMeasurementsandFluxPartitioning[J].AtmosphericEnvironment,2005,39(18):3255-3266.HORVATHL,WEIDINGERT,NAGY乙EstimationofDryDepositio