全球气候变化加剧了暴雨发生的频率与强度, 快速城市化引起地表水文循环过程改变, 两者共同作用使得城市内涝问题日益突出。城市硬质地面阻碍了雨水截留与下渗, 将其转化为可渗透性生态铺装是减少暴雨径流、缓解内涝的有效途径。建筑屋面约占城市不透水地面的40%—50%^[1], 实施屋顶绿化能够充分利用闲置空间对治暴雨内涝, 促进城市可持续发展。

国内外针对绿化屋顶水文效应的定量研究主要基于野外实验观测方法获取降雨事件中的雨量和径流量动态数据, 重现降雨-滞蓄-径流过程, 揭示产流及洪峰延缓效应, 探索关键影响因子。由于研究区降雨特征及植被土壤性质不同, 文献报导的绿化屋顶径流效应之间存在较大差异, 径流削减率从10%到100%不等^[2-9], 洪峰流量可降低18%—100%^{[1, 8-10]}, 产流出现时间延迟5—130 min^[11-12], 洪峰延迟3—55 min^{[8-9, 13]}, 雨量、雨强、土壤深度及前期土壤湿度是影响径流削减效益的关键因子^{[3, 13-15]}。因此, 单个研究结果不具有普适性, 针对特定气候区及绿化类型开展专门研究对指导当地屋顶绿化实践及雨洪管理具有重要意义。

由于建筑屋面排水系统复杂, 实时观测径流量难度较大, 大多数研究采用定时收集排水管下落雨水并根据重量或水位推算径流深度^[13-14], 或采用自制绿化屋面微模型替代现实绿化屋顶的方法以简化观测^[16-20], 前者费时费力, 难以获得长时间序列观测数据, 后者存在尺度上推等不确定性问题。水文经验方程的引入可弥补上述观测实验缺陷, 且能分拆植被、土壤截留及蒸发等各个部分对径流削减的贡献, 有助于了解水文效应形成机制, 因此, 基于现场观测的经验方程法在大面积绿化屋顶径流削减效应研究中得到广泛应用^[21-23]。

SCS-CN模型是由美国农业部水土保持局(USDA-SCS)开发的用于流域径流量估算的水文经验模型, 具有结构简单、物理概念明确、所需参数少等特点。该模型假定流域降雨-径流关系与径流曲线数(CN)有关, CN值又由植被覆盖、土壤质地及入渗能力、前期土壤湿度等下垫面特征综合决定, 取值介于0—100之间, CN值越高, 径流系数越大。SCS基于对美国植被土壤特征的大量研究建立了CN值检索表, 为SCS-CN模型的广泛应用提供了参数依据, 但由于CN值具有时空变化特征, 在实际应用过程中进行参数区域化是保证模拟结果可靠性的重要环节。国内学者对特定流域及地表覆盖的CN取值做过大量研究^[24-27], 然而, 目前针对绿化屋顶这种特殊下垫面类型的研究较少。

亚热带季风区夏热冬冷, 降雨丰富且存在明显的季节分布规律, 以南京为例, 基于全年观测数据及屋顶水量平衡经验方程, 分析简易型和花园型两种典型绿化屋顶的径流削减效应及其随季节、雨型的变化规律, 探讨关键影响因子, 基于全年降雨-径流数据反推CN值, 比较绿化屋顶与其他地表覆盖类型的径流削减效果, 并为大尺度(流域或集水区)屋顶绿化水文效应分析提供关键输入参数。

1 材料与方法 1.1 实验场地

南京气候四季分明, 夏热冬冷, 年最高气温39.7 ℃, 最低气温-13.1 ℃, 由于地处亚热带季风气候区, 降水丰富但不均匀, 年平均降雨量约1200 mm, 但70.6%的降雨量集中在5—9月汛期^[28], 期间台风暴雨频发, 城区积水内涝严重, 给城市生产生活造成诸多不利影响, 亟需有效的雨洪管理措施。本研究实验场地位于南京市栖霞区紫东国际创意园, 园区以商业办公建筑为主, 大部分建筑屋顶已实施绿化, 总面积达3.6 hm²。选择简易型(B3)和花园型(F4)两个绿化屋顶开展观测实验(图 1), 两个绿化屋顶与建筑同步规划建设, 确保屋顶荷载满足绿化要求。B3土壤基质浅, 种植南方地区屋顶绿化最常见的佛甲草；F4土壤基质厚, 能支撑小乔木、灌木和草本多层植被结构, 生物量丰富(表 1)。两个绿化屋顶自下而上包括原有建筑屋面、防水层、防根层、蓄排水层、过滤层、基质层及植被层(图 2)。

