劉 力， 李 朋,2， 胡元平， 張 雅,2， 李琳靜， 陶 良， 朱晛亭， 潘龍克

(1.湖北省地質(zhì)調(diào)查院,湖北 武漢 430034; 2.古生物與地質(zhì)環(huán)境演化湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢 430205；3.武漢市測(cè)繪研究院,湖北 武漢 430022)

濕地被譽(yù)為“地球之腎”,與森林、海洋一起并稱為全球三大生態(tài)系統(tǒng)，是地球上水陸相互作用形成的特殊自然綜合體，是各種能流、物流和信息流的交換場(chǎng)所[1]。濕地在涵養(yǎng)水源、凈化水質(zhì)、蓄洪防旱、調(diào)節(jié)氣候、維護(hù)生物多樣性和提供生產(chǎn)、生活資源等方面發(fā)揮著重要作用[2-3]。近年來(lái),有關(guān)濕地面積變化或時(shí)空分布差異的研究備受關(guān)注,周婷等[4]從地理學(xué)時(shí)空視角分析了中國(guó)大陸31個(gè)省(市、自治區(qū))濕地2003—2013年的空間格局及其關(guān)聯(lián)性；范韋瑩等[5]利用2005、2010和2017年的土地利用數(shù)據(jù),采用形態(tài)學(xué)空間格局分析方法研究了神農(nóng)架大九湖濕地空間形態(tài)演變特征；王志成等[6]從遙感影像中提取膠州灣周邊地區(qū)1987—2017年的土地利用變化信息，分析了該區(qū)潮灘濕地時(shí)空特征及演變成因；王鐵良等[7]利用1985—2017年33景Landsat遙感影像,定量分析了遼河口濱海濕地海岸線的時(shí)空變化特征。遙感技術(shù)與GIS技術(shù)是濕地宏觀研究中數(shù)據(jù)獲取與量化分析的重要手段[8-9],具有觀測(cè)范圍廣、信息量大、獲取信息快、可比性強(qiáng)等特點(diǎn)。

湯遜湖曾是武漢市最大的原生生態(tài)湖[10],隨著城市的開(kāi)發(fā)建設(shè)與經(jīng)濟(jì)社會(huì)的發(fā)展,特別是1996年以來(lái)湯遜湖周邊產(chǎn)業(yè)園與開(kāi)發(fā)區(qū)的建設(shè),對(duì)湖泊的侵占速度愈來(lái)愈快,這些侵占行為不可避免地對(duì)湖泊的生態(tài)環(huán)境產(chǎn)生負(fù)面影響。開(kāi)展湯遜湖濕地變化及驅(qū)動(dòng)因素分析,了解湯遜湖濕地1980年以來(lái)土地利用格局及其變化狀況,可為制訂湯遜湖濕地保護(hù)方案和改善湯遜湖區(qū)域人居環(huán)境提供科學(xué)依據(jù)和決策支持,更好地服務(wù)于武漢市國(guó)際濕地城市建設(shè)。

1 研究區(qū)概況

湯遜湖位于武漢市東南部,跨越江夏區(qū)、洪山區(qū)和東湖新技術(shù)開(kāi)發(fā)區(qū),水域面積47.62 km2,岸線長(zhǎng)度122.8 km,平均水深2.32 m。以江夏大道為界,湯遜湖被分為東側(cè)內(nèi)湖和西側(cè)外湖(圖1)。湯遜湖屬于湯遜湖水系,該水系以湯遜湖、南湖、野芷湖、青菱湖等為主要承載水體,通過(guò)青菱河、巡司河與長(zhǎng)江連接,屬于典型的平原型水網(wǎng)結(jié)構(gòu)。平水期、枯水期時(shí)水系匯水經(jīng)由陳家山閘、?？陂l和解放閘外排,豐水期啟動(dòng)泵站泄洪排入長(zhǎng)江[11]。

