1988—2017年洪澤湖自然湖面變化特征及其影響因素

趙玉

摘要：以L(fǎng)andsat TM/ETM+遙感影像作為數(shù)據(jù)源，基于GIS和RS技術(shù)，利用年變化率和湖泊變化強(qiáng)度指數(shù)等方法定量分析1988—2017年洪澤湖湖面變化趨勢(shì)及特征，通過(guò)收集1988—2018年洪澤湖地區(qū)的氣象數(shù)據(jù)以及人工圍養(yǎng)數(shù)據(jù)，結(jié)合當(dāng)?shù)貙?shí)際深入研究湖面變化主要原因。①1988—2017年，洪澤湖自然湖面整體上呈萎縮-平穩(wěn)-萎縮的態(tài)勢(shì)，近30年共減少125.81 km2，相對(duì)萎縮速率為9.88%，強(qiáng)度變化指數(shù)（Q）為-0.34。相比而言，第三階段（2009—2017年）減少較為顯著，第一階段（1988—2000年）次之，第二階段（2000—2009年）較穩(wěn)定，其Q分別為-0.61、-0.24和-0.13;從區(qū)域位置上來(lái)看，洪澤湖面積變化主要集中在湖區(qū)西部，其次是東北部。②氣溫與降水量因子均與湖面變化呈線(xiàn)性相關(guān)，氣溫在1988—2018年總體升高，促使蒸發(fā)量增大;其降水量的變化也導(dǎo)致湖泊補(bǔ)給量下降。③洪澤湖沿岸的養(yǎng)殖業(yè)造成自然湖面萎縮，1988—2017年人工養(yǎng)殖面積擴(kuò)大了20倍，直接導(dǎo)致自然湖面損失5.0%。④湖泊階段變化中，政策變化以及洪澇災(zāi)害也是影響洪澤湖自然湖面變化的驅(qū)動(dòng)因子。1988—2017年洪澤湖自然湖面基本呈逐步萎縮的態(tài)勢(shì)，氣象因素、政策變化以及洪澇災(zāi)害情況導(dǎo)致湖面的萎縮，但圍網(wǎng)養(yǎng)殖是導(dǎo)致其變化的主要原因。

關(guān)鍵詞：洪澤湖;自然面積變化;驅(qū)動(dòng)因子;人工圍養(yǎng)

中圖分類(lèi)號(hào)：X52? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2020）18-0048-06

DOI：10.14088/j.cnki.issn0439-8114.2020.18.010

Characteristics and influencing factors of natural lake

level change of Hongze lake in 1988—2017

ZHAO Yu

（College of Public Administration， Hohai University，Nanjing? 210098，China）

Abstract： Using Landsat TM/ETM+ remote sensing image as data source， based on GIS and RS technology， using annual change rate and lake change intensity index and other methods to quantitatively analyze the change trend and characteristics of Hongze lake in 1988—2017. By collecting the meteorological data and artificial conservation data of Hongze lake in 1988—2018， combined with the local actual situation， the main causes of area change were studied in depth. ①From 1988 to 2017， the natural lake surface of Hongze lake showed a trend of shrinking-stable-shrinking. In the past 30 years， it decreased 125.81 km2， the relative shrinking rate was 9.88%， and the intensity change index （Q） was -0.34. In comparison， the third stage （2009—2017） decreased significantly， followed by the first stage （1988—2000）， and the second stage （2000—2009）， which was relatively stable， with Q of -0.61， -0.24 and -0.13， respectively. In terms of regional location， the change of Hongze lake area was mainly concentrated in the west of the lake area， followed by the northeast. ②Both temperature and precipitation factors are linearly related to the change of lake level， and the overall increase of temperature in 1988—2018 promotes the increase of evaporation; The change of precipitation also leads to the decrease oflake supply. ③From 1988 to 2017， the area of artificial aquaculture expanded 20 times， which directly led to the loss of 5.0% of the natural lake. ④During the change of lake stage， policy change and flood disaster are also the driving factors that affect the change of the natural lake surface of Hongze lake. From 1988 to 2017， the natural lake surface of Hongze lake basically showed a trend of gradual shrinking. Meteorological factors， policy changes and flood and waterlogging disasters caused the lake surface to shrink， but purse seine farming was the main reason for its change.

