呼倫湖水面動(dòng)態(tài)變化遙感監(jiān)測及氣候因素驅(qū)動(dòng)分析

萬華偉,康　峻,高　帥,申文明*(.環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心,北京 0009；.中國科學(xué)院遙感與數(shù)字地球研究所,北京 000；3.中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院,北京 00049)

?

呼倫湖水面動(dòng)態(tài)變化遙感監(jiān)測及氣候因素驅(qū)動(dòng)分析

萬華偉1,康峻2,3,高帥2,申文明1*(1.環(huán)境保護(hù)部衛(wèi)星環(huán)境應(yīng)用中心,北京 100029；2.中國科學(xué)院遙感與數(shù)字地球研究所,北京 100101；3.中國科學(xué)院大學(xué)資源與環(huán)境學(xué)院,北京 100049)

摘要：利用長時(shí)間序列的MODIS數(shù)據(jù),采用水體指數(shù)動(dòng)態(tài)分析的方法,對(duì)2000～2013年呼倫湖的水體面積進(jìn)行了動(dòng)態(tài)變化分析,并結(jié)合區(qū)域氣候數(shù)據(jù)進(jìn)行了驅(qū)動(dòng)力分析.結(jié)果表明:2000～2012年期間,呼倫湖流域水體面積從2286k m2減少至1773km2,減少22.4%,主要減少部分分布在湖體東北部和南部;2013年水體面積尤其是南部水域水體面積略有所恢復(fù).呼倫湖主體湖區(qū)面積的變化與年均溫呈現(xiàn)不顯著的負(fù)相關(guān)關(guān)系.

關(guān)鍵詞：呼倫湖；水體面積；時(shí)間序列；氣候因素

* 責(zé)任作者, 正高級(jí)工程師, shenwenm@sepa.gov.cn

呼倫湖是我國東北部以保護(hù)草原生態(tài)、濕地系統(tǒng)和珍稀瀕危鳥類為主的最大的綜合性自然保護(hù)區(qū),在呼倫貝爾草原的生態(tài)保護(hù)和經(jīng)濟(jì)發(fā)展中發(fā)揮著不可替代的重要作用[1-2].受全球變化和人類活動(dòng)的共同影響,呼倫湖出現(xiàn)水量減少、面積萎縮、水位下降等問題,備受關(guān)注[3-6],因此,開展呼倫湖流域及湖區(qū)水體面積的動(dòng)態(tài)監(jiān)測并分析其影響因素具有重要意義.

遙感技術(shù)以及數(shù)據(jù)處理能力的迅速發(fā)展,使得大尺度全球地表水體面積變化過程的定量、動(dòng)態(tài)觀測成為可能,通過遙感手段來獲取水體面積及其動(dòng)態(tài)變化的遙感過程模型構(gòu)建亦取得了較大的進(jìn)步[7].為摸清呼倫湖近年來的水體面積變化狀況,本研究通過搜集呼倫湖流域中國境內(nèi)部分(包括新巴爾虎左旗、新巴爾虎右旗、滿洲里市)2000～2013年的MODIS數(shù)據(jù)、氣候數(shù)據(jù)和其他輔助數(shù)據(jù),對(duì)呼倫湖豐水期水體面積進(jìn)行動(dòng)態(tài)監(jiān)測,并對(duì)其驅(qū)動(dòng)力進(jìn)行了初步分析.

1　數(shù)據(jù)與方法

1.1MODIS數(shù)據(jù)和氣候數(shù)據(jù)搜集

研究選取的數(shù)據(jù)為MODIS 500m分辨率反射率8天合成產(chǎn)品(MOD09A1),在可見光和近紅外區(qū)域有7個(gè)波段,該產(chǎn)品采用了最大值合成法,選擇8日內(nèi)具有反射率最大值的像元,用以保證最佳的數(shù)據(jù)質(zhì)量[8].