实验观测时间为2016年12月8日至2017年12月7日, 期间两个屋顶植物生长状态良好, 覆盖度接近100%。首先在B3屋顶上设置自动气象站观测太阳辐射、雨量、风速风向、温湿度, 提供两个实验场地的背景气象信息；其次, 测量B3和F4的土壤湿度, 作为土壤蒸发与蓄水量的计算依据。

自动气象站和其他传感器设备均放置在屋顶中心区以减少边缘效应, 传感器的读数间隔统一设定为15 min。图 2为两个绿化屋顶纵向结构、测量环境因子的设备名称及型号。

1.2 径流量估算方法

根据屋顶水量平衡方程, 径流量等于降雨总量扣除植被、土壤、排水层截留雨水量及土壤水分蒸发量, 计算公式如下：

(1)

式中, Q为一次降雨事件的径流量(mm), R为降雨量(mm), V为植被截留量(mm), S₁为土壤截留量(mm), E为土壤蒸发量(mm), S₂为排水层结构截留量(mm)。

R通过翻斗式雨量仪测量, 其余参数通过经验公式(2)—(6)估算, 其中植被截留公式为^[29]：

(2)

式中, a为叶面上平均最大持水深度(mm), 取值为0.2；Veg为植被盖度, 本研究两场地已实施屋顶绿化4年以上, 植被覆盖度高, 故取值为1；LAI为叶面积指数, 根据两个屋顶植被类型及前人测量的经验值范围, 确定简易型绿化屋顶LAI值为2.5^[30], 花园型为4.0^[31-32]。

土壤截留量通过式(3)计算：

(3)

式中, ΔW₁为每次降雨结束后土壤湿度增加值(m³/m³), H为土壤厚度(mm)。

蒸发量估算方法为：根据非降雨天土壤湿度日变化值, 估算全年土壤日蒸发量, 并统计春、夏、秋、冬四个季节平均每小时蒸发量, 再根据降雨发生的季节及历时推算土壤蒸发量, 计算公式如下：

(4)

(5)

式中, E为一次降雨事件中绿化屋顶的蒸发量(mm), E₁为非降雨天的日蒸发量(mm), ΔW₂为非降雨天每天土壤湿度减少量(m³/m³), H为土壤厚度(mm), T为降雨历时(h)。

结构截留主要考虑绿化屋顶排水层窝状结构蓄水量, 根据单个窝状的体积及分布密度, 推算其最大蓄水量为4.3 L/m², 换算成深度为4.3 mm^[7]。当降雨量全部被植被和土壤截留时(R < V+S₁), 结构截留量为0 mm；当降雨量超过植被截留与土壤截留最大值时, 结构截留开始发挥作用；随着降雨量不断增加, 结构截留达到最大值4.3 mm后将不再变化, 计算公式如下：

(6)

本研究选择暴雨、大雨、中雨和小雨共四种典型降雨分析绿化屋顶的雨水截留效应, 一次降雨事件的雨量>50 mm为暴雨, 25—49.9 mm为大雨, 10—24.9 mm为中雨, 雨量 < 10 mm为小雨^[33]。一次降雨事件定义为连续性的、降雨停止时间间隔不超过6 h的降雨。

1.3 CN值估算方法

SCS-CN模型中CN值与降雨量(R)、径流量(Q)的关系如下：

(7)

(8)

式中, S为潜在最大蓄水能力(mm), 本研究根据全年77次降雨事件的降雨量和径流量反推S值和CN值, 再取算术平均值得到绿化屋顶的CN值^[34-37]。

2 结果与分析 2.1 全年降雨事件统计

表 2显示, 实验场地全年共降雨77次, 总降雨量为1307 mm。其中, 小雨49次, 降雨频次最高且历时最长(326.5 h), 降雨总量为188.6 mm, 占全年降雨总量的14.4%；暴雨次数最少但历时短、雨量大, 6次暴雨事件历时124 h, 降雨量达674.8 mm, 占全年降雨总量的51.6%；大雨和中雨降雨次数分别为6次和16次, 降雨量占比分别为13.4%和20.6%。

降雨季节分布规律显著：春季降雨次数最多, 共23次, 以小雨和中雨为主, 无暴雨事件, 降雨总量为207.8 mm；夏、秋两季降雨次数均为21次, 是暴雨的主要发生季节, 其中夏季降雨总量最大, 达557.4 mm, 秋季次之, 为389.2 mm；冬季降雨次数及总量均最少, 12次降雨的总雨量为152.6 mm。

2.2 径流削减效应 2.2.1 全年总况

简易型绿化屋顶全年通过蒸发及植被、土壤、结构截留的降雨量分别为42.6、38.4、282.8、184.7 mm, 径流量为758.5 mm, 径流削减率为42%；花园型绿化屋顶全年通过蒸发及植被、土壤、结构截留的降雨量分别为82.8、61.2、491.2、157.8 mm, 径流量为514 mm, 径流削减率为60.7%(表 3)。以上数据表明, 土壤基质层对雨水的吸收与储存在绿化屋顶径流削减效应中发挥主导作用；此外, 花园型绿化屋顶由于植被结构复杂、土壤基质深厚, 因而在蒸发、截留雨水方面与简易型相比表现出明显的优势。