Key words： Hongze lake; natural area change; influence factor; culture in enclosure

湖泊作為人類(lèi)發(fā)展必不可少的自然資源，其在調(diào)節(jié)徑流、保持水土及防洪灌溉等方面發(fā)揮著重大作用[1]。開(kāi)發(fā)、利用與保護(hù)湖泊資源已經(jīng)成為人類(lèi)生存與社會(huì)可持續(xù)發(fā)展的重要研究課題。隨著遙感技術(shù)的發(fā)展，利用該技術(shù)的宏觀性、動(dòng)態(tài)性和經(jīng)濟(jì)性?xún)?yōu)勢(shì)[2]，可以準(zhǔn)確客觀獲取湖面變化趨勢(shì)和特征，促進(jìn)湖泊流域的合理開(kāi)發(fā)、利用和保護(hù)。經(jīng)濟(jì)的快速發(fā)展以及人類(lèi)活動(dòng)的強(qiáng)烈干預(yù)加快了江湖阻隔，作為中國(guó)東部淺水湖較為集中區(qū)，長(zhǎng)江中下游湖泊面積持續(xù)萎縮，自然環(huán)境及水質(zhì)不斷惡化[3]。洪澤湖作為淮河中下游地區(qū)最大的攔洪蓄水平原湖泊型水庫(kù)，具有防洪灌溉、調(diào)節(jié)水量和水產(chǎn)養(yǎng)殖等綜合效能，對(duì)于流域內(nèi)經(jīng)濟(jì)發(fā)展和人民生活具有重大意義[4]。隨著人口的增長(zhǎng)和經(jīng)濟(jì)的發(fā)展，加之自然變化和人類(lèi)活動(dòng)，洪澤湖面積逐年減少。特別是由于人類(lèi)活動(dòng)導(dǎo)致自然湖面被堤壩隔絕、圍墾，用于養(yǎng)殖、造田，造成泥沙淤積，使自然湖面日益縮小，影響湖泊調(diào)節(jié)洪水的能力，加劇湖區(qū)洪澇災(zāi)害的發(fā)生，造成了湖泊生境與生物多樣性破壞，對(duì)洪澤湖水質(zhì)惡化有一定的影響[5]。

已有研究表明，20世紀(jì)30年代到20世紀(jì)末，70年間洪澤湖的湖泊面積共減少317.22 km2[6]，1979—1988年，由于人類(lèi)活動(dòng)洪澤湖西岸臨淮鎮(zhèn)附近的湖區(qū)大部分已經(jīng)變成了養(yǎng)殖池塘或者圍網(wǎng)養(yǎng)殖區(qū)[7]。人為原因?qū)е潞次s的同時(shí)，持續(xù)破壞湖泊生境，導(dǎo)致湖泊生態(tài)功能退化[8]。為了進(jìn)一步研究近30年洪澤湖面積變化現(xiàn)狀，本研究收集了洪澤湖地區(qū)1988—2017年不同時(shí)期的Landsat TM/ETM遙感數(shù)據(jù)，在遙感和地理信息系統(tǒng)的支持下，研究了近30年來(lái)自然湖面的變化狀況。并利用研究區(qū)內(nèi)人工圍養(yǎng)信息及氣象數(shù)據(jù)，結(jié)合相關(guān)文獻(xiàn)，分析了研究區(qū)變化的內(nèi)外原因，研究結(jié)果可為區(qū)域內(nèi)湖泊保護(hù)、協(xié)調(diào)區(qū)域湖泊資源開(kāi)發(fā)與生態(tài)環(huán)境保護(hù)提供一些科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

洪澤湖位于33°06′—33°40′N(xiāo)和118°10′—118°52′E，是淮河流域最大的湖泊。在江蘇省北部，京杭大運(yùn)河以西，地處淮河、沂河等水系的中下游（圖1），是中國(guó)的第四大淡水湖，也是重要的水上運(yùn)輸樞紐，具有防洪灌溉、調(diào)水運(yùn)輸、調(diào)節(jié)氣候以及保護(hù)生物多樣性等綜合利用功能[9]。洪澤湖屬于過(guò)水性湖泊，湖面隨水位而波動(dòng)。在自然區(qū)劃上，洪澤湖既屬于黃淮海平原和長(zhǎng)江中下游的過(guò)渡區(qū)，也處于中國(guó)北亞熱帶與南暖溫帶的過(guò)渡地帶，四季分明，季風(fēng)氣候顯著，植物區(qū)系豐富，區(qū)位優(yōu)勢(shì)明顯，在江淮乃至長(zhǎng)江中下游地區(qū)都是一個(gè)典型的濕地[10]。

