Mathematical Modelling of Water and Carbon Cycles | 水与碳循环的数学建模

📚 Mathematical Modelling of Water and Carbon Cycles | 水与碳循环的数学建模

Water and carbon cycles are fundamental Earth system processes that regulate climate, sustain ecosystems, and drive biogeochemical dynamics. In A-Level Mathematics, differential equation models provide a powerful tool for quantifying these cycles, enabling us to predict reservoir changes, understand feedback mechanisms, and solve real‑world environmental problems. This article builds a bridge between physical geography concepts and mathematical techniques, covering box models, linear reservoirs, coupled systems, steady‑state analysis and numerical methods, all anchored in the Edexcel specification’s emphasis on modelling with differential equations.

水循环和碳循环是地球系统的基本过程，调控气候、维持生态系统并驱动生物地球化学动态。在A-Level数学中，微分方程模型为我们量化这些循环提供了强有力的工具，使我们能够预测储库变化、理解反馈机制并解决现实环境问题。本文将自然地理概念与数学技巧连接起来，涵盖箱式模型、线性水库、耦合系统、稳态分析及数值方法，紧扣Edexcel大纲对微分方程建模的要求。

1. Overview of Water Cycle Dynamics | 水循环动力学概述

The global water cycle consists of reservoirs such as oceans, ice caps, groundwater, lakes and the atmosphere, linked by fluxes including evaporation, precipitation, runoff and infiltration. A simple mathematical statement is the conservation of mass: the rate of change of water stored in a catchment equals total input minus total output. This can be expressed as dW/dt = P − E − R, where W is water storage, P is precipitation rate, E is evapotranspiration rate and R is runoff rate.

全球水循环由各大储库（海洋、冰盖、地下水、湖泊和大气）组成，并通过蒸发、降水、径流和入渗等通量连接起来。一个简单的数学描述是质量守恒：流域内储水量的变化率等于总输入减去总输出。这可用微分方程 dW/dt = P − E − R 表示，其中 W 为储水量，P 为降水率，E 为蒸散发率，R 为径流率。

Typical units are mm day⁻¹ or km³ year⁻¹, and modellers often treat the catchment as a single control volume. The continuity equation forms the backbone of many hydrological models.

典型单位是毫米每天或立方千米每年，建模者常将流域视为一个控制体。连续性方程构成了大量水文模型的骨架。

---

2. Carbon Cycle Reservoirs and Fluxes | 碳循环储库与通量

The carbon cycle involves four active reservoirs: atmosphere, terrestrial biosphere, oceans and lithosphere. Fluxes include photosynthesis, respiration, ocean uptake, outgassing and fossil fuel combustion. For mathematical modelling, we label carbon masses C₁, C₂, C₃, C₄ and define transfer coefficients kᵢⱼ representing the fraction of carbon moving from reservoir j to reservoir i per unit time. Then dC₁/dt = k₁₂C₂ + k₁₃C₃ + … − (k₂₁ + k₃₁ + …)C₁.

碳循环涉及四个活跃储库：大气、陆地生物圈、海洋和岩石圈。通量包括光合作用、呼吸作用、海洋吸收、释放和化石燃料燃烧。在数学建模中，我们标定碳质量 C₁, C₂, C₃, C₄，并定义迁移系数 kᵢⱼ，表示单位时间从储库 j 转移到储库 i 的碳比例。那么 dC₁/dt = k₁₂C₂ + k₁₃C₃ + … − (k₂₁ + k₃₁ + …)C₁。

Carbon cycle models often span timescales from years to millennia. The pre‑industrial cycle was in a quasi‑steady state, whereas anthropogenic emissions have created a transient forcing described by a non‑autonomous differential equation.

碳循环模型通常跨越几年到几千年的时间尺度。前工业化时期循环处于准稳态，而人为排放产生了一种由非自治微分方程描述的瞬态强迫。

---

3. Mathematical Framework: Box Models | 数学框架：箱式模型

A box model assumes each reservoir is well mixed so that outgoing fluxes are proportional to the reservoir’s content. If a reservoir has mass M(t) and total outflow coefficient λ, then dM/dt = I(t) − λM(t), where I(t) is the input flux. This is a first‑order linear differential equation. Its homogeneous solution gives exponential decay: M(t) = M₀e⁻λᵗ when I(t)=0.

