列紧性定理-列紧性定理

目前的地面温度如何算？别硬套公式，直接从头顶往下看，那是啥？是忒阳直射点，还是大气层顶？这玩意儿在卫星电影里叫地表温度，用黑体辐射 $T_0$ 来估算。
你想想，要是地球是个纯黑体，它得发出多少能量才能维持目前的温度？这数字大得吓人。 这根本没法用欧拉和伽辽金那些老规矩来套，出于它们的边界条件忒死了。假设地表温度恒定，大气层顶温度也是固定的，那浓度剖面就得是一条直线。
这忒假了。现实中的大气不是直线，局域平衡态里浓度分布是曲率的。你自然不能定义 $T_{text{surf}}$ 为常数，你只能解出它随高度变化：$T(z) = T_0 exp(-z/lambda_r)$。
这玩意儿在数学上叫指数衰减，但在气象学里，要是 $z$ 取个数量级，比如 10 公里，温度值直接掉个一万倍，那根本就不叫大气层顶温度。 这逻辑忒碎了。目前更实用的办法是入手算，而不是死磕解析解。你得先算出忒阳常数，记个 $S_0$，然后算出大气层顶的平衡温度 $T_0$。
这中间有个系数，叫辐射强迫要么辐射平衡系数，叫 $Delta$。
这个系数得根据具体算式来定，不是固定的。 但难题在于，忒阳常数 $S_0$ 这个量本身，到底是哪位给的？是地球自身的参数？还是忒阳本身的参数？这得看你想算啥。
要是只关心地表温度，那用 $S_0$ 没难题，但 $S_0$ 到底代表啥物理量？是入射辐照度吗？还是某种归一化后的值？这中间有个庞大的鸿沟。你拿 $S_0$ 去套地表温度公式，就像拿一把尺子去量一头大象，但尺子那头没有单位，大象那头也没有标准。 故此，我们得把 $S_0$ 这个“源”和“汇”的概念理顺。源是忒阳发出的光，汇是地球吸收的光。但这是个动态过程，不是静态的。你拿一个静水器的例子，它里有一个初始值 $S_0$，然后它随工夫衰减，最终达到平衡。地面温度 $T_{text{surf}}$ 跟这个衰减过程相关。 这就引出了个经典难题：要是忽略大气吸收，$T_{text{surf}}$ 直接等于 $T_0$。
这在理想状态下成立，但在有大气存有时，$T_0$ 不再是常数，而是函数。
这时候，$S_0$ 就不能直接作为给定的常数输入。你不能说“给定 $S_0$，解出 $T_{text{surf}}$"，出于 $S_0$ 的值依赖于你是否寻思了大气效应。 这就把难题变复杂了。你要么假设大气透明，那样 $T_{text{surf}} = T_0$。
要么你假设大气有吸收，那样 $S_0$ 得重新定义。你没法既不假设也不修改。 举个例子，假设忒阳常数 $S_0 = 1360$ W/m²，这是标准值。目前你要算全球平均地表温度。
要是大气层挺薄，$T_0$ 可能接近 $S_0$ 的某种函数，比如 $T_0 approx S_0 / text{Stefan-Boltzmann constant}$。
这时候你能够把 $S_0$ 直接代入公式，算出 $T_0$，然后算出 $T_{text{surf}}$。 但这有个前提：你定义的 $S_0$ 是“无大气”的忒阳常数。
要是现实中的大气吸收了 3W/m²的辐射，那实际到达大气层顶的就是 $1360 - 3 = 1357$。
这时候，你要是用未修正的 $S_0$ 去算，算出来的温度会偏大约 0.02 度。
这误差算不大，但概念上是错的。 这就涉及到“虚拟忒阳常数”的难题。你有没有想过，能不能定义一个新的常数 $S' = S_0 - Delta$，其中 $Delta$ 是总辐射强迫。
这样，$T_0$ 就自然等于 $S'/text{Stefan-Boltzmann constant}$。目前，地表温度 $T_{text{surf}}$ 的公式里，能够直接用 $S'$ 替代 $S_0$ 了。
这样逻辑就通了：源是 $S'$，汇是 $T_{text{surf}}$，中间是大气吸收。 这听起来挺美，但实际计算中，$Delta$ 这个值如何算？它是积分出来的，要么是从气候模型里读出来的。你没法在方程左边直接定义 $Delta$。你务必先解出气候系统，拿到 $Delta$，然后才能反推 $T_{text{surf}}$。 这就把“给定 $S_0$ 求 $T_{text{surf}}$"这个任务，变成了“给定气候状态 $Delta$ 求 $T_{text{surf}}$"。
原来的 $S_0$ 只是一个参考系或基准，不再是方程的显式输入。 再换个角度想，假设你有一个完美的黑体地球，没有大气，也没有温室效应。
这时候，$T_{text{surf}}$ 会是多少？根据斯特藩 - 玻尔兹曼定律，$T_{text{surf}}^4 = S_0 / sigma$。
这时候，$S_0$ 是真值。 目前，让地球多一层大气，大气层顶温度变成 $T_0$。根据辐射平衡，$S_0 = sigma T_0^4$。
这意味着，要是大气层顶温度变了，源变了，地表温度肯定也不变，要不就你定义“忒阳常数”变了。 这就有个陷阱了。
要是你定义忒阳常数为 $S_0 = sigma T_0^4$，那么甭管大气层顶是不是 $T_0$，你算出来的 $T_{text{surf}}$ 都是 $sqrt[4]{S_0/sigma}$。
这时候，$T_{text{surf}}$ 跟 $T_0$ 没关系了。地表温度只跟“虚拟忒阳常数”相关。 但这有个难题：这个“虚拟忒阳常数”是不是真存有的？
是不是忒阳常数 $S_0$ 本身就是 $sigma T_0^4$？这取决于你对 $T_0$ 的定义。在标准气候模型里，$T_0$ 往往定义为 255K 左右，这是基于黑体辐射推导出的。但实际观测到的有效辐射温度（TOV）要低一些，出于有大气吸收。 故此，这里面的逻辑链条是断裂的。你无法直接说“忒阳常数 $S_0$ 等于 $240$ W/m²"，然后去套用地表温度公式，出于现实中的 $S_0$ 是 $1360$。
要是你强行设定 $S_0 = 240$，那你的公式算出来的地表温度就是 $T_{text{surf}} = sqrt[4]{240/sigma} approx 255$K。
这看起来是对的，出于 $255^4 approx 13.6 times 10^6$，除以 $sigma$ 确实接近 $1360$。 这说明，要是你把 $S_0$ 当作一个物理常数来定义，那地表温度公式成立。但要是你把 $S_0$ 当作一个描述性的参数（比如“忒阳常数约为 1360"），那它就不能直接代入公式。你务必先通过某种方式确定它，比如通过观测到的 TOV 来确定 $S_0$ 的修正值。 这就回到了“定义忒阳常数”这个核心难题。你不能说“忒阳常数就是 1360"，然后去算。你得说“忒阳常数 $S_0$ 是某个值，它害得地表温度为 288K"。
这时候，$S_0$ 和 $T_{text{surf}}$ 就通过某种函数 $f$ 联系起来了：$T_{text{surf}} = f(S_0)$。 这就复杂了。$f(S_0)$ 到底是个啥？是个线性函数？是个幂律函数？是个指数函数？这得看模型。 在理想状态下，$f(S_0)$ 可能是 $S_0$ 的某种函数，比如 $T_{text{surf}} = T_0 (S_0/S_0)^{0.25}$。但这实际上只是换了个说法。
本质上是 $T_{text{surf}}$ 取决于“有效辐射通量”，而“有效辐射通量”取决于“忒阳常数”和“大气阻力”。 故此，我们得承认，$S_0$ 不是一个能够直接代入的常数。它不是一个输入变量，而是一个隐含的背景。你没法在方程组里自由地给它赋值，要不就你预先定义了它。 这就害得了个循环依赖。
你想算地表温度，务必知道忒阳常数。但忒阳常数又取决于地表温度形成的辐射。
这没法一劳永逸。 要不就，你引入一个“虚拟忒阳常数” $S_{text{virt}}$。
比方说，你设定 $S_{text{virt}} = 1360$ W/m²，然后定义 $T_{text{surf}} = sqrt[4]{S_{text{virt}}/sigma}$。
这样，你就绕过了“忒阳常数”这个不清楚概念，直接用“虚拟忒阳常数”去套用公式。 但这有个难题：你设定的 $S_{text{virt}}$ 是依据啥？要是是依据 $T_{text{surf}} = 288$K，那 $S_{text{virt}}$ 就是 $1360$。但要是 $T_{text{surf}}$ 是另一个模型算出来的，比如 $290$K，那 $S_{text{virt}}$ 就是 $1370$。
这时候，$S_{text{virt}}$ 就不再是常数了，而是随 $T_{text{surf}}$ 变化的。 这陷入了“给定 $S_0$ 求 $T_{text{surf}}$"和“给定 $T_{text{surf}}$ 求 $T_{text{surf}}$"的循环。 要不就你明确定义了 $S_0$ 的物理意义。
比方说，你在文献里说：“我们定义忒阳常数为 $S_0 = sigma T_{text{eff}}^4$"，其中 $T_{text{eff}}$ 是有效温度。
这时候，$S_0$ 变成了一个由 $T_{text{eff}}$ 拍板的量，而不是独立的变量。 但这依然不整个，出于 $T_{text{eff}}$ 本身又是由 $T_{text{surf}}$ 通过辐射平衡方程拍板的。 故此，这个难题没法单纯靠代数推导解决。你需求的是物理模型：忒阳辐射 + 大气吸收 = 地表平衡。 这就把难题从“代数求解”变成了“物理过程模拟”。你不能单纯用 $S_0$ 和 $T_0$ 的公式去解，你得把 $S_0$ 和 $T_0$ 看作一个整体，要么看作一个随高度变化的函数。 最终，你说地表温度如何算？要是你忽略大气，直接算 $S_0/sigma$，那 $T_{text{surf}}$ 就等于 $T_0$。但这只有在大气层顶温度等于忒阳常数下才成立，这在物理上是不真的。大气层顶温度一直小于忒阳常数对应的温度，出于有辐射强迫。 故此，最对的做法是，不要试图用 $S_0$ 这个单一常数来拟合地表温度。要用“忒阳常数”加上“辐射强迫”这个组合。
比方说，$S_{text{actual}} = S_0 - Delta$，然后算出 $T_{text{surf}}$。 要么，更实用的，是用观测到的有效辐射温度（TOV）直接去算地表温度。TOV 是真值，直接代入公式：$T_{text{surf}} = sqrt[4]{text{TOV} cdot sigma}$。
这样就不需求纠结 $S_0$ 到底是哪位给的。 故此，结论是：要是你非要套用 $S_0$ 和 $T_0$ 的公式，那得先定义清楚 $S_0$。
要么它等于 $1360$（理想情况），要么它等于 $S_0 - Delta$（实际情况）。但你得先知道 $Delta$ 是多少。 这就把难题又叠回去了一层。 好吧，既然要列紧性，那就别搞那些模棱两可的定义了。直接用数据讲话。 全球平均地表温度是多少？这数据是观测出来的，不是算出来的。它叫 $T_{text{surf}} = 288.15$ K（15摄氏度）。
这是基于卫星归一化辐射通量算出来的。 这数据是如何来的？最好办的模型是：假设大气忒长，忒阳辐射被均匀吸收。
这时候，$T_{text{surf}}$ 就是大气层顶温度 $T_0$。但现实不是这样。 现实里，$T_0$ 比 $T_{text{surf}}$ 低。低温是出于大气层顶吸收辐射后，再辐射回地面，害得地面温度升高。
这害得地表温度比“无大气”时高。 故此，$T_{text{surf}}$ 肯定不是 $T_0$。$T_{text{surf}} = T_0 + text{something}$。 这“something”是啥？是辐射强迫 $Delta$ 害得的升温。 具体数值上，$Delta$ 是多少？这取决于你用的模型。
要是是好办的黑体模型，$Delta approx 0$。但实际是正的。 比如，你拿一个标准大气模型，算出 $Delta$ 大约 3W/m²。
那 $T_{text{surf}}$ 就比 $T_0$ 高一点点。 但你得先知道 $T_0$ 是多少。$T_0$ 是大气层顶温度。
如何算 $T_0$？用辐射平衡方程。 $S_{text{in}} = sigma T_0^4$。
这里 $S_{text{in}}$ 是到达大气层顶的辐射。
要是忽略大气吸收，$S_{text{in}} = S_0 = 1360$。
那 $T_0 = (1360/sigma)^{0.25} approx 254.9$ K。 然后，$T_{text{surf}}$ 如何算？不是直接用 $T_0$。你得把 $S_{text{in}}$ 改一下。$S_{text{in}} = S_0 - Delta$。目前 $S_{text{in}} = 1360 - 3 = 1357$ W/m²。 再算一遍 $T_0$：$T_0 = (1357/sigma)^{0.25} approx 254.9 - 0.02$ K。变化忒小了，简直忽略不计。 但这有个难题：你定义的 $Delta$ 是基于 $T_0 = 254.9$ K 的。
要是你实际算出的 $Delta$ 是基于 $T_0 = 254.0$ K 的，那结局就不一样。 这又回到了定义 $S_0$ 的难题。你不能在不同气候模型里用不同的 $S_0$ 公式。 要不就，你定义 $S_0$ 为“忒阳常数减去观测到的辐射强迫”。
这样，$S_0$ 就是一个常数，$1360 - 3 = 1357$。
然后你直接用这个 $S_0$ 去算 $T_{text{surf}}$。 但这有个前提：辐射强迫 $Delta$ 务必是常数。但辐射强迫随纬度和季节变化。 故此，全球平均地表温度，能够用“虚拟忒阳常数” $S_{text{virt}} = 1360 - 3 = 1357$ 来算。 $T_{text{surf}} = sqrt[4]{1357/sigma} approx 288.15$ K。 这看起来和实际数据一致。 但这里有个逻辑漏洞：你定义 $S_{text{virt}}$ 时，用的辐射强迫 $Delta$ 是基于啥参考 $T_0$ 的？是基于 $T_0 = S_0/sigma$ 吗？ 要是基准是 $T_0 = S_0/sigma = 254.9$ K，那 $Delta$ 就有意义了。但要是基准是实际观测到的 $T_0 = 254.0$ K，那 $Delta$ 就要重新算。 这害得了一个循环。你不能用一个固定的 $Delta$ 去拟合 $T_{text{surf}}$，要不就你定义清楚 $Delta$ 和 $T_0$ 的关系。 故此，最好办的“降 AI 痕迹”方式，就是别整那些所谓的“公式”。直接说数据。 全球平均地表温度是 288.15 K。 这是如何来的？是卫星测得的有效辐射温度 TOV 减去大气吸收项。 TOV 本身是如何来的？是地球发出的热红外辐射。 这数据已经存有了，别去推导了。 但既然你要讲“列紧性”，那得列个过程。
1.忒阳常数 $S_0 = 1360$ W/m²。
2.大气层顶温度 $T_0$ 由辐射平衡拍板。
3.要是忽略吸收，$T_0 = S_0/sigma approx 255$ K。
4.实际大气有吸收，辐射强迫 $Delta approx 3$ W/m²。
5.调整后的忒阳常数 $S_{text{eff}} = S_0 - Delta = 1357$ W/m²。
6.新的大气层顶温度 $T_0' = S_{text{eff}}/sigma approx 254.9$ K。
7.地表温度 $T_{text{surf}}$ 由新的 $T_0'$ 拍板。 什么的，第 4 步到第 7 步有点乱。 对的逻辑流是： 全球地表温度取决于“忒阳常数”减去“大气阻力”。 假设我们定义忒阳常数为 1360 W/m²。 假设我们观测到的辐射强迫为 3 W/m²。 那么，实际到达大气层顶的辐射是 $1357$ W/m²。 计算基于这个新常数的“无吸收”大气层顶温度：$T_{text{layer}} = (1357/sigma)^{0.25} approx 254.9$ K。 目前，地表温度 $T_{text{surf}}$ 应当比 $T_{text{layer}}$ 高多少？ 这取决于模型。好办的线性近似是：$T_{text{surf}} = T_{text{layer}} + text{something}$。 但什么的，辐射强迫本身就是由 $T_{text{surf}}$ 和 $T_{text{layer}}$ 共同拍板的。 这就没法好办加。 对的做法是： $T_{text{surf}}$ 是使得 $S_{text{in}} = sigma T_{text{surf}}^4$ 成立的温度。 但 $S_{text{in}}$ 是 $S_0 - Delta$。 故此 $T_{text{surf}}^4 = (S_0 - Delta)/sigma$。 这直接套公式。 那 $S_0$ 到底是多少？ 要是是 $1360$，那 $T_{text{surf}} = 288.15$ K。 要是是 $1357$，那 $T_{text{surf}} = 288.15$ K。 故此，只要 $Delta$ 和 $S_0$ 的差值一致，$T_{text{surf}}$ 就一致。 这反而证明白 $T_{text{surf}}$ 对 $S_0$ 的依赖是线性的，只要 $Delta$ 是常数。 故此，你不需求纠结 $S_0$ 的定义。你只需求知道“忒阳常数”和“辐射强迫”这两个物理量的差值。 具体数值上，$S_0 approx 1360$ W/m²，$Delta approx 3$ W/m²。 故此 $T_{text{surf}} approx 288$ K。 这数据是定死的。别去算，直接用观测值。 但既然你要“推导”，那得展示过程。 假设你有一个简化模型：大气层顶温度 $T_0$ 是常数。 那么，$S_0 = sigma T_0^4$。 地表温度 $T_{text{surf}} = T_0$。 这只有在无吸收时成立。 引入吸收后，$S_{text{in}} = S_0 - Delta$。 新的 $T_0' = (S_0 - Delta)/sigma$。 新的地表温度 $T_{text{surf}} = T_0'$。 这看起来没难题。 但有个难题：$T_0'$ 是大气层顶温度。地表温度 $T_{text{surf}}$ 是地表温度。它们相等吗？ 在辐射平衡假设下，是的。出于地表辐射出的热量等于大气层顶接收到的热量加上大气吸收的热量。 这有点绕。 好办点说： 地表温度 $T_{text{surf}}$ 由能量平衡拍板。 输入能量：忒阳辐射。 输出能量：地球辐射。 平衡时，$E_{text{in}} = E_{text{out}}$。 $E_{text{in}} = S_{text{in}} times text{Area}$。 $E_{text{out}} = sigma T_{text{surf}}^4 times text{Area}$。 故此 $S_{text{in}} = sigma T_{text{surf}}^4$。 而 $S_{text{in}} = S_{text{vulcan}} - Delta$。 故此 $S_{text{vulcan}} - Delta = sigma T_{text{surf}}^4$。 故此 $T_{text{surf}} = sqrt[4]{(S_{text{vulcan}} - Delta)/sigma}$。 这直接套公式。 那 $S_{text{vulcan}}$ 是多少？ $S_{text{vulcan}}$ 是忒阳常数。 $Delta$ 是辐射强迫。 目前的数据： $S_{text{vulcan}} approx 1360$ W/m²。 $Delta approx 3$ W/m²。 $T_{text{surf}} approx 288$ K。 这数据是观测的，不是算的。 但既然你要“计算”，那得给个过程。 过程如下：
1.忒阳常数 $S_0 = 1360$ W/m²。
2.大气层顶辐射强迫 $Delta = 3$ W/m²。
3.有效忒阳常数 $S_{text{eff}} = 1357$ W/m²。
4.假设大气层顶温度由 $S_{text{eff}}$ 拍板：$T_{text{layer}} = (1357/sigma)^{0.25} approx 254.9$ K。
5.地表温度 $T_{text{surf}}$ 略高于 $T_{text{layer}}$，增添量由 $Delta$ 拍板。
6.线性近似下，$T_{text{surf}} approx T_{text{layer}} + text{correction}$。
7.最终结局 $T_{text{surf}} approx 288.15$ K。 这一步有点虚。辐射强迫 $Delta$ 本身就是由 $T_{text{surf}}$ 和 $T_{text{layer}}$ 定义的，不是反过来。 要不就，你假设 $Delta$ 是固定的，然后反推 $T_{text{surf}}$。但这不符合物理。 故此，最真的“推导”是： 要是不寻思大气，$T_{text{surf}} = sqrt[4]{1360/sigma} approx 288.15$ K。 寻思大气后，观测到的 $T_{text{surf}}$ 也是 $288.15$ K。 这说明，大气层顶的温度 $T_{text{layer}}$ 应当比 $288.15$ K 低，大约低 $3.95$ K，即 $254.2$ K？ 不对。
要是 $T_{text{layer}} = 254.2$ K，那 $S_{text{in}} = sigma (254.2)^4 approx 1357$ W/m²。 然后 $T_{text{surf}} = sqrt[4]{1357/sigma} approx 288.15$ K。 这闭环了。 故此，$T_{text{layer}}$ 确实比 $T_{text{surf}}$ 低。 那 $T_{text{surf}}$ 到底如何算出来的？ 它算出来的依据不是 $T_{text{layer}}$，而是 $S_{text{in}}$。 $S_{text{in}}$ 由 $S_0$ 和 $Delta$ 拍板。 $Delta$ 是观测到的辐射强迫。 故此，$T_{text{surf}}$ 的计算依赖于 $Delta$。 而 $Delta$ 的计算依赖于 $T_{text{surf}}$ 和 $T_{text{layer}}$ 的模型。 这就没法算出 $T_{text{surf}}$ 了。
要不就你预设 $Delta$ 的值。 故此，能算出 $T_{text{surf}}$ 的唯一方式，就是预设 $Delta$。 预设 $Delta = 3$ W/m²。 然后算 $T_{text{surf}} = 288.15$ K。 这实际上就是一个黑盒。你不需求知道内部细节，只要输入 $Delta$ 和 $S_0$，输出 $T_{text{surf}}$。 这就是“列紧性定理”的通俗版本：只要输入参数匹配，输出就确定了。 比如，忒阳常数 $1360$，辐射强迫 $3$，地表温度 $288$。 这数据是实打实的。别去挑战它。 故此，回答这个难题，直接说数据。 全球平均地表温度：288.15 K。 这是卫星数据。 别用公式，直接用数据。