1.2 數(shù)據(jù)處理與方法

1.2.1 數(shù)據(jù)源 4期Landsat TM/ETM+影像均來(lái)自于美國(guó)地質(zhì)勘探局（United States geological survey，USGS）。由于湖面存在一定的季節(jié)性變化，因此本研究選用的影像均在相近月份。研究區(qū)能夠獲取的比較完整的Landsat TM數(shù)據(jù)有4期，分別為1988年10月、2000年9月、2009年10月、2017年10月的Landsat TM影像（表1）。

氣象數(shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)（http：//data.cma.cn/），選用泗洪（58135）地面基準(zhǔn)站1988—2018年氣溫和降水參數(shù)數(shù)據(jù)，以此分析影響湖泊自然面積變化的自然因素。

1.2.2 數(shù)據(jù)處理 4期遙感影像均完好無(wú)損，空間分辨率為30 m。數(shù)據(jù)已由USGS經(jīng)過(guò)正射糾正。不同時(shí)期的影像沒(méi)有明顯的位置偏移，匹配程度都較高，不需要進(jìn)行幾何校正。TM影像采用波段 5、4、3合成假彩圖像，分別以紅、綠、藍(lán)3種顏色顯示，符合人類(lèi)感官習(xí)慣。在此種波段組合上，水體呈現(xiàn)藍(lán)色或淡藍(lán)色，但建筑物呈現(xiàn)紅色，植被呈現(xiàn)綠色，各地物有明顯的特征差別[11]。自然湖面與其他地物不同，容易區(qū)分。

1.2.3 研究方法 自然湖面、圍養(yǎng)區(qū)域邊界提取采用ArcGIS 10.2軟件，以人機(jī)交互方式通過(guò)目視解譯生成洪澤湖自然湖面和圍養(yǎng)區(qū)的矢量圖層數(shù)據(jù)[12]，分別建立4個(gè)時(shí)期面積的數(shù)據(jù)庫(kù)，以此為基礎(chǔ)，利用GIS的疊加分析功能，獲取不同時(shí)期洪澤湖變化面積。Landsat TM 影像以波段5、4、3組合，湖水顯示為藍(lán)色，可以明顯看出與周?chē)匚锊煌?，易于提取。而養(yǎng)殖區(qū)一般分布在湖區(qū)周?chē)伾c水體相似，從影像上看不易區(qū)分。但是人工圍養(yǎng)區(qū)域一般是規(guī)則排列，有整齊的圍墻線(xiàn)，而圍墻線(xiàn)的顏色在影像上容易識(shí)別。因此，以圍欄為界劃分圍養(yǎng)區(qū)和自然湖面。

針對(duì)以上處理不同時(shí)期的自然湖面邊界矢量數(shù)據(jù)，利用ArcGIS 10.2，快速地對(duì)不同年代湖面多邊形進(jìn)行處理，計(jì)算面積并快速地估算湖面擴(kuò)張或萎縮速率，以及湖泊強(qiáng)度變化指數(shù)，進(jìn)而有效地研究自然湖面變化趨勢(shì)。并結(jié)合研究區(qū)實(shí)際情況，分析氣象數(shù)據(jù)、人工圍養(yǎng)數(shù)據(jù)、政策變化以及洪澇災(zāi)害等情況，深入探討其變化驅(qū)動(dòng)機(jī)制。

研究引入年變化率[13]作為衡量洪澤湖自然湖面擴(kuò)張或萎縮的一個(gè)指標(biāo)，其計(jì)算公式為：

[B=S末-S初S初×1T×100%]? ? ? （1）

式中，B表示研究時(shí)段內(nèi)湖泊年變化率，如果B>0，表示自然湖面是增加的，反之，則表示自然湖面是減少的;S末、S初分別為湖泊研究期末和研究期初的自然面積;T為研究間隔時(shí)段年數(shù)。

