基于1992-2013年中国城市遥感模拟反演碳排放数据,采用空间自相关、空间马尔科夫矩阵和动态空间面板数据模型,在同时考虑碳排放的时空滞后效应和不同地理经济空间权重矩阵的条件下,对城市碳排放的演化路径和关键影响因素进行了定量识别和减排政策探讨。研究表明,中国城市能源消费碳排放的区域差异正逐步缩小,空间上呈现出明显的高排放俱乐部集聚特征,同时碳排放类型演化具有明显的路径依赖特征;面板数据模型估计结果表明经济增长与人均碳排放呈现显著的倒“U”型曲线关系,而绝大多数城市的人均碳排放处于随经济发展而增加的阶段,二产偏重的经济结构和投资的粗放增长共同正向作用于城市碳排放,而人口的集聚效应、技术水平的提升、对外开放度和公路运输强度的增加则共同抑制城市碳排放水平的提高。因此未来要抑制促增因素和发挥促降因素的作用才能有效降低城市碳排放;优化产业结构、精简粗放投资、增加研发强度以及提升公路通达性是未来实现中国城市节能减排的有效途径。

自2000年提出生态功能保护区的概念以来,历时17年不断探索和发展完善,中国已经形成了一批国家重点生态功能区和相关配套政策,成为生态保护红线划定的重要依据和主体功能区战略实施的主要组成,对保障国家生态安全和优化国土空间格局具有十分重要的意义。海南岛中部山区在2001年成为国家生态功能保护区试点,2010年被确定为国家重点生态功能区,是较早实施生态功能保护的重点管控区域之一。本文以海南岛中部山区热带雨林重点生态功能区为研究区域,综合利用多源卫星遥感数据和地面观测统计数据,集成地统计学和生态评估模型模拟方法,以生态系统格局、植被长势、服务功能为主要评估内容,通过时空尺度拓展和多维时空对比分析,探讨了该重点生态功能区生态系统变化和保护成效。结果发现：① 2013年功能区森林面积占比为84.5%,显著高于海南岛平均水平。1990-2013年功能区内建设用地逐渐增加。2010年功能区实施后,人类活动强度仍然呈现持续增加态势,面积比例从0.5%增加到0.8%,主要表现为城镇建设和耕地开垦。② 功能区水源涵养服务功能明显优于功能区外部,1990-2013年水源涵养服务功能总体上略有增加,但是不同时间段的波动变化明显,特别是1990-2000年和2000-2013年两个时间段内均表现为先减少后增加的变化特征。功能区单位面积的水源涵养量为51.78万m²/km²,明显高于功能区外部。③ 功能区土壤保持服务功能明显要优于功能区外部,1990-2013年表现出剧烈波动、稳定性较差。功能区单位面积的土壤保持量为1.95万t/km²,明显高于功能区外部。④ 功能区人类扰动指数明显小于功能区外部,生物多样性威胁程度明显偏小,有利于生物多样性的保护。功能区平均人类扰动指数为0.3664,比功能区外部低0.1152。1990-2013年功能区外部人类扰动指数变化幅度为0.0152,约为功能区内部变化幅度的5.31倍。特别是2010-2013年,功能区内部人类扰动指数增加幅度要明显小于功能区外部。

因区域间贸易而转移的碳排放对于全球各个国家或者地区的碳排放核算及其减排责任划分具有重要而深远的影响。通过构建多区域投入产出分析模型和基于扩展的STIRPAT模型,本文核算了全球39个主要国家基于消费责任制的碳排放,并在此基础上深入探讨了全球环境压力的影响因素问题。研究发现：首先,从以净流出为主的中国和俄罗斯来看,尽管贸易给这些地区带来了大量的资源,但随之会产生严重的区域生态环境问题。而对于美国和欧盟等地区来说,这些国家或者地区通过全球贸易规避了大量的碳减排责任。此外,在全球贸易中,区域净流出的贸易隐含碳排放越小,基于消费责任制核算的区域碳排放量就越大,故而在全球减排目标分配中承担的减排任务也需相应地增加。全球贸易隐含碳排放净流出量较大的地区主要位于亚洲和东欧等地,而净流出量最小的地区主要是以西欧和北美等高度发达的经济体为主。其次,就基于消费责任制核算下的全球环境压力而言,人口因素和富裕程度是导致全球环境压力不断增加的两个重要因素,而提高生产技术水平以及逐渐提高清洁能源在总能源消费中的比重是推动全球节能减排、缓解全球环境压力的有效途径。

