随着全球气候变暖问题的日益严峻,人类社会高度重视全球气候变化问题
[1]。气候是人类生存环境的重要因素,也是生产生活的重要资源。然而当前全球气候变化却成为自然生态系统、世界经济可持续发展和人类社会面临的严峻挑战
[2-3]。相关研究
[4]表明,人类活动对气候系统有显著影响。其中,人类活动引起的土地利用变化是人类活动影响气候系统的主要途径之一
[5]。当前土地不合理利用带来的全球气候变暖问题日益加剧,给局部或更大范围内的气候造成影响,甚至给全球气候带来干扰。因此,土地利用变化对气候变化的影响研究在全球气候变化的研究领域占有重要地位,对人类社会可持续发展具有重要意义
[6-7]。随着中国社会经济的快速发展,资源需求量不断增大,导致资源环境开发强度加大,土地利用变化剧烈,深刻改变陆地表层的物理生物属性,通过地表反射率与地面粗糙度等地表通量变化影响大气环流,对区域气候产生重要影响。
早期有关土地利用变化对气候影响的研究主要集中在理论分析方面
[8-9],主要运用气象观测法
[10]和理论推导法
[11],研究重心集中在全球大尺度土地利用与气候变化的关系上
[12-13]。例如,FOLEY等
[14]的研究显示,森林减少导致碳循环和水循环的显著变化,而BONAN
[15]关于农业扩张的研究则展示了土地利用变化对气温和降水的长远影响。然而,全球尺度的研究受限于气候自然变率的掩盖,土地利用变化带来的局地效应常难以识别
[16]。随着21世纪区域气候模式的发展,数值模拟逐渐成为研究土地利用变化与气候相互作用的重要手段。区域气候模式不仅可进行长时间序列模拟,还可通过敏感性试验控制单一变量,获得土地利用变化的直接效应和反馈机制
[17],如RegCM和WRF能够有效模拟不同土地类型变化对气候的差异性影响
[18-19]。研究方法和数据分辨率的进步使得土地利用变化对气候影响的研究对象从全球扩展至区域和局地尺度,尤其是生态系统脆弱
[20]和气候变化敏感地区
[21]。有研究
[22]证实土地利用变化对气候影响显著,且不同土地类型变化对区域气候影响不同
[23-24]。然而目前对不同时期、不同区域和不同气候背景下土地利用变化对气候的影响仍缺乏系统深入的研究。
鉴于此,本研究运用最新的区域气候模式RegCM 4.7,结合数值模拟,聚焦中国大陆区域,模拟2009—2018年土地利用变化的气候效应,旨在揭示中国土地利用变化对区域气候的影响机制,并通过地表能量平衡和低层环流场的变化,探讨气温和降水的响应规律。本研究不仅可为中国区域土地利用变化气候效应研究提供案例,同时为构建以气候韧性为导向的国土空间规划提供依据,有助于制定应对气候变化的国土用途管理策略。
1 材料与方法
1.1 数据来源
研究数据包括:1)区域边界矢量图。来源于中华人民共和国自然资源部,审图号为GS(2020)4632。2)土地利用数据。来源于IGBP(国际地圈生物圈计划)全球GLCC土地覆盖数据库(edcdaac.usgs.gov),结合中国土地利用动态遥感监测数据集(CLUD)进行逐年更新与校正,生成2009—2018年土地利用序列。3)大气边界强迫数据。使用2009—2018年区域尺度初始和侧边界条件、海温(SST)等大气边界强迫数据来源于ERA-interim 气候再分析数据集。本研究选择处理后的EIN15(ERA-interim with the resolution of 1.5°)对模型进行初始条件设置,时间分辨率为6 h,空间分辨率为1.5°×1.5°,以此为基础合成逐月数据。4)气象观测数据。用于进行敏感性分析、与模型模拟结果进行对比的气象观测数据来源于英国东英格利亚大学气候研究所(CRU)高分辨率全球数据(
www.cru.uea.ac.uk),分辨率为0.5°。
1.2 RegCM 4.7区域气候模式
RegCM(regional climate model)是基于MM5发展的静力平衡区域气候模式
[25],也是本文采用的模型。该动力学模型是一个可压缩的有限差分模型,并带有流体静应力平衡和垂向σ坐标
[26]。垂直方向分为18层,地形追随坐标用于底层模拟。RegCM 4.7旨在提供对局部气候变化的高分辨率模拟,其物理过程是一个复杂而系统的体系,涉及大气、海洋和陆地表面等多个方面。本研究主要涉及辐射方案、陆面模式、行星边界层方案、积云对流方案、大尺度降水方案、海洋通量参数、海温预报等参数方案,默认情况下,根据内插的每周或每月SST(海面温度)产品,每6 h规定1次SST。
1.3 模型预处理
本研究中RegCM 4.7区域气候模式主要以土地利用数据为驱动模拟研究区气候变化情况。将模型的中心纬度设为30°21′N,中心经度设为104°21′E。输出数据分辨率设为50 km,将网格设为110×110个单元;垂直方向压力层数目设为18,最高压力为60 hPa。地图投影选择兰伯特投影。大气模型步长设为150 s,辐射模型步长设为30 min,陆面模型被调用时间间隔为450 s,各模型输出的时间频率均设为6 h。模型预处理中,土地利用数据采用BATS土地覆被分类方案进行分类,输出的土地利用文件包含被分配到每个格点的土地类型。边界条件上的横向缓冲区宽度设为12个格点,边界层方案设为Holtslag PBL行星边界层;大尺度降水方案选择次网格显示水汽方案(SUBEX);积云对流方案选择MIT-Emanuel方案以适配东亚季风区水汽输送特征。模型运行前进行为期2 a的Spin-up模拟,以消除初始场的不平衡状态,确保各物理变量(如土壤湿度、地表温度等)达到动态平衡。
2 结果与分析
2.1 RegCM 4.7模拟能力检验
模型采用连续动态模拟设计,输入2009—2018年逐年更新的土地利用数据。通过连续10 a的动态积分,模型逐日输出气温、降水等气候变量,完整捕捉土地利用变化的累积效应及气候响应的时空演变。同时,外部强迫(温室气体浓度、气溶胶排放)基于ERA-interim再分析数据集保持一致,以凸显土地利用变化的独立影响。对逐日输出的气候数据取平均值,在ArcGIS中,基于格点化的中国气候观测数据,在时间和空间尺度上对模拟值和观测值进行相关性分析,以检验研究区气温、降水变化对土地利用变化的敏感性。
图1为研究区10 a空间尺度上平均气温和降水的模拟、观测值分布情况。结果表明,RegCM 4.7模型能够基于土地利用变化较好地模拟研究区平均气温和降水的空间分布,符合中国气温“东南地区高,西南、东北地区低”、降水量“南多北少、由东南到西北逐渐递减”的空间分布特征。利用ArcGIS多元分析工具分析栅格尺度气温模拟值与观测值的空间点序列相关性,气温模拟值与观测值的决定系数(
R²)为0.9,降水
R²为0.5,且通过95%置信水平检验。因此,该区域气候模式对气温的模拟性能优于对降水的模拟性能。主要是因为该模式对气候的模拟主要基于能量平衡,而非辐射对流;同时,气温变动有一定规律性,而降水变化缺乏连续性。同时,由
图1e和
图1f可以看出,该区域气候模式对青藏高原及周边的气温和降水模拟误差相对较大,模拟的气温整体偏低于观测结果,而降水模拟值则高于观测结果。可能是因为青藏高原的气候条件特殊、地形复杂,且气象观测站点稀少,已有多位学者
[27-28]对RegCM在该区域模拟性能不足的原因进行系统分析,因此,本文不再对青藏高原区域的模拟性能展开详细讨论。由检验结果可知,气候变化在一定程度上受到土地利用变化的影响,且RegCM 4.7区域气候模式能够基于研究区土地利用变化背景较好地模拟出平均气温和降水量在空间尺度上的分布特征和时间尺度上的变化趋势
[29]。因此,该模式模拟结果可用于探索研究区2009—2018年土地利用变化对气候变化的影响。
2.2 土地利用变化特征
结合研究目的,将研究区土地利用重分类为7种类型,分别为冰川、草地、林地、水域、裸地、农用地和建设用地,对2009年和2018年土地利用数据开展统计分析,得到土地利用数量及空间变化情况。由
表1可知,研究区2009—2018年土地资源置换规模较大,用途转变方式较为复杂。其中,农用地转出规模最大,高达47 759.89 km
2,占转出总面积的32.7%,建设用地占用面积最大,占农用地流出总面积的67.9%;同时,裸地与草地、草地与林地、林地与农用地的置换规模也较大。
图2为研究区土地利用的空间总体变化情况,在空间上呈“东部变动较为剧烈、西部变动较为缓和”的特征。其中华东和华中地区,特别是山东、江苏、河南等省份农用地流失最为严重,且多转化为建设用地,主要原因是这些省份近年来城市化发展迅速,部分农用地因城市化、工业化需求被占用。西部地区裸地、草地流失较为严重,多被农用地、草地侵占;而西南地区草地、水域流失较为严重。新增裸地多集中分布在新疆、西藏、内蒙古、青海等省份的边缘地区,主要由该区域不利的自然环境要素导致。
2.3 土地利用变化对气温的影响
把研究区划分为东北地区、华北地区、华东地区等7大区域,根据模型模拟结果,在
表1与
图2的基础上对比分析7个区域主要的土地利用变化在不同季节引起的气候变化情况。
图3为研究区2009—2018年土地利用变化对夏季、冬季及全年气温影响的空间分布情况。由
图3可知,2009—2018年研究区气温变化呈显著空间分异与季节差异。夏季大部分地区气温上升,其中华东地区平均增幅达0.5~2.0 ℃,局部区域(如山东半岛)最高升温4.0 ℃;冬季则以降温为主,降幅介于0.50~5.45 ℃。该趋势与大多数区域气候模拟研究中的结论相符。该现象主要源于土地利用变化引起的地表物理属性(如反照率、粗糙度等)和水热平衡的调整,进一步影响低层大气的环流模式。
模拟结果显示,2009—2018年研究区气温变化呈显著的区域分异。东北地区夏季气温升高显著,是因为该区域农用地与林地的置换规模较大,导致区域地表蒸散发减少,地表反照率上升,从而造成增温效应,与华文剑等
[6]的模拟结果一致。东北地区的纬度更高,植被覆盖对土壤湿度的调节作用较为显著。华东地区(如山东、江苏等)大量农用地转为建设用地,驱动气温升高。西北地区由于林地和裸地被草地大量侵占,夏季存在明显增温,增幅最高为4.37 ℃,但干旱区水汽限制潜热释放,夜间长波辐射损失增加,冬季降温幅度可达1.00~5.45 ℃。西南地区气候变动较为复杂,特别是西藏、云南的部分地区在夏、冬季气温的变化完全相反。其中,西藏部分地区夏季气温降低,冬季气温升高。2处变动区基本位于同一位置,该区域存在水域和草地面积的增加,夏季蒸发增强抑制增温,冬季地表反照率下降吸收更多辐射,导致气温夏降冬升;而云南气温变动相反区域,也存在不同程度草地、农用地的增加。华中地区不同省份在夏季和全年尺度上气候变化存在差异,其中,以河南省和湖南省差异性最为显著。河南农用地转建设用地驱动夏季升温,气温增幅为0~2 ℃;而湖南省农用地面积显著增加,多为林地转入,引起该省份夏季气温的下降。冬季华中地区整体降温显著,最大下降4.26 ℃。华南地区林地面积的增加,使得区域冬季气温显著降低,降幅为0.5~4.0 ℃。华北地区林地和草地虽有大面积增加,但草地、农用地也有大面积流失,使得该区域整体气温在研究期内无明显变化。
2.4 土地利用变化对降水的影响
从
图4可以看出,2009—2018年中国降水变化呈显著的空间异质性与季节性。东部与南部地区的降水变化显著高于西部地区,且夏季降水量变化更为显著。降水变化受土地利用类型转换的直接影响,城市化通过降低地表动力属性来抑制水汽输送,而植被扩张通过潜热释放增强局地对流,与CAO等
[30]的东亚季风区降水调控机制一致。
从不同区域来看,2009—2018年研究区降水变化呈现显著的区域分异:东北地区的多种土地利用变化方式,导致该区域降水量存在两极分化,辽宁中部局部减少30~165 mm,而吉林东部局部增加30~165 mm,表明土地利用变化会对局地气候的水热平衡产生重要影响。华东地区和华中地区冬季降水量变化不明显,大部分省份主要土地变化类型为农用地转为建设用地,驱动夏季降水减少50~389 mm(如河南省和湖北省)。但湖南、江西、福建三省在大量林地转农用地后,夏季与冬季降水量增加显著,增幅为5~280 mm。华南地区林地转入较多,促进对流降水,夏季和冬季降水均有显著增加,全年尺度平均降水量变化基本与夏季降水量变化情况一致。青藏高原西部受裸地转草地影响,驱动夏季降水量增加10%~15%;云南地区在草地转林地的大部分区域,夏季、冬季及全年平均降水量均有所增加。降水量的变化不仅与土地利用变化有关,而且与地表能量通量、区域蒸散发情况和水汽辐合通量等也密切相关。
2.5 气温、降水的变化机理分析
2.5.1 地表能量平衡分析
为进一步分析模拟结果,从地表能量平衡差异与近地面风场差值来研究气温与降水的变化机理,以探究土地利用变化对区域气候变化的影响方式。土地利用变化引起的区域变量差异变化见
表2。东北地区与华北地区土地资源置换情况相似,2个区域农用地与林地的置换规模较大。该土地利用方式的转变使得区域内植被面积与密度减少,地面粗糙度下降261 μm,导致地表反照率下降1.6%,地面吸收太阳辐射增加1.84 W/m
2,植被蒸腾作用增强,使对流层大气直接增温并形成逆辐射,引起区域年平均气温升高。
华中地区是7个区域中地表反照率唯一增加的区域。区域内大面积农用地被建设用地侵占,大量林地又转为农用地,使得地面粗糙度增加167 μm,地表反照率平均增加0.5%,地面吸收的太阳辐射减少2.71 W/m2,植被蒸腾作用减弱,蒸发通量减少128.11 mm,导致区域平均气温10 a内降低0.33 ℃。华东地区建设用地大面积占用农用地后,区域地表反照率与地面粗糙度有所下降,但下降幅度不显著,气温增加幅度也相对较小。因此,建设用地侵占农用地所引起的增温效应并不显著。西北地区与西南地区裸地大面积转为草地后,区域地面粗糙度和地表反照率均有不同程度下降,地面吸收太阳辐射增加,但平均气温略有下降。因此,裸地转为草地后引起地表反照率、地面粗糙度的下降,但所引起的降温效应不及华中地区显著。
2.5.2 低层环流场分析
图5为冬季和夏季的近地面10 m的风场差值空间分布情况。模拟结果表明,土地利用变化显著改变研究区低层环流场的空间格局,且季节性差异显著。冬季,西南地区北部近地面东南风增强,西北地区(青海、新疆)西风增强0.8~1.6 m/s;夏季,华南地区西南风减弱,西藏高原北风减弱。环流变化的季节性差异显著,表明夏季风场扰动强度高于冬季。
低层环流场调整与土地利用转型的地表动力属性变化密切相关。就研究区冬季而言,西南地区及西北地区低层环流场对土地利用变化的响应较强。西南地区东南风增强,该区域草地转农用地导致地面粗糙度下降,地表摩擦阻力减小,加速东南风向内陆输送,导致北部降水减少,而喜马拉雅山脉附近受海风影响降水增多。西北地区青海省风力差值较大,主要以西风和西北风的增强为主,使得上空气流向东部及南部运动,导致东南部地区冬季降水量增多。同时,该地区以新疆为中心形成一个反气旋环流,该省份土地利用变化使得近地面气流辐散,空气以下沉运动为主,抑制云雨生成,导致冬季平均降水量下降。就研究区夏季具体而言,华南地区和西藏地区大气环流对土地利用变化的响应相对较强,东北、华东、华中地区响应较弱。西北地区仍以新疆为中心,地面土地利用变化导致该区域形成较小的反气旋,夏季降水量减少。西藏地区夏季盛行西北风和北风,但该区域近地面风向受山脉地形影响严重,季风被阻挡后风力减弱,气流被抬升,增加了区域降水量。土地利用变化通过改变近地面风场动力结构,影响水汽输送与辐合强度。建设用地扩张导致粗糙度下降,从而加速干冷空气入侵,而植被扩张会阻碍季风水汽输送。
3 结 论
1)RegCM 4.7区域气候模式能根据研究区土地利用变化情况较好地模拟气候变化情况,气温模拟值与观测值相关系数为0.9,而降水模拟值与观测值相关系数为0.5。因此,该区域气候模式对气温的模拟性能优于对降水的模拟性能,模式模拟结果可用于探索土地利用变化对气候变化的影响研究。
2)土地利用变化对气候的影响存在明显的季节差异。夏季大部分地区增温显著,冬季降温显著,且南北方降水量差异较大,有明显的空间分异。
3)不同土地利用类型变化所引起的气候效应不同。研究区农用地转为林地后引起夏季气温的升高和降水量增加,而林地转为农用地则引起降温效应,农用地转为建设用地也引起气温的升高,但降水量有所下降,裸地转为草地也带来降水量大幅的增加及降温效应,但所引起的降温效应不及林地转农用地所带来的降温效应强。
4)土地利用变化主要通过改变下垫面等方式影响地面粗糙度、地表反照率等物理通量的变化,从而影响地表水热状况。华东地区农用地转建设用地通过地表反照率降低驱动增温,而华南地区林地扩张通过潜热通量增强促进降水增加,揭示土地利用类型转换对水热平衡的差异化调控路径。
5)土地利用变化会对大气环流场造成一定程度的影响,进而导致区域降水量增多。华东地区建设用地扩张导致粗糙度下降,从而加速干冷空气入侵,抑制水汽辐合,减少降水;而西藏草地转农用地则通过增强东南风促进冬季降水增加,表明土地利用变化对环流的扰动具有显著的区域与季节依赖性。
4 不足与展望
由于时间与条件限制,本文无法尝试模型中的每种降水参数化方案,且降水过程复杂,不确定性较大。今后研究,应对降水方案的参数进行优化,以提高该模式对降水的模拟性能,使研究结果更加精确。本文研究的重点是区域土地利用变化对气候的影响分析,未针对相关影响提出针对性对策。今后应对研究区未来土地利用变化产生的气候效应进行模拟研究,并根据具体结果,基于研究区不同区域的经济、生态及土地利用实际情况,从土地利用的角度提出应对未来气候变化的具体对策。本研究通过固定外部强迫来减少干扰,但自然气候变率(如年际振荡)的影响仍需通过多模式集合模拟进一步分离。