• <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
    美麗的黃石公園
    国产在线一区二区三区精| 97在线视频观看| 91精品国产国语对白视频| 99热6这里只有精品| 日本欧美视频一区| 久久久久久久久久久免费av| 最新中文字幕久久久久| 啦啦啦中文免费视频观看日本| 一边摸一边做爽爽视频免费| 80岁老熟妇乱子伦牲交| 中文字幕人妻丝袜制服| 日韩中文字幕视频在线看片| 亚洲四区av| av在线app专区| 性高湖久久久久久久久免费观看| 草草在线视频免费看| 日韩成人av中文字幕在线观看| 午夜激情av网站| 国产成人精品一,二区| 亚洲av电影在线进入| 免费在线观看黄色视频的| 免费大片18禁| 久久精品久久久久久久性| 两个人免费观看高清视频| 纵有疾风起免费观看全集完整版| 国产国拍精品亚洲av在线观看| 中文字幕精品免费在线观看视频 | 大香蕉久久成人网| 视频中文字幕在线观看| 国产精品蜜桃在线观看| 黄色配什么色好看| 69精品国产乱码久久久| 黄色毛片三级朝国网站| 新久久久久国产一级毛片| 色94色欧美一区二区| 99久久中文字幕三级久久日本| 中文字幕av电影在线播放| 热re99久久精品国产66热6| 午夜免费男女啪啪视频观看| 欧美精品一区二区免费开放| 免费观看性生交大片5| 久久精品久久久久久久性| 免费在线观看完整版高清| 2022亚洲国产成人精品| 免费观看性生交大片5| 国产无遮挡羞羞视频在线观看| 18禁国产床啪视频网站| 亚洲av中文av极速乱| 九色成人免费人妻av| 午夜福利影视在线免费观看| 精品国产一区二区三区四区第35| 日韩一本色道免费dvd| 国产精品欧美亚洲77777| 美女国产高潮福利片在线看| 久久亚洲国产成人精品v| 日产精品乱码卡一卡2卡三| 亚洲精品国产色婷婷电影| 狂野欧美激情性bbbbbb| 国产成人一区二区在线| 日本wwww免费看| 美女福利国产在线| 高清毛片免费看| 久久精品人人爽人人爽视色| 日韩大片免费观看网站| 曰老女人黄片| 亚洲欧美成人精品一区二区| 观看av在线不卡| 久久99热这里只频精品6学生| 另类精品久久| 色94色欧美一区二区| 精品一区二区三区四区五区乱码 | 久久亚洲国产成人精品v| 欧美日韩成人在线一区二区| 免费观看无遮挡的男女| av在线老鸭窝| 国产成人精品久久久久久| 国产成人精品婷婷| 国产精品一二三区在线看| 亚洲欧美成人精品一区二区| 99热国产这里只有精品6| 国产免费视频播放在线视频| av在线app专区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧洲日产国产| 天堂俺去俺来也www色官网| 国产伦理片在线播放av一区| 制服诱惑二区| 久久午夜福利片| 亚洲少妇的诱惑av| 亚洲第一区二区三区不卡| 大码成人一级视频| 久久精品久久久久久久性| 十八禁网站网址无遮挡| 精品国产一区二区久久| av卡一久久| 涩涩av久久男人的天堂| 久久免费观看电影| 啦啦啦中文免费视频观看日本| 精品99又大又爽又粗少妇毛片| 波野结衣二区三区在线| 亚洲av福利一区| 欧美精品亚洲一区二区| 国产在线视频一区二区| 极品人妻少妇av视频| 午夜福利影视在线免费观看| 国产精品国产av在线观看| 只有这里有精品99| 大片免费播放器 马上看| 久久久a久久爽久久v久久| 日本黄大片高清| 亚洲内射少妇av| 少妇人妻精品综合一区二区| 大香蕉久久网| 国产亚洲欧美精品永久| 欧美亚洲日本最大视频资源| 你懂的网址亚洲精品在线观看| 国产男人的电影天堂91| 黄色 视频免费看| 国产白丝娇喘喷水9色精品| 桃花免费在线播放| 水蜜桃什么品种好| 国产一区二区激情短视频 | 午夜福利影视在线免费观看| 97超碰精品成人国产| 观看av在线不卡| 午夜91福利影院| 国产色爽女视频免费观看| 亚洲国产精品一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| 亚洲成人一二三区av| 国产白丝娇喘喷水9色精品| 99久久中文字幕三级久久日本| 美女国产高潮福利片在线看| 九草在线视频观看| 日韩三级伦理在线观看| a级毛色黄片| 亚洲精品美女久久av网站| 日日撸夜夜添| 国产69精品久久久久777片| 日韩熟女老妇一区二区性免费视频| 日日撸夜夜添| 黄网站色视频无遮挡免费观看| 日本黄色日本黄色录像| 国产精品国产三级国产av玫瑰| 亚洲国产色片| 国产在线视频一区二区| xxx大片免费视频| 51国产日韩欧美| 欧美激情国产日韩精品一区| 日韩三级伦理在线观看| 丝袜在线中文字幕| 精品熟女少妇av免费看| 校园人妻丝袜中文字幕| 国产xxxxx性猛交| 一区二区三区精品91| 久久国产亚洲av麻豆专区| 91成人精品电影| 亚洲精品视频女| 亚洲精品久久午夜乱码| 日日啪夜夜爽| 一级,二级,三级黄色视频| 看免费av毛片| 日韩电影二区| 最新中文字幕久久久久| 欧美xxxx性猛交bbbb| 国产一区二区三区综合在线观看 | 国产不卡av网站在线观看| 少妇被粗大猛烈的视频| a 毛片基地| 久久久久国产网址| 色哟哟·www| 欧美成人精品欧美一级黄| 九九在线视频观看精品| 综合色丁香网| 免费大片18禁| 日韩制服丝袜自拍偷拍| 视频在线观看一区二区三区| 9191精品国产免费久久| 久久这里有精品视频免费| 午夜av观看不卡| 一级毛片我不卡| 赤兔流量卡办理| 国产国语露脸激情在线看| 成年人午夜在线观看视频| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 国产爽快片一区二区三区| 日本黄大片高清| 免费大片18禁| 亚洲精品久久成人aⅴ小说| 日本色播在线视频| 国产精品99久久99久久久不卡 | 蜜桃在线观看..| 99热国产这里只有精品6| 国产片特级美女逼逼视频| 亚洲性久久影院| 日本vs欧美在线观看视频| 免费黄色在线免费观看| 99精国产麻豆久久婷婷| 亚洲av电影在线观看一区二区三区| 久久人人97超碰香蕉20202| 成人毛片60女人毛片免费| 午夜免费鲁丝| 午夜福利乱码中文字幕| 色婷婷av一区二区三区视频| 999精品在线视频| 免费观看av网站的网址| 美女福利国产在线| 国产又色又爽无遮挡免| 在线观看三级黄色| 日日啪夜夜爽| av福利片在线| 丰满饥渴人妻一区二区三| 国产av精品麻豆| 搡老乐熟女国产| 高清在线视频一区二区三区| 欧美日韩精品成人综合77777| av不卡在线播放| 亚洲,欧美精品.| 七月丁香在线播放| 国产无遮挡羞羞视频在线观看| 精品国产露脸久久av麻豆| 久久久久人妻精品一区果冻| 国产色婷婷99| 久久这里只有精品19| 侵犯人妻中文字幕一二三四区| 久久国产精品大桥未久av| 国产精品一区二区在线不卡| 国产精品久久久av美女十八| 波多野结衣一区麻豆| 成年人午夜在线观看视频| 精品酒店卫生间| 免费av中文字幕在线| 日韩制服丝袜自拍偷拍| 午夜日本视频在线| 久久久久网色| 亚洲成人av在线免费| 欧美精品av麻豆av| 天天影视国产精品| 久久精品久久久久久噜噜老黄| 亚洲欧美一区二区三区国产| 亚洲国产精品一区二区三区在线| 最后的刺客免费高清国语| 熟女av电影| 一本—道久久a久久精品蜜桃钙片| 免费看av在线观看网站| 婷婷色麻豆天堂久久| 婷婷色综合www| 久久久久精品人妻al黑| 国产精品一区二区在线不卡| 日韩欧美一区视频在线观看| 精品午夜福利在线看| 午夜精品国产一区二区电影| 国产精品无大码| 久久久久网色| 肉色欧美久久久久久久蜜桃| av国产久精品久网站免费入址| 中国美白少妇内射xxxbb| 一本—道久久a久久精品蜜桃钙片| 丝袜在线中文字幕| 91国产中文字幕| a级毛片黄视频| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| 99re6热这里在线精品视频| 午夜激情久久久久久久| 久久久a久久爽久久v久久| 中文字幕人妻熟女乱码| 999精品在线视频| 日韩一本色道免费dvd| 午夜福利影视在线免费观看| 老司机影院毛片| 人人妻人人添人人爽欧美一区卜| 欧美激情 高清一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久大尺度免费视频| 22中文网久久字幕| 免费日韩欧美在线观看| 国产有黄有色有爽视频| 亚洲婷婷狠狠爱综合网| 日韩成人伦理影院| 少妇的逼好多水| 亚洲少妇的诱惑av| 极品人妻少妇av视频| 欧美亚洲日本最大视频资源| 两性夫妻黄色片 | 欧美日韩一区二区视频在线观看视频在线| 精品少妇内射三级| 亚洲中文av在线| 久久综合国产亚洲精品| 亚洲精品一二三| 国产一区亚洲一区在线观看| 精品亚洲成国产av| 精品久久久久久电影网| 国产深夜福利视频在线观看| 成人国产麻豆网| 热re99久久精品国产66热6| 美女大奶头黄色视频| 午夜福利影视在线免费观看| 99国产综合亚洲精品| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 午夜免费鲁丝| 精品少妇黑人巨大在线播放| 国产日韩一区二区三区精品不卡| 国产精品国产av在线观看| 午夜激情久久久久久久| 国产色婷婷99| 久久人人爽人人爽人人片va| 成人综合一区亚洲| 色吧在线观看| 伊人久久国产一区二区| 少妇被粗大猛烈的视频| 午夜精品国产一区二区电影| 国产熟女午夜一区二区三区| 欧美 日韩 精品 国产| 国产精品国产三级专区第一集| 9热在线视频观看99| 精品国产露脸久久av麻豆| 91精品伊人久久大香线蕉| 精品人妻熟女毛片av久久网站| 国产精品一区二区在线不卡| 国产精品国产av在线观看| 九色成人免费人妻av| 国产成人精品一,二区| 亚洲av男天堂| 精品人妻一区二区三区麻豆| 久久久久久伊人网av| 九九爱精品视频在线观看| 亚洲天堂av无毛| 韩国高清视频一区二区三区| 街头女战士在线观看网站| 99热全是精品| 91国产中文字幕| 一边摸一边做爽爽视频免费| 一区二区三区精品91| 一级黄片播放器| 国产精品 国内视频| 亚洲精品日本国产第一区| 狠狠婷婷综合久久久久久88av| 精品卡一卡二卡四卡免费| 少妇高潮的动态图| 国产成人免费无遮挡视频| 亚洲精品美女久久av网站| 视频中文字幕在线观看| 久久鲁丝午夜福利片| 日本av手机在线免费观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 韩国av在线不卡| 成人毛片60女人毛片免费| av不卡在线播放| 日本免费在线观看一区| 国产精品熟女久久久久浪| 免费女性裸体啪啪无遮挡网站| 一本久久精品| 日本黄大片高清| 亚洲四区av| 天堂中文最新版在线下载| 又黄又粗又硬又大视频| 亚洲国产精品一区二区三区在线| 亚洲国产精品成人久久小说| 寂寞人妻少妇视频99o| 一区二区三区乱码不卡18| 亚洲av电影在线进入| 亚洲成人手机| 香蕉丝袜av| av电影中文网址| 全区人妻精品视频| 在线 av 中文字幕| 国产永久视频网站| 亚洲av免费高清在线观看| 精品熟女少妇av免费看| 在现免费观看毛片| 精品人妻偷拍中文字幕| 亚洲精品第二区| freevideosex欧美| 成人二区视频| 成人国产麻豆网| 韩国精品一区二区三区 | 亚洲精品久久成人aⅴ小说| 2018国产大陆天天弄谢| 亚洲人成网站在线观看播放| 丰满迷人的少妇在线观看| 国产不卡av网站在线观看| 国产高清国产精品国产三级| 久久精品国产亚洲av涩爱| 黄片播放在线免费| 欧美激情极品国产一区二区三区 | 亚洲欧美一区二区三区黑人 | 日韩视频在线欧美| 丰满乱子伦码专区| 伦精品一区二区三区| 人人妻人人添人人爽欧美一区卜| 伦理电影大哥的女人| 免费观看无遮挡的男女| 精品国产一区二区三区久久久樱花| 性色avwww在线观看| 日本免费在线观看一区| 丝袜脚勾引网站| 亚洲激情五月婷婷啪啪| 婷婷色av中文字幕| 亚洲av成人精品一二三区| 乱人伦中国视频| 一级毛片 在线播放| 在线观看一区二区三区激情| √禁漫天堂资源中文www| 久久毛片免费看一区二区三区| av天堂久久9| 欧美+日韩+精品| av在线播放精品| 色94色欧美一区二区| 成人无遮挡网站| 99精国产麻豆久久婷婷| 精品一区二区三区视频在线| 国产精品国产三级专区第一集| 在线观看www视频免费| 观看av在线不卡| 亚洲美女搞黄在线观看| 日韩三级伦理在线观看| 国产av一区二区精品久久| 成人18禁高潮啪啪吃奶动态图| 99re6热这里在线精品视频| 一区二区三区四区激情视频| 老司机影院毛片| 国产亚洲精品久久久com| 男的添女的下面高潮视频| 97精品久久久久久久久久精品| 大香蕉久久成人网| 熟妇人妻不卡中文字幕| 国产一区亚洲一区在线观看| 另类亚洲欧美激情| xxxhd国产人妻xxx| 人体艺术视频欧美日本| 午夜激情av网站| 亚洲精品,欧美精品| 黄色 视频免费看| 亚洲美女视频黄频| 午夜影院在线不卡| 搡老乐熟女国产| 日本黄色日本黄色录像| 精品久久久精品久久久| 国产永久视频网站| 国产一级毛片在线| 美女国产视频在线观看| 亚洲精品久久成人aⅴ小说| 亚洲av电影在线观看一区二区三区| 免费看光身美女| 国产免费一区二区三区四区乱码| 国产福利在线免费观看视频| 天堂俺去俺来也www色官网| 最黄视频免费看| 精品亚洲成a人片在线观看| 国产亚洲欧美精品永久| 欧美+日韩+精品| 精品少妇内射三级| 永久免费av网站大全| 三级国产精品片| 国产极品粉嫩免费观看在线| 国产又爽黄色视频| 9色porny在线观看| 久久国产精品大桥未久av| 中文乱码字字幕精品一区二区三区| 这个男人来自地球电影免费观看 | 蜜桃在线观看..| 久久久精品94久久精品| 国产在线视频一区二区| 一本久久精品| 亚洲av国产av综合av卡| 熟妇人妻不卡中文字幕| 一级片'在线观看视频| a级片在线免费高清观看视频| 高清不卡的av网站| 男女免费视频国产| 久久99热6这里只有精品| 人人妻人人澡人人看| 国产精品麻豆人妻色哟哟久久| 黑人猛操日本美女一级片| 亚洲精品美女久久久久99蜜臀 | 视频在线观看一区二区三区| 狠狠精品人妻久久久久久综合| 最近最新中文字幕免费大全7| 成人午夜精彩视频在线观看| 亚洲精品久久午夜乱码| 狠狠婷婷综合久久久久久88av| 少妇 在线观看| av福利片在线| 精品一区二区三卡| 久久久久久久精品精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久国产精品人妻一区二区| 日韩伦理黄色片| 亚洲第一av免费看| 人人澡人人妻人| 精品国产一区二区三区久久久樱花| 人妻少妇偷人精品九色| 午夜av观看不卡| 在线观看人妻少妇| 制服人妻中文乱码| 亚洲成人av在线免费| 九色亚洲精品在线播放| 久久精品久久久久久噜噜老黄| 日本欧美视频一区| 日韩成人av中文字幕在线观看| 欧美最新免费一区二区三区| 亚洲欧美一区二区三区国产| 日韩精品有码人妻一区| 成年动漫av网址| 日韩成人av中文字幕在线观看| 精品国产露脸久久av麻豆| 少妇人妻精品综合一区二区| 少妇人妻 视频| h视频一区二区三区| 啦啦啦中文免费视频观看日本| 免费观看在线日韩| 免费看av在线观看网站| 久久久久视频综合| 亚洲久久久国产精品| 中文字幕精品免费在线观看视频 | 97在线人人人人妻| 成人影院久久| 黑丝袜美女国产一区| 97精品久久久久久久久久精品| 视频区图区小说| 久久这里有精品视频免费| 高清不卡的av网站| 美女大奶头黄色视频| 久久久久久久国产电影| 色哟哟·www| 伦理电影大哥的女人| a 毛片基地| 日韩精品有码人妻一区| 午夜视频国产福利| 男男h啪啪无遮挡| 天天躁夜夜躁狠狠久久av| 内地一区二区视频在线| 激情五月婷婷亚洲| 亚洲经典国产精华液单| 亚洲五月色婷婷综合| 亚洲国产精品专区欧美| 蜜桃在线观看..| 性色avwww在线观看| 免费人妻精品一区二区三区视频| 日本91视频免费播放| 国产精品免费大片| 18+在线观看网站| 久久综合国产亚洲精品| 少妇高潮的动态图| 性色avwww在线观看| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 中文欧美无线码| 欧美 日韩 精品 国产| 久久久久久伊人网av| 男人添女人高潮全过程视频| 青春草亚洲视频在线观看| 久久精品国产亚洲av涩爱| 卡戴珊不雅视频在线播放| 少妇高潮的动态图| 一级毛片我不卡| 丝瓜视频免费看黄片| 最新中文字幕久久久久| 十八禁高潮呻吟视频| 午夜激情久久久久久久| 国产午夜精品一二区理论片| 最近中文字幕高清免费大全6| 国产成人精品在线电影| 国产精品熟女久久久久浪| 国产精品国产三级专区第一集| 国产欧美日韩一区二区三区在线| 日韩大片免费观看网站| 久久毛片免费看一区二区三区| 亚洲色图综合在线观看| 日本免费在线观看一区| 国产黄色视频一区二区在线观看| 大香蕉久久成人网| 熟妇人妻不卡中文字幕| 又粗又硬又长又爽又黄的视频| 免费看不卡的av| 免费观看无遮挡的男女| 午夜福利在线观看免费完整高清在| 日韩 亚洲 欧美在线| 亚洲熟女精品中文字幕| 国产深夜福利视频在线观看| 国产亚洲欧美精品永久| 国产毛片在线视频| 国精品久久久久久国模美| 永久网站在线| 国产69精品久久久久777片| 亚洲综合色惰| 国产午夜精品一二区理论片| 卡戴珊不雅视频在线播放| 亚洲精品aⅴ在线观看| 999精品在线视频| 黑人高潮一二区| 日韩成人伦理影院| av国产久精品久网站免费入址| a级毛片在线看网站| 女人被躁到高潮嗷嗷叫费观| 日韩三级伦理在线观看| 亚洲性久久影院| 亚洲国产最新在线播放| 1024视频免费在线观看| 熟妇人妻不卡中文字幕| 国产男人的电影天堂91| 亚洲久久久国产精品| 亚洲精品国产av成人精品| 免费高清在线观看视频在线观看| 国产在线一区二区三区精| 夫妻午夜视频| 久久精品人人爽人人爽视色| 曰老女人黄片| 如日韩欧美国产精品一区二区三区| 高清在线视频一区二区三区| www.熟女人妻精品国产 | 激情五月婷婷亚洲|