为研究北京市甲型H1N1流感对气象因子的时空响应规律,本文结合地理探测器及空间自相关分析方法,分析北京市甲型H1N1流感的时空分布情况,基于地理加权回归（GWR）模型,建立北京市甲型H1N1流感与平均风速、降雨量、平均湿度以及平均温度之间的回归模型。结果表明,降雨量对流感传播的影响不显著,平均风速与流感的传播主要呈现正相关,平均相对湿度与平均温度与流感的传播主要呈现负相关,其中,平均相对湿度的影响更为复杂。时间上,在流感感染人数快速上升的阶段,各气象因子的影响程度较为显著,而在流感传播速度减慢的阶段,气象因子的影响减弱;空间上,各气象要素对流感传播的影响在北京市4大功能分区上存在空间分异性,这与地区地形、风带、局地气候等因素有关。

运用ArcGIS软件建立了长江中游地区9.0~4.0 ka BP文化遗址的空间数据库,采用叠加分析和统计学方法,对文化遗址的规模、海拔分布、坡度、距河流远近等特征进行了系统分析。结果表明：① 由新石器时代早期到晚期,遗址的数量、规模和分布范围等均呈现增长或扩大趋势,叠置型遗址个数增加;② 遗址一般具有低海拔、平缓地形和距河流较近的选址倾向,晚期的文化遗址大体保持该倾向,但对海拔、坡度及距河流远近等的选择开始出现分化。③ 新石器时代晚期不同规模的文化遗址在选址上存在一定分异。小型遗址多分布在坡度为6°~15°的缓坡地,中型及以上规模的遗址更倾向于坡度为2°~6°平缓地;中小型至大型遗址的海拔分布均以30~50 m为主,特大型遗址的海拔分布以50~100 m为主;特大型遗址在距河流1000 m以内的分布比例高于其他规模遗址;④ 新石器文化在早中期受气候、水源、地形等自然环境因素影响较大,晚期文化的适应能力有所增强。

晚新生代以来青藏高原的强烈隆升造成青藏高原东南边界形成大量的构造断裂,同时构造拉伸和挤压作用也塑造了许多构造断陷湖泊。近年来很多湖泊的面积和深度逐渐减小,对当地生态环境和社会经济发展带来极大影响。选取青藏高原东南缘的西昌邛海流域为研究区,探讨邛海北岸官坝河泥石流发育特征及对邛海的泥沙淤积效应。研究发现：① 官坝河泥石流暴发频率高且规模大,崩塌滑坡和面源侵蚀是泥石流最主要的固体物质来源,动储量约为428.03万m³;② 山洪泥石流的频繁发生是邛海泥沙淤积的主要原因,而频繁地震和强降雨耦合作用控制着区域山洪泥石流活动,预测邛海流域泥石流活动仍将频繁发生,泥沙淤积问题更加严重;③ 自1950年以来,邛海流域泥沙淤积速率约为17.09 mm/a,如按当前淤积速率推算,邛海的寿命将仅约为600年左右;④ 建议对官坝河流域实施工程治理和生物防护相结合的综合治理,并加强流域管理和宣传培训等行政法治管理。该成果可为邛海流域综合规划管理和防灾减灾提供参考依据,有助于揭示青藏高原东南缘湖泊逐渐萎缩的真正原因和控制要素。

中国已建立各类自然保护区2740个,占全国陆地面积的14.8%。本文以典型国家级自然保护区为例,基于遥感反演、模型模拟、空间分析等方法,获得植被覆盖度、净初级生产力（NPP）和土地覆被时空数据集,分析2000-2015年中国不同区域、不同类型自然保护区生境状况时空变化及人类扰动的影响,进而评估自然保护区在栖息地及生物多样性保护方面的效果。结果表明：2000-2015年,国家级自然保护区植被覆盖度从36.3%提高到37.1%,各类型自然保护区均有不同程度提高,其中森林类保护区年增速0.11%,草原草甸类0.84%,内陆湿地类0.21%,荒漠生态类0.09%,野生动物类和野生植物类则分别为0.11%和0.08%。草原草甸类、内陆湿地类、荒漠生态类、野生动物类自然保护区植被NPP年增速分别为2.06 g·m^-2、1.23 g·m^-2、0.28 g·m^-2、0.4 g·m^-2,而森林类和野生植物类则分别以3.45 g·m^-2和2.35 g·m^-2的年速率减少。近15年,国家级自然保护区内人类扰动呈现微弱变化,除青藏高原区和南亚热带湿润区内保护区人类扰动略微下降以外,其他区域均有所增强,特别是北亚热带和温带湿润区,其保护区人类扰动由4.70%明显增至5.35%。

