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

    黃土高原植被凈初級生產(chǎn)力的時空變化及其與氣候因子的關(guān)系*

    2016-11-16 09:02:48史曉亮楊志勇王馨爽
    中國農(nóng)業(yè)氣象 2016年4期
    關(guān)鍵詞:氣候因子黃土高原生產(chǎn)力

    史曉亮,楊志勇,王馨爽,高 軍,胡 艷

    ?

    黃土高原植被凈初級生產(chǎn)力的時空變化及其與氣候因子的關(guān)系*

    史曉亮1,楊志勇2**,王馨爽3,高 軍4,胡 艷1

    (1.西安科技大學(xué)測繪科學(xué)與技術(shù)學(xué)院,西安 710054;2.中國水利水電科學(xué)研究院水資源研究所,北京 100038;3.國家測繪地理信息局陜西基礎(chǔ)地理信息中心,西安 710054;4.中國林業(yè)科學(xué)研究院資源信息研究所,北京 100091)

    以西北典型植被脆弱區(qū)黃土高原為研究區(qū),利用AVHRR GIMMS和MODIS兩種NDVI數(shù)據(jù)源,基于CASA模型對1982-2014年黃土高原植被凈初級生產(chǎn)力(NPP)進(jìn)行模擬,并分析其時空變化特征及其與氣候因子的關(guān)系。結(jié)果表明:黃土高原年均植被NPP為254.0gC·m-2,1982-2014年總體呈增加趨勢。不同植被類型NPP有較大差異,落葉闊葉林NPP值最高,年均NPP達(dá)513.0gC·m-2,其次為常綠針葉林、草甸、農(nóng)田、灌叢和草原。黃土高原植被NPP空間分布差異顯著,表現(xiàn)出南高北低的特點(diǎn)。從NPP年際變化的空間分布來看,在退耕還林還草生態(tài)工程實(shí)施之前(1982-1998年),黃土高原大部分區(qū)域植被NPP變化不明顯。自1999年后該區(qū)植被NPP增加趨勢顯著,增速達(dá)到5.38gC·m-2·a-1。在空間分布上,66.6%的區(qū)域植被NPP呈顯著增加趨勢,主要分布在陜北高原、山西中西部的呂梁-太行山等地。退耕還林等生態(tài)工程的實(shí)施,使該區(qū)植被狀況得到改善。黃土高原植被NPP與降水量具有顯著的相關(guān)性,但與氣溫相關(guān)性不大,說明降水是影響黃土高原植被NPP的主要因素。

    凈初級生產(chǎn)力;AVHRR GIMMS NDVI;MODIS NDVI;CASA模型;黃土高原;時空變化;氣候因子

    植被凈初級生產(chǎn)力(Net Primary Productivity,NPP)直接反映植物群落在自然環(huán)境條件下的生產(chǎn)能力,作為表征植被活力的關(guān)鍵因素,NPP不僅關(guān)系到生態(tài)系統(tǒng)固碳能力的強(qiáng)弱,同時也是判定碳匯和調(diào)節(jié)生態(tài)過程的主要因子[1]。因此,隨著全球變化和碳循環(huán)研究的不斷深入,研究NPP估算方法及其與全球變化的響應(yīng)關(guān)系,在陸地生態(tài)系統(tǒng)的碳循環(huán)研究中具有重要的地位,并被國際地圈-生物圈計劃(IGBP)、人與生物圈計劃(MAB)、全球變化與陸地生態(tài)系統(tǒng)(GCTE)和京都議定書確定為核心研究內(nèi)容之一[2-3]。

    隨著遙感技術(shù)的發(fā)展,利用其時間序列長和覆蓋范圍廣的特點(diǎn),以及在反演獲取地表植被信息和關(guān)鍵參數(shù)方面的優(yōu)勢,基于遙感數(shù)據(jù)驅(qū)動的估算模型已被廣泛應(yīng)用于長時間序列NPP時空特征及其對全球變化的響應(yīng)研究[4-5]。其中建立在植物光合作用過程和光能利用率概念基礎(chǔ)上的CASA模型,以歸一化植被指數(shù)(Normal Difference Vegetation Index,NDVI)為驅(qū)動,目前已成為全球和區(qū)域尺度較為常用的植被NPP估算模型[6],并在國內(nèi)外得到廣泛應(yīng)用[7-11]。氣溫和降水等主要?dú)夂蛞蜃訉χ脖坏纳L過程具有重要影響,并對氣候變化存在明顯響應(yīng)[12]。目前在全球氣候變暖的背景下,植被NPP對氣候變化的響應(yīng)特征表現(xiàn)出較大的空間異質(zhì)性。毛德華等[9]研究認(rèn)為,中國東北地區(qū)NPP與降水間相關(guān)性空間差異顯著,而與氣溫的相關(guān)性不明顯。穆少杰等[10]研究發(fā)現(xiàn),內(nèi)蒙古不同植被類型NPP對氣候因子的敏感性差異明顯。森林植被NPP主要受溫度限制,而荒漠、農(nóng)田和草地植被NPP則主要受降水量影響。周杰等[13]在渭河流域的研究則表明溫度和降水是影響渭河流域NPP的主要因子。因此,有必要在典型區(qū)域開展植被NPP研究,對于深入探討區(qū)域水平植被NPP時空演變的驅(qū)動機(jī)制,豐富全球碳循環(huán)研究具有重要意義。

    黃土高原地處半干旱半濕潤氣候帶,土質(zhì)疏松,水土流失嚴(yán)重,是中國的生態(tài)脆弱區(qū)和水土保持重點(diǎn)區(qū)域[14]。為改善黃土高原日益惡化的生態(tài)環(huán)境,從1999年開始國家在該地區(qū)實(shí)施退耕還林還草、封山育林等大規(guī)模植被建設(shè),植被得到一定的恢復(fù)[15]。目前已有學(xué)者在黃土高原開展了植被NPP研究。韋莉等[16]基于CASA模型估算了2003年黃土高原草地NPP;宋富強(qiáng)等[17]研究了2005年陜北黃土高原地區(qū)的植被凈初級生產(chǎn)力及其時空分布;李登科等[18]基于MOD17A3數(shù)據(jù)集的研究認(rèn)為陜西省2000-2006年大部分地區(qū)NPP呈增加趨勢,表明通過退耕還林等生態(tài)建設(shè)工程的實(shí)施,該地區(qū)植被狀況得到較好的改善。謝寶妮等[19]研究認(rèn)為黃土高原2000-2010年NPP呈增加趨勢,由大規(guī)模退耕還林還草植被建設(shè)工程引起的土地利用覆被變化是退耕區(qū)NPP增加的主要因素。但是以往研究多局限于單一遙感數(shù)據(jù)源,從而導(dǎo)致研究時間序列較短,較少涉及對黃土高原地區(qū)在大規(guī)模生態(tài)工程建設(shè)實(shí)施前后植被NPP時空演變特征的研究。鑒于此,本文綜合利用AVHRR GIMMS和MODIS兩種NDVI數(shù)據(jù)源,估算黃土高原1982-2014年長時間序列植被NPP,并以1999年為分界點(diǎn),分析探討該區(qū)大規(guī)模植被建設(shè)前后植被NPP的時空變化及其與氣候因子的關(guān)系,以期為黃土高原退耕還林還草生態(tài)工程成效評價提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    黃土高原地區(qū)(100°54′-114°33′E、33°43′- 41°16′N)為太行山以西、青海省日月山以東,秦嶺以北、長城以南的廣大地區(qū),行政區(qū)涉及陜西省、山西省、甘肅省、青海省、河南省、寧夏回族自治區(qū)以及內(nèi)蒙古自治區(qū)等?。▍^(qū)),總面積約63.2萬km2。

    研究區(qū)地形西高東低,海拔200-3000m。地處半干旱半濕潤氣候帶,屬暖溫帶大陸性季風(fēng)氣候,冬季寒冷干燥,夏季炎熱少雨,且雨季集中,年際變率較大。年均氣溫3.6~14.3℃,年均降水量300~800mm。植被類型主要包括常綠針葉林、落葉闊葉林、灌叢、草地和農(nóng)田(圖1),主要土壤類型為黃綿土、草原風(fēng)沙土和褐土。

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

    遙感數(shù)據(jù)源包括AVHRR GIMMS和MODIS兩種NDVI數(shù)據(jù)集。GIMMS NDVI數(shù)據(jù)來源于國家自然科學(xué)基金委員會“中國西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心”(http://westdc.westgis.ac.cn)??臻g分辨率為8km,為每15d最大值合成數(shù)據(jù),時間尺度為1982-2006年。數(shù)據(jù)集經(jīng)過輻射和幾何校正等預(yù)處理,為減少云、大氣、太陽高度角等的影響,采用最大合成法(Maximum Value Composites,MVC)獲取逐月NDVI值[20]。MODIS NDVI數(shù)據(jù)來自美國NASA地球觀測系統(tǒng)的2001-2014年的MOD13A3月值NDVI數(shù)據(jù),空間分辨率為1km。通過專業(yè)處理軟件MRT(MODIS Reprojection Tools)對下載獲取的MODIS NDVI數(shù)據(jù)完成鑲嵌、裁切、投影與格式轉(zhuǎn)換等處理,并基于兩種NDVI數(shù)據(jù)源重合時間段的NDVI數(shù)據(jù)集(2001-2006年),通過建立逐像元一元線性回歸模型,應(yīng)用重采樣為8km分辨率的MODIS NDVI對GIMMS NDVI數(shù)據(jù)集進(jìn)行插補(bǔ)延長[21],從而建立黃土高原33a(1982-2014年)長時間序列8km分辨率的NDVI數(shù)據(jù)集?;谥袊茖W(xué)院植被圖編輯委員會2001年編著的1:100萬中國植被圖集,經(jīng)掃描數(shù)字化得到黃土高原植被類型圖。

    對模擬結(jié)果進(jìn)行驗(yàn)證的NPP產(chǎn)品數(shù)據(jù)來自美國NASA EOS/MODIS 2000-2014年MOD17A3數(shù)據(jù),時間分辨率為年,空間分辨率為1km。MOD17A3數(shù)據(jù)是參考BIOME-BGC模型與光能利用率模型模擬得到的年NPP,并已在全球和區(qū)域碳循環(huán)研究中得到廣泛應(yīng)用[22]。

    氣象數(shù)據(jù)來源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http://cdc.nmic.cn),包括研究區(qū)45個氣象站點(diǎn)1982-2014年的逐月降水量、平均氣溫和日照時數(shù)。太陽總輻射量由大氣上界太陽輻射量和日照時數(shù)計算得到。采用IDW插值方法,獲取與NDVI數(shù)據(jù)具有相同投影和像元大小的月太陽輻射、降水量和平均氣溫柵格數(shù)據(jù)集。

    1.3 NPP估算模型

    采用CASA(Carnegie-Ames-Stanford Approach)模型對黃土高原1982-2014年的植被NPP進(jìn)行估算。該模型是一個充分考慮環(huán)境條件和植被本身特征的光能利用率模型,由遙感、氣象、植被以及土壤類型數(shù)據(jù)共同驅(qū)動,已被廣泛應(yīng)用于陸地生態(tài)系統(tǒng)NPP的估算[10,23],可以由植被吸收的光合有效輻射(APAR)和光能利用率(e)兩個因子來確定。

    式中,NPP(x,t)是像元x在t月份的植被凈初級生產(chǎn)力,以碳為單位(gC·m-2),APAR(x,t)表示像元x在t月份吸收的光合有效輻射(MJ·m-2),e(x,t)表示像元x在t月份的實(shí)際光能利用率(gC·MJ-2),與植被類型有關(guān)。植被吸收的光合有效輻射(APAR)計算方法為

    式中,S(x,t)表示像元x在t月份的太陽總輻射量(MJ·m-2),可由大氣上界太陽輻射量和日照百分比計算[24]。常數(shù)0.5表示植被利用的太陽有效輻射占太陽總輻射的比例,計算式為[25]

    式中,SRmin取值為1.08;SRmax的大小與植被類型相關(guān),可由文獻(xiàn)[26]獲取。SR由NDVI計算得到。

    光能利用率e表示植被將吸收的光合有效輻射轉(zhuǎn)化為有機(jī)碳的效率(gC·MJ-2),可由下式計算。

    式中,Te1(x,t)和Te2(x,t)分別表示低溫和高溫對光能利用率的影響系數(shù),其中,Te1為低溫和高溫情況下植物內(nèi)在生化作用對光合的限制,由下式計算:

    式中,Topt(x)表示研究區(qū)1a內(nèi)NDVI值達(dá)到最大時當(dāng)月的平均氣溫。Te2(x,t)表示環(huán)境溫度由最適宜溫度向高溫和低溫變化時植被光能利用率逐漸變小的趨勢。

    We(x,t)為水分脅迫影響系數(shù),表示水分條件對光能利用率的影響。隨著環(huán)境中有效水分的增加,其取值范圍為0.5(極端干旱條件)~1(非常濕潤條件),計算式為

    式中,EET(x,t)為實(shí)際蒸散量(mm),可根據(jù)區(qū)域?qū)嶋H蒸散模型求取;PET(x,t)為潛在蒸散量(mm),可根據(jù)Boucher提出的互補(bǔ)關(guān)系求 取[24]。

    emax表示在理想條件下植被的最大光能利用率,其取值因植被類型不同有較大差別,本文emax取值參考朱文泉等[26]在中國范圍內(nèi)的研究結(jié)果確定黃土高原地區(qū)各植被類型的最大光能利用率(表1)。

    利用以上模型和NDVI、氣象資料,估算得到黃土高原764個格網(wǎng)點(diǎn)1982-2014年逐月植被NPP值,每年1-12月值累加得到相應(yīng)的年NPP。以所有格點(diǎn)的算術(shù)平均值作為整個研究區(qū)對應(yīng)時段的NPP值。

    表1 CASA模型中emax取值

    Note: ECF is evergreen coniferous forest;DBF is deciduous broadleaf forest. The same as below.

    1.4 分析方法

    應(yīng)用一元線性回歸法在像元尺度分析黃土高原在大規(guī)模植被建設(shè)前后兩個階段(1982-1998和1999-2014年)植被NPP年值的變化趨勢[10]。

    基于像元相關(guān)分析方法對1982-2014年黃土高原植被NPP年值與主要?dú)夂蛞蜃樱昃邓亢推骄鶜鉁兀┻M(jìn)行逐像元相關(guān)分析[20],用相關(guān)系數(shù)的空間分布情況反映植被NPP對氣候因子的響應(yīng)。

    2 結(jié)果與分析

    2.1 植被NPP估算精度評價

    由于獲取研究區(qū)NPP的實(shí)測值較難,本文采用與MOD17A3產(chǎn)品數(shù)據(jù)對比的方法進(jìn)行模型模擬精度評價[23]。將本研究估算得到的不同類型的NPP均值與部分實(shí)測數(shù)據(jù)進(jìn)行對比,如果偏差不大,即認(rèn)為估算結(jié)果較合理[27]。由表2可見,對于落葉闊葉林,本研究估算的NPP值與MOD17A3數(shù)據(jù)差別較大,但在實(shí)測范圍之內(nèi)。其它植被類型差異較小,且反映的趨勢較為一致,說明本研究估算的植被NPP基本可靠。

    2.2 黃土高原植被NPP時空變化特征

    2.2.1 研究區(qū)域植被NPP平均值的年際變化

    由圖2可見,黃土高原1982-2014年植被NPP的年際波動變化較大,但整體呈增加趨勢,年增速為1.7gC·m-2(P<0.05)。年NPP值在1982年為最低(214.6gC·m-2),2012年最高(322.0gC·m-2),平均254.0gC·m-2。其中1982-1998年研究區(qū)植被NPP年均增速為1.2gC·m-2,年均NPP為244.1gC·m-2。但是自1999年實(shí)施大規(guī)模生態(tài)工程建設(shè)以來,該區(qū)植被NPP隨時間推移呈現(xiàn)明顯增加趨勢,1999-2014年年均增速達(dá)到5.4gC·m-2,年均NPP增至264.6gC·m-2。

    表2 不同植被類型NPP年均估算值與實(shí)測值的比較(gC·m-2)

    2.2.2 大規(guī)模植被建設(shè)前后不同植被類型NPP的變化

    1982-2014年不同植被類型年均NPP存在較大差異,統(tǒng)計比較表明,黃土高原生態(tài)系統(tǒng)中落葉闊葉林NPP值最高,33a間年均NPP達(dá)513.0gC·m-2,但是區(qū)域內(nèi)落葉闊葉林分布較少,僅占全區(qū)總面積的2.63%。

    常綠針葉林為346.0gC·m-2,草甸為300.1gC·m-2,農(nóng)田為261.2gC·m-2,灌叢為271.7gC·m-2,草原最小,為235.4gC·m-2。后3種是黃土高原的主要植被類型,占總面積的88.53%。

    相對于退耕還林還草工程實(shí)施之前(1982-1998年),經(jīng)過近15a的大規(guī)模植被建設(shè),黃土高原不同植被類型的NPP均呈不同程度的上升趨勢(表3)。其中,落葉闊葉林的NPP增長趨勢最大,達(dá)到32.4gC·m-2·a-1。常綠針葉林和草甸、灌叢等植被類型NPP增長也較明顯。

    表3 黃土高原大規(guī)模植被建設(shè)實(shí)施(1998年)前后不同植被類型年均NPP(gC·m-2)

    2.2.3 整個分析期植被年NPP平均值空間分布

    1982-2014年黃土高原植被年均NPP分布受到植被類型、氣候和地形等多種因素的綜合影響,空間差異顯著,總體為南高北低,分布在15.12~940.21gC·m-2。從圖3可以看出,在內(nèi)蒙古自治區(qū)的鄂爾多斯地區(qū)、河套平原,甘肅省隴中地區(qū)和寧夏回族自治區(qū)南部大部分地區(qū),以及陜北地區(qū),上述區(qū)域主要植被類型為灌叢、草原和荒漠等無植被地區(qū),植被覆蓋相對較少,自然環(huán)境惡劣,植被年均NPP多在200gC·m-2以下,相對較低。而在青海省南部、陜西省關(guān)中地區(qū)以及山西省大部分地區(qū),主要分布植被類型多為農(nóng)田和落葉闊葉林,植被年均NPP大于200gC·m-2,尤以陜西省的銅川、寶雞和渭南等地,植被年均NPP大于600gC·m-2。

    2.2.4 植被年NPP變化趨勢的空間分布

    采用逐像元線性趨勢法對黃土高原大規(guī)模生態(tài)工程建設(shè)實(shí)施前后的植被NPP變化趨勢進(jìn)行分析,從圖4a可以看出,在大規(guī)模植被建設(shè)實(shí)施前(1982-1998年),雖然黃土高原73.3%的區(qū)域植被年均NPP有所增加,但呈顯著增加趨勢的區(qū)域占整個黃土高原的面積比例為20.4%,主要分布在內(nèi)蒙古南部、山西省北部,甘肅省隴東和陜北地區(qū)等,其它大部分地區(qū)植被NPP變化不明顯。由此可見,在大規(guī)模植被建設(shè)前,黃土高原大部分區(qū)域植被年NPP無顯著變化。而自1999年開始經(jīng)過16a的大規(guī)模植被建設(shè),該區(qū)在1999-2014年88.3%的區(qū)域植被年NPP有所增加,其中呈顯著增加趨勢的區(qū)域面積比例達(dá)到了66.6%,平均增速達(dá)到7.49gC·m-2·a-1(P<0.05)(圖4b)。另有11.7%的地區(qū)植被NPP變化趨勢不明顯,主要分布在研究區(qū)西北部的干旱荒漠地區(qū),以及陜西省關(guān)中地區(qū)的西安、寶雞等地??梢钥闯?,黃土高原植被的生長狀況整體有所改善。

    2.3 黃土高原植被NPP與氣候因子的相關(guān)分析

    降水和氣溫等氣候因子的變化對陸地植被的生長發(fā)育具有重要影響,因此本文利用逐像元的相關(guān)分析方法,分別對1982-2014年的年植被NPP、年降水量和年平均氣溫的相關(guān)系數(shù)進(jìn)行逐像元分析,從而在空間尺度分析植被NPP與主要?dú)夂蛞蜃拥年P(guān)系。如圖5a所示,整個黃土高原植被NPP與降水量的相關(guān)系數(shù)為0.39,呈顯著正相關(guān)(P<0.05)。而從像元尺度來看,除了陜西關(guān)中地區(qū),以及其它?。ㄗ灾螀^(qū))的部分地區(qū)外,黃土高原大部分地區(qū)植被NPP與降水之間表現(xiàn)為顯著正相關(guān)(P<0.05),呈顯著正相關(guān)(P<0.05)的區(qū)域面積占全區(qū)總面積的比例達(dá)到65.3%。尤其是在內(nèi)蒙古自治區(qū)鄂爾多斯地區(qū),甘肅的隴東以及陜北的榆林等地植被NPP與降水量的相關(guān)性較高。

    從圖5b可以看出,整個黃土高原地區(qū)植被NPP與年平均氣溫的相關(guān)系數(shù)為0.05(P>0.05),不存在明顯的相關(guān)性。而從像元尺度來看,研究區(qū)內(nèi)植被NPP與氣溫呈顯著正相關(guān)的地區(qū)占總面積的9.8%(P<0.05),主要分布在內(nèi)蒙古自治區(qū)鄂爾多斯地區(qū)、山西省北部和陜西省銅川等地,研究范圍內(nèi)植被NPP與氣溫呈顯著負(fù)相關(guān)的地區(qū)占總面積的4.9%,主要分布在山西省中部,以及寧夏、甘肅和內(nèi)蒙古等地,其余地區(qū)(占像元面積85.3%)兩者間的相關(guān)系數(shù)未通過0.05水平的顯著性檢驗(yàn)。說明黃土高原大部分地區(qū)植被NPP與氣溫之間相關(guān)性不大。

    3 結(jié)論與討論

    3.1 討論

    黃土高原1999-2014年降水量呈緩慢增加趨勢,而氣溫略有降低[15],且大部分地區(qū)降水與NPP呈正相關(guān),因此該區(qū)降水的增加在一定程度上導(dǎo)致植被NPP增加。但是NPP與氣溫相關(guān)性不大,氣溫的增加對植被NPP影響較小,因此可以認(rèn)為,降水量的增加是黃土高原植被NPP增加的主要因素。

    由于不同NPP反演模型在機(jī)理和關(guān)鍵植被參數(shù)取值、研究時間序列長短和模型輸入數(shù)據(jù)分辨率和質(zhì)量等方面存在的差異,導(dǎo)致植被NPP反演模型結(jié)果有所不同。本文基于CASA模型估算的1982-2014年黃土高原植被年平均NPP為254.0gC·m-2,與謝寶妮等[19]研究結(jié)果基本一致。此外,對比朱文泉等[26]提供的各植被類型NPP的實(shí)測值,本研究的模擬結(jié)果也在其范圍之內(nèi)。

    以往關(guān)于黃土高原植被NPP的模擬研究[16-19],由于遙感數(shù)據(jù)源時間尺度限制,未能分析該區(qū)在大規(guī)模生態(tài)工程實(shí)施前后植被NPP的時空演變。而本文綜合利用AVHRR GIMMS和MODIS兩種NDVI數(shù)據(jù)源,并應(yīng)用逐像元線性回歸模型方法對GIMMS NDVI數(shù)據(jù)進(jìn)行插補(bǔ)延長,從而實(shí)現(xiàn)了黃土高原長時間序列植被NPP的估算,對于評價退耕還林還草等生態(tài)工程的實(shí)施效應(yīng)更具參考價值。本文分析了植被NPP與主要?dú)夂蛞蜃拥捻憫?yīng)關(guān)系,但是NPP的變化還與地形和人類活動等多種因素有關(guān),因此,對于植被NPP的驅(qū)動機(jī)制尚需進(jìn)一步深入研究。

    3.2 結(jié)論

    (1)黃土高原1982-2014年植被年平均NPP為254.0gC·m-2,總體呈增加趨勢,尤以1999年該區(qū)實(shí)施大規(guī)模生態(tài)環(huán)境工程建設(shè)以來,NPP增加趨勢更為明顯。對于不同植被類型而言,落葉闊葉林NPP值最高,其次為常綠針葉林、草甸、農(nóng)田、灌叢和草原。

    (2)黃土高原植被NPP空間分布差異顯著,總體為南高北低。而從植被NPP變化趨勢的空間分布來看,在大規(guī)模植被建設(shè)實(shí)施前(1982-1998年),該區(qū)大部分區(qū)域植被NPP無顯著變化。而自1999年開始經(jīng)過16a大規(guī)模生態(tài)工程的建設(shè)實(shí)施,黃土高原在1999-2014年66.6%的區(qū)域植被NPP呈顯著增加趨勢,年增速達(dá)到7.49gC·m-2(P<0.05)。尤其以陜北高原、山西中西部的呂梁-太行山等地,植被NPP增加趨勢更為明顯。

    (3)黃土高原植被NPP與降水量具有較強(qiáng)的正相關(guān),而與氣溫相關(guān)性不大,因此,降水是該區(qū)植被NPP的主要影響因素。

    參考文獻(xiàn)References

    [1]樸世龍,方精云,郭慶華.利用CASA模型估算我國植被凈第一性生產(chǎn)力[J].植物生態(tài)學(xué)報,2001, 25(5):603-608.

    Piao S L,Fang J Y,Guo Q H.Application of CASA model to the estimation of Chinese terrestrial net primary productivity[J].Acta Phytoecologica Sinica,2001,25(5):603- 608.(in Chinese)

    [2]Field C B,Randerson J T,Malmstrom C M.Global net primary production:combining ecology and remote sensing[J].Remote Sensing of Environment,1995,51(1):74-88.

    [3]IGBP.The terrestrial carbon cycle:implications for the Kyoto Protocol[J].Science,1998,280:1393-1394.

    [4]Ricotta C,Avena G,de Palma A.Mapping and monitoring net primary productivity with AVHRR NDVI time-series: statistical equivalence of cumulative vegetation indices[J]. ISPRS Journal of Photogrammetry & Remote Sensing,1999, 54(5/6):325-331.

    [5]Nemani R R,Keeling C D,Hashimoto H,et al.Climate-driven increases in global terrestrial net primary production from1982 to 1999[J].Science,2003,300(5625):1560-1563.

    [6]潘竟虎,李真.2001-2012年西北干旱區(qū)植被凈初級生產(chǎn)力時空變化[J].生態(tài)學(xué)雜志,2015,34(12):3333-3340.

    Pan J H,Li Z.Temporal-spatial change of vegetation net primary productivity in the arid region of Northwest China during 2001 and 2012[J].Chinese Journal of Ecology,2015, 34(12): 3333-3340.(in Chinese)

    [7]李秀芬,鄭有飛,王晨軼,等.黑龍江省森林NPP的遙感反演[J].中國農(nóng)業(yè)氣象,2009,31(1):88-92.

    Li X F,Zheng Y F,Wang C Y,et al.Retrieval of forest NPP by remote sensing in Heilongjiang Province[J].Chinese Journal of Agrometeorology,2009,31(1):88-92.(in Chinese)

    [8]陳福軍,沈彥俊,李倩,等.中國陸地生態(tài)系統(tǒng)近30年NPP時空變化研究[J].地理科學(xué),2011,31(11):1409-1414.

    Chen F J,Shen Y J,Li Q,et al.Spatio-temporal variation analysis of ecological systems NPP in China in Past 30 years[J].Scientia Geo Graphica Sinica,2011,31(11):1409-1414. (in Chinese)

    [9]毛德華,王宗明,羅玲,等.1982-2009年東北多年凍土區(qū)植被凈初級生產(chǎn)力動態(tài)及其對全球變化的響應(yīng)[J].應(yīng)用生態(tài)學(xué)報,2012,23(6):1511-1519.

    Mao D H,Wang Z M,Luo L,et al.Dynamic changes of vegetation net primary productivity in permafrost zone of Northeast China in 1982-2009 in response to global change[J].Chinese Journal of Applied Ecology,2012,23(6): 1511-1519.(in Chinese)

    [10]穆少杰,李建龍,周偉,等.2001-2010年內(nèi)蒙古植被凈初級生產(chǎn)力的時空格局及其與氣候的關(guān)系[J].生態(tài)學(xué)報,2013,33 (12):3752-3764.

    Mu S J,Li J L,Zhou W,et al.Spatial-temporal distribution of net primary productivity and its relationship with climate factors in Inner Mongolia from 2001 to 2010[J].Acta Ecologica Sinica,2013,33(12):3752-3764.(in Chinese)

    [11]石志華,劉夢云,吳健立,等.基于CASA模型的陜西省凈初級生產(chǎn)力時空分析[J].水土保持通報,2016,36(1):206-211.

    Shi Z H,Liu M Y,Wu J L,et al.Spatial-temporal analysis of vegetation net primary productivity in Shaanxi province based on CASA model[J].Bulletin of Soil and Water Conservation,2016,36(1):206-211.(in Chinese)

    [12]張禹舜,賈文雄,劉亞榮,等.近11a來祁連山凈初級生產(chǎn)力對氣候因子的響應(yīng)[J].干旱區(qū)地理,2016,39(1):77-85.

    Zhang Y S,Jia W X,Liu Y R,et al.Response of net primary productivity to climatic factors in Qilian Mountains in resent eleven years[J].Arid Land Geography,2016,39(1):77-85.(in Chinese)

    [13]周杰,張昌盛,解修平.渭河流域植被凈初級生產(chǎn)力時空格局及影響因素分析[J].水土保持學(xué)報,2015,29(2):274-277.

    Zhou J,Zhang C S,Xie X P.Vegetation net primary productivity temporal and spatial patterns and influence factors analysis in Weihe watershed[J].Journal of Soil and Water Conservation,2015,29(2):274-277.(in Chinese)

    [14]易浪,任志遠(yuǎn),張翀,等.黃土高原植被覆蓋變化與氣候和人類活動的關(guān)系[J].資源科學(xué),2014,36(1):166-174.

    Yi L,Ren Z Y,Zhang C,et al.Vegetation cover,climate and human activities on the Loess Plateau[J].Resources Science,2014,36(1):166-174.(in Chinese)

    [15]張寶慶,吳普特,趙西寧.近30a黃土高原植被覆蓋時空演變監(jiān)測與分析[J].農(nóng)業(yè)工程學(xué)報,2011,27(4):287-293.

    Zhang B Q,Wu P T,Zhao X N.Detecting and analysis of spatial and temporal variation of vegetation cover in the Loess Plateau during 1982-2009[J].Transactions of the CSAE,2011,27(4):287-293.(in Chinese)

    [16]韋莉,趙軍,潘竟虎,等.基于MODIS數(shù)據(jù)的黃土高原草地凈初級生產(chǎn)力的估算研究[J].遙感技術(shù)與應(yīng)用,2009,24(5): 660-664.

    Wei L,Zhao J,Pan J H,et al.A research of net Primary productivity model of grassland based on MODIS data in the Loess Plateau of China[J].Remote Sensing Technology and Application,2009,24(5):660-664.(in Chinese)

    [17]宋富強(qiáng),康慕誼,陳雅如,等.陜北黃土高原植被凈初級生產(chǎn)力的估算[J].生態(tài)學(xué)雜志,2009,28(11):2311-2318.

    Song F Q,Kang M Y,Chen Y R,et al.Estimation of vegetation net primary productivity on North Shaanxi Loess Plateau[J].Chinese Journal of Ecology,2009,28(11):2311- 2318.(in Chinese)

    [18]李登科,范建忠,王娟.基于MOD17A3的陜西省植被NPP變化特征[J].生態(tài)學(xué)雜志,2011,30(12):2776-2782.

    Li D K,Fan J Z,Wang J.Variation characteristics of vegetation net primary productivity in Shaanxi Province based on MO17A3[J].Chinese Journal of Ecology,2011,30 (12):2776-2782.(in Chinese)

    [19]謝寶妮,秦占飛,王洋,等.黃土高原植被凈初級生產(chǎn)力時空變化及其影響因素[J].農(nóng)業(yè)工程學(xué)報,2014, 30(11):244-253.

    Xie B N,Qin Z F,Wang Y,et al.Spatial and temporal variation in terrestrial net primary productivity on Chinese Loess Plateau and its influential factors[J].Transactions of the CSAE,2014,30(11):244-253.(in Chinese)

    [20]史曉亮,李穎,趙凱,等.諾敏河流域植被覆蓋時空演變及其與徑流的關(guān)系研究[J].干旱區(qū)資源與環(huán)境, 2013,27(6): 54-60.

    Shi X L,Li Y,Zhao K,et al.Temporal and spatial changes of vegetation cover and the relationship with runoff in Nuomin river basin[J].Journal of Arid Land Resources and Environment,2013,27(6):54-60.(in Chinese)

    [21]毛德華,王宗明,宋開山,等.東北多年凍土區(qū)植被NDVI變化及其對氣候變化和土地覆蓋變化的響應(yīng)[J].中國環(huán)境科學(xué),2011,31(2):283-292.

    Mao D H,Wang Z M,Song K S,et al.The vegetation NDVI variation and its responses to climate change and LUCC from 1982 to 2006 year in northeast permafrost region[J].China Environmental Science,2011,31(2):283-292. (in Chinese)

    [22]王新闖,王世東,張合兵.基于MOD17A3的河南省NPP時空格局[J].生態(tài)學(xué)雜志,2013,32(10):2797- 2805.

    Wang X C,Wang S D,Zhang H B.Spatiotemporal pattern of vegetation net primary productivity in Henan Province of China based on MOD17A3[J].Chinese Journal of Ecology,2013,32(10):2797-2805.(in Chinese)

    [23]沃笑,吳良才,張繼平,等.基于CASA模型的三江源地區(qū)植被凈初級生產(chǎn)力遙感估算研究[J].干旱區(qū)資源與環(huán)境,2014, 28(9):45-50.

    Wo X,Wu L C,Zhang J P,et al.Estimation of net primary production in the Three-River headwater region using CASA model[J].Journal of Arid Land Resources and Environment, 2014,28(9):45-50.(in Chinese)

    [24]羅玲.基于遙感-機(jī)理模型的松嫩西部草地凈初級生產(chǎn)力(NPP)模擬[D].長春:中國科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所,2010.

    Luo L.Estimation of grassland Net Primary Production(NPP) in the western Songnen Plain based on the remote sensing-mechanism model[D].Changchun:Northeast Institute of Geography and Agroecology,Chinese Academy of Sciences,2010.(in Chinese)

    [25]朱文泉,潘耀忠,龍中華,等.基于GIS和RS的區(qū)域陸地植被NPP估算:以中國內(nèi)蒙古為例[J].遙感學(xué)報,2005,9(3): 300-307.

    Zhu W Q,Pan Y Z,Long Z H,et al.Estimating net primary productivity of terrestrial vegetation based on GIS and RS:a case study in Inner Mongolia,China[J].Journal of Remote Sensing,2005,9(3):300-307.(in Chinese)

    [26]朱文泉,潘耀忠,張錦水.中國陸地植被凈初級生產(chǎn)力遙感估算[J].植物生態(tài)學(xué)報,2007,31(3):413-424.

    Zhu W Q,Pan Y Z,Zhang J S.Estimation of net primary productivity of Chinese terrestrial vegetation based on remote sensing[J].Journal of Plant Ecology,2007,31(3): 413-424.(in Chinese)

    [27]Yuan Q Z,Wu S H,Zhao D S,et al.Modeling net primary productivity of the terrestrial ecosystem in China from 1961-2005[J].Journal of Geographical Science,2014,24(1): 3-17.

    Spatial and Temporal Variation of Net Primary Productivity and Its Relationship with Climate Factors in the Chinese Loess Plateau

    SHI Xiao-liang1, YANG Zhi-yong2, WANG Xin-shuang3, GAO Jun4, HU Yan1

    (1.College of Geomatics, Xi’an University of Science and Technology, Xi’an 710054, China;2.Department of Water Resources, China Institute of Water Resources and Hydropower Research, Beijing 100038; 3.Shannxi Geomatics Center of National Administration of Surveying, Mapping and Geoinformation, Xi’an 710054;4.Institute of Forest Resources Information Techniques, Chinese Academy of Forestry, Beijing 100091)

    Based on AVHRR GIMMS NDVI and MODIS NDVI data, spatial-temporal variations of the net primary productivity (NPP) in the Chinese Loess Plateau from 1982 to 2014 was simulated using CASA model. Besides, the relationship between NPP and climate factors was analyzed at grid scale. The results showed that average annual NPP from 1982 to 2014 was 254.0gC·m-2in the Chinese Loess Plateau and the NPP has increased over the past three decades. The average annual NPP of deciduous broadleaf forest was 513.0gC·m-2, which was the highest for different vegetation types, followed by evergreen coniferous forest, meadow, cropland, shrub and steppe. The NPP differs greatly between different vegetation types. There was significant spatial difference of vegetation NPP, and the NPP of the south parts of the study area was higher than that of the north part. Since 1999, the government of China began to carry out the police that is “replace agriculture with forestation and conserve forest” in the study area. Before the project of returning farmland to forestland or grassland (1982-1998), there was no significant change of vegetation NPP in most regions of the study area. However, since the implementation of large-scale vegetation construction from 1999, the annual mean NPP of the study area has increased significantly at a rate of 5.38gC·m-2. About 66.6% of the study area showed an increasing NPP trend, especially in northern Shaanxi Plateau, Taihang-Lvliang mountains, where significant increase in vegetation NPP has been seen. The green for grain project has greatly improved the vegetation state. The NPP had a significant positive correlation with the precipitation, but it had no obvious relationship with temperature. Therefore, precipitation is the main factor driving vegetation NPP change.

    Net primary productivity(NPP); AVHRR GIMMS NDVI; MODIS NDVI; CASA model; The Chinese Loess Plateau; Spatial and temporal variation; Climate factors

    10.3969/j.issn.1000-6362.2016.04.009

    2016-03-20

    通訊作者。E-mail:yangzy@iwhr.com

    國家自然科學(xué)基金項(xiàng)目(51409204)“流域土地利用/覆被變化對干旱的影響評價及綜合調(diào)控研究”;陜西省教育廳科學(xué)研究計劃項(xiàng)目(15JK1473)“基于CASA模型的黃土高原植被凈初級生產(chǎn)力的遙感估算研究”;西安科技大學(xué)科研培育基金項(xiàng)目(201605)“基于分布式水文模擬和遙感技術(shù)的干旱評價方法研究”

    史曉亮(1985-),講師,主要從事資源環(huán)境遙感研究。E-mail:s_xiaoliang@126.com

    猜你喜歡
    氣候因子黃土高原生產(chǎn)力
    氣候因子對烤煙質(zhì)量風(fēng)格特色的影響
    把“風(fēng)景”轉(zhuǎn)化成“生產(chǎn)力
    紅土地(2016年3期)2017-01-15 13:45:17
    選舉 沸騰了黃土高原(下)
    公民與法治(2016年3期)2016-05-17 04:09:00
    選舉沸騰了黃土高原(上)
    公民與法治(2016年1期)2016-05-17 04:07:56
    四川水稻產(chǎn)量及其構(gòu)成要素對不同生育期氣候因子的響應(yīng)分析
    灑向黃土高原的愛
    中國火炬(2015年7期)2015-07-31 17:39:57
    人口紅利與提高生產(chǎn)力
    海峽姐妹(2015年9期)2015-02-27 15:13:07
    天山北坡NDVI對氣候因子響應(yīng)的敏感性分析
    展望2014:推動科技創(chuàng)新轉(zhuǎn)化為第一生產(chǎn)力
    杭州科技(2014年3期)2014-02-27 15:26:46
    沙漠地區(qū)微波地表發(fā)射率年內(nèi)變化規(guī)律與氣候因子的關(guān)系分析
    美女中出高潮动态图| 最近最新中文字幕大全免费视频 | 久久99蜜桃精品久久| 卡戴珊不雅视频在线播放| 高清毛片免费看| 亚洲熟女精品中文字幕| 在线免费观看不下载黄p国产| 久久毛片免费看一区二区三区| 国产精品一国产av| 婷婷成人精品国产| 久久精品国产自在天天线| 欧美日韩亚洲高清精品| 午夜av观看不卡| 婷婷成人精品国产| 久久精品久久久久久噜噜老黄| 18禁观看日本| 精品一区二区免费观看| 男女午夜视频在线观看 | 免费黄色在线免费观看| 日韩av不卡免费在线播放| 最近2019中文字幕mv第一页| 大话2 男鬼变身卡| 国产成人aa在线观看| 国产极品天堂在线| 蜜桃国产av成人99| 少妇人妻 视频| 涩涩av久久男人的天堂| 久久久久视频综合| h视频一区二区三区| 丰满乱子伦码专区| 亚洲性久久影院| 岛国毛片在线播放| 亚洲中文av在线| 国产色爽女视频免费观看| 三级国产精品片| 日本猛色少妇xxxxx猛交久久| 久久这里有精品视频免费| 国产精品一区二区在线不卡| 亚洲国产精品999| 一二三四中文在线观看免费高清| 免费av不卡在线播放| 亚洲精品国产av成人精品| 成年美女黄网站色视频大全免费| 热re99久久国产66热| 老司机影院毛片| 老司机亚洲免费影院| 一级,二级,三级黄色视频| 丰满饥渴人妻一区二区三| 国产精品一国产av| 精品亚洲乱码少妇综合久久| 一边亲一边摸免费视频| 18+在线观看网站| av国产精品久久久久影院| 欧美日韩综合久久久久久| 日韩欧美精品免费久久| 天天操日日干夜夜撸| 亚洲精品成人av观看孕妇| 99香蕉大伊视频| 国产色爽女视频免费观看| 国产片内射在线| 亚洲国产日韩一区二区| 欧美人与性动交α欧美软件 | 最近中文字幕高清免费大全6| 久久久久久久久久人人人人人人| 成年人午夜在线观看视频| 大香蕉久久成人网| 99久久人妻综合| 国国产精品蜜臀av免费| 看免费成人av毛片| 老司机影院成人| 一区在线观看完整版| 老女人水多毛片| 精品一区在线观看国产| 久久久久久久久久久免费av| 久久人妻熟女aⅴ| 国产深夜福利视频在线观看| 最近的中文字幕免费完整| av卡一久久| 久久ye,这里只有精品| 美女主播在线视频| 日韩欧美一区视频在线观看| 亚洲高清免费不卡视频| 婷婷色麻豆天堂久久| 在线观看人妻少妇| 狂野欧美激情性xxxx在线观看| 亚洲av电影在线进入| 最近2019中文字幕mv第一页| 国产一区二区在线观看av| 蜜桃在线观看..| 免费少妇av软件| 成年av动漫网址| 韩国高清视频一区二区三区| 高清欧美精品videossex| 看免费av毛片| 激情五月婷婷亚洲| 亚洲成国产人片在线观看| 国产在线免费精品| 一级黄片播放器| 又粗又硬又长又爽又黄的视频| 成人免费观看视频高清| 青春草国产在线视频| 国产男人的电影天堂91| 看十八女毛片水多多多| 三级国产精品片| 国产精品三级大全| 夫妻性生交免费视频一级片| 亚洲国产色片| 日本vs欧美在线观看视频| 人妻一区二区av| 日本欧美视频一区| 国国产精品蜜臀av免费| 狠狠婷婷综合久久久久久88av| 欧美成人精品欧美一级黄| 久久久精品94久久精品| 免费观看av网站的网址| 纯流量卡能插随身wifi吗| 搡老乐熟女国产| 国产成人精品久久久久久| 18禁观看日本| 观看美女的网站| 亚洲av国产av综合av卡| 亚洲少妇的诱惑av| 亚洲精品视频女| 夜夜爽夜夜爽视频| av播播在线观看一区| 亚洲精品久久成人aⅴ小说| 日日爽夜夜爽网站| 蜜臀久久99精品久久宅男| 亚洲第一av免费看| 国产精品欧美亚洲77777| videosex国产| 亚洲精品456在线播放app| 国产成人精品婷婷| 看十八女毛片水多多多| 麻豆精品久久久久久蜜桃| 五月玫瑰六月丁香| 精品人妻熟女毛片av久久网站| 欧美精品亚洲一区二区| 中文精品一卡2卡3卡4更新| 少妇的逼水好多| 热99国产精品久久久久久7| 免费播放大片免费观看视频在线观看| 高清视频免费观看一区二区| 秋霞在线观看毛片| 黑人猛操日本美女一级片| 咕卡用的链子| 啦啦啦啦在线视频资源| 综合色丁香网| 国产成人精品久久久久久| 欧美老熟妇乱子伦牲交| 亚洲欧美清纯卡通| 亚洲国产精品专区欧美| 国产视频首页在线观看| 一区二区三区乱码不卡18| 天天躁夜夜躁狠狠躁躁| 国产成人av激情在线播放| 午夜福利网站1000一区二区三区| 免费av不卡在线播放| 免费高清在线观看日韩| www.av在线官网国产| 尾随美女入室| 国产成人一区二区在线| 国产综合精华液| 欧美人与性动交α欧美软件 | 咕卡用的链子| 大码成人一级视频| 寂寞人妻少妇视频99o| 巨乳人妻的诱惑在线观看| 国产精品嫩草影院av在线观看| 丝袜美足系列| 丝袜脚勾引网站| 亚洲av在线观看美女高潮| 伦精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 王馨瑶露胸无遮挡在线观看| 七月丁香在线播放| 制服丝袜香蕉在线| 一区二区日韩欧美中文字幕 | 国产不卡av网站在线观看| 亚洲第一av免费看| 国产淫语在线视频| 精品视频人人做人人爽| 久久久精品免费免费高清| 国产又爽黄色视频| 一边摸一边做爽爽视频免费| 色视频在线一区二区三区| 亚洲av成人精品一二三区| 永久网站在线| 亚洲人成网站在线观看播放| 妹子高潮喷水视频| 91精品国产国语对白视频| 亚洲美女视频黄频| 寂寞人妻少妇视频99o| 亚洲人成网站在线观看播放| 亚洲精品中文字幕在线视频| 精品少妇内射三级| 精品久久蜜臀av无| 国产精品 国内视频| 亚洲精品美女久久久久99蜜臀 | 女人精品久久久久毛片| 欧美精品一区二区大全| 一级爰片在线观看| 少妇熟女欧美另类| 亚洲欧美日韩卡通动漫| 桃花免费在线播放| 精品福利永久在线观看| 久久久久久久久久成人| 欧美成人午夜免费资源| 99久国产av精品国产电影| 97在线视频观看| 免费观看a级毛片全部| 国产男女超爽视频在线观看| 男女无遮挡免费网站观看| av黄色大香蕉| 欧美+日韩+精品| 国产69精品久久久久777片| 丝袜美足系列| 日本猛色少妇xxxxx猛交久久| 久久99蜜桃精品久久| 色网站视频免费| 日韩电影二区| 国产精品蜜桃在线观看| 国内精品宾馆在线| 精品第一国产精品| 久久精品国产自在天天线| 亚洲第一区二区三区不卡| 97超碰精品成人国产| 老熟女久久久| 免费日韩欧美在线观看| 欧美性感艳星| 一级毛片 在线播放| 2022亚洲国产成人精品| 免费av不卡在线播放| 久久毛片免费看一区二区三区| 天天躁夜夜躁狠狠躁躁| 人妻少妇偷人精品九色| 亚洲,一卡二卡三卡| 2018国产大陆天天弄谢| 国产一区二区三区av在线| 欧美日韩一区二区视频在线观看视频在线| 赤兔流量卡办理| 精品人妻一区二区三区麻豆| 少妇 在线观看| 成人无遮挡网站| 久久久久人妻精品一区果冻| 男人操女人黄网站| 亚洲国产精品999| 亚洲欧美色中文字幕在线| 成人二区视频| av女优亚洲男人天堂| 国产1区2区3区精品| 精品一区在线观看国产| 免费黄色在线免费观看| a级片在线免费高清观看视频| av在线app专区| 亚洲av综合色区一区| 老司机影院成人| 男男h啪啪无遮挡| 国产亚洲午夜精品一区二区久久| 人体艺术视频欧美日本| 亚洲欧美色中文字幕在线| 咕卡用的链子| 亚洲 欧美一区二区三区| 多毛熟女@视频| 成人国语在线视频| 最黄视频免费看| 亚洲一区二区三区欧美精品| 久久久久久久久久人人人人人人| 久久av网站| 成年女人在线观看亚洲视频| av线在线观看网站| 亚洲成人一二三区av| 水蜜桃什么品种好| 人妻人人澡人人爽人人| 最近中文字幕高清免费大全6| 国产精品三级大全| 成人国语在线视频| 亚洲欧美成人综合另类久久久| 男女啪啪激烈高潮av片| 高清黄色对白视频在线免费看| 精品久久国产蜜桃| 国产成人a∨麻豆精品| 欧美日韩av久久| 精品亚洲乱码少妇综合久久| 97超碰精品成人国产| 国产熟女欧美一区二区| 波多野结衣一区麻豆| 91国产中文字幕| 蜜臀久久99精品久久宅男| 91午夜精品亚洲一区二区三区| 亚洲国产精品一区二区三区在线| 全区人妻精品视频| 欧美激情 高清一区二区三区| 成人午夜精彩视频在线观看| 国产色婷婷99| 亚洲欧美色中文字幕在线| 国产日韩欧美亚洲二区| 国产深夜福利视频在线观看| 国产精品蜜桃在线观看| 99久久人妻综合| 国产有黄有色有爽视频| 亚洲情色 制服丝袜| 久久久a久久爽久久v久久| 日韩一本色道免费dvd| 日日爽夜夜爽网站| 午夜精品国产一区二区电影| 高清黄色对白视频在线免费看| 亚洲内射少妇av| 欧美亚洲日本最大视频资源| 熟女电影av网| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品久久久久久| 国产麻豆69| 亚洲国产色片| 亚洲精品一二三| 人人妻人人澡人人看| 五月天丁香电影| a级片在线免费高清观看视频| 视频区图区小说| 国产老妇伦熟女老妇高清| 另类亚洲欧美激情| 日本与韩国留学比较| 国产一区亚洲一区在线观看| 欧美成人午夜精品| 欧美+日韩+精品| 最黄视频免费看| 国产精品一区www在线观看| 国产av一区二区精品久久| 国产一区二区在线观看日韩| 日韩欧美精品免费久久| 美女中出高潮动态图| 纵有疾风起免费观看全集完整版| 精品第一国产精品| 精品亚洲成国产av| 免费人妻精品一区二区三区视频| 最近中文字幕高清免费大全6| 丁香六月天网| 亚洲一区二区三区欧美精品| xxxhd国产人妻xxx| 青春草视频在线免费观看| 日本黄大片高清| 亚洲av.av天堂| 老司机影院成人| 免费大片黄手机在线观看| 亚洲综合色惰| 日韩成人av中文字幕在线观看| 天天躁夜夜躁狠狠躁躁| 日韩熟女老妇一区二区性免费视频| 色哟哟·www| 女的被弄到高潮叫床怎么办| 欧美精品人与动牲交sv欧美| 国产日韩一区二区三区精品不卡| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕| 丝袜美足系列| 午夜免费男女啪啪视频观看| 国产成人免费观看mmmm| 国产乱来视频区| 99视频精品全部免费 在线| 少妇精品久久久久久久| 亚洲精品久久久久久婷婷小说| 国产免费一区二区三区四区乱码| 国产又色又爽无遮挡免| 日日撸夜夜添| 在线观看www视频免费| 中文字幕人妻丝袜制服| 天美传媒精品一区二区| 国产免费又黄又爽又色| 欧美人与性动交α欧美精品济南到 | 免费av不卡在线播放| 国产精品一区二区在线不卡| 免费黄频网站在线观看国产| 国产精品国产三级专区第一集| 国产精品久久久久成人av| 在线精品无人区一区二区三| av又黄又爽大尺度在线免费看| 亚洲欧美一区二区三区国产| 五月开心婷婷网| 国产欧美另类精品又又久久亚洲欧美| 亚洲国产精品成人久久小说| 久久久国产精品麻豆| 国产精品欧美亚洲77777| 成人午夜精彩视频在线观看| 亚洲欧洲日产国产| 日本欧美视频一区| www.熟女人妻精品国产 | 亚洲 欧美一区二区三区| 免费大片黄手机在线观看| 啦啦啦视频在线资源免费观看| av有码第一页| 精品久久久精品久久久| 满18在线观看网站| 狂野欧美激情性bbbbbb| 丝袜人妻中文字幕| 国产精品一二三区在线看| 91aial.com中文字幕在线观看| 国产一级毛片在线| 国产精品秋霞免费鲁丝片| 成年av动漫网址| 亚洲国产av影院在线观看| av黄色大香蕉| 少妇的逼好多水| 少妇的逼水好多| 久久影院123| 免费看不卡的av| 亚洲av电影在线观看一区二区三区| 久久久国产欧美日韩av| 有码 亚洲区| 国产不卡av网站在线观看| 久久综合国产亚洲精品| 成年人免费黄色播放视频| 国产成人a∨麻豆精品| 大陆偷拍与自拍| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 国产精品嫩草影院av在线观看| 中文天堂在线官网| 久久精品久久久久久久性| 哪个播放器可以免费观看大片| 如何舔出高潮| 亚洲国产av影院在线观看| 免费黄频网站在线观看国产| 国产亚洲午夜精品一区二区久久| 人妻一区二区av| 国产av国产精品国产| 亚洲av电影在线进入| 99久久精品国产国产毛片| 亚洲av欧美aⅴ国产| 一边摸一边做爽爽视频免费| 欧美日韩精品成人综合77777| 亚洲av在线观看美女高潮| 各种免费的搞黄视频| h视频一区二区三区| 制服人妻中文乱码| 五月开心婷婷网| 亚洲精品国产av成人精品| 边亲边吃奶的免费视频| 久久精品久久久久久久性| 99热国产这里只有精品6| av在线app专区| 久久人人爽人人爽人人片va| 久久亚洲国产成人精品v| 看免费av毛片| 人体艺术视频欧美日本| 国产精品99久久99久久久不卡 | 你懂的网址亚洲精品在线观看| 曰老女人黄片| 午夜av观看不卡| 老熟女久久久| 免费久久久久久久精品成人欧美视频 | 亚洲精品aⅴ在线观看| 韩国av在线不卡| 日韩中字成人| 丝瓜视频免费看黄片| 国产欧美另类精品又又久久亚洲欧美| 久久毛片免费看一区二区三区| 人人妻人人添人人爽欧美一区卜| 国产精品不卡视频一区二区| 一级a做视频免费观看| 高清欧美精品videossex| 国产男女超爽视频在线观看| 免费观看无遮挡的男女| 久久精品国产亚洲av涩爱| 一级片'在线观看视频| 欧美xxxx性猛交bbbb| 99热国产这里只有精品6| 亚洲人成网站在线观看播放| 国产国语露脸激情在线看| 人妻系列 视频| 免费av不卡在线播放| 2018国产大陆天天弄谢| 在线观看一区二区三区激情| 亚洲国产av影院在线观看| 久久久久精品性色| 成人亚洲精品一区在线观看| 免费观看无遮挡的男女| 宅男免费午夜| 久久影院123| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 国产高清三级在线| 日本vs欧美在线观看视频| 国产欧美日韩一区二区三区在线| 国产熟女午夜一区二区三区| 蜜桃国产av成人99| 日本午夜av视频| 婷婷成人精品国产| 日本午夜av视频| 国产欧美日韩一区二区三区在线| 欧美bdsm另类| 香蕉精品网在线| 看非洲黑人一级黄片| 18禁观看日本| 亚洲美女视频黄频| 国产日韩欧美亚洲二区| 18禁裸乳无遮挡动漫免费视频| 丝袜喷水一区| 99久久中文字幕三级久久日本| 亚洲国产精品国产精品| 99久久中文字幕三级久久日本| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 精品视频人人做人人爽| 九色成人免费人妻av| 精品视频人人做人人爽| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级| 国产午夜精品一二区理论片| 国产亚洲午夜精品一区二区久久| 欧美老熟妇乱子伦牲交| 亚洲 欧美一区二区三区| 成年人免费黄色播放视频| 两个人看的免费小视频| 99re6热这里在线精品视频| 黄色毛片三级朝国网站| 国产高清不卡午夜福利| 春色校园在线视频观看| 欧美精品一区二区大全| 黄网站色视频无遮挡免费观看| 亚洲熟女精品中文字幕| 亚洲精品国产av蜜桃| 国产国语露脸激情在线看| 日韩制服丝袜自拍偷拍| 精品一区二区免费观看| 亚洲国产毛片av蜜桃av| 有码 亚洲区| 人妻一区二区av| av国产精品久久久久影院| 国产精品一区二区在线观看99| 大香蕉97超碰在线| 中文字幕另类日韩欧美亚洲嫩草| 国产亚洲精品久久久com| 欧美精品人与动牲交sv欧美| 日韩 亚洲 欧美在线| 永久网站在线| 亚洲国产精品成人久久小说| 男的添女的下面高潮视频| a级毛色黄片| 亚洲熟女精品中文字幕| 男人舔女人的私密视频| 美女主播在线视频| 99久久人妻综合| 亚洲一级一片aⅴ在线观看| 国产又爽黄色视频| 大片电影免费在线观看免费| 久久这里有精品视频免费| 制服人妻中文乱码| 伦理电影大哥的女人| 亚洲精品国产av成人精品| 纵有疾风起免费观看全集完整版| 韩国av在线不卡| 搡女人真爽免费视频火全软件| 天天影视国产精品| 亚洲国产精品一区三区| 日本91视频免费播放| 亚洲欧洲精品一区二区精品久久久 | 乱码一卡2卡4卡精品| 亚洲内射少妇av| 亚洲欧美成人精品一区二区| av国产精品久久久久影院| 99视频精品全部免费 在线| 国产 精品1| 日本av免费视频播放| 精品99又大又爽又粗少妇毛片| 成人国产av品久久久| 天天躁夜夜躁狠狠久久av| 久久精品夜色国产| 哪个播放器可以免费观看大片| 久久精品国产a三级三级三级| 欧美精品国产亚洲| 亚洲性久久影院| 男男h啪啪无遮挡| 免费观看在线日韩| 少妇熟女欧美另类| 日韩不卡一区二区三区视频在线| 久久久a久久爽久久v久久| 免费av中文字幕在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| 亚洲一码二码三码区别大吗| 欧美3d第一页| 日韩制服骚丝袜av| av在线老鸭窝| 国产精品人妻久久久影院| 在线观看一区二区三区激情| 久久久久视频综合| av视频免费观看在线观看| 久久青草综合色| 国产熟女欧美一区二区| 久久精品夜色国产| 观看av在线不卡| 在线观看一区二区三区激情| 成人漫画全彩无遮挡| 免费高清在线观看日韩| 五月天丁香电影| 成人手机av| 亚洲av在线观看美女高潮| 亚洲内射少妇av| 精品少妇黑人巨大在线播放| 欧美日韩国产mv在线观看视频| 日韩人妻精品一区2区三区| 在线亚洲精品国产二区图片欧美| 精品久久蜜臀av无| 深夜精品福利| 中文字幕精品免费在线观看视频 | 国产精品久久久久成人av| 国产精品偷伦视频观看了| videos熟女内射| 黑人猛操日本美女一级片| 亚洲精品美女久久久久99蜜臀 | 精品人妻在线不人妻| 亚洲精品av麻豆狂野| 国产免费一级a男人的天堂|