通過(guò)構(gòu)建湖泊變化強(qiáng)度指數(shù)[13，14]定量分析自然湖面變化相對(duì)強(qiáng)度，其計(jì)算公式為：

[Q=ΔSaΔTa×S×100]? ? ? （2）

式中，Q表示自然湖面變化強(qiáng)度指數(shù);ΔSa表示a時(shí)間段內(nèi)自然湖面變化面積;ΔTa表示研究時(shí)間差，S表示湖泊總面積（本研究選取1988年數(shù)據(jù)）。

1.2.4 不確定性及誤差評(píng)估 湖泊面積誤差可視為一種測(cè)量誤差，由系統(tǒng)誤差和偶然誤差組成。由陰影和人類(lèi)理解差異引起的誤差（或錯(cuò)誤）可視為系統(tǒng)誤差，這種誤差可以通過(guò)選擇高質(zhì)量的遙感圖像和提高湖泊解譯精度來(lái)減少，但不能準(zhǔn)確估計(jì)[15]。而分辨率、光譜特征、云層覆蓋、主觀誤差屬于偶然誤差，可用誤差理論來(lái)估算：

[δ=δ12+δ22+δ32…+δn2]? ? （3）

式中，δ1，δ2，δ3，…，δn，為互不相關(guān)的誤差變量;δ為總體誤差。

人為主觀誤差及云層覆蓋差異造成的誤判無(wú)法獨(dú)立估算，可將其視為一個(gè)獨(dú)立誤差來(lái)對(duì)待，通過(guò)多次測(cè)量的方法來(lái)估算[16]。在本研究區(qū)，湖泊表面較清晰，只有2017年影像有少量云層影響。因此，對(duì)于湖泊面積而言，主要的影響因素來(lái)自圖像分辨率和不同操作者帶來(lái)的主觀誤差。

分辨率造成的誤差主要表現(xiàn)在湖區(qū)邊緣和養(yǎng)殖區(qū)的界限有誤差。因此，由分辨率造成的誤差可通過(guò)湖泊與陸地邊緣的像元數(shù)量計(jì)算獲得[17]。在實(shí)際計(jì)算中，可利用ArcGIS緩沖區(qū)（Buffer）計(jì)算獲得。4期影像的糾正精度均為15 m，1988、2000、2009、2017年的總精度分別為15.3、15.3、15.4、15.6 m。4期的誤差均小于2.5%，即精度在97%以上。不同操作者的解譯誤差可以用多次測(cè)量的標(biāo)準(zhǔn)誤差來(lái)表示，真實(shí)值可以用多次測(cè)量值的平均值來(lái)表示。通過(guò)不同操作者解譯獲得的湖泊面積，4期的誤差均小于2%，即精度在98%以上。

分辨率、主觀誤差屬于偶然誤差，基于誤差理論估算，4期誤差分別為0.57%、0.58%、0.56%和3.15%，均控制在3.5%之內(nèi)，即精度在96%以上，表明所有被使用的影像數(shù)據(jù)和所提取的邊界數(shù)據(jù)有效可靠。

2 結(jié)果與分析

2.1 自然湖面時(shí)空變化特征

通過(guò)目視解譯獲取的湖泊面積基礎(chǔ)數(shù)據(jù)，由表2可以看出，1988年洪澤湖自然面積為1 273.93 km2，2017年面積為1148.12 km2，1988—2017年面積共減少了125.81 km2，占湖泊總面積的9.88%，年變化率為-0.34%，湖泊變化強(qiáng)度指數(shù)為-0.34，整體呈逐年萎縮趨勢(shì)。洪澤湖自然面積變化具有明顯的階段性差異，第一階段（1988—2000年），自然湖面面積減少了39.88 km2，占湖泊總面積的3.13%，年變化率為-0.26%;湖泊變化強(qiáng)度指數(shù)為-0.24;第二階段（2000—2009年），面積減少了16.09 km2，占湖泊總面積的1.3%，年變化率為-0.14%，湖泊變化強(qiáng)度指數(shù)為-0.13，與第一階段相比，第二階段依然呈減少趨勢(shì)，減速僅為第一階段的1/2;第三階段（2009—2017年），湖泊面積減少了69.84 km2，占湖泊總面積的5.73%，年變化率為-0.72%，湖泊變化強(qiáng)度指數(shù)為-0.61，與第二階段相比，第三階段萎縮速率成倍加快，約是第二階段的5倍。30年間，萎縮速率呈增大-減小-再增大的趨勢(shì)。

