艾爾肯·圖爾蓀， 玉素甫江·如素力,2， 崔一爽， 喀迪阿依·阿力木，米日阿依·買土地

（1.新疆師范大學(xué)地理科學(xué)與旅游學(xué)院/流域信息集成與生態(tài)安全實驗室，新疆 烏魯木齊 830054；2.新疆師范大學(xué)新疆干旱區(qū)湖泊環(huán)境與資源重點實驗室，新疆 烏魯木齊 830054）

湖泊作為水圈，陸地水循環(huán)和全球水平衡的重要組成部分，在全球氣候、生態(tài)環(huán)境變化過程中起著重要作用［1-3］。湖泊的水循環(huán)過程不僅對當(dāng)?shù)氐臍夂蜃兓鹬匾饔茫遗c人類活動也有緊密地關(guān)系，尤其是在干旱區(qū)域水循環(huán)和人類活動中作用更加明顯。新疆大部分湖泊均勻分布在高山與綠洲區(qū)域，這些湖泊為人類生產(chǎn)活動與生態(tài)環(huán)境平衡提供重要的水資源［4-8］。湖冰是記錄湖泊環(huán)境變化過程的重要指示因子，湖冰物候、厚度和形狀跟氣候變化過程密切相關(guān)，湖冰的生消和持續(xù)時間直接受到氣候變化的影響，所以在一些情況下湖冰物候比氣溫可以更精準(zhǔn)地記錄當(dāng)?shù)貐^(qū)域氣候的變化過程［9-11］。湖冰能夠為人類生活提供娛樂活動場所，在博斯騰湖舉辦的冰雪節(jié)每年可以接待幾萬名國內(nèi)外旅客，并且每年有幾千名旅客沿著博斯騰湖參加馬拉松比賽，能夠推動當(dāng)?shù)亟?jīng)濟發(fā)展。此外，湖冰對野生動物的棲息空間也起著重要作用，博斯騰湖是中國西北干旱區(qū)域的野生鳥類棲息地之一，在冬天隨著氣溫下降，博斯騰湖會迎來數(shù)萬只水鳥棲息越冬［12-15］。因此，湖泊結(jié)冰和消融時間的精確提取對人類活動和野生動物的棲息空間起著重要作用［16-18］。

在記錄湖冰物候的研究過程中主要研究方法是通過光譜遙感和微波遙感數(shù)據(jù)，根據(jù)湖冰和湖水的反射率特征提取湖冰。在目前國內(nèi)外研究中自動提取的湖冰物候研究方法較少，如Cai 等利用MODIS 產(chǎn)品為數(shù)據(jù)來源通過多種遙感算法監(jiān)測湖冰物候［11］；此外，波段閾值法也廣泛應(yīng)用于湖冰研究當(dāng)中［19-20］，如Latifovic 等［5］利用AVHRR（Advanced very high resolution radiometer）數(shù)據(jù)應(yīng)用波段閾值法提取了1985—2004年加拿大北部湖泊的湖冰物候，研究結(jié)果指出利用AVHRR數(shù)據(jù)提取的湖冰物候跟實測數(shù)據(jù)有較強的相關(guān)性。鞏志［21］利用MODIS 數(shù)據(jù)（MYD10A2）分析了2001—2010年青藏高原共59個湖泊的湖冰物候特征。MODIS和AVHRR等光譜遙感數(shù)據(jù)有多光譜和高時空分辨率優(yōu)勢，由于在冬天受到云、雪等因素的干擾，光譜遙感數(shù)據(jù)的湖冰提取精度受到一些影響，因此有些學(xué)者通過微波遙感來提取湖冰物候［22-23］。Dibike 等［16］通過RADAR?SAT-2數(shù)據(jù)提取了加拿大北部湖冰物候特征。諸多研究結(jié)果可以看出，雖然通過微波遙感提取湖冰時可以減少雪、云等因素的干擾，但由于微波遙感的時間分辨率和空間分辨率較低，所以對新疆一些湖泊數(shù)量稀少且面積較小的區(qū)域提取湖冰物候的結(jié)果可能不佳［24-26］。Magnuson等［1］研究結(jié)果指出，1846—1995 年之間北半球諸多湖泊和河流的開始結(jié)冰日出現(xiàn)明顯地推遲，變化趨勢為5.80 d·(100a)-1，湖冰的開始融化日也顯著提前，變化趨勢為6.50 d·(100a)-1，總體上北半球湖冰覆蓋期明顯地縮短。Du 等［9］通過微波遙感數(shù)據(jù)分析結(jié)果指出，在2002—2015年北半球71 個湖泊中共43 個湖泊（約60.6%）的湖冰覆蓋期有明顯地縮短。Latifovic 等［5］分析了從1950—2004年加拿大湖冰物候狀況，研究結(jié)果指出湖冰的開始結(jié)冰日呈現(xiàn)提前趨勢，為0.12 d·a-1；湖冰的開始融化日呈現(xiàn)提前趨勢，為0.18 d·a-1。王智穎等［26］用MODIS 數(shù)據(jù)對2000—2015 年青藏高原上面積大于100 km2的湖泊湖冰物候進(jìn)行了分析，研究結(jié)果指出青藏高原北部湖區(qū)開始凍結(jié)期較早，而完全融化期較晚，湖冰完全封凍期較長；青藏高原南部湖區(qū)開始凍結(jié)時間較晚，完全消融時間早，湖冰持續(xù)時間較短，且冰情變化趨勢主要受溫度、風(fēng)速等自然要素的影響。

