东洞庭湖水域面积变化与三峡水库水域面积变化的关系

　　洞庭湖是中国第二大淡水湖，与长江干流直 接通联，能拦蓄流域上游来水，分蓄长江干流洪水，消减干流河段的洪峰流量，缓解长江中游地区洪涝灾害，减轻下游洪水压力，维系地区的洪水蓄泄和泥沙的冲淤平衡。由于受泥沙淤积、围湖造田等因素的影响，洞庭湖被分为南洞庭湖、西洞庭湖和东洞庭湖3部分。3个湖体一方面能够相互连通，另一方面由于受各自湖底高程、河流入水量及人类活动等的影响而各自具有不同的特征。不同湖体水面变化特征与整个洞庭湖水面面积变化是否一致值得关注。长江上游三峡工程的建设是近年对洞庭湖影响最为剧烈的人类活动，对洞庭湖水域面积变化产生了深刻的影响。洞庭湖水域面积变化及其与三峡工程相关性研究一直受到国内外学者的关注，袁敏等利用中分辨率成像光谱仪 (moderateresolution imaging spectrometer，MODIS)数据研究三峡工程对洞庭湖水面的影响，认为三峡工程建成后洞庭湖水面有减少的趋势。泥沙淤积、东方田鼠暴发等洞庭湖生态环境问题也因为洞庭湖水位的变化而发生改变。

　　遥感技术具有大尺度、可回溯等特点。基于遥感技术研究洞庭湖水面或水量的时空变化成为了近年的研究热点。龟山哲等、彭定志等和龚伟等均基于 MODIS 卫星对洞庭湖区水面变动情况进行了分析。Ding Xianwen 等利用ENVI–ASAR 影像分析了 2002—2009 年洞庭湖水面变化情况，认为洞庭湖枯水期水面较 20 世纪 90 年代有所增加。虽然关于整个洞庭湖水面变化的研究较多，但关于东洞庭湖研究的报道极少。在洞庭湖水域面积变化的已有研究中，通常以 3~5 年或者1 年为研究周期，这不易反映长时间序列的变化情况。考虑到东洞庭湖直接与长江干流相连通，也是江湖水量交换的直接场所，所以本研究中重点关注东洞庭湖，利用 MODIS 影像数据分析 2002—2012年偶数年份东洞庭湖水域面积的月、季节变化以及水域消落带变化，并与三峡水库水域面积变化进行相关性分析。

　　1、 东洞庭湖概况与研究方法

　　1.1 东洞庭湖概况

　　东洞庭湖是“国际湿地公约”收录的由中国政府指定的21个国际重要湿地自然保护区之一。东洞庭湖(N 28°59″~29°38″，E 112°43″~113°15″)位于洞庭湖东北部。洞庭湖水通过东洞庭湖的城陵矶注入长江。东洞庭湖水域面积变化极大，存在明显的洪水期与枯水期。每年4月开始，东洞庭湖水位上涨，7—8月达到最高峰，9月开始下降，进入平水期，12月到翌年3月为枯水期，水位达年内最低。

　　1.2 水域面积的提取方法

　　MODIS 遥感影像来源于从 NASA 的 EarthObserving System Data Gateway (EDG)网站免费下载的空间分辨率为 250m 的 Terra/MODIS MOD09Q1数据。数据格式为 MODIS V005 版本。东洞庭湖与三峡水库各相应年份遥感数据的时间分布见表 1。利用 EDG 网站提供的免费软件 MODIS ReprojectionTool–MRT 4.0，对以上遥感数据分别进行投影转换和研究区提取。

　　现有的利用卫星遥感数据进行水体信息提取的研究中，针对湖泊水体的提取方法主要有单波段阈值法、植被指数法、水体归一化指数法等。

　　这些方法的共同特点是利用 1 个阈值作为最后区分水体与其他地物的标准。对于少量遥感影像数据而言，这些方法可能会取得较好的效果，但对较长时间序列的影像数据而言，只利用 1 个阈值可能会造成某些影像数据的提取误差。本研究中利用 MODIS传感器的 250 m 分辨率影像数据，采用单波段阈值法，以多阈值灵活设置和目视解译为辅进行水体提取，这样能有效地剔除所提取水体中混淆的其他地物信息，提高水域提取准确率，能更真实地反映水体空间分布。在获得各研究年份东洞庭湖各月水域面积的基础上，为进一步明确东洞庭湖水域面积的季节动态，取各季节不同月水域面积的均值作为东洞庭湖该季节的水域面积。分别对东洞庭湖各研究年份不同季节水域边界进行叠加分析，可获得东洞庭湖水域消落带的分布情况。利用 ENVI 4.8 与ArcGIS 10.0 软件进行水体解译与空间分析。利用SPSS 18.0 的 Pearson 进行相关性分析。

