气候和人类活动对沿海水资源影响的系统性水文模拟-GSFLOW-郑一-2018

Title:    An integrated hydrological modeling approach for detection and attribution of climatic and human impacts on coastal water resources

Authors:   Feng, Dapeng， Zheng, Yi*，et.al

Journal:   Journal of Hydrology - 2018

1. Introduction

            伴随着人口的急剧增长，良好的水资源供给变得尤为重要。已有 研究表明，区域水文过程广泛的受到人类活动和气候变化的双重影响。但人类活动和气候变化这两个因素各自对于水资源的影响研究较少。区分这两个因素各自影响的方法主要有（1）对比流域法，选取下垫面较为一致，相距较近的两个流域，通过控制变量等措施，研究某个因素改变对于流域水文水资源的影响；（2）趋势分析法，通过使用统计学方法、水文敏感性分析、水量平衡估算等方法，计算出气候变化因素的影响，剩下的影响则归因于人类活动；（3）基于物理过程的水文模拟方法。方法1经济成本高，且很难寻找到完全一致的小流域，在大中型流域研究中难以执行；方法2假定气候和人为因素是各自独立的，这一假定缺乏物理基础 ；方法3理论基础可靠，在许多区域已得到应用。但是更多的是针对地表水、地下水单一流域单元的研究，地表水地下水的耦合模拟及区分气候及人类因素各自影响的研究较少。

1. Materials and methods

      1） 研究区位置：

              滦河流域位于我国东北部，由上游山区和下游平原区两部分组成，最终东流注入渤海，总面积7770平方千米。研究区60%（2010年）土地利用类型为农田，大约55%需进行取水（地表水、地下水）灌溉。长期的农业灌溉及持续增长的人类生产生活用水，加之水库及调水工程建设，导致了自1970年以来地下水的持续消减。

      2）GSFLOW介绍：

            GSFLOW 由美国地质调查局于2008年研发，其主要由PRMs(表层水/根系水模拟）和MODFLOW（包气带水/地下水/河湖水循环模拟）耦合而成，其可对陆地主要水循环过程进行模拟。

      模型中，地表区域被划分为二维水文反馈单元（HRUs）， 地下区域依据水文地质信息纵向划分为多层，横向划分为有限差分网格。 GSFLOW引入了重力水蓄水层的概念，使地表HRUs与其相连的地下网格之间能够进行水交换。此外，含水层单元与其相连的河流或湖泊段之间的直接交换是根据达西定律（Darcy's law）计算的。本文中划分7770个地表单元（1Km*1Km ），3层垂直地下水层。

         模型所需基本参数：

              地学参数： DEM(SRTM90)/土地利用类型图、土壤类型图（黑河数据中心）；

             气象参数（来自7个气象站点）：  最高温/最低温（day）、日降雨；通过反距离加权（ IDW ）内插到每个HRU （1：研究区地势起伏不大，IDW精度满足要求；2 ： 本研究关注的是整个区域的长期水资源预算，而不是区域内水文变量的时空分布（ 说明模型无法模拟区域内水文变量，结果是一个总和 ） ） 。

            在研究区边界4个水文观测站获取每日流入地表水量，收集边界山区侧流作为地下水流入量，大约每年0.1-0.5billion立方米（后期进行校准）。月地下水位数据来自研究区24个监测井1977-1983数据，用于后续模型校准。

        模型所使用的参数没有包括对于农业灌溉用水的模拟。但是由于农业生产中对于地表水的转移和地下水的抽取显著影响着区域水文循环，因此本文通过参考已有文献，确定地表水地下水在农业生产中的总使用量，并按照土地利用类型、生长季与否、农作物产量按比例进行分配。（尽管由上文看出数据大量缺乏，但是 由于本研究侧重于流域尺度水文响应气候变化与人类活动的评价与分离，时空不确定性可能不会对本研究的主要结论产生显著影响，后续将进行研究）

       3）模型校准与验证：

           地下水部分（MODFLOW）：作为模型常规步骤， 首先应进行稳态模拟（ 预热 ），为后续模拟设定初始条件；稳态模拟后进行1975-2000年的水文模拟，其中1975-1976年作为 “spin-up”（ what？ ） ，为后续提供合理的土壤水分状态； 随后的年份又进一步划分为四个应力期:1977-1980年、1981-1990年、1991-1993年和1994-2000年（ 有啥用？ ）；通过适当更改地下水取水量（±10%）分配以及某些地下水参数来校准模型，但应力期总取水量不变； 为了检验校准的有效性，使用参数校准后 数据对第四应力期(1994-2000年)进行模型验证。

           地表水 部分 （PRMs）：由于仅有流域边界4个观测站点，缺乏校准站点，因此只能参考当地水文信息，手动调整具体参数。（ 但是 这种数据匮乏的情况实际上反映了地表水-地下水一体化建模的优势，而不是本研究的不足之处。 在最近的一项研究中( Yao et al.， 2017 )，利用MODFLOW及其径流路径包(SFR2)对青藏高原没有地下水位观测，只能参考径流观测的断块水量补给进行了量化。 ）

       4）不确定性分析：

            设计4个单一参数 极端条件 实验，与原有参数进行模拟结果对比。

      5）归因（分离）分析：

           选取80年代和90年代，两个十年进行对比分析；使用M-K分析（趋势分析、突变检验）查明研究时段内关键营力（  趋势检验前进行自相关分析。如果发现lg-1自相关，则使用Yue等人(2002)的方法对时间序列进行预白？？ 在高阶自相关情况下，M-K检验采用Hamed和Rao(1998)的方法。还进行了非参数秩和检验，以检验这20年间的任何统计上的显著差异 ） ，最终最高温/最低温，降雨（ 气候因素 ）、上游来水（ 气候+人为 ）、地下水抽取量（ 人为因素 ）被选择作为关键营力。基于此，设计5个状态情景来查明水文变化影响因素：（结果相减即得到独立因素的影响！）

      6）未来情景模拟：

           1.温度变化：+0/1/2/3/4℃、 降水变化：+0/5/10/15%、表水转移、地下水抽取、农业灌溉保持历史水平，使用上述情景配合进行模拟（2075-2100年）；

           2.使用21中地球系统模式（ESMs）输出结果进行更系统的模拟， GSFLOW模型在42种气候情景下运行（ 两种排放情景下RCP4.5/8.5的21个ESMs)，地表水引水、地下水抽水和农业灌溉保持在历史水平。 利用模型结果对气候不确定性下的水资源进行了评价。

