Water pattern changes of Tongxiang City based on multi-source remote sensing data

浙江省嘉兴市桐乡市位于浙北平原地区，杭嘉湖平原腹地，30°28′~30°48′N，120°17′~120°40′E，行政区下辖3个街道和8个乡镇。为了分析研究区内水系演变在城市化程度不同区域内的演变特征，定义市政府所在的梧桐街道为市区，与市区直接接壤的6个街道与乡镇为近郊，远离市区不直接与之相邻的4个乡镇为远郊。桐乡市境内河流属长江流域太湖运河水系，京杭大运河斜贯全境，全市有骨干河道47条，组成纵横交错的水利引排灌溉渠道和水上运输网络。桐乡市行政区划分与城市主要水系分布见图1。

Figure 1.  Overview of Tongxiang City and the distribution of major water systems

考虑到数据的可获得性以及桐乡市不同阶段发展政策的差异，选取了1998、2005、2012、2019年4个时间节点。通过新闻报道和对当地居民的调研，发现桐乡市居住区的水环境在2016年前后发生了明显变化，因此，将2016年作为2012—2019年的中间节点，并将其纳入到分析过程中。为了准确了解改革开放前的水系格局，选择经过解密的美国锁眼侦查卫星KH-9在1975年的影像作为原始对照。选取的影像数据成像时间均为当年的冬季，云量小于10%，具体参数与来源信息说明详见表1。

以不同时期土地利用调查数据为基础，通过对遥感影像的目视解译来提取不同时期的水系数量和分布数据。1975年水系面状矢量数据以可获取的最早土地利用调查数据(1988 年)为底图，通过遥感影像目视解译获取。另外5期的水系面状矢量数据以对应年份的研究区土地利用变更调查数据(1998、2005、2012、2016年)、第3次全国土地调查数据(2019年)为基础目视解译获取。在 ArcGIS 中将解译得到的面状河流水系转换为栅格，使用软件中针对栅格的中心线自动捕捉工具获取线性河流矢量数据，对自动提取的结果进行手动修正，提高数据的准确度。河流面积和河流长度分别指河流水面面积和河流中心线长度，使用河流水系面状矢量与线性矢量计算获得。

池塘的实际功能一般分为景观池塘、坑塘水面和养殖坑塘。由于土地利用调查数据中已有坑塘水面与养殖坑塘2类数据，因此本研究在分类时直接参考了这2类池塘。目视解译获取的新增池塘，通过比对同期确定功能的池塘的形态、结构及周边地类分布情况来判断其是否为养殖坑塘。同时，本研究将公园绿地和旅游景区内部的池塘水面归类为景观池塘。

各时期的水系数据在进行目视解译时均由当地政府提供的土地利用调查数据为基础，真实性高，精度可满足后续的研究要求。本研究使用的所有土地利用调查数据与村镇级行政区矢量数据均由桐乡市自然资源和规划局提供。

考虑到水系的景观生态功能与社会属性，以河道的平均宽度作为划分依据，对河流进行了分级^{[12, 20]}。其中，将河道平均宽度大于40 m的河流定为Ⅰ级河流， 20~40 m为Ⅱ级，小于20 m定为Ⅲ级。Ⅰ级与Ⅱ级河流统称为主干河流，构成了研究区水系的骨架。

针对平原水系的特征，参考国内外相关研究常用的水系特征参数与表征指标^[28−30]，选用统一的指标识别并定量评价不同时期桐乡市河网水系格局。选取的河网结构指标包括描述河网数量特征的河网密度(D_r)、河网水面率(W_p)；描述河网形态特征的支流发育系数(K)、分形维度(D₀)；描述河网连通性的水系环度(α)、实际结合度(γ)。由于池塘在空间中分散且相互之间不连通，为了分析池塘分布结构和形态的变化特征，本研究采用景观格局指数法^[31] 对池塘斑块进行了定量计算。选用的景观指数包括景观形状指数(LSI)、几何最近邻距离加权平均指数(ENN_AW)和聚合度指数(AI)。

