1.本发明属于生态水文学中的模型应用及分析评价领域，具体涉及一种基于蒸散发分量的植物水分利用效率评价方法。
背景技术：
：：2.植物水分利用效率(wateruseefficiency,wue)既植物消耗单位水分所能同化的物质量，揭示了生态系统固碳能力与水分消耗之间的动态平衡关系以及植被对气候变化的响应与适应策略，对研究陆地生态系统的碳水耦合机制具有重要意义；且是评估植物耐旱性的重要指标，既在相同条件下，wue高的植物抗旱能力更强。植物wue可用植物净光合速率(netprimaryproductivity,npp)与蒸腾速率的比值来计算，因此变化环境下的蒸腾速率与npp的改变会显著影响wue。由于不同植物和不同的环境条件会造就不同的wue，适应性更强的物种在面对恶劣环境(比如水分胁迫、气温升高等)时常具有更大的竞争力，增加其在整个群落中的优势性，进而改变植被分布格局和群落结构，对于群落演替具有重要作用。通过对wue进行系统的评价，有助于深入了解不同植物在不同环境条件下的适应能力，为提高用水效率提供理论指导，以期从根源上改善水资源开发利用中存在的不合理、低效率等问题。3.蒸散发(evapotranspiration,et)对全球水循环和能量循环至关重要。et可简单地分为湿润表面蒸发(evaporation,e)和植株蒸腾(transpiration,t)。由于蒸腾作用占陆地表面蒸散量的主要部分，并且与植被类型和动态直接相关，因此陆地植被的变化将对t/et比值产生强烈影响。例如，植被覆盖度增加往往会造成t/et增加。然而t/et的野外测量存在难度，尤其是蒸腾量很难实现大面积、长时间连续的监测，更不用说陆地生态系统t/et的长期连续监测。4.目前植被wue的测定方法主要有直接测定法、光合仪测定法(气体交换法)、稳定碳同位素技术以及替代指标法。直接测定法通过测定季节用水和生物量计算wue，这需要大量的观测试验工作，费时费力。而光合仪测定法计算的是叶片光合速率与蒸腾速率的比值，可应用于测定室内植物代谢水利用份额和实际需水量。虽然操作相对简单，研究成本低，但结果往往只能代表某特定短期时间内植物部分叶片的行为。国际公认的判定植物长期wue的最佳方法是稳定碳同位素法，其原理是利用wue与植物碳稳定同位素δ13c值之间显著的线性正相关关系。但该方法成本费用较高，且其作用机理还有待进一步了解。代替指标法具有可观的应用前景，可用于植物长期wue的测定，但该法仍需验证其所用指标与δ13c之间的相关关系，且结果的时空不稳定性较大。上述这些技术提供的信息有限，一定程度上限制了对植物用水效率的综合评价，使这方面研究较为缺乏。5.朱燕燕(2016)基于水资源“有益消耗”和“无益消耗”的概念，将农业水资源有益消耗量定义为农田蒸散发量与有效降雨量之差，利用遥感反演的et产品并结合实测资料，评价了作物灌溉水利用效率。蒋磊等(2013)提出了利用遥感蒸散发模型的灌区蒸散发模拟值，并结合灌区净引水量和降水量观测数据，对灌溉水利用效率进行评价。黄耀欢(2012)将蒸散发作为主要水资源消耗量，并将其界定为参与到经济和生态量产出两大消耗过程，提出基于遥感反演et产品的相对水资源利用效率的概念及其计算方式，对区域水资源利用效率进行综合评价。吕翠美等人(2016)则发明了一种基于用水系统能值分析的用水效率评价方法。受实际生产需求驱动，多数研究聚焦于灌溉水利用效率或农田、作物水分利用效率，缺乏对区域生态系统整体的考虑，且评价指标较为繁杂，未达到统一。在用到遥感蒸散发产品或者模型蒸散发模拟值时，未考虑到e与t两者对整体et的相对贡献率。技术实现要素：6.(一)要解决的技术问题7.本发明要解决的技术问题是如何提供一种基于蒸散发分量的植物水分利用效率评价方法，以解决现有技术中缺乏对区域生态系统整体的考虑，且评价指标较为繁杂，未达到统一，未考虑到e与t两者对整体et的相对贡献率等方面的问题。8.(二)技术方案9.为了解决上述技术问题，本发明提出一种基于蒸散发分量的植物水分利用效率评价方法，该方法包括如下步骤：10.s1、安装rhessys模型；11.s2、预处理模型所需空间分布数据集；12.s3、准备模型输入文件，模型输入文件包括“世界文件(worldfile)、流表(flowtable)、事件发生时间控制文件(tec文件)、默认参数文件和气象数据文件；13.s4、默认参数文件中部分参数的校准与验证；14.s4、基于rhessys的蒸散发模拟，获得不同植被类型下蒸腾量t和蒸发量e的模拟结果；15.s5、计算16.s6、根据s5中计算的蒸腾量t与蒸散发总量e t的比值，评价植物水分利用效率wueact。17.进一步地，所述rhessys为开源模型，支持linux和mac系统。18.进一步地，所述s2中，空间分布数据集包括地理数字信息高程(dem)图，以及根据dem图并使用grassgis软件生成的子流域图、坡度图、坡向图、水平角图、河流网络图以及单元网格图。19.进一步地，所述步骤s3中，世界文件为根据高程、植被类型、土地利用类型和相应参数生成的用于模型模拟的研究区下垫面综合信息文件，流表为河道信息，是根据dem生成的河流地理信息。20.进一步地，世界文件和流表是根据步骤s2中的空间分布数据集，使用r语言的“rhessyspreprocessing”包制作而成。21.进一步地，所述步骤s3中，所述事件发生时间控制文件tec文件用于指定模型模拟结果输出的起止日期和时空尺度。22.进一步地，所述步骤s3中，所述默认参数文件为模型参数的默认值文件paramdb-master。23.进一步地，所述步骤s3中，所述气象数据文件包括日降水量、日最高气温和日最低气温。24.进一步地，所述步骤s4中，校准与验证的参数包括：地表土壤饱和导水率(k)及其随深度的递减率(m)、土壤饱和带到地下含水层的水分传导率(gw1)和决定地下水出流量的系数(gw2)。25.进一步地，所述步骤s6中，植物水分利用效率wueact的比值取值范围为[0,1]，值越接近于1，表明植被用水效率越高[0026](三)有益效果[0027]本发明提出一种基于蒸散发分量的植物水分利用效率评价方法，rhessys可计算流域尺度上的蒸发量e、蒸腾量t这两个计算植物wue的关键因子。模型所需输入数据主要是气象数据(即降水量、最高气温、最低气温)，这在多数地区都有长期连续观测数据集，故受数据源限制较小。且模型能在区域大尺度上进行模拟，能够克服传统蒸散发和wue监测方法的时空限制性，且省时省力、成本低。将模型输出结果与模型输入的土地利用类型以及植被类型相对应，还可用于评估不同植物类型或土地利用方式下的植物水分利用效率，例如草地、常绿阔叶林、落叶阔叶林等。附图说明[0028]图1为本发明方案的主流程图；[0029]图2为几种植被类型的值示意图。具体实施方式[0030]为使本发明的目的、内容和优点更加清楚，下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。[0031]本发明拟利用一种基于物理机理的分布式生态水文模型——rhessys(theregionalhydro-ecologicalsimulationsystem)。该模型可用于综合模拟流域内水、碳和氮循环。模型由气象数据驱动，使用dem图、土地利用类型图和植被类型图来刻画下垫面情况，可分别模拟计算蒸发量(e)和蒸腾量(t)。因此，本发明拟通过模型模拟结果计算t对总蒸散发量(et＝e t)的相对贡献率(即t/et)，提出一种简洁明了的植物水分利用效率的评价方法。[0032]类似于对蒸散发的研究，我们将植物水分利用效率分为实际值(wueact)和潜在值(wuepot)。当土壤蒸发可以忽略时(例如植被覆盖度极高的生态系统)，总蒸散发量et近似于蒸腾量t值，意味着植被极大地利用了环境中的水分，此时的wueact达到最大值，即wuepot。据此，我们可以根据t与et之间的比例关系，得到植被的相对用水效率，建立如下等式关系：[0033][0034]对于给定的环境条件与植物物种类型，一定时期内潜在值wuepot类似于蒸发能力etpot，可假定为常数值。故wueact与t/et线性正相关。当土壤蒸发可忽略时，et近似于t值，既t/et≈1，植被用水效率则高。有研究指出t/et普遍表现出单峰日变化规律，当植被覆盖度高且土壤蒸发量e可以忽略时，日变化峰值将趋于统一(villeggsetal.,2014；wangandyamanaka,2014；zhuetal.,2015)。基于wueact与t/et的这种线性正相关关系，即可利用蒸散发各分量评价植物水分利用效率。[0035]图1是本发明方案的主流程图。如图1所示，本发明提出的wue评价方法具体包括下列步骤：[0036]s1、安装rhessys模型。rhessys为一款开源模型，支持linux和mac系统。模型下载及安装步骤可参见https://github.com/rhessys/rhessys/wiki。[0037]s2、预处理模型所需空间分布数据集。空间分布数据集包括地理数字信息高程(dem)图，以及根据dem图并使用grassgis软件生成的子流域图、坡度图、坡向图、水平角图、河流网络图以及单元网格图。[0038]s3、准备模型输入文件，模型输入文件包括“世界文件(worldfile)、流表(flowtable)、事件发生时间控制文件tec文件、默认参数文件和气象数据文件。[0039]其中，世界文件为根据高程、植被类型、土地利用类型和相应参数生成的用于模型模拟的研究区下垫面综合信息文件，流表为河道信息，是根据dem生成的河流地理信息；世界文件和流表是根据步骤s2中的空间分布数据集，使用r语言的“rhessyspreprocessing”包制作而成。[0040]事件发生时间控制文件tec(temporaleventcontrolfile)文件用于指定模型模拟结果输出的起止日期和时空尺度。[0041]模型的参数文件包括气候条件参数集、土壤参数集、土地利用参数集、植被参数集等，总计约240个参数，模型对大部分参数并不敏感。官方提供了模型参数的默认值文件paramdb-master，故推荐使用模型的默认参数文件。[0042]气象数据文件主要包括日降水量、日最高气温和日最低气温。[0043]s4、默认参数文件中部分参数的校准与验证。步骤s3中所提的默认参数文件中，有一少部分参数对模型模拟结果的准确性影响较大，需要结合研究区实际情况进行校准。通常，校准与验证的参数包括：地表土壤饱和导水率(k)及其随深度的递减率(m)、土壤饱和带到地下含水层的水分传导率(gw1)和决定地下水出流量的系数(gw2)。[0044]s4、基于rhessys的蒸散发模拟，获得不同植被类型下蒸腾量t和蒸发量e的模拟结果。此处以模型示例数据为例，将植被类型分别设置为落叶阔叶林、常绿阔叶林、草地。可得到不同植被类型下蒸腾量t和蒸发量e的模拟结果。[0045]s5、计算结果如图2所示。对于示例数据，几种植被类型的值均大于0.6，植被覆盖状况较好，蒸发消耗的水分不超过可利用水分的40％，植被用水效率较高。[0046]s6、根据s5中计算的蒸腾量t与蒸散发总量e t的比值，评价植物水分利用效率wueact。比值取值范围为[0,1]，值越接近于1，表明植被用水效率越高。评估标准具体可分为以下5类：[0047][0048]以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
：的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。当前第1页12当前第1页12
背景技术：
：：2.植物水分利用效率(wateruseefficiency,wue)既植物消耗单位水分所能同化的物质量，揭示了生态系统固碳能力与水分消耗之间的动态平衡关系以及植被对气候变化的响应与适应策略，对研究陆地生态系统的碳水耦合机制具有重要意义；且是评估植物耐旱性的重要指标，既在相同条件下，wue高的植物抗旱能力更强。植物wue可用植物净光合速率(netprimaryproductivity,npp)与蒸腾速率的比值来计算，因此变化环境下的蒸腾速率与npp的改变会显著影响wue。由于不同植物和不同的环境条件会造就不同的wue，适应性更强的物种在面对恶劣环境(比如水分胁迫、气温升高等)时常具有更大的竞争力，增加其在整个群落中的优势性，进而改变植被分布格局和群落结构，对于群落演替具有重要作用。通过对wue进行系统的评价，有助于深入了解不同植物在不同环境条件下的适应能力，为提高用水效率提供理论指导，以期从根源上改善水资源开发利用中存在的不合理、低效率等问题。3.蒸散发(evapotranspiration,et)对全球水循环和能量循环至关重要。et可简单地分为湿润表面蒸发(evaporation,e)和植株蒸腾(transpiration,t)。由于蒸腾作用占陆地表面蒸散量的主要部分，并且与植被类型和动态直接相关，因此陆地植被的变化将对t/et比值产生强烈影响。例如，植被覆盖度增加往往会造成t/et增加。然而t/et的野外测量存在难度，尤其是蒸腾量很难实现大面积、长时间连续的监测，更不用说陆地生态系统t/et的长期连续监测。4.目前植被wue的测定方法主要有直接测定法、光合仪测定法(气体交换法)、稳定碳同位素技术以及替代指标法。直接测定法通过测定季节用水和生物量计算wue，这需要大量的观测试验工作，费时费力。而光合仪测定法计算的是叶片光合速率与蒸腾速率的比值，可应用于测定室内植物代谢水利用份额和实际需水量。虽然操作相对简单，研究成本低，但结果往往只能代表某特定短期时间内植物部分叶片的行为。国际公认的判定植物长期wue的最佳方法是稳定碳同位素法，其原理是利用wue与植物碳稳定同位素δ13c值之间显著的线性正相关关系。但该方法成本费用较高，且其作用机理还有待进一步了解。代替指标法具有可观的应用前景，可用于植物长期wue的测定，但该法仍需验证其所用指标与δ13c之间的相关关系，且结果的时空不稳定性较大。上述这些技术提供的信息有限，一定程度上限制了对植物用水效率的综合评价，使这方面研究较为缺乏。5.朱燕燕(2016)基于水资源“有益消耗”和“无益消耗”的概念，将农业水资源有益消耗量定义为农田蒸散发量与有效降雨量之差，利用遥感反演的et产品并结合实测资料，评价了作物灌溉水利用效率。蒋磊等(2013)提出了利用遥感蒸散发模型的灌区蒸散发模拟值，并结合灌区净引水量和降水量观测数据，对灌溉水利用效率进行评价。黄耀欢(2012)将蒸散发作为主要水资源消耗量，并将其界定为参与到经济和生态量产出两大消耗过程，提出基于遥感反演et产品的相对水资源利用效率的概念及其计算方式，对区域水资源利用效率进行综合评价。吕翠美等人(2016)则发明了一种基于用水系统能值分析的用水效率评价方法。受实际生产需求驱动，多数研究聚焦于灌溉水利用效率或农田、作物水分利用效率，缺乏对区域生态系统整体的考虑，且评价指标较为繁杂，未达到统一。在用到遥感蒸散发产品或者模型蒸散发模拟值时，未考虑到e与t两者对整体et的相对贡献率。技术实现要素：6.(一)要解决的技术问题7.本发明要解决的技术问题是如何提供一种基于蒸散发分量的植物水分利用效率评价方法，以解决现有技术中缺乏对区域生态系统整体的考虑，且评价指标较为繁杂，未达到统一，未考虑到e与t两者对整体et的相对贡献率等方面的问题。8.(二)技术方案9.为了解决上述技术问题，本发明提出一种基于蒸散发分量的植物水分利用效率评价方法，该方法包括如下步骤：10.s1、安装rhessys模型；11.s2、预处理模型所需空间分布数据集；12.s3、准备模型输入文件，模型输入文件包括“世界文件(worldfile)、流表(flowtable)、事件发生时间控制文件(tec文件)、默认参数文件和气象数据文件；13.s4、默认参数文件中部分参数的校准与验证；14.s4、基于rhessys的蒸散发模拟，获得不同植被类型下蒸腾量t和蒸发量e的模拟结果；15.s5、计算16.s6、根据s5中计算的蒸腾量t与蒸散发总量e t的比值，评价植物水分利用效率wueact。17.进一步地，所述rhessys为开源模型，支持linux和mac系统。18.进一步地，所述s2中，空间分布数据集包括地理数字信息高程(dem)图，以及根据dem图并使用grassgis软件生成的子流域图、坡度图、坡向图、水平角图、河流网络图以及单元网格图。19.进一步地，所述步骤s3中，世界文件为根据高程、植被类型、土地利用类型和相应参数生成的用于模型模拟的研究区下垫面综合信息文件，流表为河道信息，是根据dem生成的河流地理信息。20.进一步地，世界文件和流表是根据步骤s2中的空间分布数据集，使用r语言的“rhessyspreprocessing”包制作而成。21.进一步地，所述步骤s3中，所述事件发生时间控制文件tec文件用于指定模型模拟结果输出的起止日期和时空尺度。22.进一步地，所述步骤s3中，所述默认参数文件为模型参数的默认值文件paramdb-master。23.进一步地，所述步骤s3中，所述气象数据文件包括日降水量、日最高气温和日最低气温。24.进一步地，所述步骤s4中，校准与验证的参数包括：地表土壤饱和导水率(k)及其随深度的递减率(m)、土壤饱和带到地下含水层的水分传导率(gw1)和决定地下水出流量的系数(gw2)。25.进一步地，所述步骤s6中，植物水分利用效率wueact的比值取值范围为[0,1]，值越接近于1，表明植被用水效率越高[0026](三)有益效果[0027]本发明提出一种基于蒸散发分量的植物水分利用效率评价方法，rhessys可计算流域尺度上的蒸发量e、蒸腾量t这两个计算植物wue的关键因子。模型所需输入数据主要是气象数据(即降水量、最高气温、最低气温)，这在多数地区都有长期连续观测数据集，故受数据源限制较小。且模型能在区域大尺度上进行模拟，能够克服传统蒸散发和wue监测方法的时空限制性，且省时省力、成本低。将模型输出结果与模型输入的土地利用类型以及植被类型相对应，还可用于评估不同植物类型或土地利用方式下的植物水分利用效率，例如草地、常绿阔叶林、落叶阔叶林等。附图说明[0028]图1为本发明方案的主流程图；[0029]图2为几种植被类型的值示意图。具体实施方式[0030]为使本发明的目的、内容和优点更加清楚，下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。[0031]本发明拟利用一种基于物理机理的分布式生态水文模型——rhessys(theregionalhydro-ecologicalsimulationsystem)。该模型可用于综合模拟流域内水、碳和氮循环。模型由气象数据驱动，使用dem图、土地利用类型图和植被类型图来刻画下垫面情况，可分别模拟计算蒸发量(e)和蒸腾量(t)。因此，本发明拟通过模型模拟结果计算t对总蒸散发量(et＝e t)的相对贡献率(即t/et)，提出一种简洁明了的植物水分利用效率的评价方法。[0032]类似于对蒸散发的研究，我们将植物水分利用效率分为实际值(wueact)和潜在值(wuepot)。当土壤蒸发可以忽略时(例如植被覆盖度极高的生态系统)，总蒸散发量et近似于蒸腾量t值，意味着植被极大地利用了环境中的水分，此时的wueact达到最大值，即wuepot。据此，我们可以根据t与et之间的比例关系，得到植被的相对用水效率，建立如下等式关系：[0033][0034]对于给定的环境条件与植物物种类型，一定时期内潜在值wuepot类似于蒸发能力etpot，可假定为常数值。故wueact与t/et线性正相关。当土壤蒸发可忽略时，et近似于t值，既t/et≈1，植被用水效率则高。有研究指出t/et普遍表现出单峰日变化规律，当植被覆盖度高且土壤蒸发量e可以忽略时，日变化峰值将趋于统一(villeggsetal.,2014；wangandyamanaka,2014；zhuetal.,2015)。基于wueact与t/et的这种线性正相关关系，即可利用蒸散发各分量评价植物水分利用效率。[0035]图1是本发明方案的主流程图。如图1所示，本发明提出的wue评价方法具体包括下列步骤：[0036]s1、安装rhessys模型。rhessys为一款开源模型，支持linux和mac系统。模型下载及安装步骤可参见https://github.com/rhessys/rhessys/wiki。[0037]s2、预处理模型所需空间分布数据集。空间分布数据集包括地理数字信息高程(dem)图，以及根据dem图并使用grassgis软件生成的子流域图、坡度图、坡向图、水平角图、河流网络图以及单元网格图。[0038]s3、准备模型输入文件，模型输入文件包括“世界文件(worldfile)、流表(flowtable)、事件发生时间控制文件tec文件、默认参数文件和气象数据文件。[0039]其中，世界文件为根据高程、植被类型、土地利用类型和相应参数生成的用于模型模拟的研究区下垫面综合信息文件，流表为河道信息，是根据dem生成的河流地理信息；世界文件和流表是根据步骤s2中的空间分布数据集，使用r语言的“rhessyspreprocessing”包制作而成。[0040]事件发生时间控制文件tec(temporaleventcontrolfile)文件用于指定模型模拟结果输出的起止日期和时空尺度。[0041]模型的参数文件包括气候条件参数集、土壤参数集、土地利用参数集、植被参数集等，总计约240个参数，模型对大部分参数并不敏感。官方提供了模型参数的默认值文件paramdb-master，故推荐使用模型的默认参数文件。[0042]气象数据文件主要包括日降水量、日最高气温和日最低气温。[0043]s4、默认参数文件中部分参数的校准与验证。步骤s3中所提的默认参数文件中，有一少部分参数对模型模拟结果的准确性影响较大，需要结合研究区实际情况进行校准。通常，校准与验证的参数包括：地表土壤饱和导水率(k)及其随深度的递减率(m)、土壤饱和带到地下含水层的水分传导率(gw1)和决定地下水出流量的系数(gw2)。[0044]s4、基于rhessys的蒸散发模拟，获得不同植被类型下蒸腾量t和蒸发量e的模拟结果。此处以模型示例数据为例，将植被类型分别设置为落叶阔叶林、常绿阔叶林、草地。可得到不同植被类型下蒸腾量t和蒸发量e的模拟结果。[0045]s5、计算结果如图2所示。对于示例数据，几种植被类型的值均大于0.6，植被覆盖状况较好，蒸发消耗的水分不超过可利用水分的40％，植被用水效率较高。[0046]s6、根据s5中计算的蒸腾量t与蒸散发总量e t的比值，评价植物水分利用效率wueact。比值取值范围为[0,1]，值越接近于1，表明植被用水效率越高。评估标准具体可分为以下5类：[0047][0048]以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
：的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。当前第1页12当前第1页12
再多了解一些
本文用于企业家、创业者技术爱好者查询，结果仅供参考。