圖2表明，1988—2017年洪澤湖面積變化主要集中于東北部及淮河入湖口的西南部，其在4個(gè)時(shí)期內(nèi)都存在明顯的變化，東南及西北部面積變化較不明顯。湖泊西部在1988—2000年變化幅度較大，呈不斷萎縮的趨勢(shì);湖泊南部和北部變化幅度相對(duì)較小，仍有部分區(qū)域，如淮河入湖口的西部部分區(qū)域呈現(xiàn)擴(kuò)大的趨勢(shì)，這與該水域?yàn)榻K泗洪洪澤湖濕地自然保護(hù)區(qū)核心區(qū)的一部分有關(guān)，由于這里建立了保護(hù)區(qū)，是洪澤湖區(qū)保護(hù)得最好的湖泊部分。從2000和2009年兩期面積對(duì)比可以看出，這10年來(lái)面積變化較穩(wěn)定，除東北部及西部一小部分萎縮外，基本處于重合狀態(tài)。2009—2017年湖泊萎縮較為顯著，湖區(qū)的南部、西部以及東北部均存在面積萎縮現(xiàn)象。

2.2 年平均氣溫與降水量變化

如圖3所示，1988—2018年洪澤湖地區(qū)氣溫變化較平緩，有輕微上升趨勢(shì)。升溫速率為0.4 ℃/10年，明顯高于全國(guó)平均水平[18]。其中，2007年溫度最高，達(dá)16.16 ℃，1991年最低，僅為14.14 ℃。將1988—2018年年均溫與洪澤湖自然面積變化進(jìn)行相關(guān)性分析（圖4），決定系數(shù)（R2）為0.888 8，在α=0.05水平上顯著相關(guān)，表明洪澤湖自然面積減少與氣溫升高密切相關(guān)。高溫會(huì)導(dǎo)致蒸發(fā)量增加，從而造成湖泊面積萎縮。

降水作為湖泊水量最直接的補(bǔ)給，近30年來(lái)，受季風(fēng)氣候影響洪澤湖地區(qū)降水較為豐沛。由圖3可知，1988—2018年洪澤湖地區(qū)降水量波動(dòng)變化，呈微弱減少趨勢(shì)，極大值為1 526 mm，出現(xiàn)在2003年，而在2004年達(dá)到了降水極小值，僅為521.6 mm，多年平均降水量為944.59 mm。其中，1988—2003年降水量呈增加趨勢(shì)，而在2003—2018年波動(dòng)最為明顯，呈現(xiàn)減少趨勢(shì)。由相關(guān)性分析（圖4）可知，R2=0.661 7在α=0.05水平上顯著，說(shuō)明降水量的變化對(duì)洪澤湖自然面積的變化有一定的影響。

2.3人工養(yǎng)殖面積的變化

基于實(shí)地考察，目前洪澤湖的周邊遍布圍養(yǎng)區(qū)，故人類(lèi)活動(dòng)的影響因素主要選取與湖區(qū)重合的圍養(yǎng)區(qū)變化進(jìn)行研究。對(duì)4期遙感影像解譯顯示，湖區(qū)周邊由于水產(chǎn)業(yè)發(fā)展，西部與東北部養(yǎng)殖范圍較大，從圖5可以看出，1988年養(yǎng)殖區(qū)較少，到2017年明顯呈倍數(shù)擴(kuò)張。此階段洪澤湖周邊漁民靠水吃水，大力發(fā)展養(yǎng)魚(yú)、養(yǎng)蟹等項(xiàng)目。