箱式模型假设每个储库充分混合，因此输出通量与储库含量成正比。若储库质量为 M(t)，总输出系数为 λ，则 dM/dt = I(t) − λM(t)，这是 一阶线性微分方程。其齐次解为指数衰减：当 I(t)=0 时，M(t)=M₀e⁻λᵗ。

The mean residence time is τ = 1/λ, one of the most informative diagnostics in cycle analysis. For the global atmospheric water reservoir (≈13×10¹² m³) with a total precipitation flux of about 5×10¹⁴ m³ year⁻¹, τ ≈ 0.026 years, or about 9.5 days.

平均滞留时间 τ = 1/λ 是循环分析中最有用的诊断量。全球大气水储库约为 13×10¹² 立方米，总降水通量约 5×10¹⁴ 立方米每年，得出 τ ≈ 0.026 年，大约 9.5 天。

---

4. Simple Water Balance Model and Solution | 简单水量平衡模型及其解

Consider a lake of volume V(t) with constant inflow Q_in and outflow Q_out = αV. The water balance is dV/dt = Q_in − αV, α > 0. This is a linear constant‑coefficient equation. Setting an initial condition V(0)=V₀, the solution is V(t) = Q_in/α + (V₀ − Q_in/α)e⁻αᵗ. The steady‑state volume is V* = Q_in/α, and the lake relaxes to this equilibrium exponentially with time constant 1/α.

考虑一个体积为 V(t) 的湖泊，有恒定流入量 Q_in 和出流量 Q_out = αV。水量平衡为 dV/dt = Q_in − αV，α > 0。这是一阶常系数线性方程。设初始条件 V(0)=V₀，解为 V(t) = Q_in/α + (V₀ − Q_in/α)e⁻αᵗ。稳态体积 V* = Q_in/α，湖泊以时间常数 1/α 指数式趋近该平衡。

Such modelling appears in Edexcel A‑Level maths as an application of integrating factors or separation of variables. The algebraic steps reinforce skills in logarithmic integration and limit evaluation.

这种建模在Edexcel A-Level数学中作为积分因子或分离变量法的应用出现。代数步骤强化了对对数积分和极限求值的掌握。

---

5. Linear Reservoir Carbon Model | 线性水库碳模型

The terrestrial biosphere carbon stock C(t) gains carbon through net primary productivity (NPP) and loses it via heterotrophic respiration, proportional to C. Thus dC/dt = NPP − k C, identical in form to the water tank model. The general solution is C(t) = NPP/k + (C₀ − NPP/k)e⁻ᵏᵗ. If NPP increases due to CO₂ fertilisation, the steady state rises, but the adjustment follows an exponential trajectory.

陆地生物圈碳储量 C(t) 通过净初级生产力 (NPP) 获得碳，并通过异养呼吸（与 C 成正比）损失碳。因此 dC/dt = NPP − k C，形式与水槽模型相同。通解为 C(t) = NPP/k + (C₀ − NPP/k)e⁻ᵏᵗ。若 NPP 因 CO₂ 施肥效应增加，稳态上升，但调整过程遵循指数轨迹。

This model is a starting point for more complex representations that include soil carbon pools with different turnover times. Multiple boxes can be cascaded, leading to systems of differential equations.

该模型是更复杂表征的起点，后者包括具有不同周转时间的土壤碳库。可将多个箱子级联，形成微分方程组。

---

6. Coupled Water and Carbon Systems | 水碳耦合系统

Water and carbon cycles interact tightly: plant stomatal opening allows CO₂ uptake for photosynthesis while losing water vapour through transpiration. A simple coupled model can link soil moisture W and biomass carbon B: dW/dt = P − αW − βB, dB/dt = γW B/(W₀+W) − δB, where the photosynthetic rate saturates with water availability following a Michaelis‑Menten function.

水循环和碳循环紧密相互作用：植物气孔张开时吸收 CO₂ 进行光合作用，同时通过蒸腾作用损失水汽。一个简单的耦合模型可将土壤湿度 W 和生物量碳 B 联系起来：dW/dt = P − αW − βB, dB/dt = γW B/(W₀+W) − δB，其中光合速率随水分可用性按米氏函数饱和。

