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

    基于MODIS NPP數(shù)據(jù)的青海湖流域產(chǎn)草量與載畜量估算研究

    2019-09-21 02:40:12王琪吳成永陳克龍張肖張樂樂丁俊霞
    生態(tài)科學 2019年4期
    關鍵詞:載畜量產(chǎn)草量溫性

    王琪, 吳成永, 陳克龍, 張肖, 張樂樂, 丁俊霞

    基于MODIS NPP數(shù)據(jù)的青海湖流域產(chǎn)草量與載畜量估算研究

    王琪1,2, 吳成永1,2, 陳克龍2,*, 張肖3, 張樂樂1,2, 丁俊霞1,2

    1. 青海師范大學地理科學學院, 青海, 西寧 810008 2. 青海省自然地理與環(huán)境過程重點實驗室, 青海, 西寧 810008 3. 四川師范大學地理與資源科學學院, 四川, 成都 610101

    快速評價區(qū)域草地生產(chǎn)力狀況, 是制定牧區(qū)草畜平衡策略, 保障草地生態(tài)系統(tǒng)健康可持續(xù)發(fā)展的前提基礎?;贛ODIS NPP產(chǎn)品數(shù)據(jù), 對青海湖流域草地生產(chǎn)力進行了估算和評價。結果表明: (1)牛羊可食性最高的高寒草甸、高寒草原和溫性草原, 其單位面積產(chǎn)草量估算值的均方根誤差為26.15 g·m-2, 表明該方法可以快速估算區(qū)域尺度的產(chǎn)草量。(2)全年干草產(chǎn)量為145.42萬t, 其中能被牛羊直接采食利用的牧草共計59.18萬t, 理論載畜量為81.07萬羊單位; (3)影響單位面積干草產(chǎn)量的主要氣候因子, 海拔3500 m以下地區(qū), 是活動積溫且與單位面積干草產(chǎn)量呈負相關關系; 3500 m以上地區(qū), 是溫度且與單位面積產(chǎn)草量呈正相關關系。研究結果可為高寒地區(qū)合理規(guī)劃牧業(yè)生產(chǎn)活動、生態(tài)補償?shù)忍峁┛茖W依據(jù)和決策支持, 也可為全國其他牧區(qū)的草原產(chǎn)草量的估算提供參考。

    產(chǎn)草量; 草地生產(chǎn)力; 植被凈初級生產(chǎn)力; 青海湖流域

    0 引言

    草地是畜牧業(yè)發(fā)展的最基本生產(chǎn)資料和基地。我國草地面積廣闊, 約占全國陸地總面積的41%。青藏高原廣泛分布著高寒草甸、高寒草原、溫性草原等草地生態(tài)系統(tǒng)[1], 總面積約為1.28×108km2[2], 是中國主要的畜牧業(yè)基地之一。青海湖流域是青藏高原東北部重要的牧區(qū), 也是青海省畜牧業(yè)生產(chǎn)的主要基地[3-4], 流域內集中了青海省 24% 的牧業(yè)人口和 29% 的草食牲畜[5]。自20世紀80年代以來, 在人類不合理活動和全球氣候變化的雙重影響下, 草地退化十分嚴重, 草地生態(tài)系統(tǒng)嚴重惡化。1977—2000 年, 流域內草地共減少 206.68 km2(其中, 大部分草地轉變?yōu)楦睾蜕车豙6]); 2000—2010 年, 呈現(xiàn)出整體性草地極輕微退化趨勢[7]。草地退化致使草地生產(chǎn)水平持續(xù)下降, 草畜供求的時空利用格局出現(xiàn)失衡[8-9], 直接影響到區(qū)域生態(tài)安全與畜牧業(yè)經(jīng)濟的可持續(xù)發(fā)展。快速準確地估算流域內的產(chǎn)草量和載畜量, 是解決上述科學與現(xiàn)實問題的首要前提和關鍵。

    草地載畜量, 是指在一定時期和一定草地面積上, 在適牧條件下, 并保證家畜正常發(fā)育和生長的狀態(tài)下, 能飼養(yǎng)放牧家畜的頭數(shù)[10]。載畜量過高則導致牧草利用過度, 草地狀況惡化, 家畜營養(yǎng)匱乏, 降低單個家畜增重。相反, 載畜量過低會造成牧草浪費, 牧草利用率降低, 減少單位草地面積上家畜的總增重[11]。因而, 放牧利用不應超過能被科學計算的載畜量域值[12]。估算理論載畜量, 能科學確定草地載畜量的臨界值, 這一結果不僅能作為指導草地管理的重要依據(jù), 而且也是科學地確定放牧強度, 避免超載過牧, 維持草地持續(xù)生產(chǎn)的有效手段。測算草地載畜量的基礎是估算草地產(chǎn)草量。

    草地產(chǎn)草量的測量估算主要有三種方法[13]: 一是直接收獲法(刈割法), 測量值是齊地面刈割所獲得的產(chǎn)草量, 通常在草地產(chǎn)量最高的時候進行測定[14]。該方法的優(yōu)點是測量準確, 缺點是耗時費力, 受人為因素的干擾大, 適合小區(qū)域草地產(chǎn)量測定, 在大尺度范圍內估產(chǎn)精度不高且它難于估算草食類動物所消耗的生物量或能量以及自養(yǎng)有機體自身的代謝生長和發(fā)育所用的能量[15]; 二是產(chǎn)量模擬模型, 通常需要輸入連續(xù)詳細的氣候、土壤等數(shù)據(jù)參數(shù)。優(yōu)點是估算精度較高, 缺點是數(shù)據(jù)參數(shù)獲取困難, 難以應用于大尺度區(qū)域; 三是遙感模型測量法, 該方法是目前區(qū)域產(chǎn)草量估算的常用方法之一, 借助于植被NDVI或者NPP(Net Primary Productivity)兩種方法來估算?;谥脖籒DVI的方法, 通過回歸模型進行估算, 適合于大區(qū)域。但不同區(qū)域, 其回歸估算模型也不同, 因此模型移植性差。 NPP估算法, 以植被凈初級生產(chǎn)力為基礎計算產(chǎn)草量, 優(yōu)點是宏觀、快速、省力、省時、受外在因素干擾少, 可為草原管理與決策提供及時的草原信息, 如光能利用率模型CASA模型(目前應用較廣泛, 輸出結果為月NPP)。如樸世龍等[16]利用CASA模型對我國植被凈第一性生產(chǎn)力進行了估算; 李剛等[17]利用改進后的CASA 模型估算了青藏高原草地載畜量。以上兩位學者在估算NPP時, 氣象參數(shù)均在插值基礎上獲得。秦泗國等[18]利用綜合模型估算了那曲高原高寒草地生態(tài)系統(tǒng)NPP并在此基礎上估算載畜量, 綜合模型估算NPP是以與植被光合作用密切相關的實際蒸散為基礎的, 關鍵參數(shù)的獲取(如PER)用氣象數(shù)據(jù)反演得到。喬凱等[19]在定量估算青海湖流域2001—2011年草地凈初級生產(chǎn)力時, 氣象數(shù)據(jù)用域內和周邊氣象站點數(shù)據(jù)內插得到, 500 m DEM數(shù)據(jù)用90 m重采樣方法獲得。鄭中等[20]在利用CASA模型估算青海湖流域草地 NPP時, 面對青海湖流域內氣象站稀少, 氣象觀測數(shù)據(jù)相當缺乏的情況下, 用WRF區(qū)域氣候模式模擬了近地表氣象數(shù)據(jù), 但模擬精度只在80%左右。因此, CASA模型模擬NPP, 雖然得到了廣泛應用, 但其缺點是模型需輸入氣溫、降水、太陽輻射以及植被指數(shù)和植被類型等眾多參數(shù)。部分參數(shù)難以獲取, 部分參數(shù)(如氣象參數(shù))需用空間插值技術得到, 而氣象站點的數(shù)目和空間分布等都會影響到插值結果, 精度難以保障。

    青海湖流域境內國家級氣象站點少, 且其地形多為山地, 地勢起伏大, 插值獲取的參數(shù)誤差大, 致使模型模擬得到的NPP精度難以保證。所以獲取精度相對較高的NPP數(shù)據(jù)是估算產(chǎn)草量與載畜量的首要前提。綜上所述, 本文選取分辨率較高(500 m)和準確性已通過驗證的MOD17A3H產(chǎn)品, 以青海湖流域為研究區(qū)域, 選取干草產(chǎn)量與理論載畜量2個指標, 進行草地生產(chǎn)力評價。其一是探討其快速估算產(chǎn)草量的可行性與準確性, 其二是通過數(shù)學方法統(tǒng)計出流域內不同草地類型及行政區(qū)域的產(chǎn)草量, 分析探討草地產(chǎn)草量的空間格局及其影響因素, 旨在為高寒地區(qū)牧業(yè)生產(chǎn)的合理規(guī)劃和畜牧業(yè)可持續(xù)發(fā)展提供科學依據(jù)。

    1 研究區(qū)概況

    青海湖流域地處青藏高原東北部, 總面積為29661 km2, 介于 36°15' — 38°20' N 與97°50' — 101°20'E之間, 地形復雜, 地勢起伏大, 西北高、東南低。地處黃土高原、西北干旱區(qū)和青藏高寒區(qū)的過渡地帶, 深居內陸, 海拔較高, 常年受西風帶、東亞季風和青藏高原季風的交匯控制, 屬典型的高寒干旱大陸性氣候, 干旱、寒冷、多風。多年平均氣溫-1.4—1.7 ℃, 年平均降水量約為340 mm, 蒸發(fā)量為800—1000 mm[21]。植被類型有草原植被(溫性草原和高寒草原)、灌叢植被(溫性河谷灌叢與高寒灌叢)、草甸(高寒草甸、鹽生草甸)和高山流石植被以及荒漠植被(高寒荒漠類和溫性荒漠類)等。其中高寒草原和高寒草甸是青海湖流域的主要草地類型, 兩者占流域內草地總面積的82.64%。2010年青海湖流域總人口為 13.23萬人, 平均人口密度4.47人·km-2[22], 農(nóng)業(yè)活動以畜牧業(yè)生產(chǎn)為主[23]。

    2 數(shù)據(jù)與方法

    2.1 數(shù)據(jù)來源與處理

    (1)遙感數(shù)據(jù)。采用2015年的MODIS數(shù)據(jù)的MOD17A3H產(chǎn)品, 來源于美國國家航空航天局(NASA) (https://ladsweb.modaps.eosdis.nasa.gov/search/)。該產(chǎn)品提供了500 m空間分辨率的NPP年數(shù)據(jù)[24], 已通過第三方驗證, 其準確性得到了評估[25]。

    (2)氣象數(shù)據(jù)。2015年青海省各氣象站點的年平均氣溫、年降水量、平均濕度等, 來自中國氣象科學數(shù)據(jù)共享服務平臺(http://data.cma.cn/)。通過整理、計算和空間插值生成相應的柵格數(shù)據(jù)。

    (3)植被數(shù)據(jù)。運用中國1: 100萬植被類型圖(2000年版)—青藏部分, 根據(jù)中國植物群落的分類方法, 分類得到青海湖流域的草地類型分布圖(圖1)。

    上述三類數(shù)據(jù), 在ArcGIS 10.2軟件平臺中轉換為Albers等積投影。

    2.2 研究方法

    2.2.1 干草產(chǎn)量計算

    干草產(chǎn)量是指某地域某一段時間內生產(chǎn)的牧草干物質總產(chǎn)量。單位面積的干草產(chǎn)量可以利用草地NPP計算得到[26], 計算公式如下:

    Figure 1 Distribution of grassland types in the Qinghai Lake Basin

    式中,B為年總單位面積干草產(chǎn)量,/(m-2·a-1);為草地年總NPP,/(m-2·a-1) ;S為草地生物量換算成NPP 的轉換系數(shù), g·g C-1, 數(shù)值為 0.45[27];S為地下與地上部分生物量比例系數(shù), 高寒草甸為8.99, 高寒草原類為3.81[28], 溫性草原類為4.25, 溫性荒漠草原類、高寒荒漠和溫性荒漠類為7.89, 暖性灌草叢為4.42, 沼澤類為15.68[29]。

    2.2.2 產(chǎn)草貢獻率計算

    產(chǎn)草貢獻率是指某草地類型的產(chǎn)草量占所有草地類型產(chǎn)草量的比重。產(chǎn)草貢獻率計算公式如下:

    其中,C代表第類型草地的產(chǎn)草貢獻率;M為第類型草地的平均產(chǎn)量;A為第類型草地的面積。

    2.2.3 理論載畜量計算

    理論載畜量是指在適度放牧的條件下, 單位草地面積上所能放牧的牲畜頭數(shù)。它是衡量草原生產(chǎn)能力的主要指標之一, 可根據(jù)牧草的生長情況(包括牧草的品質、產(chǎn)量)來評定, 取決于草地初級生產(chǎn)力[30]。計算公式如下:

    式中,Z為草地理論載畜量(羊單位);H為全年干草產(chǎn)量(t·km-2);G為不同草地類型的牧草利用率[31], 草原類為50%, 草甸類為60%, 灌(草)叢類和沼澤類為55%, 荒漠類為40%。C為草地可利用率,E為草地可食牧草比率; 在青海湖流域,C按92%計算,E按80%計算[32];L為一個羊單位家畜的日食量, 我國采用的載畜量為綿羊單位[33], 1個羊單位是活重為40 kg的綿羊及其哺乳羔羊, 每天飼草(青草)需求量為5—7.5 kg, 本研究取6 kg, 據(jù)中國農(nóng)業(yè)行業(yè)提出的天然草地載畜量計算標準[34], 牧草干重與鮮重比例按1:3計算, 換算成干草即2 kg /羊單位/;為草地放牧時間, 按全年365計算。

    3 結果與分析

    3.1 產(chǎn)草量估算結果及精度分析

    首先利用草地類型的矢量圖, 用ArcGIS 10.2軟件中的裁剪工具, 裁出NPP柵格圖, 然后運用公式(1)進行計算, 得到2015年青海湖流域單位面積干草產(chǎn)量分布圖(圖2)。最后利用單位面積干草產(chǎn)量乘以像元面積得到流域內的全年干草產(chǎn)量圖(圖3)。

    青海湖流域地形復雜, 區(qū)域小氣候、土壤和草地類型差異大, 加之人類活動影響程度不一, 因而流域內各草地類型的單位面積干草產(chǎn)量空間差異性明顯。由圖(2)可知, 環(huán)湖地區(qū)是單位面積干草產(chǎn)量的高值區(qū), 呈現(xiàn)出從東南環(huán)湖地區(qū)向西北山地遞減趨勢。環(huán)湖地區(qū)海拔相對較低, 地形平坦, 氣溫較高, 水熱組合較好, 植被以溫性草原, 高寒草原為主, 且為國家級自然保護區(qū)的核心區(qū), 因而整體上產(chǎn)草量最高。西北部, 多山地, 海拔高, 氣溫低, 植被以低矮的高寒草甸為主, 單位面積產(chǎn)草量較低。

    圖2 青海湖流域草地單位面積干草產(chǎn)量空間分布圖

    Figure 2 Hay yield per unit area of grassland in Qinghai Lake Basin

    圖 3 青海湖流域草地年干草產(chǎn)量空間分布圖

    Figure 3 Annual hay yield of grassland in the Qinghai Lake Basin

    采用均方根誤差(式4)檢驗產(chǎn)草量的估算精度。越小, 模型精度越高、預測能力越強。將青海湖流域牛羊可食性最高的高寒草甸、高寒草原和溫性草原的產(chǎn)草量估算值與文獻中產(chǎn)草量進行比較, 得出兩者間的為26.15 g·m-2(表1), 表明模型具有較好的估算精度。

    3.2 干草產(chǎn)量影響因素分析

    氣候影響牧草的生長發(fā)育, 進而影響產(chǎn)草量。青海湖流域屬高海拔地區(qū), 地勢起伏大, 氣候隨海拔梯度的變化明顯。海拔3500 m等高線是青海湖流域湖濱平原區(qū)、河谷區(qū)與高山區(qū)的大致界線。故本文將3500 m等高線作為海拔節(jié)點, 來探討此海拔節(jié)點上下的單位面積干草產(chǎn)量的水熱影響模式。為了定量厘定出氣候因子對產(chǎn)草量的影響, 利用相關系數(shù)r(式5), 對不同海拔梯度內的單位面積干草產(chǎn)量和各氣候因子(氣溫、降水量、活動積溫(≥10 ℃)和濕度)進行相關性分析。

    結果表明(表2), 海拔3500 m以下區(qū)域, 氣象因子中, 對單位面積干草產(chǎn)量影響最大的為活動積溫(≥10℃), 其次為氣溫、降水, 最小的為濕度。其中,年降水量與單位面積干草產(chǎn)量為正相關關系, 而活動積溫、氣溫與干草產(chǎn)量為負相關關系, 即活動積溫、氣溫越高, 干草產(chǎn)量越低。原因可能是: 其一、在氣候變暖背景下, 流域內年均氣溫呈現(xiàn)增加趨勢[20],活動積溫超過了適宜植被生長的溫度范圍, 從而不利于高寒草原植被的生長。其二、流域內年均氣溫的升高導致蒸發(fā)量隨之增大, 土壤水分損失加快, 抑制了植被的生長。在此情況下, 牧草的需水脅迫加重, 從而降水量越大, 單位面積干草產(chǎn)量越大, 這與龍慧靈等[37]研究結果一致。3500 m以上區(qū)域, 除了濕度與單位面積干草產(chǎn)量為負相關以外, 其余均為正相關關系。其中, 對單位面積干草產(chǎn)量影響最大的氣候因子為氣溫, 這與周偉等[38]和陸晴等[39]人的研究結果一致。原因可能是氣溫為高寒草甸等植被生長發(fā)育的最主要限制性因子, 氣溫越高, 植被生長周期越長, 單位面積產(chǎn)草量就越高。另外, 高寒草甸的植被根系發(fā)達, 能夠利用深層土壤水, 從而受降水的影響相對較弱[41]。

    表1 不同草地類型產(chǎn)草量估算結果精度分析

    表2 青海湖流域不同海拔單位面積產(chǎn)量與氣候因子的相關性

    注:*在 0.01 水平(雙側)上顯著相關

    3.3 產(chǎn)草量與理論載畜量分析

    3.3.1 不同草地類型產(chǎn)草量與理論載畜量分析

    對產(chǎn)草量估算結果進行數(shù)學統(tǒng)計(表3), 得出流域內全年干草產(chǎn)量為145.42萬t, 其中能被牛羊直接采食利用的牧草59.18萬 t, 理論載畜量為81.07萬羊單位。高寒草甸、高寒草原、亞高山灌叢、溫性草原、溫性灌叢草原和溫性荒漠草原六大類是主要產(chǎn)草的草地類型, 其年干草產(chǎn)量: 高寒草甸>高寒草原>亞高山灌叢>溫性草原>溫性灌草叢>溫性荒漠草原。高寒草甸的草地面積、理論產(chǎn)草量、可食牧草量和采食總量最大, 分別為17127.2 km2、712843.72 t、314791.79 t、431222羊單位; 其次為高寒草原, 分別為3258.64 km2、391418.61 t、144042.05 t、197318羊單位。高寒草甸與高寒草原面積占流域內草地總面積的82.64%, 產(chǎn)草量占整個流域產(chǎn)草量的75.93%, 載畜能力也最大, 占所有草地載畜量的77.53%。

    表3 青海湖流域各草地類型草地生產(chǎn)力統(tǒng)計表

    據(jù)圖4, 青海湖流域所有草地類型中, 高山草甸的產(chǎn)草貢獻率最高, 占整個流域內產(chǎn)草量的49.01%; 其次為高寒草原, 產(chǎn)草貢獻占整個流域的26.91%; 產(chǎn)草貢獻率最低的是沼澤化草甸, 僅占0.18%。

    3.3.2 行政區(qū)劃產(chǎn)草量與載畜狀況分析

    統(tǒng)計得到青海湖流域內各行政區(qū)劃的理論載畜量, 天峻縣、剛察縣、共和縣和海晏縣理論載畜量分別為307092羊單位、251267羊單位、189142羊單位和63191羊單位。根據(jù)2015統(tǒng)計年鑒和縣農(nóng)牧部門公布的4縣實際載畜量數(shù)據(jù)[41]計算4縣在流域內的實際載畜量

    (公式為: 某縣所占青海湖流域面積÷該縣面積×該縣實際載畜量):

    圖4 青海湖流域不同草地類型的產(chǎn)草貢獻率

    Figure 4 Grassland production contribution rate of different grassland types in the Qinghai Lake Basin

    天峻縣實際載畜量為497261(頭、匹), 剛察縣789743(頭、匹), 共和縣656065(頭、匹), 海晏縣191559(頭、匹)。進一步計算青海湖流域超載率, 達163.31%, 其中共和縣超載最為嚴重, 達246.86%, 次之為剛察縣與海晏縣, 超載率分別為214.30%和203.14%, 天峻縣為61.93%。

    為了制定切實可行的牧業(yè)規(guī)劃和生態(tài)補償政策的落地實施等, 本文對青海湖流域內16個主要鄉(xiāng)鎮(zhèn)的產(chǎn)草量和理論載畜量進行統(tǒng)計(表4), 結果顯示: 產(chǎn)草量和理論載畜量較高的為泉吉鄉(xiāng)、石乃亥鄉(xiāng)、哈爾蓋鎮(zhèn)和吉爾孟鄉(xiāng), 年干產(chǎn)草量共計46.04萬t, 可采食牧草總量達182552.8 t, 可放牧羊單位數(shù)達250072只, 占全流域理論載畜量的30.84%。

    4 討論

    模擬估算結果的準確性驗證是應用NPP估算區(qū)域產(chǎn)草量的前提, 這可以通過田間調查獲得的實測生物量或者已有研究結果相比較。青海湖流域牧草可食性較大的三種草地類型的產(chǎn)草量, 其單位面積內估算值的均方根誤差為26.15 g·m-2(其中, 高寒草原的誤差較大, 原因是驗證數(shù)據(jù)為2010年數(shù)據(jù), 而本文所采用的是2015年數(shù)據(jù), 根據(jù)駱成鳳等[9]人的研究, 湖濱平原和河谷地區(qū)目前仍處于過牧狀態(tài), 草地生產(chǎn)力持續(xù)下降, 因此這年際間誤差可以忽略不計), 在可接受范圍內。由于NPP的影響因素較多, 如天氣狀況, 土壤水分的多寡, 植被類型及其所處地形部位等, 導致不同草地類型的NPP模擬精度也不一致, 因而產(chǎn)草量的估算精度可能不同, 甚至有的草地類型估算精度較低, 但只要研究區(qū)中面積比重大的草地類型產(chǎn)草量估算精度高, 則估算方法就可以推廣應用。

    表4 青海湖流域各鄉(xiāng)鎮(zhèn)草地生產(chǎn)力統(tǒng)計表

    依據(jù)統(tǒng)計年鑒和縣農(nóng)牧部門公布的實際載畜量數(shù)據(jù), 然后結合理論載畜量, 計算了超載率。青海湖流域4縣, 總超載率達163.31%。整體上, 青海湖流域內各行政區(qū)劃的實際載畜量高于理論載畜量且已超過草地所能承受的最大載畜能力。需要說明的是, 理論載畜量的計算基礎是所有植被的NPP, 但并非所有的植被能被牛羊等牲畜采食, 因此, 實際上, 各地區(qū)的超載率應該更高, 這也同時提出了未來研究的一個方向, 即如何識別與剔除牲畜非食用植被, 進而提高估測載畜量的精度。草地畜牧業(yè)要得到健康、穩(wěn)定的發(fā)展, 減輕草地載畜壓力是亟待解決的問題, 農(nóng)牧部門需進行調控和管理(如加大對牲畜的補飼力度以及減少牲畜的存欄數(shù)量), 使牧民遵循草原的可更新性和季節(jié)生長規(guī)律, 通過確定科學的放牧時間、放牧強度、放牧方式及牧場輪換來穩(wěn)定流域內草地自然健康生長。

    5 結論

    基于MODIS NPP遙感產(chǎn)品數(shù)據(jù), 結合植被類型圖, 本文估算了青海湖流域草地生產(chǎn)力, 并得出以下結論:

    (1)青海湖流域草地平均單位面積干產(chǎn)草量為74.38 g·m-2, 呈現(xiàn)出自環(huán)湖區(qū)域向高山區(qū)遞減趨勢, 影響其分布的主要氣象因子是活動積溫、氣溫與降水。3500 m以下區(qū)域, 活動積溫是主要影響因子; 3500 m以上區(qū)域, 溫度是主要影響因子。降水對3500 m以上單位面積干草產(chǎn)量影響較小原因是高寒草甸的根系較發(fā)達, 能吸收深層土壤水, 從而受降水限制小。

    (2)高寒草甸、高寒草原、亞高山灌叢、溫性草原、溫性灌叢草原和溫性荒漠草原六大類草地是青海湖流域的主要產(chǎn)草類型, 年總干草產(chǎn)量為145.42萬t。其年干草產(chǎn)量: 高寒草甸>高寒草原>亞高山灌叢>溫性草原>溫性灌草叢>溫性荒漠草原。

    (3)青海湖流域理論載畜量為81.07萬羊單位, 實際載畜量為213.46頭。牲畜超載十分嚴重, 超載率達163.31%。其中共和縣超載最為嚴重, 達246.86%;次之為剛察縣與海晏縣, 超載率分別為214.30%和203.14%; 天峻縣為61.93%。

    本文探討了基于MODIS NPP遙感產(chǎn)品數(shù)據(jù)的草地生產(chǎn)力計算方法的可行性, 后續(xù)研究將結合實地采樣, 進一步產(chǎn)提高草地產(chǎn)草量估算的準確性。在此基礎上, 結合土地第二次調查資料和鄉(xiāng)鎮(zhèn)村牛羊數(shù)量, 分析地塊尺度的產(chǎn)草量, 確定合理載畜量, 以便采取合理的放牧策略, 實現(xiàn)草地畜牧業(yè)生產(chǎn)的生態(tài)和經(jīng)濟效益可持續(xù)發(fā)展。

    [1] 袁子茹, 任靈, 陳建綱, 等. 祁連山不同草地類型土壤有機質與全氮分布的關系[J]. 草原與草坪, 2016, 36(3): 12– 16.

    [2] 李以康, 冉飛, 包蘇科, 等. 高寒草甸矮嵩草對氮肥的生理響應[J]. 草業(yè)學報, 2010, 19(3): 240–244.

    [3] 魏克家, 胡瑞寧, 馬建華, 等. 青海湖環(huán)湖區(qū)草地生態(tài)環(huán)境問題與防治對策研討[J]. 四川草原, 1997, 16(2): 2–6.

    [4] 張旭萍, 郭連云, 田輝春. 環(huán)青海湖盆地氣候變化對草地生態(tài)環(huán)境的影響[J]. 草原與草坪, 2008, 26(2): 64–69.

    [5] 周薇薇, 殷青軍. 青海湖流域草地產(chǎn)草量與MODIS植被指數(shù)的相關性分析[J]. 鹽湖研究, 2012, 20(4): 8–12.

    [6] 張明, 崔軍, 曹學章. 青海湖流域草地退化時空分布特征[J]. 生態(tài)與農(nóng)村環(huán)境學報, 2017, 33(5): 426–432.

    [7] 駱成鳳, 許長軍, 游浩妍, 等. 2000—2010年青海湖流域草地退化狀況時空分析[J]. 生態(tài)學報, 2013, 33(14): 4450–4459.

    [8] 王輝, 任繼周, 袁宏波. 黃河源區(qū)天然草地沙化機理分析研究[J]. 草業(yè)學報, 2006, 15(6): 19–25.

    [9] 衛(wèi)亞星, 陳全功, 王一謀, 等. 利用TM資料調查土地利用狀況動態(tài)變化——以瑪曲縣為例[J]. 草業(yè)科學, 2002, 17(3): 6–8.

    [10] 甘肅農(nóng)業(yè)大學. 草原調查與規(guī)劃[M]. 北京: 農(nóng)業(yè)出版社, 1985.

    [11] 朝克圖, 那亞, 那松烏力吉, 金花. 載畜量的研究進展概述[J]. 內蒙古草業(yè), 2010, 22(1): 12–14.

    [12] 楊理, 侯向陽. 草畜平衡管理與草地資源可持續(xù)利用[J]. 中國農(nóng)業(yè)經(jīng)濟評論, 2005, 4(3): 453–461.

    [13] 安卯柱, 高娃, 朝魯. 內蒙古第四次草地資源調查草地生產(chǎn)力測定及計算方法簡介[J]. 內蒙古草業(yè), 2002, 14(4): 20–21.

    [14] 任繼周.草業(yè)科學研究方法[M]. 北京: 中國農(nóng)業(yè)出版社, 1998.

    [15] 李高飛, 任海, 李巖, 等. 植被凈第一性生產(chǎn)力研究回顧與發(fā)展趨勢[J]. 生態(tài)科學, 2003, 22(04): 360–365.

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

    [17] 李剛, 孫煒琳, 張華, 等. 基于秸稈補飼的青藏高原草地載畜量平衡遙感監(jiān)測[J]. 農(nóng)業(yè)工程學報, 2014, 30(17): 200–211.

    [18] 秦泗國, 鐘國輝, 王景升. 那曲草地氣候格局對草地NPP的影響及載畜量研究[J]. 干旱區(qū)資源與環(huán)境, 2010, 24(07): 159–164.

    [19] 喬凱, 郭偉. 青海湖流域植被的凈初級生產(chǎn)力估算[J]. 水土保持通報, 2016, 36(06): 204–209.

    [20] 鄭中, 祁元, 潘小多, 等. 基于WRF模式數(shù)據(jù)和CASA模型的青海湖流域草地NPP估算研究[J]. 冰川凍土, 2013, 35(02): 465–474.

    [21] 朱延龍, 韓昆, 王芳. 青海湖流域氣候變化特點及水文生態(tài)響應[J]. 中國水利水電科學研究院學報, 2012, 9 (4): 260–266.

    [22] 王琪, 陳克龍, 吳成永. 基于GIS的青海湖流域公路交通通達性分析[J]. 云南地理環(huán)境研究, 2017, 29(2): 29– 32+56+79.

    [23] 蔡國英, 尹小娟, 趙繼榮. 青海湖流域人類福祉認知及綜合評價[J]. 冰川凍土, 2014, 34(2): 469–478.

    [24] RUNNING S W, MU Q, ZHAO M. MOD17A3H MODIS/Terra Net Primary Production Yearly L4 Global 500m SIN Grid V006[DB/OL]. NASA EOSDIS Land Processes DAAC, 2015. https: //doi. org/10. 5067/MODIS/ MOD17A3H. 006.

    [25] NASA. MOD17A凈初級生產(chǎn)力產(chǎn)品的第3階段驗證準確性聲明[DB/OL]. (2014-4-1)[ 2015-11-13] . https: // landval. gsfc. nasa. gov/ProductStatus. php?ProductID= MOD17.

    [26] 陳世榮, 王世新, 周藝. 基于遙感的中國草地生產(chǎn)力初步計算[J]. 農(nóng)業(yè)工程學報, 2008, 22(1): 208–212.

    [27] 方精云, 劉國華, 徐嵩齡. 中國陸地生態(tài)系統(tǒng)的碳循環(huán)及其全球意義[A]. 王庚辰, 溫玉璞. 溫室氣體濃度和排放監(jiān)測及相關過程[D]. 北京: 中國環(huán)境科學出版社, 1996: 129–139.

    [28] ZENG Chaoxu, WU Jianshuang, ZHANG Xianzhou. Effects of Grazing on Above- vs. Below-Ground Biom-ass Allocation of Alpine Grasslands on the Northern Tibetan Plateau[J]. Plos One, 2015, 10(8): e0135173.

    [29] 樸世龍, 方精云, 賀金生, 等. 中國草地植被生物量及其空間分布格局[J]. 植物生態(tài)學報, 2004, 28(4): 491–498.

    [30] SCARNECCHIA D L. Concepts of carrying capacity and substitution ratios: a systems viewpoint[J]. Journal of Range Management, 1990, 43(6): 553–555.

    [31] 陳百明. 中國農(nóng)業(yè)資源綜合生產(chǎn)能力和人口承載能力[M]. 北京: 氣象出版社, 2001.

    [32] 樊江文, 邵全琴, 王軍邦, 等. 三江源草地載畜壓力時空動態(tài)分析[J]. 中國草地學報, 2011, 33(3): 64–72.

    [33] 任繼周. 草原調查與規(guī)劃[M]. 北京: 中國農(nóng)業(yè)出版社, 1985.

    [34] 中華人民共和國農(nóng)業(yè)行業(yè)標準. 天然草地合理載畜量的計算(NY/T635-2002 )[S].

    [35] 辛玉春, 李曉明. 海西州草地資源生產(chǎn)力評價[J]. 青海草業(yè), 2016, 25(1): 23–29.

    [36] 武彥朋, 陳克龍, 張斐, 等. 青海湖流域典型草地物種豐富度與生產(chǎn)力的關系[J]. 生態(tài)學雜志, 2011, 30(7): 1449– 1453.

    [37] 龍慧靈, 李曉兵, 黃玲梅, 等. 內蒙古草原生態(tài)系統(tǒng)凈初級生產(chǎn)力及其與氣候的關系[J]. 植物生態(tài)學報, 2010, 34(7): 781–791.

    [38] 周偉, 王倩, 章超斌, 等. 黑河中上游草地NDVI時空變化規(guī)律及其對氣候因子的響應分析[J]. 草業(yè)學報, 2013, 22(1): 138–147.

    [39] 陸晴, 吳紹洪, 趙東升. 1982~2013年青藏高原高寒草地覆蓋變化及與氣候之間的關系[J]. 地理科學, 2017, 37(2): 292–300.

    [40] 武建雙, 李曉佳, 沈振西, 等. 藏北高寒草地樣帶物種多樣性沿降水梯度的分布格局[J]. 草業(yè)學報, 2012, 21(3): 17–25.

    [41] 青海省統(tǒng)計局. 青海省統(tǒng)計年鑒(1985-2017) [M]. 北京: 中國統(tǒng)計出版社, 2017: 279–280.

    Estimating grassland yield and carrying capacity in Qinghai Lake Basin based on MODIS NPP data

    Wang Qi1,2, Wu Chengyong1,2, Chen Kelong2,*, Zhang Xiao3, ZHANG Lele1,2, Ding Junxia1,2

    1. College of Geography, Qinghai Normal University, Xining 810008, China 2. Key Laboratory of Natural Geography and Environmental Processes of Qinghai Province, Xining 810008, China 3. College of Geography and Resource Science, Sichuan Normal University, Chengdu 610101, China

    Fast evaluating grassland productivity is critically significant to make a balancing development strategy between grassland and livestock in pastoral areas, as well as promoting the healthy and sustainable development of grassland ecosystems. Based on the NPP products of MODDIS, the grassland productivity was estimated in the Qinghai Lake Basin. The results showed that: (1)Theof the estimated grass yield per unit of alpine meadow, alpine grassland and warm grassland with the highest edible degree of cattle and sheep was 26.15 g?m-2,which indicated that the method could be used to quickly estimate the grassland yield in a regional scale. (2) The annual hay yield in the Qinghai Lake Basin was 1454200 tons, of which 591800 tons could be directly eaten by cattle and sheep; theoretical stock-carrying capacity was 810,700 sheep unit. (3) In the areas below the altitude of 3500 m, the accumulated temperature was the main factor that affected hay production per unit, and their correlation was negative. While above 3500 m, temperature was the main factor; the relationship between the factor of temperature and the hay production per unit was positive. The research results can provide scientific basis and supportive strategies for the planning of livestock production activities and ecological compensation in the alpine regions.

    grass yield;grassland productivity; vegetation net primary productivity; Qinghai Lake Basin

    10.14108/j.cnki.1008-8873.2019.04.024

    S127; S283

    A

    1008-8873(2019)04-178-08

    2018-07-14;

    2019-06-26

    國家自然科學基金項目(41661023); 國家社會科學基金(14XMZ072); 青海省科技計劃項目(2017-ZJ-782)

    王琪(1992—), 男, 甘肅平?jīng)鋈? 現(xiàn)為青海師范大學在讀碩士研究生, 主要從事生物地理與自然保護方面研究, E-mail: 1186599488@qq.com

    陳克龍, 男, 教授, 博士生導師, 主要從事生物地理與濕地生態(tài)方面的研究, E-mail: ckl7813@163.com

    王琪, 吳成永, 陳克龍, 等. 基于MODIS NPP數(shù)據(jù)的青海湖流域產(chǎn)草量與載畜量估算研究[J]. 生態(tài)科學, 2019, 38(4): 178-185.

    Wang Qi, Wu Chengyong, Chen Kelong, et al. Estimating grassland yield and carrying capacity in Qinghai Lake Basin based on MODIS NPP data[J]. Ecological Science, 2019, 38(4): 178-185.

    猜你喜歡
    載畜量產(chǎn)草量溫性
    瑪曲天然草場載畜量分析及過載放牧應對策略*
    甘肅科技(2023年9期)2023-10-30 01:31:40
    研究認為溫血動物出現(xiàn)于2.33億年前
    科學24小時(2022年10期)2022-10-11 01:50:06
    露天開采背景下錫林浩特市產(chǎn)草量時空變化及對氣象因子的響應
    2016年~2020年新源縣不同草地類草原載畜量的變化及經(jīng)濟效益分析
    草食家畜(2021年3期)2021-06-08 03:18:16
    論酒性綻放與人性澄明的四重關聯(lián)
    呼倫貝爾草原載畜量偏低原因探究
    黑龍江省大慶地區(qū)20個紫花苜蓿引種試驗
    一種含轉動副間隙多連桿機構非線性動力學行為分析方法
    施肥和補播對高寒草甸草原載畜能力的影響
    草原與草坪(2018年6期)2019-01-07 01:57:10
    苜蓿新品系內農(nóng)1號農(nóng)藝性狀的研究
    草原與草坪(2013年3期)2013-10-22 07:14:02
    欧美激情极品国产一区二区三区| 久久精品aⅴ一区二区三区四区| 免费观看精品视频网站| 高清av免费在线| www.熟女人妻精品国产| 色播在线永久视频| 黄色女人牲交| 男人舔女人的私密视频| 男女下面进入的视频免费午夜 | 亚洲人成77777在线视频| 一级片免费观看大全| 日本一区二区免费在线视频| av片东京热男人的天堂| 久久影院123| 在线观看免费高清a一片| 亚洲欧美日韩另类电影网站| 男女午夜视频在线观看| 亚洲第一青青草原| 琪琪午夜伦伦电影理论片6080| 老鸭窝网址在线观看| 视频区图区小说| 精品欧美一区二区三区在线| 人妻丰满熟妇av一区二区三区| 色婷婷久久久亚洲欧美| 午夜视频精品福利| 久久亚洲精品不卡| 亚洲性夜色夜夜综合| 一区二区三区精品91| 色综合欧美亚洲国产小说| 99国产精品99久久久久| 亚洲欧美激情在线| 精品电影一区二区在线| 午夜福利一区二区在线看| 午夜亚洲福利在线播放| 大香蕉久久成人网| 两人在一起打扑克的视频| 成人18禁高潮啪啪吃奶动态图| 亚洲精品av麻豆狂野| 水蜜桃什么品种好| 神马国产精品三级电影在线观看 | 麻豆久久精品国产亚洲av | 电影成人av| bbb黄色大片| svipshipincom国产片| 黄色怎么调成土黄色| 色婷婷久久久亚洲欧美| 国产精品综合久久久久久久免费 | 久久欧美精品欧美久久欧美| 国产不卡一卡二| 露出奶头的视频| 国产精品偷伦视频观看了| 成人国语在线视频| 久久久久久久久免费视频了| 9热在线视频观看99| 亚洲 国产 在线| 欧美黄色片欧美黄色片| 久久精品国产亚洲av高清一级| 热99re8久久精品国产| 制服人妻中文乱码| 午夜福利影视在线免费观看| 少妇粗大呻吟视频| 国产主播在线观看一区二区| 成在线人永久免费视频| 别揉我奶头~嗯~啊~动态视频| 高清av免费在线| 男男h啪啪无遮挡| 国产精品国产高清国产av| 午夜精品国产一区二区电影| 日韩欧美一区二区三区在线观看| 丝袜美腿诱惑在线| 国产欧美日韩一区二区精品| 欧美中文综合在线视频| 美女福利国产在线| 亚洲熟女毛片儿| 97碰自拍视频| 国产又爽黄色视频| 91成年电影在线观看| 免费观看人在逋| 桃红色精品国产亚洲av| 三级毛片av免费| 这个男人来自地球电影免费观看| 一个人观看的视频www高清免费观看 | 中文字幕色久视频| 欧美另类亚洲清纯唯美| 欧美最黄视频在线播放免费 | 不卡一级毛片| 成人永久免费在线观看视频| 欧美国产精品va在线观看不卡| 可以免费在线观看a视频的电影网站| 可以免费在线观看a视频的电影网站| 欧美最黄视频在线播放免费 | 男女做爰动态图高潮gif福利片 | 欧美老熟妇乱子伦牲交| 长腿黑丝高跟| 久久久国产精品麻豆| 中亚洲国语对白在线视频| 高清黄色对白视频在线免费看| 久久天堂一区二区三区四区| 日本黄色日本黄色录像| 国产精品一区二区三区四区久久 | 在线看a的网站| 青草久久国产| 免费不卡黄色视频| 国产精品一区二区在线不卡| 免费一级毛片在线播放高清视频 | 19禁男女啪啪无遮挡网站| 一级毛片精品| 亚洲av熟女| 狠狠狠狠99中文字幕| 一个人观看的视频www高清免费观看 | 1024视频免费在线观看| 久久午夜亚洲精品久久| 日本三级黄在线观看| 久久国产亚洲av麻豆专区| 日韩欧美在线二视频| 国产精华一区二区三区| 日韩免费av在线播放| 80岁老熟妇乱子伦牲交| 久久久国产精品麻豆| 香蕉国产在线看| 老司机午夜福利在线观看视频| 欧美久久黑人一区二区| 久久国产精品影院| 亚洲视频免费观看视频| 色精品久久人妻99蜜桃| 日韩欧美在线二视频| 少妇 在线观看| 搡老岳熟女国产| 久久香蕉国产精品| 黑丝袜美女国产一区| 国产成人影院久久av| 亚洲成国产人片在线观看| av天堂久久9| 9热在线视频观看99| 国产不卡一卡二| 亚洲性夜色夜夜综合| 国产熟女午夜一区二区三区| 女人被狂操c到高潮| 黄色毛片三级朝国网站| 水蜜桃什么品种好| 亚洲第一av免费看| 婷婷六月久久综合丁香| 久久中文字幕人妻熟女| 母亲3免费完整高清在线观看| 很黄的视频免费| 亚洲精品久久午夜乱码| 中文字幕精品免费在线观看视频| 黄色 视频免费看| 久久九九热精品免费| 十分钟在线观看高清视频www| 欧美日韩亚洲综合一区二区三区_| 成人免费观看视频高清| av欧美777| 国产一区二区三区在线臀色熟女 | 夜夜爽天天搞| 99香蕉大伊视频| 亚洲av片天天在线观看| 婷婷六月久久综合丁香| 成人手机av| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 亚洲精品中文字幕一二三四区| 亚洲国产中文字幕在线视频| 99国产精品99久久久久| 亚洲一区二区三区不卡视频| 成人亚洲精品一区在线观看| 久久国产精品男人的天堂亚洲| 一级黄色大片毛片| 免费一级毛片在线播放高清视频 | 91国产中文字幕| 视频区图区小说| av网站在线播放免费| 国内久久婷婷六月综合欲色啪| 国产日韩一区二区三区精品不卡| 久久国产亚洲av麻豆专区| 又黄又爽又免费观看的视频| 中文亚洲av片在线观看爽| 日韩欧美国产一区二区入口| 80岁老熟妇乱子伦牲交| 国产午夜精品久久久久久| 久久久久久亚洲精品国产蜜桃av| 日韩中文字幕欧美一区二区| 国产亚洲精品久久久久5区| 成人三级做爰电影| 欧美日韩亚洲高清精品| 人成视频在线观看免费观看| 日日摸夜夜添夜夜添小说| 欧美av亚洲av综合av国产av| 国产色视频综合| 久久国产亚洲av麻豆专区| 搡老岳熟女国产| 18美女黄网站色大片免费观看| 高潮久久久久久久久久久不卡| 久久久久亚洲av毛片大全| 免费人成视频x8x8入口观看| 看片在线看免费视频| 亚洲国产精品sss在线观看 | 老熟妇乱子伦视频在线观看| 日本黄色日本黄色录像| 午夜久久久在线观看| 久久久久久免费高清国产稀缺| videosex国产| 黄色丝袜av网址大全| 纯流量卡能插随身wifi吗| 性色av乱码一区二区三区2| 国产精品永久免费网站| 国产精品影院久久| 久久午夜亚洲精品久久| a级毛片黄视频| 亚洲成av片中文字幕在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品免费视频内射| 久久中文字幕一级| 亚洲欧美日韩高清在线视频| 在线观看免费视频网站a站| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜激情av网站| 国产午夜精品久久久久久| 成人亚洲精品av一区二区 | 久久久久久亚洲精品国产蜜桃av| 国产99白浆流出| 久久中文字幕一级| 可以在线观看毛片的网站| 亚洲国产看品久久| 亚洲成a人片在线一区二区| 亚洲欧美日韩另类电影网站| 精品国产乱码久久久久久男人| 国产日韩一区二区三区精品不卡| 亚洲国产看品久久| 免费不卡黄色视频| 国产成人欧美| 中亚洲国语对白在线视频| 看片在线看免费视频| 国产高清激情床上av| 99热国产这里只有精品6| 亚洲午夜精品一区,二区,三区| 在线天堂中文资源库| 欧美成人免费av一区二区三区| 亚洲精品在线美女| 亚洲aⅴ乱码一区二区在线播放 | 香蕉丝袜av| 国产av一区在线观看免费| 嫩草影视91久久| 日韩人妻精品一区2区三区| 一级毛片高清免费大全| 88av欧美| 久久伊人香网站| 久久久久亚洲av毛片大全| www.www免费av| 亚洲专区字幕在线| 成人免费观看视频高清| 人妻丰满熟妇av一区二区三区| 欧美av亚洲av综合av国产av| 69精品国产乱码久久久| 国产精品亚洲一级av第二区| 久久99一区二区三区| 亚洲色图综合在线观看| 国产一卡二卡三卡精品| 欧美成人免费av一区二区三区| 国产精品综合久久久久久久免费 | 亚洲五月天丁香| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 日本vs欧美在线观看视频| 亚洲精品成人av观看孕妇| 亚洲九九香蕉| 啦啦啦在线免费观看视频4| 两个人免费观看高清视频| 国产国语露脸激情在线看| 少妇 在线观看| 久久人妻福利社区极品人妻图片| 99国产精品99久久久久| 99在线人妻在线中文字幕| 高清av免费在线| 91字幕亚洲| www.www免费av| 国产精品久久久久成人av| 久久精品影院6| 夫妻午夜视频| 国产真人三级小视频在线观看| 无人区码免费观看不卡| 亚洲成人久久性| 最近最新免费中文字幕在线| 女性被躁到高潮视频| 亚洲av日韩精品久久久久久密| 国产日韩一区二区三区精品不卡| 日韩欧美一区二区三区在线观看| 精品久久蜜臀av无| 午夜影院日韩av| 69精品国产乱码久久久| 不卡一级毛片| 桃红色精品国产亚洲av| 亚洲五月色婷婷综合| 91av网站免费观看| 国产亚洲精品一区二区www| 一二三四社区在线视频社区8| 熟女少妇亚洲综合色aaa.| 波多野结衣高清无吗| 黄片大片在线免费观看| 亚洲精品成人av观看孕妇| 最近最新中文字幕大全免费视频| 亚洲专区字幕在线| 久久人妻福利社区极品人妻图片| 亚洲欧美一区二区三区久久| 国产一区二区激情短视频| 国产精品一区二区在线不卡| 黄色女人牲交| 久久伊人香网站| 久久国产精品男人的天堂亚洲| 很黄的视频免费| 国产成人欧美| 久9热在线精品视频| 变态另类成人亚洲欧美熟女 | 成人三级做爰电影| 日本欧美视频一区| 亚洲欧美精品综合久久99| 国产精品影院久久| 极品人妻少妇av视频| 丝袜在线中文字幕| 久久性视频一级片| 在线观看免费视频日本深夜| 一边摸一边抽搐一进一小说| 久热爱精品视频在线9| 黄频高清免费视频| 在线十欧美十亚洲十日本专区| 国产一区二区三区在线臀色熟女 | cao死你这个sao货| 国产精品国产高清国产av| 国产精品av久久久久免费| 黄色毛片三级朝国网站| 视频区欧美日本亚洲| 免费在线观看日本一区| 黄色视频不卡| 可以在线观看毛片的网站| 亚洲久久久国产精品| 欧美大码av| 久久天堂一区二区三区四区| 国产日韩一区二区三区精品不卡| 88av欧美| 国产精品电影一区二区三区| 亚洲精品美女久久av网站| 制服人妻中文乱码| 国产91精品成人一区二区三区| 成人18禁在线播放| 在线永久观看黄色视频| 国产片内射在线| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久| 日本 av在线| av网站在线播放免费| 日本欧美视频一区| 午夜久久久在线观看| 精品一区二区三卡| 美女扒开内裤让男人捅视频| 婷婷六月久久综合丁香| 制服人妻中文乱码| 国产精品二区激情视频| 久久青草综合色| 亚洲精品美女久久av网站| 亚洲精品久久午夜乱码| 国产亚洲欧美在线一区二区| 啦啦啦在线免费观看视频4| 国产一区二区三区综合在线观看| 我的亚洲天堂| 久久香蕉国产精品| 精品高清国产在线一区| 女性生殖器流出的白浆| 丝袜美足系列| 中文字幕人妻丝袜制服| 99国产精品一区二区三区| 欧美日韩乱码在线| 亚洲va日本ⅴa欧美va伊人久久| 18禁美女被吸乳视频| 亚洲专区字幕在线| 交换朋友夫妻互换小说| 亚洲五月色婷婷综合| 18禁裸乳无遮挡免费网站照片 | 成人18禁高潮啪啪吃奶动态图| 日韩欧美一区二区三区在线观看| 两个人看的免费小视频| 亚洲成av片中文字幕在线观看| 午夜免费鲁丝| 久久久久久大精品| 中文字幕精品免费在线观看视频| 免费人成视频x8x8入口观看| 亚洲人成网站在线播放欧美日韩| 如日韩欧美国产精品一区二区三区| 精品国产亚洲在线| 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人av教育| av天堂在线播放| 性欧美人与动物交配| 国产在线观看jvid| www.精华液| 一区二区日韩欧美中文字幕| 91大片在线观看| 宅男免费午夜| 久久青草综合色| av有码第一页| 久久精品亚洲精品国产色婷小说| 精品福利永久在线观看| 日韩中文字幕欧美一区二区| 亚洲人成电影免费在线| 精品人妻1区二区| 久久久国产成人精品二区 | 黄色a级毛片大全视频| 亚洲成av片中文字幕在线观看| av国产精品久久久久影院| av有码第一页| 高清欧美精品videossex| 免费女性裸体啪啪无遮挡网站| 岛国在线观看网站| 久久久久国内视频| 亚洲精品中文字幕在线视频| 日本wwww免费看| 精品久久久久久,| 一级片'在线观看视频| 一区二区日韩欧美中文字幕| 淫妇啪啪啪对白视频| 老汉色av国产亚洲站长工具| 亚洲午夜理论影院| 国产精品久久久av美女十八| 美女扒开内裤让男人捅视频| 午夜福利,免费看| xxxhd国产人妻xxx| av在线天堂中文字幕 | 91成年电影在线观看| 咕卡用的链子| 国产成人精品无人区| 亚洲欧美精品综合一区二区三区| 波多野结衣一区麻豆| 嫩草影视91久久| 欧美中文日本在线观看视频| 80岁老熟妇乱子伦牲交| 免费在线观看日本一区| 俄罗斯特黄特色一大片| 老司机在亚洲福利影院| 久久九九热精品免费| 大型黄色视频在线免费观看| 少妇的丰满在线观看| 亚洲片人在线观看| 国产成人av激情在线播放| 久久香蕉国产精品| 欧美性长视频在线观看| 日韩大码丰满熟妇| 纯流量卡能插随身wifi吗| e午夜精品久久久久久久| 丁香欧美五月| 亚洲欧洲精品一区二区精品久久久| 美女高潮喷水抽搐中文字幕| 看黄色毛片网站| 午夜精品在线福利| 老司机靠b影院| 国产国语露脸激情在线看| 村上凉子中文字幕在线| www.精华液| 国产精品二区激情视频| 91九色精品人成在线观看| 亚洲国产精品一区二区三区在线| 亚洲欧美激情在线| 免费在线观看黄色视频的| 91成人精品电影| 免费av中文字幕在线| 十八禁网站免费在线| 黑人操中国人逼视频| 在线看a的网站| 欧美成人午夜精品| 真人一进一出gif抽搐免费| 日韩免费av在线播放| 久久精品国产99精品国产亚洲性色 | 国产激情久久老熟女| 久久久久久人人人人人| 久久人人97超碰香蕉20202| 这个男人来自地球电影免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利免费观看在线| a级毛片黄视频| 欧美精品一区二区免费开放| 日本五十路高清| 亚洲男人的天堂狠狠| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 国产精品久久视频播放| 天天躁夜夜躁狠狠躁躁| 精品久久久久久成人av| 国产欧美日韩精品亚洲av| 91九色精品人成在线观看| av超薄肉色丝袜交足视频| 亚洲成人久久性| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 欧美日韩国产mv在线观看视频| 男女下面插进去视频免费观看| 日本黄色视频三级网站网址| 久久久精品国产亚洲av高清涩受| 日韩 欧美 亚洲 中文字幕| 91大片在线观看| 午夜福利在线免费观看网站| 777久久人妻少妇嫩草av网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲第一av免费看| 人人妻人人添人人爽欧美一区卜| 俄罗斯特黄特色一大片| 国产亚洲精品久久久久久毛片| 热99re8久久精品国产| 午夜a级毛片| 一级毛片女人18水好多| 欧美日韩乱码在线| 久久久久久免费高清国产稀缺| av片东京热男人的天堂| 宅男免费午夜| av电影中文网址| 这个男人来自地球电影免费观看| 99国产极品粉嫩在线观看| 夫妻午夜视频| 成人亚洲精品一区在线观看| 久热爱精品视频在线9| 黄色片一级片一级黄色片| avwww免费| 久久青草综合色| 亚洲七黄色美女视频| 在线观看www视频免费| av片东京热男人的天堂| 亚洲第一青青草原| 韩国av一区二区三区四区| 丝袜人妻中文字幕| 涩涩av久久男人的天堂| 久久影院123| 叶爱在线成人免费视频播放| 丰满饥渴人妻一区二区三| aaaaa片日本免费| 久久精品国产清高在天天线| 女人爽到高潮嗷嗷叫在线视频| 啪啪无遮挡十八禁网站| 自线自在国产av| 免费av中文字幕在线| 午夜福利免费观看在线| 国产精品亚洲一级av第二区| 桃色一区二区三区在线观看| 午夜福利,免费看| 免费在线观看视频国产中文字幕亚洲| 制服诱惑二区| 一a级毛片在线观看| 可以在线观看毛片的网站| 精品福利永久在线观看| 韩国精品一区二区三区| www.熟女人妻精品国产| 国产一区二区三区在线臀色熟女 | 无人区码免费观看不卡| 日韩欧美一区二区三区在线观看| 国产亚洲精品综合一区在线观看 | 日韩精品青青久久久久久| 久久伊人香网站| 国产一区二区三区在线臀色熟女 | 极品人妻少妇av视频| 18禁国产床啪视频网站| 国产免费男女视频| 午夜精品久久久久久毛片777| av欧美777| 亚洲av片天天在线观看| 亚洲五月婷婷丁香| 巨乳人妻的诱惑在线观看| 脱女人内裤的视频| 精品福利永久在线观看| 波多野结衣高清无吗| 91大片在线观看| 男女之事视频高清在线观看| 国产精品综合久久久久久久免费 | 免费一级毛片在线播放高清视频 | 国产日韩一区二区三区精品不卡| 午夜亚洲福利在线播放| 麻豆久久精品国产亚洲av | 国产欧美日韩精品亚洲av| 亚洲一区中文字幕在线| 国产欧美日韩综合在线一区二区| 国产av一区二区精品久久| 亚洲av成人一区二区三| 两性夫妻黄色片| 国产不卡一卡二| 一区二区三区国产精品乱码| 在线观看午夜福利视频| 亚洲国产毛片av蜜桃av| xxxhd国产人妻xxx| 亚洲精华国产精华精| 国产野战对白在线观看| 精品国内亚洲2022精品成人| 中亚洲国语对白在线视频| 午夜影院日韩av| 国产99久久九九免费精品| 母亲3免费完整高清在线观看| 日韩精品青青久久久久久| 亚洲情色 制服丝袜| 99国产精品一区二区蜜桃av| 欧美+亚洲+日韩+国产| 久久久久精品国产欧美久久久| 精品国产乱码久久久久久男人| 12—13女人毛片做爰片一| 在线播放国产精品三级| 国产有黄有色有爽视频| 老熟妇仑乱视频hdxx| 欧美大码av| 高清黄色对白视频在线免费看| 亚洲av成人不卡在线观看播放网| 最近最新中文字幕大全电影3 | 国产又色又爽无遮挡免费看| 亚洲av片天天在线观看| 制服诱惑二区| av视频免费观看在线观看| 一本综合久久免费| 日韩有码中文字幕| 亚洲精品一二三| 欧美国产精品va在线观看不卡| 国产成人一区二区三区免费视频网站| 亚洲男人天堂网一区|