• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

天津平原区城市化地区蒸散发特征

曹润祥 李发文 李建柱 冯平

曹润祥, 李发文, 李建柱, 冯平. 天津平原区城市化地区蒸散发特征[J]. 水科学进展, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
引用本文: 曹润祥, 李发文, 李建柱, 冯平. 天津平原区城市化地区蒸散发特征[J]. 水科学进展, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
CAO Runxiang, LI Fawen, LI Jianzhu, FENG Ping. Exploring evapotranspiration characteristics in urbanized areas of Tianjin plain[J]. Advances in Water Science, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
Citation: CAO Runxiang, LI Fawen, LI Jianzhu, FENG Ping. Exploring evapotranspiration characteristics in urbanized areas of Tianjin plain[J]. Advances in Water Science, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005

天津平原区城市化地区蒸散发特征

doi: 10.14042/j.cnki.32.1309.2021.03.005
基金项目: 

国家自然科学基金资助项目 51779165

国家自然科学基金资助项目 51879181

详细信息
    作者简介:

    曹润祥(1990—), 男, 河南新乡人, 博士研究生, 主要从事水文学及水资源研究。E-mail: runx.cao@foxmail.com

    通讯作者:

    李发文, E-mail: lifawen@tju.edu.cn

  • 中图分类号: P339;TV11

Exploring evapotranspiration characteristics in urbanized areas of Tianjin plain

Funds: 

the National Natural Science Foundation of China 51779165