本文以新疆10個大型湖泊為例（薩利吉勒干南庫勒湖面積小于100 km2），通過250 m 空間分辨率的MODIS 數(shù)據(jù)研究了新疆大型湖泊長時間序列湖冰物候及其變化模式，同時分析了湖泊發(fā)育過程和湖泊要素變化。本文利用光學(xué)影像、遙感反演、算法模型、統(tǒng)計分析等方法闡明了新疆大型湖泊湖冰物候變化特征，揭示了氣候變化背景下的湖冰物候變化模式，為干旱區(qū)域水循環(huán)與氣候變化奠定了理論和研究技術(shù)基礎(chǔ)，同時也為新疆大型湖泊治理、湖泊生態(tài)環(huán)境保護(hù)提供了方法指導(dǎo)。

1 研究區(qū)概況

新疆地處中國西北干旱區(qū)，屬溫帶大陸性氣候，具有獨特的山地-綠洲-荒漠生態(tài)系統(tǒng)，面積約占全國陸地總面積的1/6，是“一帶一路”經(jīng)濟帶的核心區(qū)。新疆位于中國北部的阿爾泰山和南部的昆侖山之間，天山橫穿該地區(qū)，此外，塔里木盆地和準(zhǔn)噶爾盆地位于3 座山脈之間，形成了“三山夾兩盆”的獨特地形（圖1）。新疆邊緣由積雪覆蓋的山脈圍繞著盆地，在盆地內(nèi)廣泛分布沙漠和綠洲，這種奇特地貌類型使新疆形成大量高山冰川圍繞沙漠的獨特景觀。由于新疆氣候較干燥、降水量稀少、蒸散發(fā)高于降水量，新疆大部分區(qū)域植被覆蓋率較低，生態(tài)系統(tǒng)較薄弱，是氣候變化與人類活動較為敏感的區(qū)域。天山是塔里木河的發(fā)源地同時也是塔里木盆地的主要水資源的源地，對綠洲的生態(tài)環(huán)境和人們生產(chǎn)活動有非常重要的作用［11］。

圖1 研究區(qū)概況Fig.1 Overview of the study area

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

2.1 數(shù)據(jù)來源

本文利用MOD09GQ（Terra）和MYD9GQ（Aqua）空間分辨率為250 m 的2 種MODIS 數(shù)據(jù)集。MOD09GQ 數(shù)據(jù)共有2 個地表反射率波段，覆蓋0.620～0.876 μm的光譜范圍，MOD09GQ的可用時間范圍是2000 年8 月1 日至今一直持續(xù)觀測地球，MYD9GQ的可用時間范圍為2002年8月1日至今。

Landsat 衛(wèi)星是USGS 和NASA 的一項聯(lián)合探索計劃，從1972 年至今一直持續(xù)觀測地球，共有近50 a的遙感監(jiān)測記錄。Landsat衛(wèi)星大約每2周以30 m的分辨率對整個地球表面成像多光譜數(shù)據(jù)。本文選擇1978—2020 年云覆蓋率小于10%的Landsat 影像，獲取新疆大型湖泊的湖泊面積數(shù)據(jù)，并且提取了湖泊矢量數(shù)據(jù)（表1）。

表1 新疆主要湖泊的面積、海拔高度和礦化度Tab.1 Area,elevation and salinity of large lakes in Xinjiang