　　2、 结果与分析

　　2.1 东洞庭湖水域面积的月变化和季节变化2002、2004、2006、2008、2010、2012 年东洞庭湖各月水域面积的最大值分别为 1 633.93(8 月)、1 258.93(8 月)、1 469.82 (7 月)、1 212.46(11 月)、1 253.57(9 月)、1 401.11 km2(8 月)，各年水域面积最大值均超过了 1 200 km2。除 2008 年外，东洞庭湖各年水域面积的最大值均出现在 7、8、9 月，水域面积的变化趋势与洞庭湖区域降水量的变化趋势相同。2002、2006、2010、2012 年东洞庭湖水域面积的最小值均出现在 1 月，水域面积分别为294.64、282.25、267.86、189.47 km2，而 2004 年和2008 年东洞庭湖水域面积的最小值均出现在 2 月，分别为 300.00、261.72 km2。东洞庭湖水域面积最大值与最小值的差值在 2002 年最大，达 1 339.29km2，在 2008 年最小，仅为 950.74 km2(图 1)。

　　整体而言，2002—2012 年偶数年份，东洞庭湖全年的水域面积最小值出现在 1、2 月，随后逐渐增加，大致在 7、8、9 月达到最大值，随后降低。

　　春季东洞庭湖水域面积约为 599.17 km2，研究期内春季东洞庭湖的水域面积 2002 年的最大，为763.98 km2，2008 年的最小，仅为 343.40 km2，研究期内夏季东洞庭湖水域面积相差不大，平均面积为 1 120.71 km2，最大面积为 1 162.90 km2，最小面积为 1 030.19 km2；研究期内秋季东洞庭湖的水域面积 2006 年的最小，仅为 587.10 km2，其他年份约为 900 km2；研究期内冬季东洞庭湖的水域面积相差不大，平均面积约为 326.92 km2，是一年中水域面积最小的季节。由图 2 可见，东洞庭湖夏季的水域面积最大，春季的次之，冬季的最小。

　　2.2 东洞庭湖水域消落带的年内变化图3为2002年至2012年偶数年份东洞庭湖水域消落带的空间分布。2002、2006、2008、2010年春、秋、冬季的水域空间分布形状较为相似，夏季的水域面积较大，且增加的区域主要集中在南部区域；2004 年，东洞庭湖夏、秋季水域的空间分布形状相似，春、冬季的水域空间分布形状相似；2012 年，东洞庭湖夏季的水域面积最大，春季水域空间分布范围较小。整体而言，东洞庭湖消落带的变化呈相似的季节变化和空间分布：一方面，东洞庭湖春、秋、冬季的水域空间分布形状较为相似；另一方面，东洞庭湖水域消落带主要分布在南部湖区沿岸。

　　【图3略】

　　2.3 东洞庭湖与三峡水库水域面积变化的对比分析图 9～14 为 2002—2012 年偶数年份各月东洞庭湖和三峡水库水域面积的变化情况。2002 年，三峡水域年均面积为(17.98±4.06) km2，8 月达到最大，为 24.63 km2，较东洞庭湖水域最大面积的出现时间提前了 1 个月；4 月的面积最小，为 12.38 km2。三峡水库水域面积全年存在 2 个峰值(2 月和 8 月)和 2个谷值(4 月和 9 月)。东洞庭湖水域面积与三峡水库水域面积的相关性分析结果表明：2002 年东洞庭湖水域面积和三峡水库水域面积存在负相关关系，但相关性不显著(r = –0.035, P>0.05，图 9)。2004 年，三峡水库水域面积全年波动不大，年均水域面积为(26.59±4.27) km2，6 月最大，2 月最小。2004 年，东洞庭湖水域面积变化和三峡水库水域面积的变化存在一定的相关性，相关系数为 0.726(P<0.05)，但东洞庭湖水域面积的变化幅度比三峡水库水域面积的变化幅度大(图 10)。2006 年，东洞庭湖水域面积与三峡水库水域面积也呈正相关关系，但相关性不显著(r=0.311，P>0.05，图 11)。2008 年，东洞庭湖水域面积变化与三峡水库水域面积变化呈负相关关系，但相关性不显著(r= –0.547, P>0.05，图12)。2010 年和 2012 年三峡水库水域面积在 3 月达到最大，在 6 月最小。东洞庭湖水域面积与三峡水库水域面积的负相关关系均不显著(P>0.05)。1—3月东洞庭湖和三峡水库水域面积变化趋势相同，存在正相关关系；4—5 月二者的变化趋势相反，存在负相关关系； 6—8 月的变化趋势再次趋同(图 13、14)。