This nonlinear system cannot be solved by elementary methods; phase‑plane analysis and numerical simulations reveal steady states, oscillatory regimes or collapse depending on parameters. Such models are explored in the optional applied modules of advanced maths.

该非线性系统无法用初等解法求解；相平面分析和数值模拟可根据参数揭示稳态、振荡机制或崩溃。此类模型在进阶数学的应用选修模块中有所探究。

---

7. Steady‑State Analysis and Equilibria | 稳态分析与平衡点

For an autonomous system dX/dt = F(X), a steady state X* satisfies F(X*) = 0. In the water balance model dW/dt = P − E − R(W), setting R(W)=kW yields W* = (P−E)/k. Stability is assessed via the derivative: since d(dW/dt)/dW = −k < 0, the equilibrium is stable. In the coupled model, the Jacobian matrix determines whether small perturbations grow or decay.

对于自治系统 dX/dt = F(X)，稳态 X* 满足 F(X*) = 0。在水均衡模型 dW/dt = P − E − R(W) 中，设 R(W)=kW 得 W* = (P−E)/k。稳定性通过导数判断：因 d(dW/dt)/dW = −k < 0，平衡点稳定。在耦合模型中，雅可比矩阵决定小扰动是放大还是衰减。

Understanding equilibrium behaviour is essential for interpreting model projections, especially under climate change scenarios where forcing terms shift gradually. Students practise solving F(X)=0 for multiple linked equations, a skill also tested in pure maths.

理解平衡行为对解读模型预测至关重要，尤其是在气候变化情景下强迫项逐渐变化时。学生练习求解多个关联方程的 F(X)=0，这也是纯数学考试的一项技能。

---

8. Feedback Loops and Sensitivity Analysis | 反馈循环与灵敏度分析

Feedback occurs when a change in a variable alters a process that further modifies that variable. A positive feedback amplifies, a negative feedback dampens. In carbon cycle modelling, warming increases soil respiration rate k, reducing soil carbon. This can be expressed as k = k₀exp(βΔT), where ΔT is temperature anomaly. The system becomes dC/dt = NPP − k₀exp(βΔT) C, requiring numerical solution.

当某个变量的变化改变某一过程，进而进一步改变该变量时，就产生了反馈。正反馈起放大作用，负反馈起阻尼作用。在碳循环建模中，变暖会增大土壤呼吸速率 k，减少土壤碳。这可表示为 k = k₀exp(βΔT)，其中 ΔT 为温度异常。系统变为 dC/dt = NPP − k₀exp(βΔT) C，需要数值求解。

Sensitivity analysis quantifies how output uncertainty can be apportioned to different input parameters. Partial derivatives ∂C(t)/∂β, calculated analytically or by finite differences, identify the most influential parameters, a core element of mathematical modelling projects.

灵敏度分析量化输出不确定性如何归因于不同输入参数。偏导数 ∂C(t)/∂β 可通过解析或有限差分计算，从而确定最有影响的参数，这是数学建模课题的核心要素。

---

9. Numerical Simulation Techniques | 数值模拟技术

Euler’s method provides a straightforward way to integrate dX/dt = f(X,t). For a water reservoir model dW/dt = I(t) − λW, the update is W_{n+1} = W_n + h[I(t_n) − λW_n], with step size h. The global error is O(h), so small h is needed. Improved methods like the Runge‑Kutta 4th order (RK4) reduce error significantly and are mentioned in advanced textbooks.

欧拉法为积分 dX/dt = f(X,t) 提供了简单途径。对水库模型 dW/dt = I(t) − λW，其更新公式为 W_{n+1} = W_n + h[I(t_n) − λW_n]，步长为 h。全局误差为 O(h)，因此需要较小的 h。改进方法如四阶龙格-库塔法 (RK4) 显著减少误差，在进阶教材中有所提及。

Students may implement a simple Euler spreadsheet for seasonal precipitation inputs, observing how storage evolves. This links discrete mathematics and iteration logic with continuous modelling, a valuable revision exercise.

