本页目录
系统 I · 储库、通量与驻留时间
地球系统最小的可计算语言不是“几个圈”，而是储量、通量、边界和时间尺度。先问一份物质存在哪里，再问它以多快的速度进出；很多看似相反的地质、气候和水文直觉，都会在同一张账本上变得清楚。
学习层：让一个储库先接受审计
1. 具体情境：一座山地水库为何没有立刻追上入流
想象一个用于供水的山地储库。它当前有一笔可用水量，雨季平均入流把水送进来，蒸发、取水和下泄把水带走。工程师关心的不是“今年下了多少雨”这一句，而是：储量会增加还是减少？储库对扰动的记忆大约持续多久？如果入流突然改变，多久后才接近新的状态？
这里的数字先是教学设定，不是某个真实水库的预报，也不是观测报告。我们把所有水量写成同一套“储量单位”，把时间写成年；这样做是为了把守恒关系暴露出来，而不是制造一个虚假的现场精度。
2. 揭示前预测：先写方向，再看计算
打开实验前，先回答三件事：
- 若瞬时入流大于瞬时出流，储量的斜率是正还是负？
- 若出流与储量成正比，把出流系数减半，系统的驻留时间会怎样变化？
- 一个更大的初始储量，是否一定意味着它会更快达到新的平衡？
预测只描述模型内的方向。点击“揭示”后，实验会把预测和账本逐项对照；调参是在揭示后进行的，避免把场景输入误读成答案。
3. 公式桥：从守恒到时间响应
对储库储量 \(S(t)\)，最基本的守恒式是
若把输出近似成一阶线性过程 \(F_{\mathrm{out}}=kS\)，其中 \(k\) 的单位是 \(\mathrm{time}^{-1}\)，则
在输入恒定时，解为
这里的 \(S_*\) 是这个简化系统的稳态储量，\(S_0\) 是初始储量。最有用的时间尺度是
称为 e-folding 时间或一阶模型的驻留时间。经过 \(t=\tau\)，偏离稳态的部分还剩 \(e^{-1}\approx0.37\)；它不是“全部水分子平均待多久”的唯一答案，而是一个由模型定义的响应时间。
数量级要靠单位检查：若 \(S\) 是 \(10^3\) 个储量单位，通量是 \(10^2\) 个单位/年，\(\tau\) 常在 \(10\) 年量级；若同样的储量交换只发生在百万年尺度，\(\tau\) 就是 \(10^6\) 年量级。不能把两者的“通量很大”直接比较而忘了储量基数。
4. 互动实验：把储量、输入和记忆放进一张账本
无 JavaScript 时的静态 fallback：以下是教学设定：初始储量 \(S_0=1200\)，输入通量 \(F_{\mathrm{in}}=80\) 个单位/年，输出系数 \(k=0.05\) /年，观察窗口 \(t=10\) 年。于是 \(S_*=1600\)，\(\tau=20\) 年，窗口末储量约为 \(1357.4\)，输出通量约为 \(67.9\)，累计储量净变化约为 \(+157.4\) 个单位；窗口末瞬时净通量约为 \(+12.1\) 个单位/年。比例响应 \(1-e^{-kt}\approx0.393\)。
| 账本项 | 默认读数 | 单位 / 解释 |
|---|---|---|
| 输入通量 \(F_{\mathrm{in}}\) | 80.0 | 储量单位/年，教学输入 |
| 窗口末输出通量 \(kS(10)\) | 67.9 | 储量单位/年，随储量变化 |
| 稳态储量 \(S_*\) | 1600.0 | 储量单位，若输入保持不变 |
| 驻留时间 \(\tau\) | 20.0 | 年，一阶响应时间 |
| \(t=10\) 年储量 | 1357.4 | 储量单位 |
| 累计净变化 \(S(10)-S_0\) | +157.4 | 储量单位，10 年窗口的积分结果 |
| 窗口末瞬时净通量 | +12.1 | 储量单位/年，\(F_{\mathrm{in}}-kS(10)\)，不是累计变化 |
图中的水滴只表示储库与通量的关系，不能代表一个真实流域的面积、降水量或预报。
5. 边界：一个线性出流模型能解释到哪里
这个模型把输出写成储量的线性函数，隐含了许多条件：
- 输入在窗口内近似恒定；真实降水有季节性、脉冲和空间差异，取水也可能是阶跃或规则调度。
- \(k\) 是有效系数，不是一个跨尺度不变的自然常数。它可能随温度、植被、闸门、饱和度和管理规则改变。
- 一个储库往往需要多个子库。土壤水、深层地下水和冰雪的驻留时间不同，合成成一个 \(S\) 会隐藏快响应与慢响应的叠加。
- 稳态只在边界条件固定且系统有机会调整时有意义；持续外 forcing、相变或阈值会让“新稳态”不存在或不可达。
因此，账本能给出守恒、方向和数量级，却不能凭一条指数曲线替代水文观测、质量守恒校验和不确定度分析。
6. 迁移：把同一语法带到碳和热
碳库可以把 \(S\) 换成大气碳量，把输入写成排放，把输出写成海洋和陆地吸收；海洋热含量也可写成储量，表面净热通量是输入。公式形式相同，但 \(k\)、时滞和子库数量不同。
迁移时请回答：如果输入在 \(t=0\) 后增加一段时间再恢复，单一一阶模型会留下什么记忆？若两个储库以不同 \(\tau\) 串联，观测到的曲线为什么可能同时有快段和慢段？哪些额外观测能够区分“通量改变”和“储量估计偏差”？
1. 储库不是静态容器
“库”不是说物质永远停在那里，而是对指定边界和时间窗口做的聚合。大气中的水汽、土壤中的水、河流中的水和冰盖中的水可以被分别记账，也可以在更粗的水循环模型里合成一个总储量。边界一变，储量定义和通量定义也会一起变。
最先要写清楚的是系统边界：空间范围、物质种类、时间分辨率、是否包含相变、是否包含人为抽取。没有边界的“增加了很多”既不能和另一份数据相减，也不能判断守恒。
2. 通量是带方向的速率
通量常以质量/时间、体积/时间或能量/时间表示。对面积归一化后，则是质量/(面积·时间)或能量/(面积·时间)。同一过程在不同教材里可能采用不同符号，最稳妥的做法是先把正方向写在图上，再检查每一项的量纲。
净通量不是总通量。若两个方向各有 \(100\) 个单位/年，净通量可以是零，但交换并没有停止；这正是“稳态”与“没有过程”的区别。
3. 驻留时间与响应时间不是同义词
常用的库平均驻留时间可写成 \(\tau_{\mathrm{res}}=S/F_{\mathrm{out}}\)，但只有在近似稳态、输出定义清楚时，它才与一阶响应时间 \(1/k\) 相同。对多个子库、非线性输出或非稳态输入，平均年龄分布、脉冲响应和 e-folding 时间可能分别回答不同问题。
一个短驻留库可以有很大的通量，也可以对扰动快速响应；一个大而慢的库即使交换通量不小，也可能保存很长的系统记忆。比较系统时，必须同时报储量、通量、边界和时间尺度。
4. 证据纪律
在真实地球系统分析中，储量常由现场观测、卫星反演、再分析、质量守恒和模型同化共同约束。每个动态数字都应带指标年或时间窗、参考基准、不确定度和核验日期；没有这些元数据，不应把一列数当成可以跨年份直接比较的事实。
水循环的机制说明可优先查 NOAA 与 USGS 的水文资料，地球系统和长期气候背景可查 NASA 与 IPCC 的一手报告。这里实验的数字全是教学输入，未宣称来自某一观测产品。
5. 小结
一张合格的地球系统账本至少有四列：储量、输入、输出、时间尺度。先用 \(dS/dt=F_{\mathrm{in}}-F_{\mathrm{out}}\) 看方向，再用 \(S/F\) 或 \(1/k\) 判断记忆，最后把过程的边界和证据来源写出来。下一章把时间轴拉长：当地表记录被侵蚀、折叠或缺失时，怎样从深时证据重建先后顺序和绝对年龄？