the National Natural Science Foundation of China 51879181

  • 摘要: 为揭示城市蒸散发特征,采用考虑人为热影响的SEBS-Urban模型计算天津平原区2015—2017年逐月蒸散量。使用两步验证法对结果进行检验,首先利用MODIS MOD16产品对非建成区模拟结果合理性进行验证,其次利用城市耗水(UWD)模型得到的天津大学卫津路校区月蒸散量对建成区模拟结果验证。在SEBS-Urban模拟结果基础上,检验Budyko方程在城市化地区的适用性。结果表明:① SEBS-Urban模型精度可靠,体现了蒸散发项在城市地表能量平衡与二元水循环中的相合性;②人为热或社会侧的耗水会导致蒸散发量增加,建成区增幅达85%~115%,整个研究区增幅为7.2%~8.7%;③ Budyko方程在研究区具有较好的拟合效果,可应用于城市地区。本研究成果可为城市化地区蒸散发研究提供必要的参考。
  • 图  1  研究区位置示意

    Figure  1.  Location of the study area

    图  2  SEBS反演结果相对于MOD16 A2产品的泰勒图

    Figure  2.  Taylor diagram of SEBS estimated ET vs. MOD16 A2 product

    图  3  UWD、SEBS与SEBS-Urban结果

    Figure  3.  UWD, SEBS and SEBS-Urban simulation results

    图  4  2015—2017年SEBS与SEBS-Urban模型对天津平原地区蒸散发分布模拟结果

    Figure  4.  Annual evapotranspiration distribution modelled by SEBS and SEBS-Urban over 2015—2017

    图  5  2015—2017年Budyko模拟结果分析

    Figure  5.  Budyko monthly simulation results analysis over 2015—2017

    表  1  SEBS与SEBS-Urban模型对研究区2015—2017年ET计算结果 mm

    Table  1.   Simulation ET results of SEBS vs. SEBS-Urban for study area over 2015—2017

    月份 2015年 2016年 2017年
    SEBS SEBS-Urban SEBS SEBS-Urban SEBS SEBS-Urban
    1月 0.3 3.0 1.5 2.9 -1.5 0.4
    2月 11.8 14.6 16.7 20.0 9.8 11.3
    3月 35.7 37.5 31.9 34.5 27.5 29.9
    4月 55.6 57.5 51.3 53.0 49.9 51.3
    5月 83.1 86.2 60.6 62.2 85.2 88.2
    6月 75.7 79.2 79.3 83.8 83.1 86.8
    7月 80.1 86.0 77.4 82.2 88.2 95.0
    8月 92.0 98.3 79.6 85.2 100.2 105.4
    9月 51.8 56.7 46.5 50.4 61.4 66.3
    10月 25.3 28.3 12.3 14.5 8.9 11.1
    11月 -14.3 -12.3 -0.9 2.4 6.2 8.4
    12月 -8.6 -5.5 -8.6 -4.7 -1.4 1.0
    总量 488.5 529.6 447.4 486.4 517.7 555.1
    下载: 导出CSV
  • [1] ZHAO J F, LI C, YANG T Y, et al. Estimation of high spatiotemporal resolution actual evapotranspiration by combining the SWH model with the METRIC model[J]. Journal of Hydrology, 2020, 586: 124883. doi:  10.1016/j.jhydrol.2020.124883
    [2] 宋晓猛, 张建云, 贺瑞敏, 等. 北京城市洪涝问题与成因分析[J]. 水科学进展, 2019, 30(2): 153-165. https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201902001.htm

    SONG X M, ZHANG J Y, HE R M, et al. Urban flood and waterlogging and causes analysis in Beijing[J]. Advances in Water Science, 2019, 30(2): 153-165. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201902001.htm
    [3] 陈似蓝, 刘家宏, 王浩. 城市水资源需求场理论及应用初探[J]. 科学通报, 2016, 61(13): 1428-1435. https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201613007.htm

    CHEN S L, LIU J H, WANG H. Initial research on the theory and application of urban water demand field[J]. Chinese Science Bulletin, 2016, 61(13): 1428-1435. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201613007.htm
    [4] GAO X, SUN M, LUAN Q, et al. The spatial and temporal evolution of the actual evapotranspiration based on the remote sensing method in the Loess Plateau[J]. The Science of the Total Environment, 2020, 708: 135111. doi:  10.1016/j.scitotenv.2019.135111
    [5] LI F W, CAO R X, ZHAO Y, et al. Remote sensing Penman-Monteith model to estimate catchment evapotranspiration considering the vegetation diversity[J]. Theoretical and Applied Climatology, 2017, 127(1/2): 111-121.
    [6] KUMAR V, KHAN S, SHAKTI B, et al. Review of evapotranspiration methodologies[J]. American Journal of Earth Science and Engineering, 2018, 1(2): 72-84.
    [7] 张圣微, 张鹏, 张睿, 等. 科尔沁沙地典型区生长季蒸散发估算及其变化特征[J]. 水科学进展, 2018, 29(6): 768-778. https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201806002.htm

    ZHANG S W, ZHANG P, ZHANG R, et al. Estimation of growing season evapotranspiration and its variation in a typical area of Horqin Sandy Land[J]. Advances in Water Science, 2018, 29(6): 768-778. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201806002.htm
    [8] LIU X Y, XU J Z, ZHOU X Y, et al. Evaporative fraction and its application in estimating daily evapotranspiration of water-saving irrigated rice field[J]. Journal of Hydrology, 2020, 584: 124317. doi:  10.1016/j.jhydrol.2019.124317
    [9] 刘曼晴, 胡德勇, 于琛, 等. 辽河三角洲湿地生长季蒸散量时空格局及影响因素[J]. 生态学报, 2020, 40(2): 701-710. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB202002029.htm

    LIU M Q, HU D Y, YU C, et al. Temporal and spatial change characteristics of growing season evapotranspiration and its cause analysis in Liaohe River delta wetland, China[J]. Acta Ecologica Sinica, 2020, 40(2): 701-710. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-STXB202002029.htm
    [10] 周琳. 北京市城市蒸散发研究[D]. 北京: 清华大学, 2015.

    ZHOU L. Urban evapotranspiration in Beijing[D]. Beijing: Tsinghua University, 2015. (in Chinese)
    [11] NIE W S, ZAITCHIK B F, NI G H, et al. Impacts of anthropogenic heat on summertime rainfall in Beijing[J]. Journal of Hydrometeorology, 2017, 18(3): 693-712. doi:  10.1175/JHM-D-16-0173.1
    [12] SAILOR D J. A review of methods for estimating anthropogenic heat and moisture emissions in the urban environment[J]. International Journal of Climatology, 2011, 31(2): 189-199. doi:  10.1002/joc.2106
    [13] CONG Z T, SHEN Q N, ZHOU L, et al. Evapotranspiration estimation considering anthropogenic heat based on remote sensing in urban area[J]. Science China Earth Sciences, 2017, 60(4): 659-671. doi:  10.1007/s11430-016-0216-3
    [14] FARIDATUL M I, WU B, ZHU X L, et al. Improving remote sensing based evapotranspiration modelling in a heterogeneous urban environment[J]. Journal of Hydrology, 2020, 581: 124405. doi:  10.1016/j.jhydrol.2019.124405
    [15] JÄRVI L, GRIMMOND C S B, CHRISTEN A. The surface urban energy and water balance scheme (SUEWS): evaluation in Los Angeles and Vancouver[J]. Journal of Hydrology, 2011, 411(3/4): 219-237. http://www.sciencedirect.com/science/article/pii/S0169204611002362
    [16] 周晋军, 刘家宏, 董庆珊, 等. 城市耗水计算模型[J]. 水科学进展, 2017, 28(2): 276-284. https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201702013.htm

    ZHOU J J, LIU J H, DONG Q S, et al. Simulation model for urban water dissipation[J]. Advances in Water Science, 2017, 28(2): 276-284. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201702013.htm
    [17] ZHOU J J, LIU J H, YAN D Y, et al. Dissipation of water in urban area, mechanism and modelling with the consideration of anthropogenic impacts: a case study in Xiamen[J]. Journal of Hydrology, 2019, 570: 356-365. doi:  10.1016/j.jhydrol.2018.12.054
    [18] 郭生练, 郭家力, 侯雨坤, 等. 基于Budyko假设预测长江流域未来径流量变化[J]. 水科学进展, 2015, 26(2): 151-160. https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201502001.htm

    GUO S L, GUO J L, HOU Y K, et al. Prediction of future runoff change based on Budyko hypothesis in Yangtze River basin[J]. Advances in Water Science, 2015, 26(2): 151-160. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SKXJ201502001.htm
    [19] YU X, YANG H B, LI S E, et al. An improved conceptual model quantifying the effect of climate change and anthropogenic activities on vegetation change in arid regions[J]. Remote Sensing, 2019, 11(18): 2110. doi:  10.3390/rs11182110
    [20] YANG W J, WANG Y B, LIU X, et al. Evaluation of the rescaled complementary principle in the estimation of evaporation on the Tibetan Plateau[J]. Science of the Total Environment, 2020, 699: 134367. doi:  10.1016/j.scitotenv.2019.134367
    [21] DU C, SUN F, YU J, et al. New interpretation of the role of water balance in an extended Budyko hypothesis in arid regions[J]. Hydrology and Earth System Sciences, 2016, 20(1): 393-409. doi:  10.5194/hess-20-393-2016
    [22] HAN S J, HU H P, YANG D W, et al. Irrigation impact on annual water balance of the oases in Tarim Basin, Northwest China[J]. Hydrological Processes, 2011, 25(2): 167-174. doi:  10.1002/hyp.7830
    [23] CHEN H, HUO Z L, ZHANG L, et al. New perspective about application of extended Budyko formula in arid irrigation district with shallow groundwater[J]. Journal of Hydrology, 2020, 582: 124496. doi:  10.1016/j.jhydrol.2019.124496
    [24] 杨汉波. 流域水热耦合平衡方程推导及其应用[D]. 北京: 清华大学, 2008.

    YANG H B. Derivation and application of the coupled water-energy balance equation[D]. Beijing: Tsinghua University, 2008. (in Chinese)
    [25] WEI H B, ZHANG J E, SINGH S K, et al. Response of hydrological system to urbanization: a case study in Tianjin City, China[J]. Desalination and Water Treatment, 2019, 149: 305-314. doi:  10.5004/dwt.2019.23847
    [26] ZHANG S, HUANG G, QI Y J, et al. Impact of urbanization on summer rainfall in Beijing-Tianjin-Hebei metropolis under different climate backgrounds[J]. Theoretical and Applied Climatology, 2018, 133(3): 1093-1106. doi:  10.1007/s00704-017-2225-3
    [27] SU Z. The surface energy balance system (SEBS) for estimation of turbulent heat fluxes[J]. Hydrology and Earth System Sciences, 2002, 6(1): 85-99. doi:  10.5194/hess-6-85-2002
    [28] NIE W S, SUN T, NI G H. Spatiotemporal characteristics of anthropogenic heat in an urban environment: a case study of Tsinghua Campus[J]. Building and Environment, 2014, 82: 675-686. doi:  10.1016/j.buildenv.2014.10.011
    [29] 舒松, 余柏蒗, 吴健平, 等. 基于夜间灯光数据的城市建成区提取方法评价与应用[J]. 遥感技术与应用, 2011, 26(2): 169-176. https://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201102007.htm

    SHU S, YU B L, WU J P, et al. Methods for deriving urban built-up area using night-light data: assessment and application[J]. Remote Sensing Technology and Application, 2011, 26(2): 169-176. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YGJS201102007.htm
    [30] KARSSENBERG D, SCHMITZ O, SALAMON P, et al. A software framework for construction of process-based stochastic spatio-temporal models and data assimilation[J]. Environmental Modelling & Software, 2010, 25(4): 489-502. http://www.sciencedirect.com/science/article/pii/S1364815209002643
    [31] LANDERER F W, SWENSON S C. Accuracy of scaled GRACE terrestrial water storage estimates[J]. Water Resources Research, 2012, 48(4): W04531. http://gji.oxfordjournals.org/external-ref?access_num=10.1029/2011WR011453&link_type=DOI
    [32] 杨汉波, 杨大文, 雷志栋, 等. 任意时间尺度上的流域水热耦合平衡方程的推导及验证[J]. 水利学报, 2008, 39(5): 610-617. doi:  10.3321/j.issn:0559-9350.2008.05.015

    YANG H B, YANG D W, LEI Z D, et al. Derivation and validation of watershed coupled water-energy balance equation at arbitrary time scale[J]. Journal of Hydraulic Engineering, 2008, 39(5): 610-617. (in Chinese) doi:  10.3321/j.issn:0559-9350.2008.05.015
    [33] 姜艳阳, 王文, 周正昊. MODIS MOD16蒸散发产品在中国流域的质量评估[J]. 自然资源学报, 2017, 32(3): 517-528. https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZX201703014.htm

    JIANG Y Y, WANG W, ZHOU Z H. Evaluation of MODIS MOD16 evapotranspiration product in Chinese river basins[J]. Journal of Natural Resources, 2017, 32(3): 517-528. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZRZX201703014.htm
    [34] JAHANGIR M H, ARAST M. Remote sensing products for predicting actual evapotranspiration and water stress footprints under different land cover[J]. Journal of Cleaner Production, 2020, 266: 121818. doi:  10.1016/j.jclepro.2020.121818
    [35] 叶红, 张廷斌, 易桂花, 等. 2000-2014年黄河源区ET时空特征及其与气候因子关系[J]. 地理学报, 2018, 73(11): 2117-2134. doi:  10.11821/dlxb201811006

    YE H, ZHANG T B, YI G H, et al. Spatio-temporal characteristics of evapotranspiration and its relationship with climate factors in the source region of the Yellow River from 2000 to 2014[J]. Acta Geographica Sinica, 2018, 73(11): 2117-2134. (in Chinese) doi:  10.11821/dlxb201811006
  • [1] 孙青言, 陆垂裕, 郭辉, 严聆嘉, 何鑫, 吴初.  三江平原土地利用变化对水量平衡的影响 . 水科学进展, 2021, 32(5): 694-706. doi: 10.14042/j.cnki.32.1309.2021.05.005
    [2] 任梅芳, 徐宗学, 庞博.  变化环境下城市洪水演变驱动机理——以北京市温榆河为例 . 水科学进展, 2021, 32(3): 345-355. doi: 10.14042/j.cnki.32.1309.2021.03.003
    [3] 徐宗学, 陈浩, 任梅芳, 程涛.  中国城市洪涝致灾机理与风险评估研究进展 . 水科学进展, 2020, 31(5): 713-724. doi: 10.14042/j.cnki.32.1309.2020.05.008
    [4] 宋晓猛, 张建云, 贺瑞敏, 邹贤菊, 张春桦.  北京城市洪涝问题与成因分析 . 水科学进展, 2019, 30(2): 153-165. doi: 10.14042/j.cnki.32.1309.2019.02.001
    [5] 张圣微, 张鹏, 张睿, 刘廷玺.  科尔沁沙地典型区生长季蒸散发估算及其变化特征 . 水科学进展, 2018, 29(6): 768-778. doi: 10.14042/j.cnki.32.1309.2018.06.002
    [6] 王卫光, 李进兴, 魏建德, 邵全喜, 邓超, 余钟波.  基于蒸散发数据同化的径流过程模拟 . 水科学进展, 2018, 29(2): 159-168. doi: 10.14042/j.cnki.32.1309.2018.02.002
    [7] 胡庆芳, 张建云, 王银堂, 黄勇, 刘勇, 李伶杰.  城市化对降水影响的研究综述 . 水科学进展, 2018, 29(1): 138-150. doi: 10.14042/j.cnki.32.1309.2018.01.016
    [8] 周晋军, 刘家宏, 董庆珊, 王忠静, 王浩, 向晨瑶, 栾勇.  城市耗水计算模型 . 水科学进展, 2017, 28(2): 276-284. doi: 10.14042/j.cnki.32.1309.2017.02.013
    [9] 曾燕, 邱新法, 刘昌明.  黄河流域蒸散量分布式模拟 . 水科学进展, 2014, 25(5): 632-640.
    [10] 孙艳伟, 王文川, 魏晓妹, POMEROY C A, 王富强.  城市化生态水文效应 . 水科学进展, 2012, 23(4): 569-574. doi: CNKI:32.1309.P.20120614.2251.013
    [11] 赵军, 单福征, 杨凯, 吴阿娜.  平原河网地区河流曲度及城市化响应 . 水科学进展, 2011, 22(5): 631-637.
    [12] 田辉, 文军, 马耀明, 王介民, 吕世华, 张堂堂, 孙方林, 刘蓉.  夏季黑河流域蒸散发量卫星遥感估算研究 . 水科学进展, 2009, 20(1): 18-24.
    [13] 陈莹莹, 施建成, 杜今阳, 蒋玲梅.  基于GLDAS的中国区地表能量平衡数值试验 . 水科学进展, 2009, 20(1): 25-31.
    [14] 吴炳方, 熊隽, 闫娜娜, 杨雷东, 杜鑫.  基于遥感的区域蒸散量监测方法——ETWatch . 水科学进展, 2008, 19(5): 671-678.
    [15] 孙丽, 宋长春.  三江平原典型沼泽湿地能量平衡和蒸散发研究 . 水科学进展, 2008, 19(1): 43-48.
    [16] 许士国, 王昊.  测量芦苇沼泽蒸散发量的渗流补偿方法 . 水科学进展, 2007, 18(4): 496-503.
    [17] 姜杰, 张永强, 刘昌明.  运用高分辨率ETM+数据计算区域作物蒸散 . 水科学进展, 2005, 16(2): 274-279.
    [18] 王新平, 康尔泗, 张景光, 李新荣.  草原化荒漠带人工固沙植丛区土壤水分动态 . 水科学进展, 2004, 15(2): 216-222.
    [19] 庞治国, 付俊娥, 李纪人, 肖乾广.  基于能量平衡的蒸散发遥感反演模型研究 . 水科学进展, 2004, 15(3): 364-369.
    [20] 陈鸣, 潘之棣.  用卫星遥感热红外数据估算大面积蒸散量 . 水科学进展, 1994, 5(2): 126-133.
  • 加载中
