这里是科研人的「方法急救站」！🚑

• • 欢迎加入我们的**"学术社群"**，每日更新热点论文复现指南！

我们通过小程序科研零时差追踪到: Nature Communications 近期发表题为“A universal scaling law of intra-urban inequality”的文章。第一单位为南京农业大学。

doi: 10.1038/s41467-026-73015-1
数据(代码)链接: https://doi.org/10.6084/m9.figshare.30122284
作者邮箱：xunan2025@szu.edu.cn, owx@njau.edu.cn

标签：#城市环境公平 #气候变化 #可持续发展 #尺度分析 #空间分析 #最小二乘法回归

本文内容速览：

1. 1. 提出科学问题
2. 2. 文章的主要结论
3. 3. 分析过程和方法
4. 4. 研究的局限性

1. 提出科学问题

1.1 研究领域现状

城市是人类应对全球环境变化影响的核心区域。快速的城市化进程改变了地表形态，带来了经济的集聚效益，同时也引发了显著的社会和环境成本。当前的城市科学研究中，城市标度律（urban scaling theory）已被广泛应用于解释城市宏观属性随人口规模的变化规律，例如财富和创新等社会经济产出呈现超线性增长，而基础设施呈现亚线性增长。与此同时，环境公平研究记录了城市内部资源和风险分配的不均现象，例如低收入群体面临更高的热暴露风险和更少的绿地资源。然而，现有标度律研究通常将城市视为一个整体的单一节点，忽略了城市内部的空间异质性，而环境公平研究则多依赖于特定城市的案例分析，缺乏全球尺度的统一理论框架。

1.2 本文要解决的关键科学问题

为了将宏观标度理论与微观空间不平等现象相连接，本文重点探讨了以下问题：

• • 问题 1: 城市内部在热暴露、绿地分布和经济活动方面的不平等程度，是否也遵循随城市人口规模变化的普遍标度律？
• • 问题 2: 全球不同国家和地区的社会经济发展水平以及背景气候条件，是如何调节这种空间不平等标度关系的？

1.3 研究的理论/现实意义

本研究通过建立覆盖全球一万多个城市的空间数据库，证实了城市内部不平等具备普遍的标度特征。这一发现拓展了传统城市标度理论的适用边界，将其从描述城市整体属性延伸至描述城市内部的资源分配结构。同时，该研究为环境公平领域提供了定量化、预测性的理论基础，证明实现公平的城市发展不能仅依赖规模效应，而需要具有针对性的政策干预，这对实现全球可持续发展目标（SDGs）具有直接的指导价值。

2. 文章的主要结论

本研究利用多源卫星遥感数据对全球11,000多个城市进行分析，揭示了城市规模与内部不平等之间的系统性关联：

• • 结论 1: 城市内部的热暴露、绿地资源和经济活动均与城市人口规模呈显著的超线性标度关系（super-linear scaling laws）。具体而言，城市人口规模每增加一倍，这三个维度的不平等基尼系数（Gini coefficient）将增加约 8% 至 9%。这意味着大城市在本质上倾向于拥有更不平等的环境与经济资源分布。
• • 结论 2: 标度关系受到国家社会经济状态和自然气候条件的显著调节。在低收入国家，绿地和经济不平等随着城市规模扩张而加剧的速度远高于高收入国家；而在干旱和半干旱气候区，热暴露不平等的加剧效应最为强烈。

3. 分析过程和方法

更加详细的复现指南可加入学术社区获取!

为了系统论证城市内部不平等的标度律，作者构建了一套严谨的大尺度空间数据分析工作流。该分析流程主要包括数据获取与预处理、空间不平等指标的计算，以及标度律建模与调节效应分析三个核心阶段。

首先，在数据获取与预处理阶段，作者基于全球人类聚居地城市中心数据库（GHS-UCDB）确定了11,210个有效的高密度城市核心区样本。随后，作者利用 Google Earth Engine (GEE) 平台处理海量遥感数据，以提取三个维度的环境与经济暴露指标。
在绿地分布方面，采用 Sentinel-2 表面反射率图像计算生长期内的归一化植被指数（NDVI）。作者通过比对高分辨率影像，确定了 NDVI > 0.39 为提取有效城市功能绿地的阈值。
在热暴露方面，利用 Landsat 8/9 数据提取夏季地表温度（LST）。该步骤应用了统计单窗（SMW）算法，动态结合裸土发射率和植被覆盖度，并进行大气校正，以获取精确的地表温度。
在经济活动方面，选用可见光红外成像辐射仪（VIIRS）的夜间灯光（NTL）数据作为微观经济活动强度的代理变量。同时，结合 WorldPop 的100米分辨率网格人口数据，为后续的加权计算提供基础。