2.2.2 季节变化规律

绿化屋顶在春、冬两季的径流削减效应明显优于夏秋两季：简易型绿化屋顶春、夏、秋、冬四季的径流削减率分别为78.6%、32.9%、33.2%、47.5%, 花园型屋顶分别为98%、48.1%、49.5%、84.3%。春、冬两季以小雨和中雨为主, 降雨总量小、强度弱、历时长, 雨水缓慢进入场地, 被植被与土壤截留后不断蒸发到空气中, 通过径流排出的雨量较小；而夏秋两季暴雨事件频发, 绿化屋顶在短时间内蓄积大量雨水, 使得植被、土壤、结构层蓄水很快达到饱和, 多余的水分从场地排出, 因而径流削减率较低。

土壤截留在四季均发挥主导作用(图 3), B3土壤截留占降雨量的15%—42%, F4土壤截留占比为30%—58%, 最大值发生在春季, 最小值在夏季或秋季；结构截留占比仅次于土壤截留, 在B3和F4场地分别占9%—26%和10%—21%, 表明简易型绿化屋顶的排水结构层发挥的作用比花园型更大；土壤蒸发量占降雨量比重略高于植被截留量, 前者在两个场地的占比分别为3%—6%和5%—10%, 后者占比分别为2%—6%和3%—9%。

2.2.3 随雨型变化规律

径流削减率随降雨量增加呈递减趋势, 简易型绿化屋顶全年对暴、大、中、小雨的径流削减率分别为14.6%、38%、73.1%、99.1%, 花园型屋顶分别为30.6%、91.4%、90.1%、97.8%。花园型屋顶的径流削减优势在中雨、大雨和暴雨条件下表现突出；小雨时, 由于两类绿化屋顶几乎可以截留全部降雨, 差距较小。

径流削减各部分占比也随降雨量大小而发生变化(图 4), 其中土壤截留在大部分降雨条件下占主导作用, B3土壤截留在中雨时占比最高(43%), 暴雨时最小(9%)；F4土壤截留在大雨时占比最高(72%), 小雨时最小(19%)。植被截留与土壤蒸发量随雨量增加而减少, 小雨时B3场地两者分别占降雨总量的13%和10%, F4场地分别占21%和19%, 暴雨时两场地植被截留和土壤蒸发量均降到5%以下。结构截留随雨量增大而减小, 在小雨时占比甚至高于土壤截留量, B3和F4达到了39%, 暴雨时降到4%。

2.3 径流削减率与环境气象因子的相关性分析

采用Pearson相关分析研究径流削减率与初始土壤湿度、降雨量、降雨强度、大气温度、相对湿度、风速及太阳辐射等环境气象因子之间的关系(表 4)。

绿化屋顶径流削减率与降雨量、降雨强度及相对湿度呈显著负相关关系(P < 0.01), 植被与土壤对径流的削减是有限的, 降雨量超过阈值时, 多余雨水全部以径流形式排出, 因而雨量越大, 径流削减率越低；降雨强度越大, 短期内进入植被土壤的雨量越多, 重力作用导致大量雨水迅速排出, 整体径流削减率越小。初始土壤湿度与径流削减率在简易型绿化屋顶存在显著负相关, 而与花园型屋顶的相关性不显著, 可能的原因为, 简易型屋顶土壤深度仅为50 mm, 前期土壤湿度对其蓄水能力有决定作用, 而花园型屋顶土壤厚度达350 mm, 前期土壤湿度对截蓄雨水能力影响有限。大气温度、太阳辐射和风速等微气候因子虽影响土壤与植被的蒸发蒸腾速率, 但因蒸发量占比较小, 因而这些因子与整体径流削减率之间的相关性不显著。

2.4 CN值估算

图 5显示, 简易型绿化屋顶CN值的变化范围是75—99, 平均值为92, 其中最小值75出现在夏季的中雨事件(R=16.8 mm), 最大值99出现在夏季的小雨事件(R=0.4 mm)；花园型屋顶CN值的范围是59—99, 平均值为88, 且CN值在暴雨和大雨条件下均低于简易型屋顶, 径流削减优势突出, 而在小雨和中雨事件中可能高于后者, 优势不明显。尽管两个绿化屋顶在小雨时径流量很低, 但计算得到的CN值往往高于平均值, 说明SCS-CN模型在小雨条件下的径流计算误差相对较大, Fassman-Beck等^[38]建议在降雨量高于2 mm的条件下使用SCS-CN估算绿化屋顶径流削减效应较为适宜, 与本研究结论具有一致性。