图(5) / 表 (1)
计量
  • 文章访问数:  396
  • HTML全文浏览量:  110
  • PDF下载量:  87
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-30
  • 网络出版日期:  2021-01-06
  • 刊出日期:  2021-05-30

天津平原区城市化地区蒸散发特征

doi: 10.14042/j.cnki.32.1309.2021.03.005
    基金项目:

    国家自然科学基金资助项目 51779165

    国家自然科学基金资助项目 51879181

    作者简介:

    曹润祥(1990—), 男, 河南新乡人, 博士研究生, 主要从事水文学及水资源研究。E-mail: runx.cao@foxmail.com

    通讯作者: 李发文, E-mail: lifawen@tju.edu.cn
  • 中图分类号: P339;TV11

摘要: 为揭示城市蒸散发特征,采用考虑人为热影响的SEBS-Urban模型计算天津平原区2015—2017年逐月蒸散量。使用两步验证法对结果进行检验,首先利用MODIS MOD16产品对非建成区模拟结果合理性进行验证,其次利用城市耗水(UWD)模型得到的天津大学卫津路校区月蒸散量对建成区模拟结果验证。在SEBS-Urban模拟结果基础上,检验Budyko方程在城市化地区的适用性。结果表明:① SEBS-Urban模型精度可靠,体现了蒸散发项在城市地表能量平衡与二元水循环中的相合性;②人为热或社会侧的耗水会导致蒸散发量增加,建成区增幅达85%~115%,整个研究区增幅为7.2%~8.7%;③ Budyko方程在研究区具有较好的拟合效果,可应用于城市地区。本研究成果可为城市化地区蒸散发研究提供必要的参考。

