CN105701360B - 一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 - Google Patents
一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 Download PDFInfo
- Publication number
- CN105701360B CN105701360B CN201610117048.XA CN201610117048A CN105701360B CN 105701360 B CN105701360 B CN 105701360B CN 201610117048 A CN201610117048 A CN 201610117048A CN 105701360 B CN105701360 B CN 105701360B
- Authority
- CN
- China
- Prior art keywords
- precipitation
- carbon flux
- events
- event
- carbon
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 83
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 83
- 230000004907 flux Effects 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000011156 evaluation Methods 0.000 title abstract 2
- 238000001556 precipitation Methods 0.000 claims abstract description 92
- 230000004044 response Effects 0.000 claims abstract description 33
- 230000035945 sensitivity Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 abstract description 11
- 238000004177 carbon cycle Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 7
- 230000000241 respiratory effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000036387 respiratory rate Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000004158 soil respiration Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于箱式法和考虑降水影响的干旱区碳通量估算方法，该方法包括四个步骤：（1）碳通量对不同大小降水事件响应的敏感性观测；（2）降水事件的等级划分，将降水事件由小至大依次划分成5个等级（P 0 、P 1 、P 2 、P 3 和P 4 ），其中P 0 为零降水（无降水事件发生），P 1 、P 2 和P 3 为降水事件发生概率主体，P 4 为极端大降水事件；（3）构建响应模型及参数化，构建碳通量对降水事件的响应模型和确定相应的模型参数；（4）碳通量估算，根据研究区不同等级降水事件统计结果和响应模型，通过计算各等级降水事件碳通量的微积分并求和，得到估算时段内的碳通量估算值。本发明具有计算简单、估算准确性高和野外工作量低等优点，可广泛应用于干旱区碳通量监测和碳循环等相关研究领域。
Description
技术领域
本发明涉及一种碳通量估算方案，具体来说是一种基于箱式法和考虑降水对碳通量波动影响的碳通量估算方法。
背景技术
碳通量估算是全球变化背景下研究碳循环和碳收支的关键技术环节。利用箱式法（静态箱法）进行碳通量观测，具有操作简单、成本低、可实现多种碳通量组分观测等优势，得到国际上普遍认可和广泛采用。然而，箱式法在估算季节或年际等相对较长时间尺度碳能量时，传统的方法多是将多次观测值采用数学平均值法方法进行估算，但这种估算方法存在两个不可回避的问题，一是需要在野外密集而繁重的观测工作；二是忽视了干旱区降水事件对碳通量具有显著波动影响这一客观事实。大量研究证明，干旱区的碳交换过程对降水事件极其敏感，如在5mm小降水事件发生后，生态系统碳吸收和碳释放速率会增加数倍并持续数天，特别是土壤呼吸在几小时内即可增加几十倍。由此可知，传统的数学平均值法对碳收支的估算具有观测工作量大，特别是估算结果且具有很大的不确定性的缺陷。若观测事件多在降水前则造成偏低估算；若观测事件多在降水事件刚刚发生后则造成偏高估算。因此，如何考虑降水事件对碳通量的波动影响，而提高箱式法对碳通量估算的合理性和准确性，是干旱区碳循环和碳收支相关研究亟待解决的重要技术问题。
发明内容
本发明目的在于提供一种充分考虑降水波动影响的碳通量估算的方法，是一种基于箱式法的碳通量观测，能有效提高干旱区碳通量估算的准确性，且降低野外观测工作量的方法。
本发明的目的主要通过以下4个步骤实现（图1）：
1. 碳通量对不同大小降水事件响应的敏感性观测
采用箱式法利用碳通量观测仪器（如Li-840、Li-6400或Li-8150），观测不同大小降水事件发生后碳通量日平均交换速率的动态变化，持续数日，直至日平均交换速率恢复到干旱状态水平。
上述的不同大小的降水事件可以是自然降水事件，亦可以是人工模拟降水事件。
2. 降水事件的等级划分
根据干旱区降水事件大小分布频率特征（小于5mm的降水事件为体，占总体降水事件的60-80%），将降水事件由小至大依次划分为5个等级，即P 0 、P 1 、P 2 、P 3 和P 4 。其中，P 0 为无降水事件发生（零降水）； P 1 、P 2 和P 3 为降水事件发生概率的主体，分别对应的降水事件发生频率分别为50%、30%和20%；P 4 为极端大降水事件，发生概率较小。
3. 构建响应模型及参数化
根据实地观测获得降水事件发生后碳通量监测数据以及国内外干旱区相关的数据资料的统计结果表明，碳通量交换速率对降水事件的理论响应过程是一种非线性过程，即降水事件发生后先迅速增加直至峰值，然后逐渐降低并恢复至干旱状态水平。上述碳通量交换速率对不同等级降水事件响应的动态过程论模型如图2所示。在上述步骤1和步骤2的数据基础上，利用专业数学软件Matlab7.0构建碳通量对降水事件的响应模型，其数学表达式为：
式中f(x)为碳通量的响应量；x为响应持续时间；y 0 、x c 、w和A为模型的参数；π和e为常数，数值分别为3.142和2.718。
上述的响应模型适用于不同等级降水事件，模型的具体参数由降水事件等级和季节有关。
上述响应模型可通过观测事件对模型参数进行验证和校正，以提高模型模拟的精度。
4.碳通量估算
根据上述响应模型和研究区估算时间段内不同等级降水事件分布频率的统计结果，分别计算P 1 、P 2 、P 3 和P 4 相应等级降水事件碳通量的微积分，分别得到不同等级降水事件碳通量估算值；将干旱状态碳通量平均值与干旱时间步长（天数）乘积，得到无降水发生（P 0 ）时段的碳通量估算值；最后，将各等级降水事件碳通量估算值求和，得到估算时间段内的碳通量估算值。
上述计算过程的数学公式为：
式中，为碳通量估算值；为无降水事件发生时（P 0 ）碳通量平均值；d 0 为无降水时间步长（无降水天数）；n 1 、n 2 、n 3 和n 4 分别为估算时间段内P 1 、P 2 、P 3 和P 4 降水事件发生的次数。
上述计算公式中的d 0 、n 1 、n 2 、n 3 和n 4 是由估算时间段研究区降水资料统计获得。
上述计算的具体过程可通过编程实现或专业数学软件Matlab7.0实现。
本发明的优点及有益效果：
1）本发明所建立的响应模型为一种数学统计模型，模型的构建过程和参数简单，具体应用时不需要输入气象、植被、土壤等复杂的参数。
2）本发明的估算方法，充分考虑了干旱区碳通量对降水事件的响应极其敏感这一客观现实。避免传统数学平均值法忽视了降水事件对碳通量的波动影响而导致估算结果不可靠性，有效提高了相对长时间尺度对碳通量估算的准确性和合理性。
3）传统方法在估算长时间尺度碳通量时，需要在野外进行密集而繁重的观测工作，本发明的估算方法仅需要对典型等级的降水事件进行碳通量观测，显著降低了野外观测工作量。
4）本发明的估算方法可根据研究区降水资料和响应模型，可实现对观测空白的碳通量可靠估算。一定程度上避免了箱式法难以实现连续观测的缺陷，拓展了箱式法可应用性。
5）本发明的估算方法在干旱区的适用性具有普遍性，在应用时仅需要根据实地观测数据对模型参数进行校正即可。
附图说明
图1 是本发明碳通量估算方法的步骤流程图。
图2 是碳通量对不同等级降水事件响应的理论模型示意图。
图中：P 0 、P 1 、P 2 、P 3 和P 4 为按降水事件由大至小划分的5个降水事件等级，其中P 0 表示无降水事件发生。
具体实施方式
下面结合对典型干旱区腾格里沙漠东南缘（37°29′ N, 104°25′E）碳通量的估算，对本发明进一步详细说明。
本发明的估算方法通过以下步骤实现：
1. 碳通量对降水事件响应的敏感性观测
采用箱式法利用CO2/H2O气体分析仪（LI-840A, Li-Cor, USA），于2012-2014年7-9月份植物生长旺季，在野外观测了不同大小自然降水事件发生后净碳通量（NEE）和生态系统总呼吸（RE）日交换速率的动态变化，每次观测持续数日，直至碳通量日交换速率恢复到干旱状态常值。
2. 降水事件的等级划分
根据腾格里沙漠东南缘多年降水资料的统计分析，将降水事件由小至大依次划分为5个等级，即P 0 、P 1 、P 2 、P 3 和P 4 。其中，P 0 为无降水事件发生；P 1 、P 2 和P 3 为降水事件发生概率的主体，分别对应的降水事件发生频率为50%、30%和20%；P 4 为极端大降水事件，发生概率较小。按上述降水事件的等级划分方法，本研究区P 1 、P 2 、P 3 和P 4 代表的具体降水事件范围分别为0mm<P<3mm、3mm≤P<5mm、5mm≤P<10mm和P≥10mm。
3. 构建响应模型及参数化
根据对本研究区持续3年的碳通量观测资料，将碳通量对降水事件的理论响应过程（图2）进行量化。在上述步骤1和步骤2的数据基础上，利用专业数学软件Matlab7.0构建碳通量对降水事件的响应模型，其数学表达式为：
式中f(x)为碳通量的响应量；x为响应持续时间；y 0 、x c 、w、A和为模型的参数；π和e为常数，数值分别为3.142和2.718。
根据后续观测数据对模型进行了验证和参数校正，具体参数如下表：
腾格里沙漠东南缘碳通量对降水事件的响应模型参数表：
4.碳通量估算
根据上述响应模型和本研究区植物生长季不同等级降水事件分布频率的统计数据，分别计算P 1 、P 2 、P 3 和P 4 各等级降水事件碳通量的微积分，得到不同等级降水事件碳通量估算值；干旱状态均值与干旱时间步长（天数）之乘积，得到无降水发生（P 0 ）时段的碳通量估算值；最后，将各等级降水事件碳通量估算值求和，得到估算时间段内的碳通量估算值。
上述计算过程的数学公式为：
式中，为估算时间段内的碳通量估算值；为无降水事件发生时碳通量平均值，即干旱状态平均值；d 0 为无降水时间步长（无降水天数）n 1 、n 2 、n 3 和n 4 分别为估算时间段内P 1 、P 2 、P 3 和P 4 降水事件发生的次数。
上述计算公式中的d 0 、n 1 、n 2 、n 3 和n 4 由估算时间段内的本研究区降水资料统计获得。
上述计算的具体过程通过专业数学软件Matlab7.0实现。
根据上述计算方法，将相应模型参数和降水事件统计结果代入响应模型，最终得到本研究区2012年、2013年和2014年7-9月份的净碳交换量速率（NEE）分别为-49.46 g m-2、-37.34 g m-2和-41.99 g m-2，相应的生态系统总呼吸速率（RE）分别为106.25 g m-2、81.66 g m-2和110.22 g m-2。
以上内容是对本发明优选的实施例的说明，可以帮助本领域技术人员更充分地理解本发明的技术方案。但这些实施例仅仅是举例说明，不能认定本发明的具体实施方式仅限于这些实施例的说明。
Claims (2)
1.一种基于箱式法和考虑降水影响的干旱区碳通量估算方法，其特征在于包括以下步骤：
（1）碳通量对不同大小降水事件响应的敏感性观测；
（2）降水事件的等级划分；
（3）构建响应模型及参数化；
（4）碳通量估算；
所述碳通量对不同大小降水事件响应的敏感性观测，具体包括：
采用箱式法利用碳通量观测仪器，观测不同大小降水事件发生后碳通量日平均交换速率的动态变化，持续数日，直至日平均交换速率恢复到干旱状态水平；上述的不同大小降水事件可以是自然降水事件，亦可以是人工模拟降水事件；
所述构建的响应模型的数学表达式为：
式中f(x)为碳通量的响应量；x为响应持续时间；y0、xc、w和A为模型的参数；π和e为常数，数值分别为3.142和2.718；
所述的步骤（4）的估算时间段的碳通量估算的数学公式为：
式中，为估算时间段内的碳通量估算值；为无降水事件发生时碳通量平均值，即干旱状态平均值；d0为无降水时间步长；
所述步骤（2）的降水事件的等级划分依据是干旱区降水事件大小分布频率特征，将降水事件由小至大依次划分为5个等级，即P 0 、P 1 、P 2 、P 3 和P 4 ；其中，P 0 为零降水；P 1 、P 2 和P 3 为降水事件发生概率的主体，分别对应的降水事件频率为50%、30%和20%；P 4 为极端大降水事件；P 1 、P 2 、P 3 和P 4 代表的具体降水事件范围分别为0mm<P<3mm、3mm≤P<5mm、5mm≤P<10mm和P≥10mm；n1、n2、n3和n4分别为估算时间段内P1、P2、P3和P4降水事件发生的次数。
2.如权利要求1所述的基于箱式法和考虑降水影响的干旱区碳通量估算方法，其特征在于：所述步骤（3）的响应模型适用于不同等级降水事件，模型的具体参数与降水事件等级和季节有关。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610117048.XA CN105701360B (zh) | 2016-03-02 | 2016-03-02 | 一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610117048.XA CN105701360B (zh) | 2016-03-02 | 2016-03-02 | 一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105701360A CN105701360A (zh) | 2016-06-22 |
| CN105701360B true CN105701360B (zh) | 2018-11-13 |
Family
ID=56222657
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610117048.XA Active CN105701360B (zh) | 2016-03-02 | 2016-03-02 | 一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105701360B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106568920A (zh) * | 2015-10-13 | 2017-04-19 | 中国科学院沈阳应用生态研究所 | 一种土壤呼吸代表性监测时段确定方法 |
| CN108828158B (zh) * | 2018-07-31 | 2024-01-30 | 北京市水科学技术研究院 | 一种人工降雨条件下碳通量的测量设备 |
| CN109856362A (zh) * | 2018-11-30 | 2019-06-07 | 中国水利水电科学研究院 | 一种涝渍事件对农作物土体单元碳通量的影响识别方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102507906A (zh) * | 2011-11-15 | 2012-06-20 | 莫路锋 | 一种基于大范围土壤碳通量监测系统的wsn森林环境效益监测系统 |
| CN102506939A (zh) * | 2011-11-15 | 2012-06-20 | 周国模 | 一种大范围土壤碳通量监测系统及方法 |
| CN102721435A (zh) * | 2011-11-15 | 2012-10-10 | 周国模 | 一种基于wsn的土壤表面co2通量监测系统 |
| CN103630667A (zh) * | 2013-11-19 | 2014-03-12 | 南京信息工程大学 | 一种快速测定土壤二氧化碳通量的装置及方法 |
| CN104897873A (zh) * | 2014-12-26 | 2015-09-09 | 浙江农林大学 | 一种开放型土壤碳通量监测仪及监测方法 |
-
2016
- 2016-03-02 CN CN201610117048.XA patent/CN105701360B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102507906A (zh) * | 2011-11-15 | 2012-06-20 | 莫路锋 | 一种基于大范围土壤碳通量监测系统的wsn森林环境效益监测系统 |
| CN102506939A (zh) * | 2011-11-15 | 2012-06-20 | 周国模 | 一种大范围土壤碳通量监测系统及方法 |
| CN102721435A (zh) * | 2011-11-15 | 2012-10-10 | 周国模 | 一种基于wsn的土壤表面co2通量监测系统 |
| CN103630667A (zh) * | 2013-11-19 | 2014-03-12 | 南京信息工程大学 | 一种快速测定土壤二氧化碳通量的装置及方法 |
| CN104897873A (zh) * | 2014-12-26 | 2015-09-09 | 浙江农林大学 | 一种开放型土壤碳通量监测仪及监测方法 |
Non-Patent Citations (2)
| Title |
|---|
| The effects of extreme rainfall events on carbon release from biological soil crusts covered soil in fixed sand dunes in the Tengger Desert, northern China;Yang Zhao et al.;《Sciences in Cold and Arid Regions》;20131231;第5卷(第2期);第191-196页 * |
| 模拟降水对古尔班通古特沙漠生物结皮表观土壤碳通量的影响;吴林 等;《生态学报》;20120731;第32卷(第13期);第4105-4106页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105701360A (zh) | 2016-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Gaydon et al. | Evaluation of the APSIM model in cropping systems of Asia | |
| Lu et al. | Quantifying the impacts of small dam construction on hydrological alterations in the Jiulong River basin of Southeast China | |
| Li et al. | Impact assessment of climate change on rice productivity in the Indochinese Peninsula using a regional‐scale crop model | |
| Awotwi et al. | Climate change impact on streamflow in a tropical basin of Ghana, West Africa | |
| Pennock et al. | The role of field studies in landscape-scale applications of process models: an example of soil redistribution and soil organic carbon modeling using CENTURY | |
| CN105701360B (zh) | 一种基于箱式法和考虑降水影响的干旱区碳通量估算方法 | |
| Pandey et al. | The role of modelling in the quest for sustainable farming systems | |
| CN112715322B (zh) | 一种农业灌溉用水获取方法和装置 | |
| Cao et al. | Computing statistical indices for hydrothermal times using weed emergence data | |
| Verma et al. | Comparative analysis of CMIP5 and CMIP6 in conjunction with the hydrological processes of reservoir catchment, Chhattisgarh, India | |
| Yonaba et al. | Future climate or land use? Attribution of changes in surface runoff in a typical Sahelian landscape | |
| Kisekka et al. | Crop modeling applications in agricultural water management | |
| Tatsumi | Effects of automatic multi-objective optimization of crop models on corn yield reproducibility in the USA | |
| CN114365682A (zh) | 一种设施栽培土壤水分预测方法、装置及电子设备 | |
| CN116401882A (zh) | 基于多因子胁迫的地表臭氧对冬小麦产量影响评估方法 | |
| Uddameri et al. | A tiered stochastic framework for assessing crop yield loss risks due to water scarcity under different uncertainty levels | |
| Gemechu et al. | The current and future trend of rainfall and its variability in Adami-Tulu Jido-Kombolcha Woreda, Central Rift Valley of Ethiopia | |
| Altaweel | Simulating the effects of salinization on irrigation agriculture in southern Mesopotamia | |
| Gao | Dealing with missing data in hydrology: Data analysis of discharge and groundwater time-series in Northeast Germany | |
| Kumar et al. | Use of a decision support system to establish the best model for estimating reference evapotranspiration in sub-temperate climate: Almora, Uttarakhand. | |
| Hoang et al. | A robust parameter approach for estimating CERES-Rice model parameters for the Vietnam Mekong Delta | |
| Zhang et al. | A universal agro-hydrological model for water and nitrogen cycles in the soil–crop system SMCR_N: Critical update and further validation | |
| Amiri et al. | Comparison of the AquaCrop and CERES‐Maize models for simulating maize phenology, grain yield, evapotranspiration and water productivity under different irrigation and nitrogen levels | |
| Acharjee et al. | Yield and water productivity variation of Boro rice with irrigation strategies and transplanting dates under climate change–a case study in south-western Bangladesh | |
| Golabi et al. | Performance Evaluation of the BUDGET Model in Simulating Different Strategies of Irrigation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wang Zengru Inventor after: Li Xinrong Inventor after: Liu Lichao Inventor after: Gao Yanhong Inventor after: Yang Haotian Inventor before: Wang Zengru Inventor before: Li Xinrong Inventor before: Li Lichao Inventor before: Gao Yanhong Inventor before: Yang Haotian |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200812 Address after: 730000 No.318, Donggang West Road, Chengguan District, Lanzhou City, Gansu Province Patentee after: NORTHWEST INSTITUTE OF ECO-ENVIRONMENT AND RESOURCES, CAS Address before: 730000 Gansu city of Lanzhou province Donggang West Road No. 320 Patentee before: Institute of environment and Engineering in cold and dry areas, Chinese Academy of Sciences |
|
| TR01 | Transfer of patent right |