土地利用/覆被变化(Land Use/Cover Change, 简称LUCC)是全球变化的重要驱动因素, 是迄今人类改变地球系统最重要的方面之一。土地利用类型的变化将影响着地表反照率、地面粗糙度及土壤的热力、水分特性等地表属性, 从而对气候以及陆面过程产生极其重要的影响^[1-2]。在区域尺度上, 植被覆盖变化也能对地表能量平衡和水分平衡产生实质性的影响^[3]。因此, LUCC能够显著影响区域气候变化已经毋庸置疑^[4-5]。而且对于区域尺度而言, 多数地区气候内部变率的贡献大于外强迫, 区域尺度气温变化的机制较全球、半球尺度更为复杂^[6]。因此, 开展LUCC对区域气候变化的影响研究, 对明确区域气候变化机理具有重要的科学意义。

LUCC的区域气候效应研究的主要方法有两大类。一是利用统计学分析法, 分析观测数据与再分析资料之间的关系, 从而揭示LUCC对区域气候的影响^[7]。二是运用陆-气耦合模式, 结合模拟试验, 研究LUCC产生的区域气候效应^[8-10]。后者因具有良好的物理基础, 能够揭示区域气候变化的机制而被广泛应用^[11-14]。例如张学珍等^[15]利用WRF(Weather Research and Forecasting)模式, 研究了中国中东部耕地扩张对表面气温的影响, 结果发现农业扩张对大尺度的区域温度影响较小；王明娜等^[5]利用区域模式WRF模拟分析了21世纪初期2001—2010年中国北方半干旱区土地利用变化对地表温度的影响; 董思言等^[16]利用区域环境集成模拟系统模拟了土地利用/覆被变化对中国不同季节气温的影响；陈海山等^[17]模拟了近20年中国土地利用变化对区域气候的影响。多数研究表明, 由于正、负温度效应在区域平均的过程中相互抵消, 使得LUCC对区域平均温度变化的影响甚微^{[5, 15]}, 主要表现在同一时间上, LUCC引起某些土地利用类型的温度增加, 某些土地类型的温度却减少；而在同一空间范围内, 不同时间尺度也表现出不同的温度变化趋势。可见, LUCC对区域气候的影响主要限于局地^{[5, 16]}。然而, 多数研究主要集中在较大的区域(大于1×10⁶km²), 且空间分辨率多集中在30km及以上水平, 时间分辨率多介于季节和年尺度之间^{[16, 18]}, 这从很大程度上影响了研究结果的精度。另外, LUCC对区域气温的影响研究较多, 而对局地降水的研究相对较少。

同时, LUCC对区域气候变化产生的重要影响, 也势必会影响水文过程, 即为LUCC的间接水文效应。但是, 多数研究注重探索LUCC引起产流、汇流等直接水文过程的变化^[19-20], 鲜见有探讨LUCC对水文过程产生的间接影响的研究。因此, 本研究基于卫星观测的研究区1990年和2010年LUCC客观事实, 利用WRF模式, 以鲁南局地(沂河流域)为例, 从更高空间(3km)和时间分辨率(月和日尺度)尺度, 探析LUCC对局地气温和降水变化的影响；并采用弹性系数法, 定量评估LUCC产生的间接水文效应。研究结论对准确解读LUCC对局地气候的影响机制, 明确局地水文过程变化机理具有重要的科学指导意义。

1 材料与方法 1.1 研究区概况

沂河是淮河流域沂沭泗水系中较大河流, 位于33.5°—36.5°N, 117°—119.0°E之间(图 1), 北接海河平原, 南邻黄淮平原, 为临沂地区的母亲河。沂河源出山东省沂源县, 经沂水、沂南、临沂、蒙阴、平邑、郯城等县、市, 至江苏省邳县吴楼村入新沂河, 抵燕尾港入黄海, 流域面积1.7×10⁴km²。

1.2 试验设计与数据来源 1.2.1 基于WRF模式的LUCC区域气候效应

采用WRF 3.7研究LUCC引起的区域气候要素的变化, 从而揭示LUCC的间接水文效应。由于气温和降水是影响径流变化的主导气候因子, 因此本文主要利用WRF模式研究LUCC引起的区域气温和降水的变化, 从而分析其对径流的影响。