English Abstract

曹润祥, 李发文, 李建柱, 冯平. 天津平原区城市化地区蒸散发特征[J]. 水科学进展, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
引用本文: 曹润祥, 李发文, 李建柱, 冯平. 天津平原区城市化地区蒸散发特征[J]. 水科学进展, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
CAO Runxiang, LI Fawen, LI Jianzhu, FENG Ping. Exploring evapotranspiration characteristics in urbanized areas of Tianjin plain[J]. Advances in Water Science, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
Citation: CAO Runxiang, LI Fawen, LI Jianzhu, FENG Ping. Exploring evapotranspiration characteristics in urbanized areas of Tianjin plain[J]. Advances in Water Science, 2021, 32(3): 366-375. doi: 10.14042/j.cnki.32.1309.2021.03.005
  • 蒸散发是水文循环和能量循环的关键环节[1], 研究蒸散发对深入认识水文循环机理、正确评估和科学管理水资源具有十分重要的意义。城市蒸散发是区域蒸散发中的特殊形式。在高强度的人类活动和气候变化下, 城市水循环表现出典型的“自然-社会”二元特性[2], 城市的社会侧水循环过程可概化为“供水-用水-耗水-排水”[3]。城市耗水可以理解为城市区域所有可能的蒸散发, 包括自然下垫面蒸散发和不透水面的蒸发、建筑内部人类生产生活活动的蒸散发等。城市蒸散发对城市微气象、城市生态环境及城市的供水都有着显著影响。

    在已往的研究中, 许多方法被用来计算蒸散发量, 包括水量平衡法[4]、气象学方法[5-6]与地表能量平衡法[7]。传统的方法一般仅考虑均匀下垫面, 对城市复杂下垫面的理解不足, 结合遥感反演方法可弥补此方面的不足[8]。目前基于遥感的蒸散发计算方法如SEBS模型(Surface Energy Balance System)和SEBAL模型(Surface Energy Balance Algorithm for Land), 往往得出城市下垫面蒸散发量比森林、农田下垫面小的结论[9], 但这无法解释城市高耗水的现象[10]。从能量的角度分析, 传统基于遥感的地表能量平衡法计算区域蒸散时, 能量平衡方程中的人为热项和水平能量传输项常常被忽略。而在城市环境中, 密集的人口导致人为热大量排放[11], 使其在地表能量收支平衡中占有不容忽视的比重[12]。Cong等[13]使用基于SEBS模型改进的SEBS-Urban模型计算了北京地区的蒸散量, 模型考虑了城市人为热对模拟结果的影响, 显示人为热对城市蒸散发具有显著的影响作用。Faridatul等[14]提出了改进的适用于城市的能量平衡算法uSEBAL, 发现考虑了城市下垫面类型与人为热通量后, 计算结果优于传统SEBAL模型。从二元水循环的角度分析, 城市蒸散不仅包括自然侧的草地、树木灌丛以及其他植被覆盖地面的蒸散[15], 也包括人类活动下的室内环境中产生的蒸发[16-17], 但这部分蒸发在传统方法计算过程中没有很好体现。为涵盖城市各个环节的耗水, Zhou等[17]提出了城市耗水模型(UWD)计算框架, 并据此对厦门地区的蒸散进行了估算, 并通过区域水量平衡的方法对结果进行验证, 证明两者对应较好且数值上明显高于传统的遥感反演结果。

    SEBS-Urban及UWD模型在城市化地区的蒸散发模拟对数据有较高的要求。基于水-热耦合平衡理论的Budyko方程不需要大量的参数和计算资源, 广泛应用于流域水文过程的研究[18]。已往的研究证明, 不同Budyko模型中水热耦合控制参数受到流域下垫面特征的影响, 包括植被覆盖、土壤类型、地型要素及季节性气候指数[19-20]。很多研究也探讨了非闭合流域水量平衡关系与拓展Budyko框架假设(Extended Budyko Hypothesis)的可解释性[21], 且一些研究也证明了Budyko在灌区尺度的适用性[22-24], 但较少有研究针对城市化地区探讨其适用性。

    本研究以天津平原区为例, 用SEBS-Urban方法计算区域蒸散发, 分别在非城镇区和城市建成区内用MOD16产品和UWD对比分析, 验证模型的合理性, 以期初步揭示城市内特殊的水-热平衡关系;分析Budyko方程在城市地区的适用性, 评估现阶段研究区人为热排放对区域水热平衡关系的影响。

    • 本研究以天津市平原区为研究区域。天津是中国经济最发达也是城市化进程最为迅速的城市之一, 中心城区高强度的人类经济活动对区域的水文循环过程及局地气候产生了显著的影响[25-26]。天津市年内温差较大, 多年平均气温在12 ℃左右, 极端最高气温为42.7 ℃, 极端最低气温为-27.4 ℃;年平均相对湿度在50%~70%;区域多年平均降水在534.5~649.9 mm之间, 属半湿润半干旱地带, 年平均陆面蒸散发512.3 mm, 水面蒸发961.8 mm。研究过程中以天津大学卫津路校区为城市建成区代表特征点。天津大学卫津路校区地属天津市南开区, 位于天津市主城区西南部。研究区位置示意如图 1所示。

      图  1  研究区位置示意

      Figure 1.  Location of the study area

    • SEBS-Urban模型的数据输入包括遥感数据和气象数据两部分。本研究中, 地表反照率、地面辐射率、植被覆盖度、叶面积指数(LAI)、归一化植被参数(NDVI)以及地表温度通过MODIS遥感数据反演得到;地面参考高度处风速、大气压、相对湿度、比湿、地表参考高度处气温、向下短波辐射、向下长波辐射等输入值通过中国区域地面气象要素驱动数据集[24]获得。此外, 本研究为估算城区人为热使用了DMSP/OLS的夜间灯光V4数据, 资料来源于美国地球物理数据中心, 空间分辨率为1 km×1 km。

    • 计算天津大学卫津路校区位置城市耗水模型UWD所用到的气象数据也采用中国区域地面气象要素驱动数据集。校区2017年的土地利用及逐月供用水统计数据由天津大学后勤保障部提供。

    • 为构建研究区Budyko方程, 除上述数据集中的降水数据之外, 还需土壤湿度数据与区域入境水量数据。其中, 区域土壤湿度数据的获取采用GLDAS-Noah陆面模式的V2.1数据, 数据长度为2015—2017年, 时间分辨率为1个月, 空间分辨率为0.25°×0.25°。天津地区的入境水量由天津市水务局提供, 时间精度为1个月。

    • SEBS模型是由Su于2002年提出的蒸散单层模型[27]。SEBS模型基于能量平衡理论, 增加了地表粗糙度的计算拓展模块以估算显热传输。在未考虑人为热和水平能量传输的情况下, 地表能量平衡可用式(1)表示:

      $$ {R_{\rm{n}}} = H + {G_0} + \lambda {E_{\rm{T}}} $$ (1)

      式中: Rn为净辐射, W/m2G0为土壤通量, W/m2H为显热通量, W/m2λET为潜热通量, W/m2, 其中λ为水汽化热, ET为实际蒸散发。如果求得RnG0H的值, 便可以求得余项潜热λET, 进而可知实际蒸散发ET的值。

      高度密集的城市建筑区的人为热可高达400 W/m2[28], 若在能量平衡项中忽略此影响, 势必会对城市地区蒸散量的估计造成较大的偏差。因此, 本文采用Cong等[13]改进的SEBS-Urban模型对天津市城区的蒸散量进行计算。SEBS-Urban主要是将能量平衡方程的左侧添加了人为热项Qf, 并用区域夜间灯光指数估算人为热。其中地表能量方程可写为

      $$ {R_{\rm{n}}} + {Q_{\rm{f}}} = {G_0} + H + \lambda {E_{\rm{T}}} $$ (2)

      本研究中, 人为热影响区域的灯光值域设为52~63[29], 对人为热的赋值情况为: 冬夏季值域为50~75 W/m2, 春秋季值域为30~50 W/m2[13, 28], 其中夏冬人为热较高的原因是空调和暖气的使用。本研究中SEBS-Urban的模型实现基于开源地理及水文建模工具PCRaster[29]

    • 相比基于能量平衡的SEBS-Urban模型, UWD模型[17]从城市二元水循环的角度解析城市蒸散发: 城市环境中建筑的窗、门及其他通风口类似自然环境中的植被叶片气孔, 水汽通过这些通道逸出而重新进入水循环。依据UWD模型对天津大学卫津路校区的建筑体、硬化地面与屋顶、水面及绿地下垫面的耗水进行分项计算。

      将建筑依宿舍区、教学办公区与食堂分为3类, 对于不同类型建筑内部耗水可按式(3)计算:

      $$ {E_{\rm{B}}} = {S_{\rm{B}}}\varphi + 0.03{A_{\rm{G}}}\delta {D_{\rm{f}}} $$ (3)
      $$ \varphi = \sum\limits_{i = 1}^n {{\kappa _i}{\xi _i}} $$ (4)

      式中: EB为建筑内部耗水量, m3SB为建筑体内各月的用水量, m3φ为不同类型建筑的耗水比例;AG为室内面积;δ为湿地面比例系数;Df为室内地面的耗水强度;κi为各种用水项的比例;ξi为不同用水项的耗水比例;n为6, 分别为洗漱、做饭、洗澡、清洁、饮用与洗衣6类。

      对于发生于硬化路面与屋顶的蒸散(EH)可表示为

      $$ {E_{\rm{H}}} = {P_0} + \left( {1 - \psi } \right)({P_{\rm{Y}}} - {P_0}) + {W_{\rm{A}}} $$ (5)

      式中: P0为各月未产流的日降水的总和, mm;PY为各月份降水量总和, mm;WA为各月校园内部道路的人为洒水量, mm。

      对于水面蒸发(Ew)及绿地蒸散发(EG)分别采用Penman和Penman-Monteith公式计算。

    • 目前大多数城市地区难以获得可靠的地面通量观测数据, 对城市地区蒸散发模拟结果进行精确验证存在着客观上的困难[14]。有关城市地区蒸散发模型验证一般常采用与水量平衡法[13, 16-17]或蒸发皿数据对比的方法[14]。这些方法可以从整体上说明模型的合理性, 但难以反映改进模型在建成区的细节效果, 且当研究区内城市面积占比较小时, 水量平衡计算过程中区域蓄水量的不确定性[31]容易掩盖城市化引起的蒸散发响应。

      本研究采用了两步验证方法: 整个研究区内非建设用地占据更大比重(约85%), 首先采用MOD16数据对非城镇地区验证, 从而可从整体上把握模型的合理性;在城市化地区, UWD模型原理简洁, 模型结构的不确定性较小, 可靠性较高, 本文选用UWD模型结果进行验证, 以反映出SEBS-Urban模型在考虑人为热前后的蒸散差异的细节。

    • 杨汉波等[32]基于Budyko曲线, 将任意时段内的实际蒸散量(E)表示为该时段内的可用水量(Q)和可用能量(采用潜在蒸发量E0表示)的函数。其中时段内Q包括该时段的降雨量(P)、外部供水量(I)与时段初始时刻的土壤可供水量(S)。

      $$ E = \frac{{{E_0}\left( {P + I + S} \right)}}{{{{({{\left( {P + I + S} \right)}^m} + E_0^m)}^{1/m}}}} $$ (6)

      式中: m为水热耦合参数。

    • 在非城镇区, 模型无人为热项添加, SEBS-Urban退化为普通SEBS模型, 采用MOD16数据对非城镇地区模拟结果做验证。MOD16是目前应用最广的遥感蒸散发产品之一[33]。该数据集已经通过了全球232个流域和46个野外涡度通量站台的检验, 总体模拟精度达到86%[34]。MOD16数据集在国内被应用到多个区域的蒸散发时空变化研究中[35], 且研究[33]表明, MOD16产品在外流区与水量平衡法的一致性高于内流区, 在北方优于南方。因此, 可认为用MOD16产品来检验SEBS在非城市区域的表现有较高的可靠性。MOD16数据获取算法基于Penman-Monteith方程, 该产品不能获得建成区的蒸散发数据, 因此仅被用来检验非城镇地区SEBS/SEBS-Urban模型的准确性[31]。从2017年MOD16 A2产品中筛选了18景质量较好的产品, 并从中随机抽取30个具有连续有效值的非城区格点和对应位置处SEBS/SEBS-Urban的计算结果进行对比, 如图 2所示。从相关系数来看, 30个点中除了点22和点26外, 相关系数整体介于0.5和0.8之间。从SEBS/SEBS-Urban与MOD16 A2的标准差比值(σ)来看, 仅有5个点σ < 1, 剩余25个点σ>1, 表明SEBS-Urban计算结果的年内变化幅度多数大于MOD16 A2产品, 但相差幅度不大(除点14 σ>2.0外, 其他点均小于1.8)。从均方根偏差(RMSD)来看, RMSD值均不超过2.0倍的MOD16 A2标准差σMOD16。综上, SEBS/SEBS-Urban模型可以有效地模拟非城镇地区的蒸散发年内变化过程, 模型在整体上具有合理性。

      图  2  SEBS反演结果相对于MOD16 A2产品的泰勒图

      Figure 2.  Taylor diagram of SEBS estimated ET vs. MOD16 A2 product

    • 将SEBS-Urban计算的城市部分ET与典型地市区域(天津大学卫津路校区)的UWD作对比验证, 以反映模型添加人为热前后在建成区ET的细节效果。SEBS和SEBS-Urban模型计算结果与UWD计算结果对比如图 3所示。

      图  3  UWD、SEBS与SEBS-Urban结果

      Figure 3.  UWD, SEBS and SEBS-Urban simulation results

      图 3可以看出, 在未考虑人为热时, SEBS模型与UWD计算结果有较大的偏差, SEBS计算的ET全年均低于UWD计算结果。增加人为热之后, SEBS-Urban模型计算的ET结果整体较SEBS模型显著增加, 除7月、8月外, SEBS-Urban与UWD的计算结果表现出很好的一致性。模型7月、8月的较大差异可能是SEBS-Urban的蒸发比插值引起对ET的高估。SEBS-Urban利用遥感数据仅可反演卫星过境日的当日ET值, 但遥感模型常常因为天气状况和遥感质量的原因无法获取清晰的图像而造成数据缺失, 1年中实际可用的遥感数据最多只有80 d左右。而由于区域气候特征, 天津地区7月、8月雨水充沛, 2017年7月、8月每月可用的遥感数据仅有1幅(晴日), 因此, 据此插值不可避免地造成了对ET的高估。总体看来, SEBS-Urban模型在考虑人为热影响后, 与UWD计算结果的相关系数达0.84, Nash效率系数为0.82, 两者表现出了较好的一致性。

      UWD计算的年总蒸散发量为422.6 mm, ET最大的月份出现在6—8月, 由图 3(b)可以看出, 夏季前后社会侧耗水, 包括建筑内(宿舍、教学区、食堂)蒸发、硬化路面洒水蒸发, 对整体蒸散发的贡献比最小, 而冬季则为最大。这是因为冬季天然蒸散发很小, 相对稳定的室内人类用水活动造成水量蒸散发占据了较大的比例;而夏季校园内植物蒸腾量极大增加, 成为总蒸散发量的主导因素。

      上述结果反映了城市地区特殊水-热关系。从能量的角度来看, 一方面城市由于不透水面积的增加导致了显热的增加, 另一方面, 城市人为热贡献了额外的潜热增长。而从二元水循环的角度来看, 一方面城市不透水面的增加导致了自然侧的蒸散发减少;另一方面, 城市中人类活动用水形成了社会侧的额外蒸散发形式。无论从哪种角度分析城市地区蒸散, 忽略人类活动影响都将无法准确描述区域蒸散发特征。

    • 利用SEBS模型和SEBS-Urban模型获得研究区2015—2017年的蒸散发反演结果如表 1所示, 其空间分布如图 4所示。由表 1可以看出, 对于整个天津平原区, 由于考虑了人为热, SEBS-Urban各个月份的ET计算结果较SEBS均有不同程度的增高, 3 a增幅范围为7.2%~8.7%。表 1中冬季部分月份结果为负, 是由于出现净辐射量为负的情形, 对应时段可能出现了较大强度的凝结水[10]。从ET年内时程分配来看, 两模型计算结果年内变化趋势具有明显的一致性, 各年相关系数均达到0.99以上。从ET的空间分布来看, SEBS计算结果的ET低值区恰好与高度城市化地区相吻合。在以往蒸散模拟研究中, 这一特征常被用来解释城市地区的陆面硬化对蒸散发的抑制作用;而计算结果中ET高值集中的地方出现在水域较多的沿海和湿地较多的地区, 这些地方水分供给较好, 蒸散发相对剧烈。SEBS-Urban模型结果中, 由于考虑人为热的影响, 城区的模拟结果较SEBS显著增高, 建成区3 a增幅范围为85%~115%。中心城区考虑了最高的人为热, 但ET仍较周边略微偏小, 这反映了城区因硬化路面导致的ET减小和人为热剧烈导致潜热增加两者互相作用的结果, 综合作用的结果仍然降低了城区的ET;郊区城镇化区域在考虑人为热后, ET显著增加, 这可能是因为郊区聚集了较多的人口(表现为较高的夜间灯光指数), 但对下垫面的改造远不如城市中心剧烈。对于非城市化地区(夜间灯光指数 < 52), SEBS模型与SEBS-Urban具有相同的模拟结果。

      表 1  SEBS与SEBS-Urban模型对研究区2015—2017年ET计算结果 mm

      Table 1.  Simulation ET results of SEBS vs. SEBS-Urban for study area over 2015—2017

      月份 2015年 2016年 2017年
      SEBS SEBS-Urban SEBS SEBS-Urban SEBS SEBS-Urban
      1月 0.3 3.0 1.5 2.9 -1.5 0.4
      2月 11.8 14.6 16.7 20.0 9.8 11.3
      3月 35.7 37.5 31.9 34.5 27.5 29.9
      4月 55.6 57.5 51.3 53.0 49.9 51.3
      5月 83.1 86.2 60.6 62.2 85.2 88.2
      6月 75.7 79.2 79.3 83.8 83.1 86.8
      7月 80.1 86.0 77.4 82.2 88.2 95.0
      8月 92.0 98.3 79.6 85.2 100.2 105.4
      9月 51.8 56.7 46.5 50.4 61.4 66.3
      10月 25.3 28.3 12.3 14.5 8.9 11.1
      11月 -14.3 -12.3 -0.9 2.4 6.2 8.4
      12月 -8.6 -5.5 -8.6 -4.7 -1.4 1.0
      总量 488.5 529.6 447.4 486.4 517.7 555.1

      图  4  2015—2017年SEBS与SEBS-Urban模型对天津平原地区蒸散发分布模拟结果

      Figure 4.  Annual evapotranspiration distribution modelled by SEBS and SEBS-Urban over 2015—2017

    • 以研究区SEBS-Urban模型2015—2017年计算结果为样本进行Budyko方程拟合。由于序列较短, Budyko参数的验证借鉴了“留一法”交叉验证方法的思路, 即每次取1 a对模拟结果作验证。由于人为热在1 a中不同时段具有一定的变化, 因此在对水热耦合参数m率定时, 将m看成Hj(j=1, 2, …, 12)向量, 以率定期RMSE值最小为优化目标, 并用和声算法(Harmony Search)对其优化求解。水热耦合方程模拟得到天津地区2015—2017年逐月的蒸散发量, 结果如图 5(a)所示, 可以看出3种验证情形下, Budyko都可大致反映出月尺度下研究区蒸散规律。由图 5(b)可以看出, 验证期内两者保持了很好的线性关系, 两者线性回归的方程斜率为0.8~1.2。同时, 相关系数在0.9以上, p值均小于0.001, 表明验证期模拟的月蒸散过程与SEBS-Urban吻合较好, 在一定程度上可以说明水热耦合平衡方程仍可以适用于模拟受人为热影响剧烈的城市地区的月蒸散量。利用天津市2015—2017年数据拟合得到参数m的年内变化过程, 如图 5(c)所示。参数m反映了区域下垫面特征, 对于天津地区, 城区多样的景观植被、城市周边的作物及人为热的排放共同影响着区域水热平衡关系。由图 5(c)可以看出, 研究区Budyko参数m年内变化最主要体现的是季节变化趋势。用研究区SEBS-Urban与SEBS的各月差值表征人为热影响, 同参数m作相关分析, 并未发现显著相关(p>0.05)。将SEBS模型结果拟合得到的Budyko参数向量与SEBS-Urban的Budyko参数向量做t检验, 亦未发现显著差异。这可能是因为在整个研究区范围, 人为热影响的区域比例仅占研究区15%, 虽然在城市建成区内人为热对蒸散发影响明显, 但上升到整个区域尺度, 这种影响未能在Budyko参数上有所体现。在更小的空间尺度上(如仅建成区)的研究有待进一步探索。

      图  5  2015—2017年Budyko模拟结果分析

      Figure 5.  Budyko monthly simulation results analysis over 2015—2017

    • 本研究以天津平原区为研究区, 建立了考虑人为热影响的SEBS-Urban蒸散发模型, 采用两步验证法对SEBS-Urban模型的可靠性进行了验证;利用计算得出的研究区2015—2017年蒸散发结果探讨了Budyko方程在城市地区的适用性。主要得到以下结论:

      (1) 通过与MOD16产品和天津大学校区城市耗水模型计算结果对比分析, 证明SEBS-Urban模拟结果在研究区内具有合理性, 且能反映建成区考虑人为热后的蒸散发变化细节, 体现了蒸散发项在城市地表能量平衡与二元水循环系统中的相合性。

      (2) 考虑人为热后, 2015—2017年天津平原区总蒸散量各年增幅为7.2%~8.7%, 而建成区增幅达85%~115%。

      (3) Budyko方程对SEBS-Urban结果的拟合效果较好; 但研究区尺度上, 水热耦合参数未与人为热表现出显著的相关性。

参考文献 (35)

目录

    /

    返回文章
    返回