如圖1所示,呼倫湖主要水體的范圍為117°00′10″～117°41′40″E, 48°30′40″～49°20′40″N,流域范圍為115°31′51″～120°41′43″E, 47°17′28″～50°43′42″N.研究選取數(shù)據(jù)的條帶號(hào)為H25V04,時(shí)間范圍為2000～2013年,為消除因湖面結(jié)冰反射率發(fā)生變化,每年選取第89～273d各24景影像,其時(shí)間范圍處于4月初和9月底,地面平均氣溫高于零度(圖2).對(duì)于MODIS數(shù)據(jù),首先使用ENVI插件MCTK(MODIS conversion toolkit)工具提取第4波段(綠波段)與第6波段(中紅外波段)反射率數(shù)據(jù),采用Delaunay三角網(wǎng)最鄰近采樣法將Sinusoidal投影轉(zhuǎn)換為地理坐標(biāo)投影,再將整個(gè)影像經(jīng)投影變換由地理坐標(biāo)轉(zhuǎn)換為UTM投影(分度帶為51N),最后根據(jù)經(jīng)緯度范圍,裁剪得到呼倫湖流域與呼倫湖主體湖區(qū)的綠波段、中紅外波段反射率影像.

圖1　呼倫湖地理位置示意Fig.1　Location of study area