按照河流功能与宽度对桐乡市的河流进行分级，统计各级河流在历史各时期的长度和面积变化(表2)。桐乡市的主干河流在1998年相较1975年增加了8.29%的长度和30.95%的面积，其中Ⅰ级河流的长度增加了102.40%。尽管Ⅰ级河流数量最少，但在1998—2016年，桐乡市Ⅰ级河流长度持续增加，仅在2016—2019年略微减小。在40 a间，Ⅱ级河流的长度和面积不断波动，总体略有减小，其中长度和面积的降幅分别为6.25%和10.62%。支流是桐乡市河流水系变化的主要因素，其长度在所有河流中占比最大，各时期变化幅度也较为明显。仅有1975和2005年，支流面积大于主干河流，其余年份支流面积均小于主干河流。在2012和2019年，支流长度和面积接近，但整个过程中有明显减小和重新增加的趋势。

采用水系指数法对桐乡市各时期的河网数量、结构、连通性特征的各项指标进行定量计算。根据表3可知：桐乡市的各项河网结构指标在2019年较1975年均有所增加，其中河网密度增加了12.64%，河流水面率上升了12.12%，支流发育系数增加了1.30%，分形维度增加了3.97%，水系环度和实际结合度则基本保持一致。桐乡市的河网密度和河流水面率在1975—1998年均呈增长趋势，但同时支流发育系数和分形维度的变化并不明显。1998—2005年，桐乡市的河网密度继续增加，增长率达9.58%，但河流水面率基本保持稳定，同时支流发育系数增幅达41.60%，河网连通性也恢复到了30 a前的水平。

2005—2012年，桐乡市的河网密度、支流发育系数、分形维度、水系环度、实际结合度都急剧下降。城区河道重新规划布局，导致市内主干河流持续增长。由于主干河道宽度较大，尽管河网密度急剧下降，但桐乡市的河流水面率在2012年并未下降，这说明河网主干化和单一化的趋势更加明显。

2012—2016年，桐乡市的河网密度和水面率分别下降了16.61%和12.77%，而支流发育系数则减小了12.92%。与此同时，河流的分形维度、水系环度和实际结合度也显著增加。2016—2019年，桐乡市的河流密度和河流水面率分别增加了21.99%和18.50%，其中河流水面率到2019年已经达到历史最高值(5.83%)。在这一阶段，支流发育系数和分形维度也恢复增加，整体呈现出向好的发展态势。

定量计算桐乡市城市各区域在各个阶段的河网数量特征变化，得到了城市不同区域河网密度与河流水面率的变化情况，具体结果见表4。支流发育系数是评估人为工程对河流改造的主要指标，以城市区域作为划分标准，探究各时期该指标在空间上的数值大小差异，具体结果见表5。1975—1998年，桐乡市的河网密度与河流水面率变化情况相同。近郊区域的支流发育系数受工程影响变化显著，从4.86减小为3.66，减少了约25%。1998—2005年，各区域河网密度均发生增加，涨幅近郊＞远郊＞市区。近郊河网发育系数快速回升，达到30 a前的水平。2005年后，桐乡市开始进行城市改建与扩张，导致所有区域河网密度显著减小。但由于防涝的需要，桐乡市的主干河流在城市化建设中进一步提升宽度，各区域河流水面率不降反升。这一阶段也是河网发育系数衰减最为剧烈的时期，衰减幅度市区＞近郊＞远郊。这与城市扩张的方向相吻合，因为建设用地的扩张往往导致低等级的河流被填埋或转入地下。2012—2019年，桐乡市开始对市内水环境进行修复，各区域的支流发育系数降速得到了遏制，远郊区域甚至有一定的回升。40多年的发展过程中，市区、近郊和远郊的支流发育系数分别减小了62.68%、34.98%和21.24%。1975年，市区的支流发育系数约为远郊的2倍，并与近郊相当，但到2019年时，市区的数值仅为近郊的一半。

2005—2012年，受到快速城镇化的影响，桐乡市的河流数量明显减少。通过计算各区域单位面积内河流水面被转化为其他地类的面积占比，可以发现这种趋势。由图2可见：在市区内，建设用地的增长是导致河流转出的最主要原因，其次是园地，这两者合计占到了转出河流面积的70%。在近郊和远郊，耕地是河流转出最多的地类，占总值的一半。远郊区域是桐乡市主要的水产养殖片区，当地居民利用天然河流水面围建鱼塘进行养殖，在建设过程中有8.68%的转出河流变成了池塘。

Figure 2.  Proportion of land types transferred out of river water surface in Tongxiang City between 2005 and 2012