1. Results and discussion

       1）水量变化（1977-2000年）：降水、ET>表层水输入、输出>地下水输入、输出；如果停止农业灌溉、地下水抽取等人类活动，ET将减少17%更多的地下水地表水流入渤海，但地下水储量将更加稳定。上游水库大坝建设导致的输入水量减少将显著影响输出量，地下水抽取显著降低地下水位。

      2）气候变化及人类活动对于区域地表水地下水均有影响（共同或相反），但侧重点不同。归因分析结果表明：近十年来ET的变化、 出水量变化、 地下水位变化分别 主要归因于气候因素（降水改变所导致的蒸散发量下降）、 上游来水、人类活动。此外3者协作作用导致了地下水位的升高，其原因：气候变化导致的地表干旱使地下水补给降低，地下水头下降，地下水抽取导致地表水头升高，地下水头进一步下降，这将加强 地下水的补给和减少地下水的排放，从而减轻气候变化造成的地下水储量下降（但不显著）。

          不确定分析表明归因分析结果不会因为模拟的不确定性产生很大的差异。

      3）气候及人类影响下， 地下水比地表水 更易变化。气候变化敏感性分析表明：高温将导致更少的区域水量输出和更快的地下水位下降，而温度升高伴随的 降水增加 （21个地球系统模式中19个为此结果，且主要发生在湿季）将缓解上述趋势。但温度持续增加，降水缓解作用将下降；变化气候条件下，地下水量变化与河流输出量相比更敏感。 多个地球系统模式结果也表明地下水储量在假定未来气候下将加速下降，且下降速率RCP8.5>RCP4.5。降低地下水抽取量，不仅能够直接缓解地下水位下降的趋势，同时缓解了气候变化通过影响区域水循环对地下水储量的不利影响；反之，如果地下水加速抽取，气候对于区域水循环的不利影响也将加强。

1. Overview

       模型缺点：无法看出研究区内部单元水量变化。

       如果解决上一问题，可以尝试采用最新的曲面建模HASM方式内插，获得更为精确的气候数据插值结果，通过检查HASM和IDW插值结果来表明曲面建模优势。