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

    瑪曲高寒草甸放牧強度的遙感監(jiān)測

    2018-05-26 07:05:02馬青青柴林榮馬海玲常生華程云湘
    草業(yè)科學(xué) 2018年5期
    關(guān)鍵詞:牧戶草甸植被指數(shù)

    馬青青,柴林榮,馬海玲,常生華,程云湘

    (1.蘭州大學(xué)草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點實驗室 蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院,甘肅 蘭州 730020;2.蘭州大學(xué)農(nóng)業(yè)農(nóng)村部草牧業(yè)創(chuàng)新重點實驗室,甘肅 蘭州 730020; 3.青海省氣象災(zāi)害防御技術(shù)中心,青海 西寧 810001)

    青藏高原生態(tài)系統(tǒng)對全球氣候變化較為敏感[1],高原生態(tài)系統(tǒng)的相關(guān)研究在全球變化研究中占據(jù)極為重要的地位[2]。高寒草甸生態(tài)系統(tǒng)是青藏高原最主要的生態(tài)系統(tǒng)類型,約占青藏高原總面積的35%,在高原草地畜牧業(yè)中占有十分重要的地位[3]。放牧活動是高寒草甸的主要利用方式之一。多數(shù)研究表明,適度的放牧可以有效促進草地生產(chǎn)力的提高和生物多樣性的維持,但過度放牧往往會引起草地的大面積退化以及草地初級生產(chǎn)力和土壤理化性狀的惡化[4]。近年來,由于長期受超載過牧等因素的影響,高寒草甸生態(tài)系統(tǒng)退化顯著[5],從20世紀90年代開始,尤其青藏高原東緣草地出現(xiàn)了明顯的退化,其基本特征為區(qū)域氣候趨于干旱、地表徑流量顯著減少、水土流失日益嚴重[6],草畜矛盾日益尖銳。

    草地生物量是評價青藏高原生態(tài)脆弱性和敏感性的重要指標(biāo)及度量草地退化的重要指標(biāo)[7],也是整個草地生態(tài)系統(tǒng)運行的能量基礎(chǔ)和物質(zhì)來源,不僅是體現(xiàn)草地生態(tài)系統(tǒng)穩(wěn)定性和種群多樣性的指標(biāo),同時也被用于指導(dǎo)畜牧業(yè)生產(chǎn)[8]。對生物量的定量研究,可以及時準(zhǔn)確地獲取草地信息,為科學(xué)合理確定載畜量和有效進行草地資源保護提供參考依據(jù),同時也為草地生態(tài)環(huán)境研究提供基礎(chǔ)數(shù)據(jù)。草地生物量是草地資源遙感研究的主要指標(biāo)之一,監(jiān)測草地生物量的時空變化動態(tài)能夠較好地反映草地的生長狀況,在除放牧外無其他干擾因素的情況下能夠指示放牧強度的大小,對牧區(qū)草地退化評價和草畜平衡決策研究具有重要意義[9]。

    遙感數(shù)據(jù)已被廣泛應(yīng)用于大尺度植被活動的研究[10-12],其在實現(xiàn)草地估產(chǎn)、確定草地載畜量以及加強草地科學(xué)管理等方面都有重要的指導(dǎo)意義[13]。Paruelo等[14]在美國中部草原地區(qū)將NOAA/AVHRR遙感數(shù)據(jù)用于草地資源動態(tài)監(jiān)測,研究表明,歸一化植被指數(shù)與草地地上生物量高度相關(guān),且發(fā)現(xiàn)冪函數(shù)回歸模型是該區(qū)的最佳遙感反演模型。Bénié等[15]將遙感數(shù)據(jù)和社會經(jīng)濟因素加到模型中,以至少80%的準(zhǔn)確性模擬了薩赫勒半干旱地區(qū)十多年的生物量,因其適應(yīng)性強的模擬能力減少了數(shù)據(jù)采集成本,且有助于薩赫勒環(huán)境的管理。韓波[16]結(jié)合三江源區(qū)高寒草地地上生物量(aboveground biomass,AGB)野外實測數(shù)據(jù)和同期MODIS-NDVI(normalized difference vegetation index)和MODIS-EVI(enhanced vegetation index)數(shù)據(jù),通過回歸分析,建立研究區(qū)高寒草地AGB遙感反演模型,并對三江源區(qū)近十年高寒草地AGB時空格局進行模擬分析。周鑫等[17]分析了EVI和NDVI與實測草地地上生物量的一元線性、指數(shù)和二次多項式回歸模型,并對各種回歸模型進行分析比較,利用優(yōu)選模型的反演結(jié)果分析了伊犁地區(qū)草地地上生物量的空間分布。綜合上述研究,采用MODIS植被指數(shù)能夠監(jiān)測和預(yù)測草地產(chǎn)量,只是不同的地區(qū)會因?qū)嶋H情況采用不同的指數(shù)或模型。但是,至今為止,很少有研究直接涉及放牧強度和植被指數(shù)之間的關(guān)系。張艷楠等[18]利用NDVI和生物量建立的生物量估產(chǎn)模型計算出34家牧戶的7、8月生物量的差值,以放牧強度為應(yīng)變量,差值為自變量,通過回歸分析建立了放牧強度估測模型(R2= 0.461),發(fā)現(xiàn)二者之間具有較好的相關(guān)性。該研究中還是以生物量作為橋梁,若是直接利用NDVI和放牧強度來建立大范圍放牧強度估測模型,一方面突破了放牧強度量化的區(qū)域限制,另一方面可以對放牧行為實時監(jiān)控,調(diào)整家畜活動范圍,從一定程度上緩解草地退化壓力。

    本研究首先利用野外62個草地調(diào)查樣點實測地上生物量數(shù)據(jù),建立了研究區(qū)植被指數(shù)與草地地上生物量的回歸模型,并結(jié)合實測地上生物量和不同相對放牧強度間的相關(guān)性,分析了不同放牧強度條件下草地生物量的空間分布特征。然后利用3家牧戶的相對放牧強度,結(jié)合同期MODIS影像獲取的植被指數(shù)數(shù)據(jù),將相對放牧強度和對應(yīng)的NDVI在區(qū)域空間尺度上建立回歸模型,獲取研究區(qū)植被景觀格局和家畜放牧利用強度的時空分布關(guān)系,以期為研究區(qū)草地資源的可持續(xù)利用和草畜平衡提供依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)設(shè)在甘肅省瑪曲縣阿孜畜牧科技示范園區(qū)(圖1),地理位置為33°40′ N,101°52′ E,海拔3 650 m,年均氣溫1.2 ℃,年均降水量620 mm,主要集中在5-9月,年日照時數(shù)約2 580 h,年平均霜日大于270 d,無絕對無霜期,只有冷暖季之分。土壤為亞高山草甸土。草地類型為高寒草甸,植被以嵩草(Kobresiasp.)、苔草(Carexsp.)、針茅(Stipasp.)、羊茅(Festucasp.)、披堿草(Elymussp.)和早熟禾(Poasp.)等為主。牧草一般在4月中、下旬開始返青,8月中、下旬開始枯黃。放牧家畜主要是地方藏系牦牛。

    1.2 地面實測數(shù)據(jù)

    試驗樣地主要設(shè)置在阿孜畜牧科技示范園區(qū)的牧戶地,選擇草地類型一致(高寒草甸為主)、生產(chǎn)和利用方式相似(以牦牛放牧為主)的3家牧戶的草地作為研究樣地,而且這些草地都處于同一生長水平(即草地的理論可食牧草產(chǎn)量相等,且生長周期一致,除放牧外無其他干擾因素)。根據(jù)牧戶牦牛數(shù)量、草地面積、距離牧戶棚圈的遠近及牦牛選擇攝食的活動路線和在空間停留頻次時所表現(xiàn)出的放牧活動強度,將牧戶草地劃分為3個放牧梯度區(qū),即距離畜圈入口較近區(qū)域為第一放牧梯度區(qū),距離畜圈入口中等距離區(qū)域為第二放牧梯度區(qū),距離畜圈入口較遠區(qū)域為第三放牧梯度區(qū)(圖2)。在每個放牧梯度區(qū)隨機選擇10個0.5 m×0.5 m的樣方,且采樣點間距符合遙感數(shù)據(jù)空間分辨率的要求,每個放牧梯度區(qū)重復(fù)4次,在6、7、8、9月的下旬齊地刈割,進行草地植被調(diào)查。記錄樣方內(nèi)主要物種的名稱、蓋度、高度等植物群落信息,然后齊地面刈割,稱其鮮重,再在70 ℃烘至恒重稱其干重。

    圖1 瑪曲縣阿孜站的地理位置Fig. 1 Location of Aze station in Maqu County

    圖2 試驗樣地Fig. 2 Sample plot

    1.3 植被指數(shù)的提取和相對放牧強度的計算

    歸一化植被指數(shù)(NDVI)通常是紅波段(R)和近紅外波段(NIR)反射率之差與之和的比值,用以表征地表植被的數(shù)量分配和質(zhì)量情況,即:

    NDVI=(NIR-R)/(NIR+R)。

    (1)

    為了保持與地面實測數(shù)據(jù)的時空一致性,選取了2016年6月25日、7月27日、8月28日和9月29日的美國航空航天局(NASA)數(shù)據(jù)中心提供的16 d合成的250 m空間分辨率的MOD13Q1產(chǎn)品用于提取NDVI值。

    根據(jù)設(shè)定的不同放牧梯度,計算3家牧戶的相對放牧強度:

    相對放牧強度=每個牧戶家牦牛頭數(shù)/不同放牧梯度區(qū)域面積(不同半徑的同心圓面積)。

    (2)

    相對放牧強度不表示量化的放牧強度,而是對放牧強度強弱程度的體現(xiàn)。在本研究中,假定每個放牧梯度的區(qū)域是不同半徑的圓形,3個梯度的區(qū)域就是3個以牧場入口為圓心的同心圓,不同梯度上的區(qū)域面積就是圓的面積,且假設(shè)同一牧戶每個梯度上的牦牛頭數(shù)相同。雖然與傳統(tǒng)放牧強度的算法一樣,即家畜頭數(shù)除以草地面積,但因假設(shè)每個梯度上牦牛頭數(shù)相同,發(fā)生變化的只是面積,故將計算結(jié)果定為相對放牧強度,定性表示放牧的強弱程度。

    1.4 數(shù)據(jù)分析

    利用MODIS數(shù)據(jù)重投影工具軟件(MRT)將下載的影像進行NDVI提取,利用ArcMap軟件提取各樣方位置點對應(yīng)的NDVI值,并與實測生物量建立模型。在野外樣點實測過程中,由于GPS定點的誤差以及在稱量、烘干過程中出現(xiàn)失誤(樣品外溢或吸濕),所以在建立回歸模型前需剔除異常值。采用Microsoft Excel和SPSS軟件進行相關(guān)性的統(tǒng)計分析。

    2 結(jié)果

    2.1 相對放牧強度

    根據(jù)公式(2)計算相對放牧強度,如表1所列。

    對比3家牧戶間的相對放牧強度,可以看出,牧戶1的放牧強度比牧戶2和牧戶3都要大,牧戶3次之,牧戶2則最輕。在同一個放牧梯度水平上,除了第三個放牧梯度上牧戶2和牧戶3相對放牧強度相同外,其他兩個梯度上都是牧戶1相對放牧強度最大,牧戶3次之,牧戶2最小。

    表1 3家調(diào)查牧戶的相對放牧強度Table 1 Relative grazing intensity of the three herds under survey

    2.2 植被指數(shù)與草地地上生物量的相關(guān)性分析和模型建立

    建立調(diào)查點植被指數(shù)和地上生物量的線性回歸模型,發(fā)現(xiàn)6-9月植被指數(shù)和地上生物量間的擬合呈相同趨勢,即隨著NDVI值的增長生物量呈上升趨勢。其中,9月的擬合效果最好(圖3)。

    2.3 實測地上生物量和相對放牧強度的相關(guān)分析

    為進一步研究實測地上生物量和相對放牧強度之間的關(guān)系,建立了6-9月實測地上生物量和相對放牧強度間的回歸模型(圖4),發(fā)現(xiàn)二者擬合很好(R2=0.965 1),在某一放牧強度下,地上生物量達到了最大值,而后隨著放牧強度的增大,地上生物量呈下降趨勢。

    圖3 調(diào)查點植被指數(shù)和地上生物量線性回歸散點圖(6-9月)Fig. 3 Linear regression scatter plots of vegetation index at survey points and aboveground biomass (Jun. to Sep.)

    圖4 實測地上生物量和相對放牧強度的冪回歸模型Fig. 4 Power regression model for measured aboveground biomass and relative grazing intensity

    2.4 相對放牧強度與植被指數(shù)的相關(guān)性分析

    分析了地上生物量和植被指數(shù)、相對放牧強度間的相關(guān)關(guān)系后,又建立了實驗月相對放牧強度與對應(yīng)梯度上的NDVI值之間的冪回歸模型(圖5),結(jié)果表明兩者之間不存在相關(guān)性。在剔除了兩個特殊點之后,得到新的模型(圖6),發(fā)現(xiàn)二者相關(guān)性較好(R2=0.631 0)??梢钥闯?,隨著NDVI的變小,相對放牧強度不斷增大,當(dāng)相對放牧強度達到一定程度后,NDVI對相對放牧強度的響應(yīng)不再靈敏。

    圖5 調(diào)查點植被指數(shù)和相對放牧強度的冪回歸模型Fig. 5 Power regression model for vegetation index and relative grazing intensity

    3 討論

    3.1 植被指數(shù)與草地地上生物量之間的相關(guān)性

    在瑪曲高寒草甸基于MODIS-NDVI的線性模型能更好地模擬高寒草甸地上生物量,張艷楠等[18]發(fā)現(xiàn)當(dāng)生物量低于370 g·m-2時,建立的估產(chǎn)模型都是一元線性的,當(dāng)生物量在370~720 g·m-2時,一元線性模型和指數(shù)模型的模擬效果都很好。本研究中生物量最大值為544.12 g·m-2,比較發(fā)現(xiàn)線性模型模擬較好,這與前人研究結(jié)果一致。羅媛等[19]在基于紅光和近紅外反射光譜特征參數(shù)反演草地地上生物量的研究中發(fā)現(xiàn),紅光區(qū)域光譜特征參數(shù)與草地地上生物量之間存在顯著的相關(guān)關(guān)系,試驗表明,地上生物量的增加能夠顯著降低“紅谷”反射率,顯著升高近紅外850 nm處反射率,因此NDVI值隨著生物量的增加而增大,這就可以解釋本研究中NDVI隨著地上生物量的增加而增長的趨勢。比較分析不同月份的地上生物量和NDVI值,發(fā)現(xiàn)9月時生物量有所下降,而且NDVI值也較6月、7月和8月有所下降。一般情況下,牧草會在5月開始返青,地上生物量從生長季開始逐漸增加,到8月達到最大值,而后進入非生長季,牧草開始枯黃,生物量的補償生長滯后于放牧的影響,放牧?xí)绊懫淙郝浣M成且降低禾本科牧草的生物量[20],故到9月時地上生物量有所下降。而且,因為紅光波段是綠色植物進行光合作用吸收的主要波段,9月地上生物量枯黃會導(dǎo)致光合作用減弱,對紅光波段的反射增強,所以NDVI值較于6月、7月和8月會變小。此外,模型的相關(guān)系數(shù)比較低,可能是受研究區(qū)土壤、氣候等條件的影響,另外,選取的MODIS產(chǎn)品數(shù)據(jù)空間分辨率有限,也會影響估產(chǎn)模型的精度。

    圖6 調(diào)查點植被指數(shù)和相對放牧強度的冪回歸模型(剔除特殊點之后)Fig. 6 Power regression model for vegetation index and relative grazing intensity (after exclusion of special points)

    3.2 地上生物量與放牧強度間的相關(guān)性

    仁青吉等[21]研究結(jié)果表明,隨著放牧強度的增加,青藏高原東部瑪曲高寒嵩草草甸群落的高度、蓋度和地上生物量都呈現(xiàn)顯著降低的趨勢。張偉華等[22]研究發(fā)現(xiàn),隨著牧壓強度的增加,錫林河北岸一級臺地退化草地的植物群落高度降低,地表蓋度下降,地上生物量較大幅度減少,特別是優(yōu)質(zhì)牧草生物量迅速減少。這與本研究結(jié)果一致。李建龍等[23]在新疆石河子紫泥泉種羊場進行的不同強度的分區(qū)輪牧試驗也表明,不同輪牧強度處理對草地地上部產(chǎn)量有顯著影響。范國艷[24]研究發(fā)現(xiàn),隨著放牧強度的增加,貝加爾針茅典型草原區(qū)的貝加爾針茅(Stipabaicalensis)草原群落植物物種數(shù)、植被蓋度、地上部現(xiàn)存量均顯著降低,植物群落種類組成結(jié)構(gòu)發(fā)生顯著變化。Joneydi等[25]研究發(fā)現(xiàn),雖然重度放牧降低了干旱灌叢帶地上生物量,但中度放牧對該區(qū)植物生物量沒有影響。在本研究中,中度放牧對該地區(qū)地上生物量是有一定影響的,該差異可能與研究區(qū)域的植被類型、土壤和物候條件等因素有關(guān)。

    3.3 放牧強度與植被指數(shù)間的相關(guān)性

    不同水分條件下的不同植被群落類型對牧草生長速率有一定的影響。模型中提到的兩個特殊點分別是牧戶2和牧戶3在重牧梯度下的NDVI值,一般情況下,NDVI值隨著放牧強度的增大而變小,但此兩點處NDVI值并沒有變小,究其原因發(fā)現(xiàn),牧戶2和牧戶3的牧場有河流經(jīng)過(圖2),在沿河地區(qū),植被層保持著較高的含水量,植物種類成分屬于沼澤草甸,而且植物生產(chǎn)力比周邊草原地區(qū)的要高。因此放牧強度雖然較大,其生物量也較高,故NDVI值也較高。Kawamura等[26]也發(fā)現(xiàn),植物生物量隨著放牧強度的增加而下降,除了錫林沿河地區(qū),本研究的情況與該研究結(jié)果一致。作為家畜喜食,但又不耐放牧踐踏的植物,如鵝絨委陵菜(Potentillaanserina)、羊茅(Festucaovina)、垂穗披堿草(Elymusnutans)、草地早熟禾(Poapratensis)等受到放牧干擾影響較大,不能充分生長和發(fā)育,其蓋度和豐富度有所下降[27],即植物的生成速率小于家畜的采食速率,這可能是NDVI值隨著放牧強度的增強呈下降趨勢的原因之一。另外,從整個群落角度分析,放牧顯著降低了地上生物量,主要原因在于葉片光合面積的減少,有機物積累效率的降低[28]。曹廣民等[29]選取處于不同退化階段的具有典型代表意義的草甸為研究對象,通過對其土壤氮素礦化補給能力、牧草對氮素的需求量的研究,發(fā)現(xiàn)隨著放牧強度的增加,土壤-牧草營養(yǎng)供求失調(diào),土壤養(yǎng)分被地下根系過多吸收被生物固定,成為無效養(yǎng)分[30],這也是生物量降低的主要原因之一,也可能導(dǎo)致研究區(qū)域NDVI值的下降。分析相對放牧強度和NDVI之間的關(guān)系,可以發(fā)現(xiàn)一個相對放牧強度的臨界值,在該值下,NDVI對相對放牧強度的響應(yīng)不再靈敏,這可以作為一個衡量放牧程度強弱的警戒值,對于確定合理放牧強度具有重要的參考價值,而且生物量隨放牧強度的變化也有一個放牧強度閾值存在。草地地上生物量的直接收獲估測方法需要花費相當(dāng)大的時間和勞力,同時要求足夠多的樣本數(shù)目,以保證樣本符合統(tǒng)計學(xué)的要求(樣本數(shù)量少則估測精度有限)[31]。此外,直接收獲法還會破壞調(diào)查地點的植被,具有破壞性,對于一些生態(tài)脆弱區(qū)域的植被恢復(fù)有一定影響,難以大面積開展,再加上家畜的活動是空間動態(tài)變化的過程,如果利用放牧強度和地上生物量的相關(guān)關(guān)系來尋找放牧強度閾值,不僅工作量大,而且有空間限制,對研究結(jié)果有很大的影響。但是,與基于樣方的生物量監(jiān)測的傳統(tǒng)方法相比,遙感估測不僅有更高的時效性,而且具有傳統(tǒng)方法無法比擬的優(yōu)勢,如可以監(jiān)測高寒無人區(qū)草地的變化動態(tài)等[32]。所以,利用遙感植被指數(shù)來獲取家畜的放牧強度閾值具有可行性。

    4 結(jié)論

    本研究利用歸一化植被指數(shù)(NDVI)和地面實測數(shù)據(jù)進行相關(guān)性分析,結(jié)果發(fā)現(xiàn)瑪曲高寒草甸不同月份的地上生物量與對應(yīng)的植被指數(shù)有較好的相關(guān)性,據(jù)此建立了6-9月份調(diào)查點植被指數(shù)和地上生物量的線性回歸模型。發(fā)現(xiàn)6-9月的模擬呈現(xiàn)相同趨勢,即隨著NDVI值的增長,地上生物量呈現(xiàn)增長趨勢。實測地上生物量和相對放牧強度有很好的擬合關(guān)系(R2=0.965 1),隨著放牧強度的增加地上生物量呈下降趨勢。另外,相對放牧強度和植被指數(shù)間相關(guān)性較好,并建立了相對放牧強度估測模型(y=0.000 5x-28.21,R2=0.631 0),結(jié)果表明隨著NDVI值的變小,相對放牧強度逐漸增強,當(dāng)相對放牧強度增強到一定程度時,NDVI對其響應(yīng)不靈敏。量化這一程度上的放牧強度(即閾值),是今后研究中需要做的工作。

    參考文獻References:

    [1] 周華坤,周立,趙新全,劉偉,李英年,古松,周興民.青藏高原高寒草甸生態(tài)系統(tǒng)穩(wěn)定性研究.科學(xué)通報,2006,51(1):63-69.

    Zhou H K,Zhou L,Zhao X Q,Liu W,Li Y N,Gu S,Zhou X M.Study on the stability of alpine meadow ecosystem in Qinghai-Tibet Plateau.Chinese Science Bulletin,2006,51(1):63-69.(in Chinese)

    [2] 周華坤,趙新全,趙亮,李英年,汪詩平,徐世曉,周立.青藏高原高寒草甸生態(tài)系統(tǒng)的恢復(fù)能力.生態(tài)學(xué)雜志,2008,27(5):697-704.

    Zhou H K,Zhao X Q,Zhao L,Li Y N,Wang S P,Xu S X,Zhou L.Restoration capability of alpine meadow ecosystem on Qinghai-Tibetan Plateau.Chinese Journal of Ecology,2008,27(5):697-704.(in Chinese)

    [3] 韓道瑞,曹廣民,郭小偉,張法偉,李以康,林麗,李婧,唐艷鴻,古松.青藏高原高寒草甸生態(tài)系統(tǒng)碳增匯潛力.生態(tài)學(xué)報,2011,31(24):7408-7417.

    Han D R,Cao G M,Guo X W,Zhang F W,Li Y K,Lin L,Li J,Tang Y H,Gu S.The potential of carbon sink in alpine meadow ecosystem on the Qinghai-Tibetan Plateau.Acta Ecologica Sinica,2011,31(24):7408-7417.(in Chinese)

    [4] 寶音賀希格.放牧對草地生態(tài)系統(tǒng)的影響.畜牧與飼料科學(xué),2013,34(7-8):71-74.

    Baoyinhexige.Effect of grazing on grassland ecosystem.Animal Husbandry & Feed Science,2013,34(7-8):71-74.(in Chinese)

    [5] 周華坤,趙新全,周立,劉偉,李英年,唐艷鴻.青藏高原高寒草甸的植被退化與土壤退化特征研究.草業(yè)學(xué)報,2005,14(3):31-40.

    Zhou H K,Zhao X Q,Zhou L,Liu W,Li Y N,Tang Y H.A study on correlations between vegetation degradation and soil degradation in the ‘a(chǎn)lpine meadow’ of the Qinghai-Tibetan Plateau.Acta Pratacultural Sinica,2005,14(3):31-40.(in Chinese)

    [6] 趙雪雁.黃河首曲地區(qū)草地退化的人文因素分析:以甘肅省瑪曲縣為例.資源科學(xué),2007,29(5):50-56.

    Zhao X Y.Human dimensions of grassland degradation:A case study in Maqu County of Gansu Province.Resources Science,2007,29(5):50-56.(in Chinese)

    [7] 杜玉娥,劉寶康,郭正剛.基于MODIS的青藏高原牧草生長季草地生物量動態(tài).草業(yè)科學(xué),2011,28(6):1117-1123.

    Du Y E,Liu B K,Guo Z G.Changes of forage biomass of grasslands during the growing season in the Qinghai-Tibetan Plateau based on MODIS data.Pratacultural Science,2011,28(6):1117-1123.(in Chinese)

    [8] 渠翠平,關(guān)德新,王安志,金昌杰,倪攀,袁鳳輝,張曉靜.基于MODIS數(shù)據(jù)的草地生物量估算模型比較.生態(tài)學(xué)雜志,2008,27(11):2028-2032.

    Qu C P,Guan D X,Wang A Z,Jin C J,Ni P,Yuan F H,Zhang X J.Comparison of grassland biomass estimation models based on MODIS data.Chinese Journal of Ecology,2008,27(11):2028-2032.(in Chinese)

    [9] 方金.基于多源遙感數(shù)據(jù)的甘南牧區(qū)草地生物量遙感監(jiān)測研究.蘭州:蘭州大學(xué)碩士學(xué)位論文,2013.

    Fang J.Monitoring grassland biomass based on multi-source remote sensing data in Gannan pastoral area.Master Thesis.Lanzhou:Lanzhou University,2013.(in Chinese)

    [10] Turker C J,Sellers P J.Satellite remote sensing of primary production.International Journal of Remote Sensing,1986,7(11):1395-1416.

    [11] Moses A C,Andrew K S,Fabio C,van Sip W,Istiak S.Estimation of green grass/herb biomass from airborne hyperspectral imagery using spectral indices and partial least squares regression.International Journal of Applied Earth Observation and Geoinformation,2007,9(4):414-424.

    [12] David S,Emanuel W,Rizzo B,Thomas M S,Woodward I,Fisher H,Timothy G F K,McKeown R,Thomas H P,Rosenbloom N,Ojima D S,Parton W J,David W K,McGuire A D,Melillo J,Pan Y,Haxeltine A,Prentice I C,Sitch S,Hibbard K,Nemani R,Lars L P,Running S,Borchers J,Chaney J,Ronald P N,Bobby H B.Continental scale variability in ecosystem processes:Models,data,and the role of disturbance.Ecological Monographs,1997,67(2):251-271.

    [13] 段敏杰.放牧干擾下藏北紫花針茅高寒草地生物量遙感監(jiān)測.北京:中國農(nóng)業(yè)科學(xué)院碩士學(xué)位論文,2011.

    Duan M J.Remote sensing monitoring ofStipapurpureaalpine grassland biomass under the grazing disturbance in northern Tibet.Master Thesis.Beijing:Chinese Academy of Agricultural Sciences,2011.(in Chinese)

    [14] Paruelo J M,Epstein H E,Lauenroth W K,Burke I C.ANPP estimates from NDVI for the central grassland region of the United States.Ecology,1997,78(3):953-958.

    [15] Bénié G B,Kaboré S S,Go?ta K,Courel M F.Remote sensing-based spatio-temporal modeling to predict biomass in sahelian grazing ecosystem.Ecological Modelling,2005,184(2-4):341-354.

    [16] 韓波.三江源區(qū)高寒草地地上生物量遙感反演模型的建立.淮南:安徽理工大學(xué)碩士學(xué)位論文,2015.

    Han B.Modeling aboveground biomass of alpine grassland in the three-river headwaters region based on remote sensing data.Master Thesis.Huainan:Anhui University of Science & Technology,2015.(in Chinese)

    [17] 周鑫,盛建東,張文太,武紅旗,王新軍.基于MODIS數(shù)據(jù)的伊犁地區(qū)草地地上生物量反演.草地學(xué)報,2015,23(1):27-33.

    Zhou X,Sheng J D,Zhang W T,Wu H Q,Wang X J.Aboveground biomass inversion of grassland in Ili area using MODIS data.Acta Agrestia Sinica,2015,23(1):27-33.(in Chinese)

    [18] 張艷楠,牛建明,張慶,楊艷,董建軍.植被指數(shù)在典型草原生物量遙感估測應(yīng)用中的問題探討.草業(yè)學(xué)報,2012,21(1):229-238.

    Zhang Y N,Niu J M,Zhang Q,Yang Y,Dong J J.A discussion on applications of vegetation index for estimating aboveground biomass of typical steppe.Acta Prataculturae Sinica,2012,21(1):229-238.(in Chinese)

    [19] 羅媛,謝堂民,龍顯靜,馮樹林,陳功.基于紅光和近紅外反射光譜特征參數(shù)反演草地地上生物量.草原與草坪,2015(5):65-69.

    Luo Y,Xie T M,Long X J,Feng S L,Chen G.Pasture aboveground biomass estimating based on spectral reflectance characters of red and near-infrared bands.Grassland and Turf,2015(5):65-69.(in Chinese)

    [20] 李春剛.放牧干擾對高寒草甸群落特征的影響.貴州農(nóng)業(yè)科學(xué),2010,38(5):128-132.

    Li C G.Influence of grazing on community characteristics of alpine meadow.Guizhou Agricultural Sciences,2010,38(5):128-132.(in Chinese)

    [21] 仁青吉,武高林,任國華.放牧強度對青藏高原東部高寒草甸植物群落特征的影響.草業(yè)學(xué)報,2009,18(5):256-261.

    Ren Q J,Wu G L,Ren G H.Effect of grazing intensity on characteristics of alpine meadow communities in the eastern Qinghai-Tibetan Plateau.Acta Prataculturae Sinica,2009,18(5):256-261.(in Chinese)

    [22] 張偉華,關(guān)世英,李躍進.不同牧壓強度對草原土壤水分、養(yǎng)分及其地上生物量的影響.干旱區(qū)資源與環(huán)境,2000,14(4):61-64.

    Zhang W H,Guan S Y,Li Y J.Effect of grazing capacity on water content,nutrient and biomass of steppe soil.Journal of Arid Land Resources & Environment,2000,14(4):61-64.(in Chinese)

    [23] 李建龍,許鵬,孟林,王建華.不同輪牧強度對天山北坡低山帶蒿屬荒漠春秋場土草畜影響研究.草業(yè)學(xué)報,1993,2(2):60-65.

    Li J L,Xu P,Meng L,Wang J H.The comprehensive effects of different rotational grazing intensities (DRGI) on the soil,grassland and sheep production in a spring-autumn pasture of sagebruch desert in the northern slope of tianshan mountain.Acta Prataculturae Sinica,1993,2(2):60-65.(in Chinese)

    [24] 范國艷.放牧對貝加爾針茅草原群落地上生物量和根系分布特征的影響.蘭州:甘肅農(nóng)業(yè)大學(xué)碩士學(xué)位論文,2010.

    Fan G Y.The influence of different grazing intensities on the above-biomass and root distribution character ofStipabaicalensisgrassland.Master Thesis.Lanzhou:Gansu Agricultural University,2010.(in Chinese)

    [25] Joneydi J H,Azarnivand H,Zare C M A,Jafari M.Study of above-ground and below-ground biomass ofArtemisiasieberishrub lands with different grazing intensities in Semnan Province,Iran.Watershed Management Researches,2013,26(2):33-41.

    [26] Kawamura K,Akiyama T,Yokota H,Tsutsumi M,Yasuda T,Watanabe O,Wang S P.Quantifying grazing intensities using geographic information systems and satellite remote sensing in the Xilingol steppe region,Inner Mongolia,China.Agriculture Ecosystems and Environment,2005,107(1):83-93.

    [27] 朱紹宏.牦牛放牧強度對高寒草原植被群落及土壤養(yǎng)分影響的研究.蘭州:甘肅農(nóng)業(yè)大學(xué)碩士學(xué)位論文,2006.

    Zhu S H.Study on the effect of the plant community and soil nutrient of alpine grassland on yak grazing intensity.Master Thesis.Lanzhou:Gansu Agricultural University,2006.(in Chinese)

    [28] 譚紅妍,閆瑞瑞,閆玉春,陳寶瑞,辛?xí)云?不同放牧強度下溫性草甸草原土壤生物性狀及與地上植被的關(guān)系.中國農(nóng)業(yè)科學(xué),2014,47(23):4658-4667.

    Tan H Y,Yan R R,Yan Y C,Chen B R,Xin X P.The relationship between temperate meadow steppe soil’s biological properties and aboveground vegetation under different grazing intensities.Scientia Agricultura Sinica,2014,47(23):4658-4667.(in Chinese)

    [29] 曹廣民,吳琴,李東,胡啟武,李月梅,王溪.土壤-牧草氮素供需狀況變化對高寒草甸植被演替與草地退化的影響.生態(tài)學(xué)雜志,2004,23(6):25-28.

    Cao G M,Wu Q,Li D,Hu Q W,Li Y M,Wang X.Effects of nitrogen supply and demand status of soil and herbage system on vegetation succession and grassland degradation in alpine meadow.Chinese Journal of Ecology,2004,23(6):25-28.(in Chinese)

    [30] 杜巖功,梁東營,曹廣民,王啟蘭,王長庭.放牧強度對嵩草草甸草氈表層及草地營養(yǎng)和水分利用的影響.草業(yè)學(xué)報,2008,17(3):146-150.

    Du Y G,Liang D Y,Cao G M,Wang Q L,Wang C T.Effect on incrassation of mattic epipedon and pasture nutrient-moisture using by grazing in alpineKobresiameadow.Acta Prataculturae Sinica,2008,17(3):146-150.(in Chinese)

    [31] 梁衛(wèi)衛(wèi).草地地上生物量的簡易估測方法研究.楊凌:西北農(nóng)林科技大學(xué)碩士學(xué)位論文,2015.

    Liang W W.The study on simple estimation method of aboveground biomass of grassland.Master Thesis.Yangling:Northwest Agriculture and Forestry University,2015.(in Chinese)

    [32] 楊淑霞,張文娟,馮琦勝,孟寶平,高金龍,梁天剛.基于MODIS逐日地表反射率數(shù)據(jù)的青南地區(qū)草地生長狀況遙感監(jiān)測研究.草業(yè)學(xué)報,2016,25(8):14-26.

    Yang S X,Zhang W J,Feng Q S,Meng B P,Gao J L,Liang T G.Monitoring of grassland herbage accumulation by remote sensing using MODIS daily surface reflectance data in the Qingnan region.Acta Prataculturae Sinica,2016,25(8):14-26.(in Chinese)

    猜你喜歡
    牧戶草甸植被指數(shù)
    高山草甸
    The world's narrowest river
    山地草甸
    生計資本對牧戶牲畜養(yǎng)殖規(guī)模決策的影響
    ——以內(nèi)蒙古錫林郭勒盟為例
    飼料博覽(2021年4期)2021-06-29 05:48:38
    草原牧戶游升級記
    基于草場生存評估模型的阿巴嘎旗北部牧戶草場退化現(xiàn)狀及原因分析
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    AMSR_2微波植被指數(shù)在黃河流域的適用性對比與分析
    河南省冬小麥產(chǎn)量遙感監(jiān)測精度比較研究
    我國牧戶生計脆弱性研究進展及應(yīng)對措施
    免费观看的影片在线观看| 精品人妻一区二区三区麻豆 | av天堂中文字幕网| 少妇人妻一区二区三区视频| 90打野战视频偷拍视频| 亚洲五月婷婷丁香| 亚洲人成网站在线播| 午夜日韩欧美国产| 国内少妇人妻偷人精品xxx网站| 国产精品免费一区二区三区在线| 中国美女看黄片| 欧美黄色淫秽网站| 亚洲成av人片在线播放无| 国内少妇人妻偷人精品xxx网站| 亚洲18禁久久av| 在线观看午夜福利视频| 成人18禁在线播放| 白带黄色成豆腐渣| 国产真人三级小视频在线观看| 香蕉丝袜av| 亚洲人成伊人成综合网2020| 午夜免费男女啪啪视频观看 | 亚洲一区高清亚洲精品| 亚洲国产日韩欧美精品在线观看 | 蜜桃亚洲精品一区二区三区| 国产免费男女视频| 亚洲美女黄片视频| 麻豆国产97在线/欧美| 国产成人影院久久av| 国产精品 欧美亚洲| 欧美不卡视频在线免费观看| 免费av毛片视频| 久久亚洲真实| 久久人妻av系列| 一级黄色大片毛片| 99久久精品热视频| 国产精品久久视频播放| svipshipincom国产片| 老司机在亚洲福利影院| 嫩草影视91久久| 国产一区二区亚洲精品在线观看| av片东京热男人的天堂| 亚洲午夜理论影院| 精品人妻1区二区| 90打野战视频偷拍视频| 在线播放无遮挡| 在线观看66精品国产| 免费观看精品视频网站| 中文亚洲av片在线观看爽| 高清日韩中文字幕在线| 俺也久久电影网| 最近最新中文字幕大全电影3| 国产精品98久久久久久宅男小说| 亚洲激情在线av| 国模一区二区三区四区视频| 亚洲在线自拍视频| 精品国产超薄肉色丝袜足j| 在线免费观看的www视频| www.www免费av| 国产高清三级在线| 国产极品精品免费视频能看的| 久久中文看片网| 嫁个100分男人电影在线观看| 国产在视频线在精品| 又紧又爽又黄一区二区| 午夜免费观看网址| 日韩高清综合在线| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看| 美女黄网站色视频| 18禁裸乳无遮挡免费网站照片| av黄色大香蕉| 亚洲成人久久性| 哪里可以看免费的av片| 亚洲成人久久爱视频| 免费一级毛片在线播放高清视频| 亚洲专区中文字幕在线| 黄色视频,在线免费观看| 国产亚洲欧美98| 国产探花极品一区二区| 村上凉子中文字幕在线| 18禁美女被吸乳视频| 久久草成人影院| www日本在线高清视频| eeuss影院久久| 国产精品 国内视频| 日本在线视频免费播放| 日韩中文字幕欧美一区二区| 欧美日韩中文字幕国产精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 欧美乱色亚洲激情| 国产一区在线观看成人免费| a在线观看视频网站| 两个人看的免费小视频| 1000部很黄的大片| 国产精品日韩av在线免费观看| 国产成人av教育| 亚洲成av人片在线播放无| 国产真实伦视频高清在线观看 | 嫩草影院精品99| 欧美色欧美亚洲另类二区| 又黄又粗又硬又大视频| 欧美bdsm另类| 欧美成狂野欧美在线观看| 脱女人内裤的视频| 波多野结衣高清无吗| 免费观看人在逋| 人妻久久中文字幕网| 黄色日韩在线| 国产一区二区激情短视频| 很黄的视频免费| 国产精品亚洲美女久久久| 亚洲精品成人久久久久久| 天天躁日日操中文字幕| 两个人看的免费小视频| 亚洲专区中文字幕在线| 一个人免费在线观看电影| 亚洲av电影不卡..在线观看| 日本成人三级电影网站| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 久久精品国产综合久久久| 亚洲成人中文字幕在线播放| 99riav亚洲国产免费| 久久久精品大字幕| 99国产极品粉嫩在线观看| 狠狠狠狠99中文字幕| 国产久久久一区二区三区| 国产精品久久久久久久电影 | 国产免费男女视频| 国产精品久久久人人做人人爽| 黄片大片在线免费观看| 日日摸夜夜添夜夜添小说| 精品人妻1区二区| 99国产精品一区二区蜜桃av| 午夜精品一区二区三区免费看| 嫩草影视91久久| 搡老妇女老女人老熟妇| 国产av一区在线观看免费| 人妻夜夜爽99麻豆av| 欧美性猛交黑人性爽| 首页视频小说图片口味搜索| 中文亚洲av片在线观看爽| 女生性感内裤真人,穿戴方法视频| 国产视频内射| 免费av不卡在线播放| 亚洲av免费在线观看| 操出白浆在线播放| 久久精品国产亚洲av涩爱 | 久久久久国内视频| 国产激情欧美一区二区| 午夜影院日韩av| 99久久无色码亚洲精品果冻| 真人做人爱边吃奶动态| 国产成人系列免费观看| 成人无遮挡网站| 可以在线观看毛片的网站| 最好的美女福利视频网| 欧美成人一区二区免费高清观看| 国产亚洲精品综合一区在线观看| 国产伦精品一区二区三区视频9 | 国产一级毛片七仙女欲春2| 欧美日韩中文字幕国产精品一区二区三区| 久久伊人香网站| 国产一区二区三区视频了| 九九热线精品视视频播放| 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| 亚洲av一区综合| 中文字幕人妻熟人妻熟丝袜美 | 精品午夜福利视频在线观看一区| 麻豆久久精品国产亚洲av| 亚洲狠狠婷婷综合久久图片| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久,| 午夜福利在线观看吧| 国产精品野战在线观看| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 色老头精品视频在线观看| 岛国视频午夜一区免费看| 99热6这里只有精品| 欧美日韩瑟瑟在线播放| 男女床上黄色一级片免费看| 成人鲁丝片一二三区免费| 国产精品一区二区三区四区免费观看 | 91在线精品国自产拍蜜月 | 中文字幕熟女人妻在线| svipshipincom国产片| 国内少妇人妻偷人精品xxx网站| 亚洲七黄色美女视频| 国产又黄又爽又无遮挡在线| 啦啦啦免费观看视频1| 毛片女人毛片| 国产毛片a区久久久久| 波多野结衣高清作品| 最新在线观看一区二区三区| 悠悠久久av| 成人午夜高清在线视频| 中文字幕久久专区| 国产真人三级小视频在线观看| 一本一本综合久久| 日韩精品中文字幕看吧| 亚洲精品粉嫩美女一区| 婷婷丁香在线五月| 精品人妻偷拍中文字幕| 免费搜索国产男女视频| 一区二区三区国产精品乱码| 国产成年人精品一区二区| 性欧美人与动物交配| 精品国产亚洲在线| 搡老妇女老女人老熟妇| 窝窝影院91人妻| 午夜a级毛片| 在线观看日韩欧美| 黄色成人免费大全| 国产精品精品国产色婷婷| 精品人妻一区二区三区麻豆 | 别揉我奶头~嗯~啊~动态视频| 亚洲国产精品合色在线| 亚洲精品粉嫩美女一区| 国产av一区在线观看免费| 午夜福利在线在线| 十八禁人妻一区二区| 91九色精品人成在线观看| 香蕉av资源在线| 99久久精品国产亚洲精品| 欧美日韩综合久久久久久 | 午夜福利18| 日本一二三区视频观看| 久久精品国产亚洲av涩爱 | 在线免费观看不下载黄p国产 | 国内精品一区二区在线观看| 亚洲av电影在线进入| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产91精品成人一区二区三区| 欧美色欧美亚洲另类二区| 亚洲电影在线观看av| 天美传媒精品一区二区| 国产精品野战在线观看| 高清日韩中文字幕在线| 婷婷丁香在线五月| 此物有八面人人有两片| 午夜福利视频1000在线观看| 亚洲av五月六月丁香网| 叶爱在线成人免费视频播放| 看黄色毛片网站| 一进一出好大好爽视频| 国产成人欧美在线观看| 一区二区三区高清视频在线| 18禁裸乳无遮挡免费网站照片| 国产伦一二天堂av在线观看| 国产亚洲欧美在线一区二区| 日韩大尺度精品在线看网址| 日韩欧美三级三区| 老司机福利观看| 国产精品综合久久久久久久免费| 国产亚洲精品久久久com| 岛国在线观看网站| 亚洲av电影在线进入| 99热这里只有是精品50| 国产高清激情床上av| e午夜精品久久久久久久| 1000部很黄的大片| 一本久久中文字幕| 网址你懂的国产日韩在线| or卡值多少钱| a级一级毛片免费在线观看| 国产精品av视频在线免费观看| x7x7x7水蜜桃| 精品人妻偷拍中文字幕| 禁无遮挡网站| 国产亚洲精品久久久久久毛片| 制服丝袜大香蕉在线| 国产精品电影一区二区三区| 美女高潮喷水抽搐中文字幕| 精品日产1卡2卡| 又粗又爽又猛毛片免费看| 十八禁人妻一区二区| 欧美激情久久久久久爽电影| 天美传媒精品一区二区| 亚洲午夜理论影院| 国产一区二区激情短视频| 国产 一区 欧美 日韩| 搡老岳熟女国产| 美女 人体艺术 gogo| 欧美日韩黄片免| 99久久无色码亚洲精品果冻| 少妇高潮的动态图| 欧美日韩乱码在线| 欧美日本视频| 亚洲狠狠婷婷综合久久图片| 中文字幕高清在线视频| 国产欧美日韩精品亚洲av| 天堂影院成人在线观看| 草草在线视频免费看| 免费看日本二区| 成人国产一区最新在线观看| 国产成人aa在线观看| 午夜两性在线视频| 欧美成人a在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 最近视频中文字幕2019在线8| 精品一区二区三区av网在线观看| 欧美黄色片欧美黄色片| 天堂影院成人在线观看| 在线观看一区二区三区| 精品日产1卡2卡| 精品久久久久久,| 免费在线观看成人毛片| 午夜福利高清视频| 久久久色成人| 国产一级毛片七仙女欲春2| 色精品久久人妻99蜜桃| 免费无遮挡裸体视频| 少妇裸体淫交视频免费看高清| 小蜜桃在线观看免费完整版高清| 亚洲精品国产精品久久久不卡| 亚洲av成人精品一区久久| 草草在线视频免费看| 人妻久久中文字幕网| 无遮挡黄片免费观看| 非洲黑人性xxxx精品又粗又长| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美日韩高清专用| 999久久久精品免费观看国产| 最后的刺客免费高清国语| 欧美激情久久久久久爽电影| 91久久精品国产一区二区成人 | or卡值多少钱| 国产伦精品一区二区三区四那| 97超视频在线观看视频| 国产综合懂色| 99精品欧美一区二区三区四区| 亚洲天堂国产精品一区在线| 国产伦精品一区二区三区视频9 | 五月伊人婷婷丁香| 特大巨黑吊av在线直播| 在线视频色国产色| 成人鲁丝片一二三区免费| 国产男靠女视频免费网站| 国产美女午夜福利| 日韩免费av在线播放| av女优亚洲男人天堂| 国内精品美女久久久久久| a在线观看视频网站| 99久久精品一区二区三区| 亚洲av美国av| 天堂影院成人在线观看| 亚洲,欧美精品.| 亚洲国产精品久久男人天堂| 国产精品久久久人人做人人爽| 亚洲不卡免费看| avwww免费| 精品99又大又爽又粗少妇毛片 | 午夜福利在线在线| 久久久久性生活片| 男人舔奶头视频| 非洲黑人性xxxx精品又粗又长| 大型黄色视频在线免费观看| 在线看三级毛片| 日韩精品青青久久久久久| 少妇人妻精品综合一区二区 | 久久欧美精品欧美久久欧美| 久久天躁狠狠躁夜夜2o2o| 午夜视频国产福利| 一本综合久久免费| 91在线观看av| 99精品欧美一区二区三区四区| 成人亚洲精品av一区二区| 中文字幕av成人在线电影| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜理论影院| 在线观看av片永久免费下载| 亚洲aⅴ乱码一区二区在线播放| 18+在线观看网站| 亚洲av美国av| 亚洲国产精品久久男人天堂| 97碰自拍视频| 国产亚洲av嫩草精品影院| 国产精品乱码一区二三区的特点| 色在线成人网| 久9热在线精品视频| 一本精品99久久精品77| 国产精品亚洲一级av第二区| 听说在线观看完整版免费高清| 欧美极品一区二区三区四区| av福利片在线观看| 欧美性猛交黑人性爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成人亚洲精品av一区二区| 少妇高潮的动态图| 午夜福利欧美成人| 最近最新免费中文字幕在线| 亚洲国产欧美网| 在线观看一区二区三区| 成人av在线播放网站| 久久精品亚洲精品国产色婷小说| 18禁黄网站禁片免费观看直播| 久久午夜亚洲精品久久| 中文资源天堂在线| 97超视频在线观看视频| tocl精华| 99久久99久久久精品蜜桃| 国产野战对白在线观看| 欧美成人a在线观看| 成人欧美大片| 最新中文字幕久久久久| 在线观看日韩欧美| 少妇熟女aⅴ在线视频| 日韩欧美一区二区三区在线观看| 成年免费大片在线观看| 国产精品1区2区在线观看.| 日日干狠狠操夜夜爽| 国产熟女xx| 色吧在线观看| 中出人妻视频一区二区| h日本视频在线播放| 亚洲激情在线av| 夜夜夜夜夜久久久久| 欧美日韩瑟瑟在线播放| 久久香蕉国产精品| 国产视频内射| 亚洲一区二区三区色噜噜| 欧美激情在线99| 脱女人内裤的视频| 欧美日本亚洲视频在线播放| 亚洲国产精品成人综合色| 在线观看舔阴道视频| 国产伦精品一区二区三区视频9 | 熟妇人妻久久中文字幕3abv| 高清日韩中文字幕在线| 免费在线观看成人毛片| 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 午夜日韩欧美国产| 深夜精品福利| 免费在线观看成人毛片| 亚洲成人免费电影在线观看| 亚洲国产精品成人综合色| 国产伦精品一区二区三区四那| 久久精品国产清高在天天线| 在线免费观看不下载黄p国产 | 欧美丝袜亚洲另类 | 久久精品国产99精品国产亚洲性色| 两个人的视频大全免费| 国产成人影院久久av| av福利片在线观看| 麻豆国产97在线/欧美| e午夜精品久久久久久久| 少妇高潮的动态图| 国产成年人精品一区二区| 熟女人妻精品中文字幕| 国产亚洲精品综合一区在线观看| 国产蜜桃级精品一区二区三区| 免费观看的影片在线观看| 天堂av国产一区二区熟女人妻| 在线视频色国产色| 欧美3d第一页| 国产精品99久久99久久久不卡| 免费看美女性在线毛片视频| 国产熟女xx| 欧美性猛交╳xxx乱大交人| 美女被艹到高潮喷水动态| 琪琪午夜伦伦电影理论片6080| 日日摸夜夜添夜夜添小说| 丰满乱子伦码专区| 中文字幕av在线有码专区| av在线天堂中文字幕| 脱女人内裤的视频| 久久亚洲精品不卡| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| 老司机福利观看| 一区福利在线观看| 精品日产1卡2卡| 一区二区三区国产精品乱码| 亚洲熟妇中文字幕五十中出| 制服人妻中文乱码| 高清日韩中文字幕在线| 国产精品永久免费网站| 久久久久久久久大av| 一本精品99久久精品77| 亚洲,欧美精品.| 久久久精品大字幕| АⅤ资源中文在线天堂| 欧美最黄视频在线播放免费| 天堂√8在线中文| 天堂av国产一区二区熟女人妻| 九色成人免费人妻av| xxx96com| 欧美黄色片欧美黄色片| 日韩大尺度精品在线看网址| 免费高清视频大片| 在线播放国产精品三级| 久久精品国产亚洲av香蕉五月| 亚洲成av人片免费观看| 婷婷精品国产亚洲av| 国产高清有码在线观看视频| 色老头精品视频在线观看| 观看美女的网站| www日本在线高清视频| 啦啦啦免费观看视频1| 亚洲人成网站在线播| 国产精品av视频在线免费观看| 真人一进一出gif抽搐免费| 人人妻人人看人人澡| 欧美性猛交╳xxx乱大交人| 内射极品少妇av片p| 国产一区二区三区在线臀色熟女| 97超级碰碰碰精品色视频在线观看| 最好的美女福利视频网| 国内精品久久久久久久电影| 精品久久久久久,| 性欧美人与动物交配| 日韩 欧美 亚洲 中文字幕| tocl精华| 九九在线视频观看精品| 最近最新免费中文字幕在线| 午夜两性在线视频| 美女cb高潮喷水在线观看| 日本黄色视频三级网站网址| 国产成人影院久久av| 国产高潮美女av| 久99久视频精品免费| 国产老妇女一区| 亚洲最大成人手机在线| 成人国产一区最新在线观看| 国产黄色小视频在线观看| 国产亚洲精品av在线| 女生性感内裤真人,穿戴方法视频| 高清在线国产一区| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 欧美又色又爽又黄视频| 精品免费久久久久久久清纯| 尤物成人国产欧美一区二区三区| 97人妻精品一区二区三区麻豆| 两个人视频免费观看高清| 观看免费一级毛片| 婷婷六月久久综合丁香| av黄色大香蕉| 深夜精品福利| 日韩欧美国产在线观看| 3wmmmm亚洲av在线观看| 最新在线观看一区二区三区| 久久久久久九九精品二区国产| 18禁裸乳无遮挡免费网站照片| www日本在线高清视频| 我要搜黄色片| 午夜免费男女啪啪视频观看 | 一a级毛片在线观看| 成年免费大片在线观看| 大型黄色视频在线免费观看| 日本一二三区视频观看| h日本视频在线播放| 中国美女看黄片| 国产一级毛片七仙女欲春2| x7x7x7水蜜桃| 五月伊人婷婷丁香| 国产高清三级在线| 19禁男女啪啪无遮挡网站| 国产亚洲精品av在线| 国产高清有码在线观看视频| 又黄又爽又免费观看的视频| 日日干狠狠操夜夜爽| 国产乱人伦免费视频| 一个人观看的视频www高清免费观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久视频播放| 国产精品99久久99久久久不卡| 一本久久中文字幕| 国产高潮美女av| 久9热在线精品视频| 每晚都被弄得嗷嗷叫到高潮| 1024手机看黄色片| 99久久久亚洲精品蜜臀av| av女优亚洲男人天堂| 精品电影一区二区在线| 久久久国产成人免费| 悠悠久久av| 熟女人妻精品中文字幕| 久久国产乱子伦精品免费另类| 久久伊人香网站| 亚洲熟妇中文字幕五十中出| 9191精品国产免费久久| 啪啪无遮挡十八禁网站| 亚洲中文字幕一区二区三区有码在线看| 中文字幕精品亚洲无线码一区| 亚洲国产精品sss在线观看| 尤物成人国产欧美一区二区三区| 哪里可以看免费的av片| 特级一级黄色大片| 一边摸一边抽搐一进一小说| 日韩欧美 国产精品| 老司机午夜福利在线观看视频| 国产91精品成人一区二区三区| 亚洲av中文字字幕乱码综合| 狠狠狠狠99中文字幕| 精品国产超薄肉色丝袜足j| 老司机午夜十八禁免费视频| 亚洲美女视频黄频| 可以在线观看毛片的网站| 欧美激情在线99| 日韩欧美 国产精品| 桃色一区二区三区在线观看| 一级a爱片免费观看的视频| 国产在视频线在精品| 国产一区二区在线av高清观看| 亚洲无线观看免费| 99久久无色码亚洲精品果冻| 午夜福利高清视频| 最近最新中文字幕大全电影3| 久久6这里有精品| 国产av在哪里看|