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

    氣候變化對滇西北碧塔海流域景觀演變的影響

    2016-04-13 02:09:58喻慶國陳先剛楊宇明
    生態(tài)學(xué)報 2016年5期
    關(guān)鍵詞:演變氣候變化景觀

    姚 平,喻慶國,陳先剛,*,楊宇明

    1 西南林業(yè)大學(xué)環(huán)境科學(xué)與工程學(xué)院,昆明 650224

    2 西南林業(yè)大學(xué),國家高原濕地研究中心,昆明 650224

    3 云南省林業(yè)科學(xué)院,昆明 650201

    氣候變化對滇西北碧塔海流域景觀演變的影響

    姚平1,喻慶國2,陳先剛1,*,楊宇明3

    1 西南林業(yè)大學(xué)環(huán)境科學(xué)與工程學(xué)院,昆明650224

    2 西南林業(yè)大學(xué),國家高原濕地研究中心,昆明650224

    3 云南省林業(yè)科學(xué)院,昆明650201

    摘要:氣候變化是景觀演變的重要驅(qū)動力,高海拔地帶生態(tài)系統(tǒng)對氣候變化的響應(yīng)較為敏感和迅速。選擇地處青藏高原東南緣,受人為干擾相對較小的碧塔海流域為研究對象,利用1958—2011年氣象資料以及1955、1974、1981、1994、2005和2011年6個年代的景觀類型面積,采用氣候傾向率估計、Mann-Kendall突變檢驗方法和Pearson相關(guān)系數(shù)等統(tǒng)計方法,研究了香格里拉縣氣候變化背景下碧塔海流域景觀演變特征。結(jié)果表明:(1)1958—2011年來,滇西北香格里拉縣多年平均氣溫為5.9℃,多年平均活動積溫為2146.1℃,多年平均年極端高溫為24.3℃,多年平均年極端低溫為-18.6℃,多年平均降水量為631.7mm。氣溫、活動積溫和年極端低溫呈顯著上升趨勢,近54年來氣溫平均上升了1.94℃,升溫速率遠高于云南和全國水平。年極端高溫和降水量呈增加趨勢,但增加不明顯。滇西北高原碧塔海流域背景氣候呈現(xiàn)顯著的變暖趨勢。(2)1955—2011年來,碧塔海流域自然景觀演變規(guī)律為草甸和灌草叢向有林地演變,氣溫、活動積溫和年極端低溫升高趨勢與有林地景觀面積的增加趨勢大體一致,與草甸和灌草叢變化趨勢相反。(3)氣溫、活動積溫和年極端低溫與碧塔海流域草甸、灌草叢和有林地的景觀面積變化具有顯著的相關(guān)性,是流域景觀演變的主導(dǎo)氣候要素,在草甸和灌草叢演變成有林地的過程中發(fā)揮著重要作用。極端高溫和降水量對流域景觀演變的作用相對較弱。

    關(guān)鍵詞:氣候變化;碧塔海流域;景觀;演變

    據(jù)政府間氣候變化專門委員會(IPCC)第五次報告指出,近132年(1880—2012年)來地表平均溫度增加了0.85℃[1],中國近百年來氣候也發(fā)生了與全球氣候變化總趨勢基本一致的明顯變化[2]。氣候變化是景觀演變的重要自然驅(qū)動力之一,植被的形成是植物對特定氣候長期適應(yīng)的結(jié)果[3],氣候變化中水、熱條件的改變直接或間接地控制著各種生態(tài)過程和生態(tài)系統(tǒng)生產(chǎn)力,使植被類型改變,導(dǎo)致景觀格局的變化[4- 8]。高海拔地帶的生態(tài)系統(tǒng)對氣候變化的響應(yīng)最為敏感和迅速[9- 10],氣候變化對高原景觀格局影響的研究日益受到關(guān)注。

    地處云南省西北部 (簡稱滇西北)高原的碧塔海是滇西湖群組成之一[11],在區(qū)域獨特的地質(zhì)地貌和氣候背景下,碧塔海流域發(fā)育了以湖泊、沼澤、沼澤化草甸、草甸、森林等為主的滇西北典型高原景觀類型[12]。由于碧塔海地處金沙江流域中上游,且位于全球34個生物多樣性熱點地區(qū)之一的“中國西南山地”重要地帶[13],而備受關(guān)注。目前,對碧塔海景觀空間格局[12]、濕地生態(tài)環(huán)境[14]、濕地植被[15]、群落演替[16]等方面已有相關(guān)研究,但對在氣候變化背景下,碧塔海流域景觀格局演變的研究尚未見報道。因此,本文針對低緯高原、人為干擾較小的滇西北碧塔海流域,開展氣候變化背景下的景觀演變研究,以期從景觀層面為碧塔海流域的有效保護與合理利用,以及進一步研究低緯高原生態(tài)系統(tǒng)的變化規(guī)律和驅(qū)動因子提供科學(xué)依據(jù)。

    1研究區(qū)概況

    碧塔海流域地處青藏高原東南緣橫斷山脈三江并流區(qū)東部,位于滇西北迪慶藏族自治州香格里拉縣縣城東側(cè),距縣城約30 km,屬于金沙江水系[11]。流域介于北緯27°47′44″—27°51′19″和東經(jīng)99°57′15″—100°00′31″之間,面積1945.0 hm2,其內(nèi)最高海拔3967 m,最低海拔3538 m,高差429 m[12],為高海拔、低緯度地帶[17]。流域由碧塔海湖泊水域、山間平緩谷地和湖泊匯水面山組成[18]。碧塔海是滇西北一個人為干擾較小的[11]、保存最為完整的封閉型高原淡水湖泊[19],由于其匯水面山分布有茂密的長苞冷杉林、油麥吊云杉林、川滇高山櫟林等原始森林[20],起到了極好的水土保持和水源涵養(yǎng)作用。碧塔海湖水補給穩(wěn)定,湖面變化較小,水位較為穩(wěn)定[11]。碧塔海流域是三江并流國家級風(fēng)景名勝區(qū)、碧塔海省級自然保護區(qū)、三江并流世界自然遺產(chǎn)地、碧塔海國際重要濕地和香格里拉普達措國家公園的重要組成部分[12]。綜上可見,該區(qū)是開展氣候變化背景下景觀自然演變研究的理想場所,具有重要的研究價值。

    2資料及方法

    2.1資料

    碧塔海流域內(nèi)沒有長期氣象觀測資料,考慮到自然驅(qū)動因子常常是在較大的時空尺度上作用于景觀格局[21],本文旨在研究氣候變化背景下碧塔海流域景觀演變特征,因此采用距離碧塔海流域最近的國家基準(zhǔn)氣象站(香格里拉縣氣象觀測站),自1958年建站至2011年來近54年的年平均氣溫、年活動積溫、年極端高溫、年極端低溫、年降水量(文中簡稱為氣溫、活動積溫、年極端高溫、年極端低溫和降水量)和日平均氣溫資料作為碧塔海流域的氣候背景資料。植物通常在日平均溫度上升到5℃的時候開始發(fā)育,文中活動積溫為年內(nèi)日平均氣溫穩(wěn)定通過5℃期間的日平均氣溫之和[22]。

    在本研究中,反映碧塔海流域景觀演變的景觀類型面積數(shù)據(jù)取自本文作者的前期研究成果,1955—2011年6個年代的景觀類型面積數(shù)據(jù)和景觀類型圖參閱文獻[12],在本文中不再列出相關(guān)景觀類型面積數(shù)據(jù)和附相關(guān)景觀類型圖。

    2.2研究方法

    2.2.1氣候變化特征

    (1)氣候傾向率估計方法

    (2)Mann-Kendall突變檢驗方法(簡稱M-K方法)

    該方法是一種非參數(shù)統(tǒng)計檢驗方法,其優(yōu)點是不需要樣本遵從一定分布,也不受少數(shù)異常值的影響,更適用于類型變量和順序變量。其方法如下:對于具有n個樣本量的時間序列x,構(gòu)造一秩序列:

    (1)

    (2)

    式中,UF1=0,而E(sk)、var(sk)是累計數(shù)sk的均值和方差。在x1,x2,,…,xn相互獨立且具有相同連續(xù)分布時,有:

    (3)

    2.2.2氣候要素與景觀面積的相關(guān)性

    皮爾遜(Pearson)相關(guān)系數(shù)是描述兩個隨機變量線性相關(guān)的統(tǒng)計量,設(shè)有兩列時間變量x=x1,x2,…,xn和 y=y1,y2,…,yn,n為樣本容量,相關(guān)系數(shù)的計算公式為:

    (4)

    Pearson相關(guān)系數(shù)r的取值在-1.0—1.0之間,當(dāng)r>0時,表明兩變量呈正相關(guān),越接近1.0,正相關(guān)越顯著;當(dāng)r<0時,表明兩變量呈負(fù)相關(guān),越接近-1.0,負(fù)相關(guān)越顯著;當(dāng)r=0時,表示兩變量相互獨立。據(jù)統(tǒng)計學(xué)大樣本定理[25],樣本量大于30才有統(tǒng)計意義,由于景觀類型面積樣本容量較小,用無偏相關(guān)系數(shù)對Pearson相關(guān)系數(shù)加以矯正,無偏相關(guān)系數(shù)記為r*,計算公式為:

    (5)

    最后利用t檢驗方法對矯正后的相關(guān)系數(shù)進行顯著性檢驗[24]。采用無偏相關(guān)系數(shù)探討氣候要素與碧塔海流域各景觀類型面積之間的相關(guān)性。本文采用Matlab 7.01軟件進行計算分析。

    3結(jié)果與分析

    3.1背景氣候變化特征

    近54年來香格里拉縣多年平均氣溫為5.9℃;多年平均活動積溫為2146.1℃;多年平均年極端高溫為24.3℃,且最高值出現(xiàn)在2005年和2006年,達26.0℃;多年平均年極端低溫為-18.6℃,最低值出現(xiàn)在1982年,低至-27.4℃;多年平均降水量為631.7mm。

    利用氣候傾向率估計方法對各氣候因子的時間序列做線性傾向分析結(jié)果見圖1。從圖1中可知香格里拉縣各氣候要素隨時間變化的總趨勢均為增加,氣溫、活動積溫、年極端高溫、年極端低溫和降水量的氣候變化傾向率分別為0.36℃/10a、56.5℃/10a 、0.15℃/10a、0.98℃/10a和8.8mm/10a,近54年來氣溫平均上升了約1.94℃,與全省和全國相比升溫趨勢一致,氣溫升溫幅度和升溫速率均高于云南[26]和全國[2]的平均狀況(云南近48年平均上升了0.66℃,增溫速率為0.14℃/10a,全國近50年來平均氣溫上升了1.3℃,增溫速率為0.25℃/10a),年極端高溫上升了約0.81℃,年極端低溫上升了約5.29℃,活動積溫增加了305.1℃,降水量增加了約47.5mm。

    圖1 1958—2011年滇西北香格里拉縣氣候要素變化趨勢Fig.1 The trend of climate elements over Shangri-La County of northwestern Yunnan from 1958 to 2011直線:線性傾向擬合線

    采用M-K方法對香格里拉縣各氣候要素的變化規(guī)律和突變特征進行檢驗結(jié)果見圖2。結(jié)果發(fā)現(xiàn)氣溫M-K計量曲線中(圖2),1979年以前UF值小于零,1979年以后大于零,1988年通過α=0.05的信度檢驗,表明1958—2011年間,1979年以前氣溫有所下降,但趨勢不明顯,1979年以后氣溫呈上升趨勢,且1988年以后上升趨勢極為顯著。香格里拉縣活動積溫的M-K計量曲線中(圖2),1981年之前UF值小于零, 1981年以后UF值大于零,1993年通過α=0.05的信度檢驗,UF和UB的交點出現(xiàn)在1988和1990年,表明活動積溫在1980年以后為上升趨勢,且上升和突變趨勢顯著,突變開始的時間大致為1988年左右。年極端高溫M-K計量曲線中(圖2),1976年以后UF值大于零,在近54年中上升趨勢未通過α=0.05的信度檢驗,表明年極端高溫在1976年以后為上升趨勢,但上升趨勢不明顯。年極端低溫M-K計量曲線中(圖2),1958年以來UF值均大零,1990年通過α=0.05的信度檢驗,UF和UB的交點出現(xiàn)在1986年,表明年極端低溫自1958年以來主要為升高趨勢,突變開始的時間為1986年,且上升和突變趨勢顯著。1989年以前降水量的UF曲線在零線附近擺動,在1989年之后均位于零線以上,表明降水量在20世紀(jì)80年代以后有所增加,但增加趨勢不明顯(圖2)。

    圖2 1958—2011年滇西北香格里拉縣氣候要素M-K突變檢驗Fig.2 The M-K mutation test of climate elements over Shangri-La County of northwestern Yunnan from 1958 to 2011粗實線:UF;粗虛線:UB;細實線:檢驗臨界值;點虛線:零值線

    M-K方法檢驗表明滇西北香格里拉縣活動積溫在1988年附近存在顯著的突變趨勢,突變前1958—1988年的多年平均活動積溫為2072.6℃,突變后1989—2011年的多年平均活動積溫為2239.7℃,突變前后平均上升了167.1℃。年極端低溫突變開始的時間為1986年,突變前1958—1986年的多年平均年極端低溫為-20.2℃,突變后1987—2011年的多年平均年極端低溫為-16.8℃,平均上升了3.4℃,香格里拉縣年極端低溫突變的時間與全國年極端低溫顯著升高的時間基本一致[27]。

    綜上所述,1958—2011年近54年來滇西北香格里拉縣氣溫、活動積溫和年極端低溫大體上呈顯著上升趨勢,活動積溫和年極端低溫突變趨勢顯著,年極端高溫和降水量稍有上升或增加,但上升或增加趨勢不明顯。滇西北碧塔海流域背景氣候具有明顯的變暖趨勢。

    3.2景觀面積變化特征

    碧塔海流域景觀類型有湖泊水域、沼澤、沼澤化草甸、草甸、有林地、寒溫性稀樹灌草叢(簡稱灌草叢)和公路,其中有林地為流域的優(yōu)勢景觀類型。1955—2011年,除湖泊水域、沼澤和沼澤化草甸沒有變化,草甸和灌草叢分別減少了11.08 hm2和135.19 hm2;有林地增加了118.24 hm2,由草甸和灌草叢轉(zhuǎn)變而來;公路自2005年開始修建,增加了28.03 hm2,由草甸、灌草叢和有林地轉(zhuǎn)變而成??傮w上流域景觀演變表現(xiàn)出以自然驅(qū)動力為主,人為驅(qū)動力為輔的景觀演變驅(qū)動力格局[12]。

    圖3 滇西北碧塔海流域景觀面積距平標(biāo)準(zhǔn)化序列 Fig.3 The standardized anomalies sequence of landscape area over Bita lake watershed of northwestern Yunnan

    選取碧塔海流域內(nèi)景觀類型發(fā)生變化的草甸、有林地、灌草叢和公路景觀面積進行距平標(biāo)準(zhǔn)化,以便在同一水平上分析其隨時間的變化趨勢(圖3)。圖3可知,1955—2011年近56年來草甸和灌草叢的景觀面積隨時間推移一直呈減少的趨勢,1955—2005年間有林地的景觀面積隨時間推移大體上一直呈增加的趨勢,2005—2011年間,因修建公路使草甸、灌草叢和有林地被占用,有林地的景觀面積有所下降。

    3.3景觀演變與氣候要素的相關(guān)性

    景觀演變是自然驅(qū)動因子和人為驅(qū)動因子共同作用的結(jié)果[28],自然驅(qū)動因子包括氣候、地質(zhì)、地貌、地形、土壤、水文和生物作用。1955—2011年的56年間,碧塔海流域的地質(zhì)、地貌、地形、土壤等自然因子基本沒有變化,也沒有發(fā)生過自然災(zāi)害,且人為干擾較小[12],氣候則呈現(xiàn)顯著的變化趨勢,這必然會導(dǎo)致景觀格局發(fā)生變化,在景觀類型層面直接表現(xiàn)為植被群落演替,這主要由于群落總是朝向和當(dāng)?shù)厣?包括氣候)保持協(xié)調(diào)和平衡的方向演替。也就是說,氣候是碧塔海流域自然景觀類型演變的主要自然驅(qū)動因子。

    由于碧塔海流域遙感數(shù)據(jù)始于1955年,而氣候要素在1958年才有最早的觀測記錄,故以離1955年最近連續(xù)3a(1958、1959和1960年)的氣候要素平均值代替1955年缺測的氣候資料,其它時段則采用同期1974,1981,1994,2004和2011年的氣候資料與碧塔海流域自然景觀類型(草甸、有林地、灌草叢)面積進行相關(guān)分析(表1),篩選影響景觀演變的主導(dǎo)氣候要素。

    表1 碧塔海流域景觀類型面積與氣候要素的無偏相關(guān)系數(shù)

    *,**,***分別為通過α=0.1,0.05,0.01的顯著t檢驗

    碧塔海流域景觀類型面積與氣候要素的無偏相關(guān)系數(shù)表(表1)中加粗的相關(guān)系數(shù)值均通過顯著性檢驗,近年來氣候因子與碧塔海流域有變化的自然景觀類型草甸、有林地、灌草叢之間存在以下關(guān)系:

    草甸和灌草叢面積與氣溫、活動積溫和年極端低溫呈顯著的負(fù)相關(guān),而有林地面積與氣溫、活動積溫和年極端低溫呈顯著的正相關(guān),表明當(dāng)氣溫、活動積溫或年極端低溫上升(降低)時,草甸和灌草叢的面積減小(增加),而有林地的面積增加(減小)。草甸和灌草叢景觀面積與年極端高溫呈負(fù)相關(guān),與降水量呈正相關(guān),而有林地景觀面積與年極端高溫呈正相關(guān),與降水量呈負(fù)相關(guān),但相關(guān)性都不明顯。

    根據(jù)前文3.1的分析結(jié)果可知,1958—2011年香格里拉縣氣溫總體呈上升趨勢(圖1),1958—1979年氣溫有所下降(圖2),這主要是受到香格里拉縣1958年才有氣候觀測資料的限制,難以確知1958年以前的氣候變化情況。據(jù)萬云霞等的研究,香格里拉縣所處的云南縱向嶺谷區(qū)的氣溫,存在明顯冷、暖位相交替的年代際變化特征,在1901—2002年近百年來總體呈現(xiàn)增溫趨勢[29]。結(jié)合萬云霞等的研究以及圖1和圖2可推知,1901年以來香格里拉縣的氣溫在百年時間尺度上總體呈現(xiàn)增溫趨勢。對照圖1、圖2和圖3,香格里拉縣的氣溫、活動積溫和年極端低溫隨時間大體呈明顯上升趨勢,與有林地面積隨時間明顯增加的趨勢是一致的,而與草甸和灌草叢面積隨時間明顯減少的趨勢相反。這也充分說明氣溫、活動積溫和年極端低溫在有林地增加,以及草甸、灌草叢減少中起著重要作用。

    以上分析表明,近56年來,氣溫、活動積溫和年極端低溫與碧塔海流域草甸、灌草叢和有林地面積關(guān)系密切,降水量和年極端高溫與流域景觀面積相關(guān)性較弱,氣溫、活動積溫和年極端低溫是影響碧塔海流域自然景觀類型演變的主導(dǎo)氣候要素。

    氣溫對高海拔地區(qū)樹木的生長具有重要的作用[5],碧塔海流域有林地的主要植被類型為冷杉、云杉、高山櫟林等[20],已有研究發(fā)現(xiàn)長期增溫使高海拔地區(qū)冷杉和云杉幼苗的光合作用增強,有利于促進幼苗的生長和生物量的積累[7,30- 31];其次增溫能延長植物的生長季[32],加快土壤有機質(zhì)的分解速率和養(yǎng)分的可利用性,間接促進樹木生長[33]。熱量是影響樹木生長發(fā)育的重要氣候因子之一,樹木光合作用的下限溫度一般為5℃[34],對于受低溫限制的高寒生態(tài)系統(tǒng),較短的生長期是限制植物生物量生產(chǎn)的關(guān)鍵因素之一[35],活動積溫升高意味著植物生長期的延長[36],有利于樹木的生長發(fā)育。碧塔海流域背景氣候相對較寒冷(多年平均年極端低溫為-18.6℃,最低值低至-27.4℃),長期適應(yīng)高寒等環(huán)境的土壤微生物,在溫度稍微升高時,其生物活動便會急劇加強[37],使土壤中酶活性升高[38],有利于凋落物和土壤有機物的分解、轉(zhuǎn)換和養(yǎng)分元素的釋放,提高森林土壤肥力[39],年極端低溫升高可能通過增強土壤生物活性、土壤養(yǎng)分循環(huán)等方式,間接影響樹木的生長。氣溫和活動積溫的升高,有利于流域樹種的生長、發(fā)育,有利于流域內(nèi)草甸和灌草叢向有林地的演變。

    降水量和極端高溫對碧塔海流域景觀演變的作用相對較弱,主要原因是近年來降水量有所增加,但增加趨勢不明顯,且碧塔海湖泊水域面積沒有發(fā)生變化[12],碧塔海四周分布有郁郁蔥蔥的原始森林[20],其森林水源涵養(yǎng)和水土保持功能較好,氣溫升高能引起湖泊水體溫度的升高,導(dǎo)致蒸發(fā)量增大,使空氣濕度增大,使得降水量對流域植被的生態(tài)促進作用有限。另外,由于流域氣候背景的年極端高溫值并不高(多年平均年極端高溫為24.3℃,近54年來最高為26℃),升高趨勢不明顯,對流域植被的生長發(fā)育影響不大。

    4結(jié)論與討論

    (1)近54年來地處低緯度、高海拔的滇西北香格里拉縣氣溫上升了約1.94℃,升溫速率為0.36℃/10a,上升趨勢遠高于全國水平?;顒臃e溫和年極端低溫增溫和突變趨勢明顯,并分別于1988和1986年附近產(chǎn)生突變,突變前后活動積溫平均增加了167.1℃,年極端低溫上升了3.4℃。年極端高溫和降水量稍有上升或增加,但上升或增加趨勢不明顯。滇西北高原碧塔海流域背景氣候呈現(xiàn)顯著的變暖趨勢,且變暖程度遠高于全國水平。

    (2)在1955—2011年,碧塔海流域內(nèi)自然景觀演變規(guī)律為草甸和灌草叢向有林地演變,且有林地景觀面積隨時間的增加趨勢與氣溫、活動積溫和年極端低溫隨時間的上升趨勢大體一致,草甸和灌草叢景觀面積隨時間的變化特征則與氣溫、活動積溫和年極端低溫的變化趨勢相反。

    (3)氣溫、活動積溫和年極端低溫與碧塔海流域草甸、灌草叢和有林地的景觀面積變化具有顯著的相關(guān)性,是碧塔海流域自然景觀演變的主導(dǎo)氣候要素,降水量和極端高溫的驅(qū)動作用相對較弱。氣溫、活動積溫和年極端低溫在草甸和灌草叢演變成有林地的過程中發(fā)揮著重要作用。

    (4)景觀演變是一個極為復(fù)雜的過程,受自然驅(qū)動因子和人為驅(qū)動因子共同的影響。1955—2011年在碧塔海流域人為驅(qū)動因子主要為兩方面因素,其一為碧塔海流域周邊藏民56年間的持續(xù)性放牧習(xí)俗,由于沒有出現(xiàn)過牧現(xiàn)象,放牧沒有造成流域內(nèi)景觀格局的變化;其二為道路建設(shè),其對流域景觀格局演變的影響表現(xiàn)為2005—2011年,因修建公路使草甸、灌草叢和有林地被占用,三類景觀面積均有所減小,因此道路建設(shè)是碧塔海流域景觀演變的主要人為驅(qū)動因子[12]。綜上所述,在氣候變化的背景下,近56年來氣候變化是碧塔海流域自然景觀格局演變的重要驅(qū)動力,人為驅(qū)動因子的作用相對較弱。

    (5)據(jù)IPCC第五次評估報告,與1986—2005年相比,2016—2035年全球平均地表溫度將升高0.3—0.7℃,2081—2100年將升高0.3—4.8℃[1]。據(jù)此,滇西北碧塔海流域的背景氣候可能將進一步變暖,未來變暖趨勢是否會導(dǎo)致流域其它景觀類型之間的相互演變,以及氣候要素影響景觀類型演變的具體過程,都是值得繼續(xù)研究的內(nèi)容。另外,植被覆蓋類型的改變,將會改變地表反照率、表面粗糙度、土壤濕度等地表屬性,影響地表與大氣間的水、熱交換過程,從而影響氣候變化,也即碧塔海流域景觀格局的演變與局地氣候變化之間的正、負(fù)反饋機制也是以后需要進一步研究的內(nèi)容。

    參考文獻(References):

    [1] Intergovernmental Panel on Climate Change. Climate Change 2013: The Physical Science Basis: The Working Group I Contribution to the IPCC Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 2013.

    [2]丁一匯, 任國玉. 中國氣候變化科學(xué)概論. 北京: 氣象出版社, 2008: 61- 86.

    [3]張新時. 研究全球變化的植被-氣候分類系統(tǒng). 第四紀(jì)研究, 1993, 13(2): 157- 169.

    [4]Allen C D, Breshers D D. Drought-induced shift of a forest-woodland ecotone: rapid landscape response to climate variation. Proceedings of the National Academy of Sciences of the United States of America, 1998, 95(25): 14839- 14842.

    [5]K?mer C, Paulsen J. A world-wide study of high altitude treeline temperatures. Journal of Biogeography, 2004, 31(5): 713- 732.

    [6]劉國華, 傅伯杰. 全球氣候變化對森林生態(tài)系統(tǒng)的影響. 自然資源學(xué)報, 2001, 16(1): 71- 78.

    [7]Yin H J, Liu Q, Lai T. Warming effects on growth and physiology in the seedlings of the two conifersPiceaasperataandAbiesfaxonianaunder two contrasting light conditions. Ecological Research, 2008, 23(2): 459- 469.

    [8]李杰, 胡金明, 董云霞, 袁寒. 1994—2006年滇西北納帕海流域及其濕地景觀變化研究. 山地學(xué)報, 2010, 28(2): 247- 256.

    [9]Jorgenson M T, Racine C H, Walters J C, Osterkamp T E. Permafrost degradation and ecological changes associated with a warming climate in central Alaska. Climate Change, 2001, 48(4): 551- 579.

    [10]Nelson F E, Aniamovand O A, Shiklomanov N I. Climate change and hazard zonation in the circum-arctic permafrost regions. Natural Hazards, 2002, 26(3): 203- 225.

    [11]楊嵐, 李恒, 楊曉君. 云南濕地. 北京: 中國林業(yè)出版社, 2010: 23- 23, 34- 40.

    [12]喻慶國, 董躍宇, 楊宇明, 劉朝蓬, 袁磊, 李昊民. 滇西北碧塔海流域景觀演變規(guī)律及其驅(qū)動因子. 西北林學(xué)院學(xué)報, 2014, 29(4): 160- 168.

    [13]Mittermeier R A, Gil P R, Hofman M, Pilgrim J, Brooks T, Mittermeier C G, Lamoreux J, da Fonseca G A B, Seligmann P A, Ford H. Hotspots Revisited: Earth′s Biologically Richest and Most Endangered Terrestrial Ecoregions. USA, Chicago: The University of Chicago Press, 2005: 3- 78.

    [14]唐明艷, 楊永興. 旅游干擾下滇西北高原湖濱濕地植被及土壤變化特征. 應(yīng)用生態(tài)學(xué)報, 2014, 25(5): 1283- 1292.

    [15]尹五元, 杜凡, 張大才. 碧塔海自然保護區(qū)珍稀保護植物研究. 西南林學(xué)院學(xué)報, 2004, 24(1): 20- 22.

    [16]韓大勇, 楊永興, 楊楊. 滇西北高原碧塔海湖濱沼澤植物群落分布與演替. 生態(tài)學(xué)報, 2013, 33(4): 2236- 2247.

    [17]和麗忠. 香格里拉年鑒2010. 昆明: 云南科技出版社, 2010: 49- 49.

    [18]董躍宇, 喻慶國, 劉朝蓬, 余紅忠, 趙建林. 基于最佳空間粒度的碧塔海流域景觀格局分析. 西部林業(yè)科學(xué), 2013, 42(3): 45- 53.

    [19]楊學(xué)光. 香格里拉縣林業(yè)志. 昆明: 云南民族出版社, 2006: 128- 128.

    [20]劉德隅, 劉中天, 錢德仁, 湯家生, 王建皓, 高平龍, 施蔚云. 云南自然保護區(qū). 北京: 中國林業(yè)出版社, 1989: 260- 268.

    [21]郭程軒, 徐頌軍. 基于3S與模型方法的濕地景觀動態(tài)變化研究述評. 地理與地理信息科學(xué), 2007, 23(5): 86- 90.

    [22]賀慶棠, 陸佩玲. 氣象學(xué)(第3版). 北京: 中國林業(yè)出版社, 2010: 84- 86, 263- 265.

    [23]Jones P D. Hemispheric surface air temperature variations: recent trends and an update to 1987. Journal of Climate, 1988, 1(6): 654- 660.

    [24]魏鳳英. 現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)(第2版). 北京: 氣象出版社, 2007: 18- 19, 30- 31, 63- 66.

    [25]王梓坤. 概率論基礎(chǔ)及其應(yīng)用. 北京: 科學(xué)出版社, 1976.

    [26]段旭, 陶云. 云南近50年來的氣候變化. 熱帶氣象學(xué)報, 2012, 28(2): 243- 250.

    [27]陳隆勛, 朱文琴, 王文, 周秀驥, 李維亮. 中國近45年來氣候變化的研究. 氣象學(xué)報, 1998, 56(6): 257- 271.

    [28]李傳哲, 于福亮, 劉佳. 分水后黑河干流中游地區(qū)景觀動態(tài)變化及驅(qū)動力. 生態(tài)學(xué)報, 2009, 29(11): 5832- 5842.

    [29]萬云霞, 張萬誠, 肖子牛. 近百年云南縱向嶺谷區(qū)氣溫的時空變化特征. 自然災(zāi)害學(xué)報, 2009, 18(5): 183- 188.

    [30]楊兵, 王進闖, 張遠彬. 長期模擬增溫對岷江冷杉幼苗生長與生物量分配的影響. 生態(tài)學(xué)報, 2010, 30(21): 5994- 6000.

    [31]Wang J H, Duan B L, Zhang Y B. Effects of experimental warming on growth, biomass allocation, and needle chemistry ofAbiesfaxonianain even-aged monospecific stands. Plant Ecology, 2012, 213(1): 47- 55.

    [32]Bronson D R, Gower S T, Myron T, Van Herk I. Effect of ecosystem warming on boreal black spruce bud burst and shoot growth. Global Change Biology, 2009, 15(6): 1534- 1543.

    [33]Str?mgren M, Linder S. Effects of nutrition and soil warming on stemwood production in a boreal Norway spruce stand. Global Change Biology, 2002, 8(12): 1194- 1204.

    [34]吉林林業(yè)學(xué)校, 四川林業(yè)學(xué)校, 陜西農(nóng)林學(xué)校. 森林學(xué). 北京: 中國林業(yè)出版社, 1981.

    [35]Gugerli F, Bauert M R. Growth and reproduction ofPolygonumviviparumshow weak responses to experimentally increased temperature at a Swiss alpine site. Botanica Helvetica, 2001, 111(2): 169- 180.

    [36]徐興奎, 王小桃, 金曉青. 中國區(qū)域1960—2000年活動積溫年代變化和地表植被的適應(yīng)性調(diào)整. 生態(tài)學(xué)報, 2009, 29(11): 6042- 6050.

    [37]趙亮, 徐世曉, 李英年, 唐艷鴻, 趙新全, 古松, 杜明遠, 于貴瑞. 青藏高原矮嵩草草甸和金露梅灌叢草甸CO2通量變化與環(huán)境因子的關(guān)系. 西北植物學(xué)報, 2006, 26(1): 133- 142.

    [38]Dilly O, Munch J C. Microbial biomass content, basal respiration and enzyme activities during the course of decomposition of leaf litter in a black alder (Alnusglutinosa(L.) Gaertn.) forest. Soil Biology & Biochemistry, 1996, 28(8): 1073- 1081.

    [39]楊萬勤, 王開運. 森林土壤酶的研究進展. 林業(yè)科學(xué), 2004, 40(2):152- 159.

    Landscape evolution characteristics of the Bita Lake watershed over northwestern Yunnan province under the background of climate changes

    YAO Ping1, YU Qingguo2, CHEN Xiangang1,*, YANG Yuming3

    1DepartmentofEnvironmentalScienceandEngineering,SouthwestForestryUniversity,Kunming650224,China2SouthwestForestryUniversity,NationalPlateauWetlandsResearchCenter,Kunming650224,China3YunnanAcademyofForestry,Kunming650201,China

    Abstract:Climate change is one of the major driving forces for landscape evolution. The ecosystem in high altitude area has the most sensitive and fast response to climate change. In this study, we focused on the Bita lake watershed of northwestern Yunnan located in the Shangri-la county and in the southeast margin of Three Rivers Parallel Area over Qinghai-Tibet plateau mountains. The study area with small anthropogenic interference, is the ideal area to identify the coupling relationship between climate change and landscape change. Based on the meteorological data from 1958 to 2011, and the landscape area indices of class level in the years of 1955, 1974, 1981, 1994, 2005, 2011, the characteristics of landscape evolution over the Bita lake watershed under the background of climate change were studied using climate tendency estimation method, Mann-Kendall mutation method and Pearson correlation coefficient method. The results were as follow: (1) The values of mean annual temperature, mean annual active accumulated temperature, mean annual extreme high temperature, mean annual extreme low temperature, average annual precipitation were 5.9℃, 2146.1℃, 24.3℃, -18.6℃ and 631.7mm from 1958 to 2011, respectively. The temperature, active accumulated temperature and annual extreme low temperature increased significantly in the Shangri-La of northwestern Yunnan, while the increases of annual extreme high temperature and precipitation weren′t significant. The average temperature was elevated about 1.94 ℃ from 1958 to 2011, and the increasing rate was much higher than that in Yunnan and China. Active accumulated temperature and annual extreme low temperature were shown to be significant warming mutation trend with an increase of 167.1℃ and 3.4℃ after the mutation in the years of 1988 and 1986, respectively. (2) The mechanism of landscape evolution was in a succession from meadow and scrub-grassland to forestland in the Bita lake watershed from1955 to 2011. The rising trends of temperature, active accumulated temperature and annual extreme low temperature were consistent with the increasing trend of forestland landscape area, but were contrary to the decreasing trend of meadow and scrub-grassland landscape area. (3) Significant correlation was found between temperature, active accumulated temperature, annual extreme low temperature and natural landscape area of meadow, scrub-grassland and forestland from 1955 to 2011. When temperature,active accumulated temperature and annual extreme low temperature were increased quickly, the area of meadow and scrub-grassland was decreased, whereas the area of forestland was increased. Temperature, active accumulated temperature and annual extreme low temperature were the main driving factors of landscape evolution over the Bita lake watershed of northwestern Yunnan, while the ecological effects of annual extreme high temperature and precipitation on landscape evolution were relatively lesser. Climate warming was conducive to the growth of the watershed dominant tree, which played an important role in formation of a succession from meadow and irrigation grass to forestland in Bita lake watershed.

    Key Words:climate change; the Bita lake watershed; landscape; evolution

    基金項目:云南省優(yōu)勢特色重點學(xué)科(生態(tài)學(xué)); 國家自然科學(xué)基金(U0933601); 云南玉溪森林生態(tài)系統(tǒng)定位研究站項目

    收稿日期:2014- 07- 10; 網(wǎng)絡(luò)出版日期:2015- 07- 22

    DOI:10.5846/stxb201407101407

    *通訊作者Corresponding author.E-mail: ch-x-g@vip.sina.com

    姚平,喻慶國,陳先剛,楊宇明.氣候變化對滇西北碧塔海流域景觀演變的影響.生態(tài)學(xué)報,2016,36(5):1453- 1461.

    Yao P, Yu Q G, Chen X G, Yang Y M.Landscape evolution characteristics of the Bita Lake watershed over northwestern Yunnan province under the background of climate changes.Acta Ecologica Sinica,2016,36(5):1453- 1461.

    猜你喜歡
    演變氣候變化景觀
    《應(yīng)對氣候變化報告(2022)》發(fā)布
    景觀別墅
    火山塑造景觀
    包羅萬象的室內(nèi)景觀
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    異類女性故事在中國古代文學(xué)作品中的演變
    文教資料(2016年20期)2016-11-07 11:23:34
    美國公共住房退出管理中的兩難抉擇及啟示
    談中國戲曲藝術(shù)的發(fā)展進程
    戲劇之家(2016年19期)2016-10-31 17:41:38
    中國證券市場制度演變及制度性風(fēng)險分析
    国产精品久久电影中文字幕 | 免费在线观看日本一区| 久久精品亚洲精品国产色婷小说| 在线观看免费高清a一片| 热99国产精品久久久久久7| 国产精品久久久久久精品古装| 国产xxxxx性猛交| 亚洲精品国产精品久久久不卡| 高清黄色对白视频在线免费看| x7x7x7水蜜桃| netflix在线观看网站| 黄色 视频免费看| 国产精品国产av在线观看| 亚洲精品中文字幕一二三四区| 很黄的视频免费| 久久青草综合色| 久久香蕉精品热| 久热这里只有精品99| 国产有黄有色有爽视频| 久久精品国产99精品国产亚洲性色 | 国产男靠女视频免费网站| 日本黄色日本黄色录像| 制服诱惑二区| 丰满迷人的少妇在线观看| 国产乱人伦免费视频| 欧美精品高潮呻吟av久久| 久久国产精品大桥未久av| 香蕉国产在线看| 久久ye,这里只有精品| 激情在线观看视频在线高清 | 国产精品一区二区在线不卡| 亚洲五月色婷婷综合| 欧美黑人精品巨大| 夜夜爽天天搞| 欧美中文综合在线视频| 成年人午夜在线观看视频| 村上凉子中文字幕在线| 极品教师在线免费播放| 国产麻豆69| 99国产精品99久久久久| 亚洲成人免费av在线播放| 国产精品自产拍在线观看55亚洲 | 成年动漫av网址| 亚洲精品中文字幕一二三四区| 国产成人欧美在线观看 | 亚洲精品自拍成人| 国产成人啪精品午夜网站| 欧美激情久久久久久爽电影 | 丝瓜视频免费看黄片| 日韩三级视频一区二区三区| 如日韩欧美国产精品一区二区三区| 欧美国产精品一级二级三级| 久久国产精品影院| 欧美黑人欧美精品刺激| 男人操女人黄网站| 亚洲欧美日韩高清在线视频| 天天躁夜夜躁狠狠躁躁| 久久中文字幕一级| 欧美另类亚洲清纯唯美| 中文亚洲av片在线观看爽 | 中文字幕高清在线视频| 久久久国产欧美日韩av| 欧美精品人与动牲交sv欧美| 欧美中文综合在线视频| 久久国产精品影院| 日本黄色日本黄色录像| 国产成人免费观看mmmm| 热99国产精品久久久久久7| 国产又爽黄色视频| 1024视频免费在线观看| 久久人妻av系列| 黄色成人免费大全| 亚洲成av片中文字幕在线观看| 嫁个100分男人电影在线观看| 校园春色视频在线观看| 亚洲国产欧美网| 精品少妇一区二区三区视频日本电影| 国产一区二区三区视频了| 国产无遮挡羞羞视频在线观看| 国产精品免费视频内射| 国产99白浆流出| 男女午夜视频在线观看| 91成人精品电影| 欧美精品高潮呻吟av久久| 人妻 亚洲 视频| 男女床上黄色一级片免费看| 亚洲精品国产精品久久久不卡| 国产精品av久久久久免费| cao死你这个sao货| 一区二区三区精品91| 欧美日韩福利视频一区二区| 精品无人区乱码1区二区| 国产成人欧美在线观看 | 亚洲av美国av| 国产精品久久视频播放| 丁香六月欧美| 国产精品久久久av美女十八| 国产成人精品久久二区二区91| 午夜久久久在线观看| 9191精品国产免费久久| 校园春色视频在线观看| 极品少妇高潮喷水抽搐| 午夜激情av网站| 一区在线观看完整版| av网站免费在线观看视频| 女警被强在线播放| 亚洲国产中文字幕在线视频| 怎么达到女性高潮| 色综合婷婷激情| 久久中文字幕一级| 婷婷丁香在线五月| 9色porny在线观看| 亚洲av成人一区二区三| 美女福利国产在线| 淫妇啪啪啪对白视频| 精品高清国产在线一区| 一级黄色大片毛片| 亚洲精品一卡2卡三卡4卡5卡| 12—13女人毛片做爰片一| 身体一侧抽搐| av线在线观看网站| 精品人妻在线不人妻| 人妻一区二区av| 美女高潮到喷水免费观看| 极品教师在线免费播放| 亚洲七黄色美女视频| 日韩大码丰满熟妇| 99精品久久久久人妻精品| 老司机午夜福利在线观看视频| 国产片内射在线| 亚洲av片天天在线观看| 精品国产亚洲在线| 国产不卡一卡二| 精品福利永久在线观看| 多毛熟女@视频| 久久香蕉精品热| 成年人免费黄色播放视频| 一进一出抽搐动态| 欧美人与性动交α欧美软件| 久久中文看片网| 欧美人与性动交α欧美精品济南到| 色综合婷婷激情| 亚洲第一av免费看| 后天国语完整版免费观看| 看黄色毛片网站| 国产精品久久久av美女十八| 欧美亚洲 丝袜 人妻 在线| 我的亚洲天堂| 精品亚洲成国产av| 一级片免费观看大全| 91老司机精品| videos熟女内射| 日本a在线网址| 高清在线国产一区| 婷婷精品国产亚洲av在线 | 中文字幕av电影在线播放| 亚洲精品在线观看二区| 超色免费av| 中文字幕av电影在线播放| 波多野结衣av一区二区av| 久久亚洲真实| 看片在线看免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲熟女精品中文字幕| 下体分泌物呈黄色| 久久久久精品人妻al黑| 久久草成人影院| 精品少妇一区二区三区视频日本电影| 久久中文字幕人妻熟女| 欧美久久黑人一区二区| 国产真人三级小视频在线观看| 这个男人来自地球电影免费观看| 99精品在免费线老司机午夜| 巨乳人妻的诱惑在线观看| 又大又爽又粗| 搡老熟女国产l中国老女人| 成年动漫av网址| 麻豆国产av国片精品| 高清欧美精品videossex| 亚洲精品久久成人aⅴ小说| 欧美乱妇无乱码| 黑人操中国人逼视频| 巨乳人妻的诱惑在线观看| 视频在线观看一区二区三区| 久久国产亚洲av麻豆专区| 午夜福利欧美成人| 欧美日韩一级在线毛片| 欧美国产精品va在线观看不卡| 亚洲精品一二三| 无人区码免费观看不卡| 又紧又爽又黄一区二区| 国产av精品麻豆| 久久久国产一区二区| 婷婷成人精品国产| 99re6热这里在线精品视频| 成年人黄色毛片网站| 咕卡用的链子| 91精品三级在线观看| 少妇的丰满在线观看| 精品国产超薄肉色丝袜足j| 老司机影院毛片| 黄色丝袜av网址大全| 亚洲专区中文字幕在线| 国产精华一区二区三区| 下体分泌物呈黄色| 久久中文字幕一级| 国产成人免费无遮挡视频| 国产真人三级小视频在线观看| 亚洲熟女精品中文字幕| 一边摸一边抽搐一进一出视频| 亚洲成人国产一区在线观看| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费 | 国产亚洲精品第一综合不卡| 欧美精品一区二区免费开放| 欧美乱色亚洲激情| 十八禁网站免费在线| 久久草成人影院| 黄色 视频免费看| 久久精品国产亚洲av高清一级| 最近最新中文字幕大全电影3 | 久久久久精品人妻al黑| www.熟女人妻精品国产| 国产野战对白在线观看| 一级作爱视频免费观看| 女同久久另类99精品国产91| 久久久久久久精品吃奶| 黑人操中国人逼视频| 777米奇影视久久| 欧美人与性动交α欧美软件| 国产精品.久久久| 国产精品 欧美亚洲| 亚洲 国产 在线| 一边摸一边抽搐一进一小说 | 免费看十八禁软件| 天天躁日日躁夜夜躁夜夜| 免费在线观看完整版高清| 最新在线观看一区二区三区| ponron亚洲| 中亚洲国语对白在线视频| 91在线观看av| 国产av一区二区精品久久| 这个男人来自地球电影免费观看| 国产精品一区二区免费欧美| 在线十欧美十亚洲十日本专区| 麻豆乱淫一区二区| 亚洲精品国产区一区二| 午夜福利欧美成人| 久久99一区二区三区| 美女 人体艺术 gogo| 99久久综合精品五月天人人| 韩国av一区二区三区四区| 亚洲一区中文字幕在线| 国产在线观看jvid| 香蕉丝袜av| 最新在线观看一区二区三区| 99国产精品一区二区三区| tube8黄色片| av网站免费在线观看视频| 夫妻午夜视频| 欧美日韩黄片免| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 婷婷精品国产亚洲av在线 | 丝袜美足系列| 嫩草影视91久久| 视频区欧美日本亚洲| 精品久久久久久电影网| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久,| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av片天天在线观看| 老熟女久久久| 极品教师在线免费播放| 国产精品欧美亚洲77777| 麻豆乱淫一区二区| 天天操日日干夜夜撸| 午夜91福利影院| 亚洲成av片中文字幕在线观看| 欧美日韩亚洲国产一区二区在线观看 | 99在线人妻在线中文字幕 | 成人免费观看视频高清| 黄色a级毛片大全视频| 午夜福利乱码中文字幕| 亚洲第一欧美日韩一区二区三区| 亚洲成人免费av在线播放| 日韩有码中文字幕| 欧美亚洲日本最大视频资源| 超色免费av| 欧美日韩av久久| 少妇裸体淫交视频免费看高清 | 99国产精品一区二区三区| 亚洲av日韩精品久久久久久密| 国产高清国产精品国产三级| 一区二区日韩欧美中文字幕| 亚洲国产中文字幕在线视频| avwww免费| 日本撒尿小便嘘嘘汇集6| 成人av一区二区三区在线看| 最新美女视频免费是黄的| 在线观看一区二区三区激情| 久久精品国产99精品国产亚洲性色 | 99久久精品国产亚洲精品| 精品国产超薄肉色丝袜足j| 日韩人妻精品一区2区三区| 超碰成人久久| 丝袜美腿诱惑在线| 在线观看舔阴道视频| 日韩三级视频一区二区三区| 欧美最黄视频在线播放免费 | 欧美日韩中文字幕国产精品一区二区三区 | 成人免费观看视频高清| 91精品国产国语对白视频| 一本一本久久a久久精品综合妖精| 99久久综合精品五月天人人| 亚洲三区欧美一区| 国产精品一区二区在线观看99| 国产成人免费观看mmmm| 不卡av一区二区三区| 在线永久观看黄色视频| 国产av又大| 欧美乱色亚洲激情| 亚洲精品一二三| 免费在线观看亚洲国产| 美女国产高潮福利片在线看| 免费在线观看完整版高清| 精品久久久久久电影网| 国产乱人伦免费视频| 久久久水蜜桃国产精品网| 久久精品国产99精品国产亚洲性色 | 欧美在线黄色| 久9热在线精品视频| 高清欧美精品videossex| 中文字幕人妻丝袜一区二区| 国产日韩一区二区三区精品不卡| 日韩欧美一区二区三区在线观看 | 国产亚洲精品第一综合不卡| 久久久久久亚洲精品国产蜜桃av| 午夜福利视频在线观看免费| 久久草成人影院| 色在线成人网| 国产单亲对白刺激| 国产一区在线观看成人免费| 美女扒开内裤让男人捅视频| 久久精品国产亚洲av香蕉五月 | 国产一区二区三区综合在线观看| 50天的宝宝边吃奶边哭怎么回事| 十八禁网站免费在线| 黄片播放在线免费| 天堂√8在线中文| 在线播放国产精品三级| 午夜免费成人在线视频| 精品久久久久久久毛片微露脸| 两人在一起打扑克的视频| 国产乱人伦免费视频| 国内毛片毛片毛片毛片毛片| 久久午夜综合久久蜜桃| ponron亚洲| 久久99一区二区三区| 在线观看免费高清a一片| 日韩欧美一区二区三区在线观看 | 久久青草综合色| 亚洲第一青青草原| 久久久久视频综合| av线在线观看网站| 国产成人影院久久av| 欧美精品啪啪一区二区三区| 国产一区二区激情短视频| 久久影院123| videos熟女内射| 黄片大片在线免费观看| 午夜激情av网站| 亚洲av电影在线进入| 国产aⅴ精品一区二区三区波| 日韩制服丝袜自拍偷拍| 亚洲人成77777在线视频| 中文字幕最新亚洲高清| 国产真人三级小视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲国产精品合色在线| 亚洲熟妇熟女久久| 成熟少妇高潮喷水视频| 国产欧美日韩一区二区三区在线| 欧美在线一区亚洲| 免费一级毛片在线播放高清视频 | 在线观看免费高清a一片| 人妻 亚洲 视频| 色婷婷av一区二区三区视频| 精品少妇一区二区三区视频日本电影| 久久国产精品大桥未久av| 精品第一国产精品| 亚洲国产精品sss在线观看 | 狠狠狠狠99中文字幕| 日日摸夜夜添夜夜添小说| 80岁老熟妇乱子伦牲交| 日本vs欧美在线观看视频| 身体一侧抽搐| av网站免费在线观看视频| 一级毛片高清免费大全| 亚洲九九香蕉| 手机成人av网站| 精品福利永久在线观看| 黑人巨大精品欧美一区二区mp4| 18在线观看网站| 在线观看www视频免费| 麻豆国产av国片精品| 亚洲国产欧美日韩在线播放| 三级毛片av免费| 日韩欧美国产一区二区入口| 在线av久久热| 亚洲精品久久成人aⅴ小说| 久久久久精品国产欧美久久久| 精品人妻在线不人妻| 中文字幕最新亚洲高清| 国产精品免费一区二区三区在线 | tube8黄色片| av网站在线播放免费| 国产蜜桃级精品一区二区三区 | 成人亚洲精品一区在线观看| av电影中文网址| 亚洲av日韩在线播放| av免费在线观看网站| 国产黄色免费在线视频| 欧美 日韩 精品 国产| 久久久精品国产亚洲av高清涩受| 在线观看一区二区三区激情| 伦理电影免费视频| 国产精品自产拍在线观看55亚洲 | 国产成人精品在线电影| av不卡在线播放| 精品福利永久在线观看| 日本黄色视频三级网站网址 | 搡老岳熟女国产| 欧美中文综合在线视频| 中文字幕人妻丝袜一区二区| 精品卡一卡二卡四卡免费| 男女之事视频高清在线观看| 变态另类成人亚洲欧美熟女 | 亚洲五月天丁香| 麻豆乱淫一区二区| 欧美在线一区亚洲| 亚洲专区字幕在线| 亚洲av美国av| 变态另类成人亚洲欧美熟女 | 国产精品.久久久| 欧美大码av| avwww免费| 女人被狂操c到高潮| 国产精品久久久久久精品古装| 中国美女看黄片| 啦啦啦视频在线资源免费观看| 亚洲视频免费观看视频| 精品午夜福利视频在线观看一区| av天堂在线播放| 免费观看a级毛片全部| 女同久久另类99精品国产91| 国产成人精品无人区| 久久精品国产综合久久久| 国产欧美亚洲国产| 亚洲情色 制服丝袜| 两个人免费观看高清视频| 91麻豆av在线| 自线自在国产av| 国产成人一区二区三区免费视频网站| 啦啦啦免费观看视频1| 欧美大码av| 国产亚洲欧美在线一区二区| www日本在线高清视频| 99riav亚洲国产免费| 免费在线观看影片大全网站| 首页视频小说图片口味搜索| 国产不卡av网站在线观看| 日日夜夜操网爽| 中文字幕人妻丝袜制服| 国产亚洲av高清不卡| 久99久视频精品免费| 欧美精品高潮呻吟av久久| 中文字幕最新亚洲高清| 色94色欧美一区二区| 国产成+人综合+亚洲专区| 男女下面插进去视频免费观看| 一区二区三区激情视频| 热99re8久久精品国产| 国产在线观看jvid| 黄片播放在线免费| 久久中文字幕人妻熟女| av免费在线观看网站| 久久久国产成人免费| 91成人精品电影| 在线播放国产精品三级| 久久精品人人爽人人爽视色| 亚洲色图av天堂| 男人操女人黄网站| 国产在线观看jvid| 欧美精品高潮呻吟av久久| 亚洲avbb在线观看| 老汉色av国产亚洲站长工具| 不卡av一区二区三区| 在线观看免费日韩欧美大片| 新久久久久国产一级毛片| www.熟女人妻精品国产| 午夜福利免费观看在线| 美女国产高潮福利片在线看| 午夜精品久久久久久毛片777| 欧美av亚洲av综合av国产av| 国产精品一区二区免费欧美| 久久人妻福利社区极品人妻图片| 亚洲熟女毛片儿| 性少妇av在线| 午夜免费观看网址| 午夜91福利影院| 1024视频免费在线观看| 少妇粗大呻吟视频| 少妇的丰满在线观看| 99热只有精品国产| 韩国av一区二区三区四区| 视频在线观看一区二区三区| av线在线观看网站| 亚洲精品国产色婷婷电影| 成人亚洲精品一区在线观看| 国产精品1区2区在线观看. | 亚洲熟女毛片儿| 日本vs欧美在线观看视频| 热99久久久久精品小说推荐| 91在线观看av| 桃红色精品国产亚洲av| av中文乱码字幕在线| 夫妻午夜视频| 国产成人精品在线电影| 黑人巨大精品欧美一区二区mp4| 国产麻豆69| 他把我摸到了高潮在线观看| 精品少妇久久久久久888优播| 国产一区二区三区综合在线观看| 一级毛片精品| av超薄肉色丝袜交足视频| 伊人久久大香线蕉亚洲五| 国产精品欧美亚洲77777| 久久精品人人爽人人爽视色| 一级黄色大片毛片| 国产亚洲欧美98| 欧美激情 高清一区二区三区| 亚洲一区二区三区不卡视频| 欧美日韩成人在线一区二区| 在线视频色国产色| 精品国内亚洲2022精品成人 | 午夜福利影视在线免费观看| 成年人黄色毛片网站| 在线观看一区二区三区激情| 国产99久久九九免费精品| а√天堂www在线а√下载 | 中出人妻视频一区二区| 欧美日本中文国产一区发布| 午夜激情av网站| x7x7x7水蜜桃| 日韩欧美三级三区| 高清欧美精品videossex| 欧美日韩视频精品一区| 久久精品熟女亚洲av麻豆精品| 国精品久久久久久国模美| 人成视频在线观看免费观看| 久久人人97超碰香蕉20202| 国产片内射在线| 欧美日韩亚洲综合一区二区三区_| 欧美成人午夜精品| 国产欧美日韩综合在线一区二区| av天堂在线播放| 高清av免费在线| 99国产精品99久久久久| 美女高潮喷水抽搐中文字幕| 亚洲欧洲精品一区二区精品久久久| 在线免费观看的www视频| 亚洲 欧美一区二区三区| 亚洲中文av在线| 午夜福利影视在线免费观看| a级毛片在线看网站| 免费在线观看影片大全网站| 国产成人欧美在线观看 | 亚洲人成电影免费在线| 亚洲专区中文字幕在线| 无人区码免费观看不卡| 国产一区在线观看成人免费| 伦理电影免费视频| 99国产精品免费福利视频| 亚洲专区字幕在线| 免费看十八禁软件| 亚洲国产欧美日韩在线播放| 欧美丝袜亚洲另类 | 欧美av亚洲av综合av国产av| 十分钟在线观看高清视频www| 国产淫语在线视频| 91老司机精品| 午夜福利在线免费观看网站| 9热在线视频观看99| 色精品久久人妻99蜜桃| 在线观看66精品国产| 成人黄色视频免费在线看| 大型黄色视频在线免费观看| 两人在一起打扑克的视频| 亚洲国产看品久久| 精品第一国产精品| 久久精品人人爽人人爽视色| 欧美性长视频在线观看| a在线观看视频网站| 亚洲专区中文字幕在线| 老司机午夜福利在线观看视频| 热99久久久久精品小说推荐| 窝窝影院91人妻| 91av网站免费观看| 麻豆乱淫一区二区| 国产成人一区二区三区免费视频网站| 国产又爽黄色视频| 亚洲情色 制服丝袜|