　　虽然除 2004 年外，研究期间东洞庭湖水域面积与三峡水库水域面积的变化均不存在显著相关性，但是二者之间的关系日趋明显，即在三峡工程放水、蓄水等的调控下，在全年中的某些时间段(洪水期与枯水期)，二者之间可能呈现明显的相关性。

　　由于正相关与负相关的相互掩盖，导致以年为单位进行的相关性分析结果不明显。此外，除 2004、2006年东洞庭湖水域面积与三峡水库水域面积呈现正相关关系外，其他年份二者更多地呈负相关关系。

　　东洞庭湖水域面积与三峡水库水域面积的正相关关系反映了区域降雨的一致性，而二者负相关关系则可能反映三峡调控对东洞庭湖的影响。

　　3、 结论与讨论

　　本研究中发现东洞庭湖水域面积每年从 1 月开始逐渐增加，在 7、8、9 月达到峰值，随后开始降低，到 12 月近似达到全年最低。该变化趋势与整个洞庭湖水面的变化趋势一致。洞庭湖是典型的过水性湖泊，受流域降水和河流入水量的季节变动影响，常呈现典型的“洪水一片、枯水一线”景观。

　　在当前的研究中，洞庭湖被划分为洪水期(4 月至 10月)和枯水期(11 月至翌年 3 月)。这种区分具有较好的理论意义，但是也可能忽视了农业生产的需求。

　　本研究中关于东洞庭湖水面的月变化和季节变化分析结果，可为农业生产和生物保护提供参考。

　　洞庭湖水域面积的空间分布研究结果对洞庭湖区的防灾、减灾和农业生产具有指导意义。JixiGao 等通过分析洞庭湖水面面积的月变化，建立了洞庭湖洪水灾害指数，Yamei Wang 等利用半量化模型和模糊层次分析方法，绘制了洞庭湖洪水风险分布图，但当前还没有关于洞庭湖季节消落带的相关研究。淹–干交替是洞庭湖湿地大部分区域的典型特征，季节消落带区域对维持洞庭湖湿地生物多样性具有重要作用。本研究中发现东洞庭湖区域的季节消落带主要集中在南部湖区沿岸。该区域是东洞庭湖国家级自然保护区的核心区，具有较高的生物多样性，也是保护的重点区域。虽然东洞庭湖成立了国家级自然保护区，但当前仍然存在杨树、芦苇种植和龙虾养殖等农业生产活动。由于水面状况不稳定，该区域其实并不适宜进行龙虾养殖等活动。三峡工程、退田还湖等重大政策性工程是近年影响洞庭湖水面面积变化的重要因素，较气候变化、泥沙淤积等的影响更加剧烈，更加直接。

　　利用 BP 神经网络模拟，发现三峡工程蓄水后对洞庭湖水位的影响剧烈，其中，2009 年蓄水使城陵矶水位下降了 2.11 m。据邹邵林等的预测，三峡工程建成后洞庭湖洲滩露出时间将增加。本研究中通过分析 2002—2012 年偶数年份各月东洞庭湖和三峡水库的水域面积变化，发现三峡工程蓄水前三峡水库水域面积和东洞庭湖水域面积的相关性发生了巨大的改变，发生这种改变的原因可能是 2003年三峡大坝开始蓄水对下游洞庭湖水域面积及水位产生了重要影响，起到了调节径流的作用。

　　冬、夏季三峡水库与东洞庭湖的水域面积变化趋势相同，春、秋季水域面积的变化趋势则相反，因此，冬、夏季不论是三峡水库还是东洞庭湖均主要受自然降水的影响，而三峡工程的调控作用主要体现在春、秋季。

　　本研究结果表明：2002—2012 年偶数年份东洞庭湖水域面积的最低值出现在 1、2 月，随后增加，到 7、8、9 月达到最大值，随后又降低(个别年份有异常)；从季节变化来看，通常表现为夏季的水域面积最大，春季次之，冬季最小；东洞庭湖的水域消落带主要分布在南部湖区沿岸；三峡工程蓄水前后，三峡水库和东洞庭湖水域面积的相关性发生了改变：三峡水库蓄洪之前，其水域面积与东洞庭湖的水域面积并未出现同步变化的现象，而蓄洪之后，2个水域面积的变化呈现出明显相关性。

　　参考文献：
　　[1] 田伟国，彭嘉栋，沈军，等．基于 MODIS 影像序列的三峡截流前后洞庭湖面积变化序列分析[J]．安徽农业科学，2012(16)：9141–9145．