跨海拔的城市分布：一种最大熵方法

《Environment and Planning B: Urban Analytics and City Science》：City distributions across elevation: A maximum-entropy approach

　　

世界上大多数城市集中在海平面附近，远超仅凭低地可用性所预测的数量。研究人员通过结合全球47,916个独立城市的海拔（每个城市由单个点位置表示）与匹配的国家高程学数据，并在最大熵（ME）框架内对海拔强度进行建模，量化了这一关系。在六个候选强度模型中，数据强烈支持凹成本设定优于线性成本指数模型。最佳AIC（赤池信息准则）拟合是一个具有陡峭低海拔状态和浅平高海拔状态的分段指数模型，而首选的平滑机制模型是幂律衰减，这意味着一个凹有效成本函数C(h) ∝ ln(1 + h/h₀)。此处h表示海拔高于海平面（以米为单位），而h₀是一个拟合尺度参数，设定了成本函数曲率集中的位置。针对城市海拔，齐普夫位频标度（Zipfian rank-frequency scaling）被拒绝，且主要推论对数据扰动和国家分层具有鲁棒性。交通证据作为物理背景保留，而非核心实证检验：规模经济和海上准入为陡峭的近海岸惩罚和更平缓的高海拔梯度提供了合理的机制，但并未识别出统计模式。

研究背景与意义
以往研究指出人类聚落更多集中在低海拔地区，McGranahan et al. (2007) 揭示了低海拔沿海带（LECZ, 低海拔沿海带, <10 m）仅占全球陆地面积的2%，却容纳了10%的全球人口和13%的城市人口。由于陆地运输成本远高于水路运输，高海拔定居点被排除在海运网络之外，导致严重的交通阻抗。尽管自Auerbach于1913年发现城市规模分布的齐普夫定律以来，研究者已在城市水平分布的多种现象中记录了类似模式，Gibrat定律的比例增长也被证明能产生幂律分布，但针对城市跨海拔分布的参数形式缺乏统计模型测试。海拔作为一个全球可比的空间梯度，能够系统性地反映土地可用性、气候暴露和交通准入的变化。开展这项研究旨在填补这一空白，寻找一个能总结城市在不同国家和地区背景下如何分布的稳定统计规律。Dunne开展了这项研究，论文发表在《Environment and Planning B: Urban Analytics and City Science》。该研究的重要意义在于揭示了城市垂直分布由不同于水平分布的统计结构主导，并证明了海拔累积准入惩罚的边际递减效应，为聚落地理学、交通经济学和最大熵推理之间的联系提供了一种紧凑的方式。

技术方法
研究人员收集了SimpleMaps世界城市数据库中的47,916个城市点位，提取了哥白尼全球数字高程模型（GLO-90, 哥白尼全球数字高程模型）的海拔数据，并结合全球人类居住层城市中心数据库（GHSL UCDB, 全球人类居住层城市中心数据库）的城市足迹数据。排除了海平面以下和海拔5000 m以上的记录后，采用匹配的各国高程学数据计算了不同海拔带的土地面积。在最大熵（ME, 最大熵）框架下，将海拔作为累积交通阻抗的代理变量，构建了以土地面积为先验的泊松点过程模型，比较了六种候选强度模型（包括指数、分段指数和幂律衰减等）。使用赤池信息准则（AIC, 赤池信息准则）和贝叶斯信息准则（BIC, 贝叶斯信息准则）进行模型选择，并进行了Clauset-Shalizi-Newman（CSN, Clauset-Shalizi-Newman）检验以验证齐普夫标度，辅以非参数自举法等稳健性检验。

研究结果
描述性模式
通过比较城市数量与土地面积随海拔的分布，发现城市数量在0-50 m处达到峰值并单调下降。在约500 m处，可用土地面积下降至0-50 m水平的60%，而城市数量下降至约30%。这表明城市随海拔减少的速度快于几何土地可用性的下降，证实了存在随海拔增加的成本惩罚。

强度模型比较
通过对比47,916个城市点过程的六种模型，所有凹成本模型均明显优于线性成本指数模型。AIC拟合最好的是分段指数模型，在断点h*=124 m以下具有陡峭的低海拔状态（每100 m下降78.5%），在断点以上具有浅平的高海拔状态（每100 m下降4.1%）。而理论推导的首选平滑机制模型是幂律衰减模型，其隐含的凹有效成本函数为C(h)∝ln(1+h/h₀)，其中h₀=4.01 m，α=0.539，证实了城市密度随海拔衰减是强次指数形式。

齐普夫位频检验
通过CSN检验，最优下界截断x_min=429.28 m，估计指数α_Z=2.36，Kolmogorov-Smirnov（KS, Kolmogorov-Smirnov）距离 = 0.093，p值<0.001，表明强拒绝齐普夫位频标度。城市海拔分布虽是重尾分布，但并不具有齐普夫意义上的标度无关性，证实了海拔强度的幂律衰减与齐普夫位频标度是不同现象。

稳健性与诊断
通过敏感性测试和国家分层发现，幂律衰减模型对数据扰动和分层具有鲁棒性。区域分解揭示了热带地区海拔惩罚最弱，欧洲和北美衰减较快，内陆国家比沿海国家衰减更浅。

交通层级与海上准入
通过对比不同交通模式的燃料消耗，揭示了交通系统具有强烈规模经济和跨模式成本差异。海上运输成本最低，失去海上准入会导致集中在海平面附近的巨大离散成本惩罚，而在高海拔则是连续累积的小额增量，这一机制解释了恢复出的凹成本形式。

讨论与结论
研究人员总结出三点主要发现：第一，控制高程学土地可用性后，城市密度以强非线性方式随海拔下降；第二，凹成本模型明显优于线性成本指数模型，首选平滑表示为幂律衰减；第三，针对城市海拔拒绝了齐普夫位频标度，表明城市垂直分布由不同于城市规模水平分布的统计结构主导。研究表明，在最大熵框架下，有效成本函数的边际惩罚在近海平面最大，并随海拔升高而递减，即C(h) = -β^-1 ln f(h)所反映的对数成本设定提供了一种简约的凹成本设定。这一发现为聚落地理学、交通经济学和最大熵推理之间的联系提供了一种紧凑的方式。