由于水產(chǎn)養(yǎng)殖業(yè)的迅猛發(fā)展及經(jīng)濟(jì)利益的驅(qū)使，在引水條件較好的耕地上普遍出現(xiàn)開(kāi)挖魚(yú)塘的現(xiàn)象。特別是洪澤湖東北部地區(qū)，由于水草豐美適合養(yǎng)蟹，該區(qū)域湖泊利用強(qiáng)度相對(duì)較大。有研究表明，1979—2002年洪澤湖天然濕地內(nèi)已有約6.5 km2的濕地轉(zhuǎn)變?yōu)轲B(yǎng)蟹塘、圍網(wǎng)養(yǎng)蟹區(qū)等各類(lèi)養(yǎng)殖區(qū)，特別是20世紀(jì)80年代末以來(lái)，許多天然濕地由于養(yǎng)蟹熱的興起被圍墾[7]，這與本研究變化趨勢(shì)相似。如圖6所示，1988年人工養(yǎng)殖剛興起，面積較小，僅有3.12 km2，2017年面積達(dá)66.5 km2，1988—2017年養(yǎng)殖面積平均增長(zhǎng)速率為2.19 km2/年。其中，1988—2000年的養(yǎng)殖面積平均增長(zhǎng)速率為2.30 km2/年，2000—2009年的平均增長(zhǎng)速率僅有1.20 km2/年，2009—2017年繼續(xù)加快，平均增長(zhǎng)速率為3.12 km2/年。洪澤湖周邊養(yǎng)殖具有明顯的擴(kuò)張趨勢(shì)，不斷由湖區(qū)邊界向湖區(qū)中心擴(kuò)散，西部和東北部尤為明顯，其他養(yǎng)殖區(qū)域呈零散分布。由于湖區(qū)周邊人工圍養(yǎng)面積的逐年擴(kuò)張，占用了部分湖泊面積，從而導(dǎo)致洪澤湖水域面積的減少。

由圖7可以看出，洪澤湖養(yǎng)殖面積增加趨勢(shì)與湖泊面積減少趨勢(shì)相一致，經(jīng)Pearson相關(guān)系數(shù)分析，人工圍養(yǎng)與自然湖面變化在α=0.05水平上顯著相關(guān)（P=-0.956）。1988—2017年自然湖面共減少125.81 km2，而人工養(yǎng)殖面積增加63.38 km2，直接導(dǎo)致湖面損失5.0%，占湖面面積萎縮總量的50.38%，相較于氣溫降水自然因素的影響，人工養(yǎng)殖是導(dǎo)致湖泊變化的主導(dǎo)因素。

2.4 政策及洪澇災(zāi)害情況

將湖泊變化分為3個(gè)階段，第一階段（1988—2000年），自然湖面減少較快，人工養(yǎng)殖面積增加迅速;第二階段（2000—2009年），相比前一階段，人工圍養(yǎng)變化率呈大幅下降，自然湖面基本不變;第三階段（2009—2017年），自然湖面與養(yǎng)殖面積變化都較快。第二階段相比第一、第三階段，變化較不顯著。

研究表明， 2000—2009年湖泊開(kāi)發(fā)相關(guān)政策發(fā)生變化，尤其關(guān)于人工圍養(yǎng)。2000年十五屆五中全會(huì)著重強(qiáng)調(diào)了水利基礎(chǔ)設(shè)施的建設(shè)，這對(duì)湖區(qū)周?chē)B(yǎng)殖起到一定限制作用;2004年，江蘇省出臺(tái)了第一個(gè)針對(duì)湖泊保護(hù)的地方性法規(guī)《江蘇省湖泊保護(hù)條例》，緊接著，2006年又批準(zhǔn)了《江蘇省省管湖泊保護(hù)規(guī)劃》，這對(duì)洪澤湖的治理有了明確的目標(biāo)和要求，有效控制了人工圍養(yǎng)擴(kuò)張。2008年，為更好地管理與保護(hù)洪澤湖，設(shè)立江蘇省洪澤湖水利工程管理處，在其內(nèi)部設(shè)立專(zhuān)門(mén)的湖泊管理部門(mén)[19]。一系列政策均是圍繞保護(hù)洪澤湖濕地生態(tài)系統(tǒng)而展開(kāi)，在一定程度上抑制了人工養(yǎng)殖面積的擴(kuò)張。這一階段，洪澤湖地區(qū)在2003、2006、2007年相繼發(fā)生洪澇災(zāi)害，造成較大損失。從降水量分析中可以發(fā)現(xiàn)，洪澤湖地區(qū)降水量波動(dòng)變化，尤其在2003年出現(xiàn)降水量極大值（1 526 mm），這一年洪澤湖遭遇特大的高水位洪澇災(zāi)害[20]，湖區(qū)大部分養(yǎng)殖戶(hù)嚴(yán)重虧本，紛紛撤離該地，養(yǎng)殖積極性受到重創(chuàng)，造成很多圍網(wǎng)被撤銷(xiāo)或漁民棄塘不養(yǎng)。部分漁民的圍網(wǎng)、漁具以及船只被洪水破壞后，無(wú)法進(jìn)行大型捕撈，只有外出打工來(lái)維持生活。2006、2007年又相繼發(fā)生洪澇，這對(duì)湖區(qū)周邊養(yǎng)殖業(yè)的發(fā)展產(chǎn)生了巨大沖擊，導(dǎo)致第二階段湖面萎縮減緩。經(jīng)過(guò)一段時(shí)間的休整，在第三階段（2009—2017年），圍養(yǎng)面積又繼續(xù)擴(kuò)張，特別是隨著水產(chǎn)養(yǎng)殖業(yè)的發(fā)展，湖區(qū)面積又繼續(xù)大幅減少。