学生可以对季节性降水输入编写简单的欧拉法电子表格，观察储水量如何变化。这将离散数学和迭代逻辑与连续建模联系起来，是一项有价值的复习练习。

---

10. Parameter Estimation from Data | 基于数据的参数估计

Given observed reservoir masses and flux measurements, coefficients like λ or k can be estimated by minimising the sum of squared residuals between model output and data. If the analytical solution is C(t)=C₀e⁻ᵏᵗ, taking logarithms gives ln(C)=ln(C₀)−kt, a linear relationship. Linear regression then yields k, a technique often examined in the statistics component but applied here to differential equation models.

给定观察到的储库质量与通量测量值，可通过最小化模型输出与数据之间的残差平方和来估计 λ 或 k 等系数。若解析解为 C(t)=C₀e⁻ᵏᵗ，取对数得 ln(C)=ln(C₀)−kt，为线性关系。线性回归随即给出 k，这一技巧常在统计部分考查，但这里应用于微分方程模型。

For coupled nonlinear systems, numerical optimisation routines (e.g. gradient descent) are introduced conceptually, highlighting the synergy between calculus and statistics in environmental modelling.

对于耦合非线性系统，概念性地引入数值优化例程（如梯度下降），突出了环境建模中微积分与统计学的协同作用。

---

11. Exam‑Style Mathematical Problems | 考试风格的数学问题

Typical Edexcel A‑Level questions might present a differential equation for a carbon reservoir, ask to find the general solution via integrating factor, determine the steady state, and interpret the long‑term behaviour. For instance: “A peatland carbon stock obeys dC/dt = 2.5 − 0.02C, C in kg m⁻². Find C(t) given C(0)=80, and calculate the time for C to fall to 90% of its equilibrium value.”

典型的Edexcel A-Level题目可能给出关于碳储库的微分方程，要求使用积分因子求通解、确定稳态并解释长期行为。例如：“某泥炭地碳储量满足 dC/dt = 2.5 − 0.02C，C 单位为 kg m⁻²。已知 C(0)=80，求 C(t)，并计算 C 降至其平衡值 90% 所需的时间。”

Solution: Integrating factor e⁰·⁰²ᵗ gives C(t)=125 − 45e⁻⁰·⁰²ᵗ. Equilibrium is 125. 90% equilibrium is 112.5, solving 112.5 = 125 − 45e⁻⁰·⁰²ᵗ yields t = (ln(45/12.5))/0.02 ≈ 64.0 years. This integrates algebraic manipulation, exponentials and logarithms, core A‑Level competencies.

解答：积分因子 e⁰·⁰²ᵗ 得 C(t)=125 − 45e⁻⁰·⁰²ᵗ。平衡值为 125。90% 平衡值为 112.5，求解 112.5 = 125 − 45e⁻⁰·⁰²ᵗ 得 t = (ln(45/12.5))/0.02 ≈ 64.0 年。此题综合了代数运算、指数与对数，是A-Level核心能力。

---

12. Conclusion and Revision Tips | 结论与复习技巧

Mathematical models of water and carbon cycles consolidate techniques in differential equations, linear algebra, numerical methods and data analysis. As you revise, practise setting up conservation equations for reservoir systems, solving them analytically, and interpreting steady‑state and dynamic responses. Create flashcards linking physical processes to mathematical forms, and use technology to explore parameter sensitivity.

水循环和碳循环的数学模型整合了微分方程、线性代数、数值方法和数据分析等技巧。复习时，请练习为储库系统建立守恒方程、求解解析解并解释稳态和动态响应。制作闪卡将物理过程与数学形式联系起来，并利用技术探索参数灵敏度。

Remember that examiners reward clear structuring: state the equation, identify type, solve stepwise, and always link results back to the real‑world context. The synergy between mathematics and Earth science not only deepens understanding but also exemplifies how pure and applied mathematics unite to address global challenges.

请记住，考官青睐清晰的结构：陈述方程、判断类型、分步求解，并始终将结果联系回现实背景。数学与地球科学的协同不仅加深理解，也体现了纯数学与应用数学如何结合以应对全球挑战。

Published by TutorHao | Mathematics Revision Series | aleveler.com

更多咨询请联系16621398022（同微信）