为了揭示平原水网地区快速城市化对河流水系的影响,基于1:50000地形图和30 m空间分辨率遥感影像数据,采用格网化河网密度、支流发育系数和分形维数等水系指标,分析了20世纪80年代至21世纪10年代苏州市河流水系的时空变化特征,探讨了土地城市化水平与河网密度变化之间的定量关系。结果表明：① 苏州市河流长度衰减了10.55%,其中干流长度增加了3.24%,一级支流与二级支流的长度分别减少了18.87%和11.76%;河网密度在空间上以减小为主,其中干流河网密度在空间上以增加为主,支流河网密度在空间上以减少为主。② 支流发育系数由3.36下降至2.78,水系结构趋于主干化;水系分形维数由1.70下降至1.63,水系形态趋于简单化。③ 土地城市化率高于40%时,河流水系衰减趋于明显,且土地城市化水平越高,河流减少的越多。为满足不同的经济社会发展需要,不同等级河流对土地城市化的响应不尽相同：为提高城市防洪排涝能力,一级河流总体上有所增加;土地城市化率小于40%时,部分二级河流被侵占以增加耕地面积,并新开挖了许多三级河流以提高农田排涝与灌溉能力;土地城市化率大于40%时,大量支流被填埋并转换为建设用地,而且等级越低的河流衰减的越严重。

河湖水系连通作为水资源调配、水生态修复与改善的重要手段,已成为江河湖泊治理工作的重点之一。河湖水系连通方案的选择和连通效应的评估直接关系到连通工程对湖泊水动力和水环境的改善效果,但目前缺乏对水系连通后河湖水动力和水质的综合评估体系。以汤逊湖湖泊群为研究对象,构建二维水动力水质数学模型,选取NH₃-N、TN和TP作为水质指标,模拟不同连通方案下汤逊湖湖泊群的流场及水质变化情况,结合水动力、水质和经济社会三个方面的评价指标建立综合评价体系,对各连通方案下湖泊群的改善效果进行评估。研究结果表明：引水后湖泊水体流动性有了很大的改善,河湖连通工程能够在短时间内达到改善湖泊水质的目的,但随着引水时间的增长,改善率逐渐减小。通过对比五种连通方案发现,方案五对水动力和水质的改善都取得了最显著的效果,且能够带来最大的引水效益。

高宝诸湖位于江苏省中部,紧邻大运河西岸,是淮河重要的入江水道。近300年来受到黄淮关系演变和人类垦殖活动的影响,高宝诸湖的形态发生了较大变化。在对古今地图中高宝诸湖的形态进行数字化和空间比较的基础上,结合对历史文献和部分钻孔资料的综合分析,将高宝诸湖的演变分为三个阶段：① 清康熙至乾隆年间为近300年来水域面积最大的时期,1717年水域面积为1606.02 km²;② 清中期以后受淮河通过高家堰不定期排泄泥沙的影响,湖中沙洲开始发育;③ 1851年开始,淮河改道经由洪泽湖东南的三河闸,通过高宝诸湖进入长江,其携带的部分泥沙在高宝诸湖沉淀下来,形成三河三角洲,随着三河三角洲在在高宝诸湖内的不断淤涨,湖面萎缩速度加快,1868年湖面缩小到1494.03 km²,1916年为1071.97 km²,2011年为886.90 km²。洪泽湖泥沙是导致近300年来高宝诸湖不断萎缩的主要原因,此外,圩田的兴建进一步加剧了这一过程。

青弋江位于黄山北麓,为长江下游最长的支流。野外考察发现青弋江泾县盆地段存在溪口剖面和城北剖面等2个天然剖面,共发育了1级洪积扇台地（P）和3级河流阶地（T3、T2和T1）,并相应堆积了4级砾石层。通过对砾石层进行砾组分析,并借助电子自旋共振（ESR）测年和古地磁测年等方法,初步探讨了青弋江发育的年代、过程和成因。研究结论为：① 砾组分析表明青弋江T3阶地是青弋江的最老阶地,并且其砾石层是青弋江的最老砾石层;② 测年结果表明青弋江发育的年代区间为1300~900 ka,其中~1300 ka为青弋江发育的最早年代,而~900 ka则为青弋江发育的最晚年代;③ 青弋江发育于~1377 ka前的洪积扇辫状河,并先后经历了洪积扇及辫状河发育、辫状河下切、青弋江形成等阶段,即所谓的源于洪积扇辫状河的青弋江发育过程;④ 青弋江发育可能是降水增加和构造抬升共同作用的结果。该研究有助于为中国东部地区中小河流发育研究提供参考。