湯遜湖流域?qū)儆诒眮啛釒駶?rùn)季風(fēng)氣候區(qū),雨量充沛,四季分明,日照時(shí)間長(zhǎng),氣候條件適宜居住。地貌屬鄂東南丘陵經(jīng)江漢平原東緣向大別山南麓低山丘陵過(guò)渡類型，地形屬于殘丘型河湖沖積平原,地勢(shì)平坦低洼。湯遜湖湖泊功能為調(diào)蓄、景觀娛樂(lè)和生態(tài)調(diào)節(jié),水功能區(qū)劃為開(kāi)發(fā)利用區(qū),水質(zhì)管理目標(biāo)為Ⅲ類標(biāo)準(zhǔn)。長(zhǎng)期以來(lái),由于攔湖筑汊、管理無(wú)序、開(kāi)發(fā)過(guò)度、保護(hù)不力等因素的疊加,導(dǎo)致出現(xiàn)水體污染、環(huán)境惡化、功能退化等突出問(wèn)題[12]。

2 研究方法

2.1 數(shù)據(jù)來(lái)源

本研究采用美國(guó)陸地衛(wèi)星Landsat3 MSS、Landsat5 TM、Landsat8 OLI影像數(shù)據(jù),均下載于https://earthexplorer.usgs.gov/(美國(guó)地質(zhì)調(diào)查局官網(wǎng))。1980、1990、2000、2009和2019年影像含云量較低(云量<5%)且日期相近,避免因季節(jié)不同而導(dǎo)致雨季和非雨季存在植被類型及生長(zhǎng)情況、水田作物收割情況、灌排水情況和水體水位的差異,提高了數(shù)據(jù)的準(zhǔn)確性和有效性。1998、2016和2020年洪水年分別選擇洪水期和枯水期含云量較低的影像(表1)。

表1 遙感影像參數(shù)表Table 1 Remote sensing image parameters table

2.2 研究方法

在使用遙感影像之前,先進(jìn)行圖像的配準(zhǔn)。由于不同來(lái)源的遙感數(shù)據(jù)的軌道、成像時(shí)間、成像機(jī)理、平臺(tái)、觀測(cè)角度等不同,幾何特征差異很大，因此需要改進(jìn)原始影像的幾何變形,生成一幅滿足圖形表達(dá)要求或者地圖投影要求的新圖像[13]。本研究主要結(jié)合基于對(duì)象的思想以及機(jī)器學(xué)習(xí)算法[14]提取湯遜湖湖泊范圍。

(1) 基于對(duì)象的遙感影像分割。對(duì)象是指原始影像經(jīng)過(guò)分割得到的由相同或者相似的光譜、形狀、紋理和結(jié)構(gòu)等特征的像元集合體所構(gòu)成的多邊形單元。本課題基于多尺度分割方法,綜合考慮遙感影像的波段權(quán)重、形狀因子、光滑度、緊致度、分割尺度等因素,實(shí)現(xiàn)對(duì)象生成。對(duì)象生成后,通過(guò)定義對(duì)象的特征為對(duì)象的分類和目標(biāo)信息的提取提供基礎(chǔ)。對(duì)象的特征主要包括光譜特征、形狀特征、紋理特征和拓?fù)潢P(guān)系等。

(2) 提取方法。支持向量機(jī)(Support Vector Machines,SVM)是20世紀(jì)90年代發(fā)展起來(lái)的一種建立在統(tǒng)計(jì)學(xué)習(xí)理論基礎(chǔ)之上的機(jī)器學(xué)習(xí)方法,逐漸成為機(jī)器學(xué)習(xí)領(lǐng)域中的研究熱點(diǎn)。其基本思想是求得一個(gè)最優(yōu)超平面,通過(guò)此超平面將兩類數(shù)據(jù)進(jìn)行正確地分類,使得兩類數(shù)據(jù)樣本點(diǎn)到此超平面的最小距離之和達(dá)到一個(gè)最大值[15]。

