• <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软件| 精品酒店卫生间| 亚洲精品久久久久久婷婷小说| 国产视频首页在线观看| 永久网站在线| 亚洲国产日韩一区二区| 久久精品国产亚洲av涩爱| 美女高潮到喷水免费观看| 日本免费在线观看一区| av在线观看视频网站免费| 日韩不卡一区二区三区视频在线| 亚洲国产av新网站| 欧美精品一区二区大全| 精品一品国产午夜福利视频| 人妻一区二区av| 久久久久精品人妻al黑| 国产在线免费精品| 国产一区二区激情短视频 | 91成人精品电影| 亚洲少妇的诱惑av| 久久精品国产鲁丝片午夜精品| 一边摸一边做爽爽视频免费| 日韩 亚洲 欧美在线| 精品国产超薄肉色丝袜足j| 国产亚洲一区二区精品| 999精品在线视频| av在线app专区| 亚洲国产欧美日韩在线播放| 欧美变态另类bdsm刘玥| 久久精品夜色国产| 一级a爱视频在线免费观看| 热99国产精品久久久久久7| 欧美 亚洲 国产 日韩一| 在线 av 中文字幕| 女性被躁到高潮视频| 最近中文字幕高清免费大全6| 在线观看国产h片| 大片免费播放器 马上看| 在现免费观看毛片| 性色av一级| 女人被躁到高潮嗷嗷叫费观| 97在线人人人人妻| 一个人免费看片子| 九九爱精品视频在线观看| 一区福利在线观看| 久久国产精品男人的天堂亚洲| 国产又色又爽无遮挡免| videosex国产| 丝袜人妻中文字幕| 伦精品一区二区三区| 国产av精品麻豆| 岛国毛片在线播放| av.在线天堂| 美女大奶头黄色视频| 欧美老熟妇乱子伦牲交| tube8黄色片| 免费黄网站久久成人精品| 久久久国产欧美日韩av| 日韩电影二区| 久久这里有精品视频免费| 久久精品国产a三级三级三级| 国产深夜福利视频在线观看| 久久精品国产鲁丝片午夜精品| 国产欧美日韩一区二区三区在线| √禁漫天堂资源中文www| 秋霞在线观看毛片| av又黄又爽大尺度在线免费看| 成年女人毛片免费观看观看9 | 久久久久精品久久久久真实原创| 中文字幕色久视频| av天堂久久9| 久久国产亚洲av麻豆专区| 午夜老司机福利剧场| 天天躁日日躁夜夜躁夜夜| 久久人人爽人人片av| 99久国产av精品国产电影| 亚洲精品国产av成人精品| 欧美精品亚洲一区二区| 在线观看国产h片| 99久久综合免费| 韩国av在线不卡| 飞空精品影院首页| 成年女人在线观看亚洲视频| 99精国产麻豆久久婷婷| 中文天堂在线官网| 最黄视频免费看| 男女啪啪激烈高潮av片| 少妇人妻 视频| 久久久国产一区二区| 一级毛片电影观看| 伦理电影免费视频| 国产1区2区3区精品| 高清欧美精品videossex| 国产乱来视频区| 国产淫语在线视频| 黑丝袜美女国产一区| 亚洲精品自拍成人| 欧美黄色片欧美黄色片| 黑丝袜美女国产一区| 日韩av不卡免费在线播放| 精品国产一区二区久久| 免费观看性生交大片5| 丝袜脚勾引网站| 青春草国产在线视频| 黄色视频在线播放观看不卡| 欧美+日韩+精品| 国产成人精品久久二区二区91 | 男女午夜视频在线观看| 80岁老熟妇乱子伦牲交| 美女大奶头黄色视频| 日本午夜av视频| 亚洲欧洲日产国产| 爱豆传媒免费全集在线观看| 国产激情久久老熟女| 一本大道久久a久久精品| 日本91视频免费播放| 美女高潮到喷水免费观看| 男女啪啪激烈高潮av片| 日韩av免费高清视频| 最近手机中文字幕大全| 免费黄网站久久成人精品| 国产不卡av网站在线观看| 老司机影院毛片| 成人手机av| 久久久久国产一级毛片高清牌| 日本欧美国产在线视频| 亚洲国产精品国产精品| 亚洲,欧美,日韩| 老鸭窝网址在线观看| 国产一区亚洲一区在线观看| 日韩免费高清中文字幕av| 久久久国产欧美日韩av| 国产白丝娇喘喷水9色精品| 国产熟女午夜一区二区三区| 午夜av观看不卡| 少妇人妻久久综合中文| 中文精品一卡2卡3卡4更新| 人人妻人人添人人爽欧美一区卜| 99久国产av精品国产电影| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 黄色配什么色好看| 1024视频免费在线观看| 亚洲一级一片aⅴ在线观看| av电影中文网址| 国产无遮挡羞羞视频在线观看| av在线app专区| 亚洲精品美女久久久久99蜜臀 | 久久久久人妻精品一区果冻| 欧美国产精品一级二级三级| 热re99久久国产66热| 伊人久久大香线蕉亚洲五| 黄色配什么色好看| 搡女人真爽免费视频火全软件| 亚洲 欧美一区二区三区| 在线 av 中文字幕| 涩涩av久久男人的天堂| 亚洲视频免费观看视频| 久久亚洲国产成人精品v| 国产精品国产三级专区第一集| 日本wwww免费看| 色播在线永久视频| a 毛片基地| 国产淫语在线视频| 一级片免费观看大全| 中文字幕色久视频| 午夜91福利影院| 在线亚洲精品国产二区图片欧美| 成年动漫av网址| av.在线天堂| 欧美人与性动交α欧美精品济南到 | 国产一级毛片在线| 黑人巨大精品欧美一区二区蜜桃| 久久人人97超碰香蕉20202| 久久av网站| 熟女电影av网| av福利片在线| 91在线精品国自产拍蜜月| 韩国高清视频一区二区三区| 赤兔流量卡办理| 男女啪啪激烈高潮av片| 一级片'在线观看视频| 国产在视频线精品| 日本免费在线观看一区| 欧美精品亚洲一区二区| 丝瓜视频免费看黄片| 一区福利在线观看| 国产av码专区亚洲av| 一区在线观看完整版| 日韩中字成人| 少妇精品久久久久久久| 国产成人一区二区在线| 成人毛片60女人毛片免费| 日韩中字成人| 亚洲av成人精品一二三区| 国产亚洲欧美精品永久| 夜夜骑夜夜射夜夜干| 99久久人妻综合| 岛国毛片在线播放| 欧美日韩精品成人综合77777| 综合色丁香网| 国产亚洲一区二区精品| 人人妻人人爽人人添夜夜欢视频| 最新的欧美精品一区二区| www.熟女人妻精品国产| 水蜜桃什么品种好| 国产又爽黄色视频| 永久网站在线| 亚洲精品国产一区二区精华液| a级毛片黄视频| 国产午夜精品一二区理论片| 叶爱在线成人免费视频播放| 男女下面插进去视频免费观看| 爱豆传媒免费全集在线观看| 久久久久久人妻| 日韩av不卡免费在线播放| 丝袜美足系列| 久久av网站| 亚洲av.av天堂| 人妻系列 视频| 国产高清不卡午夜福利| 最近中文字幕高清免费大全6| 国产日韩欧美在线精品| 精品一区二区三区四区五区乱码 | 男的添女的下面高潮视频| 在线观看美女被高潮喷水网站| 精品亚洲成国产av| 久久99蜜桃精品久久| 亚洲精品国产av成人精品| 国产不卡av网站在线观看| 亚洲综合色惰| 国产精品久久久久久久久免| 国产极品天堂在线| 中文乱码字字幕精品一区二区三区| 国产精品二区激情视频| 激情五月婷婷亚洲| 久久狼人影院| 一本—道久久a久久精品蜜桃钙片| 妹子高潮喷水视频| 亚洲av欧美aⅴ国产| 啦啦啦在线免费观看视频4| 精品人妻在线不人妻| 看非洲黑人一级黄片| 水蜜桃什么品种好| 最近最新中文字幕大全免费视频 | 亚洲欧洲日产国产| 天天影视国产精品| 日本vs欧美在线观看视频| 精品少妇一区二区三区视频日本电影 | 亚洲精品中文字幕在线视频| 日本-黄色视频高清免费观看| 成人漫画全彩无遮挡| 日日啪夜夜爽| 制服诱惑二区| 免费黄色在线免费观看| 男女免费视频国产| 一区二区三区四区激情视频| 深夜精品福利| 黑人猛操日本美女一级片| 欧美成人午夜精品| 日韩av不卡免费在线播放| 国产在视频线精品| 色网站视频免费| 两个人看的免费小视频| 国产淫语在线视频| 亚洲国产欧美在线一区| 久久久国产欧美日韩av| 久久久久久伊人网av| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 99久久综合免费| 国产精品二区激情视频| 成人漫画全彩无遮挡| 国产精品香港三级国产av潘金莲 | 亚洲美女黄色视频免费看| 2021少妇久久久久久久久久久| 国产一区二区三区综合在线观看| 国产精品 欧美亚洲| 精品亚洲成a人片在线观看| av免费观看日本| 精品少妇一区二区三区视频日本电影 | www.av在线官网国产| 搡老乐熟女国产| 欧美97在线视频| 久久这里只有精品19| 一区二区三区激情视频| 亚洲精品自拍成人| 一二三四在线观看免费中文在| 肉色欧美久久久久久久蜜桃| 丰满少妇做爰视频| 国产免费现黄频在线看| 少妇猛男粗大的猛烈进出视频| 国产不卡av网站在线观看| 大片电影免费在线观看免费| 少妇精品久久久久久久| 18禁裸乳无遮挡动漫免费视频| 亚洲国产成人一精品久久久| 美女主播在线视频| 麻豆乱淫一区二区| 亚洲人成77777在线视频| 亚洲精品久久午夜乱码| 午夜福利,免费看| 亚洲欧洲日产国产| 成人午夜精彩视频在线观看| 国产97色在线日韩免费| 国产片特级美女逼逼视频| 成年av动漫网址| 丝袜在线中文字幕| 精品亚洲成国产av| 日韩中文字幕欧美一区二区 | 中国国产av一级| 日韩免费高清中文字幕av| 日韩电影二区| 成人黄色视频免费在线看| 咕卡用的链子| 国产成人精品福利久久| 午夜av观看不卡| 久久国内精品自在自线图片| a级片在线免费高清观看视频| 最新中文字幕久久久久| 美女国产高潮福利片在线看| 亚洲人成网站在线观看播放| 男人操女人黄网站| 午夜福利影视在线免费观看| 精品一品国产午夜福利视频| av福利片在线| 亚洲综合精品二区| 人妻 亚洲 视频| av免费观看日本| 亚洲成人手机| 亚洲成av片中文字幕在线观看 | 波多野结衣av一区二区av| 国产极品粉嫩免费观看在线| 国产精品久久久久久久久免| 中文字幕人妻熟女乱码| 欧美 日韩 精品 国产| 街头女战士在线观看网站| 一级a爱视频在线免费观看| 99精国产麻豆久久婷婷| av国产久精品久网站免费入址| 国产在线免费精品| 日日摸夜夜添夜夜爱| 捣出白浆h1v1| 中文字幕av电影在线播放| 91精品三级在线观看| av天堂久久9| 午夜免费男女啪啪视频观看| 激情五月婷婷亚洲| 久久亚洲国产成人精品v| 夫妻午夜视频| 1024香蕉在线观看| 在线观看美女被高潮喷水网站| 亚洲第一青青草原| 黑人欧美特级aaaaaa片| 精品人妻熟女毛片av久久网站| 搡老乐熟女国产| 精品人妻在线不人妻| 国产精品三级大全| av视频免费观看在线观看| 色吧在线观看| 精品人妻一区二区三区麻豆| 99热网站在线观看| 午夜福利视频在线观看免费| 日日撸夜夜添| 国产免费又黄又爽又色| 精品一区二区三区四区五区乱码 | 丝袜美足系列| 国产熟女欧美一区二区| 中文字幕精品免费在线观看视频| 亚洲av日韩在线播放| 不卡av一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 亚洲婷婷狠狠爱综合网| 免费女性裸体啪啪无遮挡网站| 春色校园在线视频观看| av有码第一页| 一区二区三区乱码不卡18| 精品国产露脸久久av麻豆| 亚洲综合精品二区| 成人毛片60女人毛片免费| a级毛片黄视频| av网站免费在线观看视频| 一本久久精品| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 满18在线观看网站| 最近2019中文字幕mv第一页| 国产精品国产三级专区第一集| 国产日韩欧美视频二区| 一区在线观看完整版| 在线亚洲精品国产二区图片欧美| 少妇人妻久久综合中文| 国产xxxxx性猛交| 国产男女内射视频| 成年动漫av网址| 日韩中文字幕视频在线看片| av线在线观看网站| 亚洲第一区二区三区不卡| 国产精品熟女久久久久浪| 亚洲精品aⅴ在线观看| 久久久久精品性色| 久久99精品国语久久久| 中文字幕制服av| 久久精品国产自在天天线| 伊人久久大香线蕉亚洲五| 久久久精品国产亚洲av高清涩受| 国产一区有黄有色的免费视频| 各种免费的搞黄视频| 美女午夜性视频免费| 夫妻午夜视频| 在线精品无人区一区二区三| 亚洲成人手机| 麻豆乱淫一区二区| 看免费成人av毛片| 极品人妻少妇av视频| 亚洲欧美一区二区三区久久| 久久国产精品男人的天堂亚洲| 婷婷色综合大香蕉| 一级片免费观看大全| 精品人妻一区二区三区麻豆| 蜜桃国产av成人99| 深夜精品福利| 看非洲黑人一级黄片| 一边摸一边做爽爽视频免费| 肉色欧美久久久久久久蜜桃| 国产av国产精品国产| 亚洲精华国产精华液的使用体验| 极品人妻少妇av视频| 亚洲情色 制服丝袜| 国产 精品1| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 亚洲国产看品久久| 国产精品偷伦视频观看了| 国产精品二区激情视频| 看非洲黑人一级黄片| 亚洲精品日本国产第一区| 波野结衣二区三区在线| 天美传媒精品一区二区| 日韩一区二区视频免费看| 日本午夜av视频| 少妇人妻精品综合一区二区| 国产精品av久久久久免费| 亚洲欧美一区二区三区久久| 精品一品国产午夜福利视频| 久久久久视频综合| 日日啪夜夜爽| 99久久中文字幕三级久久日本| 亚洲图色成人| 日韩中字成人| 久热久热在线精品观看| av在线观看视频网站免费| 国产亚洲一区二区精品| 免费人妻精品一区二区三区视频| 国产乱来视频区| 大香蕉久久成人网| 亚洲欧美中文字幕日韩二区| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图| 久热这里只有精品99| 国产极品粉嫩免费观看在线| 欧美成人精品欧美一级黄| 欧美精品一区二区免费开放| 亚洲欧美一区二区三区国产| 人成视频在线观看免费观看| 男人爽女人下面视频在线观看| 亚洲欧美一区二区三区黑人 | 精品国产一区二区三区四区第35| 色吧在线观看| 91成人精品电影| 久久久久久伊人网av| 日日啪夜夜爽| 亚洲av.av天堂| 免费黄频网站在线观看国产| 亚洲一级一片aⅴ在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 少妇被粗大猛烈的视频| 久久精品久久久久久久性| 久久99精品国语久久久| 热re99久久精品国产66热6| 日韩中字成人| 男人舔女人的私密视频| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 汤姆久久久久久久影院中文字幕| 激情视频va一区二区三区| 一区二区三区精品91| 少妇人妻精品综合一区二区| 2022亚洲国产成人精品| 丰满饥渴人妻一区二区三| 国产精品成人在线| 久热久热在线精品观看| 日韩不卡一区二区三区视频在线| 亚洲av在线观看美女高潮| 国产97色在线日韩免费| 国产日韩欧美在线精品| 色哟哟·www| 久久毛片免费看一区二区三区| 99香蕉大伊视频| 成年av动漫网址| 日本wwww免费看| 天天躁夜夜躁狠狠躁躁| 日韩中字成人| 久久精品国产亚洲av天美| 性少妇av在线| 新久久久久国产一级毛片| 亚洲欧洲国产日韩| 18禁观看日本| 国产精品成人在线| 九色亚洲精品在线播放| 天天躁夜夜躁狠狠久久av| 日韩伦理黄色片| 一二三四中文在线观看免费高清| 午夜福利一区二区在线看| 亚洲精品乱久久久久久| 欧美 日韩 精品 国产| 欧美日韩成人在线一区二区| 777久久人妻少妇嫩草av网站| 国产福利在线免费观看视频| 国语对白做爰xxxⅹ性视频网站| 90打野战视频偷拍视频| 一区二区三区四区激情视频| 久久婷婷青草| av网站在线播放免费| 美女中出高潮动态图| 久久久国产一区二区| 精品一区二区三卡| 大片电影免费在线观看免费| 国产野战对白在线观看| 日日撸夜夜添| 久久久a久久爽久久v久久| 成人国产av品久久久| 午夜福利在线观看免费完整高清在| av国产久精品久网站免费入址| 国产一区二区激情短视频 | 免费大片黄手机在线观看| 欧美人与性动交α欧美软件| 大香蕉久久成人网| 波多野结衣av一区二区av| 久久人人爽av亚洲精品天堂| 丝袜美腿诱惑在线| 另类精品久久| 欧美精品国产亚洲| 国产在线免费精品| 国产黄色视频一区二区在线观看| 亚洲欧洲国产日韩| 18+在线观看网站| 丰满少妇做爰视频| 国产精品成人在线| 天天影视国产精品| 高清不卡的av网站| 久久免费观看电影| 亚洲图色成人| 成人18禁高潮啪啪吃奶动态图| 国产有黄有色有爽视频| 亚洲精品中文字幕在线视频| 免费高清在线观看日韩| 亚洲精品日本国产第一区| 亚洲美女黄色视频免费看| 巨乳人妻的诱惑在线观看| 久久青草综合色| 成人毛片a级毛片在线播放| 亚洲av免费高清在线观看| 成人毛片60女人毛片免费| 国产精品免费视频内射| 久热久热在线精品观看| 一区二区三区精品91| 欧美亚洲日本最大视频资源| av免费观看日本| 九九爱精品视频在线观看| 桃花免费在线播放| 久久狼人影院| av天堂久久9| 亚洲精华国产精华液的使用体验| 久久人人爽av亚洲精品天堂| 国产一区二区三区综合在线观看| 日日摸夜夜添夜夜爱| 国产精品av久久久久免费| 老司机影院成人| 超色免费av| 亚洲精品国产av蜜桃| 中文欧美无线码| 在线亚洲精品国产二区图片欧美| 国产精品亚洲av一区麻豆 | 老女人水多毛片| 国产激情久久老熟女| 欧美黄色片欧美黄色片| 国产爽快片一区二区三区| 亚洲av成人精品一二三区| 亚洲第一av免费看| 男女边吃奶边做爰视频| 黄片小视频在线播放| 丝袜美足系列| 精品人妻偷拍中文字幕| 亚洲精品第二区| 波多野结衣一区麻豆| 伊人久久大香线蕉亚洲五| 国产欧美日韩综合在线一区二区| 综合色丁香网| 国产精品一区二区在线观看99| 久久国内精品自在自线图片| 一区二区三区精品91| 高清在线视频一区二区三区| 大片电影免费在线观看免费| 亚洲精品一二三| 免费观看无遮挡的男女| 制服丝袜香蕉在线| 在线观看www视频免费| 日韩中字成人| 欧美97在线视频|