數(shù)字高程模型（DEM）是通過一定的地形海拔數(shù)據(jù)，能夠?qū)崿F(xiàn)對全球表面地形的信息化模擬。本文所選擇的DEM 數(shù)據(jù)來自中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http://www.resdc.cn），該數(shù)據(jù)集為基于最新的SRTM V4.1 數(shù)據(jù)經(jīng)整理拼接生成的90 m 的分省數(shù)據(jù)，數(shù)據(jù)采用WGS84橢球投影。

2.2 研究方法

2.2.1 水體提取改進(jìn)的歸一化水體指數(shù)（mND?WI）、歸一化的水體指數(shù)（NDWI）和強化植被指數(shù)（EVI）等指數(shù)結(jié)合提取了地表水并且有效地排除了植被等干擾因素。水體的分類標(biāo)準(zhǔn)是：如果同時滿足mNDWI＞NDVI、EVI＜0.1 和mNDWI＞EVI 并 且EVI＜0.1 這2 種條件的時候就把地表要素分類為水體，如果無法滿足這些條件時分為非水體［26-27］。各指數(shù)計算公式如下：

式中：Red、Green、Blue、Nir 和Swir 分別為Landsat 數(shù)據(jù)的紅色波段（630～630 nm）、綠色波段（525～600 nm）、藍(lán)色波段（450～515 nm）、近紅外波段（845～885 nm）和短波紅外波段（1560～1660 nm）。

2.2.2 湖冰物候提取基于冰和水在多光譜遙感的分類特征發(fā)現(xiàn)一般情況下水和冰在近紅外波段的反射率均小于冰和水在紅外的反射率，而植被、土壤等其他地表物在紅外波段的反射率均低于其在近紅外波段的反射率，因此根據(jù)近紅外和紅外2 個波段的反射率差值可以精確的提取湖冰，并能夠有效地消除部分大氣影響和系統(tǒng)誤差。在進(jìn)行波段相減后的影像，湖冰在近紅外的反射率較大，而湖水的反射率較小，因此波段差值結(jié)合閾值法可以很好地對湖冰進(jìn)行監(jiān)測，計算公式如下［28］：

式 中：Band1、Band2 分 別 為MODIS 近 紅 外 波 段（620～670 nm）與紅外波段（841～876 nm）的反射率；a、b為根據(jù)目視解譯或者直方圖方法獲得的湖冰監(jiān)測的閾值，即2波段反射率差值大于a且Band2反射率大于b時，該像元定義為湖冰。

為便于自動提取湖冰不同狀態(tài)的時間節(jié)點信息，本文將湖泊開始凍結(jié)日期定義為下半年湖冰比例（即湖冰面積與湖泊面積的比值）為10%的日期；湖泊完全凍結(jié)日期定義為下半年湖冰比例為90%的日期；湖泊開始消融和完全消融日期則分別定義為上半年湖冰比例為90%、10%的日期，湖冰完全覆蓋期是記錄湖冰變化過程的物候要素，湖冰覆蓋期的變化能夠直接表示區(qū)域氣候變化特征，尤其對于新疆這樣的典型干旱區(qū)域來說，湖冰覆蓋期的變化趨勢能精確地反映氣候變化特征，計算方法如下［29］：

式中：FUS、FUE、BUS 和BUE 分別為湖泊開始凍結(jié)日、完全凍結(jié)日、開始消融日和完全消融日；LA 和IA分別為湖泊面積和湖冰面積。

3 結(jié)果與分析

3.1 新疆大型湖泊面積時空變化

新疆10 個湖泊在1970s 初期的總面積為5873.91 km2，到1990s初期為5263.71 km2，2000年為5952.38 km2，2011年為7446.94 km2。結(jié)合圖2可知，湖泊面積經(jīng)歷了“先大幅下降后快速上升”的過程。

圖2 1987—2020年新疆大型湖泊面積變化Fig.2 Area variation of large lakes in Xinjiang during 1987—2020

