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

    2001—2013年黃土高原植被凈初級生產(chǎn)力時空變化及其歸因

    2017-07-05 02:28:25周夏飛馬國霞曹國志賈倩於方
    安徽農(nóng)業(yè)科學(xué) 2017年14期
    關(guān)鍵詞:黃土高原時空

    周夏飛 馬國霞 曹國志 賈倩 於方

    摘要基于光能利用率模型(CASA),利用遙感數(shù)據(jù)、氣象數(shù)據(jù)和基礎(chǔ)地理數(shù)據(jù),測算了2001—2013年黃土高原植被凈初級生產(chǎn)力(NPP),并輔以一元線性回歸、Hurst指數(shù)及相關(guān)分析等方法,分析了2001—2013年黃土高原NPP時空變化特征、未來變化趨勢及其驅(qū)動因素。結(jié)果表明,2001—2013年黃土高原植被年均NPP呈顯著增加趨勢,年增速為4.9 g/(m2·a)。黃土高原植被NPP空間分布差異顯著,表現(xiàn)出由東南向西北遞減的趨勢。黃土高原植被NPP呈增加趨勢和減少趨勢的面積分別占78.0%和22.0%。Hurst指數(shù)表明研究區(qū)未來植被NPP變化的正向特征顯著,呈持續(xù)性和反持續(xù)性的比重分別為72.1%和28.9%。黃土高原植被NPP變化與降水、氣溫相關(guān)性不大,人類活動是影響植被NPP變化的重要因素,且對NPP有雙重影響。

    關(guān)鍵詞植被凈初級生產(chǎn)力;時空;黃土高原

    中圖分類號Q948.1文獻標(biāo)識碼A文章編號0517-6611(2017)14-0048-06

    AbstractBased on the CarnegieAmesStanfordApproach (CASA), the Net Primary Productivity (NPP) of vegetation in the Loess Plateau during 2001-2013 was measured and calculated with the related remote sensing data, meteorological data and basic geographic data. The temporal and spatial variation characteristics of NPP in the Loess Plateau during 2001-2013, the future trend and its driving factors were analyzed. The results showed that:The average annual NPP of vegetation in the Loess Plateau increased significantly during 2001-2013, with an annual growth rate of 4.9 g/(m2·a).The spatial distribution of vegetation NPP in the Loess Plateau showed a significant trend, with a decreasing tendency from southeast to northwest.The area of improved and degraded NPP accounts for 78.0% and 22.0% of the total study area, respectively.The NPP of the vegetation in the study area was supposed to have a significant positive trend, the area of positive and reverse characteristic accounts for 72.1% and 28.9% of the total study area, respectively.There was no significant correlation between the NPP change and the precipitation as well as the air temperature in the Loess Plateau. The human activities were the core elements affecting the change of the vegetation NPP, leading to both negative effect and positive effect.

    Key wordsNPP;Spatialtemporal;Loess plateau

    植被凈初級生產(chǎn)力(NPP)是指單位時間、單位面積上,植被通過光合作用所產(chǎn)生的有機物總量扣除自養(yǎng)呼吸后的剩余部分[1]。NPP作為地表碳循環(huán)的重要組成部分[2],不僅能反映植被群落在自然環(huán)境條件下的生產(chǎn)能力,表征陸地生態(tài)系統(tǒng)的質(zhì)量狀況[3],還是判定生態(tài)系統(tǒng)碳源、碳匯和調(diào)節(jié)生態(tài)過程的主要因子[4]。目前,全球/區(qū)域碳源、碳匯、碳平衡成為全球變化科學(xué)研究的熱點問題之一[5-7]。因此,動態(tài)監(jiān)測NPP在調(diào)節(jié)全球碳平衡、維護氣候穩(wěn)定及減緩溫室效應(yīng)等全球變化問題中具有重要意義[8]。

    黃土高原地處干旱和濕潤的過渡區(qū),是我國生態(tài)環(huán)境最為脆弱、水土流失最為嚴(yán)重的地區(qū)之一[9]。為改善黃土高原日益惡化的生態(tài)環(huán)境,從1999年開始我國在黃土高原地區(qū)實施退耕還林還草等大規(guī)模植被建設(shè)。1999—2010年黃土高原累計造林面積已達到1 890.6×104 hm2,該地區(qū)植被得到一定恢復(fù)[10-11]。然而,如此大規(guī)模的生態(tài)工程,很可能會改變黃土高原地區(qū)生態(tài)系統(tǒng)景觀格局,從而影響生態(tài)系統(tǒng)特征,而其直接體現(xiàn)在NPP變化上。因此,研究黃土高原自退耕以來NPP的變化具有重要的生態(tài)學(xué)意義[12]。

    目前,已有學(xué)者對黃土高原植被NPP開展了大量研究,并取得了一定成果。謝寶妮等[12]研究表明,2000—2010年黃土高原大部分地區(qū)NPP呈增加趨勢;史曉亮等[10]研究指出,退耕還林還草生態(tài)工程實施之前(1982—1998年),黃土高原大部分區(qū)域植被NPP變化不明顯,自1999年后該區(qū)植被NPP增加趨勢顯著;李登科等[13]基于MOD17A3數(shù)據(jù)集的研究認為,2000—2006年陜西省大部分地區(qū)NPP呈增加趨勢;許紅梅等[14]研究指出,黃土丘陵溝壑區(qū)植被NPP對溫度升高比對降雨變化的響應(yīng)更為敏感。 綜上,以往研究多局限于對植被NPP時空變化特征進行分析,并從整體上探討植被與氣候因子的相關(guān)性,而關(guān)于黃土高原地區(qū)植被NPP變化未來趨勢尚不明確。筆者基于光能利用率模型(CASA),利用遙感數(shù)據(jù)、氣象數(shù)據(jù)和基礎(chǔ)地理數(shù)據(jù)測算2001—2013年黃土高原NPP,并輔以一元線性回歸、Hurst指數(shù)、相關(guān)分析等方法,分析植被NPP時空變化特征、未來趨勢及其驅(qū)動因素,以期為黃土高原生態(tài)環(huán)境保護、建設(shè)及生態(tài)環(huán)境管理提供科學(xué)依據(jù)。

    1資料與方法

    1.1研究區(qū)概況 黃土高原位于我國中部偏北,地理坐標(biāo)為33°43′~41°16′ N,100°54′~114°33′ E,總面積64萬km2。在行政單元上,區(qū)域地跨山西、陜西、甘肅、青海、河南、寧夏以及內(nèi)蒙古等省區(qū)(圖1)。該區(qū)地勢西高東低,平均海拔在1 500~2 000 m,氣候為大陸性季風(fēng)氣候,冬季寒冷干燥、夏季炎熱少雨,年均氣溫3.6~14.3 ℃,年降水量300~800 mm。植被類型主要包括常綠針葉林、落葉闊葉林、灌叢、草地和農(nóng)田,主要土壤類型為黃綿土、草原風(fēng)沙土和褐土。

    1.2數(shù)據(jù)來源

    1.2.1遙感數(shù)據(jù)。遙感數(shù)據(jù)為2001—2013年的MOD13A1數(shù)據(jù),來源于美國國家航空航天局(NASA)的EOS/MODIS數(shù)據(jù)產(chǎn)品(http://e4ftl01.cr.usgs.gov),空間分辨率為500 m×500 m,時間分辨率為16 d。由于MODIS/NDVI產(chǎn)品經(jīng)過水、云、氣溶膠等處理,保證了數(shù)據(jù)質(zhì)量,加之較高的空間分辨率,因此被廣泛應(yīng)用于植被NPP估算的研究中。首先,利用MRT(MODIS Reprojection Tools)將4景MOD13A1數(shù)據(jù)進行合并,并將其投影為Albers(Albers Conical Equal Area);然后,采用最大合成法(MVC),消除異常值的影響,合成月NDVI數(shù)據(jù);最后,利用研究區(qū)邊界裁切NDVI影像,得到2001—2013年逐月NDVI影像。

    1.2.2土地利用類型數(shù)據(jù)。采用的土地利用類型數(shù)據(jù)為2001—2013年的MCD12Q1,來源于美國國家航空航天局(NASA)的EOS/MODIS數(shù)據(jù)產(chǎn)品(http://e4ftl01.cr.usgs.gov),空間分辨率為500 m×500 m。該數(shù)據(jù)集下共有5種土地覆蓋分類方案,采用第1種分類方案即IGBP全球植被分類體系,將其劃分為常綠針葉林、常綠闊葉林、落葉針葉林、落葉闊葉林、混交林等17類。數(shù)據(jù)處理流程同遙感數(shù)據(jù)處理流程一致。

    1.2.3氣候數(shù)據(jù)。采用的氣象數(shù)據(jù)為2001—2013年黃土高原73個氣象站點的逐日數(shù)據(jù)(平均氣溫、降水量、日照時數(shù)),來源于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)。

    1.2.4造林面積數(shù)據(jù)。來源于各省區(qū)統(tǒng)計年鑒和《中國林業(yè)統(tǒng)計年鑒》。

    1.3研究方法

    1.3.1植被NPP估算。采用基于遙感和光能利用率改進的CASA模型[4,15]來估算植被NPP。NPP主要由植被所吸收的光合有效輻射與實際光能利用率2個變量確定,詳細計算過程參考朱文泉等[4,15]的方法。

    1.3.2趨勢分析。 采用一元線性回歸分析方法分析2001—2013年黃土高原13年每個柵格植被NPP的變化趨勢,計算公式為

    式中,Slope為趨勢線的斜率,即NPP的變化趨勢;Yi為柵格第i年的NPP值;n為監(jiān)測時段的年數(shù)(n=13)。Slope>0時,表明NPP呈增加趨勢;Slope<0時,表明NPP呈下降趨勢。

    1.3.3Hurst指數(shù)。 Hurst指數(shù)是預(yù)測時間序列數(shù)據(jù)相對于過去未來發(fā)展趨勢的一個重要指數(shù),目前被廣泛應(yīng)用于水文、氣象等經(jīng)濟領(lǐng)域[16-18]?;谥貥?biāo)極差(R/S)的Hurst指數(shù)是定量描述時間序列信息長期依賴性的有效方法,其基本原理[19]是:

    (1)定義一個時間序列

    若存在R/H∝τH,則說明時間序列{ξ(t)},t = 1,2,…存在 Hurst現(xiàn)象。H值稱為Hurst指數(shù),其值可在雙對數(shù)坐標(biāo)系(lnτ,lnR/S)中用最小二乘法擬合得到。Hurst指數(shù)為0~1.0,02結(jié)果與分析

    2.1植被NPP數(shù)據(jù)驗證 對模擬結(jié)果進行驗證是后期結(jié)果分析的基礎(chǔ)。區(qū)域NPP模型模擬精度驗證方法主要有2種:一種是與實測數(shù)據(jù)進行比較;另一種是與其他模型模擬結(jié)果進行比較。由于黃土高原實測數(shù)據(jù)相對較少,筆者將模型估算的2001—2013年各植被類型平均NPP與其他模型模擬的結(jié)果及部分實測數(shù)據(jù)進行比較,結(jié)果見表1。由表1可知,該研究模擬的平均NPP值與同期的MOD17A3數(shù)據(jù)的值較為接近,基本在實測范圍值內(nèi),與其他模型估算的結(jié)果差距也不大,表明估算結(jié)果具有一定的可靠性。

    2.2植被年均NPP年際變化 由圖2可知,2001—2013年黃土高原植被NPP的年際波動較大,但總體上呈增加趨勢,年增速為4.9 g/(m2·a)(P< 0.05)。植被年均NPP在2001年最低[274.0 g/(m2·a)],2013年最高[378.0 g/(m2·a)],平均值為333.0 g/(m2·a)。值得注意的是,植被NPP在2006和2010年年均NPP下降較大,前者可能是由于2006年降水稀少,導(dǎo)致NPP有所下降,后者可能是由于新退耕還林的植被覆蓋較低所致。

    由圖3可知,2001—2013年NPP > 500 g/(m2·a)的區(qū)域面積百分比整體呈增加趨勢,而NPP < 100 g/(m2·a)的區(qū)域面積百分比整體呈減少趨勢。NPP > 500 g/(m2·a)的植被主要為林地,說明退耕還林還草使得人類對于森林干擾減少,使其得到較好的恢復(fù)。

    2.3植被年均NPP空間變化

    2.3.1空間分布特征。 從空間分布來看(圖4),黃土高原NPP格局呈現(xiàn)由東南向西北遞減的趨勢,高值區(qū)主要分布在陜西省南部、甘肅省南部、青海省、河南省等地區(qū),這主要是由于這些區(qū)域的植被類型主要為針葉林、灌叢,植被長勢也相對較好;低值區(qū)主要分布在內(nèi)蒙古、寧夏南部、甘肅北部等地區(qū),這些地區(qū)多為牧區(qū)和荒漠,植被覆蓋較少、環(huán)境較為惡劣。

    2.3.2空間變化趨勢。 從空間變化上看(圖5),黃土高原NPP整體上呈增加趨勢,呈增加趨勢和減小趨勢的面積分別占總面積的78.0%和22.0%,其中28.0%的區(qū)域呈顯著增加趨勢,主要分布在陜西北部、山西北部、寧夏南部、青海南部等區(qū)域,這可能是由于這些地區(qū)原本的植被覆蓋度相對較低,然而在“退耕還林還草”等工程實施下,該地區(qū)NPP呈增加趨勢;黃土高原僅有4.0%的區(qū)域呈顯著下降趨勢,主要分布在內(nèi)蒙古西部、寧夏北部,這主要是由于這些地區(qū)城市化進程速度加快,占用了大量林地、草地,導(dǎo)致其NPP下降。

    2001-20132.4植被NPP未來趨勢黃土高原植被NPP的Hurst指數(shù)平均值為0.551 8,其中,Hurst指數(shù)大于0.5的像元數(shù)占72.1%,說明植被NPP變化的正向特征顯著;而Hurst指數(shù)小于0.5的像元數(shù)占28.9%。Hurst指數(shù)高值區(qū)主要分布在內(nèi)蒙古、寧夏西部、甘肅省北部、青海省等地區(qū),說明這些地區(qū)植被NPP變化具有同向特征,即未來變化趨勢與過去變化趨勢一致;而Hurst指數(shù)低值區(qū)主要分布在陜西省西部及南部、寧夏東部及南部、山西省等地區(qū),說明這些地區(qū)植被NPP變化具有反向特征,即未來變化趨勢與過去變化趨勢相反(圖6)。

    將Hurst指數(shù)與一元線性回歸分析結(jié)果進行疊加分析發(fā)現(xiàn),植被NPP持續(xù)改善地區(qū)主要分布在內(nèi)蒙古南部、陜西北部、甘肅北部、寧夏南部等地區(qū),這主要是由于這些區(qū)域長期實行退耕還林工程,已取得明顯效果;植被NPP持續(xù)退化地區(qū)主要分布在陜西南部、內(nèi)蒙古北部、青海北部等地區(qū),這可能是由于這些區(qū)域未開展退耕還林工程或工程措施不到位;植被NPP由改善轉(zhuǎn)退化地區(qū)主要分布在陜西中部、甘肅南部、山西等地區(qū),這可能是由于退耕還林工程措施不到位;植被NPP由退化轉(zhuǎn)改善地區(qū)主要分布在青海省西部、內(nèi)蒙古西部、山西東部等地區(qū),這可能是由于這些區(qū)域之前未實施退耕還林工程,近年來開始開展退耕還林工程所致(圖6)。

    2.5黃土高原植被NPP影響因素

    2.5.1植被NPP與氣候因子的相關(guān)性。 降水和氣溫等氣候因子的變化對陸地植被的生長發(fā)育具有重要影響,因此筆者利用逐像元的相關(guān)分析方法,分別對2001—2013年的年植被NPP、年降水量和年平均氣溫的相關(guān)系數(shù)進行分析,從而在空間尺度分析植被NPP與主要氣候因子的關(guān)系。從圖7可以看出,黃土高原植被NPP與年降雨量和年均溫具有一定的相關(guān)性,但相關(guān)性不大。在黃土高原,73.3%的區(qū)域年降雨量與NPP呈正相關(guān),其中僅有6.9%的區(qū)域通過了0.05顯著性檢驗,主要分布在黃土高原東部及西部地區(qū),表明降水量增加,NPP增加;60.8%的區(qū)域年均溫與NPP呈負相關(guān),其中僅有9.9%的區(qū)域通過了0.05顯著性檢驗,主要分布在黃土高原北部,這些區(qū)域NPP隨著溫度的升高反而減?。▓D7)。

    2.5.2人類活動的影響。除氣溫、降水等自然因素外,人類活動也會對黃土高原NPP的時空分布格局產(chǎn)生重要影響,人類活動對其影響包括正面影響(如退耕還林)和負面影響(如城市擴張)。筆者主要分析了黃土高原退耕還林工程的歷年累積造林面積與植被NPP變化情況,結(jié)果見圖8。由圖8可知,自1999年退耕還林還草工程實施以來,黃土高原造林面積持續(xù)上升,植被NPP總體上也表現(xiàn)出升高趨勢,表明退耕還林還草工程可以在一定程度上提高植被NPP,有效改善黃土高原的生態(tài)環(huán)境;但NPP變化并未與造林面積變化嚴(yán)格一致,存在一定的波動性,尤其是2010年波動較大,這可能是由于近年來城市化、放牧等活動對植被NPP造成負面影響。

    3結(jié)論與討論

    該研究以黃土高原為研究區(qū)域,基于光能利用率模型(CASA),利用遙感數(shù)據(jù)、氣象數(shù)據(jù)和基礎(chǔ)地理數(shù)據(jù)測算了2001—2013年黃土高原植被凈初級生產(chǎn)力(NPP),并輔以一元線性回歸、Hurst指數(shù)、偏相關(guān)分析等方法,分別從3個尺度分析了植被NPP時空變化特征、未來趨勢及驅(qū)動因素,得到以下結(jié)論:

    (1)2001—2013年黃土高原植被年均NPP呈增加趨勢,年增速為4.9 g/(m2·a),植被年均NPP在2001年最低[274.0 g/(m2·a)],2013年最高[378.0 g/(m2·a)],平均值為333.0 g/(m2·a)。

    (2)黃土高原NPP格局呈現(xiàn)由東南向西北遞減的趨勢,高值區(qū)主要分布在陜西省南部、甘肅省南部、青海省、河南省等地區(qū),低值區(qū)主要分布在內(nèi)蒙古、寧夏、甘肅北部等地區(qū)。

    (3)黃土高原NPP整體上呈增加趨勢,呈增加趨勢和減小趨勢的面積分別占總面積的78.0%和22.0%,其中28.0%的區(qū)域呈顯著增加趨勢,主要分布在陜西北部、山西北部、寧夏南部、青海南部等區(qū)域。

    (4)Hurst指數(shù)分析表明,黃土高原植被NPP變化正向特征顯著,植被NPP由改善轉(zhuǎn)退化的地區(qū)主要分布在陜西中部、甘肅南部、山西等;植被NPP由退化轉(zhuǎn)改善的地區(qū)主要分布在青海西部、內(nèi)蒙古西部、山西東部等。

    (5)黃土高原植被NPP與降水量及溫度相關(guān)性不大,73.3%的區(qū)域年降雨量與NPP呈正相關(guān),其中僅有6.9%的區(qū)域通過了0.05顯著性檢驗;60.8%的區(qū)域年均溫與NPP呈負相關(guān),其中僅有9.9%的區(qū)域通過了0.05顯著性檢驗。

    (6)植被NPP變化是氣象變化和人類活動共同作用的結(jié)果,區(qū)域氣候變化不大背景下,隨著“退耕還林”等一系列生態(tài)保護工程的實施,黃土高原生態(tài)環(huán)境將進一步好轉(zhuǎn)。

    由于受遙感數(shù)據(jù)(MODIS數(shù)據(jù))時間分辨率的影響,該研究未分析黃土高原大規(guī)模生態(tài)工程實施前后的植被NPP時空演變。采用逐像元線性回歸模型融合多種數(shù)據(jù)源構(gòu)建長時間序列遙感數(shù)據(jù),并進行植被NPP估算,深入分析黃土高原退耕還林前后NPP變化及原因是未來研究的重點。

    參考文獻

    [1] FIELD C B,RANDERSON J T,MALMSTRM C M.Global net primary production:Combining ecology and remote sensing[J].Remote sensing of environment,1995,51(1):74-88.

    [2] PIAO S L,CIAIS P,LOMAS M,et al.Contribution of climate change and rising CO2 to terrestrial carbon balance in East Asia:A multimodel analysis[J].Global & planetary change,2011,75(3/4):133-142.

    [3] FIELD C B,BEHRENFELD M J,RANDERSON J T,et al.Primary production of the biosphere:Integrating terrestrial and oceanic components[J].Science,1998,281(5374):237-240.

    [4] 朱文泉,張錦水,潘耀忠,等.中國陸地生態(tài)系統(tǒng)生態(tài)資產(chǎn)測量及其動態(tài)變化分析[J].應(yīng)用生態(tài)學(xué)報,2007,18(3):586-594.

    [5] 方精云,郭兆迪,樸世龍,等.1981~2000年中國陸地植被碳匯的估算[J].中國科學(xué)(D輯:地球科學(xué)),2007,37(6):804-812.

    [6] LOBELL D B,HICKE J A,ASNER G P,et al.Satellite estimates of productivity and light use efficiency in United States agriculture,1982-98[J].Global change biology,2002,8(8):722-735.

    [7] 李晶,任志遠.基于GIS的陜北黃土高原土地生態(tài)系統(tǒng)固碳釋氧價值評價[J].中國農(nóng)業(yè)科學(xué),2011,44(14):2943-2950.

    [8] 張巖,韋振鋒,黃毅.1999-2012年三江源區(qū)植被凈初級生產(chǎn)力及固碳釋氧量測評[J].水土保持通報,2016,36(1):100-105.

    [9] 趙安周,劉憲鋒,朱秀芳,等.2000~2014年黃土高原植被覆蓋時空變化特征及其歸因[J].中國環(huán)境科學(xué),2016,36(5):1568-1578.

    [10] 史曉亮,楊志勇,王馨爽,等.黃土高原植被凈初級生產(chǎn)力的時空變化及其與氣候因子的關(guān)系[J].中國農(nóng)業(yè)氣象,2016,37(4):445-453.

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

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

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

    [14] 許紅梅,高清竹,黃永梅,等.氣候變化對黃土丘陵溝壑區(qū)植被凈第一性生產(chǎn)力的影響模擬[J].生態(tài)學(xué)報,2006,26(9):2939-2947.

    [15] 朱文泉,潘耀忠,何浩,等.中國典型植被最大光利用率模擬[J].科學(xué)通報,2006,51(6):700-706.

    [16] 張翀,任志遠.黃土高原地區(qū)植被覆蓋變化的時空差異及未來趨勢[J].資源科學(xué),2011,33(11):2143-2149.

    [17] 李雙雙,延軍平,萬佳.近10年陜甘寧黃土高原區(qū)植被覆蓋時空變化特征[J].地理學(xué)報,2012,67(7):960-970.

    [18] 劉憲鋒,任志遠,林志慧,等.2000-2011年三江源區(qū)植被覆蓋時空變化特征[J].地理學(xué)報,2013,68(7):897-908.

    [19] 江田漢,鄧蓮堂.Hurst指數(shù)估計中存在的若干問題:以在氣候變化研究中的應(yīng)用為例[J].地理科學(xué),2004,24(2):177-182.

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

    [21] 陶波.中國陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力與凈生態(tài)系統(tǒng)生產(chǎn)力模擬研究[D].北京:中國科學(xué)院地理科學(xué)與資源研究所,2003.

    [22] NI J,ZHANG X S,SCURLOCK J M O.Synthesis and analysis of biomass and net primary productivity in Chinese forests[J].Annals of forest science,2001,58(4):351-384.

    猜你喜歡
    黃土高原時空
    跨越時空的相遇
    跨越時空找自己
    童話世界(2020年31期)2021-01-04 01:41:58
    鏡中的時空穿梭
    玩一次時空大“穿越”
    “天紅2號”蘋果在豫西黃土高原的引種表現(xiàn)
    選舉 沸騰了黃土高原(下)
    公民與法治(2016年3期)2016-05-17 04:09:00
    選舉沸騰了黃土高原(上)
    公民與法治(2016年1期)2016-05-17 04:07:56
    時空之門
    灑向黃土高原的愛
    中國火炬(2015年7期)2015-07-31 17:39:57
    跨時空,再相逢
    車迷(2015年6期)2015-03-20 02:43:48
    久久精品影院6| 欧美日韩精品成人综合77777| 日本黄色视频三级网站网址| 国内精品美女久久久久久| 一进一出抽搐gif免费好疼| 国内精品美女久久久久久| 日韩av不卡免费在线播放| 在线观看av片永久免费下载| 22中文网久久字幕| 男女做爰动态图高潮gif福利片| 哪里可以看免费的av片| 男女之事视频高清在线观看| 免费观看的影片在线观看| 免费在线观看成人毛片| 亚洲不卡免费看| 久久精品国产鲁丝片午夜精品| 欧美最黄视频在线播放免费| 在线观看66精品国产| 欧美三级亚洲精品| 国产精品人妻久久久久久| 久久久久久久亚洲中文字幕| 免费av观看视频| 两个人的视频大全免费| 亚洲18禁久久av| 俺也久久电影网| 欧美国产日韩亚洲一区| 黄色日韩在线| 丰满人妻一区二区三区视频av| 国产精品久久久久久av不卡| 女的被弄到高潮叫床怎么办| 欧美最新免费一区二区三区| 亚洲欧美成人综合另类久久久 | 久久人妻av系列| .国产精品久久| 国产一区二区在线观看日韩| 亚洲国产精品成人综合色| 久久精品国产99精品国产亚洲性色| 中国美白少妇内射xxxbb| 91午夜精品亚洲一区二区三区| 少妇的逼好多水| 一本一本综合久久| 国产精品野战在线观看| 日本五十路高清| 欧美性猛交黑人性爽| 国产片特级美女逼逼视频| 亚洲av成人精品一区久久| 中文亚洲av片在线观看爽| 久久久久久久久大av| 久久天躁狠狠躁夜夜2o2o| 欧美高清性xxxxhd video| 日本熟妇午夜| 男人舔女人下体高潮全视频| 美女内射精品一级片tv| 亚洲人成网站高清观看| 久久热精品热| 偷拍熟女少妇极品色| 在线看三级毛片| 国产精品人妻久久久久久| 欧美又色又爽又黄视频| 日本免费a在线| 麻豆av噜噜一区二区三区| 国产亚洲精品av在线| 又粗又爽又猛毛片免费看| 97在线视频观看| 久久精品夜色国产| 免费看光身美女| 精华霜和精华液先用哪个| 精品免费久久久久久久清纯| 人妻久久中文字幕网| 又黄又爽又免费观看的视频| 99热只有精品国产| 国产精品一区www在线观看| 亚洲三级黄色毛片| 亚洲欧美日韩东京热| 男女啪啪激烈高潮av片| 亚洲国产欧洲综合997久久,| 国产精品久久久久久av不卡| 国产色婷婷99| 中文字幕av成人在线电影| 欧美一区二区国产精品久久精品| 能在线免费观看的黄片| a级一级毛片免费在线观看| 久久热精品热| 亚洲七黄色美女视频| 精品午夜福利视频在线观看一区| 在线观看免费视频日本深夜| 色哟哟·www| 欧美xxxx性猛交bbbb| 午夜日韩欧美国产| 国产男靠女视频免费网站| 欧美性感艳星| 亚洲第一电影网av| 麻豆久久精品国产亚洲av| 五月玫瑰六月丁香| 少妇熟女欧美另类| 大香蕉久久网| 国产国拍精品亚洲av在线观看| 在线免费十八禁| 免费一级毛片在线播放高清视频| 美女大奶头视频| 久久6这里有精品| 亚洲国产精品久久男人天堂| 三级国产精品欧美在线观看| 国产精品女同一区二区软件| 老司机午夜福利在线观看视频| 久久久久国产网址| 五月玫瑰六月丁香| 干丝袜人妻中文字幕| 日韩,欧美,国产一区二区三区 | 神马国产精品三级电影在线观看| 日本免费a在线| 色av中文字幕| 精品一区二区三区人妻视频| 一级毛片久久久久久久久女| 老女人水多毛片| 国产一级毛片七仙女欲春2| 亚洲av成人av| 欧美中文日本在线观看视频| 小说图片视频综合网站| 欧美一区二区国产精品久久精品| 成年av动漫网址| 国产视频一区二区在线看| 欧美精品国产亚洲| 中文字幕精品亚洲无线码一区| 欧美日本视频| 国产成人91sexporn| 毛片女人毛片| 熟妇人妻久久中文字幕3abv| 精品乱码久久久久久99久播| 干丝袜人妻中文字幕| 老熟妇乱子伦视频在线观看| 国产成人福利小说| 亚洲人成网站在线播| 高清午夜精品一区二区三区 | 最近的中文字幕免费完整| 91狼人影院| 男人舔女人下体高潮全视频| 简卡轻食公司| 少妇高潮的动态图| 亚洲内射少妇av| 欧美丝袜亚洲另类| 久久鲁丝午夜福利片| www.色视频.com| 亚洲四区av| 亚洲精华国产精华液的使用体验 | 69av精品久久久久久| 精品久久久久久久久久久久久| 国内精品美女久久久久久| 九九热线精品视视频播放| 国产v大片淫在线免费观看| 在线观看免费视频日本深夜| 欧美性感艳星| 内地一区二区视频在线| 此物有八面人人有两片| 校园人妻丝袜中文字幕| 亚洲在线观看片| 婷婷精品国产亚洲av在线| 国产黄色小视频在线观看| 亚洲国产精品久久男人天堂| 人人妻,人人澡人人爽秒播| 日本爱情动作片www.在线观看 | 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| 熟女电影av网| 噜噜噜噜噜久久久久久91| 女人被狂操c到高潮| 亚洲专区国产一区二区| 中国国产av一级| 校园春色视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产探花在线观看一区二区| 人人妻,人人澡人人爽秒播| 91在线精品国自产拍蜜月| 日本撒尿小便嘘嘘汇集6| 女人十人毛片免费观看3o分钟| 日本黄色片子视频| 精品一区二区三区人妻视频| 91久久精品国产一区二区成人| 午夜精品一区二区三区免费看| 国产精品久久久久久亚洲av鲁大| 一个人观看的视频www高清免费观看| 淫妇啪啪啪对白视频| 亚洲精品一区av在线观看| 欧美成人a在线观看| 亚洲丝袜综合中文字幕| 日韩制服骚丝袜av| 久久人人爽人人片av| 波野结衣二区三区在线| 美女xxoo啪啪120秒动态图| 高清毛片免费看| 人妻夜夜爽99麻豆av| 日本免费a在线| 观看免费一级毛片| 18禁黄网站禁片免费观看直播| 亚洲精品国产av成人精品 | 久久韩国三级中文字幕| 免费av观看视频| 久久久成人免费电影| 能在线免费观看的黄片| 国产亚洲精品久久久久久毛片| 97在线视频观看| 国产91av在线免费观看| 最好的美女福利视频网| 日韩亚洲欧美综合| 黄色配什么色好看| 综合色av麻豆| 国产大屁股一区二区在线视频| 男女啪啪激烈高潮av片| 99久久精品国产国产毛片| 国产伦在线观看视频一区| 久久精品国产清高在天天线| 国产v大片淫在线免费观看| 亚洲一区二区三区色噜噜| 日韩欧美免费精品| 可以在线观看毛片的网站| eeuss影院久久| 1024手机看黄色片| 成年版毛片免费区| 深爱激情五月婷婷| 欧美一级a爱片免费观看看| 午夜福利18| 久久九九热精品免费| 日韩国内少妇激情av| 久久久久久久久久黄片| 日本三级黄在线观看| 日韩精品青青久久久久久| 夜夜爽天天搞| 特级一级黄色大片| 久久精品综合一区二区三区| 日本一本二区三区精品| 九九久久精品国产亚洲av麻豆| 久久精品夜色国产| 精品一区二区三区人妻视频| 成人无遮挡网站| 亚洲无线观看免费| 有码 亚洲区| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| 好男人在线观看高清免费视频| 亚洲最大成人中文| 欧美性猛交╳xxx乱大交人| 国产高清三级在线| 婷婷六月久久综合丁香| 啦啦啦啦在线视频资源| 91午夜精品亚洲一区二区三区| 在线观看午夜福利视频| 黄色配什么色好看| 偷拍熟女少妇极品色| 亚洲av成人精品一区久久| av国产免费在线观看| 国产淫片久久久久久久久| 一个人看的www免费观看视频| 综合色av麻豆| 免费人成在线观看视频色| 热99在线观看视频| 美女高潮的动态| 国产免费一级a男人的天堂| 日日摸夜夜添夜夜爱| 中文字幕av在线有码专区| 国产精品久久久久久精品电影| 韩国av在线不卡| 国产极品精品免费视频能看的| 亚洲国产精品成人综合色| 97超碰精品成人国产| 亚洲一级一片aⅴ在线观看| 中出人妻视频一区二区| 国产乱人视频| 最近在线观看免费完整版| 人人妻人人澡人人爽人人夜夜 | 成人欧美大片| 国产一区亚洲一区在线观看| 又爽又黄a免费视频| 国产午夜精品久久久久久一区二区三区 | 亚洲国产精品久久男人天堂| eeuss影院久久| 免费在线观看成人毛片| 久久热精品热| 成人性生交大片免费视频hd| 日本黄大片高清| 国产精品久久视频播放| 特级一级黄色大片| 偷拍熟女少妇极品色| 成人综合一区亚洲| 国产精品久久久久久精品电影| 一级黄片播放器| av.在线天堂| 搡老妇女老女人老熟妇| 国产精品亚洲美女久久久| 日本精品一区二区三区蜜桃| 亚洲国产精品久久男人天堂| 人人妻人人澡欧美一区二区| 精品人妻熟女av久视频| 黄色日韩在线| 日日干狠狠操夜夜爽| 亚洲成av人片在线播放无| 精品午夜福利在线看| 内地一区二区视频在线| 最新中文字幕久久久久| 午夜福利成人在线免费观看| 亚洲精品成人久久久久久| 在线看三级毛片| 国产高潮美女av| 午夜免费激情av| 国产黄片美女视频| 毛片一级片免费看久久久久| 久久鲁丝午夜福利片| 国产精品精品国产色婷婷| 免费人成在线观看视频色| 国产精品,欧美在线| 99热这里只有精品一区| 男女下面进入的视频免费午夜| 亚洲成a人片在线一区二区| 国产高清视频在线播放一区| 在线免费十八禁| 久久精品91蜜桃| 亚洲高清免费不卡视频| 国产精品1区2区在线观看.| 国产高清有码在线观看视频| 免费看光身美女| av在线天堂中文字幕| a级一级毛片免费在线观看| 国产高清视频在线观看网站| 午夜精品一区二区三区免费看| 久久精品国产清高在天天线| 一本精品99久久精品77| 男女下面进入的视频免费午夜| 男人狂女人下面高潮的视频| 亚洲欧美日韩东京热| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 午夜精品一区二区三区免费看| 久久精品国产清高在天天线| 国产 一区 欧美 日韩| 精品国产三级普通话版| 亚洲图色成人| 国产精品久久久久久精品电影| 中文字幕免费在线视频6| 午夜影院日韩av| 韩国av在线不卡| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲网站| 人人妻,人人澡人人爽秒播| 午夜福利在线观看吧| 在线国产一区二区在线| 日韩人妻高清精品专区| 久久人人爽人人爽人人片va| 国产片特级美女逼逼视频| 精品一区二区三区av网在线观看| 日日啪夜夜撸| 69av精品久久久久久| 级片在线观看| 日本免费一区二区三区高清不卡| 美女xxoo啪啪120秒动态图| 精品久久久久久成人av| av天堂中文字幕网| 精品久久久久久久久av| 内地一区二区视频在线| 成年女人看的毛片在线观看| 精品久久久久久久久av| 亚洲图色成人| 大型黄色视频在线免费观看| 嫩草影院精品99| 麻豆av噜噜一区二区三区| 赤兔流量卡办理| 成人av一区二区三区在线看| 在线观看美女被高潮喷水网站| 乱码一卡2卡4卡精品| 亚洲av中文av极速乱| 久久综合国产亚洲精品| 免费黄网站久久成人精品| 高清午夜精品一区二区三区 | 女同久久另类99精品国产91| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 一级毛片aaaaaa免费看小| 国产高清激情床上av| 狂野欧美白嫩少妇大欣赏| 网址你懂的国产日韩在线| 亚洲国产欧美人成| 欧美国产日韩亚洲一区| av在线观看视频网站免费| 给我免费播放毛片高清在线观看| 国产探花极品一区二区| 精品午夜福利在线看| 午夜精品国产一区二区电影 | 午夜精品国产一区二区电影 | 欧美+亚洲+日韩+国产| 国产精品一区www在线观看| 久久热精品热| 久久99热6这里只有精品| 村上凉子中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 一区福利在线观看| 亚洲天堂国产精品一区在线| 亚洲国产欧洲综合997久久,| 国产精品电影一区二区三区| 国产精品爽爽va在线观看网站| 日韩成人伦理影院| av视频在线观看入口| 大香蕉久久网| 国产成人影院久久av| 亚洲第一区二区三区不卡| 国产欧美日韩一区二区精品| 男女啪啪激烈高潮av片| 免费看av在线观看网站| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 国产精品久久久久久精品电影| 老女人水多毛片| 亚洲熟妇中文字幕五十中出| 看片在线看免费视频| 中文字幕熟女人妻在线| 久久鲁丝午夜福利片| 国产成人a∨麻豆精品| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 日韩av不卡免费在线播放| 欧美色欧美亚洲另类二区| 无遮挡黄片免费观看| 干丝袜人妻中文字幕| 精品久久久久久久久亚洲| 午夜精品国产一区二区电影 | 久久精品综合一区二区三区| 国产精品99久久久久久久久| 六月丁香七月| 成人特级黄色片久久久久久久| 国产亚洲av嫩草精品影院| 色哟哟哟哟哟哟| 能在线免费观看的黄片| 日韩精品青青久久久久久| 啦啦啦啦在线视频资源| 在线观看一区二区三区| 国产精品一区www在线观看| 免费av不卡在线播放| 亚洲av熟女| 91在线观看av| 国产国拍精品亚洲av在线观看| 人人妻人人看人人澡| 深爱激情五月婷婷| 大香蕉久久网| 欧美成人a在线观看| 日本黄色视频三级网站网址| 51国产日韩欧美| 国产精品久久久久久久电影| 最近在线观看免费完整版| 国产三级中文精品| 欧美成人免费av一区二区三区| www.色视频.com| 日韩大尺度精品在线看网址| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 国产成人a∨麻豆精品| av.在线天堂| 久久精品夜色国产| 国产 一区精品| 欧美高清性xxxxhd video| 久久久久久久久久黄片| 18+在线观看网站| 搡老妇女老女人老熟妇| 插逼视频在线观看| 午夜爱爱视频在线播放| 亚洲国产精品sss在线观看| 国产精品久久电影中文字幕| 十八禁网站免费在线| 91在线观看av| 六月丁香七月| 久久久久久九九精品二区国产| 欧美3d第一页| 香蕉av资源在线| 美女被艹到高潮喷水动态| 超碰av人人做人人爽久久| 国产精品一区二区三区四区免费观看 | 国产精品电影一区二区三区| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| 亚洲精品影视一区二区三区av| 精品福利观看| 嫩草影院精品99| av专区在线播放| 长腿黑丝高跟| 精品久久久久久久人妻蜜臀av| 国国产精品蜜臀av免费| 日韩欧美在线乱码| 波多野结衣高清作品| 亚洲av电影不卡..在线观看| 国国产精品蜜臀av免费| 欧美成人免费av一区二区三区| 午夜福利高清视频| 中文字幕熟女人妻在线| 国产aⅴ精品一区二区三区波| 日产精品乱码卡一卡2卡三| 精品欧美国产一区二区三| 亚洲色图av天堂| a级毛片免费高清观看在线播放| 色视频www国产| av.在线天堂| 亚洲av第一区精品v没综合| 真实男女啪啪啪动态图| 日本 av在线| 麻豆av噜噜一区二区三区| 久久精品人妻少妇| 亚洲人与动物交配视频| 99久久久亚洲精品蜜臀av| 国产成人影院久久av| 不卡视频在线观看欧美| 成年女人永久免费观看视频| 欧美日本视频| 69av精品久久久久久| 三级男女做爰猛烈吃奶摸视频| 超碰av人人做人人爽久久| 欧美+日韩+精品| 国产国拍精品亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 日韩制服骚丝袜av| 国内久久婷婷六月综合欲色啪| 男女下面进入的视频免费午夜| 黄色欧美视频在线观看| 日韩成人av中文字幕在线观看 | 在现免费观看毛片| 97人妻精品一区二区三区麻豆| 久久久a久久爽久久v久久| 成人漫画全彩无遮挡| 亚洲成人av在线免费| 在线免费观看不下载黄p国产| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 丝袜美腿在线中文| 久久久久久九九精品二区国产| 精品日产1卡2卡| 亚洲18禁久久av| 一个人看视频在线观看www免费| 国产亚洲精品久久久com| 国产男靠女视频免费网站| 亚洲天堂国产精品一区在线| 俄罗斯特黄特色一大片| 久久精品91蜜桃| 久久99热6这里只有精品| 国产精品,欧美在线| 亚洲在线观看片| 亚洲av第一区精品v没综合| 国内精品美女久久久久久| 蜜桃亚洲精品一区二区三区| 波多野结衣高清作品| 亚洲欧美日韩东京热| 麻豆成人午夜福利视频| 午夜a级毛片| 波多野结衣巨乳人妻| 久久天躁狠狠躁夜夜2o2o| 久久久午夜欧美精品| 日韩欧美国产在线观看| 久久综合国产亚洲精品| 精品不卡国产一区二区三区| 99热全是精品| 菩萨蛮人人尽说江南好唐韦庄 | 秋霞在线观看毛片| 久久精品国产清高在天天线| 国产成人一区二区在线| 久久国内精品自在自线图片| 日韩欧美精品v在线| 黄色欧美视频在线观看| 日日摸夜夜添夜夜爱| 中文字幕av成人在线电影| 免费黄网站久久成人精品| 亚洲无线在线观看| 欧美区成人在线视频| aaaaa片日本免费| 日日撸夜夜添| 国产伦一二天堂av在线观看| 午夜精品在线福利| 久久欧美精品欧美久久欧美| 九九爱精品视频在线观看| 噜噜噜噜噜久久久久久91| 成人亚洲欧美一区二区av| 成年免费大片在线观看| 国产男人的电影天堂91| 成年版毛片免费区| 最近的中文字幕免费完整| 国产私拍福利视频在线观看| 久久精品国产亚洲av涩爱 | 中文字幕人妻熟人妻熟丝袜美| 免费看美女性在线毛片视频| 在线观看66精品国产| 欧美色视频一区免费| 国产午夜精品久久久久久一区二区三区 | 99热精品在线国产| 欧美zozozo另类| 久久精品91蜜桃| 午夜激情福利司机影院| eeuss影院久久| 中文字幕av在线有码专区| 最近手机中文字幕大全| 精品久久久久久久久亚洲| 欧美最黄视频在线播放免费| 国产精品爽爽va在线观看网站| 亚洲中文日韩欧美视频| 成人漫画全彩无遮挡| 国产探花极品一区二区| 嫩草影院入口| 国产熟女欧美一区二区| 好男人在线观看高清免费视频| 国产乱人视频| 久久99热6这里只有精品| 日本 av在线| 精品久久久久久久久久免费视频| 丰满的人妻完整版| 亚洲在线观看片| 俄罗斯特黄特色一大片| 老熟妇乱子伦视频在线观看| 亚洲四区av| 色在线成人网| 最近最新中文字幕大全电影3| 国产一区二区激情短视频|