国外多名学者对简易型绿化屋顶的CN值做过估算, 基于不同气候区与植被土壤特征的计算结果显示, 绿化屋顶的CN值介于80—96之间(表 5), 本研究计算结果也落在此区间之内。国内学者采用SCS-CN模型估算了不同土地覆盖类型的CN值, 如李常斌等^[35]、史培军等^[21]、鲁华锋^[27]估算草坪的CN值或区间分别为73.6、36—80和39—89, 均小于本研究简易型绿化屋顶的CN值；鲁华锋^[27]同时估算出乔灌草混合种植区的CN值范围为32—86, 也低于本研究花园型绿化屋顶的CN值。与地面绿化形式相比, 绿化屋顶土壤基质较浅, 且下层与建筑相连, 水文过程有较大区别, CN值计算结果表明, 在植被结构相同的条件下, 绿化屋顶的径流削减效应显著低于地面绿地。

3 结论与讨论

本研究基于野外观测数据与水量平衡方程估算径流量, 分析两个典型绿化屋顶的径流削减效应及其与环境气象因子之间的相关关系, 并采用降雨量-径流量数据反推绿化屋面的CN值。与以往基于少数降雨事件的研究相比, 全年观测数据能够消除季节、雨型、前期土壤湿度等因素造成的径流削减率估算的不确定性, 反映绿化屋顶长期滞蓄效果。

研究结果表明, 简易型、花园型绿化屋顶的全年径流削减率分别为42%和60.7%, CN值分别为92和88, 与具有相同植被属性的地面绿地相比, 绿化屋顶的CN值偏高, 径流削减效益偏低。然而, 对于用地紧凑的城市中心区来说, 增加地面绿化较为困难, 考虑到传统不透水屋面(CN=98)占城市地表面积之大, 屋顶绿化不失为缓解城市暴雨内涝的一项重要策略。以南京为例, 其主城区总面积为636 km², 建筑密度约15%^[39], 若60%的屋面实施绿化, 按公式(7)和(8)推导, 简易型和花园型两种情景可减少全年暴雨径流总量2.8 × 10⁷ m³和4.2 × 10⁷ m³, 分别相当于5.6和8.4个玄武湖的库容量, 效果非常可观。

在径流削减效应中, 土壤基质层对雨水截留起主导作用, 简易型和花园型绿化屋顶土壤截留量分别占径流削减总量的52%和62%, 表明以雨洪管理为导向的屋顶绿化实践过程中, 需重点关注土壤因素, 如厚度、成分及质地等。本研究实验场地的土壤基质成分包括园土、泥炭土、珍珠岩、陶粒、蛭石, 属沙壤土类型, 为本地易得的轻质材料, 饱和水分含量下密度约1000 kg/m³, 简易型(50 mm)和花园型(350 mm)绿化屋顶所需荷载分别为50 kg/m²和350 kg/m², 因此, 大部分现有建筑(荷载量>100 kg/m²)能满足简易型绿化的荷载要求, 新建建筑通过加固措施可达到花园型绿化荷载要求。未来研究可重点分析不同土壤配置类型的孔隙度、下渗率、饱和水分含量、田间持水量(Field capacity)等参数, 并结合不同水分含量下的荷载要求, 选择既能减轻建筑结构负担, 又能最大程度发挥径流削减效益的土壤基质形式。

需要注意的是, 本研究中“土壤蒸发量”一项只包括降雨过程中蒸发到空气中的水分, 雨过天晴后的土壤蒸发量已包含在“土壤截留量”一项中, 因此没有单独计算。土壤蒸发量与季节及背景气候条件有关, 它决定了一次暴雨后土壤从水分饱和状态降到植被凋萎点(Wilting point)之间的时间长度, 对制定合理的灌溉方案、确保植物健康生长具有重要指导意义；此外, 土壤蒸发量也与绿化屋顶的节能与微气候调节效应密切相关。因此, 未来研究可对土壤蒸发量参数做进一步探讨, 为屋顶绿化可持续设计与管理提供参考。

亚热带季风区夏热冬冷, 降雨丰富但分布不均, 由于径流削减效应与雨量、雨强、初始土壤湿度等因素密切相关, 也存在明显的季节分布规律, 因此本文研究结论不一定完全适用于其他气候区；研究选取的简易型和花园型绿化屋顶虽具代表性, 但结果仍然只能反映特定植被土壤条件下的屋面水文平衡过程。未来研究可借助发展较为成熟的水文模型(如SWMM等), 开展多尺度、多情景模拟, 深入剖析不同气候区降水条件下, 植被土壤构成各异的屋顶绿化模式对城市水文过程的影响。

致谢: 感谢刘连英、王永旗、曹绪峰、汪文对野外观测提供的帮助。