(3) 建立解譯標(biāo)志。結(jié)合濕地分類體系和湯遜湖濕地實(shí)際情況，將研究區(qū)分為湖泊濕地、人工濕地和沼澤濕地3種地類(表2)。對(duì)于5期多光譜遙感影像,使用Envi軟件進(jìn)行大氣校正和圖像裁剪,并進(jìn)行圖像拉伸和銳化；然后使用支持向量機(jī)法對(duì)影像數(shù)據(jù)進(jìn)行監(jiān)督分類,采用Majority分析對(duì)分類結(jié)果進(jìn)行小斑塊處理,通過(guò)目視解譯法對(duì)錯(cuò)分及漏分處進(jìn)行重新分類；最后采用ArcGIS 10.2制圖并提取地類信息,計(jì)算濕地面積、地類動(dòng)態(tài)度及轉(zhuǎn)移矩陣[16]。

表2 湯遜湖濕地分類表Table 2 Classification table of Tangxun Lake wetland

3 結(jié)果與討論

3.1 濕地面積時(shí)空變化

湯遜湖地區(qū)的濕地主要為人工濕地、湖泊濕地和沼澤濕地。人工濕地經(jīng)歷了先增加后減少的變化過(guò)程(表3、圖2)，由1980年的1.96 km2增加到2000年的17.55 km2,而后又減少為2019年的6.88 km2。1980—2000年主要因圍湖建庫(kù)活動(dòng)減少，使人工濕地面積不斷增加；而2000—2019年因城市化速度加快,人工濕地或被占用為建設(shè)用地，或被轉(zhuǎn)化為湖泊濕地,導(dǎo)致人工濕地面積不斷減少。湖泊濕地面積經(jīng)歷了先減少后增加的變化過(guò)程,由1980年的46.65 km2減少到2000年的37.54 km2,而后又增加到2019年的42.08 km2。1980—2000年由于東湖新技術(shù)開(kāi)發(fā)區(qū)、藏龍島產(chǎn)業(yè)園、廟山開(kāi)發(fā)區(qū)等大規(guī)模開(kāi)發(fā)建設(shè)使得湖泊濕地面積不斷減少；而2000年以后,武漢市更加重視湯遜湖自然生態(tài)系統(tǒng)保護(hù)，先后修建了湯遜湖濕地公園、當(dāng)代運(yùn)動(dòng)公園、譚鑫培濕地公園等，使環(huán)湯遜湖區(qū)域的濕地自然環(huán)境得到一定程度的改善，尤其是2002年頒布的《武漢市湖泊條例》及陸續(xù)開(kāi)展的各類綜合治理工程在很大程度上促進(jìn)了湯遜湖生態(tài)的恢復(fù)。沼澤濕地面積經(jīng)歷了整體減少的變化過(guò)程，由1998年的7.95 km2減少到2000年的4.62 km2，而后又減少到2019年的1.93 km2。湯遜湖目前已很難找到原生湖岸,大部分湖岸已被硬化，同時(shí)隨著環(huán)湯遜湖區(qū)域工業(yè)園區(qū)的擴(kuò)建及人口的不斷增多,工業(yè)和生活污水排放量也增多,致使湖水水質(zhì)惡化,打破了湯遜湖局部生態(tài)系統(tǒng)的平衡,沼澤濕地明顯減少。

表3 湯遜湖不同年份濕地面積統(tǒng)計(jì)表Table 3 Wetland area statistics of Tangxun Lake in different years

圖2 湯遜湖不同年份濕地面積變化圖Fig.2 Wetland area change map of Tangxun Lake in different years

3.2 周圍土地利用類型變化

湯遜湖地區(qū)濕地面積與周圍土地利用面積變化情況見(jiàn)表4和圖3。湯遜湖地區(qū)居民地在1980—2019年呈現(xiàn)逐年增長(zhǎng)趨勢(shì),2000年達(dá)到7.98 km2,2019年達(dá)到93.43 km2,尤其在2000—2009年為大幅增長(zhǎng)。事實(shí)上,2000年后隨著武漢市城鎮(zhèn)化速度加快，湯遜湖地區(qū)受到中心城區(qū)推力及自身吸引力的影響，逐漸由典型的郊區(qū)演變?yōu)橹匾鞘薪ㄔO(shè)區(qū)，造成建設(shè)用地逐年增加，尤其是2000—2009年增長(zhǎng)更為明顯。湯遜湖地區(qū)濕地面積總體為減少趨勢(shì),但在1980—2000年較穩(wěn)定,各濕地類型之間互相轉(zhuǎn)化，2000年后濕地面積為大幅減少趨勢(shì),主要受人工濕地面積大量減少的影響。植被覆蓋區(qū)包含了林地和草地,同樣因湯遜湖地區(qū)城市化建設(shè)影響而逐年減少。其他用地類型主要為裸地,其面積也呈逐年減少趨勢(shì),尤其在2000—2009年為大幅減少。