3 討論

洪澤湖自然湖面基本呈萎縮態(tài)勢(shì)[6-10]，近30年（1988—2017年）洪澤湖面積變化速率明顯加快。本研究發(fā)現(xiàn)，圍網(wǎng)養(yǎng)殖是導(dǎo)致湖面萎縮的主要原因，氣象因素、政策變化以及洪澇災(zāi)害也是影響湖泊變化的重要因素，如何準(zhǔn)確區(qū)分這幾者的影響有待進(jìn)一步研究。目前國(guó)家已出臺(tái)許多相關(guān)政策和措施，但實(shí)施力度與作用時(shí)間不夠持久，因此在一定區(qū)位優(yōu)勢(shì)和政策優(yōu)勢(shì)基礎(chǔ)上[21]，應(yīng)繼續(xù)推進(jìn)“退漁還湖、退圩還湖”的治理，不斷完善《防洪法》《湖泊保護(hù)條例》《江蘇省河道管理?xiàng)l例》等法律法規(guī)，充分利用現(xiàn)有的水利設(shè)施，穩(wěn)定水位，注意蓄水與泄洪保護(hù)。同時(shí)加強(qiáng)洪澤湖周?chē)鷿O民的環(huán)保意識(shí)，在保證湖面不繼續(xù)萎縮、不污染環(huán)境的條件下，適當(dāng)發(fā)展養(yǎng)殖產(chǎn)業(yè)，以實(shí)現(xiàn)該地區(qū)的資源、環(huán)境與經(jīng)濟(jì)的和諧發(fā)展。

4 小結(jié)

1）近30年（1988—2017年），洪澤湖自然面積不斷減少，2017年面積最小，為1 148.12 km2，1988年面積最大，為 1 273.93 km2，共減少了125.81 km2，湖面萎縮主要分布在湖區(qū)西部和東北部。

2）將自然湖面變化分為3個(gè)階段，其年變化率分別為-0.26%、-0.14%、-0.72%，湖泊變化強(qiáng)度指數(shù)分別為-0.24、-0.13和-0.61。

3）自然因素的作用持續(xù)施加影響，氣溫變化趨勢(shì)與自然湖面萎縮呈近線(xiàn)性相關(guān)，而年均降水量對(duì)其影響較小;1988—2017年湖區(qū)周邊人工圍養(yǎng)面積不斷增加，占用自然湖面，在一定程度上導(dǎo)致湖面萎縮，相較于氣象因素，湖區(qū)內(nèi)的人工養(yǎng)殖對(duì)其影響更大。

4） 湖泊階段變化中，政策變化以及洪澇災(zāi)害也是重要的影響因子。2000—2009年湖泊開(kāi)發(fā)相關(guān)政策發(fā)生變化，尤其關(guān)于人工圍養(yǎng);洪澤湖地區(qū)在2003、2006、2007年多次發(fā)生洪澇災(zāi)害，對(duì)養(yǎng)殖業(yè)造成巨大沖擊，導(dǎo)致湖泊利用強(qiáng)度減小。

參考文獻(xiàn)：

[1] HOU X J，F(xiàn)ENG L，CHEN X L，et al. Dynamics of the wetland vegetation in large lakes of the Yangtze Plain in response to both fertilizer consumption and climatic changes[J].Journal of photogrammetry and remote sensing，2018，141：148-160.