进入空间不平等指标的计算阶段，这是本文方法论的关键所在。作者并未简单计算整个城市的统计方差，而是计算了基于邻里暴露水平的人口加权基尼系数。
针对每一个100米分辨率的人口网格单元，作者执行了焦点统计分析（Focal Analysis），计算该单元周围500米半径缓冲区内的平均绿地覆盖率、地表温度和夜间灯光强度。这种处理方式能够真实反映特定网格内居民在日常活动范围内的环境暴露情况。随后，将每个网格的邻里平均暴露值与其人口数量配对。在 R 语言环境中，调用 dineq 包中的 gini.wtd 函数，以暴露值为向量、人口数量为权重，计算出每个城市在三个维度上的基尼系数。基尼系数介于0（绝对平等）和1（绝对不平等）之间，用以表征城市内部的资源分布离散程度。

在获得所有城市的基尼系数后，作者进行了标度律建模。由于基尼系数属于比率型强度指标（intensive metric），其基准预期是不随城市规模变化（即标度指数 ）。作者对基尼系数和城市总人口进行了对数转换，构建了普通最小二乘法（OLS）回归模型：。结果直观展示在对数散点图中，所有维度的标度指数 均显著大于0（介于 0.11 至 0.13 之间），统计学上证实了超线性关系的存在。

为了进一步探究不同背景条件对上述标度律的影响，作者开展了调节效应分析。作者依据世界银行的标准，将样本划分为四个国家收入水平，并依据干旱指数（Aridity Index）将其划分为五个主要气候类型。

在统计实现上，作者在原有的回归模型中引入了对数化人口与分类调节变量（如收入水平或气候类型）的交互项（Interaction Term）。通过交互项的显著性检验，作者能够判断不同组别之间的回归斜率（即标度指数 ）是否存在统计学差异。图表展示中，通过绘制不同颜色和线型的回归拟合线，能够清晰比对不同收入群体或气候区域下不平等随规模扩张的加剧速度。

此外，为了直观展示调节变量的定量差异，作者提取了各个分层回归模型中的标度指数 及其 95% 置信区间。通过带有误差棒的散点图形式呈现，若不同组别的置信区间不重叠，则直接证明了其统计差异的显著性。例如，低收入国家的绿地和经济不平等指数明显高于高收入国家；干旱区域的热不平等指数显著高于湿润区域。

最后，作者还进行了充分的敏感性和机制分析。考虑到城市形态可能造成的影响，除了使用人口密度作为代理变量，作者还计算了紧凑度指数（Compactness Index）、分形维数（Fractal Dimension）和伸长率等几何指标，并通过残差分析（Residual Analysis）排除了这些形态变量对标度律核心结论的干扰，确立了模型结果的稳健性。

4. 研究的局限性

尽管该研究提供了全球尺度下的定量依据，但仍存在部分局限性。首先，该标度律分析基于横截面数据（cross-sectional data），能够有效表征全球城市系统的宏观结构状态，但不应直接等同于单一城市的纵向历史发展轨迹，个体城市的演变可能因路径依赖而出现偏差。其次，研究采用人口密度作为城市形态的主要代理变量，未能完全捕捉城市三维形态的复杂性（例如街道峡谷对热量截留的微观影响）。第三，研究使用的遥感代理指标具有固有的局限性，例如地表温度（LST）并不能完全等同于人类真实的体感热舒适度，而夜间灯光也仅能反映经济活动的某一侧面。最后，基尼系数虽然量化了总体的不平等程度，但无法识别受不平等影响最严重的具体社会群体（例如按照种族或少数族裔划分），且当前的相关性分析尚未能充分揭示微观层面精确的因果机制。

如果你觉得这篇文章对你有帮助，欢迎点赞👍、收藏⭐️和分享🔗给更多的科研小伙伴们！

如果你有任何问题，欢迎加入我们的学术社群与大家讨论交流💬