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

    火災(zāi)區(qū)域林地植被變化分析—以大興安嶺和黃石公園地區(qū)為例

    2019-06-10 14:07:16達,張虎,何
    關(guān)鍵詞:黃石公園跡地大興安嶺

    錢 達,張 虎,何 龍

    (天津師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院,天津 300387)

    森林作為地球上最大的陸地生態(tài)系統(tǒng),是地球生物圈中重要的組成部分之一,在維系整個地球的生態(tài)平衡中起著至關(guān)重要的作用,同時也是人類賴以生存和發(fā)展的資源和環(huán)境.自地球出現(xiàn)森林以來,森林火災(zāi)時常發(fā)生.全世界每年平均發(fā)生森林火災(zāi)20 多萬次,燒毀森林面積約占全世界森林總面積的1%以上.中國平均每年發(fā)生森林火災(zāi)約1 萬多次,燒毀森林幾十萬至上百萬公頃,約占全國森林面積的5%~8%.火災(zāi)會給森林帶來嚴重的惡果,最直觀的危害就是燒毀大面積的植被,降低整個區(qū)域的植被覆蓋率,同時土壤中的有機質(zhì)遭到破壞容易引起土地沙化和水土流失,大量動植物和微生物被燒死,其生存環(huán)境亦遭到破壞,打破了整個森林的生態(tài)平衡,對人民生命財產(chǎn)安全造成威脅.

    由于森林火災(zāi)的毀壞性極強,近年來,森林火災(zāi)問題受到研究人員的廣泛關(guān)注,研究方向主要有火災(zāi)跡地的提取[1-3]和火災(zāi)后植被恢復(fù)等[4-6].火災(zāi)跡地的提取多圍繞燃燒進行,小范圍的火災(zāi)跡地可以利用燃燒面積指數(shù)BAI(burn area index,BAI)通過增強過火后圖像上的木炭信息進行提取[1];較大范圍的火災(zāi)跡地則需要利用歸一化燃燒指數(shù)NBR(normalized burn ratio,NBR)進行提取[2].NBR 可以將影像中的水體區(qū)域和火燒跡地區(qū)別開,并具有很好的整體效果,精度達到89%[3].但將NBR 作為評價火災(zāi)恢復(fù)情況的指數(shù)只能得到火燒跡地的燒毀面積和恢復(fù)面積,無法準確判斷出火燒跡地植被的恢復(fù)程度和燒毀程度.有關(guān)火災(zāi)后植被恢復(fù)的研究主要與植被歸一化指數(shù)NDVI(normalized difference vegetation index,NDVI) 相關(guān),如采用NDVI 像元二分線性模型估算法分析森林火災(zāi)前后火燒跡地植被覆蓋度時空動態(tài)變化[4],利用NDVI月度值監(jiān)測整個森林火災(zāi)前后火災(zāi)期間的植被行為[5],以及將混合光譜分析SMA(Spectral mixture analysis)與NDVI 的結(jié)果進行比較,以評估植被恢復(fù)[6]等.鑒于近年來火災(zāi)跡地的提取和火災(zāi)后植被恢復(fù)的研究情況,采用NBR 提取火燒跡地,然后利用火災(zāi)前后NDVI差值法評價火燒跡地植被的恢復(fù)程度以及燒毀程度,可以更加準確地獲取火燒跡地范圍,精確評估火燒跡地植被的恢復(fù)程度和燒毀程度,很好地彌補NBR 在評價恢復(fù)程度和燒毀程度方面的不足[7].

    本研究以1987年黑龍江省大興安嶺地區(qū)和黃石公園森林火災(zāi)區(qū)為研究對象,以時間序列的Landsat-TM 地表反射率數(shù)據(jù)為基礎(chǔ),利用歸一化燃燒指數(shù)(NBR)提取研究區(qū),對研究區(qū)植被指數(shù)的變化情況進行對比分析,探究大興安嶺森林火災(zāi)和黃石公園的燒毀程度以及災(zāi)后植被的恢復(fù)程度,以期為該地區(qū)林火管理提供依據(jù),同時也為火災(zāi)后林區(qū)的恢復(fù)和保護工作提供參考和理論支持.

    1 研究區(qū)域與數(shù)據(jù)

    本文研究區(qū)域為大興安嶺1987年火災(zāi)燒毀區(qū)域和1988年黃石國家公園火災(zāi)燒毀區(qū)域.圖1 為研究區(qū)及其周邊的標準假彩色影像,其中黑色區(qū)域即為火災(zāi)燒毀區(qū)域.

    圖1 大興安嶺火災(zāi)和黃石國家公園火災(zāi)區(qū)域及其周邊標準假彩色影像Fig.1 Standard false-color satellite image of the conflagration area and surrounding area in Da Hinggan Mountains and Yellowstone National Park

    圖1(a)中大興安嶺位于黑龍江省大興安嶺地區(qū)、內(nèi)蒙古自治區(qū)呼倫貝爾市東北部,是內(nèi)蒙古高原與松遼平原的分水嶺,屬大陸性季風氣候,年均溫度為-2~5°C,無霜期90~100 d,年降水量為400~500 mm.1987年5月6日大興安嶺北部發(fā)生特大森林火災(zāi),火災(zāi)范圍涵蓋了塔河、西林吉、阿爾木和圖強共4 個林業(yè)局,過火森林面積達5 600 km2,于同年6月2日被徹底撲滅.圖1(b)中黃石國家公園地處落基山脈,位于美國中西部懷俄明州的西北角,并向西北方向延伸至愛達荷州和蒙大拿州,東西距離87 km,南北寬101 km,面積達8 956 km2,屬于高原山地氣候.黃石公園的溫度和降水量變化很大,夏季期間白天最高溫度為21 ~27 ℃,冬天最高溫度為-20 ~-5 ℃.年均降水量從猛犸象溫泉附近的380 mm 到該公園西南部的2 000 mm 不等.1988年,黃石公園多個獨立的小山火先后在風勢加劇和干旱的作用下失去了控制,火勢迅速蔓延,最終形成了一場燃燒了幾個月的大火災(zāi),直到同年11月才被完全撲滅.本研究以搭載于Landsat-5 衛(wèi)星上的TM 數(shù)據(jù)為基礎(chǔ),使用美國國家航空航天局NASA 的大陸區(qū)域森林擾動系統(tǒng)LEDAPS 處理系統(tǒng)對TM 圖像進行大氣校正.LEDAPS 以6S 大氣輻射傳輸模型理論為基礎(chǔ),采用大量的實時大氣觀測數(shù)據(jù),精度較高.大興安嶺研究區(qū)對應(yīng)的Landsat 數(shù)據(jù)分幅號為h121v023 和h122v023,時間跨度為1986—1999年.黃石公園研究區(qū)對應(yīng)的Landsat 數(shù)據(jù)分幅號為h038v029,時間跨度為1988—2000年.地表反射率數(shù)據(jù)經(jīng)過波段合成、 影像拼接和裁剪后,得到覆蓋研究區(qū)的遙感影像,空間分辨率為30 m.

    受云層影響和衛(wèi)星觀測能力的限制,研究區(qū)數(shù)據(jù)可選擇性有限,難以保證研究區(qū)所有數(shù)據(jù)均為同一月份.同時,受季節(jié)氣候和植物生長規(guī)律的影響,導(dǎo)致不同季節(jié)植被覆蓋度不同,所以在研究前需對不同月份NDVI 數(shù)據(jù)的一致性進行驗證.分別在大興安嶺和黃石公園火燒跡地附近未受火災(zāi)影響的林區(qū)隨機選擇5 個50×50 像元大小的控制點,如圖1 中黃色編號所指位置,通過對比控制點處不同月份NDVI 均值與方差的變化情況檢驗TM 數(shù)據(jù)的一致性,結(jié)果如圖2 所示.

    圖2 大興安嶺和黃石公園未受火災(zāi)影響區(qū)域控制點5—10月NDVI 均值與方差變化圖Fig.2 Variation of NDVI mean value and variance of the control points in the areas which is not affected by fire in the Da Hinggan Mountains from May to October

    由圖2 可以看出,5、9 和10月份的NDVI 值較低,采取此時間段內(nèi)的數(shù)據(jù)會對研究結(jié)果的分析造成較大的不確定性;6~8月份的NDVI 值較高,基本維持在一個穩(wěn)定的數(shù)值范圍內(nèi),相互差值最大為0.03左右,數(shù)據(jù)間有較高的一致性.因此,本研究選用6~8月份的TM 地表反射率數(shù)據(jù).

    2 研究方法

    本文所選的研究區(qū)域?qū)儆诖蠓秶馂?zāi)區(qū)域,歸一化燃燒指數(shù)NBR 對燒毀面積超過200 公頃的大范圍火災(zāi)區(qū)域具有增強效果[8],因此本文采用NBR 提取火燒跡地,計算公式為

    式(1)中:ρNIR和ρWEIR和分別為近紅外波段和短波紅外波段的反射率.火災(zāi)區(qū)域植被遭到破壞,光譜不再具有健康綠色植被的特征,近紅外波段反射率下降,短波紅外波段由于缺少植被冠層覆蓋,濕度降低導(dǎo)致反射率上升[9];而健康綠色植被在短波紅外波段因為葉片水分的吸收作用常形成吸收谷.過火區(qū)域和未過火區(qū)域在近紅外和短波紅外2 個波段形成明顯差異,從而突出了過火區(qū)域[10].

    由于火燒跡地影像中裸地、云層和水體的歸一化燃燒指數(shù)值相近,需要手動確定歸一化燃燒指數(shù)的閾值范圍以進行區(qū)分,經(jīng)過多次修改閾值范圍,通過對比得到一個最佳閾值范圍[3].本研究最終選取的閾值范圍為-0.58 ~-0.1,還有極少的裸地以及云層沒有篩選出來,需要通過目視解譯進行去除.

    提取出火燒跡地后,需要通過研究區(qū)多年的地表反射率計算出與植被覆蓋度有極強相關(guān)性的歸一化植被指數(shù)NDVI,用于評價研究區(qū)多年植被覆蓋度的變化,NDVI 的計算公式為:

    式(2)中:ρNIR和ρR分別為近紅外波段和紅波段的反射率.正常植物在紅波段和近紅外波段的反射差異很大,植物光合作用吸收紅光,反射出近紅外光,因此依據(jù)吸收紅光的量可以反映植物的長勢,NDVI 值越大,植物的長勢越好.同時,由NDVI 值也可以獲知植物的數(shù)量.

    本研究利用NBR 提取的火燒跡地范圍裁剪不同年份大興安嶺和黃石公園的地表反射率數(shù)據(jù)后得到1986—1999年大興安嶺火災(zāi)區(qū)域和1988—1999年黃石公園火災(zāi)區(qū)域的地表反射率數(shù)據(jù),分別計算得到不同年份大興安嶺火災(zāi)區(qū)域和黃石公園區(qū)域的NDVI,通過不同年份的NDVI 差值評價火災(zāi)區(qū)域的受災(zāi)情況以及火災(zāi)區(qū)域和重度火災(zāi)區(qū)域災(zāi)后多年的恢復(fù)情況.同時,在火災(zāi)區(qū)域附近未受影響的林區(qū)隨機選擇5 個50×50 像元大小的控制點,通過對比控制點處不同年份NDVI 的變化情況監(jiān)測未燒毀區(qū)域植被的變化情況,用于火災(zāi)區(qū)域植被的變化情況的對照實驗.

    3 結(jié)果與討論

    3.1 火災(zāi)區(qū)域的提取及其NDVI分布圖

    利用歸一化燃燒指數(shù)(NBR)提取出的大興安嶺和黃石公園的火燒跡地范圍圖3 所示.

    圖3 大興安嶺和黃石公園火災(zāi)區(qū)域范圍的提取結(jié)果Fig.3 Extraction result of conflagration area in Da Hinggan Mountains and Yellowstone National Park

    圖3 中藍色區(qū)域為大興安嶺和黃石公園的火災(zāi)區(qū)域,即本文的研究區(qū).其中,大興安嶺火災(zāi)區(qū)域的面積為6 695.53 km2,黃石公園火災(zāi)區(qū)域的面積為2 187.95 km2.

    對大興安嶺研究區(qū)和黃石公園研究區(qū)的地表反射率數(shù)據(jù)進行計算后,分別得到二者不同年份的NDVI 結(jié)果圖,并從2 個研究區(qū)分別選取6 幅NDVI結(jié)果圖,如圖4 和圖5 所示.

    由圖4 可以看出,與1986年相比,1987年大興安嶺研究區(qū)的NDVI 大幅度降低,1988年后開始緩慢增加,至1996年,研究區(qū)的NDVI 已經(jīng)恢復(fù)到受災(zāi)前的水平,1999年研究區(qū)的NDVI 超過1986年的NDVI.

    圖4 大興安嶺研究區(qū)時間序列NDVI 結(jié)果圖Fig.4 NDVI of research area in Da Hinggan Mountains with the time series

    圖5 黃石公園研究區(qū)時間序列NDVI 結(jié)果圖Fig.5 NDVI of research area in Yellowstone National Park with the time series

    由圖5 可以看出,與1988 相比,1989年黃石公園研究區(qū)的NDVI 明顯下降,1990年后開始緩慢升高,1998年仍然沒有恢復(fù)到受災(zāi)前的水平,直到2000年NDVI 才恢復(fù)到受災(zāi)前的水平,2000年部分區(qū)域的NDVI 略高于受災(zāi)前的NDVI.2 個研究區(qū)的NDVI 變化過程都是先降低然后緩慢增加,但大興安嶺研究區(qū)的恢復(fù)速度相對較快.

    3.2 研究區(qū)受災(zāi)情況評估

    將研究區(qū)火災(zāi)前后的NDVI 進行差值計算,根據(jù)NDVI 的差值劃分林區(qū)燒毀程度,分級標準如表1 所示.

    表1 燒毀程度分級表Tab.1 Classification of burned degree

    對大興安嶺(1986年和1987年)和黃石公園(1988年和1989年)火災(zāi)前后的NDVI 進行差值分級后,得到2 個研究區(qū)的燒毀程度,結(jié)果如圖6 所示.由圖6(a)可以看出,大興安嶺研究區(qū)中輕度燒毀區(qū)域的面積為270.62 km2,中度燒毀區(qū)域的面積為1 628.96 km2,重度燒毀區(qū)域的面積為4 738.55 km2.由圖6(b)可知,黃石公園研究區(qū)中輕度燒毀區(qū)域的面積為204.78 km2,中度燒毀區(qū)域的面積為1 117.83 km2,重度燒毀區(qū)域的面積為772.72 km2.由圖6 結(jié)果可知,大興安嶺火災(zāi)在燒毀范圍和燒毀程度上都超過了黃石公園火災(zāi),受災(zāi)情況更加嚴重.

    3.3 研究區(qū)災(zāi)后恢復(fù)情況評價

    對研究區(qū)火災(zāi)前后的NDVI 差值進行分級,得到恢復(fù)程度分級標注,結(jié)果如表2 所示.

    圖7 為大興安嶺和黃石公園研究區(qū)火災(zāi)區(qū)域、重度火災(zāi)區(qū)域和未燒毀區(qū)域NDVI 隨時間的變化情況.由圖7(a)可以看出,大興安嶺火災(zāi)區(qū)域和重度火災(zāi)區(qū)域的NDVI 均從1987年開始降低,然后緩慢上升趨于平穩(wěn),未燒毀區(qū)域一直保持平穩(wěn).大興安嶺研究區(qū)火災(zāi)區(qū)域1986年NDVI 為0.69,植被生長狀況良好.1987年受火災(zāi)影響,火災(zāi)區(qū)域NDVI 均值驟降至0.23,可見火災(zāi)對林區(qū)影響程度非常大.災(zāi)后相關(guān)部門對不同燒毀程度的區(qū)域?qū)嵤┫鄳?yīng)的恢復(fù)措施,使得1988年后火災(zāi)區(qū)域NDVI 迅速回升至0.45,相比1987年火災(zāi)區(qū)域NDVI 提升很大,植被恢復(fù)速度相當可觀.1992年,火災(zāi)區(qū)域NDVI 增長至0.68,已經(jīng)接近1986年火災(zāi)區(qū)域NDVI 值,隨后幾年火災(zāi)區(qū)域NDVI 值均保持在0.7 以上.由圖7(b)可以看出,黃石公園火災(zāi)區(qū)域和重度火災(zāi)區(qū)域的NDVI 均在1989年開始降低,然后緩慢上升趨于平穩(wěn),黃石公園未燒毀區(qū)域NDVI 略有波動,但基本保持平緩.黃石公園研究區(qū)火災(zāi)區(qū)域1988年NDVI 為0.51,植被生長狀況良好.1989年受火災(zāi)影響,火災(zāi)區(qū)域NDVI 均值驟降至0.18,可見火災(zāi)對研究區(qū)域的森林造成極大的破壞.1990年火災(zāi)區(qū)域NDVI 恢復(fù)至0.24,恢復(fù)速度較為緩慢,1990—1994年火災(zāi)區(qū)域NDVI 在0.35 左右,相對于大興安嶺火災(zāi)區(qū)域來說恢復(fù)較為緩慢,1995—1998年恢復(fù)到0.45 左右,相對1989年已經(jīng)恢復(fù)較多,直到2000年才達到燒毀前的植被覆蓋度.

    圖6 大興安嶺研究區(qū)和黃石公園研究區(qū)火災(zāi)燒毀程度的分布Fig.6 Degree distribulion of fire damage of research area in Da Hinggan Mountains and Yellowstone National Park

    表2 恢復(fù)程度分級表Tab.2 Classification of recovery degree

    圖7 大興安嶺和黃石公園火災(zāi)區(qū)域、未受火災(zāi)影響區(qū)域以及重度火災(zāi)區(qū)域的多年NDVI 均值變化趨勢圖Fig.7 Tendency chart of NDVI mean over many years of burned area,unburned area and severely burned area in Da Hinggan Mountains and Yellowstone National Park

    火災(zāi)區(qū)域的少部分年份的NDVI 變化趨勢有所波動,主要原因有:①遵循云層比小于10%的原則篩選TM 數(shù)據(jù),在研究過程中盡量去除薄云,但厚云的灰度值與火燒跡相似,區(qū)別度不大[5],造成NDVI 略低,如大興安嶺火災(zāi)區(qū)域1986年TM 影像中有少量厚云(厚云的NDVI 值極小),造成1986年火災(zāi)區(qū)域的NDVI略低于其他2 個區(qū)域;②由于數(shù)據(jù)可選擇性有限,很難選擇不同年份同一月份的數(shù)據(jù),導(dǎo)致數(shù)據(jù)的連續(xù)性受到影響,造成局部結(jié)果產(chǎn)生些許誤差.同時本研究對不同年份植被覆蓋度進行比較,其中不同年份氣候因素不同造成當年植被生長狀況存在一定差異,對研究結(jié)果也會產(chǎn)生一定的影響,如黃石公園火災(zāi)區(qū)域1994年選用的是8月末的數(shù)據(jù),而1993年選用的是7月份數(shù)據(jù),導(dǎo)致1994年未燒毀區(qū)域NDVI 略高于1993年未燒毀區(qū)域NDVI;(3)本研究僅考慮了NDVI信息,沒有考慮樹種的變化對NDVI 的影響,同樣會對研究結(jié)果帶來一定不確定性.從大興安嶺火災(zāi)區(qū)域來看,火災(zāi)后落葉松林首先由雜草群落和雜草灌木群落所代替,隨后白樺幼苗大量發(fā)生,逐漸形成白樺林[10].白樺屬于喬木,葉面較大,而落葉松的葉面是針葉,白樺林代替落葉松林,會使該區(qū)域NDVI 上升.

    大興安嶺研究區(qū)1988年和1996年的NDVI 分別與1986年的NDVI 進行差值計算,分級后得到大興安嶺研究區(qū)恢復(fù)程度圖,結(jié)果如圖8 所示.

    圖8 大興安嶺研究區(qū)1988年和1996年植被恢復(fù)程度分布圖Fig.8 Distribution diagram of vegetation restoration degree of the research area in Da Hinggan Mountains in 1988 and 1996

    1988年為燒毀后第1年,火災(zāi)區(qū)域內(nèi)植被覆蓋度略有改善,與過火前的植被覆蓋度相比相差甚遠.完全恢復(fù)區(qū)域的面積為207.99 km2,基本恢復(fù)區(qū)域的面積為2008.18km2,部分恢復(fù)區(qū)域面積為3847.16km2,較少恢復(fù)區(qū)域面積為572.61 km2.完全恢復(fù)區(qū)域面積只占整個火災(zāi)區(qū)域面積的3.13%,而部分恢復(fù)區(qū)域面積和基本恢復(fù)區(qū)域面積占整個火災(zāi)區(qū)域面積的88.24%,可見整個區(qū)域的植被恢復(fù)情況不容樂觀.1996年,即燒毀后的第9年,火災(zāi)區(qū)域內(nèi)基本看不到部分恢復(fù)區(qū)域和較小恢復(fù)區(qū)域,大部分區(qū)域植被覆蓋情況已經(jīng)達到過火前的水平.完全恢復(fù)區(qū)域的面積為4 844.56 km2,基本恢復(fù)區(qū)域面積為1 769.91 km2,部分恢復(fù)區(qū)域和較少恢復(fù)區(qū)域面積僅為17.32 km2和5.93 km2.完全恢復(fù)區(qū)域面積和基本恢復(fù)面積分別占整個火災(zāi)區(qū)域面積的72.99%和26.66%,部分恢復(fù)區(qū)域和較小恢復(fù)區(qū)域基本可以忽略不計,總體上整個火災(zāi)區(qū)域的植被覆蓋度已經(jīng)和過火前相差無幾,火災(zāi)區(qū)域內(nèi)植被覆蓋情況基本恢復(fù)到火災(zāi)發(fā)生前的水平.

    黃石公園研究區(qū)1990年和1998年的NDVI 分別與1988年的NDVI 進行差值計算,分級后得到黃石公園研究區(qū)恢復(fù)程度圖,結(jié)果如圖9 所示.

    圖9 黃石公園研究區(qū)1990年和1998年植被恢復(fù)程度分布圖Fig.9 Distribution diagram of vegetation restoration degree of the research area in the Yellowstone National Park in 1990 and 1998

    1990年即燒毀后的第2年,火災(zāi)區(qū)域內(nèi)植被覆蓋度雖略有改善,但與過火前的植被覆蓋度相差甚遠.完全恢復(fù)區(qū)域的面積為176.01 km2,基本恢復(fù)區(qū)域面積為359.53 km2,部分恢復(fù)區(qū)域面積為1184.21 km2,較少恢復(fù)區(qū)域面積為465.30 km2.完全恢復(fù)區(qū)域面積只占整個火災(zāi)區(qū)域面積的8.1%,而部分恢復(fù)區(qū)域面積和基本恢復(fù)區(qū)域面積占整個火災(zāi)區(qū)域面積的70.6%,較少恢復(fù)區(qū)域面積占整個火災(zāi)區(qū)域的21.3%,可見整個區(qū)域的植被完全恢復(fù)的面積較小,大部分區(qū)域植被仍舊在恢復(fù)當中.1998年,燒毀后的第10年,火災(zāi)區(qū)域內(nèi)基本看不到較小恢復(fù)區(qū)域,部分區(qū)域植被覆蓋情況已經(jīng)達到過火前的水平.完全恢復(fù)區(qū)域的面積為469.64 km2,基本恢復(fù)區(qū)域面積為1 295.75 km2,部分恢復(fù)區(qū)域面積為419.90 km2,較少恢復(fù)區(qū)域面積為2.64 km2.完全恢復(fù)區(qū)域面積占據(jù)了整個火災(zāi)區(qū)域面積的36%,基本恢復(fù)區(qū)域面積占整個火災(zāi)區(qū)域面積的57%,部分恢復(fù)區(qū)域面積占整個火災(zāi)區(qū)域面積的7%,較小恢復(fù)區(qū)域基本可以忽略不計,總體上整個火災(zāi)區(qū)域的植被覆蓋度還沒有恢復(fù)到過火前的水平,火災(zāi)區(qū)域內(nèi)仍有一大半植被還在恢復(fù)過程中.

    1987年大火后,大興安嶺的火燒跡地根據(jù)不同的火燒強度采取了不同的管理措施.在輕度燒毀的火災(zāi)區(qū)域采取自然更新和人工促進自然更新為主的管理措施,在重度火災(zāi)燒毀區(qū)域采取人工更新的管理措施,主要以興安落葉松為主[11].1988年大火后黃石公園采取的措施是自然發(fā)生的火災(zāi)只要控制住火勢,沒有威脅人類安全,就讓火苗自生自滅,甚至要求工作人員定期把死亡的樹木燒毀,人為地制造一些隔離帶[12].同是大范圍的火災(zāi),在火災(zāi)后采取不同的管理措施,取得的效果也不相同.1987年大興安嶺火災(zāi)燒毀面積為6 638.13 km2,火災(zāi)區(qū)域的NDVI 下降了0.46,其中重度火災(zāi)區(qū)域占整個火災(zāi)區(qū)域面積的71.38%.1988年黃石公園火災(zāi)燒毀面積為2 095.33 km2,火災(zāi)區(qū)域的NDVI 下降了0.28,其中重度火災(zāi)區(qū)域占整個火災(zāi)區(qū)域面積的36.88%.上述分析表明,大興安嶺火災(zāi)的燒毀范圍和燒毀程度遠超過黃石公園,此外,大興安嶺火災(zāi)區(qū)域內(nèi)植被被燒毀的數(shù)量也遠大于黃石公園火災(zāi)區(qū)域內(nèi)植被燒毀的數(shù)量.從恢復(fù)情況看,在燒毀后的第1年大興安嶺火災(zāi)區(qū)域NDVI 上升了0.22,在燒毀后的第5年NDVI 已經(jīng)接近火災(zāi)前的NDVI,第9年時已經(jīng)超過了火災(zāi)前的NDVI;而黃石公園火災(zāi)區(qū)域在燒毀后的第2年里NDVI 上升了0.06,在燒毀后的第6年NDVI 為0.35,距離燒毀前0.46 的NDVI 還是存在一定差距,在燒毀后的第10年仍舊未達到火災(zāi)前的NDVI.由此可知,大興安嶺火災(zāi)區(qū)域的植被恢復(fù)速度和恢復(fù)程度均超過黃石公園火災(zāi)區(qū)域.短期來看,相對于黃石公園以自然更新恢復(fù)的管理措施,大興安嶺以人工促進自然更新和人工更新恢復(fù)的管理措施對火災(zāi)后森林植被的恢復(fù)具有更好的效果.

    4 結(jié)論

    本研究選擇大興安嶺和黃石公園兩地森林火災(zāi)區(qū)域為研究對象,利用NBR 提取火災(zāi)區(qū)域的范圍,以時間序列TM 地表反射率數(shù)據(jù)為基礎(chǔ),通過比較不同時間范圍內(nèi)火災(zāi)區(qū)域附近NDVI 數(shù)據(jù)的變化情況,評價數(shù)據(jù)的一致性,進而根據(jù)火災(zāi)前后火災(zāi)區(qū)域內(nèi)NDVI 的變化情況,評價研究區(qū)的受災(zāi)程度和恢復(fù)情況.研究結(jié)果表明:

    (1)大火對大興安嶺和黃石公園林區(qū)造成了巨大傷害.大興安嶺基于NBR 提取的火災(zāi)區(qū)域總面積為6 695.53 km2,其中輕度燒毀區(qū)域面積為270.62 km2,中度燒毀區(qū)域面積為1 628.96 km2,重度燒毀區(qū)域面積為4 738.55 km2,研究區(qū)NDVI 受火災(zāi)影響極為嚴重,NDVI 下降了66.32%.黃石公園基于NBR 提取的火災(zāi)區(qū)域面積為2 187.95 km2,其中,輕度燒毀區(qū)域面積為204.78 km2,中度燒毀區(qū)域面積為1117.83 km2,重度燒毀區(qū)域面積為772.72 km2,研究區(qū)NDVI 下降了64.68%.大興安嶺的植被覆蓋度整體下降幅度、火災(zāi)范圍和燒毀程度均超過黃石公園.

    (2)大興安嶺火災(zāi)區(qū)域災(zāi)后NDVI 恢復(fù)速度較快,1988年NDVI 比1987年NDVI 增長了0.22,恢復(fù)到1986年的65.65%,1993年NDVI 值超過1986年NDVI 值,隨后幾年增長率較小,趨于平衡,至1999年NDVI 值較1986年增長了0.09.黃石公園火災(zāi)區(qū)域災(zāi)后NDVI 恢復(fù)速度較為緩慢,1990年NDVI 比1989年NDVI 增長了0.06,恢復(fù)到1988年的46.42%.1990年后NDVI 緩慢恢復(fù),仍然沒有達到過火前的NDVI水平,1998年達到過火前NDVI 的84.75%,直到2000年NDVI 才恢復(fù)到火災(zāi)前的水平.

    (3) 火災(zāi)發(fā)生后火災(zāi)區(qū)域的NDVI 驟降,對火災(zāi)區(qū)域的植被造成了極大的破壞.在大興安嶺火災(zāi)輕度燒毀的火災(zāi)區(qū)域采取自然更新和人工促進自然更新為主的管理措施,在重度火災(zāi)燒毀區(qū)域采取人工更新的管理措施.人工促進恢復(fù)可以加速森林群落的更新和受災(zāi)區(qū)域森林生態(tài)系統(tǒng)的恢復(fù),但是需要投入大量的人力物力來維持.黃石公園火災(zāi)區(qū)域采取任其燒毀的管理措施.自然恢復(fù)可以減少經(jīng)濟投入,但恢復(fù)速度略慢,同時森林群落演替的時間過長,短時期內(nèi)森林群落主要被優(yōu)勢物種所占據(jù).大興安嶺和黃石公園火災(zāi)區(qū)域在2 種截然不同的管理措施下,經(jīng)過一段時間的修復(fù),基本都恢復(fù)到火災(zāi)前的植被覆蓋情況.大興安嶺火災(zāi)區(qū)域植被恢復(fù)速度快于黃石公園的植被恢復(fù)速度,說明火災(zāi)區(qū)域的森林植被依靠人工促進恢復(fù)輔助自然恢復(fù)要比自然恢復(fù)更早恢復(fù)到火災(zāi)前的植被覆蓋水平.在以后森林火災(zāi)的恢復(fù)中,可以根據(jù)火災(zāi)燒毀程度等級采取對應(yīng)措施.輕度火災(zāi)區(qū)域森林恢復(fù)可以通過自然恢復(fù),在節(jié)約成本的同時也不影響植被的恢復(fù);中度火災(zāi)區(qū)域森林恢復(fù)可以采取人工促進恢復(fù)輔助自然恢復(fù),以期在經(jīng)濟效益和生態(tài)效益中尋求平衡; 重度火災(zāi)區(qū)域森林恢復(fù)采取人工促進恢復(fù),以期加快火災(zāi)區(qū)域的恢復(fù)進度.

    猜你喜歡
    黃石公園跡地大興安嶺
    黃石公園奇遇
    封二、封三:寒蔥嶺戰(zhàn)跡地 (伏擊戰(zhàn)) 新考
    大興安嶺不會忘記你
    草原歌聲(2019年3期)2019-10-17 02:20:06
    我國科學(xué)家在國際上首發(fā)30 m分辨率全球火燒跡地數(shù)據(jù)產(chǎn)品
    大興安嶺的鐵道兵
    特別健康(2018年4期)2018-07-03 00:38:22
    2架增雨飛機為大興安嶺送雨
    伙伴(2018年7期)2018-05-14 14:19:33
    火燒跡地造林技術(shù)
    大興安嶺四季(四首)
    北極光(2017年10期)2017-11-13 23:09:01
    基于Landsat影像提取大興安嶺區(qū)域火燒跡地★
    山西建筑(2017年15期)2017-07-01 22:03:22
    美麗的黃石公園
    kizo精华| 亚洲中文av在线| 久久av网站| 在线精品无人区一区二区三| 熟女av电影| 午夜免费男女啪啪视频观看| 亚洲丝袜综合中文字幕| 久久久精品免费免费高清| 热99久久久久精品小说推荐| 国产探花极品一区二区| 日产精品乱码卡一卡2卡三| 国产欧美另类精品又又久久亚洲欧美| 日本欧美国产在线视频| av黄色大香蕉| 日本黄大片高清| 七月丁香在线播放| 99久久精品国产国产毛片| 九九爱精品视频在线观看| 91精品三级在线观看| 视频在线观看一区二区三区| 丰满迷人的少妇在线观看| 国产精品一区www在线观看| 亚洲av成人精品一区久久| 波野结衣二区三区在线| 女人精品久久久久毛片| 国产亚洲精品久久久com| 亚洲精品一区蜜桃| 全区人妻精品视频| 欧美精品国产亚洲| 99热全是精品| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 在线观看美女被高潮喷水网站| 亚洲欧美日韩另类电影网站| 国产av一区二区精品久久| 大话2 男鬼变身卡| 丰满乱子伦码专区| 国产精品国产三级国产专区5o| 日韩欧美一区视频在线观看| 汤姆久久久久久久影院中文字幕| 欧美日韩视频高清一区二区三区二| 少妇 在线观看| 日韩电影二区| 永久免费av网站大全| 亚洲精品色激情综合| 国语对白做爰xxxⅹ性视频网站| 亚洲美女搞黄在线观看| 免费高清在线观看日韩| 爱豆传媒免费全集在线观看| 亚洲国产精品成人久久小说| 久久久欧美国产精品| 美女大奶头黄色视频| 一级片'在线观看视频| 2021少妇久久久久久久久久久| 极品少妇高潮喷水抽搐| 成年美女黄网站色视频大全免费 | 国产精品一区www在线观看| 中文字幕久久专区| 亚洲一区二区三区欧美精品| 午夜免费鲁丝| av在线app专区| 老司机影院毛片| 精品国产一区二区久久| 国产在视频线精品| 精品视频人人做人人爽| 久久久午夜欧美精品| av一本久久久久| 久久久久国产精品人妻一区二区| 国产毛片在线视频| 少妇被粗大猛烈的视频| 黄色毛片三级朝国网站| 日本91视频免费播放| 全区人妻精品视频| 激情五月婷婷亚洲| 国产一级毛片在线| 波野结衣二区三区在线| 2018国产大陆天天弄谢| 日韩一区二区三区影片| 国产伦理片在线播放av一区| 五月伊人婷婷丁香| 性色avwww在线观看| 午夜福利视频在线观看免费| 伊人亚洲综合成人网| 一级毛片电影观看| 久久久精品区二区三区| 久久精品人人爽人人爽视色| 大码成人一级视频| 亚洲精品成人av观看孕妇| 少妇的逼水好多| 亚洲精品久久久久久婷婷小说| 免费大片18禁| 欧美成人午夜免费资源| 亚洲国产精品一区二区三区在线| 亚洲精品国产色婷婷电影| 一级黄片播放器| 伦精品一区二区三区| av免费观看日本| 免费观看无遮挡的男女| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| 国产视频首页在线观看| 嘟嘟电影网在线观看| 国产欧美亚洲国产| 各种免费的搞黄视频| 日韩亚洲欧美综合| 91精品三级在线观看| 日韩精品有码人妻一区| 一区二区三区乱码不卡18| 免费少妇av软件| 天天操日日干夜夜撸| 男女边吃奶边做爰视频| 晚上一个人看的免费电影| 亚洲人成77777在线视频| av免费在线看不卡| 韩国av在线不卡| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 中文字幕精品免费在线观看视频 | 日韩欧美一区视频在线观看| 一区二区三区精品91| 国产成人av激情在线播放 | 日产精品乱码卡一卡2卡三| 视频中文字幕在线观看| 自线自在国产av| 精品国产露脸久久av麻豆| 大码成人一级视频| 91精品国产九色| 人人妻人人澡人人爽人人夜夜| 99久久中文字幕三级久久日本| 国产老妇伦熟女老妇高清| 午夜免费观看性视频| 午夜日本视频在线| a级毛色黄片| 十八禁网站网址无遮挡| 亚洲一级一片aⅴ在线观看| 日韩亚洲欧美综合| 人妻制服诱惑在线中文字幕| 久久久国产一区二区| 青春草国产在线视频| 亚洲精品美女久久av网站| 成人国产麻豆网| 亚洲欧美清纯卡通| av专区在线播放| 蜜臀久久99精品久久宅男| 最近中文字幕高清免费大全6| 老熟女久久久| 天天躁夜夜躁狠狠久久av| 色视频在线一区二区三区| 免费黄网站久久成人精品| 国产av码专区亚洲av| 高清av免费在线| 日韩一区二区三区影片| 一级毛片黄色毛片免费观看视频| 国产视频首页在线观看| 老女人水多毛片| 久久久久久久久久久久大奶| 一二三四中文在线观看免费高清| 久久这里有精品视频免费| 国产精品无大码| 欧美日韩国产mv在线观看视频| 国产一区有黄有色的免费视频| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 久久久久久久亚洲中文字幕| 亚洲精品亚洲一区二区| 黄色毛片三级朝国网站| 男人添女人高潮全过程视频| 亚洲精品日韩av片在线观看| a 毛片基地| 久久久久久久久久成人| 国产精品人妻久久久影院| 国产乱来视频区| 五月开心婷婷网| 国产一区亚洲一区在线观看| 纯流量卡能插随身wifi吗| 精品熟女少妇av免费看| 黄片播放在线免费| 国产成人精品一,二区| 亚洲精品日韩在线中文字幕| 亚洲国产最新在线播放| 亚洲无线观看免费| 久久人人爽人人爽人人片va| 三级国产精品片| videosex国产| 中文字幕最新亚洲高清| 91成人精品电影| 精品人妻偷拍中文字幕| 18禁在线无遮挡免费观看视频| 男人添女人高潮全过程视频| 亚洲av成人精品一二三区| 丝袜脚勾引网站| 亚洲在久久综合| 桃花免费在线播放| 国产精品99久久99久久久不卡 | 夜夜骑夜夜射夜夜干| 少妇被粗大的猛进出69影院 | 在线观看三级黄色| 欧美成人午夜免费资源| 亚洲精品,欧美精品| 亚洲高清免费不卡视频| 午夜免费观看性视频| 精品国产乱码久久久久久小说| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 亚洲综合色网址| 国产无遮挡羞羞视频在线观看| 成人亚洲欧美一区二区av| 精品久久久久久久久av| 十八禁高潮呻吟视频| 亚洲精品亚洲一区二区| 亚洲四区av| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 婷婷成人精品国产| 母亲3免费完整高清在线观看 | 日韩制服骚丝袜av| 婷婷色av中文字幕| 欧美精品一区二区大全| 一本久久精品| 亚洲国产日韩一区二区| 欧美另类一区| 国产高清三级在线| 又大又黄又爽视频免费| 麻豆精品久久久久久蜜桃| 免费高清在线观看日韩| 啦啦啦中文免费视频观看日本| 亚洲欧美色中文字幕在线| 欧美日韩视频精品一区| 水蜜桃什么品种好| 高清午夜精品一区二区三区| 卡戴珊不雅视频在线播放| 男女无遮挡免费网站观看| 免费观看a级毛片全部| 日本欧美国产在线视频| 大陆偷拍与自拍| 成人影院久久| 精品一区二区三卡| 极品人妻少妇av视频| 日韩一区二区三区影片| 亚洲av欧美aⅴ国产| 久久国产亚洲av麻豆专区| 亚洲丝袜综合中文字幕| 亚州av有码| 18禁动态无遮挡网站| 日韩人妻高清精品专区| 高清欧美精品videossex| 日本黄色日本黄色录像| 成人漫画全彩无遮挡| 丝袜喷水一区| 国产成人91sexporn| 性色av一级| 大话2 男鬼变身卡| 日韩免费高清中文字幕av| 高清视频免费观看一区二区| 亚洲欧洲精品一区二区精品久久久 | 午夜福利,免费看| 精品午夜福利在线看| 特大巨黑吊av在线直播| 亚洲一区二区三区欧美精品| 简卡轻食公司| 亚洲成人av在线免费| 欧美性感艳星| 热99久久久久精品小说推荐| 午夜久久久在线观看| 欧美成人午夜免费资源| 亚洲熟女精品中文字幕| 免费黄色在线免费观看| 老司机影院毛片| 国产成人av激情在线播放 | 爱豆传媒免费全集在线观看| 精品人妻一区二区三区麻豆| 草草在线视频免费看| 一边摸一边做爽爽视频免费| 边亲边吃奶的免费视频| 老司机亚洲免费影院| 国产成人精品福利久久| a级毛片免费高清观看在线播放| 免费黄频网站在线观看国产| 美女国产视频在线观看| 蜜桃国产av成人99| 免费av不卡在线播放| 成人黄色视频免费在线看| 日本wwww免费看| 久久综合国产亚洲精品| 久久精品国产亚洲网站| 国产亚洲av片在线观看秒播厂| 天堂中文最新版在线下载| 国产男女内射视频| 99热6这里只有精品| 一个人看视频在线观看www免费| 亚洲av成人精品一二三区| 国内精品宾馆在线| 一级毛片aaaaaa免费看小| 久久久午夜欧美精品| 亚洲欧美日韩卡通动漫| 十八禁高潮呻吟视频| 亚洲av欧美aⅴ国产| 欧美另类一区| av线在线观看网站| 少妇丰满av| 婷婷成人精品国产| 久久性视频一级片| 激情视频va一区二区三区| av超薄肉色丝袜交足视频| 亚洲自偷自拍图片 自拍| 国产免费现黄频在线看| 青青草视频在线视频观看| 少妇被粗大的猛进出69影院| 一级毛片电影观看| 在线永久观看黄色视频| 夫妻午夜视频| 久久久精品免费免费高清| av又黄又爽大尺度在线免费看| 亚洲av成人不卡在线观看播放网| 国产欧美日韩一区二区精品| 国产精品亚洲一级av第二区| 中国美女看黄片| 色老头精品视频在线观看| 99精国产麻豆久久婷婷| 美国免费a级毛片| 美女高潮到喷水免费观看| 国产一区二区 视频在线| 考比视频在线观看| 亚洲avbb在线观看| 美女午夜性视频免费| aaaaa片日本免费| 午夜视频精品福利| 国产熟女午夜一区二区三区| 在线观看免费视频日本深夜| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美精品综合一区二区三区| 一级,二级,三级黄色视频| 久久精品亚洲熟妇少妇任你| 国产伦人伦偷精品视频| 国产高清国产精品国产三级| 老熟女久久久| 99久久精品国产亚洲精品| 午夜两性在线视频| 天堂中文最新版在线下载| 中文字幕高清在线视频| 国产免费现黄频在线看| 久久精品国产综合久久久| 五月开心婷婷网| 国精品久久久久久国模美| 99精品欧美一区二区三区四区| 在线播放国产精品三级| 宅男免费午夜| 天堂8中文在线网| 一进一出好大好爽视频| 国产一区二区 视频在线| 国产91精品成人一区二区三区 | 大香蕉久久网| 一级毛片女人18水好多| 美女高潮到喷水免费观看| 亚洲色图综合在线观看| 亚洲精品国产精品久久久不卡| 我的亚洲天堂| 国产三级黄色录像| 亚洲国产欧美网| 国产成人精品在线电影| 9191精品国产免费久久| 国产av精品麻豆| 欧美乱码精品一区二区三区| 亚洲性夜色夜夜综合| 午夜91福利影院| 久久国产精品大桥未久av| 色播在线永久视频| 黄色视频,在线免费观看| 亚洲成人手机| 欧美激情高清一区二区三区| 欧美人与性动交α欧美软件| 黄色视频,在线免费观看| 在线观看舔阴道视频| 国产av精品麻豆| 国产在视频线精品| 69av精品久久久久久 | 日本av免费视频播放| 一区二区三区激情视频| 欧美成狂野欧美在线观看| 国产一区二区三区视频了| 后天国语完整版免费观看| 黄色 视频免费看| 午夜精品久久久久久毛片777| 黄色 视频免费看| 最黄视频免费看| 国产高清视频在线播放一区| 老司机午夜十八禁免费视频| 国产在线观看jvid| 亚洲avbb在线观看| 欧美老熟妇乱子伦牲交| 怎么达到女性高潮| 91成年电影在线观看| 色综合婷婷激情| 国产成人精品在线电影| 高清欧美精品videossex| 亚洲七黄色美女视频| 色老头精品视频在线观看| 母亲3免费完整高清在线观看| 亚洲欧美精品综合一区二区三区| 国产精品偷伦视频观看了| 欧美日韩黄片免| 久久久国产成人免费| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 色在线成人网| 国产亚洲午夜精品一区二区久久| 亚洲成人国产一区在线观看| 亚洲 欧美一区二区三区| 日韩三级视频一区二区三区| 热re99久久国产66热| 午夜福利在线免费观看网站| 亚洲va日本ⅴa欧美va伊人久久| 99国产精品一区二区蜜桃av | 80岁老熟妇乱子伦牲交| 天堂8中文在线网| 亚洲精品国产精品久久久不卡| 丝袜美足系列| 亚洲七黄色美女视频| 国产精品99久久99久久久不卡| 18在线观看网站| 妹子高潮喷水视频| 女警被强在线播放| 国产亚洲精品第一综合不卡| 自拍欧美九色日韩亚洲蝌蚪91| 国产不卡一卡二| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美三级三区| 窝窝影院91人妻| 国产男女超爽视频在线观看| 亚洲自偷自拍图片 自拍| 精品福利永久在线观看| 久久久久久免费高清国产稀缺| 国产成人影院久久av| 999精品在线视频| av福利片在线| 少妇精品久久久久久久| 亚洲欧洲日产国产| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 性色av乱码一区二区三区2| 老熟妇乱子伦视频在线观看| 超色免费av| 欧美日韩亚洲高清精品| 18禁裸乳无遮挡动漫免费视频| 精品视频人人做人人爽| 天堂8中文在线网| 亚洲av日韩精品久久久久久密| 亚洲av成人一区二区三| 天天操日日干夜夜撸| 中文亚洲av片在线观看爽 | 国精品久久久久久国模美| 亚洲性夜色夜夜综合| 欧美激情久久久久久爽电影 | 久久精品国产a三级三级三级| 十分钟在线观看高清视频www| 日韩精品免费视频一区二区三区| 精品一区二区三区av网在线观看 | 色视频在线一区二区三区| 女性生殖器流出的白浆| 五月天丁香电影| 十八禁人妻一区二区| 亚洲国产成人一精品久久久| 亚洲成国产人片在线观看| 久久久精品免费免费高清| 黄色怎么调成土黄色| 一个人免费在线观看的高清视频| 欧美性长视频在线观看| 母亲3免费完整高清在线观看| 热re99久久精品国产66热6| 久久精品国产亚洲av高清一级| 在线观看66精品国产| 99香蕉大伊视频| 国产野战对白在线观看| 欧美中文综合在线视频| 日本黄色日本黄色录像| 老司机深夜福利视频在线观看| 久久亚洲真实| 久久久久久久大尺度免费视频| 日韩一卡2卡3卡4卡2021年| 免费少妇av软件| 成人18禁高潮啪啪吃奶动态图| 成人免费观看视频高清| 曰老女人黄片| 波多野结衣一区麻豆| 精品亚洲成a人片在线观看| √禁漫天堂资源中文www| 精品人妻熟女毛片av久久网站| 午夜日韩欧美国产| 一级片免费观看大全| 久久久久久免费高清国产稀缺| 极品人妻少妇av视频| 午夜福利乱码中文字幕| 高潮久久久久久久久久久不卡| 99国产综合亚洲精品| 9热在线视频观看99| 亚洲熟女精品中文字幕| 国产一区有黄有色的免费视频| 制服人妻中文乱码| videosex国产| av免费在线观看网站| 精品熟女少妇八av免费久了| 精品一区二区三卡| 免费在线观看日本一区| 亚洲七黄色美女视频| 久久久精品区二区三区| 十分钟在线观看高清视频www| 国产精品国产av在线观看| 91精品国产国语对白视频| 亚洲精品中文字幕一二三四区 | 女人精品久久久久毛片| 亚洲精品一二三| 大片免费播放器 马上看| 母亲3免费完整高清在线观看| 女人被躁到高潮嗷嗷叫费观| av网站在线播放免费| 亚洲精品美女久久久久99蜜臀| 中文字幕精品免费在线观看视频| 91国产中文字幕| 精品亚洲成国产av| 伦理电影免费视频| 国产精品自产拍在线观看55亚洲 | 久久中文看片网| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 99国产极品粉嫩在线观看| 黄色丝袜av网址大全| 午夜福利免费观看在线| 女性被躁到高潮视频| 久久久国产成人免费| 日韩 欧美 亚洲 中文字幕| 国产精品欧美亚洲77777| 一进一出好大好爽视频| 看免费av毛片| 男女免费视频国产| 女同久久另类99精品国产91| 亚洲色图 男人天堂 中文字幕| 国产成人欧美| 91麻豆av在线| 国产午夜精品久久久久久| 在线十欧美十亚洲十日本专区| 两性夫妻黄色片| 在线观看人妻少妇| 国产成人av激情在线播放| 精品免费久久久久久久清纯 | 国产精品九九99| 黄色视频在线播放观看不卡| 午夜老司机福利片| 久久久精品免费免费高清| 丝袜美足系列| 汤姆久久久久久久影院中文字幕| 在线观看人妻少妇| 亚洲欧美日韩高清在线视频 | 国产精品久久久久久精品古装| 大码成人一级视频| 国产欧美亚洲国产| 国产精品久久久久久精品电影小说| 99热国产这里只有精品6| 国产淫语在线视频| 亚洲中文日韩欧美视频| 日韩大码丰满熟妇| 少妇精品久久久久久久| 国产99久久九九免费精品| 亚洲黑人精品在线| 99久久人妻综合| 国产高清videossex| 成人亚洲精品一区在线观看| 国产精品亚洲一级av第二区| 一本色道久久久久久精品综合| 一二三四社区在线视频社区8| 国产老妇伦熟女老妇高清| 女警被强在线播放| 亚洲视频免费观看视频| 91九色精品人成在线观看| a在线观看视频网站| 少妇精品久久久久久久| 国产av国产精品国产| 国产精品久久久久久人妻精品电影 | 欧美精品人与动牲交sv欧美| 国产深夜福利视频在线观看| 最新美女视频免费是黄的| 男女之事视频高清在线观看| 美女视频免费永久观看网站| 国产精品久久电影中文字幕 | 久久亚洲真实| 97在线人人人人妻| 黄色怎么调成土黄色| 母亲3免费完整高清在线观看| 午夜免费成人在线视频| 午夜免费鲁丝| 一级片免费观看大全| 男女午夜视频在线观看| 久久久久久久久免费视频了| 男女下面插进去视频免费观看| 999久久久国产精品视频| 国产人伦9x9x在线观看| 精品免费久久久久久久清纯 | 老司机影院毛片| 精品一区二区三区视频在线观看免费 | 国产91精品成人一区二区三区 | 香蕉国产在线看| 国产主播在线观看一区二区| 欧美日韩国产mv在线观看视频| 国产av又大| 80岁老熟妇乱子伦牲交| 久久精品亚洲精品国产色婷小说| 别揉我奶头~嗯~啊~动态视频| 久久久精品免费免费高清| 亚洲专区字幕在线| 日韩有码中文字幕| 黄片小视频在线播放| 激情在线观看视频在线高清 | 中文字幕精品免费在线观看视频| 欧美日韩福利视频一区二区| 成年动漫av网址| 一本综合久久免费| 国产欧美日韩一区二区精品|