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

    基于MODIS時(shí)序數(shù)據(jù)的黑龍江流域火燒跡地提取

    2015-01-18 08:48:29張樹(shù)文姜曉麗
    生態(tài)學(xué)報(bào) 2015年17期
    關(guān)鍵詞:跡地火燒黑龍江

    楊 偉, 張樹(shù)文, 姜曉麗

    1 太原師范學(xué)院地理科學(xué)學(xué)院, 晉中 030619 2 中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所, 長(zhǎng)春 130102 3 太原師范學(xué)院城鎮(zhèn)與區(qū)域發(fā)展研究所, 晉中 030619

    基于MODIS時(shí)序數(shù)據(jù)的黑龍江流域火燒跡地提取

    楊 偉1, 張樹(shù)文2,*, 姜曉麗3

    1 太原師范學(xué)院地理科學(xué)學(xué)院, 晉中 030619 2 中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所, 長(zhǎng)春 130102 3 太原師范學(xué)院城鎮(zhèn)與區(qū)域發(fā)展研究所, 晉中 030619

    火燒跡地信息是研究火災(zāi)的重要參數(shù)和基礎(chǔ)數(shù)據(jù),也是研究全球生態(tài)系統(tǒng)和碳循環(huán)擾動(dòng)的重要依據(jù)之一。以受森林火災(zāi)影響較為嚴(yán)重的黑龍江流域?yàn)檠芯繀^(qū),以MODIS時(shí)間序列數(shù)據(jù)為數(shù)據(jù)源建立了一個(gè)分為兩階段的火燒基地提取算法(即首先設(shè)定較為嚴(yán)格的提取條件對(duì)最有可能發(fā)生火災(zāi)的像元——核心像元進(jìn)行提取,然后設(shè)定較為寬松的閾值提取距離核心像元一定范圍內(nèi)的火燒像元),對(duì)2000—2011年的火燒跡地信息進(jìn)行了提取,生成了研究區(qū)長(zhǎng)時(shí)間序列火燒跡地分布圖,并對(duì)其時(shí)空分布特征進(jìn)行了分析。選擇黑龍江省為典型驗(yàn)證區(qū)對(duì)算法精度進(jìn)行了驗(yàn)證,結(jié)果顯示算法的整體精度較之以往的算法有了一定程度的提高。

    火燒跡地; MODIS; GEMI; BAI

    火災(zāi)是影響眾多生態(tài)系統(tǒng)(森林、草地等)的一個(gè)重要的擾動(dòng)因素?;馂?zāi)對(duì)于植被的結(jié)構(gòu)和組成具有顯著的作用,被認(rèn)為是一個(gè)重要的“地表管理工具”[1]。以森林生態(tài)系統(tǒng)為例,火災(zāi)是森林生態(tài)系統(tǒng)最為重要的干擾之一,全球平均每年約有1%的森林受到火災(zāi)的影響[2]。森林生態(tài)系統(tǒng)是全球碳循環(huán)的重要組成部分,火災(zāi)的發(fā)生通過(guò)改變森林生態(tài)系統(tǒng)的格局與過(guò)程,進(jìn)而改變整個(gè)生態(tài)系統(tǒng)的碳循環(huán)以及分配過(guò)程[3]。此外,森林火災(zāi)對(duì)于氣候變暖也有著重要的響應(yīng)[4- 6]。研究表明,北方森林的火災(zāi)發(fā)生范圍對(duì)于溫度的增加非常敏感,受氣候變化的影響較為顯著[7]。此外,還有研究認(rèn)為氣候變暖很有可能影響到傳統(tǒng)的火災(zāi)循環(huán),例如縮短火災(zāi)的周期、增大火災(zāi)尺度等[8-9]。

    不論是對(duì)于碳循環(huán)的影響研究,還是森林火災(zāi)與氣候變暖的相關(guān)性研究,火燒跡地的空間信息都是一個(gè)重要的基礎(chǔ)參數(shù)。傳統(tǒng)的火燒跡地信息主要來(lái)源于統(tǒng)計(jì)數(shù)據(jù),難以覆蓋較大的區(qū)域,收集較為困難,且難以將數(shù)據(jù)進(jìn)行空間化。遙感技術(shù)的發(fā)展為解決這一問(wèn)題提供了很好的手段,特別是隨著遙感數(shù)據(jù)時(shí)空分辨率的提高,使得遙感數(shù)據(jù)能夠更為準(zhǔn)確的對(duì)地表過(guò)程進(jìn)行刻畫[10]。

    在區(qū)域或者全球尺度下,為獲取長(zhǎng)時(shí)間序列的火燒跡地信息,中空間分辨率且具有高時(shí)間分辨率的遙感數(shù)據(jù)被認(rèn)為是最好的選擇。目前,應(yīng)用最為廣泛的為AVHRR(Advanced Very High Resolution Radiometer)數(shù)據(jù)[11- 13]與MODIS(Moderate-Resolution Imaging Spectroradia-meter) 數(shù)據(jù)[14-15]。前者由于發(fā)射時(shí)間較早,時(shí)間序列較長(zhǎng)而被使用。但研究表明AVHRR數(shù)據(jù)提取火燒跡地信息存在一定的潛在誤差,主要來(lái)源于輻射的不穩(wěn)定性、云污染以及輻射傳輸問(wèn)題等方面[16- 18]。較之AVHRR數(shù)據(jù),MODIS數(shù)據(jù)在這些方面都得到了很大的改善,但其局限在于數(shù)據(jù)僅從2000年開(kāi)始。

    本文提出了一種基于MODIS時(shí)序數(shù)據(jù)的火燒跡地提取方法,以此為基礎(chǔ)對(duì)黑龍江流域2000—2011年的火燒跡地信息進(jìn)行了提取,生成研究區(qū)長(zhǎng)時(shí)間序列火燒跡地分布圖,并對(duì)其分布特征進(jìn)行了分析。

    1 材料與方法

    1.1 研究區(qū)

    研究區(qū)選擇位于西伯利亞北方森林南部的黑龍江流域(41°45′—53°33′N,115°13′—135°05′E),面積208×104km2,西起蒙古高原,包括蒙古(東方省、蘇赫巴托爾省等)、俄羅斯(阿穆?tīng)栔?、哈巴羅夫斯克、外貝加爾等)和中國(guó)(黑龍江省、吉林省、遼寧省、內(nèi)蒙古自治區(qū)等)的13個(gè)省及朝鮮的小部分,研究區(qū)植被覆蓋度較高(圖1),受森林火災(zāi)影響較為嚴(yán)重。流域的東部地區(qū)主要屬于溫帶濕潤(rùn)季風(fēng)氣候,這是全球季風(fēng)氣候的最北緣,西部主要受大陸性氣候的影響。全年平均氣溫在-8 ℃到6 ℃之間,但其時(shí)空分布差異顯著。同時(shí),流域內(nèi)降水量的時(shí)空分布也很不均衡,年平均降水量主要在250—800 mm,大約50%以上的降水量集中在最熱的夏季,而近7個(gè)月的干季(1—4月,10—12月)降水量?jī)H為25%;在空間上,降水主要集中在沿海地帶,向西逐漸遞減。

    圖1 研究區(qū)位置及其土地覆被Fig.1 Location and land cover of Heilongjiang basin

    1.2 MODIS產(chǎn)品數(shù)據(jù)

    研究選擇數(shù)據(jù)為MODIS 8d合成地表反射率數(shù)據(jù)MOD09Q1(MODIS atmospherically- correct Level 3 8-Day composite Surface Reflectance products, 空間分辨率250 m)以及MODIS 8d合成火產(chǎn)品數(shù)據(jù)MOD14A2(MODIS Level 3 8-Day composite active fire products, 空間分辨率1 km)。其中,MODIS火產(chǎn)品數(shù)據(jù)為MODIS火情監(jiān)測(cè)算法下提取的火點(diǎn)信息[19-20],其算法依據(jù)主要為火災(zāi)發(fā)生時(shí)的熱學(xué)特性,代表火災(zāi)發(fā)生時(shí)的溫度異常。

    MOD09Q1與MOD14A2產(chǎn)品均采用Sinusoidal投影系統(tǒng)發(fā)布,數(shù)據(jù)格式為HDF (Hierarchy Data Format)。研究區(qū)共涉及6景MODIS標(biāo)準(zhǔn)分幅數(shù)據(jù)。下載研究區(qū)2000—2011年的MODIS產(chǎn)品數(shù)據(jù)(每年34期×12年×6景)。時(shí)間選擇為每年的2月末至11月初,這一時(shí)間為火災(zāi)發(fā)生集中的時(shí)間段。對(duì)數(shù)據(jù)進(jìn)行投影轉(zhuǎn)換及裁剪處理。此外,由于數(shù)據(jù)空間分辨率不一致,需要將MOD14A2產(chǎn)品重采樣為250 m分辨率,與MOD09Q1產(chǎn)品相一致。

    1.3 火燒跡地提取方法

    基于遙感的火災(zāi)研究主要包括火點(diǎn)(active fire)監(jiān)測(cè)[21-22]與火燒跡地(burn scars)提取[23-24]兩類,兩者均可以產(chǎn)生火燒跡地?cái)?shù)據(jù)[25]?;瘘c(diǎn)監(jiān)測(cè)主要是基于火災(zāi)發(fā)生時(shí)的溫度異常,監(jiān)測(cè)衛(wèi)星過(guò)境時(shí)可能發(fā)生火災(zāi)的像元,對(duì)火點(diǎn)進(jìn)行實(shí)時(shí)觀測(cè)。這一過(guò)程的主要目的在于捕捉火災(zāi)發(fā)生的時(shí)間以及位置信息,雖然也可以產(chǎn)生火燒跡地信息,但結(jié)果并不可靠[26-27]?;馃E地提取通過(guò)對(duì)比火災(zāi)發(fā)生前后的光譜反射特征或者光譜指數(shù)特征變化來(lái)識(shí)別火燒跡地,從而達(dá)到提取火燒跡地的目的[28-29]。這一方法的缺陷在于,一些與火災(zāi)發(fā)生具有類似光譜特征的事件較難區(qū)分,如洪水、森林砍伐以及農(nóng)作物收獲等導(dǎo)致的光譜特征變化。本文的火燒跡地提取方法以后者為基礎(chǔ),并進(jìn)行了改進(jìn),以消除類似的混淆事件。

    1.3.1 判別指數(shù)選取

    火燒跡地識(shí)別方法是通過(guò)比較火災(zāi)發(fā)生前后的光譜特征變化來(lái)提取火燒面積。因而需要選取適合的光譜指數(shù)來(lái)進(jìn)行表征,如NDVI(Normalized Difference Vegetation Index)、BBFI(Burned Boreal Forest Index)、GEMI(Global Environmental Monitoring Index)以及BAI (Burned Area Index)等。其中,NDVI的應(yīng)用最為廣泛。NDVI能夠很好的對(duì)植被覆蓋進(jìn)行描述[30,31],但研究表明NDVI在植被覆蓋度較高的地區(qū)容易達(dá)到飽和[32],且在火燒跡地信息提取中存在較大的潛在誤差[33]。因此,選用GEMI作為識(shí)別火燒跡地的主要判別指數(shù),其計(jì)算公式如下:

    GEMI=η×(1-0.25η)-(ρred-0.125)/(1-ρred)

    η=(2(ρnir2-ρred2)+1.5ρnir+0.5ρred)/(ρnir+ρred+0.5)

    (1)

    式中,ρnir以及ρred為近紅外波段和紅光波段。火災(zāi)發(fā)生后GEMI表現(xiàn)出顯著的下降。

    為了避免采用單一光譜指數(shù)所帶來(lái)的潛在誤差,選擇了另一個(gè)光譜指數(shù)BAI來(lái)作進(jìn)一步的限定,其計(jì)算公式如下:

    BAI=1/((ρnir-ρcnir)2+(ρred-ρcred)2)

    (2)

    式中,ρcnir和ρcred分別被設(shè)定為0.06與0.1?;馂?zāi)發(fā)生后BAI值表現(xiàn)出顯著的上升。

    除此之外,在比較火災(zāi)發(fā)生前后光譜特征變化的同時(shí),為了考慮火災(zāi)發(fā)生時(shí)的熱學(xué)特性,即溫度異常,將MODIS火產(chǎn)品數(shù)據(jù)作為一個(gè)輸入波段加入到判別流程,以提高判別精度。

    1.3.2 判別流程

    火燒跡地的識(shí)別流程主要分為兩個(gè)階段:首先,設(shè)定較為嚴(yán)格的判別閾值以提取火燒的核心像元——即火災(zāi)最有可能發(fā)生的像元。這一階段的主要目標(biāo)在于盡可能的減少錯(cuò)判誤差,因而需要對(duì)火災(zāi)發(fā)生前后的光譜指數(shù)變化設(shè)定嚴(yán)格的閾值,并且同時(shí)用使用MODIS火災(zāi)產(chǎn)品進(jìn)行篩選,以表示提取像元在相關(guān)植被指數(shù)變化前曾出現(xiàn)溫度異常。這一過(guò)程同時(shí)考慮了火災(zāi)發(fā)生前后地表植被的突變以及火災(zāi)發(fā)生時(shí)的溫度異常,從而可以與其他造成地表突變的因素相區(qū)分。其次,對(duì)第一階段提取的核心像元15公里范圍內(nèi)的光譜指數(shù)變化特征進(jìn)行判別,設(shè)定較為寬松的閾值,以盡可能減少漏判誤差。判別流程及判別條件如圖2所示。

    圖2 算法流程圖Fig.2 Flowchart of algorithm

    第一階段的提取過(guò)程以GEMI、BAI、以及MOD14A2產(chǎn)品為基礎(chǔ),具體的判別條件如下所示:

    首先,火災(zāi)發(fā)生之前的GEMI值必須大于一定的閾值,以確保判別區(qū)域?yàn)橹脖桓采w。

    GEMIt-1>0.170

    (3)

    式中,t為時(shí)間(下同)。選擇數(shù)據(jù)為MODIS 8d合成數(shù)據(jù),每年共34期數(shù)據(jù),因此t的范圍為:0

    火災(zāi)發(fā)生后,GEMI值必須表現(xiàn)出顯著的下降,且這一下降過(guò)程必須持續(xù)一定的時(shí)間,以區(qū)分由云污染等造成的GEMI值的短暫下降。這一過(guò)程通過(guò)以下兩個(gè)判別條件來(lái)實(shí)現(xiàn):

    (GEMIt-GEMIt-1)/GEMIt<-0.1

    (4)

    (GEMIt+2-GEMIt-1)/GEMIt+2<-0.1

    (5)

    然后,使用BAI指數(shù)來(lái)對(duì)火燒像元做進(jìn)一步的限定?;馂?zāi)發(fā)生后,BAI值顯著增加,其判別條件如下:

    BAIt>250且 BAIt-1>200

    (6)

    最后,使用MODIS 火產(chǎn)品數(shù)據(jù)來(lái)對(duì)火燒像元進(jìn)行掩膜,以保證光譜指數(shù)變化前,所提取像元表現(xiàn)出溫度異常的特征。

    ρt>6 或者 ρt-1>6

    (7)

    式中,ρ為MODIS火產(chǎn)品數(shù)據(jù)像元值。

    第二階段的判別過(guò)程以第一階段提取的核心像元為基礎(chǔ),采用較為寬松的閾值來(lái)對(duì)鄰近像元進(jìn)行判別。在對(duì)研究區(qū)的火災(zāi)發(fā)生特征進(jìn)行分析之后,距離核心像元的最大距離被設(shè)定為15 km。第二階段的火燒跡地信息提取,僅對(duì)核心像元15 km范圍內(nèi)像元進(jìn)行判別,判別條件如下:

    GEMIt-GEMIt-1<-0.03

    (8)

    GEMIt+1-GEMIt-1<-0.02

    (9)

    GEMIt+2-GEMIt-1<0

    (10)

    GEMIt+1-GEMIt≤0

    (11)

    BAIt>250

    (12)

    最后,將兩個(gè)階段的提取結(jié)果進(jìn)行合成。采用一個(gè)3×3的變換核,對(duì)合成結(jié)果進(jìn)行濾波處理,消除提取過(guò)程中產(chǎn)生的小斑塊。

    2 結(jié)果與分析

    2.1 精度驗(yàn)證

    由于研究區(qū)涉及境外地區(qū),驗(yàn)證數(shù)據(jù)難以獲取??紤]到以上提取方法的基礎(chǔ)為火燒前后植被指數(shù)的變化以及溫度的異常,使得該方法對(duì)火燒跡地的提取具有普適性,從而可以采用選擇典型驗(yàn)證區(qū)的方法對(duì)算法進(jìn)行驗(yàn)證。故以黑龍江省為典型研究區(qū),對(duì)提取結(jié)果進(jìn)行精度驗(yàn)證。驗(yàn)證數(shù)據(jù)來(lái)源于相關(guān)林業(yè)部門2000—2005年的火災(zāi)統(tǒng)計(jì)數(shù)據(jù),包括火災(zāi)發(fā)生的時(shí)間、地點(diǎn)、經(jīng)緯度信息以及過(guò)火面積等。

    圖3 黑龍江省火災(zāi)發(fā)生位置(2000—2005年)Fig.3 Fire position of Heilongjiang province (2000—2005)

    鑒于MODIS產(chǎn)品空間分辨率以及火燒跡地信息提取后進(jìn)行去除小斑塊的濾波處理的需要,對(duì)過(guò)火面積小于60 hm2(約3×3個(gè)像元)的火災(zāi)進(jìn)行剔除,最終得到黑龍江省2000—2005年的火災(zāi)驗(yàn)證數(shù)據(jù)(圖3)。

    利用2000年—2005年火災(zāi)發(fā)生的經(jīng)緯度信息(圖3)對(duì)提取結(jié)果進(jìn)行錯(cuò)判以及漏判分析(表1)。由于驗(yàn)證數(shù)據(jù)僅提供了火災(zāi)發(fā)生的位置,因而不能對(duì)提取結(jié)果進(jìn)行空間化(逐像元)的誤差分析。以火災(zāi)發(fā)生的位置信息為參照,對(duì)提取結(jié)果進(jìn)行分析,兩者一致則認(rèn)為提取結(jié)果正確。如果在標(biāo)有火災(zāi)發(fā)生的位置沒(méi)有提取出火燒跡地信息,被認(rèn)為是漏判;相反,在沒(méi)有標(biāo)出火災(zāi)發(fā)生的位置,卻提取出火燒跡地信息,認(rèn)為是錯(cuò)判。

    表1 火燒跡地提取驗(yàn)證表

    此外,將提取結(jié)果的面積進(jìn)行匯總與驗(yàn)證數(shù)據(jù)進(jìn)行了比較(表1)。結(jié)果顯示,2000—2005年每年均有一定的漏判以及錯(cuò)判誤差存在,且提取的火燒跡地面積均小于驗(yàn)證數(shù)據(jù),總體精度為71%。其中,提取面積精度最高為2002年,達(dá)84%;提取面積精度最低為2003年,精度為61%。較之以往的研究[35],精度有所提高。

    2.2 黑龍江流域火燒跡地分布特征分析

    2.2.1 黑龍江流域火燒跡地信息提取

    圖4 黑龍江流域火燒跡地分布(2000—2011) Fig.4 Distribution of burned area in Heilongjiang basin (2000—2011)

    使用以上所驗(yàn)證的火燒跡地提取方法,以MOD09Q1數(shù)據(jù)以及MOD14A2數(shù)據(jù)為基礎(chǔ),對(duì)黑龍江流域2000—2011年的火燒跡地信息進(jìn)行了提取,得到研究區(qū)長(zhǎng)時(shí)間序列火燒跡地分布圖(圖4)。

    2.2.2 火燒跡地特征分析

    將黑龍江流域的火燒跡地面積進(jìn)行逐年匯總,從而得到流域逐年的火燒跡地面積統(tǒng)計(jì)數(shù)據(jù)(圖5)。

    圖5 黑龍江流域火燒跡地面積變化Fig.5 Burned area dynamic in Heilongjiang basin

    黑龍江流域2000—2011年受火災(zāi)影響較為嚴(yán)重,年均過(guò)火面積達(dá)53.21萬(wàn)hm2。火災(zāi)發(fā)生最嚴(yán)重的年份為2003年,面積為146.79萬(wàn)hm2;而受火災(zāi)影響最小的年份為2010年,過(guò)火面積僅有18.39萬(wàn)hm2,差距較大。火災(zāi)發(fā)生較為嚴(yán)重的年份還包括2008年,過(guò)火面積也超過(guò)了百萬(wàn)公頃,達(dá)119.41萬(wàn)hm2。其他年份受火災(zāi)影響較為平均,其中2005年相對(duì)較為嚴(yán)重,面積為62.41萬(wàn)hm2;其次為2001、2002、2004年以及2011年,火燒跡地面積分別為41.55萬(wàn)hm2、43.83萬(wàn)hm2、43.59萬(wàn)hm2以及43.23萬(wàn)hm2;最后為2000、2006、2007年以及2009年,火燒跡地面積分別為33.08萬(wàn)hm2、38.87萬(wàn)hm2、24.01萬(wàn)公頃以及23.41萬(wàn)hm2。

    從火燒跡地的空間特征來(lái)看,火燒跡地的分布與森林覆被密切相關(guān),主要分布于黑龍江流域的中、高緯度地區(qū)。俄羅斯境內(nèi)的火燒跡地分布較為均勻,原因在于其高的植被覆蓋度。中國(guó)境內(nèi)火燒跡地主要分布于大、小興安嶺以及長(zhǎng)白山地區(qū),其中以黑龍江省受影響最為嚴(yán)重。此外,蒙古和朝鮮境內(nèi)也有少部分火燒跡地分布。

    3 結(jié)論與討論

    長(zhǎng)時(shí)間序列火燒跡地?cái)?shù)據(jù)是區(qū)域或者全球尺度下森林火災(zāi)相關(guān)研究的重要基礎(chǔ)信息。本文以黑龍江流域?yàn)檠芯繀^(qū),利用MODIS時(shí)序數(shù)據(jù)對(duì)其2000—2011年的火燒跡地信息進(jìn)行了提取,主要結(jié)論如下:

    低空間分辨率高時(shí)間分辨率的遙感數(shù)據(jù)是區(qū)域或者全球尺度下火燒跡地信息提取的主要數(shù)據(jù)源,通過(guò)建立相應(yīng)的算法,可以實(shí)現(xiàn)長(zhǎng)時(shí)間序列火燒跡地信息的提取過(guò)程。

    綜合考慮火災(zāi)發(fā)生前后的植被變化(光譜指數(shù)變化)與火災(zāi)發(fā)生時(shí)的熱學(xué)特征能夠更為有效的對(duì)火燒跡地信息進(jìn)行提取,提高提取精度。

    黑龍江流域2000—2011年受火災(zāi)影響較為嚴(yán)重,年平均產(chǎn)生火燒跡地53.21萬(wàn)hm2,受火災(zāi)影響最大年份與最小年份之間的火燒跡地面積差距較大。

    本文基于MODIS數(shù)據(jù)空間分辨率以及提取結(jié)果濾波處理的需求,將過(guò)火面積大于60 hm2的森林火災(zāi)作為驗(yàn)證數(shù)據(jù),較之Emilio等使用大于200 hm2的火災(zāi)[34]作為驗(yàn)證數(shù)據(jù),提高了對(duì)算法精度的要求。結(jié)果顯示主要的漏判誤差仍來(lái)源于100 hm2左右的森林火災(zāi)。證明由于遙感數(shù)據(jù)空間分辨率的局限,算法對(duì)于面積較小的火燒跡地提取具有一定難度。而空間分辨率相對(duì)較高的遙感數(shù)據(jù),如TM數(shù)據(jù),其時(shí)間分辨率卻難以滿足火燒跡地信息提取的要求。遙感數(shù)據(jù)融合可以較好的解決這一問(wèn)題,如何采用融合之后的高空間分辨率以及時(shí)間分辨率的遙感數(shù)據(jù)進(jìn)行更為細(xì)致的火燒跡地提取有待深入研究。

    [1] Dubinin M, Potapov P, Lushchekina A, Radeloff V C. Reconstructing long time series of burned areas in arid grasslands of southern Russia by satellite remote sensing. Remote Sensing of Environment, 2010, 114(8): 1638- 1648.

    [2] Fraser R H, Li Z. Estimating fire-related parameters in boreal forest using SPOT VEGETATION. Remote Sensing of Environment, 2002, 82(1): 95- 110.

    [3] 呂愛(ài)鋒, 田漢勤, 劉永強(qiáng). 火干擾與生態(tài)系統(tǒng)碳循環(huán). 生態(tài)學(xué)報(bào), 2005, 25(10): 2734- 2743.

    [4] Amiro B D, Orchansky A L, Barr A G, Black T A, Chambers S D, Chapin III F S, Goulden M L, Litvak M, Liu H P, McCaughey J H, McMillan A, Randerson J T. The effect of post-fire stand age on the boreal forest energy balance. Agricultural and Forest Meteorology, 2006, 140(1/4): 41- 50.

    [5] Kasischke E S, Turetsky M R. Recent changes in the fire regime across the North American boreal region—Spatial and temporal patterns of burning across Canada and Alaska. Geophysical Research Letters, 2006, 33(9): 1- 5.

    [6] Morton D C, Collatz G J, Wang D, Randerson J T, Giglio L, Chen Y. Satellite-based assessment of climate controls on US burned area. Biogeosciences, 2013, 10(1): 247- 260.

    [7] Flannigan M D, Logan K A, Amiro B D, Skinner W R, Stocks B J. Future area burned in Canada. Climatic Change, 2005, 72(1/2): 1- 16.

    [8] IPCC. Climate Change 2007: The Physical Science Basis. Cambridge: Cambridge University Press, 2007.

    [9] Gillet N P, Weaver A J, Zwiers F W, Flannigan M D. Detecting the effect of climate change on Canadian forest fires. Geophysical Research Letters, 2004, 31(18), L18211, doi: 18210.11029/12004GL02 0876.

    [10] Mouillot F, Schultz M G, Yue C, Cadule P, Tensey K, Ciais P, Chuvieco E. Ten years of global burned area products from spaceborne remote sensing-A review: Analysis of user needs and recommendations for future developments. International Journal of Applied Earth Observation and Geoinformation, 2014, 26: 64- 79.

    [11] Vafeidis A T, Drake N A, Wainwright J. A proposed method for modelling the hydrologic response of catchments to burning with the use of remote sensing and GIS. Catena, 2007, 70(3): 396- 409.

    [12] Pu R L, Li Z Q, Gong P, Csiszar I, Fraser R, Hao W M, Kondragunta S, Weng F Z. Development and analysis of a 12-year daily 1-km forest fire dataset across North America form NOAA/AVHRR data. Remote Sensing of Environment, 2007, 108(2): 198- 208.

    [13] Sukhinin A I, French N H F, Kasischke E S, Hewson J H, Soja A J, Csiszar I A, Hyer E J, Loboda T, Gonrad S G, Romasko V I, Pavichenko E A, Miskiv S I, Slinkina O A. AVHRR-based mapping of fires in Russia: New products for fire management and carbon cycle studies. Remote Sensing of Environment, 2004, 93(4): 546- 564.

    [14] Oadilla M, Stehman S V, Litago J, Chuvieco E. Assessing the temporal stability of the accuracy of a time series of burned area products. Remote Sensing, 2014, 6(3): 2050- 2068.

    [15] Roy D P, Boschetti L, Justice C O, Ju J. The collection 5 MODIS burned area product—Global evaluation by comparison with the MODIS active fire product. Remote Sensing of Environment, 2008, 112(9): 3690- 3707.

    [16] Barbosa P M, Gregoire J M, Pereira J M C. An algorithm for extracting burned areas from time series of AVHRR GAC data applied at a continental scale. Remote Sensing of Environment, 1999, 69(3): 253- 263.

    [17] Martin M P, Chuvieco E. Mapping and evaluation of burned land from multi-temporal analysis of AVHRR NDVI images. EARSeL Advances in Remote Sensing, 1995, 4(3): 7- 13.

    [18] Pereira J M C. A comparative evaluation of NOAA/AVHRR vegetation indexes for burned surface detection and mapping. IEEE Transactions on Geoscience and Remote Sensing, 1999, 37(1): 217- 226.

    [19] Giglio L, Descloitres J, Justice C O, Kaufman Y J. An enhanced contextual fire detection algorithm for MODIS. Remote Sensing of Environment, 2003, 87(2/3): 273- 282.

    [20] Dozier J. A method for satellite identification of surface temperature fields of subpixel resolution. Remote Sensing of Environment, 1981, 11: 221- 229.

    [21] Movaghati S, Samadzadegan F, Azizi A. An agent-based approach for regional forest fire detection using MODIS data. Journal of Applied Sciences, 2009, 9(20): 3672- 3681.

    [22] 王彩云, 扎西頓珠, 陳濤. EOS/MODIS數(shù)據(jù)在西藏森林火災(zāi)中的應(yīng)用. 高原山地氣象研究, 2010, 30(3): 66- 69.

    [23] Ruiz J A M, Riano D, Arbelo M, French N H F, Ustin S L, Whiting M L. Burned area mapping time series in Canada (1984- 1999) from NOAA-AVHRR LTDR: A comparison with other remote sensing products and fire perimeters. Remote Sensing of Environment, 2012, 117: 407- 414.

    [24] Giglio L, Loboda T, Roy D P, Quayle B, Justice C O. An active-fire based burned area mapping algorithm for the MODIS sensor. Remote Sensing of Environment, 2009, 113(2): 408- 420.

    [25] Zhang X Y, Kondragunta S, Quayle B. Estimation of biomass burned areas using multiple-satellite-observed active fires. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(11): 4469- 4482.

    [26] Kasischke E S, Hewson J H, Stock B, van der Werf G, Randerson J. The use of ATSR active fire counts for estimating relative patterns of biomass burning-a study from the boreal forest region. Geophysical Research Letters, 2003, 30(18), doi: 10.1029/2003GL017859.

    [27] Giglio L, Van der Werf G R, Randerson J T, Collatz G J, Kasibhatla P. Global estimation of burned area using MODIS active fire observations. Atmospheric Chemistry and Physics, 2006, 6(4): 957- 974.

    [28] Key C H. Ecological and sampling constraints on defining landscape fire severity. Fire Ecology, 2006, 2(2): 34- 59.

    [29] Roy D P, Jin Y, Lewis P E, Justice C O. Prototyping a golbal algorithm for systematic fire-affected area mapping using MODIS time series data. Remote Sensing of Environment, 2005, 97(2): 137- 162.

    [30] Fraser R H, Li Z, Cihlar J. Hotspot and NDVI differencing synergy (HANDS): a new technique for burned area mapping over boreal forest. Remote Sensing of Environment, 2000, 74(3): 362- 376.

    [31] Kucera J, Yasuoka Y, Dye D G. Creating a forest fire database for the Far East Asia using NOAA/AVHRR observation. International Journal of Remote Sensing, 2005, 26(11): 2423- 2439.

    [32] Cai H Y, Zhang S W, Bu K, Yang J C, Chang L P. Intergrating geographical data and phenological characteristics derived from MODIS data for improving land over mapping. Journal of Geographical Sciences, 2011, 21(4): 705- 718.

    [33] Chuvieco E, Martín M P, Palacios A. Assessment of different spectral indices in the red-near-infrared spectral domain for burned land discrimination. International Journal of Remote Sensing, 2002, 23(23): 5103- 5110.

    [34] Chuvieco E, Englefield P, Trishchenko A P, Luo Y. Generation of long time series of burn area maps of the boreal forest from NOAA-AVHRR composite data. Remote Sensing of Environment, 2008, 112(5): 2381- 2396.

    Burned area mapping for Heilongjiang basin based on MODIS time series data

    YANG Wei1, ZHANG Shuwen2,*, JIANG Xiaoli3

    1SchoolofGeographicalSciences,TaiyuanNormalUniversity,Jinzhong030619,China2NortheastInstituteofGeographyandAgroecology,ChineseAcademyofScience,Changchun130102,China3InstituteofUrbanandDistrictDevelopment,TaiyuanNormalUniversity,Jinzhong030619,China

    Fire is an important natural disturbance that affects several ecosystems and is also one of the main factors of the terrestrial carbon cycle. As fire modifies the structure and composition of vegetation, it is considered to be an important land management tool. Burned area mapping is an essential step in forest fire research to investigate the relationship between forest fire and climate change and the effect of forest fire on carbon budgets. Traditional data collection of forest fires in field- which are statistically recorded are difficult to manipulate over a large area. The development of the remote sensing technique provides a labor-efficient method for research of land surface processes. At the regional or global scale, in order to obtain a long-time series of burned area maps, a moderate spatial resolution with high temporal resolution remote sensing data is considered as the best alternative. Currently, the most widely used remote sensing data are Advanced Very High Resolution Radiometer (AVHRR) images and Moderate-Resolution Imaging Spectroradiometer (MODIS) images. Although the AVHRR provides continuous observations for burned area analyses, some studies have identified several sources of potential errors in burned area discrimination from this sensor, mainly due to its radiometric instability, cloud obscuration, and transmission problems. Most of these problems have been notably reduced in the MODIS sensor, which offers greater spectral, spatial, and radiometric resolution than the AVHRR. This study proposes an algorithm to map areas burned by forest fire using MODIS time series data in Heilongjiang Valley, China. The algorithm is divided into two steps: First, the “core” pixels were extracted to represent the most possible burned pixels based on a comparison of the temporal change of the Global Environmental Monitoring Index (GEMI), the Burned Area Index (BAI), and the MODIS active fire products between pre- and post-fire spatial patterns. Second, a 15-km distance was set to extract the entire burned area near the “core” pixels. These more relaxed conditions were used to identify the fire pixels for reducing the omission error as much as possible. The algorithm comprehensively considered the thermal characteristics and the spectral change between pre- and post-fire spatial patterns, which were represented by the MODIS fire products and the spectral index, respectively. Heilongjiang province in China was selected as the typical study area to validate the accuracy of the algorithm. The results showed that with the use of the MODIS fire products, the accuracy of the algorithm was improved, with an overall accuracy of 71% and a highest accuracy of 84%. Consequently, the algorithm used in this study produced a long-time series of burned area maps of the study area from 2000 to 2011 with a relatively high accuracy. According to the burned area maps, the study area has been seriously affected by fire disasters on average of 0.53 million ha of burned land each year. The most affected years were 2003 and 2008 with burned areas exceeding 1 million ha. The least affected year was 2010 with a burned area of just 0.18 million ha. The relatively large disparity between the maximum and minimum values of the areas burned by forest fire indicates that there is a fluctuation in the severity of disaster during the studied period.

    burned area; MODIS; GEMI; BAI

    中國(guó)科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA05090310)

    2013- 12- 31;

    日期:2014- 11- 03

    10.5846/stxb201312313076

    *通訊作者Corresponding author.E-mail: zhangshuwen@neigae.ac.cn

    楊偉, 張樹(shù)文, 姜曉麗.基于MODIS時(shí)序數(shù)據(jù)的黑龍江流域火燒跡地提取.生態(tài)學(xué)報(bào),2015,35(17):5866- 5873.

    Yang W, Zhang S W, Jiang X L.Burned area mapping for Heilongjiang basin based on MODIS time series data.Acta Ecologica Sinica,2015,35(17):5866- 5873.

    猜你喜歡
    跡地火燒黑龍江
    封二、封三:寒蔥嶺戰(zhàn)跡地 (伏擊戰(zhàn)) 新考
    Homeless children in India publish newspaper
    火燒毒品
    南方周末(2019-06-27)2019-06-27 19:47:54
    Write a caption
    我國(guó)科學(xué)家在國(guó)際上首發(fā)30 m分辨率全球火燒跡地?cái)?shù)據(jù)產(chǎn)品
    驢肉火燒
    火燒跡地造林技術(shù)
    基于Landsat影像提取大興安嶺區(qū)域火燒跡地★
    山西建筑(2017年15期)2017-07-01 22:03:22
    不怕火燒的手帕
    莫愁(2017年15期)2017-05-11 01:21:55
    黑龍江棄風(fēng)博弈
    能源(2016年1期)2016-12-01 05:10:04
    尾随美女入室| 国产精品.久久久| 天堂av国产一区二区熟女人妻| 久久精品国产鲁丝片午夜精品| 国产一区二区三区av在线| 日韩欧美精品v在线| 日本黄大片高清| 国产精品99久久久久久久久| www.av在线官网国产| 人体艺术视频欧美日本| 国产午夜福利久久久久久| 国产永久视频网站| 亚洲精品视频女| 日韩一本色道免费dvd| 免费高清在线观看视频在线观看| 欧美日韩在线观看h| 午夜福利在线观看免费完整高清在| 国产毛片a区久久久久| 99久久九九国产精品国产免费| 国产女主播在线喷水免费视频网站 | 美女脱内裤让男人舔精品视频| 久久99精品国语久久久| 国产真实伦视频高清在线观看| 日韩制服骚丝袜av| 深爱激情五月婷婷| 在现免费观看毛片| 国产在视频线精品| 最近的中文字幕免费完整| 亚洲av在线观看美女高潮| 尤物成人国产欧美一区二区三区| 精品久久久久久久久久久久久| 国产精品一区二区三区四区免费观看| 日韩不卡一区二区三区视频在线| 69人妻影院| 亚洲一区高清亚洲精品| 伦理电影大哥的女人| 久久国产乱子免费精品| 久久久精品欧美日韩精品| 亚洲最大成人av| 美女黄网站色视频| 国产一区二区亚洲精品在线观看| 一区二区三区免费毛片| 麻豆精品久久久久久蜜桃| 国产一级毛片在线| 天堂网av新在线| 亚洲丝袜综合中文字幕| 亚洲天堂国产精品一区在线| 国产激情偷乱视频一区二区| 久久久久久久久久久免费av| 特级一级黄色大片| 午夜免费男女啪啪视频观看| 人妻一区二区av| 日韩,欧美,国产一区二区三区| 亚洲av二区三区四区| 国产免费福利视频在线观看| 国产三级在线视频| 亚洲av.av天堂| 床上黄色一级片| 偷拍熟女少妇极品色| 一个人观看的视频www高清免费观看| 18禁裸乳无遮挡免费网站照片| 国产淫片久久久久久久久| 国产精品久久久久久精品电影小说 | 国产伦在线观看视频一区| 国内精品一区二区在线观看| 国产一级毛片七仙女欲春2| 天堂网av新在线| 大又大粗又爽又黄少妇毛片口| 麻豆成人av视频| 日本免费在线观看一区| 亚洲国产高清在线一区二区三| 插逼视频在线观看| 亚洲成人精品中文字幕电影| 国产av国产精品国产| 国产精品久久久久久久久免| 三级男女做爰猛烈吃奶摸视频| 亚洲综合色惰| 午夜爱爱视频在线播放| 小蜜桃在线观看免费完整版高清| 国产亚洲5aaaaa淫片| 少妇裸体淫交视频免费看高清| 免费看日本二区| 国产一区二区三区综合在线观看 | 亚洲经典国产精华液单| av免费在线看不卡| 简卡轻食公司| 亚洲在久久综合| 26uuu在线亚洲综合色| 国产免费视频播放在线视频 | 一级片'在线观看视频| 精品欧美国产一区二区三| 国产成人精品一,二区| 国内精品美女久久久久久| 亚洲av成人av| 久久99热6这里只有精品| 免费看光身美女| 欧美日韩精品成人综合77777| 少妇猛男粗大的猛烈进出视频 | 人人妻人人澡人人爽人人夜夜 | 亚洲av不卡在线观看| 九草在线视频观看| 亚洲欧美成人综合另类久久久| 美女内射精品一级片tv| 亚洲三级黄色毛片| 22中文网久久字幕| 观看美女的网站| 国产av在哪里看| 久久精品久久精品一区二区三区| 亚洲国产精品国产精品| 国产有黄有色有爽视频| 国产 亚洲一区二区三区 | 在线观看免费高清a一片| 亚洲图色成人| 边亲边吃奶的免费视频| 综合色av麻豆| 欧美高清性xxxxhd video| 老司机影院毛片| 亚洲精品乱久久久久久| 国内精品宾馆在线| 麻豆成人av视频| 国内揄拍国产精品人妻在线| 国产老妇伦熟女老妇高清| 国产欧美日韩精品一区二区| 欧美bdsm另类| 精品一区二区三区视频在线| 欧美激情久久久久久爽电影| av免费在线看不卡| 国产69精品久久久久777片| 少妇熟女aⅴ在线视频| 久久久久精品久久久久真实原创| 日韩亚洲欧美综合| 欧美人与善性xxx| 久久这里只有精品中国| 久久国内精品自在自线图片| 尾随美女入室| 美女黄网站色视频| 亚洲一区高清亚洲精品| 精品熟女少妇av免费看| 国产精品精品国产色婷婷| 中文字幕制服av| 欧美区成人在线视频| 五月玫瑰六月丁香| 极品教师在线视频| 欧美xxⅹ黑人| 天美传媒精品一区二区| 国产成人精品一,二区| freevideosex欧美| 久久97久久精品| 又爽又黄a免费视频| 国内精品美女久久久久久| 舔av片在线| 亚洲av不卡在线观看| 国产精品久久久久久久电影| 亚洲最大成人av| 能在线免费观看的黄片| 一级毛片aaaaaa免费看小| 国产黄片美女视频| 亚洲激情五月婷婷啪啪| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产真实伦视频高清在线观看| 色哟哟·www| 禁无遮挡网站| 99热这里只有精品一区| 欧美日本视频| 国产精品一区二区性色av| 亚洲综合色惰| 能在线免费看毛片的网站| 精品一区二区三卡| av播播在线观看一区| videossex国产| 啦啦啦韩国在线观看视频| 国产精品福利在线免费观看| 99久久精品国产国产毛片| 午夜免费男女啪啪视频观看| 丰满人妻一区二区三区视频av| 能在线免费观看的黄片| 最近中文字幕2019免费版| 精品久久久久久久久av| 成年av动漫网址| 欧美三级亚洲精品| 欧美zozozo另类| 一级爰片在线观看| 五月伊人婷婷丁香| 97精品久久久久久久久久精品| 伦理电影大哥的女人| 久久久国产一区二区| 国产69精品久久久久777片| 在线a可以看的网站| 麻豆国产97在线/欧美| 一级av片app| 高清午夜精品一区二区三区| 视频中文字幕在线观看| 日韩电影二区| 国产伦精品一区二区三区视频9| 精品欧美国产一区二区三| 亚洲在久久综合| 国产精品人妻久久久影院| 国产精品国产三级国产专区5o| 波野结衣二区三区在线| 偷拍熟女少妇极品色| 老女人水多毛片| freevideosex欧美| 好男人在线观看高清免费视频| 国产伦一二天堂av在线观看| 国产一区二区亚洲精品在线观看| 国产色爽女视频免费观看| 亚洲成人精品中文字幕电影| 国产精品久久久久久精品电影小说 | 天堂俺去俺来也www色官网 | 成人无遮挡网站| 成人午夜精彩视频在线观看| av免费在线看不卡| 国产91av在线免费观看| 国产精品国产三级专区第一集| 亚洲高清免费不卡视频| 高清视频免费观看一区二区 | 国产精品人妻久久久久久| 搡女人真爽免费视频火全软件| 最近2019中文字幕mv第一页| 国产熟女欧美一区二区| 日本av手机在线免费观看| 男女那种视频在线观看| 亚洲怡红院男人天堂| 免费在线观看成人毛片| 亚洲av日韩在线播放| 国产av国产精品国产| av专区在线播放| 又大又黄又爽视频免费| 青春草国产在线视频| 欧美区成人在线视频| 噜噜噜噜噜久久久久久91| 亚洲av在线观看美女高潮| 蜜臀久久99精品久久宅男| 成人无遮挡网站| 国产伦在线观看视频一区| xxx大片免费视频| 晚上一个人看的免费电影| 午夜福利高清视频| 亚洲欧美清纯卡通| 黄色一级大片看看| 听说在线观看完整版免费高清| 乱码一卡2卡4卡精品| 日本与韩国留学比较| 精品欧美国产一区二区三| 国产精品一区二区三区四区免费观看| 青青草视频在线视频观看| 天堂网av新在线| 最新中文字幕久久久久| 国产精品一二三区在线看| 夜夜看夜夜爽夜夜摸| 日本免费a在线| 成人综合一区亚洲| 一级毛片 在线播放| 干丝袜人妻中文字幕| 大话2 男鬼变身卡| 成人无遮挡网站| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 精品一区二区三区视频在线| 黄色日韩在线| 一区二区三区四区激情视频| 一级毛片aaaaaa免费看小| 菩萨蛮人人尽说江南好唐韦庄| 久久久久性生活片| 免费观看在线日韩| videossex国产| 内射极品少妇av片p| 国产亚洲91精品色在线| 乱码一卡2卡4卡精品| 3wmmmm亚洲av在线观看| 国产一级毛片在线| 亚洲av中文av极速乱| 国产精品熟女久久久久浪| 国产白丝娇喘喷水9色精品| 色哟哟·www| 国产在视频线在精品| 深夜a级毛片| 狂野欧美白嫩少妇大欣赏| 插逼视频在线观看| 日韩欧美精品免费久久| 成年版毛片免费区| 久久精品久久精品一区二区三区| 久久久久久久久久成人| 中国国产av一级| 午夜精品一区二区三区免费看| 亚洲国产av新网站| 欧美极品一区二区三区四区| 五月伊人婷婷丁香| 91精品国产九色| 2022亚洲国产成人精品| 国产成人aa在线观看| av.在线天堂| 丝袜美腿在线中文| 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 成年女人在线观看亚洲视频 | 五月天丁香电影| 久久久久久久久久人人人人人人| 国产免费视频播放在线视频 | 男插女下体视频免费在线播放| 伦理电影大哥的女人| 观看美女的网站| 美女被艹到高潮喷水动态| 美女国产视频在线观看| 日日啪夜夜爽| 日日撸夜夜添| 国产亚洲av片在线观看秒播厂 | 国产单亲对白刺激| 国产黄频视频在线观看| 国产乱来视频区| 国产女主播在线喷水免费视频网站 | 搞女人的毛片| 免费大片18禁| 精品少妇黑人巨大在线播放| 美女xxoo啪啪120秒动态图| 国产淫语在线视频| 天堂√8在线中文| 老女人水多毛片| 精品久久久久久电影网| 久久久久久久久久久丰满| 91av网一区二区| 人妻系列 视频| 91久久精品电影网| 久久久久精品性色| 一区二区三区高清视频在线| 乱系列少妇在线播放| 免费观看在线日韩| 禁无遮挡网站| 大话2 男鬼变身卡| 国产片特级美女逼逼视频| av.在线天堂| 亚洲综合精品二区| 全区人妻精品视频| 国内精品宾馆在线| 午夜福利网站1000一区二区三区| 五月天丁香电影| 免费观看在线日韩| 插阴视频在线观看视频| av在线亚洲专区| 免费不卡的大黄色大毛片视频在线观看 | 内地一区二区视频在线| 日韩中字成人| 成人亚洲精品av一区二区| 亚洲自拍偷在线| 国产视频内射| 国产精品福利在线免费观看| 亚洲精品中文字幕在线视频 | 中文乱码字字幕精品一区二区三区 | 成人午夜高清在线视频| 国内精品一区二区在线观看| 91精品伊人久久大香线蕉| 少妇熟女aⅴ在线视频| 一夜夜www| 国产探花极品一区二区| 最近中文字幕高清免费大全6| 中国美白少妇内射xxxbb| 久久久精品欧美日韩精品| 免费看不卡的av| 成人特级av手机在线观看| 亚洲av中文字字幕乱码综合| 激情 狠狠 欧美| 亚洲精品一二三| 激情五月婷婷亚洲| 色5月婷婷丁香| 美女xxoo啪啪120秒动态图| 欧美日本视频| 久久99热这里只频精品6学生| 搞女人的毛片| 最近中文字幕高清免费大全6| 国产探花极品一区二区| 日日摸夜夜添夜夜爱| 国产一区二区三区综合在线观看 | 亚洲最大成人av| 欧美97在线视频| 最新中文字幕久久久久| 亚洲精品成人久久久久久| 午夜亚洲福利在线播放| 中文字幕人妻熟人妻熟丝袜美| 午夜日本视频在线| 日日干狠狠操夜夜爽| 成人一区二区视频在线观看| 一级毛片电影观看| 色5月婷婷丁香| 视频中文字幕在线观看| 大香蕉久久网| 好男人视频免费观看在线| 国产男女超爽视频在线观看| 国内精品宾馆在线| 日韩人妻高清精品专区| a级毛色黄片| 亚洲自拍偷在线| 91精品国产九色| 97人妻精品一区二区三区麻豆| 男人爽女人下面视频在线观看| 欧美激情久久久久久爽电影| 97热精品久久久久久| 日本熟妇午夜| 欧美性猛交╳xxx乱大交人| 波多野结衣巨乳人妻| 国产高清国产精品国产三级 | 国产69精品久久久久777片| 深夜a级毛片| 中文字幕亚洲精品专区| 99久久精品热视频| 真实男女啪啪啪动态图| 久久久久久国产a免费观看| kizo精华| 两个人视频免费观看高清| 久久久久久久久久久免费av| 一本久久精品| 两个人的视频大全免费| 非洲黑人性xxxx精品又粗又长| 婷婷色麻豆天堂久久| 亚洲精品成人av观看孕妇| 亚洲精品,欧美精品| 国产精品久久久久久久久免| 久久99热这里只有精品18| 观看免费一级毛片| 男女边吃奶边做爰视频| 中国美白少妇内射xxxbb| 国产黄片视频在线免费观看| 色网站视频免费| 美女cb高潮喷水在线观看| 日本黄大片高清| 国产精品久久视频播放| 国产成人午夜福利电影在线观看| 美女黄网站色视频| 久久久久性生活片| av天堂中文字幕网| 亚洲av成人av| 大又大粗又爽又黄少妇毛片口| 欧美另类一区| 视频中文字幕在线观看| 亚洲无线观看免费| 欧美激情在线99| 欧美最新免费一区二区三区| 成年人午夜在线观看视频 | 日韩成人av中文字幕在线观看| 大香蕉久久网| 亚洲最大成人手机在线| 人妻夜夜爽99麻豆av| 中文乱码字字幕精品一区二区三区 | 免费观看无遮挡的男女| 欧美日韩一区二区视频在线观看视频在线 | 免费大片18禁| 成人亚洲精品一区在线观看 | 日韩精品青青久久久久久| 97精品久久久久久久久久精品| 亚洲自拍偷在线| 高清毛片免费看| 一级毛片 在线播放| 亚洲三级黄色毛片| 欧美性感艳星| 老司机影院毛片| 国产片特级美女逼逼视频| 日日啪夜夜撸| 成人亚洲精品av一区二区| 我的女老师完整版在线观看| 国产淫语在线视频| 99九九线精品视频在线观看视频| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区在线观看99 | av福利片在线观看| 成人av在线播放网站| 亚洲欧美精品专区久久| 国产乱人偷精品视频| 精品午夜福利在线看| kizo精华| 美女内射精品一级片tv| 男的添女的下面高潮视频| 国产成人freesex在线| 床上黄色一级片| 夜夜爽夜夜爽视频| 欧美性猛交╳xxx乱大交人| 大香蕉久久网| 国产片特级美女逼逼视频| 日本与韩国留学比较| 99久久人妻综合| 九色成人免费人妻av| 中国美白少妇内射xxxbb| 亚洲欧美日韩卡通动漫| av.在线天堂| 少妇熟女欧美另类| 亚洲国产av新网站| 国产激情偷乱视频一区二区| 欧美zozozo另类| 永久网站在线| 女人被狂操c到高潮| 国产成人91sexporn| 可以在线观看毛片的网站| 日本av手机在线免费观看| 美女脱内裤让男人舔精品视频| 最近2019中文字幕mv第一页| 久久精品久久久久久噜噜老黄| 十八禁国产超污无遮挡网站| 噜噜噜噜噜久久久久久91| 三级经典国产精品| 国产精品一区www在线观看| 97人妻精品一区二区三区麻豆| 丝袜喷水一区| 久久久精品欧美日韩精品| 亚洲成色77777| 国产人妻一区二区三区在| 哪个播放器可以免费观看大片| 最近最新中文字幕大全电影3| 97在线视频观看| 国产精品一区二区三区四区免费观看| 国产淫片久久久久久久久| 97人妻精品一区二区三区麻豆| 十八禁国产超污无遮挡网站| 免费看美女性在线毛片视频| 最近中文字幕2019免费版| 国产中年淑女户外野战色| 草草在线视频免费看| 啦啦啦韩国在线观看视频| 最近中文字幕高清免费大全6| 精品熟女少妇av免费看| 亚洲成人一二三区av| 七月丁香在线播放| 午夜免费男女啪啪视频观看| 亚洲高清免费不卡视频| 亚洲精品视频女| av线在线观看网站| 男女边吃奶边做爰视频| 国产黄片美女视频| 97精品久久久久久久久久精品| 热99在线观看视频| 亚洲人成网站在线观看播放| 91精品伊人久久大香线蕉| 天天躁夜夜躁狠狠久久av| 特大巨黑吊av在线直播| 日韩av免费高清视频| 日韩国内少妇激情av| 久久久色成人| 国产一区二区在线观看日韩| 日韩一区二区视频免费看| 色播亚洲综合网| 国产永久视频网站| 久久精品国产亚洲av天美| 国产精品三级大全| 色综合色国产| 精品人妻熟女av久视频| 国产伦精品一区二区三区视频9| 国产在线男女| videossex国产| 国产av在哪里看| av黄色大香蕉| 乱人视频在线观看| 精品国产露脸久久av麻豆 | 国产精品综合久久久久久久免费| 国产精品蜜桃在线观看| xxx大片免费视频| 亚洲最大成人中文| 老司机影院成人| 日本色播在线视频| 夜夜爽夜夜爽视频| 日韩av在线大香蕉| 午夜日本视频在线| 乱码一卡2卡4卡精品| 超碰av人人做人人爽久久| 国产一级毛片七仙女欲春2| 欧美 日韩 精品 国产| 亚洲av成人av| 三级经典国产精品| 人妻系列 视频| 亚州av有码| 精品少妇黑人巨大在线播放| 在线免费十八禁| 国内揄拍国产精品人妻在线| 国产精品人妻久久久久久| 亚洲va在线va天堂va国产| 日韩欧美精品v在线| 青春草视频在线免费观看| 久久韩国三级中文字幕| 男女边摸边吃奶| 国产综合懂色| 亚洲av一区综合| 大话2 男鬼变身卡| 蜜桃久久精品国产亚洲av| 搡老乐熟女国产| 国产午夜福利久久久久久| 午夜福利成人在线免费观看| 91久久精品国产一区二区三区| 国产有黄有色有爽视频| 极品教师在线视频| 久久精品综合一区二区三区| 国产男女超爽视频在线观看| 亚洲成人一二三区av| 国产黄a三级三级三级人| 免费av不卡在线播放| 国产精品国产三级国产专区5o| 亚洲自偷自拍三级| 黄色配什么色好看| 亚洲精品成人av观看孕妇| 久久久久久久久久黄片| 国产一区二区三区综合在线观看 | 少妇的逼好多水| 欧美日本视频| 国产精品久久久久久av不卡| 成年女人看的毛片在线观看| 国产视频内射| 精品久久久久久电影网| 精品亚洲乱码少妇综合久久| 五月玫瑰六月丁香| 日韩欧美一区视频在线观看 | 色综合站精品国产| 精品久久久久久久人妻蜜臀av| 国产综合精华液| 在线观看人妻少妇| 亚洲内射少妇av| 中文字幕av在线有码专区| 18禁在线播放成人免费| 一个人看视频在线观看www免费| 最近中文字幕高清免费大全6| 成人亚洲欧美一区二区av|