由表6分析可知：桐乡市池塘水系面积在40多年的发展过程中迅速增长。截至2019年，池塘水系的总面积已经达到了1975年的9倍。1975年，桐乡市仅有少量的坑塘水面，到1998年时池塘水系的总面积变化不显著，但水面结构组成发生了较大的变化，养殖坑塘的面积占到了总面积的68.60%。1998—2005年，池塘水面总量增加了约300%，同时水面结构组成基本保持不变。2005—2012年，坑塘水面大量消失，面积衰减了64.06%，养殖坑塘面积的增加并不显著，桐乡市池塘水面总量整体减少，与河网的变化特征基本吻合。2012—2019年，坑塘水面面积增长了3倍，景观池塘的增长速度达到了历史高值，池塘面积总量完成了翻倍。

由表7可知：在1975年，桐乡市池塘水面的聚合度指数最低，而几何最近邻距离加权平均指数达264.04 m，这表明在自然状况下桐乡市的天然坑塘处于无序排布的状态，景观的邻近度、离散度和聚集度均处于较差的水平。在接下来的20余年间，池塘景观形状指数从49.29下降至27.12，这说明在天然坑塘转化为养殖坑塘的过程中，池塘的形状通过人工改造变得更加规整。到2005年，由于湖泊、鱼塘与稻田的紧密排布，池塘斑块的几何最近邻距离加权平均指数达到了历史最低，为66.92 m，同时景观形状指数达到最高值，为91.70。随着末级支流与池塘的大量填埋，至2012年桐乡市池塘的分布重新变得离散。由于形状规整的人工鱼塘在池塘占比中的提高，景观形状指数重新降至低值。到2019年，桐乡市池塘景观通过各项工程得到了较好的治理，斑块间连片性明显提高，呈现聚集分布的态势，同时景观形状指数也恢复至40 a前的水平。

根据表8可知：1975—1998年，桐乡市区和近郊池塘水面率快速下降，远郊池塘水面率上升。1998—2005年，所有乡镇池塘水面率均显著上升，其中市区＞近郊＞远郊。2005—2012年，随着城市的高速扩张，桐乡市区和近郊的池塘开始快速衰减，其中近郊的衰减幅度大于市区，整体衰减幅度达47.15%。远郊地区由于西部池塘聚集区的向外扩散，池塘水面率依旧保持上升，但相对上一阶段升速减缓。2012—2019年，各区域池塘水面率恢复上升，市区增加的多为景观池塘，近郊多为坑塘水面，远郊以养殖坑塘为主。由于各区域池塘基数差距明显，因此上升率依然是市区＞近郊＞远郊。

采用与河流相同的研究方式，对2005—2012年池塘水面衰减的空间规律进行了分析，计算了各区域单位面积内池塘水面转化为其他地类的面积占比。由图3可知：消失的池塘主要被转变为耕地和建设用地，而市区的转化比例则明显高于其他地区。近郊和远郊的池塘转化结构相似，各土地类型的比例差异不大。与河流不同，池塘消失后转化为园地的比例极低，在近郊和远郊却有一定数量的池塘被人工开挖成为新的河流水面。

Figure 3.  Proportion of land types transferred out of pond water surface in Tongxiang City between 2005 and 2012

桐乡市的水系变化模式具有明显的地方特色，受区域政策的影响，在不同时间段呈现不同的变化特征，与同类型平原水网城市存在异同。本研究分析总结了桐乡市的地方发展政策，从城镇化、农业结构调整和生态修复与治理三方面解释桐乡市水系演变的原因。

改革开放后，桐乡市为满足民生防涝的需求，持续开展城镇化与水利建设。自1978年起，嘉兴市开展了长山河南排工程，整个工程历时20 a，主要对桐乡市境内自然河流的河道进行了拓宽和功能改造，但并未对河网整体分布形态结构做较大改动^[32]。在此过程中，桐乡市原有的自然河流功能划分更加明确，主干河道的长度与宽度大规模增加，支流面积减小，不同等级河流形态差异增大，与嘉兴市整体变化保持一致^[25]。2005—2012年，桐乡市城镇化率从21.84%提升至39.34%，进入了快速城镇化时期^[33]，开展了大规模的重新规划以及结构改造，但同时城区的扩张也导致河道的淤积和建筑用地侵占水面。桐乡市快速城镇化时期的水系变化特征与苏州^[34]、宁波^[35]等地的变化相似，支流快速衰减，河网密度大幅减小，河网连通性下降，市区主干河道因景观和防涝需要持续拓宽，河网呈现明显的主干化和单一化趋势，大量功能不明确且经济效益较低的坑塘水面消失。

