黄河三角洲地区沿海养殖池塘空间演变阶段性划分及特征分析

摘要基于1984—2020年遥感影像数据，获取连续长时间序列的养殖池塘空间分布数据，并采用时空热点分析和PLR_FP方法对黄河三角洲地区养殖池塘空间演变的阶段性进行划分。结果表明，黄河三角洲地区养殖池塘面积呈现波动增长态势，年均增加29.3 km²，时空增长模式多样化，且具有明显的空间分异性；增长热点分布在滨州沿海、东营河口区西部、垦利黄河口南部，成为目前黄河三角洲地区池塘养殖重点分布区，而莱州湾沿岸池塘养殖规模逐渐萎缩。养殖池塘的空间演变过程可划分为起步期（1984—1992年）、平稳期（1993—1998年）、扩张期（1999—2013年）和缩减期（2014—2020年）4个阶段，不同阶段养殖池塘时空增长模式和养殖业发展特征差异明显。

关键词养殖池塘；空间演变；阶段划分；变化特征；黄河三角洲

中图分类号S967.4文献标识码A

文章编号0517-6611（2023）11-0050-06

doi：10.3969/j.issn.0517-6611.2023.11.012开放科学（资源服务）标识码（OSID）：

Stage Division and Characteristics Analysis of the Spatial Evolution of Coastal Aquaculture Ponds in Yellow River Delta

WANG Zi-yang WU Xiao-qing DOU Xiao4 et al（1. Yantai Institute of Coastal Zone Research， Chinese Academy of Sciences， Yantai，Shandong 264003;2. University of Chinese Academy of Sciences， Beijing 100049;3. Key Laboratory of Coastal Environmental Processes and Ecological Remediation， Chinese Academy of Sciences， Yantai，Shandong 264003;4. College of Public Administration， Shandong Agricultural University，Tai’an，Shandong 271000）

AbstractBased on the remote sensing image data from 1984 to 2020， this study obtained the spatial distribution data of aquaculture ponds in a long-term continuous sequence，and the spatio-temporal hotspot analysis and PLR_FP method were used to divide the stages of the spatial evolution of aquaculture ponds in Yellow River Delta.The results showed that the area of aquaculture ponds in the Yellow River Delta showed a fluctuating growth trend， with an average annual increase of 29.3 km2. The spatiotemporal growth patterns were diverse， with obvious spatial differentiation. Growth hotspots were distributed along the coast of Binzhou， the west of Dongying Estuary and the south of Kenli Yellow River Estuary， which had become the key distribution areas for pond aquaculture in the Yellow River Delta， while the scale of pond aquaculture along the Laizhou Bay was gradually shrinking. The spatial evolution process of aquaculture ponds could be divided into four stages： the initial stage （1984-1992）， the stable stage （1993-1998）， the expansion stage （1999-2013） and the shrinking stage （2014-2020）， and there were significant differences in the spatiotemporal growth patterns and aquaculture characteristics of aquaculture ponds at different stages.

Key wordsAquaculture ponds;Spatial evolution;Staged division;Change characteristics;Yellow River Delta

池塘养殖是海水养殖的一种重要模式，作为“优质粮食”的生产场所，其生态灾害风险较低，经济收益较高，对我国“蓝色粮仓”建设具有重要意义［1］。现代黄河三角洲及莱州湾沿岸为典型的三角洲海岸和淤泥质海岸，滩涂资源丰富、池塘养殖历史悠久、养殖规模大，是山东省重要的养殖基地［2］。近几十年来，随着社会经济发展水平的提高和养殖技术的进步，黄河三角洲地区池塘养殖不断向渤海扩张，占用了大量自然滩涂湿地，造成滨海湿地萎缩、生物多样性下降、土壤次生盐渍化、碳封存环境破坏等问题［3-7］。与此同时，养殖池塘也是重要的人工湿地，具有补充地下水和调节气候等生态功能［8-9］，科学合理地进行池塘养殖还能够成为践行碳中和战略的重要负排放途径［10-11］。在推动黄河流域生态保护和高质量发展战略背景下，实现池塘养殖功能优化布局，成为化解滨海湿地生态保护与高效生态化养殖之间矛盾的重要管理手段之一。而研究沿海池塘养殖时空演变特征作为其优化布局的一项重要基础性工作，有助于揭示该地区海水养殖业发展历程，获取地区池塘养殖的独特性，为黄河三角洲地区滩涂湿地池塘养殖空间优化布局和生态保护修复工程规划建设提供数据和理论支持。

当前学者多从土地利用和景观格局变化的角度揭示黄河三角洲地区养殖池塘的时空演变特征及其对滨海湿地的生态影响［12-15］，但是由于其时空分布数据不具备连续性，或仅使用经济统计数据，造成对养殖池塘时空演变特征和变化趋势的把握不够准确，不能对海水养殖业发展过程作出有效判断。该研究基于遥感、相关专项调查资料等数据，重点揭示黄河三角洲地区1984—2020年池塘海水养殖时空演变过程，划分其空间演变阶段，探讨不同阶段空间增长特征及其驱动因素，以期为池塘养殖可持续管理与开发和环境保护提供技术支持。

1资料与方法

1.1研究区概况研究区位于山东半岛北部，西起鲁冀交界的漳卫新河河口，东至烟台、潍坊交界的胶莱河口，地理位置为117°76′～119°61′E、36°92′～38°28′N，包括黄河三角洲、莱州湾自海岸线向陆延伸15 km的海岸带地区，行政区涉及滨州、东营、潍坊3个地级市。此区域为粉砂淤泥质平原海岸，地势低平，滩涂范围宽广，有黄河、小清河、潍河、胶莱河等十余条淡水河流入海，近岸海域营养盐、有机质多，饵料充盈，是山东省重要的海水养殖基地。