1987—2020 年分布在不同地理位置的10 個大型湖泊面積呈現(xiàn)顯著地空間差異特征。其中，艾比湖的面積雖有增加趨勢，但變化趨勢不一致，1990—2004年湖泊面積呈現(xiàn)波動性增加趨勢，2005年開始湖泊面積呈現(xiàn)減少趨勢，且2015年的湖泊面積最小，為1034.58 km2（圖3a）。博斯騰湖面積的變化趨勢較為復(fù)雜，1987—2000年湖泊面積呈顯著增加趨勢，由1987 年的947.48 km2擴張為2000 年的1125.34 km2，但2000 年開始湖泊面積開始退縮，至2013 年湖泊面積為925.27 km2，相比2000 年湖泊面積總共減少了200.03 km2（-17.80%）；2014湖泊面積開始增加，至2020 年達(dá)到1069.24 km2；雖然1994、2009年和2013年湖泊面積呈現(xiàn)減少趨勢，但總體上湖泊面積呈增加趨勢（圖3b）。分布在新疆北部的賽里木湖、烏倫古湖和吉力湖的面積總體上無顯著變化，賽里木湖的面積由1989 年的455.89 km2增加到2020 年的461.01 km2，1989—2020 年出現(xiàn)5.12 km2的面積波動；1992—2020 年烏倫古湖的面積增加了13.25 km2，從空間變化趨勢可以看出湖泊東南湖岸呈減少趨勢；吉力湖的面積1992—2020年沒有發(fā)生顯著變化，總體上保持穩(wěn)定趨勢（圖3c）。阿其克庫勒湖的面積呈波動性增加趨勢（圖3d），湖泊面積由1992 年的381.53 km2增加到2020 年的569.07 km2，其中2007—2014 年的面積增速最大，總體上1992—2020 年湖泊面積一共增加了187.54 km2（49.15%）。阿牙克庫木湖的面積總體上保持增加趨勢（圖3e），由1992年的600.04 km2擴張到2020年的1137.76 km2；1987—2020年阿牙克庫木湖面積一共增加了537.72 km2，其中2001—2004 年之間阿牙克庫木湖的面積增加更為突出，湖泊面積一共增加了228.77 km2。

圖3 1987—2020年新疆大型湖泊面積空間變化Fig.3 Area spatial variation of large lakes in Xinjiang during 1987—2020

3.2 湖冰物候時空變化

圖4 為2000—2019 年新疆大型湖泊湖冰物候模式的空間分布。將新疆大型湖泊的開始凍結(jié)日、湖冰完全覆蓋期、湖冰開始消融日和湖冰完全消融日的變化趨勢組合成一個完整的柱狀圖來表達(dá)，在整個柱狀圖中每個柱子的朝向表示湖冰物候變化的提前和推遲，朝上時為推遲的趨勢，朝下時為提前的趨勢。新疆大型湖泊的開始凍結(jié)日呈現(xiàn)提前和推遲的2 種變化趨勢，在研究期間開始凍結(jié)日呈現(xiàn)推遲趨勢的湖泊分別為博斯騰湖、賽里木湖、艾比湖、吉力湖、烏倫古湖、薩利吉勒干南庫勒湖和鯨魚湖，且大部分湖泊的開始凍結(jié)日推遲趨勢在0.51～1.53 d·a-1之間；開始凍結(jié)日呈現(xiàn)提前趨勢的湖泊有3個，分別為阿牙克庫木湖（變化趨勢為-1.04 d·a-1）、阿克賽欽湖（-0.41 d·a-1）、阿其克庫勒湖（-0.31 d·a-1），主要位于新疆西南部地區(qū)，占總湖泊總數(shù)的30%。近20 a新疆共有6個湖泊完全消融日呈現(xiàn)提前趨勢，主要為分布于新疆中北部的博斯騰湖、艾比湖、賽里木湖、吉力湖、薩利吉勒干南庫勒湖和烏倫古湖，其變化趨勢分別為-0.57 d·a-1、-0.61 d·a-1、-0.52 d·a-1、-0.78 d·a-1、-0.28 d·a-1和-0.36 d·a-1；完全消融日呈現(xiàn)推遲趨勢的湖泊主要分布于昆侖山高原區(qū)域，其中阿牙克庫木湖、鯨魚湖的湖冰完全消融日的推遲趨勢較顯著。湖冰完全覆蓋期是記錄湖冰變化過程的物候要素，湖冰覆蓋期的變化能夠直接表示區(qū)域氣候變化特征，尤其對于新疆這樣的典型干旱區(qū)域來說，湖冰完全覆蓋期的變化趨勢能精確地反映氣候變化特征。新疆大部分湖泊的湖冰完全覆蓋期呈現(xiàn)明顯的縮短，其中艾比湖、吉力湖、博斯騰湖等湖泊的湖冰完全覆蓋期縮短較為顯著，變化趨勢分別為-1.76 d·a-1，-2.13 d·a-1、-0.81 d·a-1，這些湖泊主要分布在新疆中北部；湖冰完全覆蓋期延長的湖泊有3個，分別為阿牙克庫木湖、阿其克庫勒湖和鯨魚湖，變化趨勢分別為3.51 d·a-1、1.54 d·a-1、1.37 d·a-1，這些湖泊均勻分布在昆侖山高原北翼。