桐乡市的农业产业结构调整是造成市区外支流与池塘变化的最主要因素。从20世纪 80 年代起在远郊的农业区中开始发展经济效益高的水产养殖业，大量坑塘水面被改造成养殖坑塘，1998 年养殖坑塘已在池塘中占比最大。在2001年后因“改水”政策，农村居民开始大规模挖掘水渠，将主干河流的地表水引入村庄和农田，由于挖掘的水渠仅为引水而宽度较窄，造成2005年前桐乡市河网密度、支流发育系数增长，但河流水面面积变化不明显，该变化特征与同时期农村地区生产方式相似的苏州^[26]和南京^[36]相同。而无锡^[20]、鄞东南^[35]、苏北里下河地区^[24]等地的农业生产方式在发展过程中未发生明显变化，故境内河网依旧呈现随时间推移持续萎缩特征。同时，在桐乡城镇化发展稳定后，随养殖塘标准化建设及老塘改造工程的开展，远郊的养殖坑塘得到整合，发展重心从大麻镇向北转移到洲泉镇中，使桐乡市后期池塘面积大幅提升且在空间上大规模聚集连片。

桐乡市的水系生态治理与修复在2012年后逐步开展。该市于2014年启动了“五水共治”工程，在初期部分末端河道被填埋，大量工程设施与材料占据了水面区域，导致河网密度、河流水面率、支流发育系数在2012—2016年有一定幅度的减小。随着工程的进展，许多被填埋的天然湖泊得到了复原，同时在常住人口聚集区与景区内增加景观池塘。各项治理措施使桐乡市支流快速恢复，河网连通度增加，各类池塘快速增长，养殖坑塘结构规整且聚集排布，水环境得到明显优化，至2019年时各项指标相较1975年均有提升，说明平原河网地区河流指标数值与城镇化水平并不是简单的负相关关系^[19]，在定量研究时还需考虑城镇化发展后期生态修复等政策的影响。

Ⅰ级河流在40余年间持续增长，低等级的支流在1998—2005年显著增长，在2005—2012年的快速城镇化时期与2012—2016年的“五水共治”初期大幅衰减达到历史低值，到2019年时支流长度与面积均恢复至较高水平。

桐乡市的河网密度在2005年前保持增加，在桐乡市进入快速城市化后持续快速衰减，于2016年后重新回增。河流水面率在2012年前持续上升，但升速显著放缓，于2012—2016年数值快速下落达到历史低值，但到2019年时数值重新回到最高点。支流发育系数和分形维度在40余年间同样经历了先增加再减小最后回增的过程，转折时间节点分别是2005和2016年。水系环度和实际结合度在整个过程中不断波动，但在其他指标均最低的2016年达到峰值。

在城市不同区域河流变化差异显著。河网密度在近郊与远郊的变化远大于市区，河流水面率在市区的变化最为显著。在2005年前河流水面率变化远郊＞近郊，但在2005年后远郊＜近郊。支流发育系数的变化市区与近郊明显高于远郊。河流消亡转入的地类在不同区域存在明显差异，人类活动的需要是造成空间变化差异的主要因素。

桐乡市池塘结构组成在发展过程中有较大的改变，养殖坑塘在40余年间从0高速增加至20.93 km²，在当前池塘水面中的占比达77%。景观池塘在2012年前变化不明显，但随后增加迅速。坑塘水面在1975年时几乎构成了桐乡市的全部池塘，但在发展过程中波动明显，与河网支流的变化规律基本相同。2005年时池塘整体形状最复杂且分布最为聚集，快速城市化进程使池塘整体结构变得更加规范但整体变得离散，经过治理后至2019年池塘间的整体连片性重新增强。

在整个发展过程中远郊池塘总量不断增加，近郊与市区变化趋势保持一致，且市区变化幅度仅在2005—2012年的快速城镇化阶段高于近郊，养殖渔业的规模化与标准化是造成池塘水面大量增加与聚集的主要因素。