• <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
    美麗的黃石公園
    天堂中文最新版在线下载| 亚洲精品日本国产第一区| 国产在视频线精品| 欧美亚洲日本最大视频资源| 精品一品国产午夜福利视频| 伦理电影大哥的女人| 香蕉丝袜av| 亚洲美女黄色视频免费看| 七月丁香在线播放| 校园人妻丝袜中文字幕| 大片免费播放器 马上看| 久久精品国产综合久久久| 欧美另类一区| 男人爽女人下面视频在线观看| 天天躁夜夜躁狠狠躁躁| 人妻一区二区av| 国产女主播在线喷水免费视频网站| 日韩中字成人| av国产精品久久久久影院| 国产成人aa在线观看| 97在线人人人人妻| 国产成人精品福利久久| 国产 一区精品| 香蕉国产在线看| 国产成人免费观看mmmm| 91精品伊人久久大香线蕉| 麻豆乱淫一区二区| 国产又色又爽无遮挡免| 男的添女的下面高潮视频| 亚洲伊人久久精品综合| 亚洲国产最新在线播放| 新久久久久国产一级毛片| 十八禁网站网址无遮挡| 丝袜人妻中文字幕| xxx大片免费视频| 亚洲成人手机| 日韩制服丝袜自拍偷拍| 看十八女毛片水多多多| 视频区图区小说| 成人18禁高潮啪啪吃奶动态图| 日本欧美国产在线视频| 亚洲第一av免费看| 最近最新中文字幕大全免费视频 | 老汉色∧v一级毛片| 天堂俺去俺来也www色官网| 国产又爽黄色视频| 国产熟女欧美一区二区| 成人国产av品久久久| 多毛熟女@视频| 咕卡用的链子| 在线观看人妻少妇| 国产av码专区亚洲av| 欧美日韩一级在线毛片| 国产1区2区3区精品| 亚洲欧美日韩另类电影网站| 男女无遮挡免费网站观看| 黑丝袜美女国产一区| 欧美日韩视频精品一区| 在线观看美女被高潮喷水网站| 狂野欧美激情性bbbbbb| 精品人妻在线不人妻| 黄片小视频在线播放| 免费高清在线观看日韩| 少妇的丰满在线观看| 久久精品国产综合久久久| 狠狠精品人妻久久久久久综合| 在线观看一区二区三区激情| kizo精华| 国产精品一国产av| 高清不卡的av网站| 亚洲美女搞黄在线观看| 免费大片黄手机在线观看| 国产黄色免费在线视频| 1024视频免费在线观看| 日本av免费视频播放| 久久久久精品久久久久真实原创| 校园人妻丝袜中文字幕| 乱人伦中国视频| 美女高潮到喷水免费观看| 男女啪啪激烈高潮av片| 久久精品国产自在天天线| 国产精品女同一区二区软件| 亚洲精品中文字幕在线视频| 久久精品人人爽人人爽视色| 亚洲av男天堂| 国产老妇伦熟女老妇高清| 国产无遮挡羞羞视频在线观看| 丝袜喷水一区| 久久99一区二区三区| 日韩三级伦理在线观看| 人妻少妇偷人精品九色| 午夜福利一区二区在线看| 叶爱在线成人免费视频播放| 国产男人的电影天堂91| 最近的中文字幕免费完整| 午夜福利乱码中文字幕| 欧美精品一区二区免费开放| 各种免费的搞黄视频| 午夜免费鲁丝| 亚洲国产欧美日韩在线播放| 老司机影院毛片| 国产麻豆69| av卡一久久| 国产爽快片一区二区三区| 99re6热这里在线精品视频| 久久精品亚洲av国产电影网| 黄片播放在线免费| 18禁观看日本| 日韩熟女老妇一区二区性免费视频| 最新中文字幕久久久久| 欧美人与性动交α欧美软件| 亚洲第一av免费看| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 建设人人有责人人尽责人人享有的| 黄片小视频在线播放| 久久精品亚洲av国产电影网| 国产成人免费无遮挡视频| 18禁动态无遮挡网站| 1024香蕉在线观看| 色94色欧美一区二区| 极品少妇高潮喷水抽搐| 丝袜人妻中文字幕| 日本免费在线观看一区| 在线精品无人区一区二区三| 亚洲男人天堂网一区| 免费高清在线观看日韩| 久久久久久久久久人人人人人人| 狠狠婷婷综合久久久久久88av| 最新的欧美精品一区二区| 日韩精品有码人妻一区| 久久久国产一区二区| 国产日韩一区二区三区精品不卡| 国产1区2区3区精品| 九色亚洲精品在线播放| 香蕉精品网在线| 国产日韩欧美亚洲二区| 久久毛片免费看一区二区三区| 久久精品国产亚洲av涩爱| 香蕉丝袜av| 69精品国产乱码久久久| 各种免费的搞黄视频| 女性被躁到高潮视频| 久久精品熟女亚洲av麻豆精品| 最近最新中文字幕大全免费视频 | 大香蕉久久网| 一级片'在线观看视频| 亚洲欧美中文字幕日韩二区| 国产欧美日韩一区二区三区在线| 亚洲精品第二区| 又大又黄又爽视频免费| 国产福利在线免费观看视频| 青草久久国产| 一本—道久久a久久精品蜜桃钙片| 中国三级夫妇交换| 国产一区二区 视频在线| 五月天丁香电影| 国产精品久久久久久精品古装| 18禁观看日本| 久久久久久久亚洲中文字幕| 久热这里只有精品99| 捣出白浆h1v1| 久久99一区二区三区| 99香蕉大伊视频| 一边亲一边摸免费视频| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 国产成人欧美| 精品视频人人做人人爽| 中文字幕人妻丝袜一区二区 | 欧美老熟妇乱子伦牲交| 制服丝袜香蕉在线| 日韩不卡一区二区三区视频在线| 精品午夜福利在线看| 欧美日韩亚洲国产一区二区在线观看 | 高清av免费在线| 国产黄色视频一区二区在线观看| 岛国毛片在线播放| 日日撸夜夜添| 日韩中文字幕欧美一区二区 | 青春草视频在线免费观看| 1024视频免费在线观看| 99九九在线精品视频| 国产日韩欧美在线精品| 91精品三级在线观看| 亚洲精品视频女| 亚洲av福利一区| 极品少妇高潮喷水抽搐| 亚洲欧美精品自产自拍| 欧美另类一区| 国产日韩欧美亚洲二区| 精品人妻偷拍中文字幕| 国产精品久久久久久av不卡| 亚洲视频免费观看视频| 亚洲伊人久久精品综合| 少妇 在线观看| 十八禁网站网址无遮挡| 国产男女内射视频| 欧美亚洲日本最大视频资源| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 亚洲五月色婷婷综合| 不卡视频在线观看欧美| www.精华液| 国产精品久久久久久久久免| 中文天堂在线官网| 在线观看美女被高潮喷水网站| 伊人久久国产一区二区| 久久精品亚洲av国产电影网| 免费日韩欧美在线观看| 免费大片黄手机在线观看| 一本色道久久久久久精品综合| 亚洲欧美精品综合一区二区三区 | 色网站视频免费| 菩萨蛮人人尽说江南好唐韦庄| av卡一久久| 丰满饥渴人妻一区二区三| 久久婷婷青草| 午夜福利在线观看免费完整高清在| 成人国产麻豆网| 国产成人a∨麻豆精品| 国产欧美日韩一区二区三区在线| 999久久久国产精品视频| 国产午夜精品一二区理论片| 欧美日韩一区二区视频在线观看视频在线| 中文欧美无线码| 大片电影免费在线观看免费| 9191精品国产免费久久| 国产爽快片一区二区三区| 亚洲男人天堂网一区| 交换朋友夫妻互换小说| 大码成人一级视频| 久久精品国产综合久久久| 建设人人有责人人尽责人人享有的| 久久精品熟女亚洲av麻豆精品| 狠狠婷婷综合久久久久久88av| 精品久久久久久电影网| 性色avwww在线观看| av在线app专区| 国产女主播在线喷水免费视频网站| av网站在线播放免费| 老司机亚洲免费影院| 在线精品无人区一区二区三| www.精华液| 成年av动漫网址| 永久免费av网站大全| 久久青草综合色| 欧美精品高潮呻吟av久久| 最新的欧美精品一区二区| 一区二区三区四区激情视频| 超碰97精品在线观看| av视频免费观看在线观看| 美女主播在线视频| 亚洲精品第二区| 精品国产一区二区三区久久久樱花| 中文字幕最新亚洲高清| videos熟女内射| 欧美亚洲 丝袜 人妻 在线| 午夜日韩欧美国产| 宅男免费午夜| av网站在线播放免费| 亚洲一区中文字幕在线| kizo精华| 熟女少妇亚洲综合色aaa.| 免费黄网站久久成人精品| 日韩av不卡免费在线播放| av卡一久久| 尾随美女入室| 黑人猛操日本美女一级片| av线在线观看网站| 69精品国产乱码久久久| 国产成人aa在线观看| 亚洲av电影在线进入| 久久影院123| 亚洲第一av免费看| 日日撸夜夜添| 日韩制服丝袜自拍偷拍| 国产女主播在线喷水免费视频网站| 久久久国产一区二区| 精品少妇一区二区三区视频日本电影 | 久久久久久久精品精品| 午夜精品国产一区二区电影| 久久国产精品大桥未久av| 极品人妻少妇av视频| 精品卡一卡二卡四卡免费| 亚洲少妇的诱惑av| 18在线观看网站| 婷婷色综合www| 国产伦理片在线播放av一区| 免费观看a级毛片全部| 欧美中文综合在线视频| 午夜久久久在线观看| 亚洲国产最新在线播放| 黄色 视频免费看| 亚洲国产日韩一区二区| 亚洲一区二区三区欧美精品| 中国三级夫妇交换| 久久精品久久精品一区二区三区| av女优亚洲男人天堂| 80岁老熟妇乱子伦牲交| 丝袜人妻中文字幕| 日本av手机在线免费观看| 久久久久久久国产电影| 叶爱在线成人免费视频播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲一区中文字幕在线| 18+在线观看网站| 2022亚洲国产成人精品| 天堂8中文在线网| www.自偷自拍.com| 最近中文字幕2019免费版| 精品久久蜜臀av无| 丝袜脚勾引网站| 久久狼人影院| 久久99一区二区三区| 亚洲国产欧美在线一区| 亚洲三区欧美一区| 蜜桃在线观看..| 国产精品成人在线| 亚洲,一卡二卡三卡| 亚洲精品一区蜜桃| 纵有疾风起免费观看全集完整版| 9色porny在线观看| 亚洲第一av免费看| 久久精品亚洲av国产电影网| 日本欧美国产在线视频| 久久人妻熟女aⅴ| 春色校园在线视频观看| 亚洲成av片中文字幕在线观看 | 亚洲欧美清纯卡通| 黄色一级大片看看| 久久午夜综合久久蜜桃| 国产深夜福利视频在线观看| 亚洲精品久久久久久婷婷小说| 欧美人与性动交α欧美软件| 精品国产露脸久久av麻豆| 久久青草综合色| 亚洲精品一区蜜桃| 一级,二级,三级黄色视频| 免费观看性生交大片5| 日日啪夜夜爽| 久久精品亚洲av国产电影网| a 毛片基地| 午夜福利在线观看免费完整高清在| 大片免费播放器 马上看| 欧美97在线视频| 色婷婷av一区二区三区视频| 高清在线视频一区二区三区| 一区二区av电影网| 国产爽快片一区二区三区| 国产 一区精品| 麻豆精品久久久久久蜜桃| 老鸭窝网址在线观看| 免费日韩欧美在线观看| 久久精品国产亚洲av高清一级| 丝袜美足系列| 夫妻性生交免费视频一级片| 精品一区在线观看国产| 国产淫语在线视频| 国产男女内射视频| 人妻 亚洲 视频| 国产一区亚洲一区在线观看| 午夜影院在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲av片在线观看秒播厂| 国产国语露脸激情在线看| 如何舔出高潮| 国产欧美日韩综合在线一区二区| 中文字幕亚洲精品专区| 美女主播在线视频| 999久久久国产精品视频| 国产熟女午夜一区二区三区| 在线亚洲精品国产二区图片欧美| 久久久久精品久久久久真实原创| 伊人久久大香线蕉亚洲五| 国产熟女午夜一区二区三区| 国产欧美日韩一区二区三区在线| 久久久久精品久久久久真实原创| 伊人久久大香线蕉亚洲五| 在线观看免费日韩欧美大片| 好男人视频免费观看在线| 国产不卡av网站在线观看| 多毛熟女@视频| 亚洲三级黄色毛片| 2018国产大陆天天弄谢| 两个人免费观看高清视频| 国产男人的电影天堂91| 亚洲精品久久成人aⅴ小说| 亚洲综合色网址| 欧美精品高潮呻吟av久久| 亚洲成国产人片在线观看| 亚洲国产精品国产精品| 亚洲国产av影院在线观看| 国产深夜福利视频在线观看| 国产成人免费观看mmmm| 欧美精品高潮呻吟av久久| 爱豆传媒免费全集在线观看| 亚洲国产看品久久| 国产精品免费大片| 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 国产成人午夜福利电影在线观看| 久久人人爽av亚洲精品天堂| 在线天堂中文资源库| 精品少妇久久久久久888优播| 亚洲男人天堂网一区| 国产一区二区 视频在线| 尾随美女入室| av国产精品久久久久影院| 看免费成人av毛片| 王馨瑶露胸无遮挡在线观看| 在线天堂最新版资源| 欧美人与性动交α欧美软件| 边亲边吃奶的免费视频| 久久精品夜色国产| 久久99蜜桃精品久久| 久久久精品区二区三区| 波多野结衣一区麻豆| 欧美亚洲日本最大视频资源| 99热国产这里只有精品6| 999久久久国产精品视频| 精品国产乱码久久久久久小说| 中文字幕制服av| 只有这里有精品99| 我要看黄色一级片免费的| 国产一区二区激情短视频 | 久久精品熟女亚洲av麻豆精品| 亚洲内射少妇av| 80岁老熟妇乱子伦牲交| 亚洲国产精品一区三区| 久久精品国产亚洲av涩爱| 久久久久久久精品精品| 亚洲中文av在线| 99久久综合免费| 波野结衣二区三区在线| 99久久精品国产国产毛片| 国产成人精品久久二区二区91 | 啦啦啦在线免费观看视频4| 国产一区二区激情短视频 | 国产成人精品在线电影| 一本久久精品| 国产成人精品久久久久久| 老司机亚洲免费影院| 电影成人av| 免费黄色在线免费观看| 日本色播在线视频| 十八禁高潮呻吟视频| 女人精品久久久久毛片| 日韩精品免费视频一区二区三区| 国产高清国产精品国产三级| 久久鲁丝午夜福利片| 色播在线永久视频| 人体艺术视频欧美日本| 丝袜脚勾引网站| 日产精品乱码卡一卡2卡三| 久久久久久免费高清国产稀缺| 欧美日韩一级在线毛片| 免费黄网站久久成人精品| 亚洲国产av影院在线观看| 国产日韩欧美亚洲二区| 有码 亚洲区| 亚洲第一青青草原| 国产亚洲午夜精品一区二区久久| 97人妻天天添夜夜摸| 国产精品久久久久久精品电影小说| tube8黄色片| 日韩熟女老妇一区二区性免费视频| 视频区图区小说| 国产精品女同一区二区软件| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| 久久久久久久久久久久大奶| 国产精品 欧美亚洲| 日韩中文字幕欧美一区二区 | 成人手机av| 亚洲婷婷狠狠爱综合网| 亚洲第一av免费看| 亚洲欧洲国产日韩| 香蕉精品网在线| 只有这里有精品99| 一级黄片播放器| 激情视频va一区二区三区| 电影成人av| 在现免费观看毛片| 黄色配什么色好看| 18禁裸乳无遮挡动漫免费视频| tube8黄色片| 欧美精品国产亚洲| 精品人妻在线不人妻| 成人影院久久| 99国产精品免费福利视频| 啦啦啦中文免费视频观看日本| 女人久久www免费人成看片| 亚洲色图综合在线观看| 国产欧美亚洲国产| 男男h啪啪无遮挡| 熟女电影av网| 在线 av 中文字幕| 国产精品av久久久久免费| 高清欧美精品videossex| 男女下面插进去视频免费观看| 狂野欧美激情性bbbbbb| 黄片无遮挡物在线观看| 免费看av在线观看网站| 精品国产乱码久久久久久小说| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 久久精品熟女亚洲av麻豆精品| 国产亚洲最大av| 宅男免费午夜| 高清不卡的av网站| 人人妻人人澡人人爽人人夜夜| 精品久久久精品久久久| 精品福利永久在线观看| 曰老女人黄片| 国产黄频视频在线观看| 天天躁夜夜躁狠狠躁躁| 免费不卡的大黄色大毛片视频在线观看| 少妇精品久久久久久久| √禁漫天堂资源中文www| 亚洲精华国产精华液的使用体验| 久久久久久久精品精品| 夫妻性生交免费视频一级片| 国产在视频线精品| 电影成人av| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| av在线老鸭窝| av网站在线播放免费| 深夜精品福利| 好男人视频免费观看在线| 国产免费福利视频在线观看| 久久久精品区二区三区| 国产精品秋霞免费鲁丝片| 国产日韩欧美在线精品| 日本猛色少妇xxxxx猛交久久| 国产av一区二区精品久久| 亚洲,一卡二卡三卡| 美女脱内裤让男人舔精品视频| 国产片特级美女逼逼视频| 交换朋友夫妻互换小说| 精品国产超薄肉色丝袜足j| 又大又黄又爽视频免费| 免费少妇av软件| 亚洲国产精品一区三区| 成年人免费黄色播放视频| 精品第一国产精品| 90打野战视频偷拍视频| 在线 av 中文字幕| 成人手机av| 999久久久国产精品视频| 日韩大片免费观看网站| 午夜免费男女啪啪视频观看| 久久国产亚洲av麻豆专区| 日本-黄色视频高清免费观看| 免费在线观看视频国产中文字幕亚洲 | 国产一区二区三区综合在线观看| 性色avwww在线观看| 男女免费视频国产| 国产精品一二三区在线看| 色婷婷av一区二区三区视频| 国产精品不卡视频一区二区| 日韩免费高清中文字幕av| 亚洲美女黄色视频免费看| 亚洲欧美色中文字幕在线| 黄片播放在线免费| 免费观看在线日韩| 啦啦啦啦在线视频资源| av又黄又爽大尺度在线免费看| 不卡视频在线观看欧美| 亚洲国产精品成人久久小说| 欧美bdsm另类| 日韩一区二区视频免费看| 91精品伊人久久大香线蕉| 亚洲av中文av极速乱| 亚洲精品一区蜜桃| 黑人欧美特级aaaaaa片| 大话2 男鬼变身卡| 狠狠精品人妻久久久久久综合| 国产男女内射视频| 国产精品一区二区在线观看99| 亚洲欧美日韩另类电影网站| 日本av免费视频播放| 久久精品国产亚洲av天美| 色吧在线观看| 波多野结衣av一区二区av| 最近最新中文字幕免费大全7| 99久久人妻综合| 国产国语露脸激情在线看| 91成人精品电影| 国产成人精品婷婷| 考比视频在线观看| 在线观看美女被高潮喷水网站| 久久国产精品大桥未久av| 成年人午夜在线观看视频| 人成视频在线观看免费观看| 一级a爱视频在线免费观看| 在线观看免费日韩欧美大片| 亚洲精品中文字幕在线视频| 日韩伦理黄色片| av在线app专区| 亚洲内射少妇av| 国产毛片在线视频| 国产精品久久久久久av不卡| 亚洲欧美精品自产自拍| 国产精品99久久99久久久不卡 | av在线观看视频网站免费| 纯流量卡能插随身wifi吗| 最近最新中文字幕免费大全7| 精品一品国产午夜福利视频| 美女xxoo啪啪120秒动态图| 黄色配什么色好看| 超碰成人久久| 欧美亚洲 丝袜 人妻 在线| 免费黄网站久久成人精品|