[2] SHAO T T，SONG K S，DU J，et al. Retrieval of CDOM and DOC using in situ hyperspectral data： A case study for potable waters in northeast China[J]. Journal of the Indian society of remote sensing，2016，44（1）：77-89.

[3] 殷立瓊，江 南，楊英寶.基于遙感技術(shù)的太湖近15年面積動(dòng)態(tài)變化[J].湖泊科學(xué)，2005，17（2）：139-142.

[4] 夏智宏，劉 敏，王 苗，等.1990s以來(lái)氣候變化和人類(lèi)活動(dòng)對(duì)洪湖流域徑流影響的定量辨識(shí)[J].湖泊科學(xué)，2014，26（4）：515-521.

[5] 李 娜，施 坤，張運(yùn)林，等.基于MODIS影像的洪澤湖水生植被覆蓋時(shí)空變化特征及影響因素分析[J].環(huán)境科學(xué)，2019，40（10）：4487-4496.

[6] 范亞民，何華春，崔云霞，等.淮河中下游洪澤湖水域動(dòng)態(tài)變化研究[J].長(zhǎng)江流域資源與環(huán)境，2010，19（12）：1397-1403.

[7] 阮仁宗，夏 雙，陳 遠(yuǎn)，等.1979～2006年洪澤湖西岸臨淮鎮(zhèn)附近湖泊變化研究[J].濕地科學(xué)，2012，10（3）：344-349.

[8] 吳 翼，戴 蓉，徐勇峰，等.洪澤湖河湖交匯區(qū)土地利用時(shí)空動(dòng)態(tài)[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，40（4）：22-28.

[9] 王德成.洪澤縣志[M].北京：中國(guó)大百科全書(shū)出版社，1999.

[10] 王 慶，陳吉余.洪澤湖和淮河入洪澤湖河口的形成與演化[J].湖泊科學(xué)，1999，11（3）：237-244.

[11] 李炳亞，潘劍君，夏 超，等.基于空間位置關(guān)系的山地湖泊水體提取方法研究[J].遙感技術(shù)與應(yīng)用，2016，31（5）：983-993.

[12] 上官冬輝，劉時(shí)銀，丁永建，等.玉龍喀什河源區(qū)32年來(lái)冰川變化遙感監(jiān)測(cè)[J].地理學(xué)報(bào)，2004，59（6）：855-862.

[13] 萬(wàn) 瑋，肖鵬峰，馮學(xué)智，等.近30年來(lái)青藏高原羌塘地區(qū)東南部湖泊變化遙感分析[J].湖泊科學(xué)，2010，22（6）：874-881.

[14] SU W Z，GU C L，YANG G S，et al. Measuring the impact of urban sprawl on natural landscape pattern of the Western Taihu Lake watershed， China[J].Landscape and urban planning，2009，95（1）：61-67.

[15] RAJASEKHAR M，RAJU G S，RAJU R S，et al. Data on comparative studies of lineaments extraction from ASTER DEM， SRTM， and Cartosat for Jilledubanderu River basin， Anantapur district， A.P， India by using remote sensing and GIS[J].Data in brief，2018，20：1676-1682.

[16] 王 媛，吳立宗，許君利，等.1964—2010年青藏高原長(zhǎng)江源各拉丹冬地區(qū)冰川變化及其不確定性分析[J].冰川凍土，2013，35（2）：255-262.

[17] 劉時(shí)銀，姚曉軍，郭萬(wàn)欽，等.基于第二次冰川編目的中國(guó)冰川現(xiàn)狀[J].地理學(xué)報(bào)，2015，70（1）：3-16.

[18] 韓翠華，郝志新，鄭景云.1951—2010年中國(guó)氣溫變化分區(qū)及其區(qū)域特征[J].地理科學(xué)進(jìn)展，2013，32（6）：887-896.

[19] 江蘇省水利廳.江蘇省洪澤湖保護(hù)規(guī)劃[Z].2006.

[20] 葉正偉.洪澤湖流域洪澇災(zāi)害的成災(zāi)機(jī)理分析與探討[J].水文，2006，26（4）：85-87，42.

[21] 趙一晗.洪澤湖綜合治理與保護(hù)的調(diào)查和思考[J].治淮，2018（11）：61-62.