本研究主要包括控制试验和敏感性试验。模拟区域的中心经纬度为35.8°N、118.1°E。试验采用3层嵌套方式, 水平分辨率分别为27km(domain1)、9km(domain2)和3km(domain3)(图 1), 对应的格点数分别为27×36、49×67和91×130。考虑到短时间内LUCC基本稳定及数据的可获取性, 采用2010年LUCC和2013年的驱动场作为控制试验, 模拟2013年不同月份的气温和降水状况。敏感性试验是将研究区(domain3)2010年LUCC数据改为1990年。除了LUCC不同之外, 所有试验采用了相同的设置和物理参数化方案, 所以控制试验与各敏感试验模拟结果之差反映了研究区LUCC对区域气候要素的影响。

由于研究区气候季节变化明显, 本文设置模拟时间(UTC)为2013年1月1日00时—2013年2月1日00时；2013年4月1日00时—2013年5月1日00时；2013年7月1日00时—2013年8月1日00时；2013年10月1日00时—2013年11月1日00时。这样一方面可以反映WRF模式对不同季节天气特征的模拟状况, 同时可以展示LUCC对不同天气动力状况下的气候要素的影响。另外, 参考前人研究成果^[15]及敏感性试验, 本试验采用的微物理过程方案为Lin等的方案, 长波辐射方案为rrtm方案, 短波辐射方案为Dudhia方案, 近地面层方案为Monin-Obukhov方案, 陆面过程方案为Noah方案, 边界层方案为YSU方案, 积云参数化方案为浅对流Kain-Fritsch(new Eta)方案。初始场和侧边界使用的是NCAR提供的1°×1°的NCEP-FNL再分析数据。

本文假设LUCC引起气候要素变化量与实际气温或降水量大小呈正比, 那么可以有：

(1)

式中, δ是研究区1990—2010年LUCC引起的2013年气温或降水变化量；ΔX₄为WRF模式模拟的LUCC引起研究区2013年1、4、7和10月气温或降水的变化量之和；X_ms2013为用1990年LUCC模拟的2013年1、4、7和10月平均气温或降水量之和；X₂₀₁₃为研究区2013年实测平均气温或降水总量。

