1. 课程导入与流域问题界定
流域土地利用变化与污染负荷关联，研究的是一条从“人怎么用地”到“水里有多少污染物”的因果链。它既不是单纯的土地覆被遥感解译，也不是孤立的水质评价，而是把空间格局、水文过程与污染输移放在同一框架下，回答三个问题：土地怎么变、负荷怎么来、两者如何定量关联。本章作为全课程起点，先厘清概念、搭好框架，再示范如何从一个真实的水质问题倒推分析边界。
三个核心概念
土地利用变化指人类活动在地表覆盖上的改变，包括类型转换（如林地转园地、水田转建设用地）和强度变化（如复种指数提高、施肥量增加）。在污染负荷研究中，真正起作用的往往不是“面积占比”本身，而是变化的位置、方向与速率——坡地上的开发与河漫滩上的开发，对泥沙输出的影响可能相差一个数量级。
污染负荷指一定时段内进入受纳水体的污染物总量，常用单位是 kg/a 或 t/a，区别于浓度（mg/L）。负荷 = 浓度 × 水量，因此它同时受“源强”和“输移水量”控制。按来源可分为点源（排污口、集中养殖）与面源（农田氮磷、畜禽散养、城镇地表径流、水土流失）。面源负荷的关键特征是时空分散、受降雨驱动、难以直接计量，这正是需要借助模型与关联分析的原因。
关联分析是在土地利用格局与负荷响应之间建立可解释、可检验的定量关系。常用手段包括：在子流域尺度上做相关与回归（如建设用地占比与总氮负荷）、用输出系数模型做负荷估算、用 SWAT 等分布式模型做情景模拟，以及用地理加权回归识别空间异质性。关联不等于因果，但好的关联分析应能回答“若某类用地增加 10%，负荷大致变化多少”。
整体分析框架
本课程采用“问题界定—数据与尺度—负荷估算—关联建模—情景与决策”的主线。下图展示了各环节的输入、输出与反馈关系，其中“尺度选择”和“源汇识别”是两个容易被忽视却决定结论可靠性的枢纽环节。
从实际问题出发界定边界
假设某水库型水源地连续三年夏季总磷超标，藻华面积扩大。面对这样的问题，第一步不是立刻收集数据，而是把模糊现象转成可分析的问题。可以按以下顺序推进：
- 锁定污染物与响应指标。总磷是藻华的限制因子，故以总磷负荷为主指标，叶绿素 a 为响应指标；若同时存在氨氮超标，则需区分生活源与农业源。
- 确定时间窗口。藻华暴发在 6—9 月，但负荷积累贯穿全年。建议以“水文年”为基本单元，并单独分析汛期（5—9 月）与非汛期，因为汛期面源贡献占比通常显著上升。
- 划定空间边界。以水库直接汇水区为分析单元，按 DEM 提取子流域。若上游有跨区来水，需明确“区内产生”与“入境输入”的区分，否则负荷归因会失真。
- 识别关键问题。把现象转成可检验的假设，例如“近十年库周建设用地增加导致汛期磷负荷上升”“坡耕地占比高的子流域贡献率显著更高”。
一个可操作的起步示例
以某 320 km² 汇水区为例（示例假设数据）：土地利用为林地 52%、耕地 31%、建设用地 9%、其他 8%；设有 1 个入库水质站、3 个雨量站，2015—2023 年逐月总磷监测。起步阶段可先做两件事：一是用输出系数法粗算年负荷，林地取 0.15、耕地取 1.8、建设用地取 1.2（单位 kg/(hm²·a)，均为示例取值，实际须按本地文献或监测率定）；二是计算各子流域建设用地与耕地占比，与汛期总磷浓度做散点相关。若相关系数偏低，优先检查子流域划分是否合理、监测点是否受点源干扰，而不是急于换更复杂的模型。