中国人工林面积居世界首位,在全球气候变暖背景下,人工林在提供木材和固定CO₂等方面同时具有重要意义。然而木材生产与碳储量之间存在此消彼长的权衡关系,如何综合考虑二者权衡过程,实现区域森林综合效益最大化,是中国森林生态建设中亟待解决的问题。选择会同森林生态实验站磨哨林场作为研究区,采用人工林固定样地多年观测数据以及InVEST模型,定量评估木材生产与碳储量,构建了不同管理模式下的权衡分析方法,提出了研究区森林生态系统综合效益最佳的权衡对策。研究发现：① 在采伐活动的影响下,木材生产与碳储量显著负相关（R = -0.907, P < 0.001）,表现出强权衡关系;② 随着采伐强度增加,木材生产与碳储量的综合效益逐渐增加,权衡值先下降后上升;每10年采伐总面积10%~20%的管理模式下,收获木材量和碳储量的综合效益相对较高且权衡值最低,是最佳森林管理模式;③ 研究区当前采伐强度略高于最佳采伐强度,适当降低采伐强度能够使森林获取更多的综合效益。因此,在明确森林管理目标的基础上,可借助基于标准差的权衡分析方法,提出中国南方人工林生态系统服务提升与可持续经营的权衡对策。

如何充分利用离散的观测数据,通过对维管植物物种分布丰富度及其与生境因子之间的相互作用和影响机理的定量分析,实现维管植物物种丰富度的空间分布及其情景模拟,是目前生物多样性研究前沿和核心内容之一。针对这一问题,在实现青藏高原37个国家自然保护区的维管植物物种数量收集和边界数据矢量化的基础上,分别进行维管植物物种数量与土地覆盖类型、环境因子和景观生态指数等三大类生境因子之间的相关关系的定量计算和对比分析,筛选和确定最佳相关分析方程,进而构建青藏高原维管植物物种丰富度的空间模拟分析模型。该模型中,维管植物物种丰富度与生境因子之间的复相关系数为0.94,模型验证结果表明,青藏高原的维管植物物种的平均丰富度为496.79种/100 km²,其空间分布格局整体上呈东南向西北逐渐减少趋势;另外,除柴达木盆地荒漠区域以外,维管植物物种的空间分布随海拔的升高而减少。基于CMIP5 RCP 2.6、RCP 4.5和RCP 8.5三种气候情景模拟获得的青藏高原维管植物物种丰富度未来情景结果显示,在T0-T4（2010-2100）时段内,青藏高原维管植物物种丰富度整体将呈减少趋势。RCP 8.5情景下青藏高原维管植物物种丰富度的变化幅度最大,而RCP 2.6情景下的维管植物物种丰富度的变化幅度最小。研究表明,本文构建的模型能够对青藏高原维管植物物种丰富度的空间分布格局及其未来情景进行模拟分析,模拟结果可为青藏高原生物多样性及其对气候变化响应的综合评估和情景模拟提供方法和技术支持。

京津冀城市群是中国雾霾最严重的区域,在京津冀协同发展背景下,探究该地区大气污染的时空分布和影响因素具有重要意义。运用空间自相关分析和三种空间计量模型,分析了京津冀202个区县PM_2.5的时空分异特征,创新性地对自然与人文影响因素贡献及其空间溢出效应进行系统地甄别和量化。结果表明：2000-2014年来京津冀城市群PM_2.5浓度整体呈上升趋势,季节上呈秋冬高、春夏低,空间上呈东南高、西北低的特点,且城市建成区PM_2.5浓度比周围郊区和农村平均高10~20 μg/m³;2014年仅有13.9%的区县空气质量达标,PM_2.5浓度存在显著的空间集聚性与扩散性,城市间交互影响距离平均为200 km,邻近地区的PM_2.5每升高1%,将导致本地PM_2.5至少升高0.5%;社会经济内因对PM_2.5主要是正向影响,自然外因主要是负向影响;影响因素中对本地大气污染的直接效应贡献强度依次是：年均风速>年均气温>人口密度>地形起伏度>第二产业占比>能源消费>植被覆盖度,人均GDP、年降水量和相对湿度对本地PM_2.5没有显著影响;对邻近地区大气污染具有显著空间溢出效应的因素排序是：植被覆盖度>地形起伏度>能源消费>人口密度;对于自然和人文影响因素应分别采取针对性的适应策略和调控策略,加强区域间联防联控与合作治理,在城市群规划中注重环保规划与立法。

