• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    植被下墊面Z0的估算及其改進(jìn)影響評(píng)估

    2022-09-20 08:42:36田春艷崔寅平陳修治
    中國(guó)環(huán)境科學(xué) 2022年9期
    關(guān)鍵詞:下墊面粗糙度風(fēng)速

    田春艷,崔寅平,申 沖,陳修治,沈 傲,樊 琦*

    植被下墊面0的估算及其改進(jìn)影響評(píng)估

    田春艷1,2,崔寅平1,申 沖3,陳修治1,沈 傲1,樊 琦1*

    (1.中山大學(xué)大氣科學(xué)學(xué)院,廣東省氣候變化與自然災(zāi)害研究重點(diǎn)實(shí)驗(yàn)室,南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(珠海),廣東 珠海 519082;2.中國(guó)民用航空中南地區(qū)空中交通管理局氣象中心, 廣東 廣州 510000;3.廣州市氣候與農(nóng)業(yè)氣象中心,廣東 廣州 511430)

    利用衛(wèi)星遙感反演數(shù)據(jù)和多種觀測(cè)資料,基于形態(tài)學(xué)-遙感反演的方法,估算我國(guó)主要植被下墊面的空氣動(dòng)力學(xué)粗糙度長(zhǎng)度(0),將其應(yīng)用于區(qū)域大氣化學(xué)模式(WRF-Chem)中改進(jìn)模式默認(rèn)0值,并結(jié)合全國(guó)氣象和污染物觀測(cè)數(shù)據(jù)進(jìn)行改進(jìn)效果評(píng)估,探究0改進(jìn)對(duì)模式氣象場(chǎng)和化學(xué)場(chǎng)模擬結(jié)果的影響.結(jié)果表明:(1)0高值主要位于森林類(lèi)下墊面,可超過(guò)1m;農(nóng)田、草地類(lèi)下墊面0值較小,低于0.5m,其余植被下墊面的0值基本介于兩者之間.(2)0改進(jìn)后,10m風(fēng)速(WS10)和地表溫度(TSK)的改進(jìn)效果較好,而2m溫度(2)和相對(duì)濕度(RH)的改進(jìn)效果不明顯;從空間分布和不同下墊面結(jié)果來(lái)看,0對(duì)TSK主要是增溫作用,對(duì)WS10是減小作用;就改進(jìn)比(各要素改進(jìn)后的模擬值與改進(jìn)前的模擬值之差比上改進(jìn)前的模擬值,下同)而言,0改進(jìn)對(duì)WS10的影響大于TSK.從對(duì)污染物濃度效果評(píng)估看,0對(duì)PM2.5和PM10模擬改進(jìn)效果較好,對(duì)其它污染物(SO2、NO2、O3)的改進(jìn)效果不明顯.(3)對(duì)比氣象要素和污染物濃度效果評(píng)估可知,0改進(jìn)對(duì)氣象要素的影響大于污染物濃度,其主要是通過(guò)影響氣象場(chǎng)來(lái)間接影響污染物濃度.總的來(lái)看,0改進(jìn)影響最大的是氣象場(chǎng)的10m風(fēng)速,考慮其原因在于改進(jìn)后的0普遍高于改進(jìn)前模式默認(rèn)0值,而0是表征地表粗糙度的變量,由于地表粗糙度的增加,增強(qiáng)了對(duì)氣流的阻礙作用,使得近地層風(fēng)速減小.

    空氣動(dòng)力學(xué)粗糙度長(zhǎng)度;植被下墊面;數(shù)值模擬;WRF-Chem模式;PM2.5

    陸面過(guò)程是指陸地表面與大氣之間進(jìn)行物質(zhì)、能量等傳輸和交換的過(guò)程,是影響區(qū)域氣候和環(huán)境變化的重要原因[1],主要通過(guò)湍流運(yùn)動(dòng)實(shí)現(xiàn),而地表粗糙是近地面層湍流運(yùn)動(dòng)的根本原因.空氣動(dòng)力學(xué)粗糙度長(zhǎng)度(0)是陸地表面風(fēng)速降低到零時(shí)對(duì)應(yīng)的高度,它表征地表的空氣動(dòng)力學(xué)特征,是邊界層湍流參數(shù)化方案中的重要參數(shù)[2-3].準(zhǔn)確得出0值及其變化規(guī)律是數(shù)值模式準(zhǔn)確模擬地表通量的前提,同時(shí)亦直接影響數(shù)值模式的模擬性能[4-5].然而,在數(shù)值模式中,0往往被處理為一個(gè)常數(shù),賦予相同的固定值,沒(méi)有體現(xiàn)時(shí)空的差異性,這無(wú)疑會(huì)對(duì)模式的模擬結(jié)果造成一定程度的影響.

    近年來(lái),一些學(xué)者嘗試將不同方法計(jì)算得到的0應(yīng)用到模式中以更真實(shí)地反映0的影響,進(jìn)而提升模式的模擬性能.王凱嘉[6]通過(guò)對(duì)夏季風(fēng)過(guò)渡區(qū)的農(nóng)田、混合型、草地3種下墊面0估算發(fā)現(xiàn),陸面模式中對(duì)應(yīng)的3種下墊面默認(rèn)0值(0.065, 0.08, 0.1m)與估算0值(0.318, 0.218, 0.049m)相差甚遠(yuǎn).同時(shí),其對(duì)3種下墊面0的動(dòng)態(tài)變化研究發(fā)現(xiàn):草地下墊面0的動(dòng)態(tài)變化范圍為0.012~0.25m,混合型下墊面0變化為 0.12~0.35m,農(nóng)田下墊面0在9月后相對(duì)于前期會(huì)迅速增加,這表明與模式將0作為幾何常量不同,植被下墊面0的動(dòng)態(tài)變化十分顯著.胡文超等[7]對(duì)河西走廊地區(qū)0實(shí)測(cè)值與模擬值的差異性分析發(fā)現(xiàn):模式默認(rèn)值不能準(zhǔn)確地反映下墊面0的非均一性,其差異對(duì)于草原下墊面最大可達(dá)375%,隨著植被下墊面復(fù)雜程度增加,差異程度也增加.劉野等[8]利用我國(guó)典型半干旱區(qū)野外試驗(yàn)站的觀測(cè)資料,結(jié)合多種方法估算試驗(yàn)站0,結(jié)果表明0在半干旱區(qū)具有明顯的季節(jié)和年際變化特征,且植被低矮的下墊面0與陸面模式中默認(rèn)值相差較大,其將修正后0代入模式,發(fā)現(xiàn)明顯改善了模式對(duì)該地區(qū)地表感熱通量的模擬.鮑艷等[9]通過(guò)對(duì)陸面模式BATS(Biosphere-Atmosphere Transfer Scheme)和LSM(Land Surface Model)設(shè)置敏感性試驗(yàn)檢驗(yàn)陸面模式參數(shù),結(jié)果表明模式默認(rèn)值對(duì)地表輻射、熱通量模擬存在偏差,地表溫度模擬偏低,通過(guò)對(duì)0參數(shù)重新取值,并考慮不同文獻(xiàn)所用方案改進(jìn)模式后,模式對(duì)地表熱過(guò)程模擬效果明顯提高,模擬的地表熱通量、地表溫度與觀測(cè)值的相關(guān)可達(dá)0.92以上.

    吳曉波[10]基于WRF(The Weather Research and Forecasting Model)和CFD(Computational Fluid Dynamic)耦合的復(fù)雜地形風(fēng)場(chǎng)數(shù)值研究表明,0對(duì)實(shí)際地形的CFD數(shù)值模擬結(jié)果影響比較大,使用精細(xì)化0模擬的風(fēng)場(chǎng)誤差較小.Jee等[11]利用高分辨率數(shù)值模式(WISE-WRF)對(duì)首爾地區(qū)城市下墊面0改善的靈敏度分析發(fā)現(xiàn),城市下墊面0對(duì)摩擦速度和風(fēng)速有一定的影響,與默認(rèn)結(jié)果相比,改進(jìn)后的模型能更準(zhǔn)確地模擬溫度和風(fēng)速.孫行知等[12]利用WRF模式設(shè)計(jì)新的有效0參數(shù)方案對(duì)2003年?yáng)|亞夏季風(fēng)過(guò)程進(jìn)行區(qū)域氣候模擬,結(jié)果表明,新方案改進(jìn)了對(duì)東亞夏季風(fēng)降水分布、強(qiáng)度和雨帶推進(jìn)等的模擬結(jié)果.Varqueze等[13]利用幾何參數(shù)估算日本范圍內(nèi)的0分布,應(yīng)用到模式中替換默認(rèn)的0,更新后的WRF提高了東京、名古屋和大阪地區(qū)溫度和風(fēng)速的模擬結(jié)果,對(duì)東京地區(qū)的研究發(fā)現(xiàn),0的影響是比較主要的,尤其是對(duì)風(fēng)場(chǎng).Shen等[14]利用更新0后的WRF-Chem模式對(duì)珠江三角洲地區(qū)的熱、動(dòng)力環(huán)境和CO濃度的模擬研究發(fā)現(xiàn),更新后的模式提高了對(duì)風(fēng)速的模擬效果,同時(shí)對(duì)地表溫度、熱通量、摩擦速度和邊界層高度等變量都有明顯影響;另外,通過(guò)影響風(fēng)速,更新后的模式使得CO夜間模擬濃度增加.

    已有模式中0的改進(jìn)工作多針對(duì)某一地區(qū)或某一類(lèi)型的下墊面進(jìn)行,在區(qū)域范圍內(nèi)多種不同下墊面方面的研究工作還需進(jìn)一步加強(qiáng)完善.我國(guó)幅員遼闊,下墊面類(lèi)型豐富,模式中默認(rèn)的0與實(shí)際存在較大的偏差,無(wú)法體現(xiàn)0真實(shí)的時(shí)空分布特征.因此,本研究利用衛(wèi)星遙感資料估算區(qū)域范圍內(nèi)不同植被下墊面類(lèi)型的0值并將其應(yīng)用于數(shù)值模式,評(píng)估0改進(jìn)對(duì)模式模擬結(jié)果的影響.

    1 研究資料與方法

    1.1 研究資料

    模式中使用的土地覆蓋類(lèi)型資料為2015年MODIS (Moderate Resolution Imaging Spectroradiometer)第6版數(shù)據(jù)產(chǎn)品(MCD12Q1),空間分辨率為500m,該土地覆蓋分類(lèi)采用國(guó)際地圈生物圈計(jì)劃(IGBP)的土地分類(lèi)方案,共分為17類(lèi)(表1).本文中估算空氣動(dòng)力學(xué)粗糙度長(zhǎng)度時(shí),所研究的范圍為中國(guó)大陸,包括海南島和臺(tái)灣島(圖1a),該區(qū)域內(nèi)土地覆蓋類(lèi)型豐富,包含了IGBP所含的17種分類(lèi),但由于估算0涉及的葉面積指數(shù)數(shù)據(jù)針對(duì)的是植被下墊面,因此,本文后續(xù)改進(jìn)的0主要為表1中所列1~10類(lèi)自然植被下墊面和12、14類(lèi)受人為活動(dòng)影響的農(nóng)田類(lèi)下墊面.由各種類(lèi)型下墊面占比(表1)可知,我國(guó)植被下墊面占比最高是草地,其次是農(nóng)田,草地和農(nóng)田共占我國(guó)近一半的陸面面積,占比最低為郁閉灌木和稀疏灌木,不到0.1%.從各植被下墊面在我國(guó)的地理位置來(lái)看,占比較高的草地主要位于我國(guó)內(nèi)蒙古中、東部地區(qū)到青藏高原東、南部地區(qū);農(nóng)田主要位于我國(guó)山東、河南省及其周邊部分地區(qū)和黑龍江、吉林和遼寧的西部地區(qū);其它下墊面主要位于我國(guó)南方地區(qū)和黑龍江、吉林和遼寧的東部地區(qū)以及內(nèi)蒙古黑龍江交界處最北端附近.

    葉面積指數(shù)(LAI)數(shù)據(jù)采用MODIS的第6版數(shù)據(jù)產(chǎn)品(MCD15A2H),是美國(guó)國(guó)家航空航天局(NASA)2015年最新發(fā)布的LAI產(chǎn)品,為衛(wèi)星Terra和Aqua上的MODIS儀器組合產(chǎn)品[15].該產(chǎn)品空間分辨率為500m,時(shí)間分辨率為8d.LAI定義為單位地表面積上的植被單側(cè)葉片面積,即單側(cè)葉片面積占土地面積的比[16-17].它與植被的密度、結(jié)構(gòu)等生物學(xué)特性和環(huán)境條件有關(guān),是反映植被較好的動(dòng)態(tài)指標(biāo).

    樹(shù)木類(lèi)植被,由于其一年四季只存在葉落葉生的過(guò)程,植被高度基本不變,而草地和農(nóng)田類(lèi)植被,植被高度隨植被生長(zhǎng)變化而變化[18].因此,本研究對(duì)研究區(qū)域內(nèi)所包含的樹(shù)木類(lèi)植被高度選取衛(wèi)星反演的森林樹(shù)冠高度數(shù)據(jù)(圖1b)進(jìn)行計(jì)算.森林樹(shù)冠高度數(shù)據(jù)[19]是從2005年的GLAS(Geoscience Laser Altimeter System)激光雷達(dá)信號(hào)獲取,其空間分辨率為1km,本研究將其插值到模式對(duì)應(yīng)的27km網(wǎng)格中使用.而農(nóng)田草地下墊面則根據(jù)文獻(xiàn)數(shù)據(jù)整理,取草地冬、夏季的植被冠層高度分別為0.06, 0.2m,農(nóng)田下墊面夏季為1.05m,冬季由于沒(méi)有統(tǒng)計(jì)高度結(jié)果,采取森林樹(shù)冠高度數(shù)據(jù)中網(wǎng)格內(nèi)的默認(rèn)值.其余未作更改的植被下墊面網(wǎng)格,認(rèn)為其下墊面為樹(shù)木類(lèi)植被,其植被高度不隨四季變化.

    全國(guó)氣象觀測(cè)數(shù)據(jù)來(lái)自國(guó)家氣象科學(xué)數(shù)據(jù)共享服務(wù)平臺(tái)上的674個(gè)觀測(cè)站點(diǎn)的氣象數(shù)據(jù),選取了2m溫度(2)、地表溫度(TSK)和10m風(fēng)速(WS10)以及相對(duì)濕度(RH)4個(gè)變量進(jìn)行研究,該數(shù)據(jù)為2015年1月和7月的日均值數(shù)據(jù).全國(guó)空氣質(zhì)量觀測(cè)數(shù)據(jù)來(lái)自中國(guó)環(huán)境監(jiān)測(cè)總站1479個(gè)觀測(cè)站點(diǎn)的污染物濃度逐時(shí)數(shù)據(jù).數(shù)據(jù)使用時(shí)間段為2015年1月和7月,使用的污染物包括SO2、NO2、O3、PM2.5、PM10.

    1.常綠針葉林;2.常綠闊葉林;3.落葉針葉林;4.落葉闊葉林;5.混交林;6.郁閉灌木;7.稀疏灌木;8.多樹(shù)草原;9.稀樹(shù)草原;10.草地;11.永久濕地;12.農(nóng)田;13.城市和建成區(qū);14.農(nóng)田和自然植被的鑲嵌體;15.雪、冰;16.裸地或植被覆蓋地;17.水體;審圖號(hào):GS(2021)5875號(hào)|GS(2021)8810號(hào),下同

    1.2 Z0估算方法

    本研究使用的0估算方法參考Raupach[20-21]通過(guò)風(fēng)洞試驗(yàn)研究提出的一種估算0的半經(jīng)驗(yàn)?zāi)M,并結(jié)合Jasinski等[22]的研究,0具體計(jì)算公式如下:

    表1 Z0估算涉及的IGBP土地覆蓋類(lèi)型參數(shù)

    注:R為粗糙元阻力系數(shù);為粗糙元遮陰系數(shù);h為風(fēng)速;為摩擦速度;為植被冠層面積指數(shù);是綠色葉面積指數(shù)估算系數(shù);最大歸一化植被指數(shù)與該植被覆蓋率為100%時(shí)的最大歸一化植被指數(shù)值的比值;L,min為不同植被類(lèi)型的月最小莖部和死葉面積指數(shù);1-為枯葉的脫落比率.

    1.3 模式介紹及其設(shè)置

    WRF-Chem(Weather Research and Forecasting- Chemistry)模式是美國(guó)國(guó)家大氣研究中心、美國(guó)大氣海洋局及美國(guó)西部太平洋國(guó)家實(shí)驗(yàn)室等單位共同開(kāi)發(fā)的新一代大氣化學(xué)傳輸模式[24].

    本研究采用WRF-Chem 3.6.1版本,模擬研究區(qū)域?yàn)?6.2758°~49.4169°N,78.9415°~141.0585°E,采用LAMERT投影方式,水平格距為27km,格點(diǎn)數(shù)為180×130,垂直方向設(shè)置35層,層頂?shù)臍鈮簽?0hPa.模式采用的物理化學(xué)參數(shù)化方案包括: Morrison (2moments)微物理方案,RRTM長(zhǎng)波輻射方案,Goddard短波輻射方案,Monin-Obukhov (Janjic)近地層方案,Grell 3D ensemble積云對(duì)流方案, unified Noah land-surface model陸面過(guò)程, Mellor-Yamada-Janjic TKE邊界層方案,UCM城市冠層方案; MOSAIC氣溶膠機(jī)制,CBMZ大氣化學(xué)機(jī)制, Fast-J光解機(jī)制.使用1°×1°的NCEP 全球再分析資料作為模式模擬的初始條件和邊界條件.模擬時(shí)間段為世界時(shí)2014年12月20日00:00~2015年02月01日00:00,2015年06月20日00:00時(shí)~ 2015年08月01日00:00時(shí).模式中采用的排放源為清華大學(xué)制作的2015年MEIC(中國(guó)大陸)和2010的MIX(其它區(qū)域)清單[25-26].土地覆蓋類(lèi)型為2015年MODIS第6版數(shù)據(jù)產(chǎn)品(MCD12Q1)插值到對(duì)應(yīng)模式網(wǎng)格數(shù)據(jù).本次模擬分為Base、Case試驗(yàn),Base試驗(yàn)采用模式默認(rèn)0值,Case試驗(yàn)采用本文估算的0值.

    2 結(jié)果與分析

    2.1 Z0估算及特征分析

    圖3 不同文獻(xiàn)Z0調(diào)研值

    由圖2可知,估算0高值區(qū)主要位于我國(guó)的西南、東南沿海以及東北的東南地區(qū),其中,西南地區(qū)的0值最高,在1月可達(dá)2m;其次是東南沿海地區(qū),約1.5m;而0值最低的地方,主要位于華北平原和內(nèi)蒙古西北部至青藏高原地區(qū).本文估算0量級(jí)與文獻(xiàn)[2,27]基本一致.進(jìn)一步結(jié)合圖1a,從下墊面類(lèi)型來(lái)看,高值區(qū)對(duì)應(yīng)的位置主要在森林樹(shù)木類(lèi)下墊面,而低值區(qū)主要位于農(nóng)田草地類(lèi)下墊面.對(duì)比文獻(xiàn)調(diào)研得到的0值[4,6-8,16,18,28-70],對(duì)于同一類(lèi)植被下墊面,受研究方法和地形區(qū)域等影響,不同文獻(xiàn)調(diào)研結(jié)果具體值存在一定差異,但是就量級(jí)而言,同一類(lèi)下墊面類(lèi)型的量級(jí)基本保持一致,本文的研究結(jié)果與文獻(xiàn)調(diào)研結(jié)果結(jié)論基本一致(圖3).由圖2可知,0呈現(xiàn)動(dòng)態(tài)變化特征,1、7月的0結(jié)果存在明顯差異.將7月0值減去1月0值得0空間差值.對(duì)比不同下墊面類(lèi)型,其中0減小的地區(qū)主要在森林類(lèi)茂盛的植被下墊面.研究表明[27,30-31],對(duì)于一些植被,由于其冠層在生長(zhǎng)季葉片密度大,郁閉度高,相鄰的葉片互相遮蔽使得表面動(dòng)力學(xué)光滑,從而導(dǎo)致其零平面位移高度抬升,而植被0值降低.而對(duì)于農(nóng)田類(lèi)密度不大的植被下墊面,則主要由植被高度決定,隨著生長(zhǎng)季的變化,植被越高,氣流流經(jīng)下墊面時(shí),受植被高度的影響就越大,垂直方向的風(fēng)速分布隨高度發(fā)生變化,0增加[71].

    2.2 模式中Z0改進(jìn)及模擬結(jié)果評(píng)估

    用本文估算0值減去模式默認(rèn)0值得到1、7月0差值空間分布(圖4).除了華北平原和內(nèi)蒙古等部分地區(qū)外,模式默認(rèn)0值明顯低于估算0值,最大差值主要出現(xiàn)在我國(guó)西南和東南沿海地區(qū),可達(dá)1m以上,對(duì)比不同的季節(jié),冬季的差值大于夏季.由于本文主要是針對(duì)植被下墊面進(jìn)行0估算,所以對(duì)于第11類(lèi)永久濕地、15類(lèi)雪和冰、16類(lèi)裸地或稀疏植被、17類(lèi)水體下墊面采取模式默認(rèn)值,未作更改,而對(duì)于第3類(lèi)落葉林和第6類(lèi)郁閉灌木,研究區(qū)域內(nèi)不存在此類(lèi)土地覆蓋類(lèi)型,因此也未對(duì)模式進(jìn)行修改.由表2可知,對(duì)于模式默認(rèn)0值,其所對(duì)應(yīng)的第1~5類(lèi)森林類(lèi)下墊面在模式中采用了同樣的值,均為0.5m,而從估算0來(lái)看,不同的森林下墊面由于植被本身特性不同,其所對(duì)應(yīng)的0不同,森林類(lèi)下墊面的0值都在1m以上,結(jié)合圖3可見(jiàn),本文對(duì)此類(lèi)型下墊面0的改進(jìn)效果最為顯著.而對(duì)于其它下墊面類(lèi)型來(lái)說(shuō),估算的0值也明顯高于模式默認(rèn)值,其中,差值最大的是第5類(lèi)混交林,在冬季,相差可達(dá)1.17m,其次是第8類(lèi),最大可差0.97m.

    通過(guò)用估算值和模式默認(rèn)值對(duì)比分析可知,模式默認(rèn)值不能很好地體現(xiàn)0時(shí)空變化特征,模式默認(rèn)0值存在明顯低估,不能很好地反映出不同下墊面類(lèi)型的特征,這對(duì)于陸氣相互作用的模擬存在一定的影響.

    表2 IGBP土地覆蓋類(lèi)型對(duì)應(yīng)的模式默認(rèn)(Base)及估算(Case)Z0值

    注:SFZ0為不同下墊面模式默認(rèn)和本論文估算的0平均值,均作為初值代入模式計(jì)算.

    2.3 模式模擬結(jié)果評(píng)估

    評(píng)估模擬效果的統(tǒng)計(jì)參數(shù)[72-73]主要有:平均偏差(BIAS)、平均絕對(duì)誤差(MAE)、均方根誤差(RMSE)、相關(guān)系數(shù)()、符合指數(shù)(IOA).

    2.3.1 氣象場(chǎng)模擬結(jié)果評(píng)估 通過(guò)WRF-Chem模式模擬得到0改進(jìn)前后氣象要素的模擬值,利用1.1節(jié)中介紹的中國(guó)地面氣象站的氣象觀測(cè)數(shù)據(jù)(2、TSK、WS10、RH)進(jìn)行氣象場(chǎng)結(jié)果評(píng)估.如圖5所示,對(duì)于2和TSK,從1月和7月模擬與觀測(cè)的對(duì)比來(lái)看,1月模式模擬的溫度變化趨勢(shì)和范圍與觀測(cè)基本接,變化范圍在-5~5℃,TSK的模擬結(jié)果略低于觀測(cè);7月模擬和觀測(cè)的TSK在20~35℃之間變化,2的變化范圍在20~30℃之間,略低于TSK,TSK的模擬較觀測(cè)低,2則相反.結(jié)合0改進(jìn)前后的溫度箱型圖(圖略)對(duì)比來(lái)看,1月2Base方案中位數(shù)、均值相對(duì)于Case而言更接近觀測(cè),但從上下四分位數(shù)來(lái)看,Case方案更接近觀測(cè);對(duì)于TSK而言,Case和Base方案的結(jié)果較為接近,但就中位數(shù)、均值、上下四分位數(shù)及最大最小值來(lái)看,Case方案的結(jié)果整體而言更好;7月T2和TSK的中位數(shù)和均值較為接近.整體來(lái)看,2Base的結(jié)果優(yōu)于Case方案,而TSK則是Case方案的結(jié)果優(yōu)于Base.對(duì)于WS10,1月和7月模式模擬結(jié)果都存在高估的現(xiàn)象,從Base、Case與觀測(cè)的對(duì)比來(lái)看,Base和Case兩者的變化趨勢(shì)與觀測(cè)基本一致,但是就變化范圍來(lái)看,可以明顯地發(fā)現(xiàn),Case的模擬結(jié)果低于Base,更接近觀測(cè).對(duì)于RH,和WS10相反,1月和7月模式模擬結(jié)果相對(duì)于觀測(cè)有所偏低,Base的模擬結(jié)果相對(duì)Case而言,更接近觀測(cè).總的來(lái)看,0改進(jìn)效果最明顯的是WS10,其次是TSK.考慮其原因在于0是表征風(fēng)速降低到零時(shí)對(duì)應(yīng)的高度,在本文中0改進(jìn)后的值相對(duì)模式默認(rèn)值有了較大的提高,其高度的增加使得同等高度對(duì)應(yīng)的風(fēng)速明顯降低,從而使得改進(jìn)后模式模擬風(fēng)速明顯較低,更接近觀測(cè)值,這進(jìn)一步反映了模式模擬0值被明顯低估.而TSK作為地表變量,更容易受到地表參數(shù)0的影響,0值的變化,使得地表粗糙度發(fā)生變化,進(jìn)一步導(dǎo)致地表湍流運(yùn)動(dòng)的變化,使得熱量的輸送交換產(chǎn)生較大的影響,從而影響TSK的模擬結(jié)果.

    由以上分析可知,0改進(jìn)對(duì)TSK和WS10兩個(gè)變量的模擬效果有較大的提升,為了解兩者在不同下墊面類(lèi)型上的提升差異,特定義“改進(jìn)比”來(lái)進(jìn)一步分析.“改進(jìn)比”定義:100′(Case-Base)/Obs和100′(Case-Base)/Base,以表示0改進(jìn)相對(duì)于觀測(cè)值和Base方案的改進(jìn)比例,其中,由于1月TSK可能為負(fù)值,為了排除其負(fù)號(hào)對(duì)結(jié)果的影響,對(duì)其分母取絕對(duì)值,即100′(Case-Base)/|Obs| ;100′(Case-Base)/ |Base|.

    由圖6可知,1月0改進(jìn)對(duì)TSK主要表現(xiàn)為增溫作用,改進(jìn)比較大的地區(qū)主要位于京津冀沿山西、陜西到四川云南一帶,最大超過(guò)50%,對(duì)比TSK的Base空間分布圖(圖略),其溫度零值帶是影響其改進(jìn)比超過(guò)50%的主要原因,其次較大的改進(jìn)比主要在長(zhǎng)江以南地區(qū),約20%左右,其余地區(qū)的改進(jìn)較小.7月由于溫度相對(duì)1月升高較大,其改進(jìn)比相對(duì)1月而言所有降低,改進(jìn)較大的地區(qū)主要位于四川盆地西北方,約15%左右,其余地區(qū)改進(jìn)效果較小,在±5%之間變化.結(jié)合圖4對(duì)比分析可知0增加的區(qū)域與1月溫度增加的區(qū)域基本一致,排除溫度零值帶影響后,兩者的改進(jìn)程度也基本一致,7月相關(guān)性則不如1月明顯.就WS10來(lái)看,除了山東及其附近省市、內(nèi)蒙古北部等部分地區(qū)風(fēng)速有所增加外,0改進(jìn)主要使風(fēng)速減小,風(fēng)速改進(jìn)最大值可達(dá)50%以上,且1月影響大于7月,考慮主要原因在于1月大氣層結(jié)穩(wěn)定,天氣變化等其他因素影響較小.結(jié)合圖4,可發(fā)現(xiàn)0改進(jìn)與WS10呈現(xiàn)較強(qiáng)的負(fù)相關(guān),考慮其主要由于0值增加,地表粗糙度增大,阻礙了氣流流動(dòng),使得WS10減小.

    考慮到灌木等下墊面類(lèi)型占比過(guò)低,進(jìn)一步針對(duì)占比大于2%的下墊面類(lèi)型了解0改進(jìn)比的影響.由圖7可知,對(duì)于1月和7月的TSK,0改進(jìn)對(duì)其影響主要是呈現(xiàn)增溫效果,其對(duì)不同下墊面類(lèi)型的改進(jìn)程度不一,在1月,相對(duì)于觀測(cè),4類(lèi)落葉闊葉林和12類(lèi)農(nóng)田下墊面的改進(jìn)比可達(dá)16%左右,其它下墊面類(lèi)型的改進(jìn)比則主要在4%左右;相對(duì)于Base方案,改進(jìn)比較高的是5類(lèi)混交林和12類(lèi)農(nóng)田下墊面,其值可達(dá)12%左右,其它下墊面類(lèi)型的改進(jìn)效果和觀測(cè)較為接近,在4%左右.7月,不同的下墊面相對(duì)于Base方案的改進(jìn)比都略高于觀測(cè),但數(shù)值較為接近,且相對(duì)于1月而言,改進(jìn)比較低,在3%以下.對(duì)于WS10,1月和7月的改進(jìn)結(jié)果較為一致,對(duì)于不同的下墊面類(lèi)型,其主要起到降低風(fēng)速的作用,相對(duì)于Base方案,觀測(cè)的改進(jìn)比大于Base,Base的改進(jìn)比在20%左右,而相對(duì)于觀測(cè)的改進(jìn)比除了農(nóng)田下墊面和7月的草地,其改進(jìn)比基本在20%以上,其中,改進(jìn)比最大的8類(lèi)多樹(shù)草原下墊面,在1月可達(dá)100%左右.綜上分析可知,0對(duì)溫度有增溫作用,對(duì)風(fēng)速有降速作用,考慮其原因在于粗糙度長(zhǎng)度增加,導(dǎo)致同等高度風(fēng)速降低,湍流運(yùn)動(dòng)減弱,對(duì)熱量的傳送交換減小.就改進(jìn)比而言,0改進(jìn)對(duì)風(fēng)速的影響大于溫度的影響,考慮其原因在于0和風(fēng)速是直接的動(dòng)力影響,而對(duì)TSK的影響則是通過(guò)熱力的間接影響.0改進(jìn)對(duì)1月的TSK影響大于7月,這可能與夏季的溫度高于冬季有關(guān),而冬夏季風(fēng)速大小差異不大有關(guān).對(duì)比可知,盡管是對(duì)Base方案做的改進(jìn),但是其模擬結(jié)果相對(duì)于觀測(cè)而言,有較大的改善,其中風(fēng)速相對(duì)于觀測(cè)而言,改進(jìn)較大.

    圖7 Case方案模式模擬的TSK和WS10相對(duì)于觀測(cè)和Base方案的改進(jìn)比

    Fig.7 Improvement ratios of the simulated TSK and WS10 for the Case scheme compared with the observations and Base scheme

    L2.常綠闊葉林; L4.落葉闊葉林;L5.混交林; L8.多樹(shù)草原;L9.稀樹(shù)草原;L10.草地;L12.農(nóng)田

    2.3.2 污染物模擬結(jié)果評(píng)估 通過(guò)WRF-Chem模式模擬得到0改進(jìn)前后污染物濃度的模擬結(jié)果,結(jié)合1.1節(jié)中介紹的中國(guó)環(huán)境監(jiān)測(cè)站的常規(guī)污染物(SO2、NO2、O3、PM2.5、PM10)數(shù)據(jù)進(jìn)行污染物濃度模擬結(jié)果評(píng)估.如圖8所示,對(duì)于PM10和PM2.5,1月和7月模式模擬的結(jié)果和觀測(cè)值的變化趨勢(shì)基本一致,但是模式模擬結(jié)果較觀測(cè)偏低.結(jié)合箱型圖(圖略)來(lái)看,1月PM2.5和PM10的濃度變化幅度大于7月,而對(duì)比改進(jìn)前后,雖然Case和Base的結(jié)果基本一致,但是仍可看出,就均值和中位數(shù)而言,Case方案相對(duì)于Base方案而言更接近觀測(cè).對(duì)于SO2和NO2,1月的模式模擬結(jié)果較觀測(cè)偏低,7月結(jié)果基本一致,就其中位數(shù)和均值而言,Case方案更接近觀測(cè).對(duì)于O3,1月和7月模式的模擬結(jié)果和觀測(cè)基本一致,1月份的中位數(shù)和均值Case方案相對(duì)Base更接近觀測(cè),而7月份均值Case優(yōu)于Base,但是中位數(shù)Base更接近觀測(cè).整體來(lái)看,除了7月O3濃度高于1月外,其余污染物濃度皆是1月高于7月.就0改進(jìn)前后的對(duì)比來(lái)看,Case和Base的結(jié)果基本一致,整體而言, Case方案的均值和中位數(shù)更接近觀測(cè).

    實(shí)線:1月;虛線:7月

    由圖9可知,0改進(jìn)對(duì)PM2.5和NO2空間分布在不同區(qū)域的影響有所不同,除了1月NO2濃度改進(jìn)比和0改進(jìn)量(圖4)空間分布之間有較強(qiáng)的相關(guān)性,其余的相關(guān)性均不明顯.從兩者改進(jìn)前濃度分布圖(圖略)來(lái)看,不同區(qū)域污染物濃度分布可能是影響其不同下墊面類(lèi)型相關(guān)性不同的主要原因,尤其是對(duì)于森林類(lèi)下墊面,由于其污染物濃度較低,容易受其它因素影響大,從而與0的相關(guān)性不如其它下墊面.

    進(jìn)一步利用“改進(jìn)比”分析不同下墊面上0對(duì)PM2.5和NO2的改進(jìn)程度,由圖10可知,對(duì)于PM2.5,1月0改進(jìn)主要是使得污染物濃度增加考慮其原因在于風(fēng)速降低,不利于污染物的輸送,Base的改進(jìn)比大于觀測(cè),約1.5%,觀測(cè)的改進(jìn)比約為1%左右,改進(jìn)比最低的8類(lèi)下墊面,低于0.5%;7月,除2類(lèi)常綠闊葉林和9類(lèi)稀樹(shù)草原污染物濃度改進(jìn)比增加外,其余下墊面均降低,和1月一樣,對(duì)相對(duì)于Base的改進(jìn)比大于觀測(cè),同時(shí),可發(fā)現(xiàn)7月的改進(jìn)比相對(duì)于1月的改進(jìn)比大,考慮主要原因在于1月污染物濃度較大,使得分母大,整體改進(jìn)比減小.就NO2來(lái)看,Base改進(jìn)比大于觀測(cè),1月稀樹(shù)草原和7月8類(lèi)多樹(shù)樹(shù)草原下墊面改進(jìn)效果不明顯,約0%.就增減情況來(lái)看,除1月8類(lèi)多樹(shù)草原下墊面和7月10類(lèi)草地下墊面外,其余下墊面在1月和7月的改進(jìn)比均是增加.7月NO2的改進(jìn)比大于1月.總得來(lái)看,0改進(jìn)對(duì)污染物的模擬結(jié)果有一點(diǎn)的改進(jìn),但是對(duì)于氣象場(chǎng)而言,其改進(jìn)效果不明顯,考慮其原因在于,0的改進(jìn)主要是作用于氣象場(chǎng),通過(guò)氣象場(chǎng)來(lái)間接影響污染物濃度,而0改進(jìn)的氣象場(chǎng)的影響存在一定的差異,同時(shí)影響污染物濃度存在較多物理化學(xué)環(huán)境等因素,從而使得污染物濃度的改進(jìn)效果不明顯.

    L2.常綠闊葉林; L4.落葉闊葉林;L5.混交林; L8.多樹(shù)草原;L9.稀樹(shù)草原;L10.草地;L12.農(nóng)田

    本研究采用觀測(cè)資料進(jìn)行模式模擬結(jié)果進(jìn)行驗(yàn)證對(duì)比,但是觀測(cè)資料在植被下墊面類(lèi)型上的數(shù)據(jù)量相對(duì)較少,因此模式在這類(lèi)型下墊面上結(jié)果的驗(yàn)證評(píng)估還可進(jìn)一步加強(qiáng).另外,本論文植被高度數(shù)據(jù)采用的是衛(wèi)星反演數(shù)據(jù)及文獻(xiàn)調(diào)研數(shù)據(jù),缺乏區(qū)域上更精細(xì)化的不同下墊面高度數(shù)據(jù),因此只針對(duì)冬季和夏季的典型月份1月和7月開(kāi)展了0改進(jìn)對(duì)模式模擬結(jié)果的影響研究,未能對(duì)多個(gè)月份開(kāi)展更多時(shí)次的模擬工作.同時(shí),只針對(duì)植被下墊面進(jìn)行了0估算和分析研究,沒(méi)有對(duì)城市等其它下墊面類(lèi)型的0進(jìn)行估算及模式模擬研究,這也是未來(lái)可進(jìn)一步拓展的部分.

    3 結(jié)論

    3.10高值主要位于森林類(lèi)下墊面,數(shù)值可超過(guò)1m.估算0值與模式默認(rèn)值之間差異最大的是5類(lèi)混交林,農(nóng)田草地類(lèi)的改進(jìn)值較小,低于0.5m.

    3.2 從對(duì)氣象要素站點(diǎn)平均模擬結(jié)果評(píng)估看:2和TSK改進(jìn)前后模擬與觀測(cè)基本一致.TSK的改進(jìn)效果優(yōu)于2.WS10模式模擬結(jié)果改進(jìn)前后都存在高估的現(xiàn)象,但0改進(jìn)后更接近觀測(cè),0改進(jìn)很大程度上提升風(fēng)速的模擬效果.0對(duì)RH的改進(jìn)效果則不明顯.從空間分布和不同下墊面類(lèi)型的平均看,0對(duì)TSK有增加作用,對(duì)WS10有減小作用.就改進(jìn)比而言,0改進(jìn)對(duì)WS10的影響大于TSK;對(duì)于WS10,其相對(duì)于觀測(cè)的改進(jìn)比大于Base,而對(duì)于TSK,7月各下墊面類(lèi)型主要是相對(duì)于Base的改進(jìn)比大于相對(duì)于觀測(cè),1月則針對(duì)不同下墊面類(lèi)型結(jié)果不同.0改進(jìn)對(duì)1月的TSK影響大于7月.

    3.3 從對(duì)污染物濃度站點(diǎn)平均模擬結(jié)果評(píng)估看,1月除O3模擬結(jié)果和觀測(cè)接近外,其余污染物模擬結(jié)果較觀測(cè)偏低,7月除PM2.5和PM10模擬較觀測(cè)偏低外,其余污染物濃度模擬和觀測(cè)基本接近;就改進(jìn)前后的對(duì)比來(lái)看,不同污染物Case和Base的結(jié)果基本一致,但Case的均值和中位數(shù)更接近觀測(cè);對(duì)于PM2.5和PM10,Case的結(jié)果要優(yōu)于Base,而對(duì)于SO2、NO2和O3,其1月改進(jìn)效果優(yōu)于7月.0改進(jìn)對(duì)污染物的影響小于污染物濃度.

    [1] 曾 劍.中國(guó)北方地區(qū)陸面過(guò)程特征和參數(shù)化及其與氣候關(guān)系[D]. 北京:中國(guó)氣象科學(xué)研究院, 2011.

    Zeng J. The characteristics and parameterization of land surface processes and its relationship with climate over northern China [D]. Beijing: Chinese Academy of Meteorological Sciences, 2011.

    [2] 劉小平,董治寶.空氣動(dòng)力學(xué)粗糙度的物理與實(shí)踐意義[J]. 中國(guó)沙漠, 2003,23(4):3-12.

    Liu X P, Dong Z B. Review of aerodynamic roughness length [J]. Journal of Desert Research, 2003,23(4):3-12.

    [3] 張雅靜.陰山北麓地表空氣動(dòng)力學(xué)粗糙度性質(zhì)的試驗(yàn)研究[D]. 呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué), 2007.

    Zhang Y J. Experiment Study on topsoil aerodynamics roughness characteristics in the north of mount Yin [D]. Huhhot: Inner Mongolia Agricultural University, 2007.

    [4] 楊耀先,李茂善,胡澤勇,等.藏北高原高寒草甸地表粗糙度對(duì)地氣通量的影響[J]. 高原氣象, 2014,33(3):626-636.

    Yang Y X, L i M S, Hu Z Y, et al. Influence of surface roughness on Surface-Air Fluxes in alpine meadow over the northern Qinghai-Xizang Plateau [J]. Plateau Meteorology, 2014,33(3):626- 636.

    [5] 張 強(qiáng),姚 彤,岳 平.一個(gè)平坦低矮植被陸面動(dòng)力學(xué)粗糙度多因子參數(shù)化方案及其檢驗(yàn)[J]. 中國(guó)科學(xué):地球科學(xué), 2015,45(11): 1713-1727.

    Zhang Q, Yao T, Yue P.Development and test of a multifactorial parameterization scheme of land surface aerodynamic roughness length for flat land surfaces with short vegetation [J]. Science China: Earth Sciences, 2015,45(11):1713-1727.

    [6] 王凱嘉.夏季風(fēng)影響過(guò)渡區(qū)動(dòng)力學(xué)粗糙度變化特征和影響機(jī)制研究[D]. 蘭州:蘭大學(xué), 2018.

    Wang K J. The change features and influence mechanism research on dynamic roughness of summer monsoon transition zone [D]. Lanzhou: Lanzhou University, 2018.

    [7] 胡文超,張文煜,張 宇,等.河西走廊下墊面粗糙度實(shí)測(cè)值與模擬值的差異性分析[J]. 高原氣象, 2010,29(1):51-55.

    Hu W C, Zhang W Y, Zhang Y, et al. Variance analysis on the simulated and observed values of underlaying surface roughness in Gansu Corridor [J]. Plateau Meteorology, 2010,29(1):51-55.

    [8] 劉 野,郭維棟,宋耀明.半干旱區(qū)關(guān)鍵地表參數(shù)的估算及其對(duì)地氣通量模擬的改進(jìn)[J]. 中國(guó)科學(xué):地球科學(xué), 2015,45(10):1524-1536.

    Liu Y, Guo W D, Song Y M. Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes Simulation [J]. Science China: Earth Sciences, 2015,45(10):1524- 1536.

    [9] 鮑 艷,呂世華.戈壁陸面參數(shù)改進(jìn)對(duì)地表熱狀況影響的數(shù)值模擬[J]. 地球科學(xué)進(jìn)展, 2007,22(特刊):185-193.

    Bao Y, Lü S H. Simulations of Impacts of Land Surface Parameters Improvements on Surface Thermal Condition over Gobi [J]. Advances in Earth Science, 2007,22(Suppl):185-193.

    [10] 吳曉波.基于WRF和CFD耦合的復(fù)雜地形風(fēng)場(chǎng)數(shù)值研究[D]. 武漢:華中科技大學(xué), 2018.

    Wu X B. Numerical Simulation of Wind Field on Complex Terrain Based on Coupling WRF And CFD [D]. Wuhan: Huazhong University of Science & Technology, 2018.

    [11] Jee J B, Min J, Yi C, et al. Sensitivity analysis of the high-resolution WISE-WRF model with the use of surface roughness length in Seoul metropolitan areas [J]. Atmosphere, 2016, 26(1):111-126.

    [12] 孫行知,鐘中,孫源.基于有效空氣動(dòng)力學(xué)參數(shù)新方案對(duì)2003年?yáng)|亞夏季降水模擬的改進(jìn)[C]. 中國(guó)氣象學(xué)會(huì).創(chuàng)新驅(qū)動(dòng)發(fā)展提高氣象災(zāi)害防御能力——S18第四屆研究生年會(huì).中國(guó)氣象學(xué)會(huì):中國(guó)氣象學(xué)會(huì), 2013:408-412.

    Sun X Z, Zhong Z, Sun Y. Parameters Improvement in new effective aerodynamic scheme and its impact on East Asian Summer precipitation simulation in 2003 [C]. Chinese Meteorological Society. Innovation-driven development Improving the Ability of Meteorological Disaster Defense- S18 The Fourth Annual Conference of Graduate Students. Chinese Meteorological Society: Chinese Meteorological Society, 2013:408-412.

    [13] Varquez A C G, Nakayoshi M, Makabe T, et al. WRF application of high resolution urban surface parameters on some major cities of Japan [J]. Journal of Japan Society of Civil Engineers Ser B1 (Hydraulic Engineering), 2014,70(4):I_175-I_180.

    [14] Shen C, Shen A, Tian C Y, et al. Evaluating the impacts of updated aerodynamic roughness length in the WRF/Chem model over Pearl River Delta [J]. Meteorology and Atmospheric Physics, 2020,132: 427–440.

    [15] Myneni R B, Hoffman S, Knyazikhin Y, et al. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data [J]. Remote Sensing of Environment, 2002,83(1/2): 214-231.

    [16] Van Der Graaf S C, Kranenburg R, Segers A J, et al. Satellite-derived leaf area index and roughness length information for surface- atmosphere exchange modelling: a case study for reactive nitrogen deposition in north-western Europe using LOTOS-EUROS v2.0 [J]. Geoscientific Model Development, 2020,13(5):2451–2474.

    [17] 方紅亮.我國(guó)葉面積指數(shù)衛(wèi)星遙感產(chǎn)品生產(chǎn)及驗(yàn)證[J]. 遙感技術(shù)與應(yīng)用, 2020,35(5):990-1003.

    Fang H L. Development and validation of satellite Leaf Area Index (LAI) products in China [J]. Remote Sensing Technology and Application, 2020,35(5):990-1003.

    [18] 劉勇洪,房小怡,欒慶祖.基于衛(wèi)星數(shù)據(jù)與GIS技術(shù)的北京地區(qū)粗糙度長(zhǎng)度估算研究[J]. 高原氣象, 2016,35(6):1625-1638.

    Liu Y H, Fang X Y, Luan Q Z. Estimation of roughness length of Beijing area based on satellite data and GIS technique [J]. Plateau Meteorology, 2016,35(6):1625-1638.

    [19] Simard M, Pinto N, Fisher J B, et al. Mapping forest canopy height globally with spaceborne lidar [J]. Journal of Geophysical Research, 2011,116(G4).

    [20] Raupach M R. Drag and drag partition on rough surfaces [J]. Boundary-Layer Meteorol, 1992,60:375–395.

    [21] Raupach M R. Simplified expressions for vegetation roughness length and zero-plane displacement as functions of canopy height and area index [J]. Boundary-Layer Meteorol, 1994,71:211–216.

    [22] Jasinski M F, Borak J, Crago R. Bulk surface momentum parameters for satellite-derived vegetation fields [J]. Agricultural and Forest Meteorology, 2005,133(1-4):55-68.

    [23] Zeng X, Shaikh M, Dai Y, et al. Coupling of the common land model to the NCAR community climate model [J]. Journal of Climate, 2002,15:1832-1854.

    [24] Grell G A, Peckham S E, Schmitz R. Fully coupled "online" chemistry within the WRF model [J]. Atmospheric Environment, 2005,39(37): 6957–6975.

    [25] Seinfeld J H, Carmichael G R, Arimoto R, et al. Ace-Asia: regional climatic and atmospheric chemical effects of Asian dust and pollution [J]. Bulletin of the American Meteorological Society, 2004,85(3):367- 380.

    [26] Li M, Zhang Q, Kurokawa J I, et al. MIX: a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP [J]. Atmospheric Chemistry and Physics, 2015,17(2):935–963.

    [27] 于名召.空氣動(dòng)力學(xué)粗糙度的遙感方法及其在蒸散發(fā)計(jì)算中的應(yīng)用研究[D]. 北京:中國(guó)科學(xué)院大學(xué)(中國(guó)科學(xué)院遙感與數(shù)字地球研究所), 2018.

    Yu M Z. Research on remote sensing methods for the aerodynamic roughness length and its application in evapotranspiration calculation [D]. Beijing: University of Chinese Academy of Sciences(Institute of Remote Sensing and Digital Earth, Chinese ,Academy of Sciences), 2018.

    [28] 王國(guó)印.半干旱區(qū)陸氣相互作用的觀測(cè)與研究[D]. 蘭州:蘭州大學(xué), 2013.

    Wang G Y. Observation of land-atmosphere interaction over semi-arid regions [D]. Lanzhou: Lanzhou University, 2013.

    [29] 姚 彤.我國(guó)北方地區(qū)地表粗糙度和反照率參數(shù)化特征研究[D]. 蘭州:蘭州大學(xué), 2014.

    Yao T. A study on parametric feature of the surface roughness and albedo in north China [D]. Lanzhou: Lanzhou University, 2014.

    [30] 趙曉松,關(guān)德新,吳家兵,等.長(zhǎng)白山闊葉紅松林的零平面位移和粗糙度[J]. 生態(tài)學(xué)雜志, 2004,23(5):84-88.

    Zhao X S, Guan D X, Wu J B, et al. Zero-plane displacement and roughness length of the mixed forest of broad-leaved and Korean-pine in Changbai Mountain [J]. Chinese Journal of Ecology, 2004,23(5): 84-88.

    [31] 鐘 中,韓士杰.長(zhǎng)白山闊葉紅松林冠層空氣動(dòng)力學(xué)參數(shù)的計(jì)算[J]. 南京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2002,38(4):565-571.

    Zhong Z, Han S J. The calculation of aerodynamic parameters of Korea Pine canopy in Changbai Mountains [J]. Journal of Nanjing University(Natural Sciences), 2002,38(4):565-571.

    [32] 邢麗瑋.基于MODIS和GLAS數(shù)據(jù)反演多時(shí)序中國(guó)陸表植被空氣動(dòng)力學(xué)粗糙度[D]. 北京:首都師范大學(xué), 2012.

    Xing L W. Retrieving the multi-temporal aerodynamic roughness length of surface vegetation in China based on MODIS and GLAS data [D]. Beijing: Capital Normal University, 2012.

    [33] 曾 劍,張 強(qiáng).2008年夏季中國(guó)干旱-半干旱區(qū)陸面主要物理參數(shù)的平均特征[J]. 高原氣象, 2012,31(6):1539-1550.

    Zeng J, Zhang Q. Mean characteristics of land surface key parameters in semi-arid and regions of China in summers of 2008. Plateau meteorology, 2012,31(6):1539-1550.

    [34] 鞠英芹.不同下墊面類(lèi)型動(dòng)力學(xué)粗糙度與熱力學(xué)粗糙度的研究[D]. 南京:南京信息工程大學(xué), 2012.

    Ju Y Q. Study of aerodynamic roughness length and thermal roughness length on different underlying surfaces [D]. Nanjing: Nanjing University of Information Science & Technology, 2012.

    [35] 周艷蓮,孫曉敏,朱治林,等.幾種典型地表粗糙度計(jì)算方法的比較研究[J]. 地理研究, 2007,26(5):887-896.

    Zhou Y L, Sun X M, Zhu Z L, et al. Comparative research on four typical surface roughness length calculation methods [J]. Geographical Research, 2007,26(5):887-896.

    [36] De Bruin H A R, Moore C J. Zero-plane displacement and roughness length for tall vegetation, derived from a simple mass conservation hypothesis [J]. Boundary-Layer Meteorol, 1985,31:39–49.

    [37] 劉 偉,魏 信,石 文,等.復(fù)雜地形條件下零平面位移和空氣動(dòng)力學(xué)粗糙度的計(jì)算——以珠海南亞熱帶常綠闊葉林地區(qū)為例[J]. 熱帶氣象學(xué)報(bào), 2016,32(4):524-532.

    Liu W, Wei X, Shi W, et al. The calculation of zero-plane displacement and aerodynamic roughness of heterogenous terrain-- With the Evergreen broad-leaf forest in Zhuhai as an example [J]. Journal of tropical meteorology, 2016,32(4):524-532.

    [38] Jasinski M F, Crago R D. Estimation of vegetation aerodynamic roughness of natural regions using frontal area density determined from satellite imagery [J]. Agricultural and Forest Meteorology, 1999, 94(1):65-77.

    [39] Gallagher M W, Nemitz E, Dorsey J R, et al. Measurements and parameterizations of small aerosol deposition velocities to grassland, arable crops, and forest: Influence of surface roughness length on deposition [J]. Journal of Geophysical Research, 2002,107(D12):4154.

    [40] Gao Z Q, Wang J M, Ma Y M, et al. Study of roughness lengths and drag coefficients over Nansha sea region, Gobi, desert, Oasis and Tibetan plateau, Physics and Chemistry of the Earth, Part B: Hydrology [J]. Oceans and Atmosphere, 2000,25(2):141-145.

    [41] Lu L, Liu S, Xu Z. et al. The characteristics and parameterization of aerodynamic roughness length over heterogeneous surfaces [J]. Advances in Atmospheric Sciences, 2009,26:180-190.

    [42] Mailler S, Menut L, Khvorostyanov D, et al. CHIMERE-2017: From urban to hemispheric chemistry-transport modeling [J]. Geoscientific Model Development, 2017,10:2397-2423.

    [43] Silva J, Ribeiro C, Guedes R, et al. Roughness length classification of Corine land cover classes [J]. Proceedings of European Wind Energy Conference, 2007.

    [44] Simpson D, Benedictow A, Berge H, et al. The EMEP MSC-W chemical transport model – technical description [J]. Atmospheric Chemistry and Physics, 2012,12(16):7825–7865.

    [45] Wiernga J. Representative roughness parameters for homogeneous terrain [J]. Boundary-Layer Meteorology, 1993,63:323-363.

    [46] 馮健武,劉輝志,王 雷,等.半干旱區(qū)不同下墊面地表粗糙度和湍流通量整體輸送系數(shù)變化特征[J]. 中國(guó)科學(xué):地球科學(xué), 2012,42(1): 24-33.

    Feng J W, Liu H Z, Wang L, et al. Seasonal and inter-annual variation of surface roughness length and bulk transfer coefficients in a semiarid area. Sci China Earth Sci, 2012,42(1):24-33.

    [47] 高志球,王介民,馬耀明,等.不同下墊面的粗糙度和中性曳力系數(shù)研究[J]. 高原氣象, 2000,19(1):17-24.

    Gao Z Q, Wang J M, Ma Y M, et al. Study on roughness lengths and drag coefficients over the different underlying surfaces [J]. Plateau meteorology, 2000,19(1):17-24.

    [48] 李宏宇,張 強(qiáng),史晉森,等.黃土高原自然植被下墊面陸面過(guò)程參數(shù)研究[J]. 氣象學(xué)報(bào), 2012,70(5):1137-1148.

    Li H Y, Zhang Q, Shi J S, et al. A study of the parameterization of land-surface processes over the natural vegetation surface of Loess Plateau [J]. Acta Meteorologica Sinica, 2012,70(5):1137-1148.

    [49] 李鎖鎖,呂世華,柳媛普,等.黃河上游瑪曲地區(qū)空氣動(dòng)力學(xué)參數(shù)的確定及其在陸面過(guò)程模式中的應(yīng)用[J]. 高原氣象, 2010,29(6):1408- 1413.

    Li S S, Lü S H, Liu Y P, et al. Determination of aerodynamical parameter in Maqu Area in the upper reach of Yellow River and Its application in land surface process model [J]. Plateau meteorology, 2010,29(6):1408-1413.

    [50] 劉嘯然,李茂善,胡文斌.藏北高原那曲地區(qū)不同下墊面地表粗糙度的變化特征研究[J]. 高原氣象, 2019,38(2):428-438.

    Liu X R, Li M S, Hu W B. Variations of surface roughness on different underlying surface at Nagqu Area over the Qinghai-Tibetan Plateau [J]. Plateau Meteorology, 2019,38(2):428-438.

    [51] 馬耀明,塚本修,王介民,等.青藏高原草甸下墊面上的動(dòng)力學(xué)和熱力學(xué)參數(shù)分析[J]. 自然科學(xué)進(jìn)展, 2001,11(8):42-46.

    Ma Y M, Zhong B X, Wang J M, et al. Analysis of dynamic and thermodynamic parameters on the underlying surface of a meadow on the Tibetan Plateau [J]. Progress in Natural Science, 2001,11(8):42-46.

    [52] 孫 俊,胡澤勇,陳學(xué)龍,等.黑河中上游不同下墊面動(dòng)量總體輸送系數(shù)和地表粗糙度對(duì)比分析[J]. 高原氣象, 2012,31(4):920-926.

    Sun J, Hu Z Y, Chen X L, et al. Comparative analysis on momentum bulk transfer coefficients and roughness length under the different underlying surfaces in uppers and middle reaches of Heihe river basin [J]. Plateau Meteorology, 2012,31(4):920-926.

    [53] 張 強(qiáng),曾 劍,姚 桐.植被下墊面近地層大氣動(dòng)力狀態(tài)與動(dòng)力學(xué)粗糙度長(zhǎng)度的相互作用及其參數(shù)化關(guān)系[J]. 科學(xué)通報(bào), 2012,57(8): 647-655.

    Zhang Q, Zeng J, Yao T. Interaction of aerodynamic roughness length and windflow conditions and its parameterization over vegetation surface [J]. Chin Sci Bull, 2012,57(8):647-655.

    [54] 張文煜,張 宇,陸曉靜,等.黃土高原半干旱區(qū)非均一下墊面粗糙度分析[J]. 高原氣象, 2009,28(4):763-768.

    Zhang W Y, Zhang Y, Lu X J, et al. Analysis on heterogeneous underlying surface roughness length in semi-arid mountains area of Loess Plateau, China [J]. Plateau Meteorology, 2009,28(4):763-768.

    [55] Tsai J L, Tsuang B J. Aerodynamic roughness over an urban area and over two farmlands in a populated area as determined by wind profiles and surface energy flux measurements [J]. Agricultural and Forest Meteorology, 2005,132(1/2):154-170.

    [56] Tsai J L, Tsuang B J, Lu P S, et al. Measurements of aerodynamic roughness, Bowen ratio, and atmospheric surface layer height by eddy covariance and tethersonde systems simultaneously over a heterogeneous rice paddy [J]. Journal of Hydrometeorology, 2010, 11(2):452-466.

    [57] 賈 立,王介民.綠洲-沙漠復(fù)合地表?xiàng)l件下的局地和有效粗糙度[J]. 氣象學(xué)報(bào), 1999,57(3):91-102.

    Jia Li, Wang J M. The local and effective aerodynamic roughness length of a complex landscape of oasis and desert [J]. Acta meteorological sincia, 1999,57(3):91-102.

    [58] Lankreijer H J M, Hendriks M J, Klaassen W. A comparison of models simulating rainfall interception of forests [J]. Agricultural and Forest Meteorology, 1993,64(3/4):187-199.

    [59] Lo A K F. Determination of zero-plane displacement and roughness length of a forest canopy using profiles of limited height [J]. Boundary-Layer Meteorology, 1995,75:381-402.

    [60] Mikami M, Toya T, Yasuda N. An analytical method for the determination of the roughness parameters over complex regions [J]. Boundary-Layer Meteorology, 1996,79:23-33.

    [61] Parlange M B, Brutsaert W. Regional roughness of the Landes forest and surface shear stress under neutral conditions [J]. Boundary-Layer Meteorol, 1989,48:69–81.

    [62] Toda M, Sugita M. Single level turbulence measurements to determine roughness parameters of complex terrain [J]. Journal of Geophysical Research, 2003,108(D12):4363.

    [63] Weligepolage K, Gieske A S M, Van Der Tol C, et al. Effect of sub-layer corrections on the roughness parameterization of a Douglas fir forest [J]. Agricultural and Forest Meteorology, 2012,162-163: 115-126.

    [64] Yang R, Friedl M A. Determination of roughness lengths for heat and momentum over boreal forests [J]. Boundary-Layer Meteorology, 2003,107:581-603.

    [65] 劉和平,劉樹(shù)華,朱廷曜,等.森林冠層空氣動(dòng)力學(xué)參數(shù)的確定[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 1997,33(4):116-122.

    Liu H P, Liu S H, Zhu T Y, et al. Determination of aerodynamic parameters of Changbai mountain forest [J]. Acta Scicntiarum Naturalium, 1997,33(4):116-122.

    [66] 張學(xué)仕,蔣 琰,薛建輝,等.江蘇句容下蜀次生櫟林的空氣動(dòng)力學(xué)參數(shù)研究[J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010,34(6):61-65.

    Zhang X S, Jiang Y, Xue J H, et al. Aerodynamic parameters of secondary oak forest in Xiashu Jurong county of Jiangsu province [J]. Journal of Nanjing Forestry University (Natural Science Edition), 2020,34(6):61-65.

    [67] 張一平,宋清海,于貴瑞,等.西雙版納熱帶季節(jié)雨林風(fēng)時(shí)空變化特征初步分析[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2006,17(1):11-16.

    Zhang Y P, Song Q H, Yu G R, et al. Temporal-spatial characteristics of wind in tropical seasonal rainf orest in Xishuangbanna of Yunnan province [J]. Chinese Journal of Applied Ecology, 2006,17(1):11-16.

    [68] 周艷蓮,孫曉敏,朱治林,等.幾種不同下墊面地表粗糙度動(dòng)態(tài)變化及其對(duì)通量機(jī)理模型模擬的影響[J]. 中國(guó)科學(xué).D輯:地球科學(xué), 2006, 36(S1):244-254.

    Zhou Y L, Sun X M, Zhu Z L, et al.Dynamic variation of surface roughness on several different underlying surfaces and its influence on flux mechanism model simulation [J]. Science in China Series D:Earth Sciences,2006,36(S1):244-254.

    [69] Driese K L, Reiners W A. Aerodynamic roughness parameters for semi-arid natural shrub communities of Wyoming [J]. USA, Agricultural and Forest Meteorology, 1997,88(1-4):1-14.

    [70] 楊興華,艾力·買(mǎi)買(mǎi)提明,張瑞軍,等.塔克拉瑪干沙漠荒漠過(guò)渡帶空氣動(dòng)力學(xué)粗糙度分析[J]. 干旱區(qū)研究, 2012,29(3):524-528.

    Yang X H, Ali M, Zhang R J, et al. Analysis on aerodynamic roughness length in a desert transitional zone of the Taklimakan desert [J]. Arid Zone Research, 2012,29(3):524-528.

    [71] 石雪峰,夏建新,吉祖穩(wěn).空氣動(dòng)力學(xué)粗糙度與植被特征關(guān)系的研究進(jìn)展[J]. 中央民族大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006,15(3):218-225.

    Shi X F, Xia J X, Ji Z W. A review:the study in relationship between aerodynamic roughness length and the characteristics of vegetation [J]. 2006,15(3):218-225.

    [72] 藍(lán) 靜.珠三角區(qū)域氣溶膠污染及低能見(jiàn)度過(guò)程的數(shù)值模擬[D]. 廣州:中山大學(xué), 2013.

    Lan J. Numerical simulation of aerosol contamination and low visibility processes in the Pearl River Delta region [D]. Guangzhou: Sun Yat-sen University, 2013.

    [73] 劉新峰,袁 惠,楊 軍,等.統(tǒng)計(jì)分析方法在大氣環(huán)境監(jiān)測(cè)數(shù)據(jù)符合性分析中的應(yīng)用探討[J]. 四川環(huán)境, 2012,31(2):36-39.

    Liu X F, Yuan H, Yang J, et al. Application of statistical analysis for data compliance analysis of atmosphere environmental monitoring [J]. Sichuan Environment, 2012,31(2):36-39.

    Estimation of vegetation underlying surface0and assessment on the impacts of its improvement.

    TIAN Chun-yan1,2, CUI Yin-ping1, SHEN Chong3, CHEN Xiu-zhi1, SHEN Ao1, FAN Qi1*

    (1.Guangdong Province Key Laboratory for Climate Change and Natural Disaster Studies, Southern Laboratory of Ocean Science and Engineering (Zhuhai), School of Atmospheric Sciences, Sun Yat-sen University, Zhuhai 519082, China;2.Meteorological Center of Middle South Regional Air Traffic Management Bureau of CAAC, Guangzhou 510000, China;3.Guangzhou Climate and Agrometeorology Center, Guangzhou 511430, China)., 2022,42(9):3969~3982

    The morphology-remote sensing inversion with satellite remote sensing and multi-observation data were utilized to estimate the aerodynamic roughness length (0) of main vegetation underlying surfaces in China. This new0dataset was then used to replace the default values in the WRF-Chem model. Combined with the national meteorological and pollutant observation data, the impacts of0improvement on the meteorological field and chemical field simulations were evaluated. The results show that the high value of0is mainly located over the forest underlying surface (higher than 1m); the0values over the farmland and grassland underlying surfaces are generally lower (lower than 0.5m). The0values over other underlying surfaces fell in between 0.5m and 1m; and the modification of0has a better improvement in the simulations of 10m wind speed (WS10) and surface temperature (TSK) even though its effects on 2m temperature (2) and relative humidity (RH) are not obvious. We also observed that an increase in0mainly has a warming effect on TSK and a weakening effect on WS10 for different types of underlying surfaces. In terms of the improvement ratio (the difference between the simulated value after the improvement of each element and the simulated value before the improvement divided by the simulated value before the improvement, the same below), the improved0demonstrates a greater impact on WS10 than on TSK. For the pollutant concentrations, the modification of0has larger effects on the simulations for PM2.5and PM10than for other pollutants (SO2, NO2, O3). Furthermore, the impact of the improved0was larger on meteorological elements than on pollutant concentrations. In general, the improvement of0has the largest impact on the 10m wind speed; because the improved0is generally higher than the model-default0value. As0is a variable for characterizing the surface roughness, the obstruction to the air flow increases with an increase in surface roughness and thus results in a reduction in the near surface wind speed.

    aerodynamic roughness length;vegetation underlying surface;numerical simulation;WRF-Chem model;PM2.5

    X171

    A

    1000-6923(2022)09-3969-14

    2022-02-21

    國(guó)家重點(diǎn)研發(fā)計(jì)劃(2019YFC0214605);國(guó)家自然科學(xué)基金資助項(xiàng)目(42075181);廣東省科技計(jì)劃項(xiàng)目(科技創(chuàng)新平臺(tái)類(lèi))(2019B121201002);廣東省氣候變化與自然災(zāi)害研究重點(diǎn)實(shí)驗(yàn)室(2020B1212060025)

    *責(zé)任作者, 教授, eesfq@mail.sysu.edu.cn

    田春艷(1994-),女,重慶人,中山大學(xué)碩士研究生,主要大氣環(huán)境數(shù)值模擬研究.發(fā)表論文1篇.

    猜你喜歡
    下墊面粗糙度風(fēng)速
    不同下墊面對(duì)氣溫的影響
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測(cè)
    基于最優(yōu)TS評(píng)分和頻率匹配的江蘇近海風(fēng)速訂正
    基于無(wú)人機(jī)影像的巖體結(jié)構(gòu)面粗糙度獲取
    甘肅科技(2020年20期)2020-04-13 00:30:18
    冷沖模磨削表面粗糙度的加工試驗(yàn)與應(yīng)用
    模具制造(2019年4期)2019-06-24 03:36:48
    北京與成都城市下墊面閃電時(shí)空分布特征對(duì)比研究
    流域下墊面變化對(duì)潮白河密云水庫(kù)上游徑流影響分析
    下墊面變化對(duì)徑流及洪水影響分析
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    鋼材銹蝕率與表面三維粗糙度參數(shù)的關(guān)系
    亚洲内射少妇av| 一级毛片我不卡| 久久人人爽人人片av| 亚洲在线自拍视频| 国内精品美女久久久久久| 26uuu在线亚洲综合色| 欧美性猛交╳xxx乱大交人| 亚洲欧美中文字幕日韩二区| 乱码一卡2卡4卡精品| 亚洲在久久综合| 亚洲精品乱码久久久久久按摩| 99久久精品热视频| 97热精品久久久久久| 极品教师在线视频| 黄色一级大片看看| 婷婷色综合大香蕉| 99国产精品一区二区蜜桃av| 免费黄网站久久成人精品| 成人永久免费在线观看视频| 国产三级中文精品| 狂野欧美激情性xxxx在线观看| av女优亚洲男人天堂| 亚洲人成网站高清观看| 麻豆久久精品国产亚洲av| 在线免费观看的www视频| 嫩草影院入口| 天堂影院成人在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产极品天堂在线| 欧美性猛交黑人性爽| 国产淫片久久久久久久久| 国产色爽女视频免费观看| 变态另类成人亚洲欧美熟女| 国产私拍福利视频在线观看| 给我免费播放毛片高清在线观看| 啦啦啦观看免费观看视频高清| 日本成人三级电影网站| 久久亚洲精品不卡| 综合色av麻豆| 卡戴珊不雅视频在线播放| 亚洲va在线va天堂va国产| 欧美xxxx黑人xx丫x性爽| 夜夜爽天天搞| 99riav亚洲国产免费| 一级毛片我不卡| 美女高潮的动态| 亚洲精品456在线播放app| 亚洲天堂国产精品一区在线| 高清毛片免费观看视频网站| 乱人视频在线观看| 亚洲欧美成人综合另类久久久 | 身体一侧抽搐| 亚洲真实伦在线观看| 26uuu在线亚洲综合色| 一区福利在线观看| av又黄又爽大尺度在线免费看 | 国产国拍精品亚洲av在线观看| 久久久精品欧美日韩精品| 六月丁香七月| 最近2019中文字幕mv第一页| 91精品国产九色| 国产色婷婷99| 国产成人freesex在线| 国产成人精品久久久久久| 黄片无遮挡物在线观看| 五月伊人婷婷丁香| 亚洲欧美成人综合另类久久久 | 日本在线视频免费播放| 免费观看在线日韩| 亚洲国产欧美人成| 国产一级毛片在线| 国产高清不卡午夜福利| 长腿黑丝高跟| 三级毛片av免费| 午夜免费男女啪啪视频观看| 久久国产乱子免费精品| 中国美女看黄片| 国产乱人偷精品视频| 欧美区成人在线视频| 免费电影在线观看免费观看| 狂野欧美白嫩少妇大欣赏| 小蜜桃在线观看免费完整版高清| 18禁黄网站禁片免费观看直播| 搡老妇女老女人老熟妇| 精品人妻一区二区三区麻豆| 成人毛片a级毛片在线播放| 色哟哟·www| 熟女人妻精品中文字幕| 日本熟妇午夜| 高清日韩中文字幕在线| 可以在线观看的亚洲视频| 欧美变态另类bdsm刘玥| av国产免费在线观看| 两个人视频免费观看高清| 亚洲av电影不卡..在线观看| 午夜福利在线观看免费完整高清在 | 最近2019中文字幕mv第一页| 天天一区二区日本电影三级| 色吧在线观看| 久久精品夜色国产| 色综合色国产| 精品午夜福利在线看| 成人性生交大片免费视频hd| 一个人看视频在线观看www免费| 97人妻精品一区二区三区麻豆| 青青草视频在线视频观看| 国产精品久久久久久久久免| 国产熟女欧美一区二区| a级毛片免费高清观看在线播放| 少妇熟女欧美另类| 国语自产精品视频在线第100页| 久久久久久久久大av| 男人和女人高潮做爰伦理| 美女被艹到高潮喷水动态| 日日啪夜夜撸| 国产精品久久久久久av不卡| 99精品在免费线老司机午夜| 91av网一区二区| 99久久无色码亚洲精品果冻| 欧美最黄视频在线播放免费| 99国产极品粉嫩在线观看| 在线观看美女被高潮喷水网站| 国产一级毛片在线| 亚洲第一区二区三区不卡| 成年女人永久免费观看视频| 女的被弄到高潮叫床怎么办| 国内精品美女久久久久久| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 免费av观看视频| 身体一侧抽搐| 少妇被粗大猛烈的视频| 精品一区二区免费观看| 亚洲av熟女| 免费在线观看成人毛片| 中出人妻视频一区二区| 黄色欧美视频在线观看| 国产黄片美女视频| 亚洲精品日韩在线中文字幕 | 成年女人永久免费观看视频| 麻豆国产97在线/欧美| 欧美日韩乱码在线| 精品无人区乱码1区二区| 中文在线观看免费www的网站| 成人高潮视频无遮挡免费网站| 久久欧美精品欧美久久欧美| 我的女老师完整版在线观看| 色哟哟·www| 好男人在线观看高清免费视频| 51国产日韩欧美| 亚洲av电影不卡..在线观看| 观看免费一级毛片| 国产精品野战在线观看| 亚洲色图av天堂| 日韩欧美精品v在线| 99在线人妻在线中文字幕| 91久久精品国产一区二区成人| 亚洲四区av| 91精品国产九色| 成人欧美大片| 国产精品久久久久久久电影| 少妇被粗大猛烈的视频| 噜噜噜噜噜久久久久久91| av女优亚洲男人天堂| 男女那种视频在线观看| 婷婷精品国产亚洲av| 3wmmmm亚洲av在线观看| 亚洲av一区综合| 久久精品国产亚洲网站| 欧美最黄视频在线播放免费| 1000部很黄的大片| 国产精品一区二区三区四区久久| 国产黄色视频一区二区在线观看 | 欧美+日韩+精品| 悠悠久久av| 97在线视频观看| 99久国产av精品| 久久精品91蜜桃| 综合色av麻豆| 久久久久久久亚洲中文字幕| 国产午夜精品久久久久久一区二区三区| 亚洲婷婷狠狠爱综合网| 亚洲性久久影院| 欧美一区二区国产精品久久精品| 精品人妻熟女av久视频| 久久久久国产网址| av又黄又爽大尺度在线免费看 | 国产一区亚洲一区在线观看| 97在线视频观看| 成人午夜高清在线视频| 国产老妇伦熟女老妇高清| 亚洲图色成人| 成年版毛片免费区| 在线观看一区二区三区| 91在线精品国自产拍蜜月| 成人鲁丝片一二三区免费| 国产伦一二天堂av在线观看| 2021天堂中文幕一二区在线观| 淫秽高清视频在线观看| 国产麻豆成人av免费视频| 亚洲av中文字字幕乱码综合| 一进一出抽搐gif免费好疼| 精品久久国产蜜桃| 中文字幕精品亚洲无线码一区| 国产伦在线观看视频一区| 99久久中文字幕三级久久日本| 51国产日韩欧美| kizo精华| 亚洲精品自拍成人| 国产精品人妻久久久影院| 岛国在线免费视频观看| 亚洲一区高清亚洲精品| 亚洲人成网站在线观看播放| 亚洲一区二区三区色噜噜| 久久亚洲精品不卡| 91久久精品国产一区二区三区| 国产片特级美女逼逼视频| 人妻久久中文字幕网| 亚洲在线观看片| 白带黄色成豆腐渣| 热99在线观看视频| 精品免费久久久久久久清纯| 国产视频内射| 在线播放国产精品三级| 久久久久久久久久成人| 亚洲精品久久久久久婷婷小说 | 欧美日韩在线观看h| 一级二级三级毛片免费看| 亚洲一区高清亚洲精品| 亚洲欧美日韩东京热| 少妇被粗大猛烈的视频| 免费搜索国产男女视频| 亚洲美女搞黄在线观看| 国产精品一区二区三区四区免费观看| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月| 99久久人妻综合| 久久精品人妻少妇| 亚洲最大成人手机在线| 夫妻性生交免费视频一级片| av视频在线观看入口| 成人一区二区视频在线观看| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 成人毛片a级毛片在线播放| 免费人成视频x8x8入口观看| 日韩在线高清观看一区二区三区| 长腿黑丝高跟| 日本一本二区三区精品| 日韩制服骚丝袜av| 欧美极品一区二区三区四区| 国产精品福利在线免费观看| 国产中年淑女户外野战色| 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩卡通动漫| 亚洲,欧美,日韩| 精品久久国产蜜桃| 国产色婷婷99| 亚州av有码| 国产v大片淫在线免费观看| 日本免费一区二区三区高清不卡| 国产亚洲精品久久久com| 看免费成人av毛片| 在线观看美女被高潮喷水网站| 欧美极品一区二区三区四区| 亚洲国产高清在线一区二区三| av又黄又爽大尺度在线免费看 | 丝袜美腿在线中文| 如何舔出高潮| 免费av不卡在线播放| 亚洲av.av天堂| 免费观看的影片在线观看| 亚洲第一电影网av| 国产探花极品一区二区| 日韩中字成人| 亚洲高清免费不卡视频| 成人美女网站在线观看视频| 在线a可以看的网站| 亚洲欧美成人综合另类久久久 | 日韩一区二区三区影片| 精品一区二区三区人妻视频| 性色avwww在线观看| 亚洲最大成人中文| 亚洲国产欧美人成| 美女被艹到高潮喷水动态| 久久中文看片网| 成人亚洲精品av一区二区| 日产精品乱码卡一卡2卡三| 男人和女人高潮做爰伦理| 波野结衣二区三区在线| 久久韩国三级中文字幕| 欧美人与善性xxx| 丰满乱子伦码专区| 国产欧美日韩精品一区二区| 99在线人妻在线中文字幕| 国产亚洲av片在线观看秒播厂 | 一个人看的www免费观看视频| 欧美性感艳星| 日韩成人伦理影院| 国产一区二区三区在线臀色熟女| 久久久欧美国产精品| 亚洲国产欧洲综合997久久,| 99热这里只有是精品50| 免费不卡的大黄色大毛片视频在线观看 | 嫩草影院精品99| 成人综合一区亚洲| 国产亚洲91精品色在线| 大又大粗又爽又黄少妇毛片口| 久久久成人免费电影| 亚洲欧美精品综合久久99| 日本五十路高清| 一进一出抽搐动态| 中文字幕av成人在线电影| 99久国产av精品国产电影| 国产精品1区2区在线观看.| 一区二区三区四区激情视频 | 成人无遮挡网站| 成人高潮视频无遮挡免费网站| 99九九线精品视频在线观看视频| 中文资源天堂在线| 欧美性感艳星| 搞女人的毛片| av在线播放精品| 亚洲性久久影院| 激情 狠狠 欧美| 日韩一区二区视频免费看| 嫩草影院入口| 岛国毛片在线播放| 99久久无色码亚洲精品果冻| 亚洲经典国产精华液单| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品合色在线| 久久6这里有精品| 国产成人freesex在线| 久久综合国产亚洲精品| 特级一级黄色大片| 嫩草影院新地址| av在线亚洲专区| 久久久久久久午夜电影| 赤兔流量卡办理| 有码 亚洲区| 欧美最黄视频在线播放免费| 亚洲国产欧美人成| 97超碰精品成人国产| 精品久久国产蜜桃| 久久久久久久久中文| 麻豆国产97在线/欧美| 亚洲三级黄色毛片| 国产精品日韩av在线免费观看| 青春草视频在线免费观看| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 国产私拍福利视频在线观看| 亚洲av成人精品一区久久| 日韩制服骚丝袜av| 老女人水多毛片| 久久精品国产清高在天天线| 亚洲人成网站高清观看| 少妇的逼好多水| 日韩视频在线欧美| 97在线视频观看| 久久这里有精品视频免费| 一边亲一边摸免费视频| 日韩一区二区视频免费看| 精品久久国产蜜桃| 免费av观看视频| 激情 狠狠 欧美| 97超碰精品成人国产| 亚洲最大成人av| 亚洲欧美成人精品一区二区| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 午夜亚洲福利在线播放| 亚洲欧美日韩东京热| 狠狠狠狠99中文字幕| 国内精品一区二区在线观看| 亚洲高清免费不卡视频| 成人特级黄色片久久久久久久| 日韩高清综合在线| 日本-黄色视频高清免费观看| 国产精华一区二区三区| 国产女主播在线喷水免费视频网站 | 亚洲自偷自拍三级| 可以在线观看毛片的网站| 有码 亚洲区| 波多野结衣巨乳人妻| 久久久久久大精品| 搞女人的毛片| 美女内射精品一级片tv| 日本-黄色视频高清免费观看| 亚洲av免费高清在线观看| 成人毛片60女人毛片免费| 国产人妻一区二区三区在| 久久久久久九九精品二区国产| 国产成人aa在线观看| 国产精品久久久久久亚洲av鲁大| 久久久国产成人精品二区| 深夜a级毛片| 亚洲av电影不卡..在线观看| 久久精品国产清高在天天线| 国产探花在线观看一区二区| 久久久精品94久久精品| 久久亚洲国产成人精品v| 日本免费a在线| 一个人看视频在线观看www免费| 成人午夜高清在线视频| 在线观看免费视频日本深夜| 亚洲欧美日韩高清专用| 波多野结衣高清无吗| 麻豆av噜噜一区二区三区| 精华霜和精华液先用哪个| 欧美高清性xxxxhd video| 久久久久久伊人网av| 国产成人精品一,二区 | 网址你懂的国产日韩在线| 午夜老司机福利剧场| 国产亚洲精品久久久久久毛片| 给我免费播放毛片高清在线观看| 欧美+亚洲+日韩+国产| 日本黄色片子视频| 亚洲美女视频黄频| 成人美女网站在线观看视频| 淫秽高清视频在线观看| 国内精品一区二区在线观看| 亚洲va在线va天堂va国产| 成人毛片60女人毛片免费| 美女 人体艺术 gogo| 欧美日本亚洲视频在线播放| 国产91av在线免费观看| 国产精品国产三级国产av玫瑰| 精品久久久噜噜| 一级毛片电影观看 | av福利片在线观看| 好男人在线观看高清免费视频| 日韩av不卡免费在线播放| 美女高潮的动态| 日韩精品有码人妻一区| 国产在线精品亚洲第一网站| 色哟哟·www| 国产伦一二天堂av在线观看| 亚洲婷婷狠狠爱综合网| 国产高潮美女av| 麻豆乱淫一区二区| 欧美日韩在线观看h| 日韩高清综合在线| 亚洲第一区二区三区不卡| 国产精品久久久久久精品电影小说 | 不卡视频在线观看欧美| 免费观看人在逋| 亚洲经典国产精华液单| 日本欧美国产在线视频| 亚洲精品乱码久久久久久按摩| 性欧美人与动物交配| 26uuu在线亚洲综合色| 国产又黄又爽又无遮挡在线| 成人特级av手机在线观看| 色5月婷婷丁香| 在线观看美女被高潮喷水网站| 欧美高清成人免费视频www| 精品久久国产蜜桃| 欧美一级a爱片免费观看看| 亚洲一级一片aⅴ在线观看| 高清毛片免费看| 26uuu在线亚洲综合色| 极品教师在线视频| 人人妻人人澡欧美一区二区| 秋霞在线观看毛片| 国产精品综合久久久久久久免费| 国产片特级美女逼逼视频| 人体艺术视频欧美日本| 卡戴珊不雅视频在线播放| 少妇熟女aⅴ在线视频| 久久精品综合一区二区三区| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| 寂寞人妻少妇视频99o| 色播亚洲综合网| 欧美人与善性xxx| 在线免费观看的www视频| 国产精品一区www在线观看| 婷婷精品国产亚洲av| 国产视频首页在线观看| 欧美色欧美亚洲另类二区| av.在线天堂| 亚洲四区av| 日韩成人伦理影院| 18禁在线播放成人免费| 日韩欧美在线乱码| 深夜精品福利| 亚洲在线自拍视频| 中文精品一卡2卡3卡4更新| 精品人妻一区二区三区麻豆| 少妇被粗大猛烈的视频| 国产单亲对白刺激| 黄片无遮挡物在线观看| 日本与韩国留学比较| 国产麻豆成人av免费视频| 天天躁夜夜躁狠狠久久av| 能在线免费观看的黄片| 永久网站在线| 久久久色成人| 国产免费一级a男人的天堂| 国产淫片久久久久久久久| 我的老师免费观看完整版| 一进一出抽搐动态| 亚洲欧美成人综合另类久久久 | 在线观看免费视频日本深夜| 免费在线观看成人毛片| 国产成人91sexporn| 亚洲av中文字字幕乱码综合| 免费人成视频x8x8入口观看| 久久久久久九九精品二区国产| 老司机福利观看| 国产成人精品婷婷| 成人高潮视频无遮挡免费网站| 亚洲av免费高清在线观看| 亚洲成人精品中文字幕电影| 亚洲图色成人| a级毛片a级免费在线| 国语自产精品视频在线第100页| 在线免费十八禁| 欧美日韩在线观看h| 丝袜喷水一区| 最近视频中文字幕2019在线8| .国产精品久久| 在线免费观看不下载黄p国产| 搡女人真爽免费视频火全软件| 国产精品电影一区二区三区| 亚洲国产精品成人综合色| 国产高潮美女av| 午夜福利成人在线免费观看| 国产黄a三级三级三级人| 成人高潮视频无遮挡免费网站| 好男人在线观看高清免费视频| 免费观看a级毛片全部| 又黄又爽又刺激的免费视频.| 99九九线精品视频在线观看视频| 中文亚洲av片在线观看爽| 麻豆乱淫一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲综合色惰| 尾随美女入室| 亚洲四区av| 亚洲精品成人久久久久久| 久久久精品大字幕| 欧美zozozo另类| 亚洲第一区二区三区不卡| 久久精品国产99精品国产亚洲性色| 99热只有精品国产| 亚洲,欧美,日韩| 午夜精品在线福利| 日韩中字成人| 精品久久久久久久人妻蜜臀av| 精品久久久久久久末码| av天堂在线播放| 国产大屁股一区二区在线视频| 国产男人的电影天堂91| 成人午夜精彩视频在线观看| 欧美一区二区亚洲| 大香蕉久久网| 能在线免费看毛片的网站| 亚洲中文字幕一区二区三区有码在线看| 最近的中文字幕免费完整| 日日啪夜夜撸| 极品教师在线视频| 热99在线观看视频| 亚洲图色成人| 国产精品一区二区三区四区久久| 国产高清三级在线| 国产成人福利小说| 久久亚洲国产成人精品v| 免费看美女性在线毛片视频| 麻豆成人av视频| 卡戴珊不雅视频在线播放| 欧美日韩一区二区视频在线观看视频在线 | 国产午夜精品久久久久久一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 男女下面进入的视频免费午夜| 久久精品国产鲁丝片午夜精品| 国产精品人妻久久久影院| 欧美变态另类bdsm刘玥| 成人漫画全彩无遮挡| 欧美日韩综合久久久久久| 女人被狂操c到高潮| 亚洲国产欧洲综合997久久,| 欧美不卡视频在线免费观看| 99热6这里只有精品| 亚州av有码| 国产精品国产高清国产av| 精品无人区乱码1区二区| 亚洲av免费在线观看| 婷婷亚洲欧美| 非洲黑人性xxxx精品又粗又长| 在线免费观看的www视频| 日韩国内少妇激情av| 熟女人妻精品中文字幕| 国产精品一及| 黄色配什么色好看| 少妇熟女欧美另类| 国产精品一区二区三区四区久久| 亚洲色图av天堂| 麻豆精品久久久久久蜜桃| 欧美精品国产亚洲| 人人妻人人澡人人爽人人夜夜 | 精品免费久久久久久久清纯| 久久热精品热| 国产综合懂色| av视频在线观看入口| 日日啪夜夜撸| 中文字幕久久专区| 欧美xxxx黑人xx丫x性爽| 成人三级黄色视频| 麻豆一二三区av精品| 国产亚洲精品久久久com|