氣候變化和人類活動(dòng)是區(qū)域生態(tài)環(huán)境變化的主要驅(qū)動(dòng)力[9-10],本研究中使用的氣象數(shù)據(jù)來自于中國氣象科學(xué)數(shù)據(jù)共享網(wǎng)(http://cdc. cma.gov.cn/)發(fā)布的中國地面氣候資料月值數(shù)據(jù)集及中國高分辨率地面氣象要素驅(qū)動(dòng)數(shù)據(jù)集[11-12].研究選取月平均氣溫與降水量兩個(gè)參數(shù),累加得到平均氣溫與降水量的逐年數(shù)據(jù),并通過MicroMet模型[13]進(jìn)行站點(diǎn)數(shù)據(jù)插值得到分辨率為1km×1km的空間柵格數(shù)據(jù),裁剪研究區(qū)域得到呼倫湖流域范圍的2000～2013年氣溫和降水量資料.模型首先采用距離加權(quán)插值方法將源數(shù)據(jù)插值到網(wǎng)格,并且將輸入的源數(shù)據(jù)的高程信息也按距離加權(quán)插值出一個(gè)虛擬的高程水準(zhǔn)面,再對(duì)參考網(wǎng)格與虛擬的高程水準(zhǔn)面的高程差進(jìn)行調(diào)整得到插值結(jié)果.氣溫的插值結(jié)果是用高程差與氣溫直減率來計(jì)算;降水量的插值則釆用了非線性的降雨—高程方程[14].

圖2　2000～2013年呼倫湖流域地表溫度逐月變化Fig.2　Monthly variation of temperature in the Hulun Lake watershed from the year 2000 to 2013

1.2水體面積提取方法

為克服利用分類后提取水體面積的方法和利用構(gòu)建水體特征指數(shù)的方法的不足,本研究采用水體指數(shù)動(dòng)態(tài)分析的方法通過將時(shí)間連續(xù)的水體指數(shù)計(jì)算結(jié)果進(jìn)行累計(jì),對(duì)累計(jì)水面觀測達(dá)到一定閾值后才認(rèn)為是有效豐水期水體面積,減小了因偶然波動(dòng)帶來的誤差[15].

呼倫湖位于中高緯度地區(qū),在冬季時(shí)易受積雪的影響,湖體難以與周圍地面明顯區(qū)分,且冬季一般為該流域的枯水期,故豐水期水體面積方法計(jì)算采用一年內(nèi)非結(jié)冰期的遙感影像進(jìn)行如下的動(dòng)態(tài)分析方法:在裁剪過后的呼倫湖主體湖區(qū)和流域的遙感影像范圍內(nèi),對(duì)每一景圖像利用第4波段(綠波段)和第6波段(中紅外波段)計(jì)算NDWI.公式如下:

取NDWI＞0的點(diǎn)即視為該點(diǎn)在該時(shí)間被淹沒一次,將24景圖像被淹沒的次數(shù)進(jìn)行累加,即取得每個(gè)像元在一年豐水期當(dāng)中被水覆蓋的次數(shù),淹沒次數(shù)≥10次的區(qū)域,視為當(dāng)年呼倫湖的豐水期水體范圍,根據(jù)得到的淹沒像元數(shù)目和像元大小計(jì)算水體面積.最后得到呼倫湖流域和呼倫湖主體湖區(qū)豐水期水體面積圖(圖3)和呼倫湖流域豐水期水體面積提取結(jié)果時(shí)序圖(圖4).

1.3水體面積與氣候要素相關(guān)性分析方法

呼倫湖流域豐水期水體面積與氣候要素相關(guān)分析采用偏相關(guān)分析的方法,通過去除另外的隨機(jī)變量的影響,來度量兩組隨機(jī)變量之間的關(guān)聯(lián)程度的方法.若隨機(jī)變量只有3組,則去除第3組變量的影響,來考慮1、2組變量之間的偏相關(guān)關(guān)系r12,3,可由公式(2)表示:

式中:r12表示變量1與變量2的相關(guān)性,如公式(3)所示,并以此類推.

式中:xi與yi分別表示變量1與變量2的個(gè)體.研究假設(shè)年均溫、年總降水量對(duì)豐水期水體面積變化均產(chǎn)生一定的作用,在分別研究年均溫、年總降水量這兩個(gè)變量與水體面積變化的驅(qū)動(dòng)力影響關(guān)系時(shí),偏相關(guān)分析可以只分析一個(gè)變量與水體面積變化之間的相關(guān)程度,而將另一個(gè)變量的影響剔除.

2　結(jié)果與分析

2.1呼倫湖水體面積變化分析

從圖3、圖4可以看出,2000～2012年期間,呼倫湖流域豐水期水體面積呈現(xiàn)出加速減小的趨勢.2012年呼倫湖流域豐水期水體面積比2000年減少了539km2.其中呼倫湖東北部和南部水域減少非常顯著,至2005年呼倫湖東北部水域已經(jīng)完全消失.2006年以后呼倫湖南部開始逐步縮小, 至2012年南部突出水域也趨于消失.相關(guān)研究綜合各種資料得到的水體面積結(jié)果也表明,呼倫湖的水體面積21世紀(jì)初(2000～2006年)由于干旱的影響,出現(xiàn)了快速的萎縮[5].同時(shí),本研究遙感監(jiān)測結(jié)果表明在呼倫湖周邊地區(qū),出現(xiàn)了一批面積較小的新增水域,主要集中在新巴爾虎左旗中南部地區(qū).呼倫湖主體湖區(qū)豐水期面積也呈現(xiàn)加速減小的趨勢,與流域豐水期水體面積減小相趨勢一致,主要體現(xiàn)在呼倫湖南部突出水域.2008年起,呼倫湖南部突出水域與主體湖區(qū)分離,形成南北兩小片獨(dú)立水域,并繼續(xù)縮小.2010年南部獨(dú)立水域的北片水域完全消失.

圖3　2000～2013年呼倫湖流域和湖區(qū)水體面積逐年變化Fig.3　Yearly variation of water area in the Hulun Lake watershed from the year 2000 to 2013

圖4　2000～2013年呼倫湖流域水體面積變化影像Fig.4　The images of water area variation in the Hulun Lake watershed from 2000 to 2013

2013年由于呼倫貝爾地區(qū)百年一遇的暴雨,呼倫湖流域水體和主體湖區(qū)豐水期面積均有所增長,體現(xiàn)在呼倫湖南部突出水域的面積部分恢復(fù),以及主體湖區(qū)周圍地區(qū)新巴爾虎左旗中部和新巴爾虎右旗中西部地區(qū)水體面積的恢復(fù).流域和主體湖區(qū)豐水期面積均恢復(fù)到2008年左右的水平.

2.2呼倫湖流域氣候要素變化分析

從圖5可以看出,2000～2013年期間,呼倫湖流域年均溫平均為2.1℃,呈現(xiàn)波動(dòng)中降低的趨勢,年均溫約以3年為一個(gè)波動(dòng)周期;2003、2006、2009、2012和2013年年均溫低于平均水平,分別為1.8、1.9、1.6、1.3和0.9℃;2007年年均溫達(dá)到14年來最高值(3.2℃),2013年年均溫達(dá)到14年來最低值(0.9℃);自2010～2013年,年均溫連續(xù)降低,降溫幅度平均為0.63℃/a.

2000～2013年期間,呼倫湖流域中國境內(nèi)部分年總降水量平均為321.4mm,呈現(xiàn)波動(dòng)中增加的趨勢,2001年總降水量達(dá)到14年來最低值(205.3mm),2013年總降水量達(dá)到14年來最高值(581.7mm);自2010～2013年,年總降水量連續(xù)增加,增加幅度平均為93.7mm/a,與相關(guān)研究得到的結(jié)論相一致[6].

圖5　2000～2013年呼倫湖流域年均氣溫、總降水量逐年變化Fig.5　Yearly variation of average annual temperature and precipitation in the Hulun Lake watershed from the year 2000 to 2013

2.3呼倫湖流域水體面積與氣候要素相關(guān)分析

從圖5可以看出,年均溫和年總降水量呈現(xiàn)較為明顯的負(fù)相關(guān)關(guān)系,在2013年年均溫達(dá)到14年來最低值的同時(shí),年總降水量達(dá)到14年來的最高值.

圖6　2000～2013年呼倫湖流域水體面積變化與年總降水量相關(guān)性分析Fig.6　Relationship between water area change and annual precipitation in the Hulun Lake watershed from 2000 to 2013

圖7　2000～2013年呼倫湖流域水體面積變化與溫度相關(guān)性分析Fig.7　Relationship between water area change and average annual temperature in the Hulun Lake watershed from 2000 to 2013

對(duì)年均溫(變量1)、年總降水量(變量2)和呼倫湖流域豐水期水體面積變化量(變量3)進(jìn)行回歸分析(圖6～),r2分別是:=0.51,=0.08,= 0.34.對(duì)年均溫(變量1)、年總降水量(變量2)和呼倫湖流域豐水期水體面積變化量(變量3)進(jìn)行偏相關(guān)分析,偏相關(guān)系數(shù)r分別如下:r12,3= -0.69, r13,2=-0.22,r23,1=0.56;由此可知:年均溫與年總降水量有明顯的負(fù)相關(guān)關(guān)系(r12,3=-0.69),年均溫較低的年份年總降水量往往較高,反之亦然;年均溫與呼倫湖流域豐水期水體面積變化有一定負(fù)相關(guān)關(guān)系,年均溫越高的年份呼倫湖流域豐水期水體面積減少越明顯,但關(guān)系并不顯著(r13,2=-0.22);年總降水量與呼倫湖豐水期流域水體面積變化有明顯的正相關(guān)關(guān)系,年總降水量越多的年份呼倫湖豐水期水體面積減少越不明顯,甚至有所增大;年總降水量是決定呼倫湖豐水期水體面積變化的主要影響因素(r23,1=0.56),其他研究綜合近50年氣候數(shù)據(jù)分析也表明,降水、溫度等氣象因子對(duì)流域豐水期水面積會(huì)產(chǎn)生影響,分析表明降水減少、溫度增加等氣候變化會(huì)造成相對(duì)蒸散提高,從而帶來流域水體減少.

圖8　2000～2013年呼倫湖流域年總降水量與年均氣溫相關(guān)性分析Fig.8　Relationship between annual precipitation and average annual temperature in the Hulun Lake watershed from 2000 to 2013

3　結(jié)論

監(jiān)測表明,2000～2012年期間,呼倫湖流域水體面積呈現(xiàn)出加速減小的趨勢,其中,2000～2012年期間呼倫湖主體湖區(qū)面積從2286km2減少至1773km2,減幅22.4%,2013年面積回升至1883km2,約為2008年左右的水平.呼倫湖主體湖區(qū)水體面積減少部分主要分布在湖體東北部和南部水域,2005年呼倫湖東北部水域已經(jīng)完全消失,2008年起,呼倫湖南部突出水域與主體湖區(qū)分離,形成南北兩小片獨(dú)立水域,并繼續(xù)縮小,至2012年趨于消失.2013年,主要恢復(fù)區(qū)域?yàn)楹魝惡喜克?分析呼倫湖主體湖區(qū)面積的變化與流域內(nèi)氣溫和降水的關(guān)系,年均溫與呼倫湖流域水體面積變化有一定負(fù)相關(guān)關(guān)系,但并不顯著;年總降水量與呼倫湖流域水體面積變化有明顯的正相關(guān)關(guān)系.因此,年總降水量是決定呼倫湖水體面積變化的主要影響因素,尤其是2013年的湖體面積增加主要原因?yàn)榻邓蛩?