高浓度PM_2.5是形成雾霾的主要原因之一,科学识别PM_2.5浓度的空间异质性与驱动因素对区域大气联动治理意义重大。研究采用2000-2011年NASA大气遥感影像反演的PM_2.5数据,结合地统计、地理探测器及GIS空间分析等方法,系统分析了中国2000-2011年PM_2.5浓度时空演化格局特征与其驱动因素。结果表明：① 2000-2011年中国PM_2.5污染平均浓度一直保持在22.47~28.26 μg/m³区间,总体呈现先快速增加后趋于稳定的演化态势,2006年是PM_2.5浓度值变化的拐点（峰值）。② 空间上PM_2.5浓度整体呈现北方高于南方,东部高于西部趋势,污染浓度高值区集中分布在黄淮海平原、长三角下游平原、四川盆地与塔克拉玛干沙漠四大区域,其中京津冀地区污染最为严重。③ 污染浓度重心研究表明PM_2.5重心总体呈现快速东移趋势,污染高值区持续向东移动,低值区向西移动,两者重心背向而行,表明东部雾霾污染程度在进一步加剧。④ 空间自相关分析表明PM_2.5年均浓度呈现强烈的局部空间正自相关特性,PM_2.5“高—高”集聚区连片分布在黄淮海平原、汾渭盆地、四川盆地及江汉平原地区,PM_2.5“低—低”集聚区分布在长城以北的内蒙古、黑龙江、青藏高原、以及台湾、海南与福建等东南沿海及岛屿地区。⑤ 地理探测分析表明气候等自然因素与人类活动共同对PM_2.5浓度空间变化产生巨大影响,其中自然地理区位、人口密度、汽车数量、工业烟尘、秸秆燃烧等因子是中国PM_2.5浓度空间变化的主要驱动因素。

多种生态系统服务之间的相互关系研究是生态系统服务研究的热点问题。目前,有关多种生态系统服务权衡与协同的研究多停留在定性描述阶段,缺乏具备空间位置信息的定量研究。综合土地利用、NDVI、土壤类型、气象等多源数据,对汉江上游流域2000-2013年的土壤保持服务、产水服务、植被碳固定（NPP）服务进行空间制图,并基于逐像元偏相关的时空统计制图方法,对三种生态系统服务之间的权衡与协同关系时空变化进行分析。结果表明：① 2000-2013年,汉江上游多年平均土壤保持量为434.20 t·hm^-2·yr^-1,且多年来以16.10 t·hm^-2·yr^-1速率递增。多年平均产水量为250.34 mm,多年递增速率为3.79 mm·yr^-1。NPP多年平均值为854.11 gC·m^-2·yr^-1,多年递减速率为8.54 gC·m^-2·yr^-1。② 汉江上游土壤保持量在空间上呈现出“南北山地高,河谷盆地低”的分布格局;植被NPP在空间上表现为汉江河谷地带低,其余地方高的态势;产水服务在空间上呈现由北向南递增的分布规律。③ 整个汉江上游流域土壤保持与NPP、土壤保持与产水服务均以权衡关系为主,权衡关系区域的面积占比分别达到62.77%和71.60%;NPP与产水服务以协同关系为主,协同关系的面积占比达到62.89%。④ 林地、湿地、耕地、人工表面及裸地中,土壤保持与NPP服务、土壤保持与产水服务在空间上以权衡关系为主,产水与NPP服务在空间上以协同关系为主;而草地表现为三种服务两两之间均以权衡关系为主。定量评估生态系统服务相互关系的时空特征,有利于研究生态系统服务之间相互关系在时间上的非线性变化以及空间上的分布异质性,对区域土地管理与生态系统服务优化具有重要的指导意义。