1.2数据源与预处理遥感数据源主要为Landsat系列卫星，时间跨度为1984—2020年，通过地理空间数据云（http：//www.gscloud.cn/）平台下载获取，共计74景影像。在ENVI中进行预处理，包括辐射校正、大气校正、全色波段和多光谱的融合等。辅助数据包括地形图、历史年份的高清影像和海岸带专项调查数据等。海水养殖相关统计数据来源于省、市统计局历年统计年鉴，部分年份缺失数据采用移动平均法进行补充。

1.3研究方法

1.3.1遥感信息提取。养殖池塘的提取首先借助水体指数并参照明显的解译标志［16］识别其潜在区域；然后，结合地物遥感解译判别标志、专家经验知识、海岸带土地利用调查资料以及实地调研，采用人工目视解译方法对养殖池塘和盐田、水库坑塘进行区分，处理获得1984—2020年黄河三角洲地区连续37年的养殖池塘分布数据。利用同时期高清谷歌影像、高分一号等高分辨率遥感数据及研究团队近几年实地调查点位进行解译精度验证，结果显示37年养殖池塘总体解译精度均在79%以上，满足研究需要。

1.3.2阶段性划分方法。该研究基于养殖池塘空间分布连续多年的遥感观测数据，采用基于特征点的分段线性表示法（piecewise linear representation based feature points， PLR_FP）提取时间序列数据分段特征点［17］，进行池塘养殖空间扩展阶段性划分。PLR_FP方法可以概化时间序列数据的主要状态特征，进行有效的数据压缩和局部噪声过滤。具体做法：首先，识别时间序列数据的局部极值点；然后，采用阈值R限定前后极值点数据的垂直差比，将符合要求的点识别为不同尺度下的分段特征点［17-18］；最后，借助双样本Wilcoxon秩和检验将不具有显著性的特征点进行排除［19］，依据最终分段线性表示的拟合系数与常数0之间的关系判断每阶段的具体趋势。

其中，秩和检验用于检验2个独立的养殖序列样本是否来自相同的总体。从2个独立总体抽取样本并将所有观测值由大到小排列后，每一个观测值在混合排列中所占位置的名次称为它的秩，当其中某样本秩和（statistic）较小，混合样本为一个随机混合样本，秩和检验中P值小于0.05，可以认为该分段点前后时间序列数据具有明显差异，可对前后2个时段内养殖池塘演变特征相同这一假设提出质疑。

1.3.3空间演变模式判别方法。对于连续性面状空间分布数据，采用时空热点分析方法量化分析近37年养殖池塘的空间格局类型及其出现频率和时空位置［20］，揭示池塘养殖空间演变分异规律。

该方法首先基于ArcGIS热点分析方法计算各时间截面Getis-Ord G*统计量并进行显著性检验；空间权重矩阵采用邻接标准中的“Queen”，即面要素有公共边或结点即为相邻，采用“CONTIGUITY_EDGES_CORNERS”工具计算，邻域步长设置为1年；然后，分析各区域截面Getis-Ord G*统计量的Mann-Kendall趋势；其中，区域截面的构造采用覆盖研究区范围的1 km格网；最后，结合养殖池塘面积遥感观测数据，按照模式分类规则（表1）对比归类为不同时空模式［20］。

2结果与分析

2.1养殖池塘时空演变总体特征从图1可以看出，1984—2020年黄河三角洲地区沿海养殖池塘规模呈现波动增长的态势，从期初的66.7 km²增长至1 150.7 km²，增长约16.3倍，年均增速为29.3 km²。其中，1986—1989、2000—2005、2010—2013年养殖池塘面积增速较快，并于2013年达到最高峰（1 359.9 km²）；2016—2018年养殖池塘规模出现缩减，3年共减少231.1 km²；而1990—1999年养殖池塘面积保持在相对稳定的水平，可见养殖池塘规模长时间序列年度增长呈现明显的阶段性变化特点。

黄河三角洲地区沿海养殖池塘空间演变也呈现出明显的时空分异特征，不同时间段养殖池塘增长的冷热点发生迁移变化，时空演变模式呈现多样化（图2）。研究期初，滨州的无棣、沾化和潍坊市的昌邑、寒亭等县（市、区）的Gi*统计量较高，为区域内沿海池塘养殖的起步区；随后一段时间，各县（市、区）养殖池塘增长变化的Gi*统计量分布相对均匀，无明显的热点区，仅在寿光、河口、利津等县（市、区）保持较低水平。但是在2000年之后，养殖池塘区域增长的冷热点分布逐渐变得两级分明，沾化、无棣的Gi*统计量骤增，出现了新一轮的池塘养殖高潮，并带动与其毗邻的东营河口区池塘养殖的扩张，成为“逐渐增强的热点”，也是目前池塘养殖的重点分布区；而莱州湾西南部沿岸县（市、区）Gi*统计量明显降低，归为“逐渐减弱的热点”，因盐化工、临港产业发展和城镇化建设，该区域养殖池塘规模大面积萎缩。养殖池塘空间增长的冷点分布较为分散，其中“逐渐增强的冷点”和“分散的冷点”位于黄河三角洲国家级自然保护区，说明自然保护地建设管理对养殖池塘空间扩展具有明显的约束作用。