參考文獻(xiàn)：

[1] 陳靜生,周家義.中國水環(huán)境重金屬研究 [M]. 北京:中國環(huán)境科學(xué)出版社, 1992.

[2] 顧潤源,李思慧,趙慧穎,等.呼倫湖流域徑流對(duì)氣候變化的響應(yīng)[J]. 生態(tài)學(xué)雜志, 2012,31(6):1517-1524.

[3] 王蘇民,竇鴻身.中國湖泊志 [M]. 北京:科學(xué)出版社, 1998.

[4] 孫標(biāo).基于空間信息技術(shù)的呼倫湖水量動(dòng)態(tài)演化研究 [D].內(nèi)蒙古農(nóng)業(yè)大學(xué), 2010.

[5] 趙慧穎,烏力吉,郝文俊.氣候變化對(duì)呼倫湖濕地及其周邊地區(qū)生態(tài)環(huán)境演變的影響 [J]. 生態(tài)學(xué)報(bào), 2008,28(3):1064-1071.

[6] 白美蘭,郝潤全,沈建國.近46a氣候變化對(duì)呼倫湖區(qū)域生態(tài)環(huán)境的影響 [J]. 中國沙漠, 2008(1):101-107.

[7] 宋平,劉元波,劉燕春.陸地水體參數(shù)的衛(wèi)星遙感反演研究進(jìn)展 [J]. 地球科學(xué)進(jìn)展, 2011,26(7):731-740.

