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

    內(nèi)蒙古草原植被最大光能利用率取值優(yōu)化研究

    2016-07-12 12:56:08辛?xí)云?/span>包玉海王牧蘭元志輝烏蘭吐雅
    光譜學(xué)與光譜分析 2016年10期
    關(guān)鍵詞:內(nèi)蒙古草原利用率生物量

    包 剛,辛?xí)云?,包玉海,王牧蘭, 3,元志輝, 3,烏蘭吐雅

    1. 呼倫貝爾草原生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站/中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081

    2. 內(nèi)蒙古師范大學(xué)內(nèi)蒙古自治區(qū)遙感與地理信息系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,內(nèi)蒙古 呼和浩特 010022

    3. 內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010022

    4. 內(nèi)蒙古自治區(qū)農(nóng)牧業(yè)科學(xué)院,內(nèi)蒙古 呼和浩特 010031

    內(nèi)蒙古草原植被最大光能利用率取值優(yōu)化研究

    包 剛1, 2,辛?xí)云?*,包玉海2,王牧蘭2, 3,元志輝2, 3,烏蘭吐雅4

    1. 呼倫貝爾草原生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站/中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,北京 100081

    2. 內(nèi)蒙古師范大學(xué)內(nèi)蒙古自治區(qū)遙感與地理信息系統(tǒng)重點(diǎn)實(shí)驗(yàn)室,內(nèi)蒙古 呼和浩特 010022

    3. 內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院,內(nèi)蒙古 呼和浩特 010022

    4. 內(nèi)蒙古自治區(qū)農(nóng)牧業(yè)科學(xué)院,內(nèi)蒙古 呼和浩特 010031

    針對(duì)目前CASA (carnegie-ames-stanford approach)模型等植被生產(chǎn)力模型植被最大光能利用率的取值未對(duì)草原進(jìn)行區(qū)分的問(wèn)題,以內(nèi)蒙古草甸草原、典型草原和荒漠草原為研究對(duì)象,結(jié)合野外實(shí)測(cè)NPP(net primary productivity)數(shù)據(jù)和CASA模型的建模思路優(yōu)化了三大草原類(lèi)型植被最大光能利用率,并以此為基礎(chǔ)模擬分析了其植被光能利用率和NPP時(shí)空格局。結(jié)果表明,基于99個(gè)地面采樣點(diǎn)所建立的一元二次方程模擬的草甸草原、典型草原和荒漠草原最大光能利用率分別為0.654,0.553和0.511 gC·MJ-1,平均為0.573 gC·MJ-1。與未對(duì)草原類(lèi)型進(jìn)行區(qū)分而統(tǒng)一取值為0.541 gC·MJ-1的結(jié)果相比,實(shí)測(cè)NPP與模擬NPP之間的決定系數(shù)和和均方根誤差分別提高了0.024和2.62 gC·(m2·month-1)-1。受水熱組合和草原類(lèi)型的空間格局的影響,內(nèi)蒙古草原植被光能利用率和NPP總體上由東北向西南逐漸下降趨勢(shì),呈明顯的單峰季節(jié)變化特征。但光能利用率和NPP的最大值出現(xiàn)的月份有所不同,分別出現(xiàn)在8月份和7月份,這可能與植被吸收的光合有效輻射和光能利用率的最高值出現(xiàn)的月份不同有關(guān)。光能利用率和NPP平均值按草甸草原>典型草原>荒漠草原的順序依次降低。

    最大光能利用率;草原類(lèi)型;CASA模型;內(nèi)蒙古;時(shí)空格局

    引 言

    光能利用率(light use efficiency, LUE,單位:gC·M·jol-1)是指植被把所吸收的光合有效輻射轉(zhuǎn)化為有機(jī)碳的效率,是估算植被凈初級(jí)生產(chǎn)力(net primary productivity,NPP)的一個(gè)關(guān)鍵變量。在估算植被生產(chǎn)力的諸多光能利用率模型中,通過(guò)對(duì)最大光能利用率(εmax)產(chǎn)生影響的環(huán)境綜合脅迫因子來(lái)計(jì)算植被光能利用率[1-3]。植被最大光能利用率是指植被在沒(méi)有任何限制條件的理想條件下通過(guò)光合作用將所截獲并吸收的光合有效輻射轉(zhuǎn)化為有機(jī)干物質(zhì)的效率,其取值因不同的植被類(lèi)型而有所不同[4]。其取值大小直接影響植被實(shí)際光能利率的計(jì)算和最終的凈初級(jí)生產(chǎn)力的模擬結(jié)果,是影響初級(jí)生產(chǎn)力估算模型精度或各模型模擬結(jié)果之間存在一定偏差的主要原因之一[4-5]。綜上所述,不同類(lèi)型植被最大光能利用率取值優(yōu)化對(duì)改進(jìn)生產(chǎn)力模型,提高模擬精度和減小不確定性等具有重要意義。

    作為光能利用率模型的代表,CASA(carnegie-ames-stanford approach)模型已經(jīng)在不同尺度植被凈初級(jí)生產(chǎn)力模擬及陸地生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)研究中得到較廣泛應(yīng)用[1, 3, 6-7]。但由于所采用的實(shí)測(cè)數(shù)據(jù)不同以及最大光能利用率本身因不同的植被類(lèi)型或同一種植被類(lèi)型在亞類(lèi)的不同或生態(tài)環(huán)境條件的不同而有所差異[4-5],在CASA模型的最大光能利用率取值上一直存在偏差。這將直接導(dǎo)致最終凈初級(jí)生產(chǎn)力模擬結(jié)果的不確定。Potter和Field在全球尺度沒(méi)有區(qū)分植被類(lèi)型,統(tǒng)一取值為0.389 gC·MJ-1[1-2]。彭少麟將廣東植被最大光能利用率取值為1.25 gC·MJ-1[8]。Mao等將東北地區(qū)沼澤濕地植被最大光能利用率取值為0.542 gC·MJ-1[3]。朱文泉利用全國(guó)690個(gè)觀測(cè)點(diǎn)的實(shí)測(cè)NPP數(shù)據(jù),模擬獲得中國(guó)10種典型植被最大光能利用率介于0.389~0.985之間[4]。Yu等在東亞地區(qū)模擬得到的13種主要植被類(lèi)型的最大光能利用率介于0.389~0.978之間[5]。盡管上述取值在其對(duì)應(yīng)尺度上具有應(yīng)用價(jià)值,但直接應(yīng)用于具體植被類(lèi)型或研究區(qū)可能帶來(lái)較大誤差。需要指出,CASA模型所取最大光能利用率值未對(duì)草甸草原、典型草原、荒漠草原等草原類(lèi)型進(jìn)行進(jìn)一步的細(xì)分,而所有草原類(lèi)型統(tǒng)一取值為0.389[1]或0.541[5],這可能對(duì)草原植被生產(chǎn)力模擬帶來(lái)不確定性。因此針對(duì)不同草原類(lèi)型,結(jié)合相應(yīng)實(shí)測(cè)NPP數(shù)據(jù),確定不同類(lèi)型草原植被最大光能利用率,基于此更準(zhǔn)確反映實(shí)際光能利用率空間格局和更高精度地估算草原植被NPP更為重要。

    鑒于此,本研究以內(nèi)蒙古草原植被為研究對(duì)象,結(jié)合內(nèi)蒙古草原植被空間數(shù)據(jù)、地面實(shí)測(cè)NPP數(shù)據(jù)、氣象數(shù)據(jù)和MODIS-NDVI等數(shù)據(jù),根據(jù)CASA模型的建模思路開(kāi)展草甸草原、典型草原和荒漠草原的最大光能利用率取值優(yōu)化研究,探討優(yōu)化不同草原類(lèi)型植被最大光能利用率對(duì)NPP模擬精度的影響。

    1 實(shí)驗(yàn)部分

    CASA模型的驅(qū)動(dòng)變量主要包括每月的NDVI數(shù)據(jù)、柵格化的氣象數(shù)據(jù)(月平均溫度、月總降水量、月太陽(yáng)總輻射)和植被類(lèi)型數(shù)據(jù),其來(lái)源及處理過(guò)程分別如下所述。

    1.1 數(shù)據(jù)處理

    1.1.1 NDVI數(shù)據(jù)

    本文選取2009年4—10月(地面樣方數(shù)據(jù)為2009年)空間分辨率為500 m×500 m,時(shí)間分辨率為16 d的MODIS NDVI產(chǎn)品MOD13A1作為CASA模型輸入。每月NDVI數(shù)據(jù)通過(guò)最大值合成法(maximum value composite, MVC)獲得[9]。利用MODIS Reprojection Tools(MRT)工具對(duì)NDVI數(shù)據(jù)進(jìn)行拼接、投影和數(shù)據(jù)格式轉(zhuǎn)換等預(yù)處理,采用Albers投影。

    1.1.2 草地生物量數(shù)據(jù)及NPP轉(zhuǎn)換

    草地生物量數(shù)據(jù):草地生物量數(shù)據(jù)為2009年和2012年野外調(diào)查的樣方數(shù)據(jù)。結(jié)合內(nèi)蒙古植被類(lèi)型圖,在草甸草原、典型草原和荒漠草原范圍內(nèi)共設(shè)計(jì)110個(gè)大小為1 m×1 m樣方(圖1),采用收割法和挖掘法獲取草地地上和地下生物量。其中,2009年樣方為99個(gè),不包含地下生物量;2012年樣方為11個(gè),包括地上和地下生物量,在此主要用于確定典型草原地下和地上生物量比例[10]。在每個(gè)樣方內(nèi)進(jìn)行植物高度、種類(lèi)和覆蓋度等調(diào)查之后,收割其所有地上綠色植物裝入土樣品袋,編號(hào),帶回實(shí)驗(yàn)室稱(chēng)其鮮重,然后在70 ℃恒溫條件下烘干至恒重,稱(chēng)量地上干重。地下生物量挖掘法:由于地下生物量的采集相對(duì)復(fù)雜和工作量大等現(xiàn)實(shí),地下生物量的采集只在2012年11個(gè)樣方內(nèi)開(kāi)展。完成地上生物量的調(diào)查之后,采用挖掘法在樣方正中間挖取大小為0.25 m×0.25 m的根樣(即樣方的1/4),將挖取的根樣連土一并帶回實(shí)驗(yàn)室晾干后,置于孔徑為0.150 mm紗網(wǎng)中沖洗,晾干和編碼等,將根系和土壤完全分離。取樣深度視植物根系分布深淺而定,一般為40 cm~70 cm之間。將由土壤分離后的根系放入烘箱在70 ℃恒溫條件下烘干至恒重,稱(chēng)量干重并乘以16獲得1 m×1 m樣方地下生物量。由于2012年11個(gè)樣方全部分布在典型草原區(qū),根據(jù)地下生物量干重和地上生物量干重,得到內(nèi)蒙古典型草原地下/地上生物量比值平均值4.127。而草甸草原和荒漠草原的比值來(lái)源于Piao等研究結(jié)果,分別取5.26和7.89[10]。

    圖1 研究區(qū)草原類(lèi)型空間格局及采樣點(diǎn)分布

    生物量數(shù)據(jù)的NPP轉(zhuǎn)換:生物量是生態(tài)系統(tǒng)一定時(shí)間內(nèi)單位面積所包含的活的生物有機(jī)體的總干物質(zhì)量,單位為g·m-2,是凈生產(chǎn)力的存留部分。而凈初級(jí)生產(chǎn)力是指在單位時(shí)間內(nèi)在單位面積上所累積的有機(jī)物質(zhì)量,一般以碳為單位[即碳通量單位:gC·(m2·t-1)-1],是生物量形成的基礎(chǔ),在一定程度上代表生物量形成的速度。因此,基于地面實(shí)測(cè)NPP數(shù)據(jù)優(yōu)化草原植被最大光能利用和驗(yàn)證模型模擬結(jié)果時(shí)必須先將地面實(shí)測(cè)生物量轉(zhuǎn)換為NPP,以保證模型輸出NPP和地面實(shí)測(cè)NPP的一致。在此,利用Gill等提出的方法將實(shí)測(cè)生物量轉(zhuǎn)換為NPP[11]

    NPP=ANPP+BNPP

    (1)

    BNPP=BGB×(liveBGB/BGB)×turnover

    (2)

    turnover=0.000 9(g·m-2)×ANPP+0.25

    (3)

    式中,ANPP和BNPP分別表示地上和地下NPP;BGB為累計(jì)的地下生物量,通過(guò)野外采樣或不同類(lèi)型草原地上地下生物量比值來(lái)計(jì)算獲取[10]。liveBGB為當(dāng)年地下活生物量,liveBGB/BGB表示活根系生物量占總根系生物量的比例,取Gill等確定的值liveBGB/BGB=0.6;turnover為草地根系周轉(zhuǎn)值。將地上生物量按照碳素比例(36.98%)換算成含碳量(ANPP)[12]。

    1.1.3 氣象數(shù)據(jù)

    氣候數(shù)據(jù)采用文獻(xiàn)[13]蒙古高原柵格化氣象數(shù)據(jù)產(chǎn)品。此數(shù)據(jù)產(chǎn)品來(lái)自于蒙古高原108個(gè)標(biāo)準(zhǔn)氣象觀測(cè)站(內(nèi)蒙古48個(gè)站點(diǎn),蒙古國(guó)60個(gè)站點(diǎn))的月平均氣溫、月總降水量和22個(gè)太陽(yáng)輻射站的太陽(yáng)輻射數(shù)據(jù)(內(nèi)蒙古8個(gè)站點(diǎn),蒙古國(guó)14個(gè)站點(diǎn))。利用內(nèi)蒙古行政界線對(duì)蒙古高原柵格化氣象數(shù)據(jù)進(jìn)行掩膜處理獲得內(nèi)蒙古范圍的氣象數(shù)據(jù)。

    1.1.4 植被類(lèi)型數(shù)據(jù)

    植被類(lèi)型數(shù)據(jù)來(lái)自數(shù)字化的1∶100萬(wàn)內(nèi)蒙古植被圖。在ArcGIS10.0環(huán)境下對(duì)植被圖進(jìn)行投影轉(zhuǎn)換、數(shù)字化和柵格化等處理,獲得投影和分辨率與遙感數(shù)據(jù)一致的植被數(shù)字化產(chǎn)品,用于不同類(lèi)型草原空間格局的提取。

    1.2 研究方法

    根據(jù)CASA模型的建模思路[式(4)和式(5)][1-2],采用最小二乘法對(duì)內(nèi)蒙古三大草原植被最大光能利用率進(jìn)行優(yōu)化。從CASA模型的表達(dá)式[式(4)和式(5)]可以看出,如果先確定某種草原類(lèi)型的一組NPP,APAR,溫度脅迫因子和水分脅迫因子,可根據(jù)誤差最小原則建立開(kāi)口向上的一元二次方程[式(6)和式(7)],模擬其最小值所對(duì)應(yīng)的橫坐標(biāo)值,即最大光能利用率[4-5]。

    (4)

    (5)

    式中,APAR(x,t)表示像元x在t月吸收的光合有效輻射,ε(x,t)表示像元x在t月的實(shí)際光能利用;εmax表示在理想條件植被最大光能利用;Tε1(x,t)和Tε2(x,t)表示低溫和高溫對(duì)光能利用率的脅迫作用,而Wε(x,t)表示水分對(duì)光能利用率的脅迫作用。

    (6)

    (7)

    式中,i表示某一草原類(lèi)型的樣本數(shù),j為某一草原類(lèi)型的最大樣本數(shù),m為NPP實(shí)測(cè)數(shù)據(jù),來(lái)自野外樣本數(shù)據(jù)[式(1)—式(3)];n來(lái)自CASA模型的基本表達(dá)式,即n=APAR×Tε1×Tε2×Wε,ε*為待定的某一草原類(lèi)型的最大光能利用率。

    在優(yōu)化不同草原類(lèi)型最大光能利用率的基礎(chǔ)上,通過(guò)式(4)和式(5)分別模擬內(nèi)蒙古草原植被實(shí)際光能利用和植被凈初級(jí)生產(chǎn)力。為評(píng)價(jià)區(qū)分草原類(lèi)型優(yōu)化最大光能利用率對(duì)NPP模擬結(jié)果的影響,分別將Yu等[5]在東亞所有草原類(lèi)型統(tǒng)一使用的最大光能利用率值0.541和本文優(yōu)化結(jié)果作為變量輸入模型,分別模擬出NPP,并與實(shí)測(cè)NPP進(jìn)行比較。

    2 結(jié)果與討論

    2.1 三大草原植被最大光能利用率

    基于式(4)—式(7),結(jié)合野外考察獲取的草甸草原、典型草原、荒漠草原實(shí)測(cè)NPP數(shù)據(jù),建立了優(yōu)化三大草原植被最大光能利用的一元二次方程[式(8)],并通過(guò)對(duì)其進(jìn)行一階求導(dǎo)得到最大光能利用率值,分別為0.654,0.553和0.511 gC·MJ-1。三類(lèi)草原最大光能利用率平均值為0.573 gC·MJ-1。表1為用于建立一元二次方程的99個(gè)觀測(cè)點(diǎn)實(shí)測(cè)NPP和APAR×Tε1×Tε2×Wε值。

    (8)

    表1 計(jì)算不同類(lèi)型草原植被最大光能利用率的模型輸入數(shù)據(jù)

    續(xù)表1

    HT472 378 97HT37156 27118 85ED464 3939 78HT593 3136 67HT38252 72166 90ED568 2651 37HT683 0948 34HT39109 7654 77ED646 2131 73HT784 2110 99HT40169 13122 37ED791 4333 32HT8105 9850 40HT41125 4295 68ED8129 2576 57HT948 5420 12HT42154 74125 91ED9111 3569 64HT1079 4067 82HT43226 20117 50ED10182 9190 12HT1161 6522 14HT44219 73151 08ED11115 1784 96HT12148 6045 14XT1207 8075 61ED12174 1079 82HT13194 5291 22XT2158 2081 38ED13107 1359 57HT14205 5674 43XT392 0764 50ED14144 7378 19HT1562 9612 98XT4152 75119 77ED15175 0461 36

    注:n和m對(duì)應(yīng)于式(7),分別表示APAR×Tε1×Tε2×Wε和實(shí)測(cè)NPP,H,X,E代表呼倫貝爾、錫林郭勒草原和鄂爾多斯高原,M,T,D分別代表草甸草原、典型草原和荒漠草原

    Note:nandmcorrespond to equation (7), represent APAR×Tε1×Tε2×Wεand field surveyed NPP, respectively, H, X, E represent Hulun Buir, Xilingol grassland and ordos pluteau, M, T, D respresent meadow, typical and desert steppe, respectively

    2.2 NPP模擬精度驗(yàn)證

    為評(píng)價(jià)分草原類(lèi)型優(yōu)化最大光能利用率后對(duì)草地NPP模擬的影響,將優(yōu)化后的最大光能利用率值和其他相關(guān)變量輸入CASA模型模擬獲得了樣方觀測(cè)點(diǎn)的NPP。圖2(a)顯示了分草原類(lèi)型(草甸草原、典型草原、荒漠草原)確定最大光能利用率后模擬得到的NPP與實(shí)測(cè)值之間的關(guān)系??傮w來(lái)看,NPP實(shí)測(cè)值與模擬值之間有良好的線性關(guān)系,兩者之間的決定系數(shù)達(dá)到0.78,通過(guò)0.001的顯著性水平[圖2(a)]。綜合考慮模擬誤差和模型模擬常用精度范圍得出[6],本文模擬結(jié)果能夠反映研究區(qū)草原植被NPP的變化狀況,可以用來(lái)分析NPP時(shí)空變化格局。與Yu等[5]未對(duì)草原進(jìn)行區(qū)分而統(tǒng)一取0.541 gC·MJ-1的結(jié)果相比[圖2(b)],本研究的模擬精度有所提升,模擬NPP和實(shí)測(cè)值之間的決定系數(shù)和均方根誤差分別提高了0.024和2.62 gC·(m2·month-1)-1,且優(yōu)化不同草原類(lèi)型最大光能利用率后,擬合線更接近于1∶1線。

    圖2 NPP模擬值與實(shí)測(cè)值的比較

    圖3 2009年內(nèi)蒙古草原植被生長(zhǎng)季各月份光能利用率空間格局

    2.3 內(nèi)蒙古草原植被光能利用率空間分布

    圖3為根據(jù)CASA模型的光能利用率子模型[式(5)]計(jì)算得到的內(nèi)蒙古草原生長(zhǎng)季光能利用率空間分布格局。由圖3可以看出內(nèi)蒙古草原植被生長(zhǎng)季光能利用率總體分布與其草原類(lèi)型的分布基本一致,高值區(qū)主要分布在降水相對(duì)豐富的東北部草甸草原區(qū),而低值區(qū)主要分布在西南部荒漠草原區(qū)。中部典型草原的光能利用率處于中間水平,各月份平均值約為0.45 gC·MJ-1。從月際變化看,隨著時(shí)間的推移,內(nèi)蒙古草原植被光能利用率呈先增加后下降的單峰趨勢(shì),與其溫度和降水量的變化高度一致(圖3),統(tǒng)計(jì)結(jié)果進(jìn)一步證實(shí)了這一點(diǎn)(圖4)。總體看,所有類(lèi)型草原光能利用率在8月份達(dá)到最大值,分別為0.605,0.431和0.355 gC·MJ-1,其次為7月份為0.598,0.430和0.310 gC·MJ-1,最低值主要出現(xiàn)在4月份,分別為0.228,0.242和0.219 gC·MJ-1。三大草原平均光能利用率曲線季節(jié)變化規(guī)律與典型草原比較一致,說(shuō)明典型草原光能利用率最為接近內(nèi)蒙古所有類(lèi)型草原植被光能利用率的平均水平(圖4)。

    圖4 內(nèi)蒙古草地生長(zhǎng)季光能利用率變化

    2.4 內(nèi)蒙古草原植被NPP空間格局

    由圖5可以看出,與植被光能利用率空間格局相似,內(nèi)蒙古草原植被NPP分布同樣表現(xiàn)為由東北向西南呈逐漸遞減趨勢(shì),該分布特征與研究區(qū)水熱梯度和植被類(lèi)型變化保持了較好的一致性,但各月份NPP大小有明顯區(qū)別,如6,7和8月份研究東北部草甸草原區(qū)NPP達(dá)到了130 gC·(m2·month-1)-1以上,而中部30~90 gC·(m2·month-1)-1之間;4和10月份所有草原類(lèi)型NPP均低于30 gC·(m2·month-1)-1(圖5)。與光能利用率的月變化不同,草甸草原、典型草原和荒漠化草原植被NPP的最高值均出現(xiàn)在7月,分別為164.2,69.1和22.9 gC·(m2·month-1)-1,其次為8月份,為156.6,62.2和19.7 gC·(m2·月-1)-1,這可能草原植被在7月份吸收的光合有效輻射比8月份大有關(guān)。最小值主要集中在10月份(圖6)。

    圖5 2009年內(nèi)蒙古草原植被生長(zhǎng)季各月份NPP空間格局

    圖6 內(nèi)蒙古草地生長(zhǎng)季NPP變化

    3 結(jié) 論

    以內(nèi)蒙古草原類(lèi)型數(shù)據(jù)、氣象數(shù)據(jù)和MODIS NDVI數(shù)據(jù)為數(shù)據(jù)源,結(jié)合地面實(shí)測(cè)NPP數(shù)據(jù)和CASA模型的建模思路,建立計(jì)算草甸草原、典型草原和荒漠草原最大光能利用率的一元二次多項(xiàng)式,開(kāi)展了三大草原植被最大光能利用率取值優(yōu)化研究。通過(guò)研究得出以下結(jié)論:

    (1)經(jīng)過(guò)一元二次多項(xiàng)式模擬優(yōu)化的內(nèi)蒙古草甸草原、典型草原和荒漠草原的植被最大光能利用率分別為0.654,0.553和0.511 gC·MJ-1,平均為0.573 gC·MJ-1。

    (2)與最大光能利用率的取值未對(duì)草原(草甸草原、典型草原、荒漠草原)類(lèi)型進(jìn)行區(qū)分的研究結(jié)果相比,優(yōu)化不同類(lèi)型草原最大光能利用率能夠提高CASA模型的NPP模擬精度。與統(tǒng)一取值為0.541 gC·MJ-1的結(jié)果相比,實(shí)測(cè)NPP和模擬NPP之間的決定系數(shù)和均方根誤差分別提升了0.024和2.62 gC·(m2·month-1)-1。這說(shuō)明對(duì)植被類(lèi)型進(jìn)行更進(jìn)一步的優(yōu)化,如優(yōu)勢(shì)種群等的優(yōu)化并確定其最大光能利用可能進(jìn)一步提高NPP模擬精度。

    (3)內(nèi)蒙古草原植被光能利用率和NPP總體上呈由高緯度的東北部地區(qū)向低緯度的西南部地區(qū)逐漸下降的趨勢(shì),與其水熱組合和草原類(lèi)型空間格局基本一致。生長(zhǎng)季所有草原類(lèi)型平均光能利用率和NPP季節(jié)變化呈明顯的單峰變化,但最大值出現(xiàn)的月份有所不同。光能利用率的最大值出現(xiàn)在8月份,而NPP最大值出現(xiàn)在7月份,可能與7月份植被吸收的光合有效輻射比8月份相對(duì)高有關(guān)。光能利用率和NPP平均值按草甸草原>典型草原>荒漠草原的順序依次降低。

    [1] Potter C S, Randerson J T, Field C B, et al. Global Biogeochemical Cycles, 1993, 7(4): 811.

    [2] Field C B, Randerson J T, Malmstr?m C M. Remote Sensing of Environment, 1995, 51(1): 74.

    [3] Mao D, Wang Z, Li L, et al. Ecological Engineering, 2014, 68: 97.

    [4] Zhu W, Pan Y, He H, et al. Chinese Science Bulletin, 2006, 51(4): 457.

    [5] Yu D, Shi P, Shao H, et al. International Journal of Remote Sensing, 2009, 30(18): 4851.

    [6] Mu S, Zhou S, Chen Y, et al. Global and Planetary Change, 2013, 108: 29.

    [7] MAO De-hua, WANG Zong-ming, LUO Ling, et al(毛德華, 王宗明, 羅 玲). Chinese Journal of Applied Ecology(應(yīng)用生態(tài)學(xué)報(bào)), 2012, 23(6): 1511.

    [8] PENG Shao-lin, GUO Zhi-hua, WANG Bo-sun(彭少麟, 郭志華, 王伯蓀). Acta Ecologica Sinica (生態(tài)學(xué)報(bào)), 2000, 20(6): 903.

    [9] Holben B N. International Journal of Remote Sensing, 1986, 7 (11): 1417.

    [10] Piao S, Fang J, Zhou L, et al. Global Biogeochemical Cycles, 2007, 21(2): GB2002.

    [11] Gill R, Kelly R, Parton W, et al. Global Ecology and Biogeography, 2002, 11(1): 79.

    [12] Chen Z, Shao Q, Liu J, et al. Science China Earth Sciences, 2012, 55(8): 1306.

    [13] Bao G, Qin Z, Bao Y, et al. Remote Sensing, 2014, 6: 8337.

    (Received Sep. 7, 2015; accepted Jan. 11, 2016)

    *Corresponding author

    Optimization of Maximum Light Use Efficiency in Inner Mongolian Steppe

    BAO Gang1, 2, XIN Xiao-ping1*, BAO Yu-hai2, WANG Mu-lan2, 3, YUAN Zhi-hui2, 3, Wulantuya4

    1. Hulunber Grassland Ecosystem Observation and Research Station/Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences,Beijing 100081,China

    2. Inner Mongolian Key Laboratory of Remote Sensing and Geographic Information System of Inner Mongolia Normal University,Huhhot 010022,China

    3. College of Geographical Science, Inner Mongolia Normal University,Huhhot 010022, China

    4. Inner Mongolia Academy of Agricultural &Animal Husbandry Sciences, Huhhot 010031, China

    For the case that the value of the maximum light use efficiency (MLUE) is not optimized for different steppes, we simulated the MLUE for meadow steppe, typical steppe and desert steppe in Inner Mongolia based on the field observed NPP and CASA ecosystem model, and analyzed the spatial and temporal pattern of the LUE and net primary productivity (NPP) in Inner Mongolia. The result indicate that the MLUE is optimized to be 0.654,0.553 and 0.511 gC·MJ-1for meadow steppe, typical steppe and desert steppe in Inner Mongolia, respectively, with an average of 0.573 gC·MJ-1. Compared to the result that used same value of 0.541 gC·MJ-1for MLUE, the correlation coefficient and relative mean square error was improved 0.024 and 2.62 gC·(m2·month-1)-1, respectively after optimization. Affected by the hydrothermal condition and distribution of grassland types, the LUE and NPP in Inner Mongolia decreased from northeast to southwest, and showed one crest shape. However, the maximum value of LUE and NPP was appeared in August and July. This difference could be attributed to the difference in the maximum value between absorbed photosynthetically active radiation and LUE. The LUE and NPP decreased by meadow steppe, typical steppe and desert steppe.

    Maximum light use efficiency;Steppe;CASA model;Inner Mongolia;Spatiotemporal pattern

    2015-09-07,

    2016-01-11

    國(guó)家自然科學(xué)基金項(xiàng)目(41301456),國(guó)家科技支撐計(jì)劃項(xiàng)目(2013BAC03B02),現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專(zhuān)項(xiàng)資金(CARS-35)和內(nèi)蒙古自治區(qū)高等學(xué)??茖W(xué)研究項(xiàng)目(NJZY029)資助

    包 剛,1978年生,呼倫貝爾草原生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站博士研究生 e-mail: baogang@imnv.edu.cn *通訊聯(lián)系人 e-mail: xinxiaoping@caas.cn

    S123

    A

    10.3964/j.issn.1000-0593(2016)10-3280-07

    猜你喜歡
    內(nèi)蒙古草原利用率生物量
    輪牧能有效促進(jìn)高寒草地生物量和穩(wěn)定性
    化肥利用率穩(wěn)步增長(zhǎng)
    做好農(nóng)村土地流轉(zhuǎn) 提高土地利用率
    第十五屆中國(guó)·內(nèi)蒙古草原文化節(jié)“守望相助在美好家園”蒙古語(yǔ)歌曲創(chuàng)作筆會(huì)在呼舉辦
    草原歌聲(2018年2期)2018-12-03 08:14:06
    第十五屆中國(guó)·內(nèi)蒙古草原文化節(jié)“永遠(yuǎn)的春天”——內(nèi)蒙古自治區(qū)紀(jì)念改革開(kāi)放四十周年合唱展演圓滿落幕
    草原歌聲(2018年2期)2018-12-03 08:13:56
    淺議如何提高涉煙信息的利用率
    內(nèi)蒙古草原上織起幸福的網(wǎng)
    內(nèi)蒙古草原旅游區(qū)觀光型家庭牧場(chǎng)的發(fā)展研究
    板材利用率提高之研究
    生物量高的富鋅酵母的開(kāi)發(fā)應(yīng)用
    人妻丰满熟妇av一区二区三区| 国产精品久久视频播放| 久久久久久大精品| 91在线精品国自产拍蜜月 | 久久中文看片网| 天天添夜夜摸| 午夜精品在线福利| 欧美色欧美亚洲另类二区| 制服人妻中文乱码| 精品国产亚洲在线| 草草在线视频免费看| 国产欧美日韩精品亚洲av| 久久精品综合一区二区三区| 欧美日本亚洲视频在线播放| 精品福利观看| 又黄又爽又免费观看的视频| 亚洲无线观看免费| 香蕉久久夜色| 亚洲一区二区三区不卡视频| 国产亚洲精品av在线| 99久久99久久久精品蜜桃| 丁香六月欧美| 亚洲在线观看片| 99久久成人亚洲精品观看| 国内精品久久久久精免费| 国产精品野战在线观看| 日韩欧美在线乱码| 亚洲av中文字字幕乱码综合| netflix在线观看网站| 亚洲自偷自拍图片 自拍| 日韩免费av在线播放| 99国产极品粉嫩在线观看| 深夜精品福利| 精品久久久久久久末码| 久久人妻av系列| 成人永久免费在线观看视频| av国产免费在线观看| 国产精品一区二区三区四区久久| 91久久精品国产一区二区成人 | 高清在线国产一区| 麻豆成人av在线观看| 啦啦啦免费观看视频1| 欧美黄色片欧美黄色片| 国产美女午夜福利| 亚洲中文日韩欧美视频| www.自偷自拍.com| 一个人免费在线观看的高清视频| 日本撒尿小便嘘嘘汇集6| 久久久久久久久中文| 久久香蕉精品热| av欧美777| 亚洲精品中文字幕一二三四区| 久久久久国产一级毛片高清牌| 91在线观看av| 国产乱人视频| 久久久国产精品麻豆| 国产美女午夜福利| 欧美高清成人免费视频www| 97碰自拍视频| 日韩欧美免费精品| 99热这里只有精品一区 | 久久香蕉国产精品| 亚洲色图 男人天堂 中文字幕| 毛片女人毛片| 久久精品国产清高在天天线| 午夜免费成人在线视频| 黄色女人牲交| 国产亚洲欧美在线一区二区| 两个人视频免费观看高清| 久久午夜亚洲精品久久| 国产精品久久久久久人妻精品电影| 18美女黄网站色大片免费观看| 免费在线观看成人毛片| 青草久久国产| cao死你这个sao货| 狠狠狠狠99中文字幕| 51午夜福利影视在线观看| 国产精品九九99| 在线观看美女被高潮喷水网站 | 日韩欧美在线二视频| 淫妇啪啪啪对白视频| 国产伦精品一区二区三区四那| 91老司机精品| 窝窝影院91人妻| 99国产精品99久久久久| 色综合婷婷激情| 一本综合久久免费| 婷婷精品国产亚洲av| 97碰自拍视频| 国产精品影院久久| 免费av毛片视频| 免费大片18禁| 国产成人福利小说| 我的老师免费观看完整版| 色噜噜av男人的天堂激情| 国产精品av久久久久免费| 三级毛片av免费| 一级毛片女人18水好多| 亚洲中文字幕一区二区三区有码在线看 | 中文字幕高清在线视频| 夜夜看夜夜爽夜夜摸| 少妇丰满av| 一夜夜www| 国产97色在线日韩免费| 国产成人精品久久二区二区免费| 动漫黄色视频在线观看| 国产99白浆流出| 欧美黄色片欧美黄色片| 国产亚洲av嫩草精品影院| 法律面前人人平等表现在哪些方面| 亚洲午夜理论影院| 老司机在亚洲福利影院| 欧美日本视频| 精品国产亚洲在线| 国产又色又爽无遮挡免费看| 国产精品女同一区二区软件 | 男女午夜视频在线观看| 手机成人av网站| 国产成人啪精品午夜网站| 99热这里只有是精品50| 亚洲国产欧洲综合997久久,| 舔av片在线| 久9热在线精品视频| 成人性生交大片免费视频hd| 特大巨黑吊av在线直播| 国产精品野战在线观看| 亚洲电影在线观看av| 精品久久久久久,| 亚洲国产欧美人成| 国产精品久久久久久人妻精品电影| 精品久久蜜臀av无| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美国产在线观看| 激情在线观看视频在线高清| 亚洲专区中文字幕在线| 91麻豆精品激情在线观看国产| 中亚洲国语对白在线视频| 美女高潮的动态| 国内久久婷婷六月综合欲色啪| 国产综合懂色| 国产免费男女视频| 日本成人三级电影网站| 黄色日韩在线| 欧美日韩瑟瑟在线播放| 亚洲熟女毛片儿| 亚洲av中文字字幕乱码综合| 激情在线观看视频在线高清| 亚洲国产高清在线一区二区三| 日韩欧美国产一区二区入口| 国产野战对白在线观看| 欧美黄色片欧美黄色片| 黄色视频,在线免费观看| 国产三级黄色录像| 亚洲精品中文字幕一二三四区| 好男人在线观看高清免费视频| 精品一区二区三区四区五区乱码| 成人无遮挡网站| 国产淫片久久久久久久久 | 一二三四在线观看免费中文在| 久久久久亚洲av毛片大全| 天堂动漫精品| 国产视频内射| 国产精品一区二区精品视频观看| 国产高清视频在线观看网站| 国产日本99.免费观看| netflix在线观看网站| 国产精品香港三级国产av潘金莲| 香蕉丝袜av| 亚洲成a人片在线一区二区| 亚洲色图av天堂| 热99在线观看视频| 欧美黑人巨大hd| 一区二区三区国产精品乱码| 真人做人爱边吃奶动态| 少妇丰满av| 久久久色成人| 国产伦一二天堂av在线观看| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 91麻豆av在线| 亚洲一区二区三区色噜噜| 亚洲av电影在线进入| 性欧美人与动物交配| 三级男女做爰猛烈吃奶摸视频| 一二三四社区在线视频社区8| 18禁黄网站禁片午夜丰满| www日本在线高清视频| 国产精品久久久人人做人人爽| 好男人在线观看高清免费视频| 国产亚洲精品一区二区www| 少妇人妻一区二区三区视频| 又大又爽又粗| 亚洲一区二区三区不卡视频| 美女高潮喷水抽搐中文字幕| 在线免费观看不下载黄p国产 | 又大又爽又粗| 一个人看视频在线观看www免费 | 在线观看日韩欧美| 成人高潮视频无遮挡免费网站| 久久久成人免费电影| 俺也久久电影网| 亚洲午夜理论影院| 一二三四社区在线视频社区8| 成年女人毛片免费观看观看9| 亚洲国产高清在线一区二区三| 听说在线观看完整版免费高清| 亚洲avbb在线观看| 又黄又爽又免费观看的视频| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 99视频精品全部免费 在线 | 成年女人看的毛片在线观看| 搡老妇女老女人老熟妇| 日韩大尺度精品在线看网址| 国产欧美日韩一区二区三| 1024手机看黄色片| 国产精品野战在线观看| 18禁国产床啪视频网站| 夜夜夜夜夜久久久久| 亚洲国产高清在线一区二区三| 亚洲狠狠婷婷综合久久图片| 丰满人妻熟妇乱又伦精品不卡| 久久精品人妻少妇| 久久天躁狠狠躁夜夜2o2o| 国产伦精品一区二区三区四那| 午夜a级毛片| 18禁黄网站禁片午夜丰满| 他把我摸到了高潮在线观看| 国产单亲对白刺激| 神马国产精品三级电影在线观看| 嫩草影院精品99| 色综合欧美亚洲国产小说| 黑人操中国人逼视频| 国产亚洲av高清不卡| 午夜日韩欧美国产| 久久国产精品影院| 日韩欧美一区二区三区在线观看| 色综合亚洲欧美另类图片| 爱豆传媒免费全集在线观看| 少妇高潮的动态图| av免费在线看不卡| 欧美性感艳星| 一区二区三区免费毛片| 97在线视频观看| 99久久无色码亚洲精品果冻| 亚洲不卡免费看| 在线观看66精品国产| 黑人高潮一二区| 久久国内精品自在自线图片| 啦啦啦观看免费观看视频高清| 欧美日韩综合久久久久久| 美女黄网站色视频| 免费观看的影片在线观看| 青春草视频在线免费观看| 久久热精品热| 变态另类丝袜制服| 非洲黑人性xxxx精品又粗又长| 尾随美女入室| 91aial.com中文字幕在线观看| 亚洲国产欧美人成| 桃色一区二区三区在线观看| 亚洲无线观看免费| 国内揄拍国产精品人妻在线| 国内精品宾馆在线| 国产色婷婷99| 精品一区二区三区人妻视频| 国产精品三级大全| 秋霞在线观看毛片| 国产探花在线观看一区二区| 国产三级中文精品| 视频中文字幕在线观看| 亚洲国产成人一精品久久久| 亚洲在久久综合| 国产三级在线视频| 嫩草影院精品99| 国产精品蜜桃在线观看| 久久久久久久久久久免费av| 五月伊人婷婷丁香| 久久久久久久久久黄片| 99国产精品一区二区蜜桃av| 国产成人aa在线观看| 免费观看精品视频网站| 国产精品野战在线观看| 久久久久久久久久久免费av| 国产一区二区在线观看日韩| 综合色av麻豆| 99久久精品国产国产毛片| 午夜久久久久精精品| 欧美日韩综合久久久久久| 国产不卡一卡二| 丰满人妻一区二区三区视频av| 99国产精品一区二区蜜桃av| 亚洲不卡免费看| 国产老妇女一区| 黑人高潮一二区| 成人av在线播放网站| 免费观看精品视频网站| 成人毛片a级毛片在线播放| 国产av不卡久久| 国产真实乱freesex| 国产精品不卡视频一区二区| 精品一区二区免费观看| 久久6这里有精品| 国产精品蜜桃在线观看| 人人妻人人澡欧美一区二区| 免费观看的影片在线观看| 伊人久久精品亚洲午夜| 一级二级三级毛片免费看| 日韩欧美国产在线观看| 26uuu在线亚洲综合色| 久久国内精品自在自线图片| 人人妻人人澡人人爽人人夜夜 | 国产成年人精品一区二区| 国产精品伦人一区二区| 欧美日韩在线观看h| 最近的中文字幕免费完整| 看免费成人av毛片| 我要看日韩黄色一级片| 久久亚洲国产成人精品v| 日韩三级伦理在线观看| 十八禁国产超污无遮挡网站| 亚洲精品,欧美精品| 男女视频在线观看网站免费| 午夜福利成人在线免费观看| 九九在线视频观看精品| 亚洲av熟女| 午夜久久久久精精品| 国产精品综合久久久久久久免费| 1024手机看黄色片| 国产免费福利视频在线观看| 国产在视频线精品| 高清午夜精品一区二区三区| 一本久久精品| 久久精品91蜜桃| 色综合站精品国产| 99热精品在线国产| 男女边吃奶边做爰视频| 欧美变态另类bdsm刘玥| 国产欧美另类精品又又久久亚洲欧美| 成人美女网站在线观看视频| 男女下面进入的视频免费午夜| 精品无人区乱码1区二区| 纵有疾风起免费观看全集完整版 | 国产女主播在线喷水免费视频网站 | 欧美精品一区二区大全| 日韩成人av中文字幕在线观看| 国产淫语在线视频| 精品午夜福利在线看| 18禁动态无遮挡网站| 日本免费在线观看一区| 国产精华一区二区三区| 熟女电影av网| 国产亚洲av嫩草精品影院| 九九在线视频观看精品| 一本久久精品| 亚洲av二区三区四区| 日本免费a在线| 青春草视频在线免费观看| 最近最新中文字幕免费大全7| 亚洲av不卡在线观看| 国产精品福利在线免费观看| 建设人人有责人人尽责人人享有的 | 国国产精品蜜臀av免费| 99热精品在线国产| 精品久久国产蜜桃| 日韩在线高清观看一区二区三区| 淫秽高清视频在线观看| 亚洲性久久影院| 美女xxoo啪啪120秒动态图| 成人一区二区视频在线观看| 国产一区亚洲一区在线观看| 老女人水多毛片| 国产精品一及| 国产免费又黄又爽又色| 国产av不卡久久| 啦啦啦观看免费观看视频高清| 亚洲国产欧美在线一区| 可以在线观看毛片的网站| 变态另类丝袜制服| 中文字幕av在线有码专区| 精品久久国产蜜桃| 18+在线观看网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 69av精品久久久久久| 人妻制服诱惑在线中文字幕| 色播亚洲综合网| 亚洲成色77777| 久久午夜福利片| 九九在线视频观看精品| 好男人在线观看高清免费视频| 永久网站在线| 国产黄色视频一区二区在线观看 | 亚洲久久久久久中文字幕| 少妇人妻一区二区三区视频| 欧美另类亚洲清纯唯美| 最新中文字幕久久久久| 精品人妻一区二区三区麻豆| 国产精品av视频在线免费观看| av线在线观看网站| 亚洲av.av天堂| 国产av一区在线观看免费| 国产精品久久久久久精品电影| 国产在视频线在精品| 亚洲美女视频黄频| 级片在线观看| 麻豆国产97在线/欧美| 精品久久久久久久人妻蜜臀av| 国产高清国产精品国产三级 | 精品久久久噜噜| 小蜜桃在线观看免费完整版高清| 国产伦精品一区二区三区视频9| 久久99热这里只频精品6学生 | 国产伦一二天堂av在线观看| 国产精品三级大全| 成人二区视频| 大又大粗又爽又黄少妇毛片口| 国产av码专区亚洲av| 欧美一区二区亚洲| 亚洲国产成人一精品久久久| 欧美日韩综合久久久久久| 一级毛片电影观看 | 亚洲欧美日韩东京热| 日韩欧美三级三区| 99在线视频只有这里精品首页| 成人高潮视频无遮挡免费网站| 长腿黑丝高跟| 午夜福利成人在线免费观看| 中文字幕av成人在线电影| 日韩av在线大香蕉| 欧美变态另类bdsm刘玥| 成人国产麻豆网| 亚洲欧美日韩高清专用| 亚洲最大成人中文| 国产黄a三级三级三级人| 最近最新中文字幕大全电影3| 午夜a级毛片| 乱码一卡2卡4卡精品| 久久精品影院6| 午夜久久久久精精品| 亚洲国产欧美在线一区| 久久亚洲国产成人精品v| 日本爱情动作片www.在线观看| 看十八女毛片水多多多| 天天躁日日操中文字幕| 国产黄片美女视频| 国产老妇伦熟女老妇高清| 人妻夜夜爽99麻豆av| 麻豆乱淫一区二区| 国产 一区精品| 精品久久久久久成人av| 国产成年人精品一区二区| 久久精品熟女亚洲av麻豆精品 | 最近中文字幕高清免费大全6| 国语自产精品视频在线第100页| 午夜老司机福利剧场| 1000部很黄的大片| 免费在线观看成人毛片| 午夜亚洲福利在线播放| 老师上课跳d突然被开到最大视频| 精品午夜福利在线看| 看黄色毛片网站| 免费大片18禁| 女人久久www免费人成看片 | 日日啪夜夜撸| 中文字幕免费在线视频6| 国产精品无大码| 亚洲人成网站在线播| 久久久久久伊人网av| 青春草视频在线免费观看| 舔av片在线| 国产精品国产高清国产av| 亚洲综合精品二区| 波多野结衣高清无吗| 建设人人有责人人尽责人人享有的 | 免费观看性生交大片5| 久久久久久久久久久免费av| 日韩亚洲欧美综合| 亚洲成人精品中文字幕电影| 精品久久久久久久久久久久久| 婷婷六月久久综合丁香| 欧美丝袜亚洲另类| 亚洲国产日韩欧美精品在线观看| h日本视频在线播放| 男人舔女人下体高潮全视频| 午夜精品国产一区二区电影 | 色吧在线观看| 热99re8久久精品国产| 男的添女的下面高潮视频| 狠狠狠狠99中文字幕| 日本三级黄在线观看| 精品久久久久久成人av| av播播在线观看一区| 日韩 亚洲 欧美在线| 欧美日韩综合久久久久久| 一二三四中文在线观看免费高清| 亚洲欧美成人精品一区二区| 九九热线精品视视频播放| 亚洲av中文字字幕乱码综合| 国产成人福利小说| 亚洲丝袜综合中文字幕| 亚洲精品色激情综合| 乱码一卡2卡4卡精品| 亚洲av男天堂| 国产色爽女视频免费观看| 在线免费观看的www视频| 乱码一卡2卡4卡精品| 国产成人免费观看mmmm| 欧美另类亚洲清纯唯美| 色噜噜av男人的天堂激情| 干丝袜人妻中文字幕| 亚洲精品久久久久久婷婷小说 | 99久久精品一区二区三区| 色哟哟·www| 免费一级毛片在线播放高清视频| 97超视频在线观看视频| 亚洲精品一区蜜桃| 日韩大片免费观看网站 | 亚洲三级黄色毛片| 久久精品影院6| 麻豆一二三区av精品| 国产高清有码在线观看视频| av.在线天堂| 国产视频内射| 免费电影在线观看免费观看| 色5月婷婷丁香| 一级毛片我不卡| 日本五十路高清| 中文字幕制服av| 国产精品,欧美在线| 91精品一卡2卡3卡4卡| 天天躁日日操中文字幕| 22中文网久久字幕| 国产精品美女特级片免费视频播放器| 亚洲成人中文字幕在线播放| 黄色欧美视频在线观看| 嫩草影院新地址| 小说图片视频综合网站| 久久久久性生活片| 精品久久久久久久久亚洲| 午夜精品一区二区三区免费看| 一级二级三级毛片免费看| 人体艺术视频欧美日本| 少妇人妻精品综合一区二区| 久久久成人免费电影| 午夜免费男女啪啪视频观看| 国模一区二区三区四区视频| 国产精品国产三级专区第一集| 亚洲精品乱码久久久v下载方式| 久久久久久国产a免费观看| 看十八女毛片水多多多| 亚洲人与动物交配视频| 久久亚洲精品不卡| 18+在线观看网站| 一区二区三区免费毛片| 自拍偷自拍亚洲精品老妇| 亚洲人成网站在线播| 成人特级av手机在线观看| .国产精品久久| 99热6这里只有精品| 国产精品久久久久久精品电影| 亚洲精华国产精华液的使用体验| 一本久久精品| 91aial.com中文字幕在线观看| 亚洲人成网站高清观看| 美女xxoo啪啪120秒动态图| 天堂中文最新版在线下载 | 女人十人毛片免费观看3o分钟| 你懂的网址亚洲精品在线观看 | 久久久成人免费电影| 国产精品不卡视频一区二区| 亚洲va在线va天堂va国产| 91久久精品电影网| 亚洲国产精品久久男人天堂| 成人特级av手机在线观看| av天堂中文字幕网| 一级黄片播放器| 美女黄网站色视频| 日韩国内少妇激情av| 午夜爱爱视频在线播放| 亚洲怡红院男人天堂| 色5月婷婷丁香| 日韩亚洲欧美综合| 久久精品久久久久久噜噜老黄 | 毛片一级片免费看久久久久| 神马国产精品三级电影在线观看| 国产亚洲av片在线观看秒播厂 | 插阴视频在线观看视频| 国产一区二区三区av在线| 亚洲在线观看片| 91久久精品国产一区二区成人| 日韩av在线免费看完整版不卡| 国产乱人偷精品视频| 舔av片在线| 国产免费福利视频在线观看| 国产一级毛片在线| 精品久久国产蜜桃| 国产av码专区亚洲av| 男人狂女人下面高潮的视频| 男女边吃奶边做爰视频| 精品不卡国产一区二区三区| 久久99热6这里只有精品| 精品午夜福利在线看| 国产黄片美女视频| 韩国高清视频一区二区三区| 国产亚洲91精品色在线| av福利片在线观看| 99热全是精品| 亚洲国产欧美在线一区| 人人妻人人澡人人爽人人夜夜 | 又爽又黄a免费视频| 国产精品爽爽va在线观看网站| 久久精品91蜜桃| 99热这里只有精品一区| 我要看日韩黄色一级片| 成人毛片60女人毛片免费|