圖3 湯遜湖濕地及土地利用類型面積變化圖Fig.3 Changes of wetland and land use types in Tangxun Lake

表4 湯遜湖濕地及土地利用類型面積統(tǒng)計(jì)表Table 4 Statistics of wetland area and land use type area in Tangxun Lake

3.3 洪水期與枯水期濕地類型對(duì)比

湯遜湖洪水期與枯水期濕地類型面積變化情況見(jiàn)表5和圖4，相比枯水期,人工濕地在洪水期減少,湖泊濕地在洪水期相應(yīng)增加。

表5 湯遜湖洪水期與枯水期濕地面積統(tǒng)計(jì)表Table 5 Statistical table of wetland area of Tangxun Lake in flood seasons and dry seasons

圖4 湯遜湖洪水期與枯水期濕地面積變化圖Fig.4 Changes of wetland area of Tangxun Lake in flood seasons and dry seasons

3.4 驅(qū)動(dòng)因素分析

3.4.1自然因素

武漢市于1998、2016和2020年發(fā)生了三次大洪水。根據(jù)收集上述三個(gè)洪水年的月平均降雨量(表6),發(fā)現(xiàn)2016年降雨量最為豐富(1 816.1 mm),其次為1998年(1 729.2 mm)和2020年(1 722.0 mm)。從單個(gè)月份來(lái)看，7月份的降雨量最大(圖5)。

表6 武漢市三次洪水年月平均降雨量統(tǒng)計(jì)表Table 6 Monthly average rainfall table of three floods in Wuhan City

圖5 武漢市三次洪水年月平均降雨量直方圖Fig.5 The monthly average rainfall histogram of three floods in Wuhan City

為探討洪水期降雨量與濕地面積的變化關(guān)系,本文分析了降雨量及各類濕地面積的Pearson相關(guān)性，結(jié)果(表7)表明湯遜湖湖泊濕地面積與降雨量之間存在弱的負(fù)相關(guān)關(guān)系,這主要是由于降雨量增加會(huì)導(dǎo)致人為圈閉的人工濕地與湖泊水域連成一片而形成更廣闊的湖泊水域,人工濕地面積因此減少。

表7 湯遜湖濕地面積與洪水期降雨量的Pearson相關(guān)性統(tǒng)計(jì)表Table 7 Pearson correlation between wetland area and precipitation during flood seasons in Tangxun Lake

3.4.2人類活動(dòng)因素

武漢市濕地多分布于相對(duì)低洼地帶,其分布格局主要受干流河道及平原內(nèi)部起伏的微地貌控制。近年來(lái),隨著人口劇增、城市擴(kuò)展、盲目開(kāi)墾土地及經(jīng)濟(jì)社會(huì)的發(fā)展,天然濕地尤其是湖泊濕地面積大幅度減少,濕地保護(hù)已成為武漢城市規(guī)劃建設(shè)的一項(xiàng)重要課題[18]。

