4.1北京老年人公园绿地可达性
基于ArcGIS10.2的网络分析模块，使用两步移动法根据公式(1)和(2)计算得到了2010—2019年各街道(乡镇)的老年人公园绿地可达性。依据自然间断点分级法(Jenks)将计算得到的2010与2019年各街道(乡镇)的老年人公园绿地可达性分为5类，并描绘其空间格局(图4-1、4-2)。
由图4-1、4-2可知，北京各街道的老年人公园绿地可达性表现出明显的空间聚集特征，高可达性街道都集中于中心城区，郊区则存在许多公园可达性为0的街道。相较于2010年，2019年老年人公园绿地可达性有所提升，可达性为0的街道减少了16个。
北京公园绿地社会公平性的变动
根据公式(3)计算出了2010—2019年北京市整体和各区的老年人公园绿地可达性基尼系数，并据此对北京市整体和各区的公园绿地社会公平性变动进行了分析。
4.2.1北京市整体的老年人公园绿地社会公平性及其变动
从图5可以看出，2010—2019年北京市老年人公园绿地可达性的基尼系数总体呈下降趋势，从2010年的0.857下降到了2019年的0.838，下降幅度达到了2.22%。这表明，2010—2019年北京公园绿地社会公平性在总体上有所改善。具体而言，2010—2013年基尼系数呈上升趋势，2014—2019年基尼系数呈下降趋势。2019年，北京的老年人公园绿地可达性基尼系数等于0.857，说明现阶段从北京市整体来看公园绿地社会公平性仍然不理想，有很大的提升空间。
4.2.2北京市各区的老年人公园绿地社会公平性及其变动
图6是北京市各区的老年人公园绿地可达性基尼系数在2010—2019年的均值。由图6可知，北京市各区的公园绿地社会公平性在空间上呈现出以下特征。1)明显的圈层特征。位于最中心的首都功能核心区(东城、西城)公平性最好，首都功能核心区外围的城市功能拓展区(朝阳、丰台、石景山、海淀)公平性较好，而城市功能拓展区外围的城市发展新区和生态涵养发展区公平性则较差。2)区域差异显著。基尼系数最高的房山区是基尼系数最低的东城区的近3倍。
图7是北京市各区2019年的老年人公园绿地可达性基尼系数相较于2010年的变化率。
由图7可知，北京市各区的公园绿地社会公平性变动在空间上呈现出以下特征。1)公平性改善最为明显的区域是城市功能拓展区(朝阳、丰台、石景山、海淀)和北京城市副中心(通州)。2)公平性变差最为明显3个区分别位于北京的西北、东北和东南部，呈离散分布。
由图7可知，北京市各区的公园绿地社会公平性变动在空间上呈现出以下特征。1)公平性改善最为明显的区域是城市功能拓展区(朝阳、丰台、石景山、海淀)和北京城市副中心(通州)。2)公平性变差最为明显3个区分别位于北京的西北、东北和东南部，呈离散分布。
4.3北京老年人公园绿地社会公平性变动的驱动因素
本文通过回归分析来研究社会经济因素对公园绿地社会公平性的影响。为在研究各因素对公园绿地社会公平性影响的同时比较多个因素的相对影响差异，同时给出标准化和非标准化2种回归结果(表3)。
回归模型的Wald检验统计量等于14210.12，P值接近于0，说明模型通过了整体显著性检验。
由非标准化的边际效应可知：非首都功能疏解政策的边际效应为-0.0190，在5%水平上显著。这表明，非首都功能疏解政策的实施对公园绿地社会公平性改善有正向驱动作用。第三产业产值与第二产业产值之比每上升1%，基尼系数就下降0.0000442，这说明产业结构优化升级是公园绿地社会公平性改善的正向驱动因素。
公园绿地覆盖率每增加1%，基尼系数就下降-0.0211，说明增加公园绿地的供给可以改善公园绿地社会公平性。但以老年人口密度衡量的公园绿地需求对基尼系数有正向影响，老年人口密度每增加1%，基尼系数就增加0.180，说明公园绿地需求的增加对公园绿地的社会公平性有负面影响。此外，人均GDP和一般公共预算收入占地区生产总值的比重对公园绿地社会公平性并无显著影响。
公园绿地覆盖率每增加1%，基尼系数就下降-0.0211，说明增加公园绿地的供给可以改善公园绿地社会公平性。但以老年人口密度衡量的公园绿地需求对基尼系数有正向影响，老年人口密度每增加1%，基尼系数就增加0.180，说明公园绿地需求的增加对公园绿地的社会公平性有负面影响。此外，人均GDP和一般公共预算收入占地区生产总值的比重对公园绿地社会公平性并无显著影响。
由标准化边际效应可知：各社会经济因素对公园绿地社会公平性的影响由大到小依次为：公园绿地需求、公园绿地供给、产业结构和非首都功能疏解政策。
