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

    阿爾金山草地鼠害發(fā)生區(qū)及鼠荒地面積遙感估算

    2016-12-02 06:31:16李培先鄭江華倪亦非吳建國(guó)吾買(mǎi)爾吾守艾合買(mǎi)江阿吉那松曹克圖
    新疆農(nóng)業(yè)科學(xué) 2016年7期
    關(guān)鍵詞:阿爾金山鼠害荒地

    李培先 ,鄭江華,3,倪亦非 ,吳建國(guó) ,吾買(mǎi)爾·吾守 ,艾合買(mǎi)江·阿吉 ,那松曹克圖

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046; 2.新疆大學(xué)智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046;3.綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046; 4.新疆維吾爾自治區(qū)治蝗滅鼠指揮部辦公室,烏魯木齊 830001;5.新疆巴音郭楞蒙古自治州草原工作站,新疆庫(kù)爾勒 841000)

    ?

    阿爾金山草地鼠害發(fā)生區(qū)及鼠荒地面積遙感估算

    李培先1,2,鄭江華1,2,3,倪亦非4,吳建國(guó)4,吾買(mǎi)爾·吾守5,艾合買(mǎi)江·阿吉5,那松曹克圖5

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,烏魯木齊 830046; 2.新疆大學(xué)智慧城市與環(huán)境建模普通高校重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046;3.綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,烏魯木齊 830046; 4.新疆維吾爾自治區(qū)治蝗滅鼠指揮部辦公室,烏魯木齊 830001;5.新疆巴音郭楞蒙古自治州草原工作站,新疆庫(kù)爾勒 841000)

    【目的】運(yùn)用“3S”技術(shù)對(duì)阿爾金山草地鼠害發(fā)生區(qū)和鼠荒地進(jìn)行監(jiān)測(cè)和解譯,獲得鼠害發(fā)生面積和分布,為研究區(qū)鼠害蔓延的防治及鼠荒地的治理提供參考。【方法】研究區(qū)Landsat8遙感影像,采用最大似然監(jiān)督分類(lèi)方法,初步分出鼠害發(fā)生區(qū)、草地、裸地、水體等地物,再添加海拔高度和植被覆蓋度的訓(xùn)練規(guī)則,運(yùn)用決策樹(shù)分類(lèi)方法,進(jìn)一步從鼠害發(fā)生區(qū)中提取鼠荒地、鼠害發(fā)生區(qū),鼠害潛在發(fā)生區(qū),最后生成專(zhuān)題圖?!窘Y(jié)果】遙感估算得阿爾金山地區(qū)鼠害發(fā)生區(qū)面積為3 419.621 km2,其中鼠荒地面積1 376.476 km2,分別占到研究區(qū)的6.92%和2.78%,與2014年11月新疆巴州草原站實(shí)地調(diào)查估算結(jié)果相近,且有擴(kuò)大趨勢(shì)?!窘Y(jié)論】鼠荒地或鼠害潛在發(fā)生區(qū)有沿水源邊緣分布特征,鼠荒地與草地,裸地多呈鑲嵌分布。鼠荒地和裸地、草地共存易形成混合像元,造成解譯誤差;最大似然法和決策樹(shù)方法的結(jié)合使結(jié)果更符合實(shí)際調(diào)查。

    “3S”技術(shù); 阿爾金山;監(jiān)督分類(lèi);決策樹(shù);鼠荒地

    0 引 言

    【研究意義】草原鼠害是有害嚙齒類(lèi)動(dòng)物對(duì)草原資源及草原畜牧業(yè)產(chǎn)生的危害性影響[1]。草原鼠荒地指主要因嚙齒類(lèi)動(dòng)物活動(dòng)和超載過(guò)牧等原因引起草原嚴(yán)重退化的次生裸地,其植被覆蓋度低于20%,地上植物生物量低于原生草原的20%[2]。草原鼠害監(jiān)測(cè)預(yù)報(bào)僅依靠人工的采樣調(diào)查雖然有較高的準(zhǔn)確性,但費(fèi)時(shí)耗力,存在時(shí)效性差等缺點(diǎn),難以適應(yīng)大范圍實(shí)時(shí)監(jiān)測(cè)的需求。結(jié)合多種現(xiàn)代技術(shù)有效監(jiān)測(cè)鼠害發(fā)生區(qū)及鼠荒地面積與分布對(duì)研究區(qū)鼠害蔓延的防治、鼠荒地的治理以及草原植被生態(tài)恢復(fù)等具有重要意義?!厩叭搜芯窟M(jìn)展】一些學(xué)者已運(yùn)用“3S”技術(shù)(Geographic Information System,GIS;Remote Sensing,RS;Global Positional System,GPS)進(jìn)行鼠害監(jiān)測(cè):黃建文等[3]利用TM影像,分析了天然梭梭林大沙鼠鼠害防治前后的土壤調(diào)節(jié)植被指數(shù)和植被生長(zhǎng)狀況,為鼠害防治管理和動(dòng)態(tài)監(jiān)測(cè)提供科學(xué)依據(jù)。何詠琪等[4]運(yùn)用“3S”技術(shù),確定了高程、坡度、坡向、草地類(lèi)型、土壤類(lèi)型、5~10月增強(qiáng)型植被指數(shù)6個(gè)鼠害監(jiān)測(cè)模型因子,計(jì)算因子權(quán)重值和量化值,建立鼠害監(jiān)測(cè)模型,提取了不同鼠害發(fā)生區(qū)和鼠害危害區(qū)的閾值,最終生成青海省草原鼠害模擬分布區(qū),結(jié)果與實(shí)地調(diào)查的鼠害區(qū)域基本一致。徐正剛等[5]基于MODIS數(shù)據(jù),提出通過(guò)分析洞庭湖區(qū)植被指數(shù)變動(dòng)規(guī)律以確定評(píng)價(jià)時(shí)間的方法,選擇增強(qiáng)型植被指數(shù)EVI進(jìn)行鼠害的敏感性分析,確定了災(zāi)害閾值,最終采用植被指數(shù)差值法獲得了研究區(qū)作物危害分布圖與危害等級(jí)分布圖。王瑋等[6]指出滿(mǎn)足草地鼠害監(jiān)測(cè)要求的有TM和SPOT、Quick Bird等衛(wèi)星遙感數(shù)據(jù), TM影像是進(jìn)行草原鼠蟲(chóng)害的研究主要數(shù)據(jù)。阿爾金山鼠害發(fā)生區(qū)屬于新疆巴音郭楞蒙古自治州四大鼠害發(fā)生區(qū)之一[7]。該發(fā)生地分布有5個(gè)鼠害區(qū),害鼠有灰尾兔(Lepusoiostolus),大耳鼠兔(Ochotonamacrotis)、高原鼠兔(O.curzoniae)、柯氏鼠兔(O.koslowi)、喜馬拉雅旱獺(Marmotahimalayana)、白尾松田鼠(Microtus(Phaiomys)leucurus Blyth(1863))等 17 種,均為寒旱型種類(lèi)[8]。其中優(yōu)勢(shì)種白尾松田鼠有效洞口數(shù) 187.5 個(gè)/hm2,高原鼠兔有效洞口數(shù) 260個(gè)/hm2。相關(guān)高原鼠兔的生態(tài)學(xué)研究涉及與草甸植物群落的關(guān)系[9-12]、食物選擇[13]、冬季死亡率[14]和行為模式[15]等,以及阿爾金山自然保護(hù)區(qū)高原鼠兔夏季微生境選擇的主導(dǎo)因子分析[16]與食性顯微組織分析[17]。賈婷婷等[18]得出:隨高原鼠兔有效洞穴密度增加,樣方內(nèi)出現(xiàn)的植物種數(shù)逐漸增多,說(shuō)明高原鼠兔活動(dòng)增加了高寒草甸植物群落的物種豐富度,但增加物種主要是有毒的乳漿大戟、圓葉筋骨草、藍(lán)花棘豆,以及適口性差的矮火絨草和蘭石草,草甸質(zhì)量變劣。馬波等[19]得出:NDVI與洞穴數(shù)有顯著正相關(guān)的線(xiàn)性回歸關(guān)系:在無(wú)家畜競(jìng)爭(zhēng)的自然環(huán)境中, 高原鼠兔有選擇地利用植被較好的生境而避免使用植被過(guò)度退化的生境?!颈狙芯壳腥朦c(diǎn)】對(duì)高原鼠兔生態(tài)學(xué)特性的研究可為更準(zhǔn)確分析鼠害的空間分布提供參考依據(jù),然而對(duì)阿爾金山鼠害發(fā)生區(qū)及鼠荒地的發(fā)生面積以及空間分布的研究卻少有涉及。研究運(yùn)用“3S技術(shù)”對(duì)阿爾金山草地鼠害發(fā)生區(qū)和鼠荒地進(jìn)行監(jiān)測(cè)和解譯,獲得鼠害發(fā)生面積和分布?!緮M解決的關(guān)鍵問(wèn)題】研究結(jié)合實(shí)地采集的GPS數(shù)據(jù)等,參照鼠害發(fā)生海拔、植被覆蓋度的生境特點(diǎn),利用研究區(qū)的Landsat8 遙感影像,采用最大似然監(jiān)督分類(lèi)方法和決策樹(shù)分類(lèi)方法提取鼠荒地、鼠害發(fā)生區(qū),鼠害潛在發(fā)生區(qū)并統(tǒng)計(jì)其面積,分析鼠荒地和鼠害發(fā)生區(qū)的空間分布特征。

    1 材料與方法

    1.1 研究區(qū)概況

    阿爾金山鼠害發(fā)生區(qū)主要位于阿爾金山自然保護(hù)區(qū)境內(nèi)低于海拔4 800 m的區(qū)域,行政區(qū)劃隸屬于巴音郭楞蒙古自治州若羌縣和且末縣境內(nèi)。在保護(hù)區(qū)中是以高原有蹄類(lèi)野生動(dòng)物為主的高原荒漠生態(tài)系統(tǒng)[20]。年平均降雨量 100~300 mm,主要集中在 6~8 月,形成明顯的雨季。年平均氣溫0~1 ℃,最熱月平均氣溫8 ℃,最低溫度在-30 ℃以下;土壤類(lèi)型主要有山地棕漠土、高寒荒漠土、高寒荒漠草原土、高寒草原土和高寒草甸土等[16]。該發(fā)生區(qū)海拔在3 500 m以上,植物異常稀疏。4 000 m 以上的高山帶為墊狀優(yōu)若黎、藏艾菊構(gòu)成的寒漠,4 000 m以下的溝谷有水柏枝和小蔚構(gòu)成的稀疏灌叢[21]。山地荒漠帶的上限高達(dá) 3 100~3 300 m,以合頭草荒漠為主,蒿屬荒漠次之[7]。主要植被類(lèi)型包括紫花針茅、墊狀駝絨藜、硬葉苔草,以及棘豆、蒿子、冰草、草地早熟禾、芨芨草、合頭草、黑枸杞、沙拐棗、琵琶柴、駝絨藜、鳶尾等。隨著海拔升高和植被類(lèi)型變化草地類(lèi)型依次為、溫性荒漠草原、荒漠草原、高山草原、高山草原和高寒草原的過(guò)渡帶、高寒草原。表1,圖1

    表1 阿爾金山鼠害發(fā)生研究區(qū)坐標(biāo)范圍

    注:研究區(qū)確定依據(jù):1)包含實(shí)地調(diào)查路線(xiàn)和采樣點(diǎn),包含2015年新疆畜牧廳治蝗辦提供的阿爾金山鼠害發(fā)生區(qū)。2)由高原鼠兔生境多發(fā)生在3 500~4 000 m海拔高度,研究區(qū)需覆蓋阿爾金山該海拔區(qū)間

    Note: the study range was determined by two reasons:(1) The study area need to contain the routes of field survey, sampling points, and the data of rodent damage rangeland in Altun Mountain grassland in 2015 provided by Locust and Rodent Control Headquarters of Xinjiang. (2) The study area need to cover the altitude interval from 3500m to 4000m in Altun Mountain grassland where the habitats of plateau pika(Ochotonacurzoniae) exist

    注:研究區(qū)范圍(綠色邊框);鼠害發(fā)生區(qū)估算范圍,由巴州草原站和治蝗辦提供(淡藍(lán)色面);野外調(diào)查路線(xiàn)和采樣點(diǎn)(紅線(xiàn)和黃色標(biāo)記);高原鼠兔和鼠害發(fā)生區(qū)照片(紅色邊框)

    Note: the scope of the study area (green border).The estimation range of rodent damage rangeland was provided by Grassland Workstation of Bayinguoleng Mongolian Autonomous Prefecture (light blue area). Routes of field investigation and sampling points were shown with red lines and yellow marks separately. The photos of Plateau pika (Ochotonacurzoniae) and rodent damage rangeland were presented with red border

    圖1 研究區(qū)位置及采樣點(diǎn)分布

    Fig.1 The location and sampling point of study area

    1.2 材 料

    1.2.1 遙感數(shù)據(jù)及預(yù)處理

    遙感數(shù)據(jù)為L(zhǎng)andsat8 OLI數(shù)據(jù),DEM數(shù)字高程數(shù)據(jù)為ASTER GDEM 30 m數(shù)據(jù)。表2

    表2 研究區(qū)影像數(shù)據(jù)信息

    遙感影像空間分辨率為30 m,時(shí)間分辨率為16 d,地圖投影坐標(biāo)系為 UTM-WGS84。影像選擇考慮到其適用性(研究區(qū)云量<10%,無(wú)陰影條帶),與調(diào)查時(shí)間2015年9月13~15日最接近的影像。選擇日期在9月14日和9月5日的影像覆蓋研究區(qū)的主體部分,與調(diào)查時(shí)間很接近。而8月22日與8月27日影像僅覆蓋在研究區(qū)的兩個(gè)邊沿,與調(diào)查時(shí)間間隔小于1個(gè)月,對(duì)研究區(qū)鼠荒地解譯的影響可以接受。

    歷史野外數(shù)據(jù):新疆畜牧廳治蝗辦提供的2014年阿爾金山鼠害發(fā)生區(qū)和鼠荒地GPS數(shù)據(jù)。

    圖像處理軟件為ENVI5.1,ArcGIS10.0和Google Earth Free 7.1.4。

    遙感數(shù)據(jù)預(yù)處理:對(duì)已進(jìn)行系統(tǒng)輻射校正和地面控制點(diǎn)幾何校正的Landsat8 OLI的影像進(jìn)行輻射定標(biāo),圖像鑲嵌與裁剪,F(xiàn)LAASH大氣校正;對(duì)14景DEM數(shù)據(jù)進(jìn)行鑲嵌與裁剪。

    1.2.2 實(shí)地調(diào)查數(shù)據(jù)

    經(jīng)實(shí)地調(diào)查,鼠害發(fā)生區(qū)的優(yōu)勢(shì)鼠種為高原鼠兔,處于高寒草原和高山草原的過(guò)渡帶,海拔基本在3 500 m以上,鼠害發(fā)生區(qū)主要植被類(lèi)型:針茅、硬葉苔草、蒿子、冰草、草地早熟禾。徐瑤[22]得出在草場(chǎng)發(fā)生輕度退化的樣方中,鼠害的發(fā)生率為 100%。草場(chǎng)退化加劇,鼠害發(fā)生的概率則隨之減少。與此次調(diào)查的情況相同:在此次調(diào)查中發(fā)現(xiàn)的鼠害發(fā)生區(qū),植被覆蓋度都較高,部分達(dá)到70%,對(duì)鼠害發(fā)生區(qū)不同草種光譜數(shù)據(jù)的采集與分析,發(fā)現(xiàn)有較明顯退化跡象。采用GPS采集研究區(qū)鼠害發(fā)生地、草地、鼠荒地、裸地等地物樣本,作為訓(xùn)練樣本輔助后期遙感影像解譯。

    1.3 方 法

    1.3.1 最大似然監(jiān)督分類(lèi)(MLC)

    基于參數(shù)化密度分布函數(shù)判別的最大似然方法MLC(Maximum Like-hood Classification)是遙感影像監(jiān)督分類(lèi)最常用的統(tǒng)計(jì)方法之一[23-24]。最大似然法與其他非參數(shù)方法比較具有易于與先驗(yàn)知識(shí)融合和算法簡(jiǎn)單而便于實(shí)施等優(yōu)點(diǎn)[25],且具有較高的分類(lèi)精度。

    根據(jù)此次調(diào)查的實(shí)際情況和影像特點(diǎn),將研究區(qū)地物分為:鼠荒地(粗分)、草地、裸地、水體、雪地、云和陰影。通過(guò)GPS數(shù)據(jù)、結(jié)合Google Earth進(jìn)行訓(xùn)練樣本的選擇。具體為將部分實(shí)地采集的草地、裸地的GPS面作為驗(yàn)證樣本,用于精度評(píng)價(jià);將實(shí)地采集的鼠害發(fā)生區(qū)樣本,結(jié)合巴州治蝗辦提供的鼠荒地坐標(biāo)點(diǎn)在Google Earth繪制60 m×60 m鼠荒地樣本, 隨機(jī)均勻選擇一半為訓(xùn)練樣本,另一半為驗(yàn)證樣本,運(yùn)用最大似然監(jiān)督分類(lèi)初步進(jìn)行阿爾金山鼠荒地(粗分)解譯和精度驗(yàn)證。

    1.3.2 決策樹(shù)分類(lèi)

    決策樹(shù)方法主要是決策樹(shù)學(xué)習(xí)和決策樹(shù)分類(lèi)兩個(gè)過(guò)程。決策樹(shù)學(xué)習(xí)過(guò)程是通過(guò)對(duì)訓(xùn)練樣本進(jìn)行歸納學(xué)習(xí)(Inductive learning) ,生成以決策樹(shù)形式表示的分類(lèi)規(guī)則的機(jī)器學(xué)習(xí)(Machine learning)過(guò)程[26]。通過(guò)對(duì)訓(xùn)練樣本進(jìn)行決策樹(shù)學(xué)習(xí)生成決策樹(shù),決策樹(shù)可以根據(jù)屬性的取值對(duì)一個(gè)未知樣本集進(jìn)行分類(lèi),就是決策樹(shù)分類(lèi)[27-28]。

    巴州復(fù)雜的地理自然環(huán)境,為各種鼠類(lèi)提供了有利的生存條件,隨著海拔高度和草地類(lèi)型的垂直變化,鼠類(lèi)的分布呈現(xiàn)出明顯的垂直分異[7]。李葉等[16]指出對(duì)阿爾金山高寒荒漠草原高原鼠兔生境選擇的主成分分析表明影響高原鼠兔生境選擇的主要因子依次是海拔、蓋度、植物種數(shù)和距道路距離。說(shuō)明海拔高度是提取阿爾金山鼠害發(fā)生區(qū)的最重要因子。分析GPS數(shù)據(jù)對(duì)海拔高度的訓(xùn)練規(guī)則如下。

    (1)鼠害發(fā)生區(qū)(3 500~4 000 m):由此次調(diào)查,以及新疆畜牧廳治蝗辦、巴州治蝗辦提供的鼠害發(fā)生區(qū)GPS坐標(biāo)數(shù)據(jù)確定。(2)高海拔潛在發(fā)生區(qū) (4 000~4 400 m):指發(fā)現(xiàn)鼠害有局部發(fā)生的高海拔地區(qū)。由此次實(shí)際調(diào)查有鼠害局部發(fā)生(4 026~4 069 m),結(jié)合巴州草原站與治蝗辦提供的鼠害局部發(fā)生GPS坐標(biāo)(海拔4 409 m)確定。以及李葉等[16]于2012年7~8月在阿爾金山海拔高度為(4 122.25±139.48)m處做高原鼠兔穴居地樣本實(shí)驗(yàn)。鼠害潛在發(fā)生區(qū)說(shuō)明鼠害發(fā)生有局部擴(kuò)大趨勢(shì)。(3)低海拔潛在發(fā)生區(qū)(3 250~3 500 m): 指發(fā)現(xiàn)鼠害有局部發(fā)生的低海拔地區(qū)。由巴州草原站與治蝗辦提供的鼠害局部發(fā)生GPS坐標(biāo)確定(最低海拔3 238 m)。(4)未確定發(fā)生區(qū)(>4 250 m or <3 250 m):高原鼠兔生境選擇在植被高度、植被密度、坡度、海拔高度、植被類(lèi)型和土壤硬度等生態(tài)因子方面差異都有統(tǒng)計(jì)學(xué)意義(P<0.01)[16]。所以?xún)H依靠海拔高度無(wú)法確定是否發(fā)生鼠害,定為不確定發(fā)生區(qū)。在該海拔高度內(nèi)有適合高原鼠兔等鼠種適宜的生境,也可能變成鼠害潛在發(fā)生區(qū)。

    由草原鼠荒地“植被覆蓋度低于20%”的定義[2]作為鼠荒地的提取規(guī)則。

    2 結(jié)果與分析

    2.1 最大似然監(jiān)督分類(lèi)方法粗分鼠害發(fā)生區(qū)

    鼠荒地(粗分)包括鼠害發(fā)生區(qū),鼠害潛在(局部)發(fā)生區(qū),鼠荒地等,因?yàn)橛?xùn)練樣本中包含鼠害發(fā)生區(qū),鼠害潛在(局部)發(fā)生區(qū),鼠荒地信息;另外分類(lèi)中可能造成的鼠荒地與裸地混分,鼠害發(fā)生區(qū)與草地之間的混分等。作為初步分類(lèi)。

    從混淆矩陣中得出總體分類(lèi)精度為86.67%,Kappa系數(shù)為0.81。鼠荒地的生產(chǎn)者精度為85.37%;用戶(hù)精度低,為31.25%,錯(cuò)分誤差高?;旌舷裨斐闪耸蠡牡嘏c裸地、草地的錯(cuò)分現(xiàn)象。從影像解譯圖中鼠荒地空間分布上看,也說(shuō)明鼠荒地和草地,鼠荒地與裸地都有交集,鑲嵌分布,符合實(shí)際的調(diào)查情況。圖2

    2.2 決策樹(shù)分類(lèi)方法細(xì)分鼠害發(fā)生區(qū)和提取鼠荒地

    對(duì)鼠荒地(粗分)結(jié)果進(jìn)行基于海拔訓(xùn)練規(guī)則的決策樹(shù)分類(lèi),分離具體鼠害發(fā)生區(qū)和未確定發(fā)生區(qū),減小錯(cuò)分影響。進(jìn)一步對(duì)結(jié)果運(yùn)用決策樹(shù)分類(lèi),添加植被覆蓋度小于20%的訓(xùn)練規(guī)則,減小草地對(duì)鼠荒地造成的錯(cuò)分誤差,從鼠害發(fā)生區(qū)中提取鼠荒地。圖3

    圖2 2015年新疆阿爾金山地區(qū)鼠荒地潛在分布衛(wèi)星影像解譯

    圖3 2015新疆阿爾金山地區(qū)鼠害發(fā)生區(qū)分布衛(wèi)星影像解譯

    首先運(yùn)用基于NDVI的像元二分法選擇不同的閾值計(jì)算植被覆蓋度,用ArcGIS計(jì)算出從0%~10%到90%~100% 10個(gè)等級(jí);再結(jié)合實(shí)地調(diào)查目測(cè)法,照片比對(duì),進(jìn)行實(shí)際調(diào)查結(jié)果和不同閾值下計(jì)算結(jié)果的誤差分析。選擇的NDVIsoil和NDVIveg的閾值分別為累計(jì)貢獻(xiàn)率的5.2%和99.5%時(shí)各誤差率相對(duì)最低。表3

    然后繪出植被覆蓋度等級(jí)圖。圖4

    表3 不同閾值下基于NDVI像元二分法計(jì)算的草地覆蓋度及其誤差

    圖4 2015年新疆阿爾金山地區(qū)植被覆蓋等級(jí)

    最后,按決策樹(shù)方法進(jìn)一步得出鼠荒地解譯圖并進(jìn)行面積統(tǒng)計(jì)。得出2015年新疆阿爾金山地區(qū)鼠荒地空間分布圖,阿爾金山地區(qū)鼠害發(fā)生區(qū)面積為3 419.621 km2(512.943 17萬(wàn)畝),其中鼠荒地面積1 376.476 km2(206.471 43萬(wàn)畝),以及潛在鼠害發(fā)生區(qū)面積。將圖2、圖5與Google Earth影像進(jìn)行比對(duì)發(fā)現(xiàn)鼠荒地及鼠害潛在發(fā)生區(qū)有沿阿牙克庫(kù)木湖與阿其克庫(kù)勒湖周邊分布以及沿山地河流邊緣分布的特征。 圖5,表4

    2014年巴州草原站由實(shí)地調(diào)查GPS點(diǎn)勾畫(huà)初步計(jì)算的阿爾金山草原鼠害面積已達(dá)3 466.667 km2(520萬(wàn)畝),嚴(yán)重危害面積(鼠荒地)約1 333.333 km2(200萬(wàn)畝),且有繼續(xù)擴(kuò)大的趨勢(shì)。這與2015年遙感估算得到的面積相近,但鼠害發(fā)生區(qū)及鼠荒地在空間分布上并非成片的面狀分布(圖1),而是由生境等多因素決定的點(diǎn)狀或斑塊狀分布,呈現(xiàn)出與草地,裸地互有交集,鑲嵌分布特點(diǎn)。圖5

    2015年鼠荒地遙感估算面積比2014年多出3.24%,為43.143 km2(6.471 43萬(wàn)畝),尤其鼠害潛在(局部)發(fā)生或活動(dòng)區(qū)的大面積出現(xiàn),說(shuō)明鼠害有一定擴(kuò)大蔓延的趨勢(shì)。

    圖5 2015年新疆阿爾金山地區(qū)鼠荒地及鼠害發(fā)生區(qū)分布衛(wèi)星影像解譯

    鼠害地類(lèi)型Thetypeofrodentdamagerangelands面積(km2)Area面積(萬(wàn)畝)Area(104mu)占研究區(qū)比(%)Ratio鼠荒地137647620647143278鼠害發(fā)生區(qū)(3500~4000m)204314530647174413高海拔潛在發(fā)生區(qū)(4000~4400m)6334533950179951281低海拔潛在發(fā)生區(qū)(3250~3500m)87092313063842176未確定發(fā)生區(qū)(>4400mor<3250m)471203170680465953其它341078365116175466898

    注:研究區(qū)總面積為49 444.944 km2( 7 416.741 65萬(wàn)畝)。Landsat8 OLI影像單個(gè)像元的實(shí)際面積為900 m2≈0.001 km2故面積(km2)保留三位小數(shù)。實(shí)際工作中多使用(萬(wàn)畝)單位,故進(jìn)行換算;0.001 km2=0.000 15萬(wàn)畝,面積(萬(wàn)畝)保留五位小數(shù)

    Note: the whole study area is 49,444.944 km2(7,416.741,65×104mu).The single pixel area of Landsat8 OLI image is 900 m2≈0.001 km2, so the area/ km2keep three decimal. The unit of mu is used more often in practical work because of 0.001 km2equal to 0.000,15×104mu. So the area/ 104mu keep five decimal

    3 討 論

    此次監(jiān)測(cè)完成了研究區(qū)鼠害發(fā)生區(qū)和鼠荒地分布及面積的提取。但從影像方面看,30 m空間分辨率下只能通過(guò)適當(dāng)波段組合后不同地物的顏色和深淺來(lái)區(qū)別,造成訓(xùn)練樣本選擇的目視誤判。區(qū)分草地、裸地、鼠荒地的四景鑲嵌影像時(shí)間不同也對(duì)影像的自動(dòng)解譯造成誤差。由于阿爾金山鼠害研究區(qū)屬無(wú)人區(qū),海拔高,空氣較稀薄,山路崎嶇或無(wú)法通行,給樣本的采集帶來(lái)很大困難,造成嚴(yán)重鼠荒地訓(xùn)練樣本少,對(duì)鼠荒地的解譯精度造成一定影響??煽紤]引進(jìn)無(wú)人機(jī)航拍技術(shù)對(duì)鼠害區(qū)進(jìn)行監(jiān)測(cè),提高解譯精度。

    4 結(jié) 論

    4.1 從影像解譯效果和方法上:30 m中空間分辨率的landsat8 OLI數(shù)據(jù)針對(duì)鼠荒地和裸地的用戶(hù)精讀依然較低,最主要原因是鼠害發(fā)生區(qū)與草地、鼠荒地與裸地存在共存現(xiàn)象。30 m中空間分辨率的影像對(duì)阿爾金山草地和裸地區(qū)分基本可行,但共存現(xiàn)象和混合像元,造成較多錯(cuò)分。針對(duì)鼠荒地高原鼠兔為主要鼠害,進(jìn)行生境因子分析,將海拔高度以及植被覆蓋度作為訓(xùn)練規(guī)則,運(yùn)用最大似然監(jiān)督分類(lèi)和決策樹(shù)方法相結(jié)合提取鼠害發(fā)生區(qū)、鼠害潛在發(fā)生區(qū)以及鼠荒地。使結(jié)果更符合實(shí)際調(diào)查。

    4.2 從鼠害發(fā)生區(qū)及鼠荒地的空間分布來(lái)看,鼠荒地或鼠害潛在發(fā)生區(qū)有沿水源邊緣分布的特征,且鼠荒地與草地,裸地互相有交集。與黃倩等[1]得出的高原鼠兔和大耳鼠兔主要沿河谷階地,陽(yáng)坡或避風(fēng)的洼地分布,鼠害草原多呈鑲嵌分布的結(jié)論相近。鼠害潛在分布區(qū)(4 000~4 400 m)、(3 250~3 500 m)的局部鼠害發(fā)生,其危害程度和面積估算等需要進(jìn)一步實(shí)地驗(yàn)證工作,為有效預(yù)防鼠害蔓延擴(kuò)散提供依據(jù)。

    4.3 從鼠害防治上看,阿爾金山草地鼠害發(fā)生區(qū)和鼠荒地面積大并有一定擴(kuò)大趨勢(shì),對(duì)當(dāng)?shù)夭菰鷳B(tài)環(huán)境勢(shì)必造成負(fù)面影響,需進(jìn)行及時(shí)有效的鼠害監(jiān)測(cè)和相應(yīng)防治工作,同時(shí)鼠荒地宜進(jìn)行合理的草地植被生態(tài)修復(fù),以實(shí)現(xiàn)阿爾金山草原畜牧業(yè)的可持續(xù)發(fā)展。

    致謝:感謝新疆維吾爾自治區(qū)產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生基地—新疆維吾爾自治區(qū)畜牧廳治蝗滅鼠指揮部辦公室和新疆巴州草原站與治蝗辦的大力支持與協(xié)助!

    References)

    [1]黃倩,花立民,曹慧,等.甘肅草原鼠害區(qū)劃研究[J]. 草業(yè)科學(xué),2009,26(2):91-99.

    HUANG Qian, HUA Li-ming, CAO Hui, et al. (2009). Study on grassland rodent pest division in Gansu Province [J].PrataculturalScience, 26(2):91-99. (in Chinese)

    [2]唐川江,周俗,張新躍,等.四川省草原鼠荒地調(diào)查報(bào)告[J].草業(yè)科學(xué),2007,24(7):58-61.

    TANG Chuan-jiang, ZHOU Su, ZHANG Xin-yue, et al.(2007). Survey On rodent damaged rangeland in Sichuan Province [J].PrataculturalScience, 24(7):58-61. (in Chinese)

    [3]黃建文,鞠洪波,特木欽,等. 阿拉善左旗天然梭梭林鼠害防治的遙感監(jiān)測(cè)[J].林業(yè)科學(xué), 2004,40(3):107-110.

    HUANG Jian-wen, JU Hong-bo, TE Mu-qin ,et al.(2004).Monitoring on Controlling Effect of Rhombomys opimus in Haloxylon Ammodendron Bunge Stands Using Remote TM Imagery in Alashan [J].ScientiaSilvaeSinicae, 40(3):107-110. (in Chinese)

    [4]何詠琪, 黃曉東, 侯秀敏,等. 基于3S技術(shù)的草原鼠害監(jiān)測(cè)方法研究[J]. 草業(yè)學(xué)報(bào), 2013,22(3):33-40.

    HE Yong-qi,HUANG Xiao-dong,HOU Xiu-min, et al. (2013). Monitoring grassland rodents with 3S technologies [J].ActaPrataculturaeSinica, 22(3):33-40. (in Chinese)

    [5]徐正剛, 趙運(yùn)林, 李波,等. 基于MODIS植被指數(shù)評(píng)估洞庭湖區(qū)東方田鼠大暴發(fā)的危害[J]. 生態(tài)學(xué)報(bào), 2014,34(23):7 101-7 109.

    XU Zheng-gang, ZHAO Yun-lin, LI Bo, et al.(2014).Assessing Yangtze vole damage in Dongting Lake region of outbreak year based on MODIS imagery [J].AetaEeologieaSinica, 34(23):7,101-7,109. (in Chinese)

    [6]王瑋, 馮琦勝, 于惠,等. "3S"技術(shù)在草地鼠蟲(chóng)害監(jiān)測(cè)與預(yù)測(cè)中的應(yīng)用[J]. 草業(yè)科學(xué), 2010,27(3):31-39.

    WANG Wei, FENG Qi-sheng, YU Hui, et al.(2010). Application of "3S" technology on monitoring and prediction of the grassland rodents and insects [J].PrataculturalScience, 27(3):31-39.(in Chinese)

    [7]李文利. 新疆巴音郭楞蒙古自治州草地優(yōu)勢(shì)種害鼠發(fā)生危害區(qū)域劃分研究[J]. 草食家畜, 2013,(1):58-61.

    LI Wen-li.(2013).Distribution of Dominant Rodents on Grassland in the Bayinguoleng Mongulia Autonomous Prefecture of Xinjiang [J].Grass-FeedingLivestock, (1):58-61. (in Chinese)

    [8]鄭智民,姜志寬,陳安國(guó).嚙齒動(dòng)物學(xué)[M].上海:上海交通大學(xué)出版社,2008:83-114.

    ZHENG Zhi-min, JIANG Zhi-kuan, CHEN An-guo. (2008).RodentsZoology[M]. Shanghai: Shanghai Jiaotong University Press:83-114. (in Chinese)

    [9]孫飛達(dá), 龍瑞軍, 路承香. 高原鼠兔不同洞穴密度對(duì)高寒草地植物群落組成及多樣性的影響[J]. 干旱區(qū)資源與環(huán)境, 2010,24(7):181-186.

    SUN Fei-da, LONG Rui-jun, LU Chen-xiang. (2010). Effects of plateau pikas (Ochotonacurzoniae) burrow densities on plant community composition and population diversity in alpine meadow [J].JournalofAridLandResourcesandEnvironment, 24(7):181-186.(in Chinese)

    [10] 劉偉, 張毓, 王溪,等. 高原鼠兔刈割行為與棲息地植物群落的關(guān)系[J]. 獸類(lèi)學(xué)報(bào), 2009,29(1):40-49.

    LIU Wei, ZHANG Yu, WANG Xi, et al.(2009).The relationship of the harvesting behavior of plateau pikas with the plant community[J].ActaTheriologicaSinica, 29(1):40-49.(in Chinese)

    [11] 劉偉, 李里, 嚴(yán)紅宇,等. 高原鼠兔挖掘活動(dòng)對(duì)植物種的豐富度和地上生物量的影響[J]. 獸類(lèi)學(xué)報(bào), 2012,32(3):216-220.

    LIU Wei, LI Li, YAN Hong-yu, et al.(2012).Effects of burrowing activity of plateau pika (Ochotonacurzoniae) on plant species richness and aboveground biomass[J].ActaTheriologicaSinica,32(3):216-220.(in Chinese)

    [12] 周雪榮, 郭正剛, 郭興華. 高原鼠兔和高原鼢鼠在高寒草甸中的作用[J]. 草業(yè)科學(xué), 2010,27(5):38-44.

    ZHOU Xue-rong, GUO Zheng-gang, GUO Xing-hua, et al.(2010).The role of plateau pika and plateall zokor in alpine meadow [J].PrataculturalScience, 27(5):38-44.(in Chinese)

    [13] 劉偉, 張毓, 王溪,等. 高原鼠兔冬季的食物選擇[J]. 獸類(lèi)學(xué)報(bào), 2009,29(1):12-19.

    LIU Wei, ZHANG Yu, WANG Xi, et al. (2009). Food selection pattern for plateau pika in winter [J].ActaTheriologicaSinica, 29(1):12-19.(in Chinese)

    [14] 楊虎虎, 曲家鵬, 劉力華,等. 食物及洞道密度對(duì)高原鼠兔冬季死亡率的影響[J].獸類(lèi)學(xué)報(bào),2011,31(3):257-264.

    YANG Hu-hu, QU Jia-peng,LIU Li-hua, et al.(2011).Effects of food supply and burrow density on mortality of the plateau pika (Ochotonacurzoniae) in winter [J].ActaTheriologicaSinica, 31(3):257-264.(in Chinese)

    [15] 張衛(wèi)國(guó), 劉蓉, 江小雷. 風(fēng)險(xiǎn)性聲訊信號(hào)對(duì)高原鼠兔行為模式的影響[J]. 草地學(xué)報(bào), 2010,18(1):115-120.

    ZHANG Wei-guo,LIU Rong,JIANG Xiao-lei. (2010). Influence of risk sound signal on behavior pattern of Pika [J].ActaAgrestiaSinica, 18(1):115-120.(in Chinese)

    [16] 李葉, 王振宇, 張翔,等. 阿爾金山自然保護(hù)區(qū)高原鼠兔夏季微生境選擇的主導(dǎo)因子分析[J]. 中國(guó)媒介生物學(xué)及控制雜志, 2014,25(1):28-31.

    Li Ye, WANG Zhen-yu, Zhang Xiang, et al. (2014). Analysis of dominant factors affecting microhabitat selection of plateau pika (Ochotonacurzoniae) during summer in Altun Mountain National Nature Reserve, Xinjiang Uygur Autonomous Region, China [J].ChineseJournalofVectorBiologyandControl, 25(1):28-31. (in Chinese)

    [17] 熱娜古麗·艾合麥提, 李葉, 張翔,等. 阿爾金山國(guó)家級(jí)自然保護(hù)區(qū)高原鼠兔的食性顯微組織分析[J]. 中國(guó)媒介生物學(xué)及控制雜志,2015,(2):164-167.

    Renagül Exmet, LI Ye, ZHANG Xiang, et al. (2015). Microscopic histological analysis of diet ofOchotonacurzoniaein the Altun Mountain National Nature Reserve [J].ChineseJournalofVectorBiologyandControl, 26(2):164-167.(in Chinese)

    [18] 賈婷婷, 毛亮, 郭正剛. 高原鼠兔有效洞穴密度對(duì)青藏高原高寒草甸群落植物生態(tài)位的影響[J]. 生態(tài)學(xué)報(bào), 2014, 34(4):869-877.

    JIA Ting-ting, MAO Liang, GUO Zheng-gang.(2014).Effect of available burrow densities of plateau pika (Ochotonacurzoniae) on plant niche of alpine meadow communities in the Qinghai-Tibet Plateau [J].ActaEcologicaSinica, 34(4):869-877.(in Chinese)

    [19]馬波, 王小明, 劉曉慶,等. 高原鼠兔洞穴數(shù)量與其棲息地植被分布格局的GIS分析[J]. 生物多樣性, 2011,19(1):71-78.

    MA Bo, WANG Xiao-ming, LIU Xiao-qin, et al.(2011).GIS analysis of the spatial relationship between plateau pika burrow distribution and vegetation distributional patterns [J].BiodiversityScience, 19(1):71-78.(in Chinese)

    [20] 劉世梁, 趙海迪, 董世魁,等. 基于SPOT NDVI的阿爾金山自然保護(hù)區(qū)植被動(dòng)態(tài)變化研究[J]. 干旱區(qū)研究, 2014,31(5):832-837.

    LIU Shi-liang, ZHAO Hai-di, DONG Shi-kui, et al.(2014).Dynamic of vegetation in the Altun Mountain Nature Reserve based on SPOT NDVI [J].AridZoneResearch, 31(5):832-837.(in Chinese)

    [21]許鵬. 新疆草地資源及其利用[M].烏魯木齊:新疆科技衛(wèi)生出版社,1993:349-357.

    XU Peng.(1993).GrasslandResourcesandUtilizationinXinjiang[M].Urumqi: Health Science and Technology Publishing House in Xinjiang:349-357.(in Chinese)

    [22]徐瑤. 藏北草地退化遙感監(jiān)測(cè)與生態(tài)安全評(píng)價(jià)[D]. 成都:成都理工大學(xué)博士論文,2014.

    XU Yao.(2014).RemotesensingdynamicmonitoringandEco-securityAssessmentofGrasslandDegradationinNorthernTibet[D]. PhD Thesis. Chengdu University of Technology,Chengdu. (in Chinese)

    [23] Richards, John A. (1986). Remote Sensing Digital Image Analysis: An Introduction. Remote sensing digital image analysis:Springer-Verlag.

    [24] Islam Abou EL-Magd Corresponding author, & T. W. Tanton. (2003). Improvements in land use mapping for irrigated agriculture from satellite sensor data using a multi-stage maximum likelihood classification.InternationalJournalofRemoteSensing, 24(21):4,197-4,206.

    [25]陳亮, 劉希, 張?jiān)? 結(jié)合光譜角的最大似然法遙感影像分類(lèi)[J]. 測(cè)繪工程, 2007,16(3):40-42.

    CHEN Liang, LIU Xi, ZHANG Yuan. (2007). A study of image classification based on MLC combined with spectral angle [J].EngineeringofSurveyingandMapping, 16(3):40-42. (in Chinese)

    [26]李德仁, 王樹(shù)良, 李德毅,等.論空間數(shù)據(jù)挖掘和知識(shí)發(fā)現(xiàn)的理論與方法[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2002,27(3):221-233.

    LI De-ren, WANG Shu-liang, LI De-yi, et al.(2002).Theories and Technologies of Spatial Data Mining and Knowledge Discovery [J].GeomaticsandInformationScienceofWuhanUniversity, 27(3):221-233. (in Chinese)

    [27] Hansen, M., Dubayah, R., & Defries, R. (2007). Classification trees: an alternative to traditional land cover classifiers.InternationalJournalofRemoteSensing, 17(5):1,075-1,081.

    Fund project:Supported by the project commissioned by Locust Rodent Headquarters Office of Xinjiang Uygur Autonomous Region "Remote sensing monitoring of grassland biological disasters in Xinjiang" (2014-2015)and High level talent plan program of Xinjiang Uygur Autonomous Region (104-40002)

    Estimating Area of Grassland Rodent Damage Rangeland and Rat Wastelands Based on Remote Sensing in Altun Mountain, Xinjiang, China

    LI Pei-xian1,2, ZHENG Jiang-hua1,2,3, NI Yi-fei4, WU Jian-guo4, Wumaier Wushou5,Aihemaijiang Aji5,Nasongcaoketu5

    (1. College of Resources & Environmental Sciences,Xinjiang University, Urumqi 830046, China;2.KeyLaboratoryofCityIntellectualizingandEnvironmentModelling,XinjiangUniversity,Urumqi830046,China; 3.KeyLaboratoryofOasisEcology,Urumqi830046,China; 4.LocustandRodentControlHeadquartersofXinjiangUygurAutonomousRegion,Urumqi830004,China;5.GrasslandWorkstationofBayingolinMongolAutonomousPrefecture,KorlaXinjiang841000,China)

    【Objective】 Monitoring and interpreting grassland rodent damage rangeland and rat wastelands in Altun Mountain with 3S technologies can acquire area and distribution of grassland rodent in order to provide a reference in grassland rodents control and management.【Method】Landsat 8 remote sensing image of study area were classified preliminarily into pest occurrence area, grassland, bare land, body of water and others by supervising classification of maximum likelihood classification. Then rat wastelands, rodent damage rangeland and rodent damage potential (local) occurrence area were extracted from preliminary pest occurrence area by decision tree classifier based on rules about altitude and vegetation coverage further. Finally, the thematic map of grassland rodent damage rangeland and rat wastelands was made and occurrence area was calculated.【Result】The area of grassland rodent damage rangeland and rat wastelands in Altun Mountain were 3,419.621 km2and 1,376.476 km2, accounting for 6.92% and 2.78% of the study area, which is similar with the estimated result of Bazhou Prairie Station in November, 2014. And the area of grassland rodent damage rangeland has the trend of expansion.【Conclusion】It is a characteristic that rat wastelands, rodent damage rangeland and rodent damage potential(local) occurrence area are along with the edge of the water, and the rat wasteland, the bare land and the grassland are mostly of mosaic distribution. The coexistence of the wasteland and bare land, grassland is easy to lead to a mixed pixel interpretation. The combination of the maximum likelihood supervised classification and the decision tree method to extract rat wastelands, rodent damage rangeland and rodent damage potential (local) occurrence area makes the results more consistent with the actual survey

    "3S" technologies; Altun Mountain; supervised classification; decision tree; rat wasteland

    10.6048/j.issn.1001-4330.2016.07.023

    2016-03-20

    新疆維吾爾自治區(qū)治蝗滅鼠指揮部辦公室委托項(xiàng)目“新疆草原生物災(zāi)害遙感監(jiān)測(cè)”(2014~2015);新疆維吾爾自治區(qū)高層次人才計(jì)劃(104-40002)。

    李培先(1991-),男,新疆伊寧人,碩士研究生,研究方向?yàn)椴菰餅?zāi)害遙感監(jiān)測(cè),(E-mail)15894623722@163.com

    鄭江華(1973-),男,浙江江山人,教授,研究方向?yàn)檫b感與地理信息系統(tǒng)應(yīng)用,(E-mail)zheng_jianghua@126.com

    S443

    A

    1001-4330(2016)07-1346-10

    猜你喜歡
    阿爾金山鼠害荒地
    南望阿爾金山
    都市(2023年5期)2023-09-13 02:17:19
    獨(dú)登南山
    草原鼠害的危害及綜合防治措施
    阿爾金山西部區(qū)域化探數(shù)據(jù)處理方法對(duì)比研究
    荒 地
    春到阿爾金山
    中亞信息(2018年4期)2018-05-28 07:44:25
    玉樹(shù)草原鼠害防治建議
    玉樹(shù)地區(qū)草原鼠害防治的建設(shè)與管理
    對(duì)外發(fā)包的荒地為何被判無(wú)效
    巴里坤縣草原鼠害防治現(xiàn)狀與對(duì)策
    99热精品在线国产| 黄色配什么色好看| 在现免费观看毛片| 国产一区二区亚洲精品在线观看| 麻豆久久精品国产亚洲av| 天堂网av新在线| 国产亚洲av嫩草精品影院| 亚洲国产精品成人久久小说 | 18禁黄网站禁片免费观看直播| 一个人看视频在线观看www免费| 久久久国产成人免费| 亚洲内射少妇av| 两个人的视频大全免费| 淫妇啪啪啪对白视频| 国产片特级美女逼逼视频| 永久网站在线| 精品人妻熟女av久视频| а√天堂www在线а√下载| 国产亚洲精品久久久久久毛片| 好男人在线观看高清免费视频| 男人舔女人下体高潮全视频| 国产精品电影一区二区三区| 91久久精品国产一区二区成人| 国产乱人视频| 久久久久久久久久久丰满| 中文字幕熟女人妻在线| 亚洲最大成人中文| 午夜免费激情av| 在线观看一区二区三区| 婷婷亚洲欧美| 别揉我奶头~嗯~啊~动态视频| 久久久成人免费电影| 亚洲精华国产精华液的使用体验 | 国产高清视频在线播放一区| 欧美在线一区亚洲| 日韩成人伦理影院| 最近手机中文字幕大全| 国产高清激情床上av| 日日摸夜夜添夜夜添av毛片| 久久久久九九精品影院| 97人妻精品一区二区三区麻豆| 亚洲欧美日韩高清在线视频| 国产精华一区二区三区| 哪里可以看免费的av片| 黄色一级大片看看| 色视频www国产| 日本精品一区二区三区蜜桃| 亚洲电影在线观看av| 麻豆精品久久久久久蜜桃| 色哟哟·www| 国产一区二区三区av在线 | 有码 亚洲区| 人妻夜夜爽99麻豆av| 午夜视频国产福利| 男女之事视频高清在线观看| av在线观看视频网站免费| 精品一区二区三区av网在线观看| 免费观看人在逋| a级一级毛片免费在线观看| 国产精品av视频在线免费观看| 桃色一区二区三区在线观看| 久久午夜福利片| 黄色视频,在线免费观看| 久久久久久伊人网av| 国产午夜精品久久久久久一区二区三区 | 综合色av麻豆| 男女下面进入的视频免费午夜| 看十八女毛片水多多多| 成年av动漫网址| 高清毛片免费观看视频网站| 国产精品国产三级国产av玫瑰| 国产视频一区二区在线看| 老师上课跳d突然被开到最大视频| 亚洲av免费高清在线观看| 久久热精品热| 国产高清三级在线| 午夜激情福利司机影院| 露出奶头的视频| 亚洲国产精品sss在线观看| 久久久久久久久中文| 久久精品国产亚洲网站| 我要看日韩黄色一级片| 悠悠久久av| 成人综合一区亚洲| 极品教师在线视频| 亚洲婷婷狠狠爱综合网| 搡老妇女老女人老熟妇| 嫩草影院新地址| 嫩草影视91久久| 国产在线男女| 在线看三级毛片| 深爱激情五月婷婷| 男人舔奶头视频| 特级一级黄色大片| 亚洲精品456在线播放app| 深爱激情五月婷婷| av国产免费在线观看| h日本视频在线播放| 韩国av在线不卡| 91在线精品国自产拍蜜月| 搡老岳熟女国产| 国语自产精品视频在线第100页| 久久人人爽人人片av| 久久这里只有精品中国| 国产91av在线免费观看| 久久韩国三级中文字幕| 欧美不卡视频在线免费观看| 欧美日韩在线观看h| 午夜老司机福利剧场| 久久这里只有精品中国| 在线观看免费视频日本深夜| 别揉我奶头 嗯啊视频| 中国美女看黄片| 午夜福利18| 亚洲色图av天堂| 99热这里只有是精品50| 亚洲欧美日韩高清专用| 国产精品女同一区二区软件| 2021天堂中文幕一二区在线观| 白带黄色成豆腐渣| 精品人妻熟女av久视频| 99在线人妻在线中文字幕| 大香蕉久久网| 一进一出抽搐动态| 女生性感内裤真人,穿戴方法视频| avwww免费| 插逼视频在线观看| 亚洲成人精品中文字幕电影| 在线a可以看的网站| 伊人久久精品亚洲午夜| 午夜日韩欧美国产| 国产色爽女视频免费观看| 极品教师在线视频| 欧美3d第一页| 51国产日韩欧美| 久久国内精品自在自线图片| 精品久久久久久久末码| 一个人免费在线观看电影| 有码 亚洲区| 天堂动漫精品| 黄片wwwwww| 国产精品久久久久久亚洲av鲁大| 国内精品美女久久久久久| 十八禁国产超污无遮挡网站| 99热全是精品| 欧美性猛交╳xxx乱大交人| 97在线视频观看| 99在线人妻在线中文字幕| 日韩三级伦理在线观看| 精品一区二区免费观看| 精品国产三级普通话版| 亚洲自拍偷在线| 国产三级中文精品| 一进一出抽搐gif免费好疼| 国产午夜福利久久久久久| 在线观看免费视频日本深夜| 网址你懂的国产日韩在线| 中文字幕av在线有码专区| 级片在线观看| 日本欧美国产在线视频| 午夜老司机福利剧场| 欧美成人a在线观看| 日本a在线网址| 免费大片18禁| 国产精品女同一区二区软件| 国产成年人精品一区二区| avwww免费| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕av成人在线电影| 九九在线视频观看精品| 不卡视频在线观看欧美| 熟妇人妻久久中文字幕3abv| 中文亚洲av片在线观看爽| 欧美bdsm另类| 色吧在线观看| 精品一区二区免费观看| 日日撸夜夜添| 国产成人福利小说| 国产一区二区激情短视频| 男女边吃奶边做爰视频| 国产av不卡久久| 国产一区亚洲一区在线观看| 久久中文看片网| 丰满人妻一区二区三区视频av| 国产精品国产三级国产av玫瑰| 免费av毛片视频| 99热只有精品国产| 国产精品,欧美在线| 国内揄拍国产精品人妻在线| eeuss影院久久| 成人国产麻豆网| 久久久成人免费电影| 国产成人91sexporn| 欧美性猛交黑人性爽| 欧美性感艳星| 中国美女看黄片| 简卡轻食公司| 成年女人毛片免费观看观看9| 麻豆乱淫一区二区| 美女高潮的动态| 欧美3d第一页| 国产精品一二三区在线看| 麻豆成人午夜福利视频| 精品欧美国产一区二区三| 久久久久九九精品影院| 亚洲精品久久国产高清桃花| 久久久久国内视频| 99riav亚洲国产免费| 麻豆乱淫一区二区| 丝袜喷水一区| 欧美激情久久久久久爽电影| 亚洲乱码一区二区免费版| 露出奶头的视频| 在现免费观看毛片| 夜夜夜夜夜久久久久| 色av中文字幕| 久久久a久久爽久久v久久| 亚洲成av人片在线播放无| 成人高潮视频无遮挡免费网站| 一级黄色大片毛片| 一a级毛片在线观看| 亚洲精品国产成人久久av| 男插女下体视频免费在线播放| 色视频www国产| 免费黄网站久久成人精品| 精品国内亚洲2022精品成人| 看十八女毛片水多多多| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频 | 亚洲精品成人久久久久久| 嫩草影院新地址| 男人的好看免费观看在线视频| 干丝袜人妻中文字幕| 免费人成在线观看视频色| 一区二区三区高清视频在线| 免费电影在线观看免费观看| 午夜福利在线观看吧| 日本在线视频免费播放| 国产精品人妻久久久久久| 国产成人精品久久久久久| 九色成人免费人妻av| 51国产日韩欧美| 国产成人a区在线观看| 亚洲欧美成人精品一区二区| 亚洲经典国产精华液单| 在线天堂最新版资源| 国产高潮美女av| 亚洲精品色激情综合| 日本一二三区视频观看| 国产精品永久免费网站| 国产精品一区二区免费欧美| 国产色婷婷99| 日韩在线高清观看一区二区三区| 人人妻人人澡欧美一区二区| 亚洲七黄色美女视频| 成熟少妇高潮喷水视频| 人妻制服诱惑在线中文字幕| 97人妻精品一区二区三区麻豆| 狠狠狠狠99中文字幕| 非洲黑人性xxxx精品又粗又长| 成人综合一区亚洲| 亚洲欧美精品自产自拍| 九色成人免费人妻av| 一个人观看的视频www高清免费观看| 国产单亲对白刺激| 久久精品国产亚洲网站| 一卡2卡三卡四卡精品乱码亚洲| 99精品在免费线老司机午夜| 99久久九九国产精品国产免费| 99久久久亚洲精品蜜臀av| 亚洲综合色惰| 18禁在线播放成人免费| 午夜免费激情av| 免费看av在线观看网站| 一级毛片aaaaaa免费看小| 亚洲国产精品合色在线| 亚洲成人久久性| 久久这里只有精品中国| 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 少妇熟女欧美另类| 99热这里只有是精品50| av天堂中文字幕网| 丰满的人妻完整版| 高清毛片免费观看视频网站| 久久精品夜色国产| 久久久久精品国产欧美久久久| 日日摸夜夜添夜夜添av毛片| 中文亚洲av片在线观看爽| 一a级毛片在线观看| 成人漫画全彩无遮挡| 给我免费播放毛片高清在线观看| 亚洲一区高清亚洲精品| av专区在线播放| 国产老妇女一区| 亚洲成人中文字幕在线播放| 久久精品综合一区二区三区| 日韩高清综合在线| 成年免费大片在线观看| 成人美女网站在线观看视频| 一个人看的www免费观看视频| 日本免费一区二区三区高清不卡| 久久久久久久久大av| 欧美日韩乱码在线| 狂野欧美白嫩少妇大欣赏| 最近最新中文字幕大全电影3| 久久久久国产精品人妻aⅴ院| 天天躁夜夜躁狠狠久久av| avwww免费| 乱人视频在线观看| 亚洲av免费高清在线观看| 老司机福利观看| 18+在线观看网站| 国内精品美女久久久久久| 成人国产麻豆网| 免费大片18禁| 国产成人aa在线观看| 国产aⅴ精品一区二区三区波| 亚洲av五月六月丁香网| 久久人人精品亚洲av| 国产精品三级大全| 日韩高清综合在线| 欧美成人精品欧美一级黄| 99热精品在线国产| 熟女电影av网| 亚洲欧美成人综合另类久久久 | 老司机午夜福利在线观看视频| 18+在线观看网站| 国产高清视频在线播放一区| 色综合亚洲欧美另类图片| 男女之事视频高清在线观看| 欧美xxxx性猛交bbbb| 国产精品电影一区二区三区| 中文亚洲av片在线观看爽| 国内揄拍国产精品人妻在线| 超碰av人人做人人爽久久| 三级毛片av免费| av专区在线播放| 色哟哟哟哟哟哟| 日韩欧美三级三区| 亚洲aⅴ乱码一区二区在线播放| 麻豆乱淫一区二区| 麻豆一二三区av精品| 你懂的网址亚洲精品在线观看 | 国产亚洲欧美98| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区色噜噜| 国产精品,欧美在线| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线播| 人妻制服诱惑在线中文字幕| 午夜福利18| 偷拍熟女少妇极品色| 麻豆乱淫一区二区| 真人做人爱边吃奶动态| 我的老师免费观看完整版| 丰满乱子伦码专区| 99久国产av精品| 网址你懂的国产日韩在线| 日韩欧美一区二区三区在线观看| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 欧美日韩综合久久久久久| 床上黄色一级片| 亚洲av电影不卡..在线观看| 亚洲电影在线观看av| 99九九线精品视频在线观看视频| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| 综合色av麻豆| 中文字幕免费在线视频6| 在线天堂最新版资源| 国产精品爽爽va在线观看网站| 成人鲁丝片一二三区免费| 亚洲专区国产一区二区| 国产av不卡久久| 特大巨黑吊av在线直播| 不卡一级毛片| 18禁在线播放成人免费| 亚洲第一电影网av| 永久网站在线| 中文字幕熟女人妻在线| 久久久久精品国产欧美久久久| 悠悠久久av| 我要搜黄色片| 在线播放国产精品三级| 亚洲国产高清在线一区二区三| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 99riav亚洲国产免费| 黄色欧美视频在线观看| 亚洲第一区二区三区不卡| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 国产毛片a区久久久久| 国产精品国产三级国产av玫瑰| 亚洲最大成人手机在线| 看黄色毛片网站| 少妇人妻一区二区三区视频| 中国美女看黄片| 亚洲在线自拍视频| 波多野结衣高清无吗| 观看美女的网站| 久久久久久久久久久丰满| 亚洲人成网站在线播放欧美日韩| 国产男人的电影天堂91| 久久久国产成人精品二区| 国产亚洲精品久久久久久毛片| 亚洲无线在线观看| 美女高潮的动态| 伊人久久精品亚洲午夜| 亚洲内射少妇av| 亚洲国产精品国产精品| 午夜福利18| 精品久久久久久久末码| 日韩大尺度精品在线看网址| 亚洲精品一区av在线观看| av在线播放精品| 高清毛片免费看| 国产高清视频在线播放一区| 18禁在线播放成人免费| 又爽又黄a免费视频| 天天躁夜夜躁狠狠久久av| 岛国在线免费视频观看| 蜜桃久久精品国产亚洲av| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区国产精品久久精品| 欧美色欧美亚洲另类二区| 三级男女做爰猛烈吃奶摸视频| 亚洲自偷自拍三级| 最新在线观看一区二区三区| 成年版毛片免费区| 久久精品影院6| 成人亚洲精品av一区二区| 成人特级黄色片久久久久久久| 国产男人的电影天堂91| 韩国av在线不卡| 看十八女毛片水多多多| 成年女人永久免费观看视频| 日本在线视频免费播放| 亚洲一区二区三区色噜噜| 亚洲最大成人av| 人妻久久中文字幕网| 成人av一区二区三区在线看| 搡老妇女老女人老熟妇| 精品久久久久久久久久免费视频| 国产av在哪里看| 俺也久久电影网| 婷婷亚洲欧美| 偷拍熟女少妇极品色| 在线国产一区二区在线| 天堂√8在线中文| 久久精品国产鲁丝片午夜精品| 1000部很黄的大片| 久久天躁狠狠躁夜夜2o2o| 久久精品国产亚洲av涩爱 | 亚洲高清免费不卡视频| 性欧美人与动物交配| 国产精品综合久久久久久久免费| 97热精品久久久久久| 国产伦精品一区二区三区四那| 黑人高潮一二区| 成年av动漫网址| 日本色播在线视频| 最新在线观看一区二区三区| 国产高潮美女av| 久久午夜亚洲精品久久| 精品不卡国产一区二区三区| av在线播放精品| 人人妻人人看人人澡| 国产男人的电影天堂91| 欧美成人免费av一区二区三区| 日日啪夜夜撸| 亚洲av免费高清在线观看| 亚洲经典国产精华液单| 久久精品综合一区二区三区| 国产乱人偷精品视频| av天堂在线播放| 一本久久中文字幕| 色5月婷婷丁香| 日本精品一区二区三区蜜桃| 久久久国产成人精品二区| 日本黄色视频三级网站网址| 欧美一区二区亚洲| 老女人水多毛片| 日日啪夜夜撸| 久久久久九九精品影院| 中文资源天堂在线| 91麻豆精品激情在线观看国产| 国产中年淑女户外野战色| 国产熟女欧美一区二区| 国产一区二区三区在线臀色熟女| 亚洲国产精品国产精品| av中文乱码字幕在线| 一本久久中文字幕| 国产伦在线观看视频一区| 免费不卡的大黄色大毛片视频在线观看 | 22中文网久久字幕| 97热精品久久久久久| 欧美人与善性xxx| 黄色视频,在线免费观看| 男女那种视频在线观看| 看非洲黑人一级黄片| 最近2019中文字幕mv第一页| 一个人看的www免费观看视频| 淫秽高清视频在线观看| 精品久久久久久久久久免费视频| 日韩中字成人| 少妇熟女欧美另类| 亚洲欧美精品自产自拍| 成人无遮挡网站| 日韩国内少妇激情av| 欧美一区二区国产精品久久精品| 久久久精品大字幕| avwww免费| 特大巨黑吊av在线直播| 97人妻精品一区二区三区麻豆| 桃色一区二区三区在线观看| 亚洲真实伦在线观看| 好男人在线观看高清免费视频| 日本黄色视频三级网站网址| 亚洲av熟女| 欧美高清成人免费视频www| 变态另类丝袜制服| 一级黄片播放器| 久久精品夜夜夜夜夜久久蜜豆| 最近在线观看免费完整版| 网址你懂的国产日韩在线| 最近视频中文字幕2019在线8| 国产成人a区在线观看| 日本 av在线| 免费看av在线观看网站| 成人毛片a级毛片在线播放| 深爱激情五月婷婷| 国产色爽女视频免费观看| 看非洲黑人一级黄片| 日本撒尿小便嘘嘘汇集6| 不卡一级毛片| 美女内射精品一级片tv| 日本熟妇午夜| 中文字幕熟女人妻在线| 青春草视频在线免费观看| 国产精品久久久久久精品电影| 老司机午夜福利在线观看视频| 中文字幕av成人在线电影| 99久国产av精品国产电影| 日韩精品中文字幕看吧| av在线天堂中文字幕| 一级毛片aaaaaa免费看小| 国产成年人精品一区二区| 精品久久久噜噜| 精品午夜福利在线看| 老女人水多毛片| 丰满的人妻完整版| 亚洲人成网站在线播| 婷婷色综合大香蕉| 久久99热6这里只有精品| 欧美+亚洲+日韩+国产| 国产麻豆成人av免费视频| 99久久久亚洲精品蜜臀av| 亚洲欧美清纯卡通| 在线观看免费视频日本深夜| 亚洲婷婷狠狠爱综合网| 看十八女毛片水多多多| av视频在线观看入口| 在线播放无遮挡| 深爱激情五月婷婷| 国产午夜精品论理片| 观看美女的网站| 日韩精品中文字幕看吧| 少妇猛男粗大的猛烈进出视频 | 亚洲性夜色夜夜综合| 国产av麻豆久久久久久久| 国产精品人妻久久久影院| 日韩大尺度精品在线看网址| 国国产精品蜜臀av免费| 久久久欧美国产精品| 免费看美女性在线毛片视频| 欧美人与善性xxx| 女人十人毛片免费观看3o分钟| 欧美不卡视频在线免费观看| 在线观看午夜福利视频| 日本-黄色视频高清免费观看| 久久午夜亚洲精品久久| 欧美三级亚洲精品| 最近在线观看免费完整版| 国产乱人视频| 秋霞在线观看毛片| 99国产极品粉嫩在线观看| 亚洲人与动物交配视频| 精品人妻偷拍中文字幕| 国产精品美女特级片免费视频播放器| 一区二区三区高清视频在线| 国内精品宾馆在线| 欧美日韩一区二区视频在线观看视频在线 | 女生性感内裤真人,穿戴方法视频| 久久人妻av系列| 国产一区亚洲一区在线观看| 国产成人a区在线观看| 亚洲精品色激情综合| 亚洲人成网站高清观看| 久久久久久久久久黄片| 校园春色视频在线观看| 一区福利在线观看| 欧美国产日韩亚洲一区| 国产精品久久久久久精品电影| 久久久久久久久久久丰满| 久久精品国产亚洲av香蕉五月| 国内久久婷婷六月综合欲色啪| 国产蜜桃级精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 国产精品一及| 国内精品久久久久精免费| 色哟哟·www| 人人妻人人看人人澡|