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

    北京市平原百萬畝大造林工程降溫效應(yīng)及其價(jià)值的遙感分析

    2017-04-10 12:28:06賈寶全仇寬彪
    生態(tài)學(xué)報(bào) 2017年3期
    關(guān)鍵詞:冷島亮溫平原

    賈寶全,仇寬彪

    1 中國林業(yè)科學(xué)研究院林業(yè)研究所, 北京 100091 2 國家林業(yè)局林木培育重點(diǎn)實(shí)驗(yàn)室, 北京 100091 3 國家林業(yè)局城市森林研究中心, 北京 100091 4 北京林業(yè)大學(xué)水土保持與荒漠化防治教育部重點(diǎn)實(shí)驗(yàn)室,水土保持學(xué)院, 北京 100083

    北京市平原百萬畝大造林工程降溫效應(yīng)及其價(jià)值的遙感分析

    賈寶全1,2,3,*,仇寬彪1,4

    1 中國林業(yè)科學(xué)研究院林業(yè)研究所, 北京 100091 2 國家林業(yè)局林木培育重點(diǎn)實(shí)驗(yàn)室, 北京 100091 3 國家林業(yè)局城市森林研究中心, 北京 100091 4 北京林業(yè)大學(xué)水土保持與荒漠化防治教育部重點(diǎn)實(shí)驗(yàn)室,水土保持學(xué)院, 北京 100083

    隨著城市化進(jìn)程的加快,城市熱島效應(yīng)已經(jīng)嚴(yán)重影響到了城市居民的生產(chǎn)與生活。大量的研究成果證明,以植被和水體為主導(dǎo)的冷島效應(yīng)是減弱城市熱島效應(yīng)的最有效、最持久、最經(jīng)濟(jì)手段,但過往的研究成果主要以城市公園作為研究對象,樣本數(shù)量少,研究結(jié)果的不確定性大,對于大范圍人工林地的冷島效應(yīng)關(guān)注較少。北京市在2012—2014年間實(shí)施了“百萬畝平原大造林”工程,累計(jì)形成10648個(gè)林地斑塊,這為人工林地冷島效應(yīng)問題的進(jìn)一步研究提供了絕好的研究樣本。利用2014年9月4日的landsat- 8衛(wèi)星影像數(shù)據(jù),通過定量遙感手段反演了北京市平原區(qū)的地表亮溫,并以造林地斑塊的GIS數(shù)據(jù)為基礎(chǔ),對造林工程的降溫效應(yīng)進(jìn)行了分析。結(jié)果表明,雖然造林時(shí)間短,但由于采用了大苗造林措施、公司實(shí)施模式,這些造林地的平均降溫幅度可達(dá)1.023℃,其冷島效應(yīng)的輻射范圍可到達(dá)林地邊緣外350m,其中0—100m距離范圍內(nèi)的降溫效果最為顯著,達(dá)到了0.392℃。從降溫效應(yīng)的區(qū)域差異來看,延慶盆地造林地塊的降溫效果最顯著,達(dá)到了3.519℃,6環(huán)以內(nèi)造林地塊的降溫效果最小,只有0.111℃。通過不同年份造林地塊的降溫效應(yīng)分析發(fā)現(xiàn),造林地塊的景觀格局和濕地保護(hù)與建設(shè)工程類型的規(guī)模對林地斑塊的降溫效果有顯著影響。從全部林地斑塊的大小看,其呈現(xiàn)了數(shù)量上以小斑塊占優(yōu)、面積上以中、小型斑塊為主的數(shù)量特征,降溫效果則呈現(xiàn)了斑塊規(guī)模愈大降溫效果愈顯著的特點(diǎn)。初步的價(jià)值估算表明,平原區(qū)大造林工程的降溫效益達(dá)到了4.8882×108元,其中林地本身的降溫效應(yīng)占到了53.73%,林地外圍輻射降溫的效益占到了46.38%。在今后的造林工程實(shí)踐中,要注意造林地塊大小的設(shè)計(jì)問題,應(yīng)以超大斑塊和巨斑塊為主體,同時(shí)要注意將林地斑塊的空間鄰接距離保持在100—500m的范圍內(nèi),以充分利用林地斑塊冷島效應(yīng)空間上的輻射擴(kuò)散特點(diǎn)。

    造林工程;降溫效應(yīng);LANDSAT- 8;價(jià)值評估;平原區(qū);北京市

    城市是人類文明的標(biāo)志,是一個(gè)時(shí)代經(jīng)濟(jì)、政治、社會(huì)、科學(xué)、文化、生態(tài)環(huán)境發(fā)展和變化的焦點(diǎn)和結(jié)晶體[1]。目前全球性的城市化進(jìn)程,既影響著生物多樣性和生態(tài)系統(tǒng)的功能,同時(shí)也影響著地方與區(qū)域氣候,以及城市居民的生活質(zhì)量[2]。城市熱島效應(yīng)作為影響城市生態(tài)環(huán)境質(zhì)量的重要因素,正吸引了越來越多學(xué)者的關(guān)注[3]。作為城市有生命的基礎(chǔ)設(shè)施,綠地以其特殊的“冷島效應(yīng)”,在減緩熱島效應(yīng)方面發(fā)揮著不可替代的作用[4- 5]。長久以來,圍繞公園綠地的“冷島效應(yīng)”,國內(nèi)外學(xué)者借助基于地面小氣候觀測到現(xiàn)代先進(jìn)的遙感技術(shù)等不同的分析研究手段,探討了公園面積、周長、幾何形狀等外部特征[6- 8],以及公園植被組成特征與垂直空間配置[5,9- 10]、公園內(nèi)部綠地與水體的結(jié)構(gòu)比例[11]等對其降溫效應(yīng)的影響,同時(shí)也對冷島效應(yīng)的季節(jié)規(guī)律[12- 14]、日變化規(guī)律[15- 16]等進(jìn)行了探討。近年來,也有學(xué)者利用景觀生態(tài)學(xué)的理論和方法,對公園的景觀空間結(jié)構(gòu)參數(shù)與城市熱場的關(guān)系進(jìn)行了卓有成效的分析[17- 20],另外也有部分學(xué)者對城市水域景觀的熱環(huán)境效應(yīng)、城市公園中植被與水體結(jié)合度對公園“冷島效應(yīng)”的影響[2,21- 23]等開展了一些研究工作。從已有的研究成果看,過去的相關(guān)工作都聚焦在城市內(nèi)部城市公園或其他類型的城市綠地上,在城市不斷擴(kuò)展,生態(tài)用地與其它用地矛盾日漸加劇的現(xiàn)實(shí)狀況下,這些研究工作對于揭示植被冷島效應(yīng)的形成機(jī)理、影響因素等科學(xué)問題是很有成效的,但在我國目前城市化加速、城市外圍鄉(xiāng)村土地轉(zhuǎn)化加快的背景下,這些研究成果能否有效指導(dǎo)城市外圍地帶的生態(tài)建設(shè),目前還沒有更多的科學(xué)根據(jù)來支撐;同時(shí),這些工作的研究結(jié)論多是以個(gè)別公園或不多的公園組合為研究對象得出的,目前所見的最多公園樣本數(shù)為24個(gè)[19],由于受樣本數(shù)量的制約,研究結(jié)果具有一定的不確定性。

    北京是我國乃至世界的著名大都市,占全市面積約1/3的平原區(qū)域承載了絕大部分的首都功能,但隨著人口與經(jīng)濟(jì)的高密度聚集,以及中心城區(qū)的持續(xù)性外擴(kuò),該區(qū)域面臨的生態(tài)與環(huán)境問題日趨嚴(yán)重。為了加強(qiáng)北京市的大氣治理、改善區(qū)域生態(tài)環(huán)境,根據(jù)平原區(qū)林少、生態(tài)與環(huán)境問題嚴(yán)重的現(xiàn)實(shí),北京市委、市政府規(guī)劃、啟動(dòng)了平原區(qū)造林工程,規(guī)劃在2012—2014年的3a時(shí)間內(nèi),在平原區(qū)新增森林面積66666.7hm2。截止2014年年底,全市已完成平原造林68050.98hm2,植樹5000余萬株,平原地區(qū)的森林覆蓋率凈增9.65個(gè)百分點(diǎn)。因其造林規(guī)模巨大、分布區(qū)又位于中心城區(qū)以外的廣大城鄉(xiāng)平原,這給我們深入探討城市森林綠地“冷島效應(yīng)”的相關(guān)問題,提供了絕好的研究樣本。

    1 研究數(shù)據(jù)與研究方法

    1.1 衛(wèi)星影像數(shù)據(jù)

    北京市域面積為1.641×104km2,因此選擇中尺度的Landsat TM衛(wèi)星影像作為本次研究工作的唯一信息源(衛(wèi)星影像數(shù)據(jù)軌道號(hào)分別為123/32、123/33)。由于北京市的造林工程以春、秋兩季造林施工為主,所以在衛(wèi)星影像的時(shí)間上,選擇了2014年9月4日的Landsat- 8 衛(wèi)星影像作為基本的分析數(shù)據(jù)源。自美國地質(zhì)調(diào)查局(USGS)網(wǎng)站(http://glovis.usgs.gov/)下載了相關(guān)分幅的L1T級(jí)影像數(shù)據(jù)之后,再從中國科學(xué)院遙感與數(shù)字地球研究所的對地觀測數(shù)據(jù)共享服務(wù)網(wǎng)(http://www.geodata.cn/)上下載同期同景數(shù)據(jù)(該網(wǎng)上數(shù)據(jù)的都是進(jìn)行過正射校正的L4級(jí)產(chǎn)品,但因缺乏band10和band11故無法直接用來做熱場反演)對其進(jìn)行了幾何校正,校正精度保持在1個(gè)像元之內(nèi),之后再在ERDAS2014軟件平臺(tái)上的ATCOR2模塊下對影像進(jìn)行大氣校正。

    1.2 平原造林?jǐn)?shù)據(jù)

    平原造林地塊數(shù)據(jù)來源于北京市林業(yè)勘察設(shè)計(jì)研究院,該數(shù)據(jù)是在航片基礎(chǔ)上通過野外調(diào)繪而成(圖1)。根據(jù)GIS統(tǒng)計(jì),3a來新造林圖斑共10648個(gè),最小的造林地斑塊面積5.25m2,最大的林地斑塊面積323.08hm2。造林過程均采用大苗造林方式,喬木的平均胸徑8.13cm、樹高5.09m。造林地主要來源于農(nóng)耕地、騰退的建設(shè)用地、沙荒地等類型。

    圖1 北京市平原造林分布與熱場分布Fig.1 The spatial distribution of the plain forestation patches and land surface temperature

    1.3 地面亮溫反演

    地面亮溫?cái)?shù)據(jù)以Landsat- 8的熱紅外波段數(shù)據(jù)反演得到,在Landsat- 8中,有兩個(gè)熱紅外波段:波段10和波段11,有關(guān)研究研究表明[24],波段10反演的的地面亮溫能更好地突出高溫中心和低溫區(qū)域,有利于城市熱環(huán)境的分析,因此采用波段10來反演相關(guān)的地面亮溫。

    首先利用下式將波段10數(shù)據(jù)的DN值轉(zhuǎn)換為大氣頂部輻射亮度Lλ:

    Lλ=ML×Qcal+AL

    然后利用下式將輻射亮度轉(zhuǎn)化為地表亮度溫度Tλ:

    Tλ=K2/ ln(K1/Lλ+1)- 273. 15

    式中,Lλ為波段10的大氣頂部輻射亮度(W m-2sr-1um-1);ML為波段乘積調(diào)整因子,從遙感影像頭文件獲?。籄L為附加因子,也從頭文件獲?。籕cal為像元亮度值DN;Tλ為星上亮度溫度(℃);K1、K2為常數(shù),從頭文件獲取。

    通過上述步驟反演的2014年9月4日北京市地面亮溫分布狀況見圖1。

    1.4 景觀斑塊分級(jí)

    景觀斑塊的大小不同,其生態(tài)學(xué)意義差別很大。大量的國內(nèi)外研究結(jié)果表明,植被與水體斑塊的大小對于城市熱島效應(yīng)的減弱或冷島效應(yīng)的發(fā)揮具有重要意義[10,25-26]。目前對于斑塊大小的劃分不同學(xué)者之間因?yàn)檠芯繀^(qū)域、研究對象、研究目標(biāo)的不同而差異很大??紤]到平原大造林的目標(biāo)不同于傳統(tǒng)的城市園林綠化,其營造與后續(xù)經(jīng)營均是按照森林生態(tài)系統(tǒng)的經(jīng)營目標(biāo)與規(guī)范進(jìn)行的,故這里的斑塊大小規(guī)模的劃分,參照郭晉平在研究山西省關(guān)帝山的森林群落時(shí)所提出的標(biāo)準(zhǔn)[27](表1)。

    表1 林地斑塊規(guī)模劃分等級(jí)標(biāo)準(zhǔn)

    1.5 研究尺度的選擇

    北京市域面積16410km2,其中平原面積6338km2,占全市面積的38.6%,主要集中分布于市域范圍的東南部(圖1),為了突出平原區(qū)熱場的時(shí)空分異特征,同時(shí)也為了更好地研究林地降溫效應(yīng)的空間分布特點(diǎn),根據(jù)平原區(qū)地貌分異差異、人類活動(dòng)強(qiáng)度以及城市建設(shè)用地的空間擴(kuò)展特點(diǎn)等,將整個(gè)平原區(qū)劃分為延慶盆地、六環(huán)以內(nèi)、六環(huán)以北平原、六環(huán)以南平原(六環(huán)以北和六環(huán)以南以通州和順義的行政邊界為限)等4個(gè)區(qū)域單元。

    2 結(jié)果與討論

    2.1 總體亮溫變化

    林地斑塊的冷島效應(yīng)空間包括了兩個(gè)部分:林地斑塊本身所占據(jù)的地表區(qū)域以及緊靠林地斑塊外圍一定距離范圍內(nèi)的非綠色植被空間范圍。文獻(xiàn)資料表明,綠地降溫效應(yīng)的最大外圍邊界距離在200—500m之間,超過500m之后就沒有直接關(guān)聯(lián)了[14],基于此,利用GIS的緩沖區(qū)分析功能,以現(xiàn)有的百萬畝平原造林地斑塊為主體,在其外圍0—500m范圍內(nèi),以50m為基礎(chǔ)做緩沖區(qū),通過比較不同緩沖區(qū)林地斑塊的平均溫度,一方面顯示林地斑塊的冷島效應(yīng)強(qiáng)度大小,另一方面也可以反映林地斑塊的降溫效應(yīng)隨林地斑塊外圍距離逐漸變化的變化特征,相關(guān)的統(tǒng)計(jì)結(jié)果見圖2。

    從圖2可以看出,無論是全部還是不同等級(jí)的造林地斑塊,在其邊界之外,隨著與林地斑塊距離的逐漸增大,林地斑塊的降溫效應(yīng)呈現(xiàn)了逐步衰減的變化過程。以全部斑塊的平均狀況而論,至350m以外的緩沖區(qū)范圍內(nèi),其溫度變化幾近飽和,據(jù)此可以認(rèn)為,林地斑塊降溫效應(yīng)的最大邊界距離在其邊界之外350m左右,若以350—400m緩沖區(qū)距離內(nèi)的平均溫度28.237℃作為不受林地影響的本底背景溫度,將其與林地斑塊本身的平均溫度27.213℃相比,則全部林地斑塊本身的降溫幅度可達(dá)1.023℃,如果以2004年9月8日作為造林前的林地斑塊范圍內(nèi)的平均溫度26.922℃做背景參考,則新造林林地斑塊的降溫效應(yīng)可以達(dá)到1.954℃;在林地斑塊外圍的有效降溫距離內(nèi),以最靠近林地斑塊的0—100m的緩沖區(qū)距離范圍內(nèi)的降溫效果最大,達(dá)到0.392—0.577℃,隨著距離林地邊緣距離的逐步遞增,其降溫效果在逐漸減弱,在林地斑塊外150—200m的范圍內(nèi)降溫效果尚可達(dá)到0.123℃,而到了300—350m的距離范圍內(nèi),降溫幅度僅有0.014℃。

    另外從圖2還可以看出,不同級(jí)別林地斑塊其降溫效應(yīng)也存在很大差異。首先斑塊規(guī)模越大,相同外圍距離處的溫度就越低,以林地斑塊外0—50m距離緩沖區(qū)為例,小斑塊、中斑塊、大中斑塊、大斑塊、超大斑塊和巨斑塊的溫度分別為27.81℃、27.48℃、27.51℃、27.59℃、26.94℃和27.21℃,巨型斑塊的影響溫度比小斑塊整整低了0.6℃,其他緩沖區(qū)也有相同的變化趨勢;其次,斑塊規(guī)模越大,其溫度影響的距離效應(yīng)愈大,從圖3可以看出,小斑塊和中斑塊對其外圍溫度影響的最大距離在林地斑塊外圍350—400m間,大中斑塊、大斑塊、超大斑塊和巨斑塊的外圍影響距離都在在400—450m之間。

    圖2 平原造林地斑塊降溫效應(yīng)的緩沖區(qū)分析Fig.2 The buffer analysis of cooling effect on forest patches outside

    由于林地斑塊外500m緩沖區(qū)范圍內(nèi)的土地覆蓋情況差異較大,這種環(huán)境背景會(huì)對研究結(jié)果造成一定的影響。2010年土地利用的分析結(jié)果表明(表2),500m緩沖區(qū)范圍內(nèi)的土地利用類型以耕地(占69.2%)和農(nóng)村居民點(diǎn)建設(shè)用地(占24.5%)為主,其他地類的占比都很小,因此這兩種土地利用類型對研究結(jié)果影響較大,由于由于百萬畝平原造林地大多是在耕地、未利用地、水域邊緣、草地等的基礎(chǔ)上建成的,因此,這些地類對研究結(jié)果的影響可以作為背景忽略不予考慮,但農(nóng)村居民點(diǎn)建設(shè)用地的影響較大,由于居民點(diǎn)建設(shè)用地的熱效應(yīng)為正向效應(yīng),其對本文研究的冷島效應(yīng)起消減作用,目前從技術(shù)上很難將這種影響完全去除。但由于農(nóng)村居民點(diǎn)一般單個(gè)斑塊面積不大、且空間分布比較零散,因此目前的研究結(jié)果可以作為冷島效應(yīng)的低限來看待。

    表2 林地斑塊外500m緩沖區(qū)內(nèi)土地利用

    2.2 不同區(qū)域造林地斑塊的亮溫差異

    圖3 北京市平原造林不同區(qū)域的林地斑塊亮溫比較 Fig.3 The land surface temperature comparison among different plain region

    根據(jù)前面的平原分區(qū)方案,對各區(qū)域內(nèi)林地斑塊2014年9月4日的亮溫所做統(tǒng)計(jì)結(jié)果見圖3。從圖3可以看出,不同區(qū)域的平原造林地塊,其溫度差異明顯,其中以延慶盆地的造林地斑塊的溫度最低,平均為24.718℃,比全市造林地斑塊的平均溫度整整低了2.5℃,冷島效應(yīng)幅度可以達(dá)到3.519℃,這可能與延慶盆地地處燕山山脈包圍之中,且其海拔平均較高有關(guān)。而在最大的北京平原區(qū),以6環(huán)以內(nèi)的新造林地塊的溫度最高,達(dá)到了28.126℃,比全市造林地的平均溫度還高出了0.913℃,受強(qiáng)烈的城市化過程的影響,其冷島效應(yīng)強(qiáng)度只有0.111℃;而六環(huán)以南和六環(huán)以北兩個(gè)區(qū)域的的造林地斑塊平均溫度相差不大,但均比全部造林地斑塊的平均溫度略高,其冷島效應(yīng)強(qiáng)度分別為0.967℃和0.883℃。而從亮溫溫度極差來看,其區(qū)域差異要比平均亮溫的變化明顯許多,總體呈現(xiàn)了沿延慶盆地—6環(huán)以北-六環(huán)以南-六環(huán)以內(nèi)這一梯度逐步擴(kuò)大的規(guī)律。

    2.3 不同造林年份的森林景觀斑塊的亮溫差異

    北京市的平原造林工程開始實(shí)施于2012年,截止2014年,共完成了平原造林任務(wù)63333hm2,由于不同造林地塊的造林年份不同,這勢必會(huì)在造林林地的降溫效應(yīng)上有所反映,為量化這種差異,也對不同造林年份的林地斑塊的平均亮溫進(jìn)行了統(tǒng)計(jì)(圖4)。

    圖4 不同造林年份林地斑塊的平均亮溫Fig.4 The mean land surface temperature of the forest patches in different year

    從圖4可以看出,不同年份的造林地斑塊的亮溫溫度差異還是客觀存在的,但其絕對差異的幅度較小,介于0.05—0.13℃之間,其中2012和2013年造林地斑塊的平均亮溫要大于區(qū)域林地整體的平均亮溫。尤其值得注意的是,統(tǒng)計(jì)結(jié)果顯示,越是造林晚的林地斑塊,其降溫效果似乎越明顯,這與一般的認(rèn)識(shí)有所差異。一般而言,造林地的時(shí)間愈長,其系統(tǒng)的穩(wěn)定性相對而言愈高、植物的生長發(fā)育狀況愈好,因此其降溫效果愈明顯。之所以會(huì)出現(xiàn)這種反常變化,主要有兩方面的原因。首先是林地斑塊的景觀格局的差異。國內(nèi)外的相關(guān)研究結(jié)果表明,景觀格局對于綠色植被的降溫效果有很大的影響,其中,綠色植被斑塊的周長-面積比率與其表面亮溫成正比,而斑塊的總面積和平均面積與其表面亮溫成反比[23, 28]。從相關(guān)年份林地斑塊的相關(guān)格局指數(shù)來看(表3),按照2012年到2014年的時(shí)間序列,其周長-面積比率呈現(xiàn)出了逐步縮小的變化過程,而平均面積與林地總面積則都呈現(xiàn)逐步增大的變化趨勢,這一切都預(yù)示著其林地斑塊的亮溫會(huì)逐年降低,這一實(shí)際變化結(jié)果既進(jìn)一步印證了相關(guān)文獻(xiàn)中的結(jié)論,也從景觀格局方面說明了上述不同年限造林地塊溫度隨時(shí)間序列逐步降低的內(nèi)在原因。

    除了景觀格局的原因之外,上述情況的出現(xiàn)可能還與工程造林類型中的濕地保護(hù)與建設(shè)有一定的關(guān)系。在平原造林工程的實(shí)施過程中,共包括了景觀生態(tài)林、綠色通道和濕地保護(hù)與建設(shè)等三大類型(表4),以三大類型工程的總體情況看,平原造林以景觀生態(tài)林建設(shè)和綠色通道建設(shè)為核心,兩者合計(jì)的工程量占到了工程總量的98%以上,從年度推進(jìn)情況來看,這兩類工程的推進(jìn)情況雖然2012年與其它兩個(gè)年份有一定的差異,但2013和2014兩年的推進(jìn)幅度與比例相差不大,年度推進(jìn)中幅度變化最顯著的是濕地保護(hù)與建設(shè)工程,2012年實(shí)施的該工程面積僅占3年來該類工程實(shí)施總量的7.94%,2013年的占22.88%,2014年的實(shí)施面積占到了工程總實(shí)施量的69.17%。Chen等對廣州公園綠地降溫效應(yīng)的研究中發(fā)現(xiàn)[29],當(dāng)公園綠地中的水體面積大于12.89hm2時(shí),公園的降溫效果會(huì)更加明顯,也就是說林水的有效結(jié)合可以增強(qiáng)綠地的降溫效果,從不同工程類別土地的地表亮溫情況看,濕地保護(hù)與建設(shè)工程造林地塊的平均亮溫為26.93℃,分別比綠色通道和景觀生態(tài)林工程造林區(qū)的平均亮溫低了0.41℃和0.31℃。2012—2014年平原造林地塊的溫度變化的年際差異也可能與此有很大的關(guān)系。

    表3 平原造林地斑塊的景觀格局指數(shù)

    表4 平原造林地不同工程類型統(tǒng)計(jì)

    *MLST:Mean land surface temperature

    2.4 造林地斑塊尺度大小與降溫效應(yīng)

    2.4.1 斑塊大小的構(gòu)成分析

    根據(jù)景觀生態(tài)學(xué)理論,景觀斑塊的大小不同,其內(nèi)部的包含的物質(zhì)與能量有差異,因此會(huì)影響到景觀斑塊的一些表觀功能特征。根據(jù)郭晉平的林地斑塊劃分標(biāo)準(zhǔn)對3年來北京平原造林地塊的斑塊尺度所做的統(tǒng)計(jì)結(jié)果見圖5。

    從圖5可以看出,在林地景觀斑塊數(shù)量上,以小斑塊占絕對優(yōu)勢,其數(shù)量比例占到了全部斑塊數(shù)量的83.64%,其次為中斑塊類型,但其數(shù)量比例只有11.82%,其他的斑塊類型的數(shù)量比例均在3%以下。

    從林地景觀斑塊的分級(jí)面積來看,與斑塊數(shù)量變化具有相同的變化趨勢,但變化的劇烈程度要有很大緩和。大致可以分為3個(gè)量級(jí):中、小斑塊為第1級(jí),其所占的面積比例都在25%以上;中大斑塊與大斑塊為第2級(jí),所占面積比例在10%—15%之間;超大斑塊與巨斑塊為第3級(jí),面積比例均在10%以下。

    圖5 平原造林地斑塊大小的尺度分析Fig.5 The scale analysis of the forest patches

    2.4.2 不同等級(jí)斑塊的亮溫分析

    不同規(guī)模大小的林地斑塊的亮溫統(tǒng)計(jì)結(jié)果見圖6。

    圖6 不同等級(jí)林地斑塊的亮溫分布Fig.6 The mean land surface temperature in different forest patch classes

    從圖6可以看出,斑塊大小對亮溫的影響總體來說是斑塊面積尺度越大,降溫效應(yīng)愈明顯。例如,面積小于10hm2規(guī)模的小斑塊,其平均亮溫為27.48℃,比林地斑塊的平均溫度27.21℃高出了0.27℃,中斑塊、大中斑塊和大斑塊類型的亮溫情況與小斑塊類似,其平均溫度也都高于全部林地斑塊的平均亮溫,只有超大斑塊和巨斑塊的平均亮溫分別比全部林地斑塊的平均亮溫分別低了0.43℃和0.42℃,如果以全部林地斑塊外圍350—400m緩沖區(qū)距離內(nèi)的平均溫度28.237℃作為不受林地影響的參考背景溫度,則其冷島效應(yīng)強(qiáng)度分別達(dá)到了1.457℃和1.449℃。

    2.5 平原造林工程降溫效應(yīng)的價(jià)值評價(jià)

    對于植被降溫效應(yīng)的價(jià)值評估,國內(nèi)外都做了大量的工作,其基本流程是:首先計(jì)算植被蒸騰所吸收的熱量,之后再在溫度降低的能量被全部用于植被蒸騰作用的假設(shè)前提下,將溫度降低的數(shù)值轉(zhuǎn)換為植被蒸騰消耗的熱量值,之后通過電能節(jié)約環(huán)節(jié),再將熱量值轉(zhuǎn)換成電能;最后通過居民用電電價(jià)就可以將夏季林地的降溫功能轉(zhuǎn)化成以貨幣量化的生態(tài)價(jià)值。

    參照相關(guān)案例的研究方法及其在北京計(jì)算時(shí)的有關(guān)參數(shù)[30-32],按照每年90d的高溫期計(jì)算,計(jì)算結(jié)果見表5。

    從表5可以看出,平原造林地塊除了其自身68050.98hm2的降溫面積之外,通過冷島效應(yīng)向周邊的輻射作用,形成的降溫面積總計(jì)可達(dá)250212.88hm2,其中降溫輻射較強(qiáng)的0—100m邊界外圍范圍即達(dá)到了86544.66hm2,已經(jīng)高于造林工程本身所覆蓋的地表面積。從其蒸騰降溫所消耗的熱能來看,每年林地本身降溫節(jié)能209912.48×108J;外圍間接降溫節(jié)能總計(jì)180833×108J,其中0—100m外圍邊界范圍內(nèi)降溫消耗的熱量占到了70.97%。

    按照居民用電價(jià)格0.5元/kWH計(jì)算,平原大造林形成的降溫效應(yīng)總價(jià)值為4.8882×108元,其中林地本身的降溫價(jià)值達(dá)到了2.626×108元,占總價(jià)值的53.72%,通過本身冷島向周邊輻射引起的間接降溫效應(yīng)的價(jià)值為2.2622×108元,這其中緊靠林地斑塊外圍0—100m范圍內(nèi)的間接降溫價(jià)值達(dá)到了1.6055×108元。

    表5 北京市平原造林降溫功能評估

    3 結(jié)論

    (1) 從2012—2014年的3a間,北京市的平原大造林工作累計(jì)造林68050.98hm2,形成新的林地景觀斑塊10648個(gè)。從其斑塊構(gòu)成看,在林地景觀斑塊數(shù)量上,以小斑塊占絕對優(yōu)勢,其數(shù)量比例占到了全部斑塊數(shù)量的83.64%,其次為中斑塊類型,但其數(shù)量比例只有11.82%,其他的斑塊類型的數(shù)量比例均在3%以下;從林地景觀斑塊的分級(jí)面積看,與斑塊數(shù)量變化具有相同的變化趨勢,但變化的劇烈程度要有很大緩和,其中中、小斑塊所占的面積比例都在25%以上;中大斑塊與大斑塊所占面積比例在10%—15%之間,而超大斑塊與巨斑塊的面積比例均在10%以下。

    (2)利用2014年9月4日landsat- 8所反演的地表亮溫看,北京市平原區(qū)百萬畝大造林工程,其區(qū)域降溫效果與效益都是非常顯著的。平均而言,林地本身的降溫效果可以達(dá)到1.023℃,林地的降溫輻射范圍可以到達(dá)林地外圍350m的距離。從斑塊尺度對亮溫降低的效應(yīng)看,不同級(jí)別林地斑塊其降溫效應(yīng)差異很大,其總體趨勢是,斑塊規(guī)模越大,相同外圍距離處的溫度就越低,同時(shí)斑塊規(guī)模越大,其溫度影響的距離效應(yīng)愈大,但其降溫效應(yīng)的影響范圍與已有文獻(xiàn)的結(jié)論類似,均未超出林地斑塊邊緣外500m的距離范圍。

    (3)平原造林地塊除了其自身68050.98hm2的降溫面積之外,通過冷島效應(yīng)向周邊的輻射作用,形成的降溫面積總計(jì)可達(dá)250212.88 hm2。初步的價(jià)值估算表明,平原大造林形成的降溫效應(yīng)總價(jià)值為4.8882×108元,其中林地本身的降溫價(jià)值占總價(jià)值的53.72%,通過本身冷島向周邊輻射引起的間接降溫效應(yīng)的價(jià)值為2.2622×108元,這其中緊靠林地斑塊外圍0—100m范圍內(nèi)的間接降溫價(jià)值達(dá)到了1.6055×108元。

    (4)有關(guān)研究的結(jié)果表明,當(dāng)公園面積大于12hm2之后,隨著公園面積的進(jìn)一步擴(kuò)大,公園的冷島效應(yīng)會(huì)越強(qiáng)[33],本次研究的結(jié)果也顯示,面積100hm2以上的超大林地斑塊和巨型林地斑塊的降溫效應(yīng)最為突出。這啟示在今后的平原人工造林與平原人工林的后續(xù)經(jīng)營中,在工程設(shè)計(jì)之初就應(yīng)該注意到林地斑塊大小的設(shè)計(jì)問題,盡量利用地形與地勢,建立大的林地斑塊,這樣一方面有助于增強(qiáng)所建林地的降溫效應(yīng),另外,面積設(shè)計(jì)大了之后,也更有利于后續(xù)的以郊野公園等形式為主體的林地生態(tài)效益的深度開發(fā)利用。

    [1] 胡嘉驄, 魏信, 陳聲海. 北京城市熱場時(shí)空分布及景觀生態(tài)因子研究. 北京: 北京師范大學(xué)出版社, 2014: 1- 6.

    [2] Li J X, Song C H,Gao L, Zhu F G, Meng X L, Wu J G. Impacts of landscape structure on surface urban heat islands: A case study of Shanghai, China. Remote Sensing of Environment, 2011, 115(12): 3249- 3263.

    [3] 黃良美, 鄧超冰, 黎寧. 城市熱島效應(yīng)熱點(diǎn)問題研究進(jìn)展.氣象與環(huán)境學(xué)報(bào), 2011, 27(4): 54- 58.

    [4] Rizwan A M, Dennis L Y C, Liu C H. A review on the generation, determination and mitigation of urban heat island. Journal of Environmental Sciences, 2008, 20(1): 120- 128.

    [5] Zhang Z, Lu Y M, Pan H T. Cooling and humidifying effect of plant communities in subtropical urban parks. Urban Forestry & Urban Greening, 2013, 12(3): 323- 329.

    [6] Chang C R, Li M H, Chang S D. A preliminary study on the local cool-island intensity of Taipei city parks. Landscape and Urban Planning, 2007, 80(4): 386- 395.

    [7] 周東穎, 張麗娟, 張利, 范懷欣, 劉棟. 城市景觀公園對城市熱島調(diào)控效應(yīng)分析——以哈爾濱市為例.地域研究與開發(fā), 2011, 30(3): 73- 78.

    [8] 王劍強(qiáng), 王志泰. 基于緩解熱島效應(yīng)的山地城市生態(tài)斑塊研究——以貴州省貞豐縣為例. 西北林學(xué)院學(xué)報(bào), 2014, 29(2): 232- 236.

    [9] 劉嬌妹, 李樹華, 楊志峰. 北京公園綠地夏季溫濕效應(yīng). 生態(tài)學(xué)雜志, 2008, 27(11): 1972- 1978.

    [10] 陳朱, 陳方敏, 朱飛鴿, 曹璐, 沈興華, 李俊祥. 面積與植物群落結(jié)構(gòu)對城市公園氣溫的影響. 生態(tài)學(xué)雜志, 2011, 30(11): 2590- 2596.

    [11] 馮曉剛, 石輝. 基于遙感的夏季西安城市公園“冷效應(yīng)”研究. 生態(tài)學(xué)報(bào), 2012, 32(23): 7355- 7363.

    [12] Onishi A, Cao X, Ito T, Shi F, Imura H. Evaluating the potential for urban heat-island mitigation by greening parking lots. Urban Forestry & Urban Greening, 2010, 9(4): 323- 332.

    [13] Ren Z B, He X Y,Zheng H F, Zhang D, Yu X Y, Shen G Q, Guo R C. Estimation of the relationship between urban park characteristics and park cool island intensity by remote sensing data and field measurement. Forests, 2013, 4(4): 868- 886.

    [14] Hamada S, Ohta T. Seasonal variations in the cooling effect of urban green areas on surrounding urban area. Urban Forestry & Urban Greening, 2010, 9(1): 15- 24.

    [15] 高玉福, 李樹華, 朱春陽. 城市帶狀綠地林型與溫濕效益的關(guān)系. 中國園林, 2012, 28(1): 94- 97.

    [16] 閆偉姣, 孔繁花, 尹海偉, 孫常峰, 許峰,李文超,張嘯天. 紫金山森林公園降溫效應(yīng)影響因素. 生態(tài)學(xué)報(bào), 2014, 34(12): 3169- 3178.

    [17] 馬雪梅, 張友靜, 黃浩. 城市熱場與綠地景觀相關(guān)性定量分析. 國土資源遙感, 2005, 17(3): 10- 13.

    [18] Li X M, Zhou W Q, Ouyang Z Y, Zheng H. Spatial pattern of greenspace affects land surface temperature: evidence from the heavily urbanized Beijing metropolitan area, China. Landscape Ecology, 2012, 27(6): 887- 898.

    [19] 馮悅怡, 胡潭高,張力小. 城市公園景觀空間結(jié)構(gòu)對其熱環(huán)境效應(yīng)的影響. 生態(tài)學(xué)報(bào), 2014, 34(12): 3179- 3187.

    [20] 周雅星, 劉茂松, 徐馳, 方芳, 鐘晶晶, 張明娟. 南京市市域熱場分布與景觀格局的關(guān)聯(lián)分析. 生態(tài)學(xué)雜志, 2014, 33(8): 2199- 2206.

    [21] 李東海, 艾彬, 黎夏. 基于遙感和GIS的城市水體緩解熱島效應(yīng)的研究——以東莞市為例. 熱帶地理, 2008, 28(5): 414- 418.

    [22] 岳文澤, 徐麗華. 城市典型水域景觀的熱環(huán)境效應(yīng). 生態(tài)學(xué)報(bào), 2013, 33(6): 1852- 1859.

    [23] Chen A L, Yao X A, Sun R H, Chen L D. Effect of urban green patterns on surface urban cool islands and its seasonal variations. Urban Forestry & Urban Greening, 2014, 13(4): 646- 654.

    [24] 王鵬龍, 張建明, 呂榮芳. 基于空間自相關(guān)的蘭州市熱環(huán)境. 生態(tài)學(xué)雜志, 2014, 33(4): 1089- 1095.

    [25] 葛偉強(qiáng), 周紅妹, 楊引明, 丁金才. 基于遙感和GIS的城市綠地緩解熱島效應(yīng)作用研究. 遙感技術(shù)與應(yīng)用, 2006, 21(5): 432- 435.

    [26] Kong F H, Yin H W, James P, Hutyra L R, He H S. Effects of spatial pattern of greenspace on urban cooling in a large metropolitan area of eastern China. Landscape and Urban planning, 2014, 128: 35- 47.

    [27] 郭晉平. 森林景觀生態(tài)研究. 北京: 北京大學(xué)出版社, 2001:170- 178.

    [28] 孟丹, 李小娟, 宮輝力, 趙文吉. 北京地區(qū)熱力景觀格局及典型城市景觀的熱環(huán)境效應(yīng). 生態(tài)學(xué)報(bào), 2010, 30(13): 3491- 3500.

    [29] Chen X Z, Su Y X, Li D, Huang G Q, Chen W Q, Chen S S. Study on the cooling effects of urban parks on surrounding environments using Landsat TM data: a case study in Guangzhou, southern China. International Journal of Remote Sensing, 2012, 33(18): 5889- 5914.

    [30] 楊士弘. 城市綠化樹木的降溫增濕效應(yīng)研究. 地理研究, 1994, 13(4): 74- 80.

    [31] Lin W Q, Wu T H, Zhang C G, Yu T. Carbon savings resulting from the cooling effect of green areas: A case study in Beijing. Environmental Pollution, 2011, 159(8/9): 2148- 2154.

    [32] Zhang B,Xie G D, Gao J X, Yang Y. The cooling effect of urban green spaces as a contribution to energy-saving and emission-reduction: A case study in Beijing, China. Building and Environment, 2014, 76: 37- 43.

    [33] Lu J, Li C D, Yang Y C, Zhang X H, Jin M. Quantitative evaluation of urban park cool island factors in mountain city. Journal of Central South University of Technology, 2012, 19(6): 1657- 1662.

    The cooling effect of plain afforestation in the Beijing Project and its remote sensing-based valuation

    JIA Baoquan1,2,3,*, QIU Kuanbiao1,4

    1ResearchInstituteofForestry,ChineseAcademyofForestry,Beijing100091,China2KeyLaboratoryofTreeBreedingandCultivation,StateForestryAdministration,Beijing100091,China3ResearchCentreofUrbanForestry,StateForestryAdministration,Beijing100091,China4KeyLaboratoryofSoilandWaterConservationandDesertificationCombating,MinistryofEducation,CollegeofSoilandWaterConservation,BeijingForestryUniversity,Beijing100083,China

    With the rapid development of urbanization, the urban heat island effect has exerted a great influence on the lives and productivity of urban residents. Many research has suggested that the cooling effect caused by vegetation and water bodies is an economic way to efficiently and substantially mitigate the urban heat island effects at a local scale. However, because the previous studies mainly concentrated on a limited number of urban parks, the past result might have a large uncertainty. Since plain afforestationwas implemented in the Beijing Project from 2012 to 2014, 10,648 forest patches have been constructed, therefore providing a great sample for examining the cooling effects of artificial forests.We analyzed the cooling effects of these artificial forests using the land surface brightness temperature as an indicator of urban heat islands. For this, we used the LANDSAT- 8 images of 4 September 2014 and additional information about the originally afforested patches. The artificial forests in this project decreased inbrightness temperature by 1.023℃,despitethe short afforestation period. The latter can be explained by the large young plants that were used for the afforestation, as well as the special implementation mode used. A buffer zone analysis showed that the influence of this cooling effect reached as far as 350m from the edge of the artificial forests. The strongest cooling effect was found within 100m from the forest edges, with a decrease in brightness temperature of 0.392℃.In addition, a zonal analysis showed that the artificial forests had the most significant cooling effects in the Yanqing Basin, with a cooling of 3.519℃. The smallest significant decrease, 0.111℃, was found within the 6th Ring Road. An analysis between different years showed that the landscape pattern of artificial forest patches, as well as wetland protection and construction,were important for establishing the cooling effects. According to a patch size frequency analysis, the small patches had the highest frequency, whereas the middle and small patches covered a large part of the total area. Nevertheless, the large artificial forest patches had a larger significant cooling effect compared with the small ones. According to the primary valuation,the cooling effect of this project had a value as high as 488 millionCNY. Approximately 54% of this total value was caused by cooling within the forest patches, whereas cooling within the influential range from the forest patch edges caused approximately 46%.For further afforestation intended for effective mitigation of the urban heat island effects, several factors should receive careful attention. First, large artificial forest patches are preferablein terms of forest patch design. Second, for the configuration of the forest patches,two adjacent forest patches should be between 100 and 500m apart, in order to benefit from the cooling effects of the forest patches′influential range.

    afforestation project; cooling effect; LANDSAT- 8; valuation; plain; Beijing

    林業(yè)公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201404301);北京市自然科學(xué)基金資助項(xiàng)目(8152031)

    2015- 08- 22;

    日期:2016- 06- 14

    10.5846/stxb201508221755

    *通訊作者Corresponding author.E-mail: jiabaoquan2006@163.com

    賈寶全,仇寬彪.北京市平原百萬畝大造林工程降溫效應(yīng)及其價(jià)值的遙感分析.生態(tài)學(xué)報(bào),2017,37(3):726- 735.

    Jia B Q, Qiu K B.The cooling effect of plain afforestation in the Beijing Project and its remote sensing-based valuation.Acta Ecologica Sinica,2017,37(3):726- 735.

    猜你喜歡
    冷島亮溫平原
    那一片平原
    黃河之聲(2022年6期)2022-08-26 06:48:50
    極熱條件下空冷島永磁電動(dòng)機(jī)散熱能力研究
    高寒條件下空冷島風(fēng)機(jī)集群永磁直驅(qū)系統(tǒng)批量改造分析
    高寒條件下空冷島變頻永磁直驅(qū)改造及應(yīng)用
    霰譜分布特征對強(qiáng)對流云高頻微波亮溫影響的模擬研究
    平原的草
    基于南太平洋的AMSR2 L1R亮溫?cái)?shù)據(jù)質(zhì)量評估
    浪起山走
    滇池(2019年1期)2019-02-14 02:36:54
    空冷島的熱態(tài)沖洗及防凍措施
    2014年2月12日新疆于田MS7.3地震熱紅外亮溫異常分析
    中國地震(2015年1期)2015-11-08 11:11:32
    国产久久久一区二区三区| 韩国av在线不卡| 91久久精品国产一区二区成人| 国产精品一区二区三区四区久久| 国产主播在线观看一区二区| 国产真实伦视频高清在线观看 | 日韩欧美精品v在线| 人人妻人人澡欧美一区二区| 国产精品综合久久久久久久免费| 国产三级中文精品| 美女大奶头视频| av中文乱码字幕在线| 欧美一区二区精品小视频在线| 亚洲精品亚洲一区二区| 人人妻,人人澡人人爽秒播| 久久精品91蜜桃| 精品人妻1区二区| 精品一区二区免费观看| www.色视频.com| 国产一区二区三区在线臀色熟女| videossex国产| 亚洲精品久久国产高清桃花| 99九九线精品视频在线观看视频| 欧美一区二区国产精品久久精品| 嫩草影院精品99| а√天堂www在线а√下载| videossex国产| 亚洲欧美日韩高清专用| 三级毛片av免费| 欧美丝袜亚洲另类 | 亚洲精品日韩av片在线观看| 精品人妻视频免费看| 色噜噜av男人的天堂激情| 亚洲av美国av| 淫秽高清视频在线观看| 日韩av在线大香蕉| 中文字幕精品亚洲无线码一区| 高清在线国产一区| 日韩一区二区视频免费看| 国产一区二区三区视频了| 久久精品国产亚洲av香蕉五月| 亚洲av熟女| 国产v大片淫在线免费观看| 特级一级黄色大片| avwww免费| 99国产精品一区二区蜜桃av| 精品人妻一区二区三区麻豆 | av在线蜜桃| 十八禁网站免费在线| 亚洲四区av| 观看免费一级毛片| 老熟妇仑乱视频hdxx| 波多野结衣高清作品| 久久国内精品自在自线图片| 成人欧美大片| 村上凉子中文字幕在线| 欧美日韩乱码在线| 久久久国产成人精品二区| 国语自产精品视频在线第100页| 亚洲欧美清纯卡通| 亚洲七黄色美女视频| 久久久精品欧美日韩精品| 国产国拍精品亚洲av在线观看| 乱人视频在线观看| 天天躁日日操中文字幕| 国产精华一区二区三区| 在线看三级毛片| 极品教师在线视频| 国产色爽女视频免费观看| 欧美日本视频| 国产黄片美女视频| 大型黄色视频在线免费观看| 一级a爱片免费观看的视频| 97超视频在线观看视频| 亚洲av一区综合| 亚洲欧美清纯卡通| 免费看av在线观看网站| 91麻豆av在线| 亚洲图色成人| 国产精品女同一区二区软件 | 国产成人a区在线观看| 麻豆一二三区av精品| 亚洲av五月六月丁香网| 精品人妻偷拍中文字幕| 亚洲国产精品sss在线观看| 我的女老师完整版在线观看| 噜噜噜噜噜久久久久久91| 亚洲七黄色美女视频| 久久久色成人| 亚洲国产精品久久男人天堂| 一区二区三区激情视频| 免费不卡的大黄色大毛片视频在线观看 | 丰满乱子伦码专区| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看| 黄色女人牲交| 成人av在线播放网站| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 久久6这里有精品| 日本三级黄在线观看| 一级黄色大片毛片| 国产一区二区激情短视频| 99热6这里只有精品| 亚洲无线观看免费| 日韩欧美一区二区三区在线观看| 精品一区二区三区视频在线| av在线亚洲专区| 少妇被粗大猛烈的视频| 性色avwww在线观看| 中文字幕久久专区| 亚洲专区中文字幕在线| 校园人妻丝袜中文字幕| 色尼玛亚洲综合影院| 又紧又爽又黄一区二区| 成年女人永久免费观看视频| 干丝袜人妻中文字幕| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜添av毛片 | 免费在线观看成人毛片| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 久久午夜福利片| 蜜桃亚洲精品一区二区三区| 人人妻人人看人人澡| 精品久久久久久久人妻蜜臀av| 日韩av在线大香蕉| 色精品久久人妻99蜜桃| 国产日本99.免费观看| 十八禁国产超污无遮挡网站| 18+在线观看网站| 身体一侧抽搐| 99久久精品热视频| 精品免费久久久久久久清纯| 国产成人aa在线观看| 天堂网av新在线| 久久久久久九九精品二区国产| 淫秽高清视频在线观看| 精品99又大又爽又粗少妇毛片 | 欧美zozozo另类| 乱码一卡2卡4卡精品| 一进一出抽搐gif免费好疼| 亚洲最大成人av| 亚洲av不卡在线观看| 国产高清不卡午夜福利| 久久精品综合一区二区三区| 精品久久久久久成人av| 亚洲 国产 在线| 最近中文字幕高清免费大全6 | 日本色播在线视频| 亚洲精品久久国产高清桃花| 在线观看舔阴道视频| 美女免费视频网站| 一区二区三区免费毛片| 美女被艹到高潮喷水动态| 亚洲国产精品sss在线观看| 中国美女看黄片| 国产欧美日韩精品一区二区| 国产精品无大码| 亚洲精品影视一区二区三区av| 亚洲三级黄色毛片| 日本 av在线| 男人和女人高潮做爰伦理| 午夜日韩欧美国产| 久久亚洲精品不卡| 亚洲成人久久爱视频| 日韩国内少妇激情av| 免费电影在线观看免费观看| 亚洲欧美激情综合另类| 成人特级黄色片久久久久久久| 久久久久久久午夜电影| 国内久久婷婷六月综合欲色啪| 一区二区三区高清视频在线| videossex国产| 天堂av国产一区二区熟女人妻| 91av网一区二区| 国产成人影院久久av| 赤兔流量卡办理| 精品久久久久久久久亚洲 | 午夜福利在线在线| 国产亚洲精品久久久com| 看黄色毛片网站| 亚洲美女视频黄频| 欧美黑人欧美精品刺激| 色av中文字幕| 窝窝影院91人妻| 亚洲国产高清在线一区二区三| 亚洲一级一片aⅴ在线观看| 99久久无色码亚洲精品果冻| 免费在线观看成人毛片| 午夜福利成人在线免费观看| 老司机深夜福利视频在线观看| 日本黄大片高清| 男女做爰动态图高潮gif福利片| 99九九线精品视频在线观看视频| x7x7x7水蜜桃| 国产精品人妻久久久久久| 欧美+亚洲+日韩+国产| 国产精品国产三级国产av玫瑰| 国产精品久久视频播放| 女的被弄到高潮叫床怎么办 | 99在线人妻在线中文字幕| 最后的刺客免费高清国语| 亚洲欧美精品综合久久99| 人人妻人人看人人澡| 欧美区成人在线视频| 久久精品国产亚洲av香蕉五月| 天天躁日日操中文字幕| 九色国产91popny在线| 国产三级中文精品| 一边摸一边抽搐一进一小说| 国产探花极品一区二区| 国产一区二区激情短视频| 国产精品久久电影中文字幕| ponron亚洲| 九色成人免费人妻av| 一级黄片播放器| 国产成人一区二区在线| 国产极品精品免费视频能看的| 他把我摸到了高潮在线观看| 一个人观看的视频www高清免费观看| 女的被弄到高潮叫床怎么办 | 日本爱情动作片www.在线观看 | 男人狂女人下面高潮的视频| 免费看光身美女| 精品久久久久久久久亚洲 | 日本熟妇午夜| 少妇被粗大猛烈的视频| 亚洲av第一区精品v没综合| 一级av片app| 99久久九九国产精品国产免费| 一级a爱片免费观看的视频| 日本欧美国产在线视频| 黄色日韩在线| 美女xxoo啪啪120秒动态图| 又黄又爽又刺激的免费视频.| 国产真实伦视频高清在线观看 | 日本与韩国留学比较| 国产精华一区二区三区| .国产精品久久| 国产一区二区三区在线臀色熟女| 欧美最黄视频在线播放免费| 日本色播在线视频| 乱系列少妇在线播放| .国产精品久久| 俄罗斯特黄特色一大片| 日韩一区二区视频免费看| 久久精品国产亚洲av香蕉五月| 热99re8久久精品国产| 精品不卡国产一区二区三区| 岛国在线免费视频观看| 波多野结衣巨乳人妻| 一进一出好大好爽视频| 九九爱精品视频在线观看| 婷婷精品国产亚洲av在线| 亚洲欧美清纯卡通| 一级黄色大片毛片| 观看美女的网站| 免费在线观看日本一区| 丰满人妻一区二区三区视频av| 91精品国产九色| 国产乱人视频| 亚洲国产色片| 欧美三级亚洲精品| 最近最新免费中文字幕在线| 亚洲最大成人中文| 日本免费a在线| 成人欧美大片| 免费电影在线观看免费观看| 成熟少妇高潮喷水视频| 美女大奶头视频| 18禁黄网站禁片午夜丰满| 精品不卡国产一区二区三区| 毛片一级片免费看久久久久 | 国产黄a三级三级三级人| 亚洲不卡免费看| 国产男靠女视频免费网站| 国内精品宾馆在线| 日本黄色视频三级网站网址| videossex国产| 亚洲va日本ⅴa欧美va伊人久久| 久久精品综合一区二区三区| 大又大粗又爽又黄少妇毛片口| 国产伦精品一区二区三区视频9| 日韩 亚洲 欧美在线| 熟女电影av网| 成人无遮挡网站| 久久人妻av系列| 男女啪啪激烈高潮av片| 国产人妻一区二区三区在| 欧美日本视频| 丰满人妻一区二区三区视频av| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品青青久久久久久| 免费观看的影片在线观看| 免费无遮挡裸体视频| 免费观看人在逋| 真人做人爱边吃奶动态| 男女边吃奶边做爰视频| 国产精品野战在线观看| 少妇的逼好多水| 日韩欧美国产在线观看| 久久久久久久久大av| 欧美极品一区二区三区四区| 天堂√8在线中文| 中出人妻视频一区二区| 亚洲美女黄片视频| 一本久久中文字幕| 中国美白少妇内射xxxbb| 午夜影院日韩av| 久久精品综合一区二区三区| 免费看日本二区| 亚洲av二区三区四区| 国产精品一及| 99精品久久久久人妻精品| 小蜜桃在线观看免费完整版高清| 亚洲自拍偷在线| 日韩欧美精品免费久久| 啦啦啦韩国在线观看视频| 精品久久久久久,| 99久久精品国产国产毛片| 国产精品98久久久久久宅男小说| 国产高清有码在线观看视频| 欧美激情久久久久久爽电影| 一区二区三区四区激情视频 | 午夜精品一区二区三区免费看| 精品一区二区免费观看| 日本五十路高清| 成人性生交大片免费视频hd| 久久热精品热| 别揉我奶头~嗯~啊~动态视频| 日本-黄色视频高清免费观看| 亚洲精品久久国产高清桃花| 熟女电影av网| 久久国内精品自在自线图片| 国产一区二区三区视频了| www.www免费av| 一个人看的www免费观看视频| 69人妻影院| 国产伦一二天堂av在线观看| 超碰av人人做人人爽久久| 亚洲午夜理论影院| 亚洲性夜色夜夜综合| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 亚洲国产日韩欧美精品在线观看| 天天躁日日操中文字幕| 成人国产综合亚洲| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 国产av在哪里看| 俺也久久电影网| 日韩在线高清观看一区二区三区 | 久久精品国产亚洲av天美| 亚洲成人中文字幕在线播放| 成熟少妇高潮喷水视频| 热99在线观看视频| 久久久久久久久久久丰满 | 成人二区视频| 3wmmmm亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区免费观看| 国产色爽女视频免费观看| 丰满人妻一区二区三区视频av| 亚洲天堂国产精品一区在线| 免费不卡的大黄色大毛片视频在线观看 | 观看免费一级毛片| 国产精品一区二区三区四区久久| 可以在线观看毛片的网站| 国产国拍精品亚洲av在线观看| 人妻丰满熟妇av一区二区三区| 天堂动漫精品| 久久久精品大字幕| 国产成人福利小说| 亚洲精华国产精华液的使用体验 | 亚洲熟妇中文字幕五十中出| 美女被艹到高潮喷水动态| 国产成人福利小说| 日韩在线高清观看一区二区三区 | 99视频精品全部免费 在线| 免费大片18禁| 不卡视频在线观看欧美| 久久6这里有精品| 伦理电影大哥的女人| 亚洲无线在线观看| 欧美一级a爱片免费观看看| 国产高清视频在线播放一区| 黄色视频,在线免费观看| 91精品国产九色| 亚洲成人免费电影在线观看| 俄罗斯特黄特色一大片| 精品久久久久久久久久久久久| 国产精品综合久久久久久久免费| 亚洲五月天丁香| 看黄色毛片网站| 国产精品免费一区二区三区在线| 精品无人区乱码1区二区| 国产视频一区二区在线看| 嫩草影院新地址| 大又大粗又爽又黄少妇毛片口| 国产av不卡久久| 免费一级毛片在线播放高清视频| 舔av片在线| 熟妇人妻久久中文字幕3abv| 欧美黑人巨大hd| 变态另类成人亚洲欧美熟女| 午夜精品在线福利| 22中文网久久字幕| 97超级碰碰碰精品色视频在线观看| av在线观看视频网站免费| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 午夜久久久久精精品| 亚州av有码| 成人一区二区视频在线观看| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲av涩爱 | 91久久精品国产一区二区三区| 韩国av一区二区三区四区| 一夜夜www| 国产精品一区二区三区四区免费观看 | 禁无遮挡网站| 久久精品人妻少妇| 成人二区视频| 色综合亚洲欧美另类图片| 老女人水多毛片| 精品午夜福利在线看| 美女被艹到高潮喷水动态| 噜噜噜噜噜久久久久久91| 欧美人与善性xxx| 亚洲国产精品成人综合色| 九九在线视频观看精品| 熟女人妻精品中文字幕| 精品人妻视频免费看| 亚洲avbb在线观看| 国产成年人精品一区二区| 日日夜夜操网爽| 欧美性猛交╳xxx乱大交人| 欧美精品啪啪一区二区三区| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| 熟女人妻精品中文字幕| 成人三级黄色视频| 亚洲图色成人| 变态另类丝袜制服| 特级一级黄色大片| 最后的刺客免费高清国语| 国产免费一级a男人的天堂| 国产中年淑女户外野战色| 村上凉子中文字幕在线| 久久精品国产鲁丝片午夜精品 | 国产精品99久久久久久久久| 亚洲五月天丁香| 美女 人体艺术 gogo| 午夜免费成人在线视频| 99热这里只有精品一区| 欧美黑人欧美精品刺激| 91在线观看av| 丝袜美腿在线中文| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 亚洲国产精品sss在线观看| 久久精品国产鲁丝片午夜精品 | 特级一级黄色大片| 午夜老司机福利剧场| 国产爱豆传媒在线观看| 性欧美人与动物交配| 久久这里只有精品中国| 日日夜夜操网爽| 超碰av人人做人人爽久久| 婷婷精品国产亚洲av| 1024手机看黄色片| 国产真实伦视频高清在线观看 | 成年女人永久免费观看视频| 日韩,欧美,国产一区二区三区 | 嫩草影视91久久| 精品午夜福利视频在线观看一区| 高清日韩中文字幕在线| 小蜜桃在线观看免费完整版高清| 丝袜美腿在线中文| 两个人的视频大全免费| 免费看av在线观看网站| 欧美日韩精品成人综合77777| 国产在视频线在精品| 亚洲内射少妇av| 久久久久免费精品人妻一区二区| 嫩草影院入口| 日韩人妻高清精品专区| 亚洲avbb在线观看| 99在线人妻在线中文字幕| 国产精品99久久久久久久久| 桃红色精品国产亚洲av| 亚洲精品一区av在线观看| 亚洲图色成人| 变态另类丝袜制服| 亚洲中文日韩欧美视频| 99久久无色码亚洲精品果冻| 国产精品综合久久久久久久免费| 热99re8久久精品国产| 18禁黄网站禁片免费观看直播| 精品久久久久久久久久久久久| 男女边吃奶边做爰视频| 99久久九九国产精品国产免费| 日韩亚洲欧美综合| 色在线成人网| 乱码一卡2卡4卡精品| 日韩强制内射视频| 国产高清不卡午夜福利| 很黄的视频免费| 久久精品夜夜夜夜夜久久蜜豆| 免费观看精品视频网站| а√天堂www在线а√下载| 亚洲av五月六月丁香网| 午夜福利在线观看吧| av.在线天堂| 国产欧美日韩一区二区精品| 在线a可以看的网站| 国产一区二区三区在线臀色熟女| 午夜免费男女啪啪视频观看 | 可以在线观看毛片的网站| 久久中文看片网| 九色国产91popny在线| 国产精品国产三级国产av玫瑰| 国产大屁股一区二区在线视频| 国产精品久久久久久精品电影| 99热这里只有精品一区| 久久这里只有精品中国| 别揉我奶头 嗯啊视频| 一区福利在线观看| 欧美国产日韩亚洲一区| 他把我摸到了高潮在线观看| 一本精品99久久精品77| 99久久精品热视频| 午夜久久久久精精品| 亚洲av免费在线观看| 亚洲欧美日韩高清在线视频| 国产高清视频在线观看网站| 亚洲三级黄色毛片| 女生性感内裤真人,穿戴方法视频| 亚洲av中文字字幕乱码综合| 国产精品自产拍在线观看55亚洲| 伦精品一区二区三区| 丝袜美腿在线中文| av福利片在线观看| 男人的好看免费观看在线视频| 日本-黄色视频高清免费观看| 国内精品宾馆在线| 1000部很黄的大片| 三级毛片av免费| 午夜影院日韩av| 91麻豆av在线| 黄色女人牲交| 亚洲三级黄色毛片| 91久久精品国产一区二区成人| 精品人妻偷拍中文字幕| 国产精品永久免费网站| 露出奶头的视频| 欧美区成人在线视频| 成年免费大片在线观看| 亚洲18禁久久av| 又黄又爽又免费观看的视频| 国产成人aa在线观看| 久久久成人免费电影| 中文资源天堂在线| 日本欧美国产在线视频| 久久精品国产自在天天线| 亚洲av中文字字幕乱码综合| 一边摸一边抽搐一进一小说| 色在线成人网| 啦啦啦韩国在线观看视频| 欧美黑人巨大hd| 日本五十路高清| 精品一区二区三区av网在线观看| 69人妻影院| 22中文网久久字幕| 国产在视频线在精品| 观看美女的网站| 国产高清视频在线播放一区| 日本在线视频免费播放| 听说在线观看完整版免费高清| 国产熟女欧美一区二区| 欧美中文日本在线观看视频| 久久精品影院6| 国产黄片美女视频| 桃色一区二区三区在线观看| 十八禁网站免费在线| 精品福利观看| 无人区码免费观看不卡| 成人毛片a级毛片在线播放| 性欧美人与动物交配| 国产黄a三级三级三级人| 变态另类成人亚洲欧美熟女| 最新在线观看一区二区三区| 亚洲18禁久久av| 国产乱人伦免费视频| 国产日本99.免费观看| 精品国产三级普通话版| 免费观看精品视频网站| 制服丝袜大香蕉在线| 全区人妻精品视频| 午夜爱爱视频在线播放| 婷婷精品国产亚洲av在线| 国产精品久久久久久久久免| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 国产高清激情床上av| 久久久精品大字幕| 嫩草影院精品99| 亚洲经典国产精华液单| 别揉我奶头~嗯~啊~动态视频| 国产精品无大码| 亚洲av成人精品一区久久| 欧美黑人欧美精品刺激| 最近在线观看免费完整版| .国产精品久久|