[论文解读] Analysis of a two-layer energy balance model: long time behaviour and greenhouse effect

Piermarco Cannarsa, Valerio Lucarini+3 more|arXiv (Cornell University)|Nov 28, 2022

Earth Systems and Cosmic Evolution参考文献 56被引用 4

一句话总结

本文分析了一个双层能量平衡模型，通过大气吸收率 $\varepsilon_a$ 描述温室效应。研究证明：当 $\varepsilon_a \in (0,2)$ 时，解全局存在；当 $\varepsilon_a > 2$ 时，解在有限时间内爆破；当 $\varepsilon_a \in (0,2)$ 时，解收敛至平衡态，且地表温度随 $\varepsilon_a$ 单调递增，从数学上确认了温室效应。

ABSTRACT

We study a two-layer energy balance model, that allows for vertical exchanges between a surface layer and the atmosphere. The evolution equations of the surface temperature and the atmospheric temperature are coupled by the emission of infrared radiation by one level, that emission being captured by the other layer, and the effect of all non radiative vertical exchanges of energy. Therefore, an essential parameter is the absorptivity of the atmosphere, denoted $ε_a$. The value of $ε_a$ depends critically on greenhouse gases: increasing concentrations of $CO_2$ and $CH_4$ lead to a more opaque atmosphere with higher values of $ε_a$. First we prove that global existence of solutions of the system holds if and only if $ε_a \in (0, 2)$, and blow up in finite time occurs if $ε_a > 2$. (Note that the physical range of values for $ε_a$ is $(0, 1]$.) Next, we explain the long time dynamics for $ε_a \in (0, 2)$, and we prove that all solutions converge to some equilibrium point. Finally, motivated by the physical context, we study the dependence of the equilibrium points with respect to the involved parameters, and we prove in particular that the surface temperature increases with respect to $ε_a$. This is the key mathematical manifestation of the greenhouse effect.

研究动机与目标

• 严格分析具有耦合地表与大气温度的双层能量平衡模型的长期行为。
• 确定解全局存在或在有限时间内爆破的条件。
• 建立平衡态对大气吸收率 $\varepsilon_a$ 的依赖关系，$\varepsilon_a$ 是温室效应的关键参数。
• 从数学上确认地表温度随 $\varepsilon_a$ 单调递增，从而体现温室效应。

提出的方法

• 构建包含地表与大气温度耦合演化方程的双层能量平衡模型。
• 引入长波辐射交换与非辐射能量通量，由大气吸收率 $\varepsilon_a$ 参数化。
• 利用能量平衡原理与斯特藩-玻尔兹曼辐射定律建模向外长波辐射。
• 应用动力系统理论研究渐近行为及收敛至平衡点的性质。
• 通过严格凸性与罗尔定理分析平衡方程 $\Phi(T_s) = q\beta_s(T_s)$，以界定解的个数。
• 对函数 $\Phi$ 进行数值与解析研究，包括其导数，以确定不同 $\varepsilon_a$ 值下的平衡态个数。

实验结果

研究问题

• RQ1 在何种 $\varepsilon_a$ 条件下，双层能量平衡模型的解全局存在？
• RQ2 当 $\varepsilon_a \in (0,2)$ 时，系统的长期行为如何？所有解是否都收敛至平衡点？
• RQ3 平衡解的个数如何依赖于 $\varepsilon_a$ 的取值，特别是在临界阈值 $\varepsilon_a = 2$ 附近？
• RQ4 地表温度是否随 $\varepsilon_a$ 单调递增？这是否为温室效应提供了数学基础？
• RQ5 该模型是否能支持多个平衡态？若能，其参数条件为何？

主要发现

• 双层能量平衡模型的解全局存在当且仅当 $\varepsilon_a \in (0,2)$；当 $\varepsilon_a > 2$ 时，解在有限时间内爆破，尽管物理上 $\varepsilon_a \in (0,1]$。
• 对所有 $\varepsilon_a \in (0,2)$，所有解均收敛至某个平衡点，表明在物理相关条件下系统具有稳定长期行为。
• 地表温度随 $\varepsilon_a$ 单调递增，为温室效应提供了严格的数学确认。
• 平衡解的个数最多为六个，且在物理相关参数下，数值证据表明最多仅有三个平衡态存在，即使在 $\varepsilon_a \to 2^-$ 附近亦如此。
• 当 $\varepsilon_a \in (0,1]$ 时，函数 $\Phi$ 严格凸且严格递增，确保最多存在三个平衡解。
• 当 $\varepsilon_a \in (\varepsilon_{a,0}, 2)$ 且 $\varepsilon_{a,0} \approx 1.99$ 时，$\Phi$ 呈现凸-凹-凸形状，最多允许五个平衡解，但数值结果表明仅有最多三个在物理上相关。