• <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
    免费av中文字幕在线| 女同久久另类99精品国产91| 天堂8中文在线网| svipshipincom国产片| 成人18禁在线播放| 9色porny在线观看| 久久精品国产99精品国产亚洲性色 | 国产成人系列免费观看| 亚洲午夜精品一区,二区,三区| 91大片在线观看| av福利片在线| cao死你这个sao货| 91麻豆精品激情在线观看国产 | 日韩熟女老妇一区二区性免费视频| 精品一区二区三卡| 欧美国产精品va在线观看不卡| 一级黄色大片毛片| 国产成人啪精品午夜网站| 黄片大片在线免费观看| www.999成人在线观看| 丰满饥渴人妻一区二区三| 国产一区二区激情短视频| 久久国产精品影院| 99国产精品99久久久久| 黄色怎么调成土黄色| 久久国产精品人妻蜜桃| 午夜激情久久久久久久| 高清在线国产一区| 国产成人啪精品午夜网站| 国产亚洲欧美精品永久| 久久精品国产99精品国产亚洲性色 | av国产精品久久久久影院| 国产精品亚洲av一区麻豆| 三级毛片av免费| 首页视频小说图片口味搜索| 日本wwww免费看| 国产麻豆69| 女警被强在线播放| 日韩有码中文字幕| 操出白浆在线播放| 国产精品成人在线| 青草久久国产| 欧美成狂野欧美在线观看| 久久九九热精品免费| 久久久久国内视频| 老鸭窝网址在线观看| 久久亚洲真实| 久久中文字幕人妻熟女| 五月开心婷婷网| 欧美黑人欧美精品刺激| 精品一品国产午夜福利视频| 下体分泌物呈黄色| 午夜激情久久久久久久| 成在线人永久免费视频| 丝袜在线中文字幕| 国精品久久久久久国模美| 亚洲熟女精品中文字幕| 欧美黑人精品巨大| 国产老妇伦熟女老妇高清| 国产区一区二久久| 在线天堂中文资源库| 在线观看66精品国产| 国产精品久久电影中文字幕 | 国产在线精品亚洲第一网站| cao死你这个sao货| 大陆偷拍与自拍| 中文字幕另类日韩欧美亚洲嫩草| 12—13女人毛片做爰片一| 女人爽到高潮嗷嗷叫在线视频| 久久国产亚洲av麻豆专区| 黄片播放在线免费| av天堂久久9| 亚洲精品国产色婷婷电影| 法律面前人人平等表现在哪些方面| 777米奇影视久久| 国产97色在线日韩免费| 别揉我奶头~嗯~啊~动态视频| 在线观看免费日韩欧美大片| 欧美亚洲 丝袜 人妻 在线| 日本vs欧美在线观看视频| 日韩欧美免费精品| 免费观看人在逋| 在线看a的网站| 黄片大片在线免费观看| 十八禁人妻一区二区| 亚洲黑人精品在线| 老鸭窝网址在线观看| 18在线观看网站| 在线观看免费日韩欧美大片| 亚洲一卡2卡3卡4卡5卡精品中文| 精品熟女少妇八av免费久了| 麻豆国产av国片精品| 视频在线观看一区二区三区| a级毛片在线看网站| 午夜福利在线观看吧| 久久精品亚洲av国产电影网| 欧美日韩成人在线一区二区| 亚洲人成电影观看| 久久性视频一级片| 午夜精品久久久久久毛片777| 后天国语完整版免费观看| 在线观看免费午夜福利视频| 视频区图区小说| av福利片在线| 夜夜夜夜夜久久久久| 最新在线观看一区二区三区| 亚洲avbb在线观看| av天堂久久9| 中文字幕av电影在线播放| 中文字幕av电影在线播放| 色视频在线一区二区三区| 久久99热这里只频精品6学生| 亚洲欧美日韩高清在线视频 | 少妇被粗大的猛进出69影院| 人妻一区二区av| 99在线人妻在线中文字幕 | 亚洲精品一二三| 一进一出好大好爽视频| 亚洲av美国av| 19禁男女啪啪无遮挡网站| 美女高潮到喷水免费观看| 日韩视频一区二区在线观看| 精品国产乱码久久久久久男人| 亚洲一码二码三码区别大吗| 欧美在线一区亚洲| 久久性视频一级片| 国产97色在线日韩免费| 午夜激情久久久久久久| 国产精品欧美亚洲77777| 精品一区二区三区av网在线观看 | 精品一区二区三区视频在线观看免费 | 五月天丁香电影| 一区二区三区激情视频| 91字幕亚洲| 精品亚洲乱码少妇综合久久| 亚洲国产成人一精品久久久| 如日韩欧美国产精品一区二区三区| 国产黄频视频在线观看| 久久久国产成人免费| 国产福利在线免费观看视频| 黑人巨大精品欧美一区二区mp4| 一二三四在线观看免费中文在| 国产精品久久久久成人av| 国产精品一区二区免费欧美| 菩萨蛮人人尽说江南好唐韦庄| 麻豆成人av在线观看| 亚洲欧洲精品一区二区精品久久久| 男女边摸边吃奶| 99精品在免费线老司机午夜| 美女扒开内裤让男人捅视频| 欧美日韩亚洲国产一区二区在线观看 | 久久午夜综合久久蜜桃| 啦啦啦视频在线资源免费观看| 欧美日韩国产mv在线观看视频| 一本综合久久免费| 国产精品亚洲一级av第二区| 黑人巨大精品欧美一区二区蜜桃| 精品国产一区二区三区四区第35| 99re6热这里在线精品视频| 人妻一区二区av| 成年人免费黄色播放视频| 免费高清在线观看日韩| 男女边摸边吃奶| 久9热在线精品视频| 狂野欧美激情性xxxx| 亚洲精品国产色婷婷电影| 真人做人爱边吃奶动态| 亚洲欧美一区二区三区久久| 手机成人av网站| 蜜桃国产av成人99| 精品亚洲成国产av| 欧美激情 高清一区二区三区| 露出奶头的视频| 9色porny在线观看| 大片免费播放器 马上看| 激情在线观看视频在线高清 | 欧美人与性动交α欧美软件| 亚洲av片天天在线观看| 国产精品欧美亚洲77777| 国产国语露脸激情在线看| 亚洲一区二区三区欧美精品| 大香蕉久久网| 日本欧美视频一区| 18禁观看日本| 国产精品1区2区在线观看. | 一区福利在线观看| 国产精品久久久久久精品古装| 精品国产国语对白av| 欧美精品av麻豆av| 国产成人免费观看mmmm| 国产熟女午夜一区二区三区| 精品国产一区二区三区久久久樱花| 日韩精品免费视频一区二区三区| 1024香蕉在线观看| 国产成人精品在线电影| 国产精品 国内视频| 精品国产乱码久久久久久小说| 一级毛片电影观看| 国产亚洲精品一区二区www | 亚洲精品成人av观看孕妇| 老熟妇仑乱视频hdxx| 色在线成人网| 国产欧美亚洲国产| 视频在线观看一区二区三区| 午夜老司机福利片| 正在播放国产对白刺激| 女人爽到高潮嗷嗷叫在线视频| 激情视频va一区二区三区| 亚洲成人免费av在线播放| 亚洲精品美女久久久久99蜜臀| 久久精品91无色码中文字幕| 性少妇av在线| 精品国内亚洲2022精品成人 | 美女扒开内裤让男人捅视频| 亚洲国产欧美网| 亚洲专区国产一区二区| 亚洲全国av大片| 一本一本久久a久久精品综合妖精| 在线看a的网站| 亚洲国产av影院在线观看| 最新美女视频免费是黄的| 岛国在线观看网站| 99久久人妻综合| 精品少妇黑人巨大在线播放| 亚洲精品自拍成人| 最新在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 国产无遮挡羞羞视频在线观看| 久久狼人影院| 免费观看av网站的网址| 国产av又大| 多毛熟女@视频| 热99国产精品久久久久久7| 成人国产一区最新在线观看| 老汉色av国产亚洲站长工具| a级毛片在线看网站| 国产精品久久久久久人妻精品电影 | 天天躁狠狠躁夜夜躁狠狠躁| 国产精品麻豆人妻色哟哟久久| 免费看a级黄色片| 久久人人97超碰香蕉20202| 老司机在亚洲福利影院| 水蜜桃什么品种好| av不卡在线播放| 手机成人av网站| 国产亚洲av高清不卡| 在线观看66精品国产| 下体分泌物呈黄色| 中文字幕精品免费在线观看视频| 亚洲国产欧美网| 啦啦啦视频在线资源免费观看| 久久99热这里只频精品6学生| 久久久久精品国产欧美久久久| 中文字幕最新亚洲高清| 国产亚洲欧美精品永久| 在线av久久热| av又黄又爽大尺度在线免费看| 亚洲国产中文字幕在线视频| 操出白浆在线播放| 人人妻人人澡人人看| 91麻豆av在线| 国产在视频线精品| 大型黄色视频在线免费观看| 欧美激情 高清一区二区三区| 亚洲精品自拍成人| 精品久久蜜臀av无| 欧美乱码精品一区二区三区| 成人国产一区最新在线观看| 啦啦啦 在线观看视频| 成年人黄色毛片网站| 久久人人爽av亚洲精品天堂| 成年版毛片免费区| 国产成人系列免费观看| videosex国产| 狠狠婷婷综合久久久久久88av| 久9热在线精品视频| av一本久久久久| 久久国产精品影院| 国产一区有黄有色的免费视频| 99re6热这里在线精品视频| 免费看十八禁软件| 国产亚洲欧美精品永久| 高清欧美精品videossex| 国产精品久久久av美女十八| 国产三级黄色录像| 免费看十八禁软件| 91大片在线观看| 欧美精品亚洲一区二区| 午夜福利乱码中文字幕| 满18在线观看网站| 成人国产av品久久久| 美女视频免费永久观看网站| 精品乱码久久久久久99久播| 高清黄色对白视频在线免费看| 成年女人毛片免费观看观看9 | 在线av久久热| netflix在线观看网站| 在线亚洲精品国产二区图片欧美| 美女高潮喷水抽搐中文字幕| 肉色欧美久久久久久久蜜桃| 久久久精品国产亚洲av高清涩受| 精品高清国产在线一区| 韩国精品一区二区三区| av有码第一页| 精品久久久久久电影网| 热re99久久国产66热| 自线自在国产av| 十八禁网站网址无遮挡| 国产精品成人在线| 久久午夜综合久久蜜桃| 黄色视频在线播放观看不卡| 国产不卡一卡二| 欧美精品人与动牲交sv欧美| 欧美国产精品一级二级三级| 亚洲午夜精品一区,二区,三区| 国产成人系列免费观看| 国产有黄有色有爽视频| 菩萨蛮人人尽说江南好唐韦庄| 热re99久久国产66热| 极品少妇高潮喷水抽搐| 极品教师在线免费播放| 久久精品国产99精品国产亚洲性色 | 99国产综合亚洲精品| 一进一出抽搐动态| 99国产精品99久久久久| 伊人久久大香线蕉亚洲五| 99国产综合亚洲精品| 久久精品国产99精品国产亚洲性色 | 叶爱在线成人免费视频播放| 老司机午夜十八禁免费视频| 中文字幕精品免费在线观看视频| 91成人精品电影| 免费观看av网站的网址| av网站在线播放免费| 欧美日韩精品网址| 1024香蕉在线观看| 亚洲色图 男人天堂 中文字幕| 黄片播放在线免费| 亚洲中文字幕日韩| 蜜桃在线观看..| 十八禁人妻一区二区| 搡老乐熟女国产| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女 | 久久久久久久国产电影| 国产精品.久久久| 首页视频小说图片口味搜索| www日本在线高清视频| 日韩欧美一区视频在线观看| a级毛片在线看网站| 欧美大码av| 欧美日韩亚洲高清精品| 午夜福利免费观看在线| 日本vs欧美在线观看视频| 久久精品亚洲精品国产色婷小说| 久久精品熟女亚洲av麻豆精品| 夜夜爽天天搞| 搡老乐熟女国产| 老熟妇仑乱视频hdxx| 欧美激情高清一区二区三区| 91麻豆av在线| 亚洲三区欧美一区| 久久精品人人爽人人爽视色| 久久这里只有精品19| av电影中文网址| 久久久国产欧美日韩av| 男女边摸边吃奶| 91老司机精品| 男人舔女人的私密视频| 91av网站免费观看| 国产精品秋霞免费鲁丝片| www日本在线高清视频| 久久毛片免费看一区二区三区| 日本av免费视频播放| 精品卡一卡二卡四卡免费| 一个人免费看片子| 老鸭窝网址在线观看| 久久99热这里只频精品6学生| 丰满饥渴人妻一区二区三| 在线观看免费视频网站a站| 亚洲av日韩精品久久久久久密| bbb黄色大片| 9191精品国产免费久久| 天天影视国产精品| 欧美日韩成人在线一区二区| 怎么达到女性高潮| 91成人精品电影| 老鸭窝网址在线观看| 色视频在线一区二区三区| 精品亚洲成a人片在线观看| 最近最新中文字幕大全免费视频| 久久ye,这里只有精品| 国产欧美日韩一区二区三区在线| 大陆偷拍与自拍| 成人18禁在线播放| 免费在线观看影片大全网站| 成人亚洲精品一区在线观看| 午夜激情久久久久久久| 欧美性长视频在线观看| 99国产精品99久久久久| 国产xxxxx性猛交| 黄色片一级片一级黄色片| 亚洲熟女毛片儿| 久久久国产精品麻豆| 高清在线国产一区| 另类精品久久| 80岁老熟妇乱子伦牲交| 久久久久国内视频| 精品第一国产精品| 中文亚洲av片在线观看爽 | 久久久久精品国产欧美久久久| 一二三四在线观看免费中文在| 国产主播在线观看一区二区| 一边摸一边抽搐一进一小说 | 在线观看免费午夜福利视频| 精品国产国语对白av| 在线观看人妻少妇| 高清黄色对白视频在线免费看| 美女扒开内裤让男人捅视频| 99国产精品一区二区三区| 亚洲专区中文字幕在线| 亚洲av美国av| 亚洲av日韩精品久久久久久密| 无限看片的www在线观看| 侵犯人妻中文字幕一二三四区| 国产在线一区二区三区精| 日本wwww免费看| 捣出白浆h1v1| 国产成人av教育| 下体分泌物呈黄色| 狠狠婷婷综合久久久久久88av| 在线观看免费高清a一片| 亚洲精品美女久久av网站| 一级毛片电影观看| avwww免费| 在线av久久热| 日韩欧美国产一区二区入口| 别揉我奶头~嗯~啊~动态视频| 一级a爱视频在线免费观看| 亚洲av日韩在线播放| 午夜91福利影院| 国产欧美日韩精品亚洲av| 国产成人啪精品午夜网站| 日本五十路高清| 最近最新中文字幕大全免费视频| 99国产精品一区二区蜜桃av | 满18在线观看网站| 日本av手机在线免费观看| 无遮挡黄片免费观看| 午夜精品久久久久久毛片777| 国产精品久久久久久精品电影小说| 女人久久www免费人成看片| 男女高潮啪啪啪动态图| 韩国精品一区二区三区| 久久久久久久久久久久大奶| 我的亚洲天堂| 人妻 亚洲 视频| 免费不卡黄色视频| 精品国产亚洲在线| 狠狠婷婷综合久久久久久88av| 国产aⅴ精品一区二区三区波| 丁香六月天网| 免费观看a级毛片全部| 国产精品98久久久久久宅男小说| 亚洲av日韩精品久久久久久密| 欧美精品av麻豆av| 欧美精品高潮呻吟av久久| 成人黄色视频免费在线看| 亚洲精品一卡2卡三卡4卡5卡| 久久国产精品大桥未久av| 国产aⅴ精品一区二区三区波| 啦啦啦免费观看视频1| 亚洲一卡2卡3卡4卡5卡精品中文| 激情在线观看视频在线高清 | 久久久久精品人妻al黑| 成人18禁在线播放| 亚洲少妇的诱惑av| 亚洲精品国产精品久久久不卡| 精品国产乱码久久久久久小说| 国产成人精品久久二区二区免费| 乱人伦中国视频| 国产精品一区二区在线观看99| 人人妻人人澡人人看| 窝窝影院91人妻| 不卡一级毛片| 亚洲成av片中文字幕在线观看| 麻豆国产av国片精品| 啦啦啦中文免费视频观看日本| 黄色 视频免费看| 高清视频免费观看一区二区| 国产欧美日韩一区二区精品| 精品少妇久久久久久888优播| 天天躁夜夜躁狠狠躁躁| 午夜福利在线免费观看网站| 最新在线观看一区二区三区| 欧美性长视频在线观看| 欧美黑人精品巨大| 超色免费av| 亚洲 欧美一区二区三区| 亚洲全国av大片| 视频在线观看一区二区三区| 99国产精品99久久久久| 高清视频免费观看一区二区| 欧美中文综合在线视频| 乱人伦中国视频| 99精品久久久久人妻精品| 欧美成人午夜精品| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| 亚洲中文av在线| 50天的宝宝边吃奶边哭怎么回事| 精品视频人人做人人爽| 国产极品粉嫩免费观看在线| 热re99久久精品国产66热6| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 久久久国产精品麻豆| 看免费av毛片| 日韩欧美一区二区三区在线观看 | 午夜两性在线视频| 国产有黄有色有爽视频| 国产1区2区3区精品| 精品少妇黑人巨大在线播放| 欧美中文综合在线视频| 国产真人三级小视频在线观看| 99热网站在线观看| 嫩草影视91久久| 国产黄色免费在线视频| 在线十欧美十亚洲十日本专区| www.999成人在线观看| svipshipincom国产片| 国产又色又爽无遮挡免费看| 色婷婷av一区二区三区视频| 男人舔女人的私密视频| 如日韩欧美国产精品一区二区三区| 亚洲五月色婷婷综合| 精品福利永久在线观看| 757午夜福利合集在线观看| 老司机影院毛片| 一区二区三区国产精品乱码| 精品乱码久久久久久99久播| 男女高潮啪啪啪动态图| 男女无遮挡免费网站观看| 亚洲国产欧美在线一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品亚洲av国产电影网| 在线观看66精品国产| 亚洲av日韩在线播放| 国产精品久久久久成人av| 男女床上黄色一级片免费看| 成人18禁高潮啪啪吃奶动态图| 国产免费av片在线观看野外av| 亚洲欧洲日产国产| 国产淫语在线视频| 亚洲第一欧美日韩一区二区三区 | 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 精品一品国产午夜福利视频| 婷婷丁香在线五月| 亚洲欧美色中文字幕在线| 精品国产一区二区久久| 伊人久久大香线蕉亚洲五| 中文字幕人妻熟女乱码| 亚洲国产欧美日韩在线播放| 午夜精品久久久久久毛片777| 国产日韩欧美视频二区| 这个男人来自地球电影免费观看| 欧美日韩精品网址| 久久天躁狠狠躁夜夜2o2o| 国产精品免费视频内射| 男女边摸边吃奶| 在线观看66精品国产| 啪啪无遮挡十八禁网站| 下体分泌物呈黄色| 中亚洲国语对白在线视频| 亚洲色图av天堂| 啪啪无遮挡十八禁网站| 成年版毛片免费区| av超薄肉色丝袜交足视频| 97人妻天天添夜夜摸| 亚洲七黄色美女视频| 成人永久免费在线观看视频 | 日日夜夜操网爽| 日韩欧美三级三区| 亚洲欧洲精品一区二区精品久久久| 亚洲成人免费av在线播放| av福利片在线| 国产一区二区三区视频了| 精品福利永久在线观看| 国产在视频线精品| 91精品三级在线观看| 在线播放国产精品三级| 久久久国产欧美日韩av| 黄色片一级片一级黄色片| 日韩一区二区三区影片| 久久久久久久精品吃奶| 一级毛片精品| 蜜桃在线观看..| 欧美日韩视频精品一区| av国产精品久久久久影院| 老汉色∧v一级毛片| 成人亚洲精品一区在线观看| 色视频在线一区二区三区| 精品一区二区三卡| 国产一区二区在线观看av| 亚洲成人免费电影在线观看| 国产亚洲精品久久久久5区| av国产精品久久久久影院| 美女午夜性视频免费| 日本vs欧美在线观看视频| 美女视频免费永久观看网站| 久久天躁狠狠躁夜夜2o2o| 夜夜骑夜夜射夜夜干| 90打野战视频偷拍视频| 50天的宝宝边吃奶边哭怎么回事|