(1) 人口因素。城市存在的基本條件就是一定數(shù)量的人口聚集。人在促進(jìn)城市發(fā)展的同時(shí)，也給濕地生態(tài)環(huán)境帶來(lái)一定影響。從武漢市總?cè)丝谠鲩L(zhǎng)曲線(圖6)來(lái)看,從1980年至今,武漢市總?cè)丝谠鲩L(zhǎng)具四個(gè)階段：①1980—1989年,人口增長(zhǎng)較快,年平均增長(zhǎng)率為2%；②1990—1999年,人口增幅有所回落,年平均增長(zhǎng)率為1.1%；③2000—2010年,人口增幅明顯提高，2000—2005年平均增長(zhǎng)率為1.55%,2006—2010年平均增長(zhǎng)率4.17%,尤其是2006—2007年、2008—2009年增長(zhǎng)率達(dá)到7.5%；④2011—2019年,2011年常住人口總數(shù)突破1 000萬(wàn),人口增幅有所放緩，年平均增長(zhǎng)率為1.42%。從人口增長(zhǎng)過(guò)程來(lái)看，2000年以來(lái)人口增長(zhǎng)加快，與湯遜湖濕地面積開(kāi)始減少的趨勢(shì)相符合。

圖6 武漢市總?cè)丝谠鲩L(zhǎng)曲線圖(據(jù)參考文獻(xiàn)[17])Fig.6 Growth curve of Wuhan population

(2) 土地利用因素。改革開(kāi)放以前,湯遜湖及周邊是重要的農(nóng)業(yè)生產(chǎn)區(qū),水產(chǎn)養(yǎng)殖是臨湯遜湖地區(qū)的主要經(jīng)濟(jì)發(fā)展方式，湖泊面積總體較穩(wěn)定。隨著經(jīng)濟(jì)社會(huì)的發(fā)展、交通條件的改善和消費(fèi)觀念的改變,湯遜湖地區(qū)開(kāi)始建設(shè)不同級(jí)別的開(kāi)發(fā)區(qū)，極大地改變了原有土地利用狀況，促使“向湖泊要土地”，導(dǎo)致湖泊面積不斷減少。特別是2010年以后，開(kāi)發(fā)區(qū)建設(shè)進(jìn)入快速發(fā)展期，土地利用類型大幅變化，與2009—2019年湯遜湖濕地面積加快減少的趨勢(shì)相符合。

(3) 工程活動(dòng)因素。湯遜湖地區(qū)自1990年開(kāi)啟了房地產(chǎn)建設(shè)的浪潮[19],各類住宅、寫字樓、別墅等主要沿湖岸和主干道分布。與此同時(shí)，湯遜湖大橋和江夏大道建成后，分割了內(nèi)湖與外湖水體,后續(xù)修建的大量水上道路設(shè)施使得湯遜湖水體進(jìn)一步被人為分割。造成的結(jié)果是，湯遜湖已很難找到原生的湖岸,大部分湖岸都被硬化,湯遜湖面積的自然恢復(fù)幾無(wú)可能，工程建設(shè)活動(dòng)人為固化了湯遜湖的邊界。

4 結(jié)論

(1) 應(yīng)用Landsat衛(wèi)星影像數(shù)據(jù)解析了五個(gè)不同年份(1980、1990、2000、2009、2019年)湯遜湖濕地的濕地類型和土地利用情況，發(fā)現(xiàn)湯遜湖主要分為人工濕地、湖泊濕地和沼澤濕地三種濕地類型，其周邊土地利用情況以居民地、濕地、植被覆蓋區(qū)和其他用地類型為主。

(2) 1980—2000年，湯遜湖人工濕地、湖泊濕地、沼澤濕地的面積分別為增加、減少、減少；2000—2019年，湯遜湖人工濕地、湖泊濕地、沼澤濕地的面積分別為減少、增加、減少；湯遜湖濕地面積整體為1980—2000年基本穩(wěn)定，2000—2019年加速減少。湯遜湖在洪水期和枯水期的濕地類型比較單一，人工濕地、湖泊濕地的面積在洪水期分別為減少和增加。

(3) 引起湯遜湖濕地格局變遷的因素主要為降雨和人口增長(zhǎng)、土地利用類型改變、工程活動(dòng)等人類活動(dòng)。

致謝：感謝湖北省地質(zhì)調(diào)查院李定遠(yuǎn)、劉成新等教授級(jí)高級(jí)工程師在論文撰寫過(guò)程中提出的寶貴建議！