东部沿海地区城市群是国家经济发展的战略核心区和国家新型城镇化的主体区之一,其经济—生态可持续发展问题是近年来地理研究领域的前沿主题。首先对城市群生态效率进行定义并将其作为城市群可持续发展测度指标;随后运用传统DEA模型及非期望产出SBM模型,对东部沿海地区四大城市群在2005年、2011年、2014年三个时间截面中的经济效率、生态效率进行比较评价并分析了城市群生态效率时空演变特征;最后通过松弛度分析,为改善沿海地区四大城市群生态效率提供了对应建议。结果表明：① 山东半岛城市群、长三角城市群、珠三角城市群经济综合效率呈先下降后上升的V字型演变格局,而京津冀城市群经济综合效率始终呈下降趋势,2011年前其受纯技术效率影响较大,2011年后受规模效率影响较大;② 京津冀城市群受污染影响带来的效率损失最大,山东半岛城市群相对较小,总体而言,四大城市群生态效率在2005-2011年为下行期,2011-2014为修复期;③ 沿海四大城市群内城市生态效率时空格局具有不同的演化特征,核心城市周边以及沿海沿江地区城市生态效率相对较优,内陆城市则较低;④ 京津冀城市群、长三角城市群、珠三角城市群核心城市资源消耗强度高、经济效益产出高、生态效率高,从整体来看,城市群内多数城市呈污染物排放量下降趋势,污染物治理效果较为明显,环境影响问题有所改善;⑤ 影响东部四大城市群经济生态效率的原因不同,本文从城市群角度进行了差异性分析,并提出了合理化建议。

淤涨型沿海滩涂作为中国重要湿地和耕地后备资源,承担着生态保护、社会保障、经济发展等多重服务功能。选取江苏中部淤涨型滨海滩涂开展实证分析,分析围垦年限为0~40年滩涂垦区耕地质量因子空间差异,采用“时空代换”方法揭示滨海滩涂围垦后耕地质量的演变趋势,为沿海滩涂围垦区耕地资源可持续利用提供决策参考。结果表明：江苏中部沿海滩涂土壤总体呈现盐碱度高、养分贫瘠等特征且土壤理化性质变异程度高,土壤质地呈现显著南北和陆海梯度特征。滩涂垦区围垦初期地形、植被状态和围垦后滨海滩涂转变为耕地路径对耕地质量演变趋势均具有重要影响。沿海滩涂转为耕地路径由垦区自然条件、开发主体、开发规模等因素共同作用,江苏中部沿海滩涂转为耕地路径已从“盐生植被→养殖水体→耕地”转变为“盐生植被→耕地”,该转变缩短了滩涂转为耕地的周期,但同时也带来围垦滩涂撂荒闲置现象。围垦后沿海滩涂土壤逐步向耕地土壤转变,耕地质量因子演变趋势可分为四种：上升型,如土壤全磷、有效磷、粘粒含量;下降型,如土壤砂粒含量;下降—上升型,如有机质、全氮、碱解氮和CEC;上升—下降型,如土壤盐分和pH。围垦初期,滩涂土壤脱离海水环境,土壤有机质发生矿化过程呈现下降趋势,土壤盐分和pH受到蒸发和蒸腾作用急剧上升;围垦中期,受到人类耕作作物秸秆等有机物输入以及施肥的影响,土壤有机质、养分等呈现上升趋势,土壤盐分和pH呈现下降趋势;围垦后期,耕地质量因子受到耕作管理影响存在一定幅度的波动,由于沿海滩涂垦区耕地形成时间短,耕地质量对耕作管理等因素敏感度高。

目前有关北京市不透水地表的相关研究多数是从整体层面开展,忽略了其内部功能区的作用及差异。运用分类回归树（CART）及系列变化检测模型得到北京市1991年、2001年、2011年和2015年四期不透水地表分布,并运用标准差椭圆、洛伦兹曲线、贡献指数及景观格局理论对各功能区不透水地表的时空变化进行分析。结果显示：1991-2015年北京市不透水地表的总面积增加了约144.18%,分布的主导方向由早期的东北—西南趋向于当前的正北—正南。各功能区间不透水地表的空间分布异质性逐渐减弱,但贡献指数值存在很大差异：功能拓展区的贡献指数最高,其四年中的最低值（1.79）高于其他功能区四年最高值,是北京市不透水地表增长最主要的贡献区;功能核心区的蔓延度指数值最高,约为其他功能区的2倍,为不透水地表的优势聚集区;发展新区的贡献值由负值变为正值并成倍增长,成为北京市不透水地表增长的主要源区;生态涵养发展区的贡献指数始终为负,并逐年减小。不同类型不透水地表的景观指数和质心偏移均存在差异,高盖度不透水地表的形状指数和斑块密度值最小,分布最为集中,对生态环境影响较大,北京市在未来发展过程中应合理规划控制其空间格局及增长模式,尽量减缓其增长速度及团聚程度。