《Ecologies》：Effects of Riparian Land Use and Land Cover on Water Quality Along the Kansas River: Seasonal and Spatial Dynamics
编辑推荐：
研究人员结合季节t检验、单因素方差分析和冗余分析（RDA），在三个河岸缓冲尺度（500 m、1000 m和2000 m）下，研究了堪萨斯河干流四个美国地质调查局（USGS）站点（2019–2026年）的流量、溶解氧（DO）、温度、浊度和pH值。DO和温度表现出
研究人员结合季节t检验、单因素方差分析和冗余分析（RDA），在三个河岸缓冲尺度（500 m、1000 m和2000 m）下，研究了堪萨斯河干流四个美国地质调查局（USGS）站点（2019–2026年）的流量、溶解氧（DO）、温度、浊度和pH值。DO和温度表现出最强的季节对比：旱季DO高出3.1–3.7 mg/L，湿季温度高出13–15 °C，这种耦合对水生栖息地适宜性至关重要。浊度在湿季显著升高，与农业径流和泥沙迁移一致，而四个站点中有三个站点的流量无显著季节差异，反映了上游水库的调节作用。空间ANOVA仅在湿季DO（F = 3.12, p = 0.04）检测到站点间差异，该指标在下游城市化站点最低。RDA将农业覆盖与浊度、城市覆盖与湿季DO降低联系起来，尽管置换检验在n = 4个重复时未达到显著水平（p > 0.05）。季节性和河岸土地利用与土地覆盖（LULC）共同塑造了这条受调节河流的水质，而500 m缓冲带是土地覆盖评估中最具空间分辨力的尺度。
**论文解读文章**
**研究背景与问题**
淡水资源的可持续性受到土地利用/土地覆盖（LULC）变化、城市化和气候变率的加剧压力[1,2,3,4]。河岸带作为水陆交错区，对水质具有不成比例的影响，其植被能拦截地表径流、稳定河岸并调节水温[12,13,14]。然而，在受上游水库调节的农业-城市河流系统中，河岸LULC对水质的尺度依赖性和季节依赖性控制机制尚不明确[20,30]。现有研究多聚焦于未调节河流或采用单一缓冲宽度，缺乏整合多尺度缓冲带与季节分异的统一分析框架[2,14,31]。堪萨斯河作为大都会区的主要饮用水源（服务人口超过70万），其径流受Milford、Tuttle Creek和Perry等水库的调节，导致季节性流量变幅被压缩，这与典型未调节河流中LULC通过季节性径流传输污染物的机制不同[30]。因此，研究人员旨在（1）量化堪萨斯河干流四个USGS站点（2019–2026年）五个水质参数的季节与空间模式；（2）基于2023年国家土地覆盖数据库（NLCD）在三个缓冲尺度（500 m、1000 m、2000 m）表征河岸LULC组成；（3）通过冗余分析（RDA）探索多尺度LULC与水质变异的关系。该研究发表在《Ecologies》上，旨在为水库调节型河流的河岸缓冲带管理提供依据。
**主要技术方法**
研究人员采用以下关键方法：（1）从USGSWater Data for the Nation门户获取四个干流站点的连续传感器数据（2019.1–2026.6），并汇总为月均值；（2）利用2023年NLCD（30 m分辨率）提取三个缓冲尺度（500 m、1000 m、2000 m）的河岸LULC，将原始20类合并为7类（开发、荒地、森林、灌木、草本、农业、湿地）；（3）进行Lilliefors修正的Kolmogorov–Smirnov正态性检验后，采用独立样本t检验比较季节差异，单因素方差分析（ANOVA）及Tukey HSD事后检验比较站点间空间差异；（4）以四个站点为空间重复（n=4），对每个季节-缓冲尺度组合分别进行RDA，通过方差膨胀因子（VIF）筛选预测变量，并用置换检验评估模型显著性。样本来源为堪萨斯河干流四个USGS连续监测站（Wamego、Topeka、De Soto、NR Lake Quivira）。
**研究结果**
**3.1 河岸土地利用/土地覆盖（LULC）组成**
通过NLCD 2023提取三个缓冲尺度的LULC比例，发现站点间差异显著：Wamego站在500 m缓冲带中农业占主导（41.7%），De Soto站类似（58.6%），而Topeka站城市覆盖占72.3%，NR Lake Quivira站在所有尺度下城市覆盖最高（500 m处61.0%）。随着缓冲带拓宽至2000 m，模式总体一致但略有减弱：Topeka站城市比例降至67.5%，出现农业边缘带。500 m缓冲带捕捉到最显著的空间分化，而1000 m和2000 m缓冲带内各站点土地覆盖排序趋于一致。
**3.2 水质的季节变化**
基于月均值时间序列和独立样本t检验，温度和DO在所有站点均表现出最显著且一致的季节对比：湿季平均温度22.2–23.1 °C，旱季8.4–8.5 °C（p<0.001）；旱季DO平均11.95–12.51 mg/L，湿季8.32–8.83 mg/L（p<0.001）。浊度在湿季显著升高（p<0.001），而流量仅在NR Lake Quivira站点有显著季节差异（p=0.04），其余三个站点均无显著差异，反映上游水库调节抑制了流量季节变幅。pH在旱季略高，但仅De Soto和NR Lake Quivira达到显著水平。
**3.3 监测站点间的空间变异**
单因素ANOVA显示，在10个参数-季节组合中，仅湿季DO存在显著站点间差异（F=3.12, p=0.04）。Tukey HSD事后检验表明，NR Lake Quivira站的湿季DO显著低于Wamego站（均值差0.32 mg/L, p=0.03）和De Soto站（均值差0.30 mg/L, p=0.04），反映最下游城市化站点因不透水表面径流热加载导致的DO抑制。其他参数（流量、温度、浊度、pH）在旱季和湿季均无显著空间差异，表明四个干流站点共享相似的水质背景，空间变异有限。
**3.4 土地覆盖模式对水质的影响**
RDA分别对每个季节-缓冲尺度组合进行分析，以四个站点为空间重复（n=4）。经VIF筛选后，500 m缓冲带保留农业和湿地两个预测变量，1000 m和2000 m缓冲带保留草本、农业和湿地三个变量。置换检验所有组合均未达到显著水平（p>0.05），这是由空间重复少导致的统计功效不足所致。在500 m缓冲带，农业和湿地共同解释60%的水质总方差（RDA1: 42–47%, RDA2: 14–18%）；在1000 m和2000 m缓冲带，三个变量解释87%的方差，且两个尺度结果相同，表明超过500 m后土地覆盖组成在站点间趋于同质。RDA双序图显示：农业箭头指向浊度，尤其在湿季和2000 m缓冲带；城市（开发）预测变量指向湿季DO降低（Topeka和NR Lake Quivira）；草本和湿地预测变量指向DO升高并远离浊度。这些趋势虽未通过显著性检验，但与非点源径流和热加载的物理机制一致。
**3.5 时间序列概览**
月均值时间序列（2019–2026）显示，流量和浊度在所有站点呈下降趋势，pH逐渐上升，DO和温度则表现为稳定的季节性循环，无方向性变化。趋势统计（Mann–Kendall检验）仅供上下文参考，因时间序列存在自相关，研究结论未依赖这些趋势。
**讨论与结论**
讨论部分指出，温度和DO的季节对比是生态显著性的核心：湿季高温降低DO溶解度，导致夏季氧亏，限制鱼类和敏感大型无脊椎动物栖息地。浊度升高主要由农业径流中的细颗粒物驱动，而流量季节差异的缺失（三站未显著）是水库调节系统的独特特征，这与未调节河流中强烈的季节流量峰值形成对比。空间上，仅湿季DO存在显著差异，证实了城市热加载的局部影响。多尺度RDA的探索性结果支持农业-浊度关联和城市-DO抑制关联，且500 m缓冲带是最具空间分辨力的管理尺度，更宽的缓冲带未提供额外信息。研究局限性包括：仅四个干流站点导致统计功效低且空间自相关；NLCD单一年份快照无法捕捉年内季节性地表条件变化（如农业裸土期与生长期）；未纳入营养盐和生物指标。结论强调：在水库调节河流中，季节性和流量调节主导水质变异，河岸LULC信号较弱；湿季DO是下游城市段最需关注的管理指标；500 m河岸缓冲带是最具操作性的管理目标，应加强城市段植被缓冲以缓解热加载，并在上游农业段维持保护性耕作措施。未来应增设支流站点并采用空间自相关感知模型（如混合效应模型或空间流网模型）以提高统计推断能力。
