• <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| 99热网站在线观看| 精品熟女少妇av免费看| 久久久久精品性色| 肉色欧美久久久久久久蜜桃| 99视频精品全部免费 在线| 80岁老熟妇乱子伦牲交| 久久亚洲国产成人精品v| 熟女av电影| 女人久久www免费人成看片| 国产av码专区亚洲av| 欧美精品一区二区大全| 在线观看免费日韩欧美大片 | 一二三四中文在线观看免费高清| 18禁在线无遮挡免费观看视频| 亚洲,一卡二卡三卡| 欧美日韩国产mv在线观看视频| av女优亚洲男人天堂| 狂野欧美激情性bbbbbb| 国产爽快片一区二区三区| 99久久精品一区二区三区| 熟女电影av网| 亚洲av男天堂| 亚洲国产最新在线播放| 久久久久久久亚洲中文字幕| 女人久久www免费人成看片| 一区二区三区乱码不卡18| 少妇 在线观看| 欧美一级a爱片免费观看看| 国产在线视频一区二区| 草草在线视频免费看| av在线播放精品| 久久久久国产精品人妻一区二区| 亚洲人成网站在线观看播放| 日韩av在线免费看完整版不卡| 亚州av有码| 日本爱情动作片www.在线观看| 久久精品夜色国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一本大道久久a久久精品| 老司机影院成人| 国产在线一区二区三区精| 久久久a久久爽久久v久久| 免费高清在线观看日韩| 精品午夜福利在线看| 日韩欧美精品免费久久| 久久久a久久爽久久v久久| 夜夜看夜夜爽夜夜摸| 桃花免费在线播放| 久久久久网色| 欧美激情 高清一区二区三区| 久久热精品热| 中国国产av一级| 国产视频内射| 黄色配什么色好看| 日韩精品有码人妻一区| 国产欧美另类精品又又久久亚洲欧美| 午夜久久久在线观看| 久久人人爽av亚洲精品天堂| 国产成人av激情在线播放 | 亚洲精品中文字幕在线视频| 热re99久久国产66热| videossex国产| a级毛片免费高清观看在线播放| 少妇的逼水好多| 精品国产乱码久久久久久小说| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 国产成人aa在线观看| 嘟嘟电影网在线观看| 国产精品麻豆人妻色哟哟久久| 97超碰精品成人国产| 国产永久视频网站| 在线观看三级黄色| 婷婷色综合www| 麻豆精品久久久久久蜜桃| 最近2019中文字幕mv第一页| 在线观看一区二区三区激情| 亚洲中文av在线| 成人手机av| 国产精品人妻久久久影院| 观看美女的网站| 人妻少妇偷人精品九色| 高清在线视频一区二区三区| 国产毛片在线视频| 亚洲精品亚洲一区二区| 色吧在线观看| 观看美女的网站| 国产精品久久久久成人av| 少妇熟女欧美另类| 在线观看一区二区三区激情| 看十八女毛片水多多多| 欧美亚洲 丝袜 人妻 在线| 亚洲精品乱久久久久久| 亚洲欧美日韩卡通动漫| 免费观看的影片在线观看| 中文字幕免费在线视频6| 在线播放无遮挡| 免费观看av网站的网址| 又大又黄又爽视频免费| 日韩大片免费观看网站| 日日摸夜夜添夜夜爱| 九色亚洲精品在线播放| 九九爱精品视频在线观看| 一级毛片aaaaaa免费看小| 日本黄色日本黄色录像| 亚洲精品久久久久久婷婷小说| 26uuu在线亚洲综合色| 成人漫画全彩无遮挡| av福利片在线| 中文字幕免费在线视频6| 天堂俺去俺来也www色官网| 国产乱人偷精品视频| 日日摸夜夜添夜夜添av毛片| 男女边摸边吃奶| 亚洲av男天堂| 夫妻午夜视频| 男女高潮啪啪啪动态图| 国产精品久久久久成人av| 欧美日韩亚洲高清精品| 日韩在线高清观看一区二区三区| 欧美最新免费一区二区三区| 视频在线观看一区二区三区| www.色视频.com| 国产精品欧美亚洲77777| 国产综合精华液| 美女cb高潮喷水在线观看| 日韩亚洲欧美综合| 三级国产精品片| 丝袜喷水一区| 国产精品99久久久久久久久| 国产熟女欧美一区二区| 永久网站在线| kizo精华| 91久久精品电影网| 国产熟女午夜一区二区三区 | 国产精品熟女久久久久浪| 99热这里只有精品一区| 少妇被粗大的猛进出69影院 | 日韩精品免费视频一区二区三区 | 十八禁高潮呻吟视频| 久久久久国产网址| 国产一区二区三区av在线| 久久久精品94久久精品| 国产免费一区二区三区四区乱码| 麻豆精品久久久久久蜜桃| h视频一区二区三区| 亚洲第一区二区三区不卡| 制服诱惑二区| 亚洲欧洲精品一区二区精品久久久 | 人成视频在线观看免费观看| 欧美人与性动交α欧美精品济南到 | 中文字幕久久专区| 一级毛片黄色毛片免费观看视频| 亚洲欧美一区二区三区国产| 国产伦理片在线播放av一区| 日本与韩国留学比较| 九九久久精品国产亚洲av麻豆| 9色porny在线观看| 在线免费观看不下载黄p国产| 国产亚洲av片在线观看秒播厂| 另类亚洲欧美激情| 亚洲熟女精品中文字幕| 91久久精品国产一区二区成人| 国产欧美日韩一区二区三区在线 | 亚洲欧洲日产国产| 我要看黄色一级片免费的| 青春草视频在线免费观看| 日本-黄色视频高清免费观看| 国产成人91sexporn| 啦啦啦中文免费视频观看日本| 在线观看免费高清a一片| 美女视频免费永久观看网站| 国产精品.久久久| freevideosex欧美| 国产男人的电影天堂91| 久久精品久久精品一区二区三区| av女优亚洲男人天堂| 男人添女人高潮全过程视频| 永久网站在线| 亚洲精品日韩在线中文字幕| 黑人高潮一二区| av国产久精品久网站免费入址| a级毛色黄片| 高清视频免费观看一区二区| 日本午夜av视频| 久久久久人妻精品一区果冻| 熟妇人妻不卡中文字幕| 热re99久久精品国产66热6| 午夜免费男女啪啪视频观看| 人妻夜夜爽99麻豆av| 欧美精品国产亚洲| 午夜免费鲁丝| a级毛色黄片| 少妇被粗大的猛进出69影院 | 精品视频人人做人人爽| 夜夜看夜夜爽夜夜摸| 亚洲综合色惰| 国产av国产精品国产| 国产精品久久久久成人av| 女的被弄到高潮叫床怎么办| 久久久久久久久久成人| 纵有疾风起免费观看全集完整版| av免费在线看不卡| 精品久久久精品久久久| 国产成人精品一,二区| 中文字幕人妻丝袜制服| 女人精品久久久久毛片| 一区二区三区免费毛片| 精品99又大又爽又粗少妇毛片| 国产成人一区二区在线| 永久免费av网站大全| 国产欧美另类精品又又久久亚洲欧美| av卡一久久| 亚洲av综合色区一区| 黄色配什么色好看| 亚洲av欧美aⅴ国产| 成人影院久久| 香蕉精品网在线| 国产成人午夜福利电影在线观看| 大陆偷拍与自拍| 日产精品乱码卡一卡2卡三| 五月玫瑰六月丁香| 亚洲av免费高清在线观看| 高清毛片免费看| 99国产精品免费福利视频| 久久99精品国语久久久| 一级毛片电影观看| 亚洲国产精品成人久久小说| 国产淫语在线视频| 亚洲精品aⅴ在线观看| 亚洲精品色激情综合| 亚洲精品av麻豆狂野| 免费看不卡的av| 能在线免费看毛片的网站| 少妇人妻精品综合一区二区| 新久久久久国产一级毛片| 99国产综合亚洲精品| 亚洲一区二区三区欧美精品| 草草在线视频免费看| 久久韩国三级中文字幕| 狠狠精品人妻久久久久久综合| 日韩电影二区| 亚洲精品日韩在线中文字幕| 亚洲精品亚洲一区二区| 国产一区亚洲一区在线观看| 亚洲综合色网址| 日本爱情动作片www.在线观看| 在线观看美女被高潮喷水网站| 99久久综合免费| 午夜91福利影院| 亚洲精品第二区| 亚洲图色成人| 高清毛片免费看| 成年女人在线观看亚洲视频| 久久国产精品大桥未久av| 五月伊人婷婷丁香| 国产精品女同一区二区软件| 亚洲人与动物交配视频| 亚洲精品久久成人aⅴ小说 | 久久国产精品大桥未久av| 久久精品国产自在天天线| 多毛熟女@视频| 老司机亚洲免费影院| 国产高清三级在线| 伊人久久精品亚洲午夜| 亚洲精品aⅴ在线观看| 日日撸夜夜添| 久久人人爽av亚洲精品天堂| 国产精品久久久久久精品古装| 波野结衣二区三区在线| 看非洲黑人一级黄片| 最近手机中文字幕大全| 青春草视频在线免费观看| 国产精品无大码| 国产免费视频播放在线视频| av在线老鸭窝| 狠狠婷婷综合久久久久久88av| 国产男女内射视频| 久久精品久久久久久久性| 青春草视频在线免费观看| 日产精品乱码卡一卡2卡三| 99九九在线精品视频| 91精品一卡2卡3卡4卡| 日韩大片免费观看网站| a级毛色黄片| 免费av不卡在线播放| 女人精品久久久久毛片| 亚洲欧美日韩另类电影网站| www.av在线官网国产| 黑人高潮一二区| 久久99热6这里只有精品| 亚洲av二区三区四区| 在线观看三级黄色| 丰满饥渴人妻一区二区三| 人体艺术视频欧美日本| 国产精品免费大片| 亚洲国产精品成人久久小说| 精品人妻在线不人妻| av播播在线观看一区| 在线观看免费高清a一片| 蜜桃久久精品国产亚洲av| 狠狠婷婷综合久久久久久88av| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 亚洲av成人精品一区久久| 人体艺术视频欧美日本| 黑人高潮一二区| videossex国产| 少妇人妻精品综合一区二区| 午夜影院在线不卡| 亚洲伊人久久精品综合| 国产午夜精品一二区理论片| 日韩av不卡免费在线播放| 久久精品久久精品一区二区三区| 少妇人妻精品综合一区二区| av视频免费观看在线观看| 亚洲精品久久久久久婷婷小说| 精品久久国产蜜桃| 人妻少妇偷人精品九色| videos熟女内射| 人体艺术视频欧美日本| 七月丁香在线播放| 久久国内精品自在自线图片| 午夜影院在线不卡| 18在线观看网站| 性色avwww在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲精品久久久久久婷婷小说| 欧美性感艳星| 国产亚洲一区二区精品| videosex国产| 亚州av有码| 热99国产精品久久久久久7| 国产不卡av网站在线观看| 亚洲怡红院男人天堂| 99久国产av精品国产电影| 国产不卡av网站在线观看| 日韩精品有码人妻一区| 精品卡一卡二卡四卡免费| 国产色爽女视频免费观看| 看十八女毛片水多多多| 亚洲av成人精品一二三区| 2021少妇久久久久久久久久久| 国产一区有黄有色的免费视频| 国产av精品麻豆| 国产欧美另类精品又又久久亚洲欧美| 日韩亚洲欧美综合| 中文字幕久久专区| 日韩人妻高清精品专区| 欧美另类一区| 在线观看www视频免费| 国产精品一区二区三区四区免费观看| 日韩一本色道免费dvd| 精品久久蜜臀av无| 精品少妇内射三级| 纯流量卡能插随身wifi吗| 久久av网站| 日本av免费视频播放| 久久97久久精品| 观看美女的网站| 国产精品嫩草影院av在线观看| 欧美97在线视频| 高清黄色对白视频在线免费看| 大码成人一级视频| 亚洲精品日本国产第一区| 少妇被粗大的猛进出69影院 | 建设人人有责人人尽责人人享有的| 尾随美女入室| 国产午夜精品久久久久久一区二区三区| 国产成人av激情在线播放 | 午夜激情福利司机影院| 亚洲国产毛片av蜜桃av| 十八禁网站网址无遮挡| 久久久久视频综合| 少妇猛男粗大的猛烈进出视频| 午夜精品国产一区二区电影| 久久久欧美国产精品| 天天操日日干夜夜撸| 亚洲国产欧美在线一区| 91精品国产国语对白视频| 美女视频免费永久观看网站| 黑丝袜美女国产一区| 一级毛片 在线播放| 黄色视频在线播放观看不卡| 欧美变态另类bdsm刘玥| 一区在线观看完整版| 久久国内精品自在自线图片| 天天影视国产精品| 视频在线观看一区二区三区| 美女福利国产在线| 欧美精品高潮呻吟av久久| freevideosex欧美| 日日爽夜夜爽网站| 日韩精品免费视频一区二区三区 | 久久精品国产亚洲av天美| 99精国产麻豆久久婷婷| av一本久久久久| 久久精品国产亚洲av天美| 九九在线视频观看精品| 69精品国产乱码久久久| 老司机影院成人| 日韩欧美一区视频在线观看| 男女啪啪激烈高潮av片| 国产精品三级大全| 99久国产av精品国产电影| 精品人妻一区二区三区麻豆| 永久网站在线| xxxhd国产人妻xxx| 18禁在线播放成人免费| 99热6这里只有精品| 99re6热这里在线精品视频| 日韩中字成人| 精品少妇内射三级| 亚洲第一av免费看| 视频在线观看一区二区三区| 天美传媒精品一区二区| 99热这里只有是精品在线观看| 亚洲,欧美,日韩| 欧美日韩视频精品一区| 最近中文字幕高清免费大全6| 亚洲精品国产av蜜桃| 久久久国产精品麻豆| 一级毛片 在线播放| 久久久精品区二区三区| 九草在线视频观看| 搡女人真爽免费视频火全软件| 人人妻人人添人人爽欧美一区卜| 亚洲欧美日韩另类电影网站| 久久精品国产亚洲av涩爱| 成人国语在线视频| 性色av一级| 高清午夜精品一区二区三区| 黑人欧美特级aaaaaa片| 精品人妻一区二区三区麻豆| 国产极品天堂在线| 精品国产国语对白av| 亚洲无线观看免费| 少妇人妻久久综合中文| 一边摸一边做爽爽视频免费| videosex国产| xxx大片免费视频| 麻豆乱淫一区二区| 免费大片黄手机在线观看| 我的女老师完整版在线观看| 亚洲欧美中文字幕日韩二区| 国产成人freesex在线| 亚洲怡红院男人天堂| 欧美 亚洲 国产 日韩一| 99国产精品免费福利视频| 桃花免费在线播放| 免费观看的影片在线观看| 国产日韩欧美亚洲二区| 97精品久久久久久久久久精品| 日韩欧美一区视频在线观看| 一级毛片 在线播放| 青春草国产在线视频| 国产毛片在线视频| 国产精品国产三级国产av玫瑰| 婷婷色综合www| 日韩精品免费视频一区二区三区 | 九九爱精品视频在线观看| 亚洲成色77777| 精品久久蜜臀av无| av黄色大香蕉| 国产有黄有色有爽视频| 国产精品免费大片| 亚洲精品中文字幕在线视频| 日本-黄色视频高清免费观看| 亚洲av电影在线观看一区二区三区| 18禁观看日本| 久久久亚洲精品成人影院| 91久久精品国产一区二区成人| 中文字幕久久专区| 最新中文字幕久久久久| 亚洲国产色片| 成年人免费黄色播放视频| 精品久久久久久久久av| 久久久久久久国产电影| 美女大奶头黄色视频| 丝袜美足系列| 久久精品国产亚洲av涩爱| 有码 亚洲区| 伊人久久国产一区二区| 少妇人妻 视频| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区三区在线 | 午夜福利视频精品| av视频免费观看在线观看| 亚洲欧美成人综合另类久久久| 精品一区二区三区视频在线| 人妻一区二区av| 爱豆传媒免费全集在线观看| 久久精品久久久久久久性| 日本-黄色视频高清免费观看| 久久久久人妻精品一区果冻| 精品亚洲成国产av| 制服诱惑二区| 欧美日韩av久久| 亚洲成人av在线免费| 亚洲情色 制服丝袜| 国产免费又黄又爽又色| 视频区图区小说| 韩国av在线不卡| 精品国产一区二区三区久久久樱花| 伦理电影免费视频| 国产日韩欧美在线精品| 国产国语露脸激情在线看| av女优亚洲男人天堂| 国产精品.久久久| 美女xxoo啪啪120秒动态图| 在线观看人妻少妇| 又大又黄又爽视频免费| a级毛色黄片| 成人18禁高潮啪啪吃奶动态图 | 哪个播放器可以免费观看大片| 国模一区二区三区四区视频| 91精品一卡2卡3卡4卡| 少妇的逼好多水| 91成人精品电影| 看非洲黑人一级黄片| 久久精品久久久久久久性| 汤姆久久久久久久影院中文字幕| 成年人午夜在线观看视频| 亚洲国产精品一区二区三区在线| av卡一久久| 亚洲av二区三区四区| av在线播放精品| 国产成人精品福利久久| 黄色毛片三级朝国网站| 国产白丝娇喘喷水9色精品| 卡戴珊不雅视频在线播放| 国产男女超爽视频在线观看| 插阴视频在线观看视频| 亚洲av日韩在线播放| 亚洲人与动物交配视频| 国产黄色免费在线视频| 久久久精品区二区三区| 最黄视频免费看| 伊人久久国产一区二区| 精品久久蜜臀av无| 久久这里有精品视频免费| 久久久久人妻精品一区果冻| 免费大片黄手机在线观看| 欧美日韩国产mv在线观看视频| 肉色欧美久久久久久久蜜桃| 午夜激情av网站| 91国产中文字幕| 熟妇人妻不卡中文字幕| 99久久中文字幕三级久久日本| 乱人伦中国视频| 日本vs欧美在线观看视频| 夜夜爽夜夜爽视频| 国产精品免费大片| 日韩欧美一区视频在线观看| 精品一区在线观看国产| 赤兔流量卡办理| 各种免费的搞黄视频| 免费黄频网站在线观看国产| 日韩,欧美,国产一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 亚洲av在线观看美女高潮| 国产成人91sexporn| 欧美97在线视频| 黄片播放在线免费| 久久久久久久久久成人| 91久久精品国产一区二区三区| 美女中出高潮动态图| 久久人人爽人人爽人人片va| 97在线人人人人妻| 日本vs欧美在线观看视频| 黄色欧美视频在线观看| 欧美少妇被猛烈插入视频| 考比视频在线观看| 午夜免费鲁丝| 99视频精品全部免费 在线| 国产老妇伦熟女老妇高清| 一级,二级,三级黄色视频| 大片免费播放器 马上看| 满18在线观看网站| a 毛片基地| 日本av免费视频播放| 精品亚洲成a人片在线观看| 色哟哟·www| 另类亚洲欧美激情| 人人妻人人添人人爽欧美一区卜| 亚洲性久久影院| 人人妻人人添人人爽欧美一区卜| 熟女av电影| 国产精品女同一区二区软件| 久久久欧美国产精品| 国产亚洲午夜精品一区二区久久| 国产一区有黄有色的免费视频| 亚洲综合色惰| 黑人巨大精品欧美一区二区蜜桃 | 国产成人精品久久久久久| 国产精品久久久久久精品电影小说| 交换朋友夫妻互换小说| 日本黄色片子视频| 久久99精品国语久久久| 国产精品蜜桃在线观看| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 大香蕉久久成人网| 天堂俺去俺来也www色官网| 中文欧美无线码| 欧美 日韩 精品 国产| 久久精品久久久久久噜噜老黄| 国产精品一二三区在线看| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 色婷婷av一区二区三区视频| 91成人精品电影| 七月丁香在线播放| 免费少妇av软件|