圖4 2000—2019年新疆湖冰物候變化模式的空間分布Fig.4 Spatial distribution of lake ice phenology pattern in Xinjiang during 2000—2019

3.3 湖冰完全覆蓋時期年際變化

湖冰完全覆蓋時期是湖冰物候的重要參數(shù)，能夠直接表示湖冰物候的變化過程。圖5 為2000—2019 年新疆大型湖泊湖冰完全覆蓋期的年際變化。從圖中可知，近20 a 雖然一部分湖泊的湖冰完全覆蓋期呈現(xiàn)延長趨勢，但大部分湖泊的湖冰完全覆蓋期呈縮短趨勢，其中阿牙克庫木湖在2001、2002、2004年和2005年沒有完全的結(jié)冰。湖冰完全覆蓋期延長的湖泊主要有阿牙克庫木湖、阿其克庫勒湖和鯨魚湖，變化趨勢分別為3.51 d·a-1、1.54 d·a-1和1.37 d·a-1；除此之外，湖冰完全覆蓋期縮短的湖泊分別為博斯騰湖、艾比湖、烏倫古湖、吉力湖、賽里木湖、阿克賽欽湖和薩利吉勒干南庫勒湖，完全覆蓋期變化趨勢分別為-0.81 d·a-1、-1.76 d·a-1、-0.61 d·a-1、-2.31 d·a-1、-1.84 d·a-1、-0.17 d·a-1和-0.11 d·a-1。根據(jù)分析結(jié)果可知，地理位置分布在昆侖山高原的湖泊湖冰完全覆蓋期總體上呈增加趨勢，其中增加幅度最大的是阿牙克庫木湖；分布在平原區(qū)域的湖泊湖冰完全覆蓋期總體上呈減少趨勢，其中減少幅度最明顯的是吉力湖。

圖5 2000—2019年新疆大型湖泊湖冰完全覆蓋期年際變化Fig.5 Interannual variation of lake ice complete coverage of large lakes in Xinjiang during 2000—2019

4 討論

本文分析了2000—2019 年新疆大型湖泊的湖冰物候特征，湖冰物候特征是區(qū)域氣候變化靈敏指示器，張音等［30］研究指出近60 a新疆年均氣溫、年均降水量總體上均呈上升趨勢，但區(qū)域空間變化規(guī)律不明顯，因此本文對區(qū)域氣候變化研究具有重要的參考意義。在全球氣候變化背景下，很多研究結(jié)果得出在北半球大部分湖泊出現(xiàn)開始凍結(jié)日推遲和開始消融日提前的趨勢。Dibike 等［16］研究結(jié)果指出在1855—2004 年分布在北半球湖泊的完全凍結(jié)日和完全消融日變化趨勢分別為0.03～0.16 d·a-1、-0.05～0.19 d·a-1，湖冰完全覆蓋期變化趨勢為-0.07～-0.43 d·a-1。