本研究的模拟运算采用WRF 3.7并行版本, 在曙光天阔服务器1620-G15(32CPU)上完成。控制试验和敏感性试验的外层(domain1和domain2)均利用USGS(美国地质勘测局)10′和2′土地利用数据, 内层(domain3)用中国科学院资源环境科学数据中心(http://www.resdc.cn) 2010年1km×1km土地利用数据(重分类、重采样和格式转换后)和1990年土地利用数据。运用研究区实测气象站数据验证控制试验模拟结果。实测气象资料包括兖州、沂源、莒县和徐州气象站的逐日气温和降水数据。

1.2.2 径流对气候要素变化的弹性系数

在气候变化对水文过程影响研究中, 通常采用弹性系数来表征水文要素对气象要素的敏感性。本文采用Zheng等^[21]提出的弹性计算方法分析1960—2013年径流量对降水或气温变化的弹性系数。计算公式如下：

(2)

式中, X_i为气象要素, Q_i为径流量, ε为弹性系数, X和Q分别为气象要素和径流量的多年平均值。ε的物理意义是指气象要素变化1%, 引起径流的变化率ε%。

因此, LUCC引起2013年径流变化率(φ)为：

(3)

1.2.3 Mann-Kendall单调趋势检验

本文运用Mann-Kendall单调趋势检验判定径流和气候要素在某一时间段的变化趋势是否显著。本文取P= 0.05的显著性检验水平, 若时间序列在此置信水平下存在显著变化趋势, 则|Z|=1.96, Z为正值, 表明时间序列具有上升或增加趋势, Z为负值, 则表示具有下降或减少趋势。

2 结果与分析 2.1 控制试验模拟结果检验

基于控制试验模拟的不同月份逐日气温与实测数据之间的关系可知, WRF模式能够较好的反映研究区各月气温状况, 模拟温度与实测温度之间的相关性均达到了P < 0.001显著性水平(图 2)。当然, 模拟值与实测值之间也具有一定的误差。从月尺度上, WRF模式对2013年1月和7月模拟的平均气温略低于实测值, 误差率分别为9.1%和6.2%；而对4月和10月模拟的平均气温略高于实测值, 误差率分别为1.8%和3.3%。总体来看, WRF模式对研究区气温模拟结果可靠, 可以较好的反映气温的时空变化特征。

图 3展示的是控制试验模拟的不同月份逐日降水与实测数据之间的关系。结果表明, 虽然WRF模式对降水的模拟精度低于气温, 但模拟降水与实测降水之间的相关性均达到了P < 0.05显著性水平。从误差率来看, WRF模式对2013年1月和4月模拟的总降水量略低于实际值, 误差率分别为4.8%和14.3%；而对7月和10月模拟的总降水量略高于实测值, 误差率分别为6.7%和17.1%。可以看出, WRF模式对研究区降水模拟结果可信。

2.2 土地利用/覆被变化分析

从各土地利用类型所占面积比例来看, 1990和2010年, 研究区土地利用类型均是以旱地为主, 分别占到总面积的70.3%和68.6%, 其次是水田, 分别占总面积的11.4%和11.1%, 草地和有林地所占面积比例介于5.4%—7.2%之间, 而其他土地利用类型所占面积比例较小(图 4)。

从土地利用类型的变化量来说, 1990—2010年, 旱地面积减少量最大, 为1089.8km²。其次是水田和草地, 分别减少了198.3km²和62.2km²。而城市建设用地增加量最大, 为821.8km², 其次为裸地和水体, 分别增加了207.8km²和214.2km²。另外, 永久湿地(75.0km²)和疏林地(55.6km²)面积略有增加, 而其他土地利用类型的面积变化量较小。由此可以看出, 近20多年来, 研究区土地利用类型主要是以旱地向建设用地(747.3km²)和裸地(132.4km²)转化的过程。

2.3 LUCC导致的气象要素变化 2.3.1 气温变化

图 5展示的是由控制试验与敏感性试验模拟气温的差值, 其结果表明由近20年的LUCC引起的不同月份气温空间变化。可以看出, 由于LUCC引起的1月和10月气温均增加了0.2℃；而7月气温呈现降低态势(0.1℃), 另外4月气温基本稳定。从各土地利用类型的平均温度来说, 1月仅有城市建设用地平均温度略有降低, 其他各土地利用类型平均温度均表现为上升态势, 其中混交林气温上升幅度最大, 为0.8℃, 其次是水体, 温度升高了0.4℃, 而落叶阔叶林和湿地气温上升幅度较小, 均小于0.1℃。4月仅有旱地温度略有增加, 但幅度很小, 其他各土地利用类型平均温度均降低, 其中混交林温度减少幅度最大, 为1.0℃, 其次是水体, 温度降低了0.3℃, 而湿地、落叶阔叶林、水田和草地温度降幅较小, 小于0.1℃。7月所有的土地利用类型温度均呈减少态势, 其中仍是混交林和水体降温幅度较大, 分别为1.3℃和0.3℃, 其次是城市建设用地和灌木林地, 降温幅度介于0.2—0.3℃, 而湿地气温减少量最小, 小于0.1℃, 其他旱地、草地、水田等土地利用类型温度降幅在0.1—0.2℃之间。10月所有的土地利用类型的温度均表现为上升态势, 其中混交林和水体升温幅度最大, 分别为0.7℃和0.3℃, 其次是灌木林地, 温度增幅在0.2℃左右, 而其他土地利用类型的温度增加幅度较小, 小于0.1℃。

2.3.2 降水变化

控制试验与敏感性试验模拟降水的差值结果表明(图 6), 近20年的LUCC引起的1月和4月降水变化量很小, 可以忽略不计；10月降水略有增加, 增量为0.3mm；而7月由LUCC引起的平均降水减少幅度较大, 降低了23.7mm。可见, LUCC对夏季降水影响较大, 这是因为研究区夏季降水量相对丰富, LUCC变化可以从一定程度上影响水汽通量的变化。

2.4 LUCC的间接水文效应

基于1960—2013年研究区年径流和气象资料, 利用公式(2)计算得到沂河流域年径流量对年降水变化的弹性系数为2.4, 意即研究区年降水量变化1%, 将引起年径流量变化2.4%。而研究区年径流量对年均气温变化的弹性系数小于降水, 为-1.8, 这表明气温升高1%, 将导致径流减少1.8%。这也说明研究区降水对径流的影响程度大于气温。

同时, 基于公式(3)可以得出1990—2010年LUCC引起2013年降水和气温变化导致的沂河年径流量的变化率分别为18.4%和1.7%。可以看出, LUCC的间接水文效应主要表现在对降水变化的影响。

3 讨论

WRF模式模拟结果表明, 近20年来研究区LUCC导致局地气候发生了明显变化, 而且这种变化对径流量产生了重要影响。因此, LUCC的水文效应不仅表现在改变地表植被的截留量、土壤水分的入渗能力等产流、汇流方面, 而且体现在通过影响局地降水和蒸发等气象要素而间接影响径流。因此, LUCC导致的间接水文效应值得进一步关注, 这也是我们未来工作的重点。在研究方法上, 可考虑基于气候模式与具有良好物理机制的水文模型耦合的方法, 全面解析LUCC产生的气候和水文效应。

虽然WRF模式模拟结果通过了显著性检验, 但本文仅模拟了各季节的代表性月份LUCC对局地气候的影响, 这给模拟结果带来一定的不确定性, 因此要加强更长时间尺度的模拟与研究, 进一步准确揭示LUCC产生的气候效应。同时, 研究结果表明, 不仅温度在区域平均的过程中会因正、负效应而相互抵消, 而且降水也存在类似现象, 这与前人研究结果一致^{[5, 15]}。因此, 应从更高空间和时间分辨率角度科学揭示LUCC对局地气候的影响。

另外, 模拟结果表明, 1990—2010土地利用类型变化使得研究区城市建设用地在个别月份的温度略有下降, 可以有以下3个方面的探讨：一是尽管城市建设用地温度略降, 但是LUCC导致的区域温度总体还是以升高为主, 而城市建设用地仅是城市土地利用的一种类型, 因此, 至于LUCC是否导致研究区城市“热岛”整体减弱的还有待进一步研究；二是热岛效应不仅表现在强度上, 也体现在空间分布上, 虽然城市建设用地热岛减弱, 但在空间上热岛范围在扩张；三是研究表明, 在2005年以后我国东部许多地区的城市热岛出现减弱的趋势^[22], 例如葛荣凤等^[23]研究发现北京二环区域的热环境呈现一定程度的好转, 张晓莉等^[24]指出1987—2013年西宁城区热岛强度小幅减弱, 韩贵锋等^[25]研究显示近20年来重庆部分城区的热岛强度明显减弱。这主要是因为近年来我国城市空间格局优化、城市绿化面积增加等现象使得人为因素对城市环境的不利影响在减弱^[22], 人居环境改善^[25]。因此, 近20年来研究区LUCC导致城市建设用地温度在个别月份略有降低也是可以理解的。

4 结论

(1) 1990—2010年, 研究区土地利用变化明显, 主要是以旱地向建设用地(747.3km²)和裸地(132.4km²)转化的过程。

(2) 近20年来, 土地利用类型变化的区域气候效应明显且具有一定的时空差异性。WRF模式模拟结果表明, LUCC引起2013年1月和10月气温都增加了0.2℃, 使得7月气温却减小了0.2℃, 而4月气温基本稳定。其中, LUCC导致混交林和水体的温度变化幅度较大, 在0.4—1.3℃之间, 而对落叶阔叶林和湿地气温的影响较小, 其温度变化幅度均小于0.1℃。降水方面, 1990—2010年LUCC引起1月和4月的降水变化几乎可以忽略不计, 导致10月降水增加了0.3mm, 而引起7月降水变化较大, 即减少23.7mm。

(3) 1960—2013年, 沂河径流对流域气候要素变化响应敏感。流域降水和气温变化1%, 将会导致径流分别改变2.4%和1.8%。同时, 1990—2010年, LUCC的间接水文效应显著, 主要表现在LUCC使得2013年径流量减少了20.1%。