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

    秦巴山區(qū)植被覆蓋變化及氣候因子驅(qū)動分析

    2019-06-04 07:31:46陳超男朱連奇李新鴿
    生態(tài)學(xué)報(bào) 2019年9期
    關(guān)鍵詞:秦巴山區(qū)氣候因子海拔

    陳超男,朱連奇,田 莉,李新鴿

    1 河南大學(xué)環(huán)境與規(guī)劃學(xué)院, 開封 475004 2 中國科學(xué)院地理科學(xué)與資源研究所, 北京 100101 >3 千煙洲生態(tài)試驗(yàn)站,生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點(diǎn)實(shí)驗(yàn)室, 北京 100101

    植被作為陸地生態(tài)系統(tǒng)的主體,在物質(zhì)循環(huán)和能量流動中起著重要作用,是聯(lián)結(jié)大氣、水分和土壤關(guān)系的自然紐帶,被認(rèn)為是全球變化研究中的“指示器”[1-3],其動態(tài)變化受土壤、氣候、人類活動等多重因素的共同影響,氣候因子中以溫度和降水的影響最為直接[4-7]。歸一化植被指數(shù)(Normalized Difference Vegetation Index,NDVI)能較好的反映區(qū)域植被覆蓋與植被生長狀況,在全球范圍內(nèi)得到廣泛應(yīng)用[8-11]。

    近百年來,氣候系統(tǒng)正經(jīng)歷以全球變暖為主要特征的顯著變化,尤其是在北半球高海拔、高緯度(如蒙古高原、北極)等自然過渡帶與生境脆弱區(qū),溫度升高顯著高于全球溫度升高的平均水平[12-13],由這種變化所引起的陸地生態(tài)系統(tǒng)效應(yīng)已受到人類的高度重視。山地因在較小的空間具有濃縮的環(huán)境梯度和高度異質(zhì)化的生境,對氣候變化的敏感性僅次于極地,亦成為全球變化背景下研究植被動態(tài)對氣候響應(yīng)的熱點(diǎn)區(qū)域[14-16]。

    秦巴山區(qū)(102°24′—112°40′E和30°43′—35°29′N)西接青藏高原東緣,東臨華北平原,漢江以北為秦嶺,南為大巴山,中部分布有漢中盆地和安康盆地。該區(qū)是中國地理格局的重要樞紐、大尺度東西向的生態(tài)廊道,又是最高級別的生態(tài)安全保障區(qū)、扶貧開發(fā)重點(diǎn)區(qū),以復(fù)雜、多樣、獨(dú)特的自然環(huán)境孕育了豐富的天然動植物資源,對其生態(tài)環(huán)境變化的研究具有重要意義。作為氣候變化的敏感區(qū)和生境脆弱區(qū),秦巴山區(qū)植被覆蓋變化及其驅(qū)動力研究是全球變化研究領(lǐng)域的重點(diǎn)內(nèi)容之一,國內(nèi)學(xué)者在這方面已做了大量工作且取得一定成果[14,17-19]。雒新萍[20]研究秦巴山區(qū)1982—2003年植被年均NDVI與氣候因子的相關(guān)性,認(rèn)為秦巴山區(qū)植被覆蓋與溫度的相關(guān)性高于降水,這與劉憲鋒等[21]的研究結(jié)果相反。崔曉臨等[22-23]分析了2000—2009年秦嶺山地植被覆蓋與溫度的響應(yīng),發(fā)現(xiàn)研究時段內(nèi)秦嶺植被覆蓋整體呈增加趨勢,高海拔植被NDVI對溫度變化更為敏感。鄧晨輝等[24]則指出2000—2015年間秦嶺覆蓋植被變化主要受人類活動影響,其次是降水與溫度。張清雨等[25]對秦嶺山地林地類型時空變化格局進(jìn)行研究,發(fā)現(xiàn)1978—2006年研究區(qū)林地面積增加但破碎化日益嚴(yán)重,經(jīng)濟(jì)發(fā)展、適耕地面積等是林地類型轉(zhuǎn)換的主要驅(qū)動力。然而,目前關(guān)于秦巴山區(qū)植被覆蓋變化及其驅(qū)動力研究主要集中于陜西省境內(nèi)[21-22,26-28],對西秦嶺、伏牛山和大巴山研究較少,且前人對秦巴山區(qū)植被覆蓋變化的主要?dú)夂蝌?qū)動因子爭議較大,大多僅用一種植被數(shù)據(jù)對植被覆蓋變化進(jìn)行探究,全面、準(zhǔn)確評價整個秦巴山區(qū)植被覆蓋時空變化特征及其對氣候因子的響應(yīng)機(jī)制、整合碎片化、局域化資料顯得尤為重要與緊迫。鑒于此,本文使用GIMMS3g(1982—2015年)、SPOT VEG(1998—2015年)、MODIS(2000—2017年)3種不同時空分辨率的NDVI數(shù)據(jù)集,對秦巴山區(qū)1982—2017年植被覆蓋的時空變化特征及其與主要?dú)夂蛞蜃拥南嚓P(guān)性進(jìn)行研究,并以秦嶺山地(陜西境內(nèi))為重點(diǎn),分析了植被覆蓋的海拔梯度差異及其與氣候因子的垂直響應(yīng)模式,以期在全球氣候變化的背景下為山地植被保護(hù)與恢復(fù)提供科學(xué)依據(jù)。

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

    1.1 研究區(qū)概況

    秦巴山區(qū)(圖1)地跨甘、川、陜、渝、鄂和豫等6個省市,東西長約1000 km,南北寬約300 km,總面積約為30萬km2。以秦嶺-大巴山為主體的秦巴山區(qū)是長江-黃河流域的分水嶺,也是我國南北過渡帶的主體區(qū)。過渡帶的地理位置使得該區(qū)植被區(qū)系復(fù)雜,秦嶺北坡主要為暖溫帶落葉闊葉林,南坡主要為混有常綠樹種的落葉闊葉林,大巴山北坡與秦嶺南坡植被相似,南坡則為典型的常綠闊葉林,隨海拔升高,氣候、植被等呈明顯的垂直地帶性分布[29-30]。該地區(qū)年平均氣溫介于12—16℃,年平均降水量介于709—1500 mm。因其地理環(huán)境的高度復(fù)雜性、生物多樣性和氣候敏感性,該區(qū)已成為我國地學(xué)、生物學(xué)研究的天然實(shí)驗(yàn)室。

    圖1 秦巴山區(qū)地理位置及其地形示意圖Fig.1 Location and topography of Qinling-Daba Mountains

    1.2 數(shù)據(jù)來源與預(yù)處理

    為了延長對區(qū)域植被覆蓋變化的研究時間,避免遙感數(shù)據(jù)本身的機(jī)械誤差,本研究采用了三套不同時空分辨率的植被遙感數(shù)據(jù)集。其中GIMMS3gV1.0(1982—2015年)數(shù)據(jù)來源于美國全球檢測與模型組(Global Inventor Modeling and Mapping Studies,GIMMS) 15 d合成的空間分辨率為8 km(https://ecocast.arc.nsds.gov)的NDVI數(shù)據(jù)集;SPOT VEG(1998—2015年)數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心提供的逐月NDVI數(shù)據(jù)集(http://www.resdc.cn),空間分辨率為1 km;MODIS(2000—2017年)數(shù)據(jù)來源于美國國家航空航天局(NASA)16 d合成的空間分辨率為250 m MOD13Q1產(chǎn)品。利用MRT(MODIS Reprojection Tool)工具對MODIS數(shù)據(jù)集進(jìn)行拼接、剪裁等預(yù)處理,提取NDVI影像。為減少云、氣溶膠、視角等影響,GIMMS3g和MODIS數(shù)據(jù)均采用最大值合成法MVC(Maximum Value Composites)[31]得到月NDVI。

    氣象數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn)提供的中國年降水量、年平均氣溫空間插值數(shù)據(jù)集,空間分辨率為1 km。

    DEM數(shù)據(jù)(空間分辨率為30 m)來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http://datamirror.csdb.cn)。

    1.3 研究方法

    采用均值法得到年均NDVI和多年累計(jì)平均NDVI,為了避免遙感數(shù)據(jù)本身帶來的誤差,將年均NDVI進(jìn)行z-score標(biāo)準(zhǔn)化處理,利用線性趨勢分析法探討研究區(qū)植被覆蓋年際變化特征,并基于像元尺度的線性趨勢分析法對多年年均NDVI進(jìn)行處理,以回歸方程斜率表征不同區(qū)域植被覆蓋空間變化特征,斜率為正時,植被呈改善趨勢(NDVI隨時間推移而增加),斜率為負(fù)時,植被呈退化趨勢(NDVI隨時間推移而降低),且斜率絕對值越大,植被覆蓋變化趨勢越明顯;最后采用相關(guān)分析法來探明區(qū)域植被覆蓋變化對氣候因子的響應(yīng)狀況。另外,水熱組合、土壤養(yǎng)分等條件隨海拔變化導(dǎo)致植被覆蓋有很大的差異[32-35],本文以秦嶺山地(陜西境內(nèi))為重點(diǎn)區(qū),參考《秦嶺山地自然地理》[36],根據(jù)秦嶺植被分布的垂直帶性將海拔劃分為6個梯度(<500 m,500—1000 m,1000—1600 m,1600—2700 m,2700—3200 m,>3200 m)以觀察植被覆蓋的垂直變化特征。

    由于本文采用了3種不同時空分辨率的植被遙感數(shù)據(jù)集,根據(jù)研究需要,在分析秦巴山區(qū)植被覆蓋長時間序列(1982—2017年)年際變化特征時采用3種數(shù)據(jù),這樣既能對各種遙感數(shù)據(jù)質(zhì)量互相驗(yàn)證又延長了對區(qū)域植被變化的研究時段;分析植被覆蓋空間分布、空間變化(包括不同海拔帶)特征時采用高分辨率的MODIS數(shù)據(jù);分析植被覆蓋變化對氣候因素的響應(yīng)時采用和氣象數(shù)據(jù)具有相同空間分辨率的SPOT VEG數(shù)據(jù)。

    2 結(jié)果與分析

    2.1 植被NDVI(GIMMS3g 1982—2017、SPOT VEG1998—2015和MODIS2000—2017)年際變化特征

    將秦巴山區(qū)1982—2017年(GIMMS3g、SPOT VEG和MODIS數(shù)據(jù))年均NDVI進(jìn)行標(biāo)準(zhǔn)化處理,用其年際變化趨勢代表研究區(qū)植被覆蓋年際變化趨勢(圖2)。

    圖2 1982—2017年秦巴山區(qū)NDVI變化趨勢Fig.2 The trends of the GIMMS3g 1982—2015, SPOT VEG1998—2015 and MODIS2000—2017 in Qinling-Daba Mountains during 1982—2017CV: 變異系數(shù),Coefficient of Variation;G1982—2015: 1982—2015年GIMMS3g NDVI數(shù)據(jù),GIMMS3g NDVI during 1982—2015;G1982—2000: 1982—2000年GIMMS3g NDVI數(shù)據(jù),GIMMS3g NDVI during 1982—2000;G2000—2015: 2000—2015年GIMMS3g NDVI數(shù)據(jù),GIMMS3g NDVI during 2000—2015;S2000—2015: 2000—2015年SPOT VEG NDVI數(shù)據(jù),SOPT VEG NDVI during 2000—2015;M2000—2017: 2000—2017年MODIS NDVI數(shù)據(jù),MODIS NDVI during 2000—2017

    結(jié)果表明,基于不同分辨率的3種遙感數(shù)據(jù)在重疊年份的年際變化趨勢一致,充分證明了該區(qū)域遙感數(shù)據(jù)的可靠性,適用于分析秦巴山區(qū)植被覆蓋的時空變化格局。2000—2017年間,三種遙感數(shù)據(jù)NDVI均呈現(xiàn)穩(wěn)定增加趨勢,且年際波動表現(xiàn)為較好的一致性,GIMMS3g 2000—2015、SPOT VEG2000—2015和MODIS2000—2017增速分別為2.2%/10a(P<0.001)、4.1%/10a(P=0.001)和4.5%/10a(P<0.001),由于傳感器自身差異,GIMMS3g數(shù)據(jù)增速略小于SPOT VEG和MODIS數(shù)據(jù)增速,但整體上可以準(zhǔn)確顯示研究區(qū)植被覆蓋變化趨勢。由GIMMS3g數(shù)據(jù)可知,秦巴山區(qū)1982—2015年間植被覆蓋整體為增加趨勢,增速為1.4%/10a,其中1982—2000年增速為1.6%/10a(P=0.098),植被覆蓋變化主要受氣候因子影響,NDVI水平較低且年際波動較大,僅1990、1991、1995、1997和1998年NDVI高于平均值,其余年份均低于平均水平,在1998—2000年間NDVI迅速下降(與SPOT VEG數(shù)據(jù)表現(xiàn)一致),該時間段為研究區(qū)實(shí)施退耕還林工程初期,耕地面積急劇減少[37],導(dǎo)致植被覆蓋度降低;2000年之后,生態(tài)保護(hù)工程成效顯著,在人為作用的影響下,植被覆蓋呈穩(wěn)定增加趨勢(1.7%/10a,P<0.001),除2005、2011、2012年有較大下降,其余年份NDVI值均處于較高水平且穩(wěn)定上升。計(jì)算GIMMS3g1982—2000NDVI年際變異系數(shù)(CV)為0.04,2000—2015年為0.02,SPOT VEG2000—2015和MODIS2000—2017CV值分別為0.04和0.05,對比不同時期GIMMS3g數(shù)據(jù)的CV值可知,2000年之前,秦巴山區(qū)植被覆蓋變化受氣候因素影響年際波動較大,而2000年之后,在人為因素(如植樹造林工程)影響下穩(wěn)定增加。

    2.2 植被NDVI(MODIS2000—2017)空間分布特征

    由于植被類型、土地利用、氣候因素等不同,秦巴山區(qū)2000—2017年植被覆蓋呈現(xiàn)明顯的區(qū)域差異(圖3a),NDVI介于0.27—0.74之間,平均值為0.59。對比地形圖(圖1)可以發(fā)現(xiàn),研究區(qū)植被覆蓋最好的區(qū)域主要為太白山、米倉山、雪寶頂和神農(nóng)架等山地,NDVI達(dá)到0.60—0.74,這些區(qū)域植被覆蓋以森林為主,植被組成主要有常綠闊葉林、落葉闊葉林以及暖溫帶針闊葉林等;研究區(qū)中部的漢中盆地土地利用類型主要為農(nóng)耕地及人工園林,與林地相比,NDVI有所下降,集中分布于0.4—0.5之間;研究區(qū)東北部海拔較低,人口密集,植被覆蓋受人為干擾嚴(yán)重相對較低,NDVI多分布于0.30—0.40之間,另外,西秦嶺(研究區(qū)西北、西南邊緣的岷縣、武都縣以及西和縣)以及研究區(qū)海拔最高的西南部等地區(qū),因受水熱條件制約,植被類型以亞高山落葉闊葉灌叢、禾草、雜草草甸為主,植被覆蓋度相對較差,NDVI小于0.35。

    圖3 2000—2017年秦巴山區(qū)NDVI空間分布及垂直分布特征Fig.3 Spatial distribution and vertical characteristics of NDVI in Qinling-Daba Mountains during 2000—2017

    由圖3可以看出,秦嶺山地隨著海拔的升高,NDVI值呈先增大后減小的變化特征,海拔<500 m時,NDVI相對較小且分布分散,平均值為0.52,主要介于0.50—0.60之間,占該海拔帶總面積的55.50%,其次為0.40—0.50和0.60—0.70,所占比例分別為23.62%、13.52%,小于0.30和0.30—0.40之間分布較少;500—1000 m,NDVI有所上升且分布較為集中,平均值為0.59,主要介于0.50—0.6與0.6—0.7之間,所占比例分別為41.69%、49.21%,0.30—0.40分布較少;1000—1600 m,NDVI繼續(xù)上升,平均值為0.63,介于0.60—0.70之間的比例為71.13%;1600—2700 m,平均值為0.67,植被覆蓋狀況最好且分布最為集中,NDVI主要介于0.60—0.70之間,所占比例為71.47%,大于0.7的分布比例次之;2700—3200 m,NDVI略有下降,平均值為0.59, 0.6—0.7之間所占比例為39.00%;>3200 m時,NDVI最小,平均值為0.34,小于0.4的區(qū)域占78.11%。

    植被分布受多種因素的影響,其中人為因素與氣候因素最為重要[38-41]。海拔小于500 m時,植被類型主要為農(nóng)耕栽培植被,500—1600 m海拔帶植被覆蓋有所升高,但由于受人類活動影響,不少山坡已被開墾,原始森林植被遭到破壞,植被類型主要為人工林地,植被組成單一[42],且山麓氣候干燥,故植被覆蓋較差;海拔升高至1600—2700 m時,人為干擾較小,森林植被保存完好,水熱條件適宜,植被覆蓋狀況最好,植被組成主要為針闊葉混交林[36];2700—3200 m為亞高山植被帶,植被成分單純,NDVI值迅速下降,海拔高于3200 m時,受溫度限制,植被主要由灌叢和草甸組成,NDVI值最小。

    2.3 植被NDVI(MODIS2000—2017)空間變化特征

    利用線性趨勢分析法對2000—2017年各年均NDVI值與其相對應(yīng)的時間序列進(jìn)行逐項(xiàng)元回歸分析,將其回歸斜率劃分為五類,用以表征該時間段內(nèi)植被覆蓋變化趨勢,劃分標(biāo)準(zhǔn)及結(jié)果見圖4、表1。結(jié)果表明,18 a間秦巴山區(qū)植被覆蓋改善的區(qū)域面積約20.07萬km2,占總面積的96.90%,退化的區(qū)域面積約0.93萬km2,所占比例為3.10%,即研究區(qū)植被覆蓋整體呈增加趨勢。其中輕度改善的面積最大,為13.79萬km2(45.97%),中度改善的面積次之,為13.42萬km2(44.73%),明顯改善的區(qū)域面積為1.86萬km2(6.19%),之后為輕度退化和中度退化,面積分別為0.71萬km2(2.35%)和0.14萬km2(0.47%),明顯退化的面積最小,為0.09萬km2(2.9%)。從空間分布來看,植被明顯改善的區(qū)域主要集中在西秦嶺(甘南的禮縣、武都縣、西和縣),另外,安康中部和十堰北部等低海拔區(qū)域植被覆蓋也以明顯改善為主;中度改善的區(qū)域較為零散,遍布整個研究區(qū);輕度改善的區(qū)域主要集中在研究區(qū)北部高海拔區(qū)域,太白山、米倉山和神農(nóng)架等山地也多表現(xiàn)為輕度改善;植被退化的面積較少,主要集中在研究區(qū)西南部的少數(shù)地區(qū),如成都北部和德陽北部的交界處,漢中盆地及南陽西部邊緣區(qū)域也有零星分布。由上述分析可得,近18 a秦巴山區(qū)植被覆蓋整體呈改善趨勢,以輕度改善和中度改善為主體,變化斜率較大的區(qū)域主要集中在海拔相對較低、人類活動較為頻繁的地區(qū),高海拔區(qū)域主要受氣候因素影響,植被覆蓋變化主要表現(xiàn)為基本不變(輕度退化、輕度改善)。

    圖4 2000—2017年秦巴山區(qū)NDVI空間及垂直變化特征Fig.4 Spatial pattern (a) and vertical characteristics (b) of NDVI change trends in Qinling-Daba Mountains during 2000—2017*和**分別為通過置信水平為95%和99%的顯著性檢驗(yàn)

    類型 Type分類標(biāo)準(zhǔn) Classification protocol面積/萬km2Area面積比例/%Area percent明顯退化 Significant degradationSlope < -0.0020.09 0.29中度退化 Moderate degradation-0.002 < Slope < -0.0010.14 0.47輕度退化 Mild degradation-0.001< Slope < 00.71 2.35輕度改善 Mild improvement0 < Slope < 0.00113.79 45.97中度改善 Moderate improvement0.001< Slope < 0.00213.42 44.73明顯改善 Significant improvementSlope > 0.0021.86 6.19

    為進(jìn)一步探究海拔對植被覆蓋變化的影響,統(tǒng)計(jì)秦嶺地區(qū)不同海拔帶多年來植被覆蓋的變化情況,結(jié)果顯示,不同海拔帶植被覆蓋均呈增加趨勢,但增速有所不同(圖4)。海拔<500 m、500—1000 m、1000—1600 m時,植被類型以農(nóng)作物與人工園林為主,植被覆蓋變化主要受人類活動(耕作制度、灌溉及施肥方式、退耕還林政策等)影響,NDVI增速較大,分別為4.8%/10a、5.7%/10a、5.1%/10a,且均通過了置信水平為99%的顯著性檢驗(yàn);海拔在1600—2700 m與2700—3200 m時,植被覆蓋變化受人為干擾較小,主要受溫度和降水等氣候因子影響,在全球變暖的大背景下,該區(qū)植被覆蓋有所增加,但與低海拔相比增速減小,分別為3.2%/10a(P<0.001)、2.6%/10a(P<0.05);隨著海拔的繼續(xù)升高(>3200 m),NDVI變化不顯著,增速為1.2%/10a(P>0.05)。進(jìn)一步說明人為因素對植被覆蓋影響的速度大于氣候因素對植被覆蓋影響的速度,人類活動可在短時間內(nèi)改變區(qū)域植被覆蓋格局。

    2.4 植被NDVI(SPOT VEG1982—2015)變化的氣候因子驅(qū)動力分析

    溫度和降水是影響植被覆蓋最直接、重要的氣候因子[7],在全球氣候變化的大背景下,1982—2015年秦巴山區(qū)以增溫效應(yīng)為主,溫度增加速率為0.49℃/10a(R=0.826,P<0.001,n=34),增濕效應(yīng)不顯著。在這種情況下,本文以年為單位,逐項(xiàng)元分別計(jì)算2000—2015年植被NDVI與年均溫和年降水的相關(guān)系數(shù),得到其空間分布圖,用以分析研究區(qū)植被覆蓋對氣候因子的響應(yīng)模式,并對相關(guān)系數(shù)進(jìn)行統(tǒng)計(jì)分類(表2)。結(jié)果顯示,研究區(qū)植被覆蓋變化與溫度主要為正相關(guān),與降水為正負(fù)相關(guān)并存,與溫度的相關(guān)性較降水高。

    圖5 秦巴山區(qū)氣溫距平、降水距平百分率變化趨勢,NDVI與溫度、降水與的相關(guān)系數(shù)分布以及對溫度、降水的垂直響應(yīng)模式 Fig.5 The change trends of temperature anomaly, precipitation anomaly percentage; the spatial distribution of correlation among NDVI with temperature and precipitation; the correlation among NDVI and temperature and precipitation at different altitudes in Qinling-Daba Mountains during 2000—2015*為通過置信水平為90%的顯著性檢驗(yàn)

    2000—2015年秦巴山區(qū)植被覆蓋變化與溫度的平均相關(guān)系數(shù)為0.22,區(qū)域正相關(guān)系數(shù)和負(fù)相關(guān)系數(shù)所占比例分別為88.81%和11.19%(圖5,表2),即在研究時段內(nèi),研究區(qū)植被覆蓋變化與溫度主要呈正相關(guān)。在正相關(guān)區(qū)域中,低度正相關(guān)的區(qū)域占總面積的57.12%,分布于整個研究區(qū);中度正相關(guān)和高度正相關(guān)比例分別為10.61%、0.58%,主要分布在太白山、米倉山和神農(nóng)架等森林面積較大的高海拔區(qū)域,這些地區(qū)降水充沛但溫度較低,熱量為該區(qū)域植被生長的主要限制性因子,溫度升高有利于植被覆蓋的增加。負(fù)相關(guān)的區(qū)域面積較小,主要分布在四川北部的松藩縣、平武縣和茂縣等地區(qū),另外,豫西地區(qū)人口密集、人為干擾較大的洛寧縣、欒川縣、內(nèi)鄉(xiāng)縣一帶也與溫度呈負(fù)相關(guān),但相關(guān)性均不顯著。

    相比之下,植被覆蓋變化與降水的正負(fù)相關(guān)性較為平衡,平均相關(guān)系數(shù)為0.03,正相關(guān)與負(fù)相關(guān)所占比例分別為41.53%和58.47%。其中呈低度正相關(guān)的區(qū)域比例為49.14%,主要分布在研究區(qū)北部和南部邊緣地區(qū),中度和高度正相關(guān)的區(qū)域比例分別為7.86%和1.66%,主要位于研究區(qū)西北部的臨潭縣、岷縣和溫縣等降雨較少的地區(qū),這些地區(qū)溫度較高,土壤水分蒸發(fā)嚴(yán)重,降水增加將促進(jìn)植被生長;呈低度負(fù)相關(guān)和中度負(fù)相關(guān)的區(qū)域比例為34.98%和6.55%,主要集中在研究區(qū)東部城鎮(zhèn)化水平較高的大面積低海拔區(qū)域,另外,西南部的高海拔區(qū)域也呈負(fù)相關(guān),這些區(qū)域降水量相對較多,植物生長所需的水分足夠,降水增加將導(dǎo)致光照強(qiáng)度下降、光合作用減弱,從而抑制植被的生長[43]。

    表2 2000—2015年秦巴山區(qū)NDVI與溫度、降水相關(guān)系數(shù)分類標(biāo)準(zhǔn)及結(jié)果統(tǒng)計(jì)

    Table 2 The classification standard for correlation coefficient among NDVI, temperature and precipitation statistical result in Qinling-Daba Mountains during 2000—2015

    相關(guān)性Correlation相關(guān)系數(shù)Correlation coefficient面積比例/% Area percent 溫度 Temperature降水 Precipitation中度負(fù)相關(guān) Moderate negative correlation< -0.30.586.55低度負(fù)相關(guān) Mild negative correlation-0.3—0.010.6134.98低度正相關(guān) Mild positive correlation0.0—0.357.1249.14中度正相關(guān)Moderate positive correlation0.3—0.526.627.68高度正相關(guān) Significant positive correlation> 0.55.071.66

    不同海拔帶植被覆蓋變化對氣候因子的響應(yīng)具有明顯的差異性(圖5)。結(jié)果表明,隨著海拔的升高,植被覆蓋與溫度的相關(guān)性出現(xiàn)由負(fù)到正的轉(zhuǎn)變。海拔小于1600 m時,植被NDVI與溫度主要呈負(fù)相關(guān),但均未通過顯著性檢驗(yàn),1000—1600 m海拔帶與溫度相關(guān)系數(shù)最大;海拔在1600—3200 m時,相關(guān)性轉(zhuǎn)為正相關(guān),2700—3200 m時,相關(guān)系數(shù)略有上升,>3200 m時,相關(guān)系數(shù)最大為0.35,且通過了置信水平為90%的顯著性檢驗(yàn),即高海拔區(qū)域植被對溫度變化更為敏感;與溫度相反,隨著海拔的升高,NDVI與降水的相關(guān)性出現(xiàn)由正到負(fù)的轉(zhuǎn)變,海拔在低于2700 m時兩者為負(fù)相關(guān),高于2700 m時兩者為正相關(guān),低海拔區(qū)域相關(guān)系數(shù)更大,但不同海拔帶相關(guān)性均未通過顯著性檢驗(yàn)。造成這種現(xiàn)象的原因主要為:在低海拔區(qū)域,植被覆蓋類型主要為農(nóng)作物和人工園林,植被覆蓋與人類活動的相關(guān)性高于與氣候因子的相關(guān)性,且栽培植被對水分要求較高,因此與降水呈正相關(guān),隨著海拔的升高,降水增加而溫度降低,植被生長主要受熱量條件限制,因此與溫度呈正相關(guān)。

    3 結(jié)論

    (1)時間變化上,1982—2017年3種NDVI遙感數(shù)據(jù)在秦巴山區(qū)均呈顯著上升趨勢。1982—2015年NDVI-GIMMS3g增速為1.4%/10a,其中1982—2000年NDVI-GIMMS3g整體較低且波動幅度較大,增速為1.6%/10a,年際間變異系數(shù)(CV)為0.04,2000年后 NDVI較高且呈穩(wěn)定增加狀態(tài),增速為1.7%/10a,CV為0.02;2000—2015年NDVI-SPOT VEG增速為4.1%/10a,CV為0.04;2000—2017年NDVI-MODIS增速為4.5%/10a,CV為0.05。

    (2)空間分布上,2000—2017年(MODIS數(shù)據(jù))秦巴山區(qū)高植被覆蓋區(qū)主要為太白山、米倉山、神農(nóng)架等山地,低植被覆蓋區(qū)主要為研究區(qū)西部的高海拔區(qū)和東部受人為干擾嚴(yán)重的低海拔區(qū);隨海拔升高,秦嶺山地NDVI呈先增加后減少的變化趨勢。

    (3)空間變化上,2000—2017年(MODIS數(shù)據(jù))秦巴山區(qū)植被覆蓋增加的區(qū)域大于減少的區(qū)域,增加和減少的區(qū)域分別占總面積的96.90%和3.10%,低海拔區(qū)域變化顯著,高海拔區(qū)域基本不變;秦嶺山地不同海拔梯度帶植被覆蓋均呈增加趨勢,低海拔增速較高海拔大。

    (4)1982—2015年(SPOT VEG數(shù)據(jù))秦巴山區(qū)增溫效應(yīng)顯著,增速為0.49℃/10a,增濕效應(yīng)不顯著;2000—2015年研究區(qū)植被覆蓋與溫度以正相關(guān)為主,主要分布在山地等降水較多的地區(qū);與降水為正負(fù)相關(guān)性并存,正相關(guān)主要分布在研究區(qū)西北部降水較少的半干旱區(qū),負(fù)相關(guān)主要分布在東部城鎮(zhèn)化水平較高的低海拔區(qū)域;秦嶺山地隨著海拔的升高,植被覆蓋與溫度的相關(guān)系數(shù)呈由負(fù)到正的轉(zhuǎn)變,高海拔區(qū)域?qū)囟雀舾?與降水的相關(guān)系數(shù)呈由正到負(fù)的轉(zhuǎn)變,低海拔區(qū)對降水更敏感。

    研究結(jié)果表明,1982—2017年秦巴山區(qū)植被覆蓋呈顯著增加趨勢,1982—2000年植被NDVI較低且年際波動大,2000年后植被NDVI較高且穩(wěn)定增加。植被動態(tài)變化受氣候、人為和歷史等多種因素共同影響,在研究時段內(nèi),研究區(qū)溫度和降水均呈增加趨勢,增溫效應(yīng)尤為明顯,但對2000—2015年植被覆蓋與氣溫和降水進(jìn)行相關(guān)分析,發(fā)現(xiàn)植被覆蓋與溫度以正相關(guān)為主,與降水為正負(fù)相關(guān)并存,且高海拔區(qū)域植被對溫度變化更為敏感,低海拔區(qū)域?qū)邓鼮槊舾?這與以往的研究結(jié)果相似[13,20,22,27,36],但無論是氣溫還是降水,呈明顯相關(guān)的區(qū)域面積并不是很大[44],造成這種現(xiàn)象的可能性主要有以下幾點(diǎn):首先,氣候雖為植被覆蓋變化的決定性因子,但兩者的相關(guān)性受研究時段和空間尺度影響[5],研究的時空尺度不同,相關(guān)系數(shù)將存在較大的差異性;其次,植被覆蓋變化對氣候的響應(yīng)模式存在顯著的空間差異性,并可能存在一定的滯后性[14,45-47],本文僅研究了年均溫、年降水量對植被NDVI的影響,未對其滯后性進(jìn)行分析;另外,有研究表明,近年來植被覆蓋變化的主要驅(qū)動力已由氣候因素轉(zhuǎn)為人為因素[48],人類活動能在短時間內(nèi)改變區(qū)域環(huán)境,退耕還林、植樹造林等生態(tài)保護(hù)工程可快速增加區(qū)域植被覆蓋[21,49-50],自2000年后,研究區(qū)全面實(shí)施退耕還林政策,林地增加面積為2.54×106hm2, 約占全國退耕還林總面積的10%[37],植被覆蓋變化受人為因素干擾較大。因此,秦巴山區(qū)植被覆蓋變化對不同時段氣候因子的響應(yīng)以及定量分析人為因素(社會經(jīng)濟(jì)、城市擴(kuò)張等)對植被覆蓋的影響是以后工作的重點(diǎn)內(nèi)容,以期更為全面的揭示秦巴山區(qū)植被覆蓋對氣候因素與人為因素的響應(yīng)機(jī)制,為山地植被保護(hù)和恢復(fù)提供科學(xué)依據(jù)。

    猜你喜歡
    秦巴山區(qū)氣候因子海拔
    青海省果洛州高原鼠兔宜生區(qū)分布及其對氣候變化的響應(yīng)
    高海拔的甘南
    散文詩(2021年22期)2022-01-12 06:14:04
    南美三大高海拔球場
    海拔4000米的父愛
    秦巴山區(qū)牟牛溝滑坡體治理施工技術(shù)
    近16年秦巴山區(qū)TRMM降水資料的降尺度研究
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    我校協(xié)辦“放飛夢想的翅膀·秦巴山區(qū)公益行動”
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    国产精品日韩av在线免费观看| 日韩国内少妇激情av| 国产高清三级在线| 干丝袜人妻中文字幕| 久久久久性生活片| 久久99热6这里只有精品| 青春草国产在线视频| 不卡视频在线观看欧美| 91狼人影院| 日日摸夜夜添夜夜添av毛片| 免费无遮挡裸体视频| 久久精品91蜜桃| 亚洲国产高清在线一区二区三| 插逼视频在线观看| 成年女人看的毛片在线观看| 亚洲av日韩在线播放| 国产午夜精品久久久久久一区二区三区| 日本午夜av视频| 热99在线观看视频| 日韩中字成人| 99热网站在线观看| 五月玫瑰六月丁香| av又黄又爽大尺度在线免费看 | 一级毛片久久久久久久久女| 日本爱情动作片www.在线观看| 六月丁香七月| 久久久久九九精品影院| 午夜激情欧美在线| 老司机福利观看| 国产高清三级在线| 日韩国内少妇激情av| 天天躁夜夜躁狠狠久久av| 国产精品一区二区三区四区久久| 午夜爱爱视频在线播放| av专区在线播放| 国产免费一级a男人的天堂| 亚洲av.av天堂| 最后的刺客免费高清国语| 国产91av在线免费观看| 美女被艹到高潮喷水动态| 少妇熟女aⅴ在线视频| 精品国产露脸久久av麻豆 | 男人舔女人下体高潮全视频| 69av精品久久久久久| 久热久热在线精品观看| 老司机影院毛片| 亚洲欧美成人综合另类久久久 | 天堂中文最新版在线下载 | 我的老师免费观看完整版| 97超视频在线观看视频| 乱系列少妇在线播放| 久久久久久久亚洲中文字幕| 看免费成人av毛片| 欧美区成人在线视频| 日本-黄色视频高清免费观看| 亚洲av中文av极速乱| 国产精品三级大全| 欧美丝袜亚洲另类| 国产精品不卡视频一区二区| 欧美日韩在线观看h| 亚洲av一区综合| 色视频www国产| 91狼人影院| av播播在线观看一区| 色哟哟·www| 亚洲成人久久爱视频| 亚洲av中文字字幕乱码综合| 国产午夜精品久久久久久一区二区三区| 免费大片18禁| 国产精品久久电影中文字幕| av.在线天堂| 免费不卡的大黄色大毛片视频在线观看 | 91久久精品国产一区二区三区| 国产美女午夜福利| 亚洲,欧美,日韩| 亚洲欧美日韩无卡精品| 精品免费久久久久久久清纯| 中文字幕制服av| 女人久久www免费人成看片 | 欧美人与善性xxx| 亚洲国产精品sss在线观看| 日本免费a在线| 国产一级毛片在线| 老司机影院毛片| 最近最新中文字幕免费大全7| 美女cb高潮喷水在线观看| 真实男女啪啪啪动态图| 男人舔奶头视频| 久久久精品欧美日韩精品| 国产乱来视频区| av在线亚洲专区| 亚洲国产精品sss在线观看| 亚洲av熟女| 日本一二三区视频观看| 亚洲最大成人手机在线| 午夜激情福利司机影院| 精品人妻一区二区三区麻豆| a级一级毛片免费在线观看| 国产伦在线观看视频一区| 九九在线视频观看精品| 日韩精品青青久久久久久| 日本wwww免费看| 一级毛片电影观看 | 国产91av在线免费观看| 啦啦啦韩国在线观看视频| 九草在线视频观看| 91午夜精品亚洲一区二区三区| 七月丁香在线播放| 一区二区三区四区激情视频| 最近2019中文字幕mv第一页| 久久人妻av系列| 网址你懂的国产日韩在线| 国产精品福利在线免费观看| 最近最新中文字幕免费大全7| 国产精品女同一区二区软件| 性色avwww在线观看| 日韩一区二区三区影片| 99久久成人亚洲精品观看| 中文字幕av在线有码专区| 亚洲色图av天堂| 日产精品乱码卡一卡2卡三| 我的女老师完整版在线观看| av免费在线看不卡| 久久久久九九精品影院| 18+在线观看网站| 欧美日本亚洲视频在线播放| 欧美xxxx黑人xx丫x性爽| 真实男女啪啪啪动态图| 久久热精品热| 女人十人毛片免费观看3o分钟| 九九爱精品视频在线观看| 全区人妻精品视频| 我要看日韩黄色一级片| 亚洲av电影不卡..在线观看| 看黄色毛片网站| 久久热精品热| АⅤ资源中文在线天堂| 亚洲综合精品二区| 蜜桃亚洲精品一区二区三区| 啦啦啦观看免费观看视频高清| 日本熟妇午夜| 中文字幕av成人在线电影| 亚洲怡红院男人天堂| 国产高清国产精品国产三级 | 亚洲人成网站在线播| 不卡视频在线观看欧美| 日韩av不卡免费在线播放| 色哟哟·www| 女人被狂操c到高潮| 一个人看的www免费观看视频| 亚洲av福利一区| 在线免费观看不下载黄p国产| 国产精品99久久久久久久久| 久久久久久久久久黄片| 国产单亲对白刺激| 精品一区二区三区人妻视频| 国产精品女同一区二区软件| 久久久精品大字幕| 亚洲成人中文字幕在线播放| 亚洲国产色片| 伦理电影大哥的女人| 日韩大片免费观看网站 | 久久精品91蜜桃| 久久久久久大精品| 久久精品国产鲁丝片午夜精品| 两个人的视频大全免费| 91久久精品国产一区二区成人| 久久人妻av系列| 亚洲精品日韩在线中文字幕| 一级毛片我不卡| 天堂av国产一区二区熟女人妻| 国产 一区精品| 亚洲国产最新在线播放| 亚洲性久久影院| 久久久成人免费电影| 亚洲欧美日韩高清专用| 欧美激情久久久久久爽电影| 国产一区二区亚洲精品在线观看| 国产精品久久电影中文字幕| 午夜福利在线在线| 亚洲欧美清纯卡通| 亚洲国产精品合色在线| 国产国拍精品亚洲av在线观看| 成人鲁丝片一二三区免费| 夫妻性生交免费视频一级片| 成人漫画全彩无遮挡| 99久久人妻综合| 久久热精品热| 久久久久九九精品影院| 尾随美女入室| 欧美精品国产亚洲| 国产在线一区二区三区精 | 国产亚洲精品久久久com| 免费av观看视频| 精品一区二区免费观看| 秋霞在线观看毛片| 天堂中文最新版在线下载 | 亚洲精品日韩av片在线观看| 亚洲综合精品二区| www.av在线官网国产| 亚洲av成人av| 国产精品伦人一区二区| 麻豆久久精品国产亚洲av| 青青草视频在线视频观看| 国产在视频线在精品| 亚洲精品自拍成人| 午夜福利视频1000在线观看| 欧美日本视频| 免费不卡的大黄色大毛片视频在线观看 | 免费看日本二区| 精品熟女少妇av免费看| 久久6这里有精品| 亚洲丝袜综合中文字幕| 色综合亚洲欧美另类图片| 久久人人爽人人片av| 国产高清不卡午夜福利| 久久久久久大精品| 黄色欧美视频在线观看| 婷婷六月久久综合丁香| 男人狂女人下面高潮的视频| 亚洲av中文av极速乱| 成人无遮挡网站| 尾随美女入室| 亚洲国产欧美人成| 观看免费一级毛片| 精品人妻视频免费看| 久久久久网色| 欧美不卡视频在线免费观看| 99九九线精品视频在线观看视频| 99久久精品一区二区三区| 男女下面进入的视频免费午夜| 免费黄色在线免费观看| 亚洲欧美一区二区三区黑人 | 亚洲av综合色区一区| 国产成人精品久久久久久| 久久久久久久精品精品| av线在线观看网站| 纵有疾风起免费观看全集完整版| 午夜91福利影院| 亚洲av在线观看美女高潮| 香蕉丝袜av| 久久人人爽av亚洲精品天堂| 久久这里有精品视频免费| 精品久久久久久电影网| 激情视频va一区二区三区| 日本免费在线观看一区| 精品少妇久久久久久888优播| 亚洲国产看品久久| 最近最新中文字幕免费大全7| 婷婷成人精品国产| av在线观看视频网站免费| 亚洲欧美清纯卡通| 免费高清在线观看日韩| 久久99热这里只频精品6学生| 午夜免费观看性视频| 亚洲国产av新网站| 欧美丝袜亚洲另类| 人成视频在线观看免费观看| 婷婷色综合www| 久久久久精品人妻al黑| 亚洲国产精品一区三区| 另类亚洲欧美激情| 亚洲成国产人片在线观看| 久久久久久久大尺度免费视频| 久久免费观看电影| 在现免费观看毛片| 免费观看av网站的网址| 人妻人人澡人人爽人人| 欧美日本中文国产一区发布| 精品第一国产精品| 蜜桃在线观看..| 视频在线观看一区二区三区| 中文字幕精品免费在线观看视频 | 一本大道久久a久久精品| 国产欧美亚洲国产| 国产综合精华液| 精品视频人人做人人爽| 国产男人的电影天堂91| 成年av动漫网址| 国产极品天堂在线| 一二三四在线观看免费中文在 | 久久97久久精品| 热99久久久久精品小说推荐| 国产男女内射视频| 亚洲欧美色中文字幕在线| 欧美bdsm另类| 婷婷成人精品国产| 日本av免费视频播放| 最近2019中文字幕mv第一页| 欧美精品人与动牲交sv欧美| 在线 av 中文字幕| 大话2 男鬼变身卡| 久久国内精品自在自线图片| 欧美日韩综合久久久久久| 久久久久久久精品精品| 亚洲成av片中文字幕在线观看 | 男男h啪啪无遮挡| 久久亚洲国产成人精品v| 国产精品免费大片| 成年动漫av网址| 桃花免费在线播放| 亚洲伊人色综图| 观看av在线不卡| 一级黄片播放器| 纵有疾风起免费观看全集完整版| 亚洲精品日本国产第一区| 久久国产精品大桥未久av| 日本av手机在线免费观看| 国产不卡av网站在线观看| 日韩三级伦理在线观看| 在线观看一区二区三区激情| 精品午夜福利在线看| 丝瓜视频免费看黄片| 午夜av观看不卡| 韩国av在线不卡| 欧美精品一区二区大全| a 毛片基地| 亚洲国产精品一区二区三区在线| 在现免费观看毛片| 满18在线观看网站| 亚洲精品日本国产第一区| 一级毛片黄色毛片免费观看视频| 一本大道久久a久久精品| 夫妻性生交免费视频一级片| 亚洲国产最新在线播放| 国产熟女欧美一区二区| 性色avwww在线观看| 欧美最新免费一区二区三区| 成人无遮挡网站| 国产免费现黄频在线看| 精品卡一卡二卡四卡免费| 亚洲精品美女久久av网站| 久久人人爽av亚洲精品天堂| 国产探花极品一区二区| 午夜免费鲁丝| 亚洲精品美女久久久久99蜜臀 | 肉色欧美久久久久久久蜜桃| 高清欧美精品videossex| 欧美日韩综合久久久久久| 丝袜美足系列| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 国产欧美亚洲国产| 美女福利国产在线| 色94色欧美一区二区| 卡戴珊不雅视频在线播放| 亚洲精品美女久久久久99蜜臀 | 你懂的网址亚洲精品在线观看| 久久精品国产综合久久久 | 97人妻天天添夜夜摸| 久久99一区二区三区| 亚洲精品国产av成人精品| 插逼视频在线观看| 大香蕉97超碰在线| 亚洲国产精品一区二区三区在线| 欧美日韩av久久| 纵有疾风起免费观看全集完整版| 国产黄色视频一区二区在线观看| 男女边摸边吃奶| 水蜜桃什么品种好| 婷婷色麻豆天堂久久| 最近的中文字幕免费完整| 久久99一区二区三区| 老司机影院成人| 天堂中文最新版在线下载| 夜夜骑夜夜射夜夜干| 国产av精品麻豆| 午夜老司机福利剧场| 大陆偷拍与自拍| 美女内射精品一级片tv| 99精国产麻豆久久婷婷| 热99国产精品久久久久久7| 亚洲欧美精品自产自拍| 国产 精品1| 22中文网久久字幕| 少妇人妻精品综合一区二区| 日产精品乱码卡一卡2卡三| 麻豆乱淫一区二区| 波多野结衣一区麻豆| 丝袜在线中文字幕| 亚洲av成人精品一二三区| 午夜精品国产一区二区电影| 99热网站在线观看| 日韩视频在线欧美| 国产视频首页在线观看| 午夜福利,免费看| 综合色丁香网| 成人免费观看视频高清| 亚洲精品国产av成人精品| 精品熟女少妇av免费看| a 毛片基地| 男女边摸边吃奶| 欧美成人午夜精品| 69精品国产乱码久久久| 国产爽快片一区二区三区| 国产国语露脸激情在线看| 内地一区二区视频在线| 十分钟在线观看高清视频www| 乱码一卡2卡4卡精品| 女人精品久久久久毛片| 夜夜爽夜夜爽视频| 97在线视频观看| 国产精品国产av在线观看| 亚洲欧洲精品一区二区精品久久久 | xxxhd国产人妻xxx| 亚洲经典国产精华液单| 欧美人与性动交α欧美软件 | 另类精品久久| 成人国产麻豆网| 极品少妇高潮喷水抽搐| xxx大片免费视频| 青春草视频在线免费观看| 夫妻午夜视频| 午夜福利,免费看| av在线播放精品| 日韩熟女老妇一区二区性免费视频| 好男人视频免费观看在线| 国产av码专区亚洲av| 国产精品嫩草影院av在线观看| 国产永久视频网站| 超碰97精品在线观看| 国产一区二区在线观看日韩| 国产亚洲av片在线观看秒播厂| 日韩三级伦理在线观看| 国产午夜精品一二区理论片| 久久久久久久精品精品| 99精国产麻豆久久婷婷| a 毛片基地| 免费人妻精品一区二区三区视频| 狠狠精品人妻久久久久久综合| 美女国产视频在线观看| 在线观看免费视频网站a站| 免费av中文字幕在线| a级片在线免费高清观看视频| 看免费av毛片| 国产av码专区亚洲av| a级毛片黄视频| 亚洲精品日韩在线中文字幕| 如何舔出高潮| 色婷婷av一区二区三区视频| 日韩视频在线欧美| 国产日韩一区二区三区精品不卡| 久久鲁丝午夜福利片| 精品国产国语对白av| 国产免费福利视频在线观看| 国产精品国产三级国产专区5o| 日本av免费视频播放| 韩国精品一区二区三区 | kizo精华| 啦啦啦视频在线资源免费观看| 最近中文字幕2019免费版| 在线 av 中文字幕| 国产精品成人在线| 亚洲美女搞黄在线观看| 777米奇影视久久| 少妇的逼水好多| 免费观看性生交大片5| 老司机影院毛片| 2021少妇久久久久久久久久久| av网站免费在线观看视频| 天堂俺去俺来也www色官网| 最新中文字幕久久久久| 午夜影院在线不卡| av在线播放精品| 色网站视频免费| 99久国产av精品国产电影| 日韩一本色道免费dvd| 亚洲精品456在线播放app| 精品久久久久久电影网| 秋霞伦理黄片| 校园人妻丝袜中文字幕| 99精国产麻豆久久婷婷| 97精品久久久久久久久久精品| 国产视频首页在线观看| 久久精品国产亚洲av涩爱| 99久久综合免费| 天堂中文最新版在线下载| 亚洲欧美日韩另类电影网站| 免费播放大片免费观看视频在线观看| 国产免费又黄又爽又色| 亚洲av电影在线观看一区二区三区| av有码第一页| 男人爽女人下面视频在线观看| 91国产中文字幕| av线在线观看网站| 国国产精品蜜臀av免费| 亚洲国产日韩一区二区| 免费日韩欧美在线观看| 大香蕉久久成人网| 乱人伦中国视频| 一级a做视频免费观看| 亚洲精品美女久久久久99蜜臀 | 免费人妻精品一区二区三区视频| 国产在线一区二区三区精| 有码 亚洲区| 久久精品国产亚洲av涩爱| 一区二区三区乱码不卡18| a级毛片在线看网站| 久久久久久人人人人人| 精品久久蜜臀av无| 久久久久久久久久成人| 最近手机中文字幕大全| freevideosex欧美| 亚洲av在线观看美女高潮| 色哟哟·www| 国产欧美亚洲国产| 精品一区二区三区视频在线| av免费观看日本| 亚洲国产毛片av蜜桃av| 黄片播放在线免费| av.在线天堂| 国产亚洲一区二区精品| 日日撸夜夜添| 尾随美女入室| 春色校园在线视频观看| 蜜臀久久99精品久久宅男| 久久精品国产a三级三级三级| 有码 亚洲区| 黄色一级大片看看| 中文乱码字字幕精品一区二区三区| 国产一区二区三区av在线| 国产乱人偷精品视频| 国产亚洲最大av| videosex国产| 少妇 在线观看| 精品一区二区三区四区五区乱码 | 大码成人一级视频| 十八禁网站网址无遮挡| 女性被躁到高潮视频| 香蕉国产在线看| 久久久久久久亚洲中文字幕| 校园人妻丝袜中文字幕| 最近最新中文字幕大全免费视频 | xxx大片免费视频| 91精品三级在线观看| 最近手机中文字幕大全| 欧美精品一区二区免费开放| 国产成人精品一,二区| 青青草视频在线视频观看| 狠狠婷婷综合久久久久久88av| 性色av一级| 最近最新中文字幕免费大全7| 色5月婷婷丁香| 一本大道久久a久久精品| 成人手机av| 在线观看一区二区三区激情| 日本爱情动作片www.在线观看| 久久99一区二区三区| 女的被弄到高潮叫床怎么办| 国产激情久久老熟女| 久久久久久久久久久免费av| 最近的中文字幕免费完整| 极品人妻少妇av视频| 日韩在线高清观看一区二区三区| 亚洲精品国产av成人精品| 两性夫妻黄色片 | 极品人妻少妇av视频| 亚洲欧美成人精品一区二区| 一本大道久久a久久精品| 国产一区二区三区av在线| 国产色婷婷99| 日日爽夜夜爽网站| 久久免费观看电影| 精品人妻一区二区三区麻豆| 国产精品女同一区二区软件| 满18在线观看网站| 熟妇人妻不卡中文字幕| 狂野欧美激情性xxxx在线观看| 在线天堂最新版资源| 999精品在线视频| 日韩av免费高清视频| 九草在线视频观看| 97精品久久久久久久久久精品| 久久久久久伊人网av| 中国美白少妇内射xxxbb| 欧美人与性动交α欧美软件 | 巨乳人妻的诱惑在线观看| 亚洲,一卡二卡三卡| 日韩av不卡免费在线播放| videos熟女内射| 久久久久久久精品精品| 国产永久视频网站| av在线app专区| 日本av手机在线免费观看| 欧美精品人与动牲交sv欧美| 十八禁高潮呻吟视频| 亚洲情色 制服丝袜| 午夜福利网站1000一区二区三区| 韩国精品一区二区三区 | 色网站视频免费| 成人二区视频| 国产精品三级大全| 99国产综合亚洲精品| 亚洲精品456在线播放app| 国国产精品蜜臀av免费| 十八禁网站网址无遮挡| 欧美日韩国产mv在线观看视频| 午夜福利网站1000一区二区三区| 国产日韩欧美视频二区| 日本午夜av视频| 男女午夜视频在线观看 | 亚洲三级黄色毛片| av天堂久久9| 岛国毛片在线播放| 9热在线视频观看99| 精品久久久久久电影网| 欧美人与善性xxx| 水蜜桃什么品种好| 最近最新中文字幕免费大全7| 久久久久国产网址| 欧美成人午夜精品| 午夜激情久久久久久久| 人人妻人人澡人人爽人人夜夜| 韩国高清视频一区二区三区| 国产男女内射视频|