圖6 為4 個新疆大型湖泊湖冰凍融的變化趨勢。由圖6a～c 所示，在2001 年12 月16 日阿牙克庫木湖湖面沒有出現(xiàn)冰蓋，但2010 年12 月16 日湖泊出現(xiàn)完全結(jié)冰，2020年12月2日湖泊大部分區(qū)域已出現(xiàn)冰層，從湖泊結(jié)冰情況可以看出，近20 a阿牙克庫木湖的開始凍結(jié)日表現(xiàn)為明顯的提前趨勢，這種變化趨勢引起了湖冰完全覆蓋期的延長。Cai 等［11］對新疆湖冰物候變化研究發(fā)現(xiàn)，新疆大部分湖泊湖冰覆蓋期出現(xiàn)顯著的縮短趨勢，但阿牙克庫木湖等有些湖泊的湖冰完全覆蓋期出現(xiàn)延長趨勢，這種趨勢與本文研究結(jié)果一致。從圖6d～f中可知，2000年6月10日鯨魚湖的湖冰完全消融，2012年6月15日鯨魚湖的湖面有湖冰覆蓋，2020年6月25日鯨魚湖大部分區(qū)域仍出現(xiàn)湖冰，從這種變化趨勢可以看出，近20 a 鯨魚湖的湖冰完全覆蓋期呈延長趨勢。圖6g～i為賽里木湖的開始凍結(jié)日變化特征，2001年1月2日賽里木湖開始出現(xiàn)湖冰，而2009年1月9日湖面開始凍結(jié)，但2020年1月13日賽里木湖才出現(xiàn)冰層，從這種變化趨勢可以看出，在研究期間賽里木湖的開始結(jié)冰日表現(xiàn)為推遲趨勢，因此這種變化趨勢縮短了賽里木湖的湖冰完全覆蓋期。圖6j～l顯示了吉力湖的湖冰凍融過程，2001 年6 月3 日吉力湖的湖冰完全消融，而2010年5月26日和2020年5月7 日湖冰消融日呈提前趨勢，通過3 期的湖冰消融日可以得知，近20 a吉力湖的湖冰完全消融日表現(xiàn)為提前趨勢，而這種變化趨勢引起了吉力湖的湖冰覆蓋期縮短。很多研究結(jié)果指出，湖泊的凍融情況主要受到湖泊自身特征變化的影響［26,28］，例如，2000—2019 年阿牙克庫木湖的開始凍結(jié)日呈顯著提前趨勢，主要原因與湖泊面積顯著擴張有關(guān)，而湖泊面積的擴張導(dǎo)致了湖泊礦化度下降，因此這種響應(yīng)關(guān)系引起了開始凍結(jié)日提前和湖冰覆蓋期的延長。秦啟勇等［31］研究結(jié)果指出，2000—2019年賽里木湖的湖冰開始凍結(jié)日呈提前趨勢，與本文的研究結(jié)果存在差異；而Cai 等［11］研究發(fā)現(xiàn)賽里木湖的凍結(jié)日呈推遲趨勢，湖冰完全覆蓋期呈縮短趨勢，與本文研究結(jié)果一致。

圖6 新疆湖冰物候變化模式Fig.6 Lake ice phenology pattern in Xinjiang

5 結(jié)論

本文基于MOD09GQ、MYD09GQ 等光學(xué)遙感數(shù)據(jù)，通過多波段閾值法提取了湖冰物候，在本研究提取的湖冰物候結(jié)果通過Sentinel-2數(shù)據(jù)進(jìn)行了精度驗證，此外還總結(jié)了氣候變化背景下新疆湖泊湖冰凍融過程并分析了湖冰物候變化模式。主要結(jié)論如下：

（1）通 過250 m 空 間 分 辨 率MOD09GQ、MYD09GQ數(shù)據(jù)集對新疆大型湖泊進(jìn)行了湖冰物候提取，過去研究通過8 d合成500 m分辨率的MODIS冰雪數(shù)據(jù)集對新疆湖冰物候進(jìn)行研究，而本研究基于更高時空分辨率的MODIS 數(shù)據(jù)集應(yīng)用波段閾值法精準(zhǔn)地提取了多種湖冰物候參數(shù)。

（2）湖冰完全覆蓋期是重要的湖冰參數(shù)，湖冰覆蓋期的延長或者縮短能夠直接表示區(qū)域氣候變化過程，新疆大部分湖泊湖冰完全覆蓋期表現(xiàn)為縮短趨勢，其中分布在新疆中北部的艾比湖、吉力湖和博斯騰湖等湖泊的湖冰完全覆蓋期縮短較為明顯，變化趨勢分別為-1.76 d·a-1，-2.13 d·a-1和-0.81 d·a-1。

（3）湖泊面積、湖泊礦化度等湖泊自身要素能夠直接影響湖泊凍融過程。因此在研究期間阿牙克庫木湖、鯨魚湖等咸水湖的湖泊面積表現(xiàn)為明顯擴張趨勢，而這種變化引起了這些湖泊礦化度的下降，因此阿牙克庫木湖、鯨魚湖等湖泊的開始凍結(jié)日表現(xiàn)為提前，而湖冰完全覆蓋期呈延長趨勢。

由于本研究缺乏實測湖冰物候數(shù)據(jù)，因此在研究過程中無法對遙感監(jiān)測的湖冰物候進(jìn)行精度評價，除此之外，為了近一步系統(tǒng)地探討湖冰物候的驅(qū)動因子機制，未來可以考慮凍土變化趨勢和湖泊所在區(qū)域的冰川面積等因素。