[8] Vermote E F, Vermeulen A. MODIS Algorithm Technical Background Document - Atmospheric Correction Algorithm: Spectral Reflectance (MOD09) Version4.0 [EB/OL].

[9] 周沙,黃躍飛,王光謙.黑河流域中游地區(qū)生態(tài)環(huán)境變化特征及驅(qū)動(dòng)力 [J]. 中國環(huán)境科學(xué), 2014,34(3):766-773.

[10] 王強(qiáng),張勃,戴聲佩,等.三北防護(hù)林工程區(qū)植被覆蓋變化與影響因子分析 [J]. 中國環(huán)境科學(xué), 2012,32(7):1302-1308.

[11] Yang K, He J, Tang W, et al. On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan Plateau [J]. Agricultural and Forest Meteorology, 2010,150:38-46.

[12] Chen Y, Yang K, He J, et al. Improving land surface temperature modeling for dry land of China [J]. Journal of Geophysical Research: Atmospheres, 2011,116(D20):1-15.

[13] Glen E L, Kelly E. A meteorological distribution system for high-resolution terrestrial modeling (MicroMet) [J]. Journal of Hydrometeorology, 2006,7(2):217-234.

[14] Peter E T, Steven W R, White M A. Generating surfaces of daily meteorological variables over large regions of complex terrain [J]. Journal of Hydrology, 1997,190(3):214-251.

[15] Jain S K, Singh R D, Jain M K, et al. Delineation of flood-prone areas using remote sensing techniques [J]. Water Resources Management, 2005,19(4):333-347.

Study on dynamic change of hulun lake water area and climate driving force analysis.

WAN Hua-wei1, KANG Jun2,3, GAO Shuai2, SHEN Wen-ming1(1.Satellite Environmental Application Center, Ministry of Environmental Protection, Bejing 100029, China；2.Institute of Remote sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China；3.College of Resources and Environment, University of Chinese Academy of Science, Beijing 100049, China). China Environmental Science, 2016,36(3)：894～898

Abstract：Water area dynamic change of the Hulun lake was analyzed from 2000 to 2013 using long time-series MODIS data and water index dynamic analysis method. The driving force of the change was also analyzed by combining meteorological data. The preliminary results showed that the water area of the lake decreased from 2286km2in 2000 to 1773km2in 2012 and the decreasing rate was 22.4%. The drastic changes mostly happened in the northeast and south of the lake. Due to the great increase of precipitation in 2013, the water area restored to 1883km2and the main growth happened in the south of the lake. The driving force analysis showed the variation of water area correlated negatively with an average annual temperature and positively with the annual total precipitation. However, the level of significance was for the precipitation higher than temperature.

Key words：Hulun Lake；water area；time series；climate

作者簡介：萬華偉(1981-),女,山東東阿人,正高級(jí)工程師,博士,主要從事定量遙感應(yīng)用、生態(tài)系統(tǒng)與生物多樣性遙感監(jiān)測與評(píng)價(jià)相關(guān)研究.發(fā)表論文20余篇.

基金項(xiàng)目：國家“863”項(xiàng)目(2012AA12A310);國家高分重大科技專項(xiàng)(05-Y30B02-9001-13/15-9)

收稿日期：2015-06-23

中圖分類號(hào)：X524,TP751

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-6923(2016)03-0894-05