基于通量观测的攀西烟区烟田水分利用效率特征研究

徐同庆，徐宜民，王程栋，刘光亮，王松峰，董建新，王树声，陶健

1 中国农业科学院烟草研究所，青岛266101；2 中国农业科学院研究生院，北京100081

农艺与调制

基于通量观测的攀西烟区烟田水分利用效率特征研究

徐同庆1,2，徐宜民1，王程栋1，刘光亮1，王松峰1，董建新1，王树声1，陶健1

1 中国农业科学院烟草研究所，青岛266101；2 中国农业科学院研究生院，北京100081

为揭示烟田生态系统水、碳通量和水分利用效率（WUE）的变化特征。基于2015—2016年攀西烟区涡度相关通量观测数据，分析烟田生态系统水、碳通量和WUE的时间动态特征，阐释烟田生态因子对水、碳通量和WUE的影响规律。结果表明：烟田WUE的日变化为日出后迅速升高，8:00～9:00达到最大值后逐渐降低，到14:00左右趋于稳定，傍晚前又有小幅度的增加；烟田WUE的季节变化呈“双峰”趋势，旺长期和成熟采收期前期较高，伸根期最低。烟田WUE与气温、净辐射均呈二次曲线关系，而与饱和水汽压差（VPD）呈负指数关系。气温、VPD和净辐射控制WUE变化的时间节点不同，其中VPD对WUE限制节点出现的时间较早，而气温、净辐射相对较晚，这主要与总初级生产力（GPP）和蒸散量（ET）对气温、VPD和净辐射变化的敏感度差异有关。根据不同生育期WUE变化规律的研究结果，降低土壤无效蒸发，增强叶片光合能力，是提高烟田生态系统水资源利用水平的关键途径。

水分利用效率；水碳通量；通量观测；烟田生态系统

烟草是重要的经济作物，提升优质烟叶的产量及比例是提高烟草种植收益的关键途径[1]。烟株通过光合作用吸收水和CO2并将其转化为有机物，实现自身生长发育[2-3]，水、碳循环过程是烟株生长发育的关键过程，烟株不同生育期水、碳循环过程动态协调是实现优质烟叶生产的重要基础[4]。烟田生态系统水分利用效率（water use efficiency, WUE）是指烟株消耗单位质量水分所同化的CO2的量，是表征烟田生态系统水、碳循环耦合关系的重要指标[5-6]。WUE不仅反映了烟株在不同生育期水分利用能力，还反映了烟株在环境因子动态影响下的响应与适应策略[7-9]。

涡度相关通量观测（eddy covariance fl ux observation）是一种基于大气湍流理论和数据统计分析相结合的技术，通过快速测定大气的物理量（如温度、湿度、CO2浓度等）与其垂直风速的协方差来计算湍流通量[10]。涡度相关技术是当前研究生态系统水、碳通量的主要手段，它不仅可以用来估算碳平衡组分[11-12]，还可以用来研究生态系统水、碳通量对环境变化的响应[13]。利用水、碳通量的观测数据来研究环境因子对生态系统WUE的驱动机理，能够深入理解气候变化背景下生态系统水分利用策略的变化[14-15]。然而，目前基于通量观测数据对生态系统水、碳通量和WUE变化特征及影响机理的研究主要集中于森林和草地生态系统，对农田生态系统尤其是烟田生态系统的研究较少[16-17]。战领等[18]基于通量数据研究了农田生态系统水、碳通量及WUE的变化规律，分析了影响WUE变化的相关因素；Tang等[19]分析了农田生态系统WUE的变化特征，并利用通量数据对MODIS的生产力产品数据进行了精度检验；Zhang等[20]利用通量数据并结合模型模拟指出了作物叶面积指数（leaf area index, LAI）对生态系统WUE的影响。当前，有关烟田生态系统水、碳通量及WUE方面的研究较少[21]，开展烟田生态系统水、碳循环及WUE的研究，有利于根据不同烟草种植区生态条件制定科学的水分管理措施，同时根据烟株不同生育期水分需求和消耗状况合理分配水资源，提高烟田生态系统的生产稳定性和应对气候变化的能力。

本研究基于2015—2016年的涡度相关通量观测数据，研究攀枝花和西昌烟区（以下简称“攀西烟区”）烟田生态系统大田生育期水、碳通量及WUE在不同时间尺度上的动态变化特征，分析了主要的生态因子对WUE的驱动机理，为进一步深入开展烟田生态系统WUE研究提供理论和数据支持。

1 材料与方法

1.1 研究区概况

攀西烟区地处四川西南的河谷农业区，光照资源丰富，昼夜温差大，年均温较高，是四川省主要植烟区之一，也是我国主要的烤烟适宜生态区之一，烟叶品质存在显著的区域化特征，烟叶风格以“清甜香”而独具特色[22]。本试验设置在四川省凉山州西昌市大兴乡的中国农业科学院烟草研究所西南基地（北纬27°49′，东经 102°22′，海拔 1700 m），属于亚热带高原季风气候。试验地烟田生态因子在整个攀西烟区具有典型代表性，年降水量在1000～1500 mm，且集中在夏季，多年平均气温17.2℃；土壤多为红壤和紫色土（酸性），种植制度为烤烟/大蒜或烤烟/水稻轮作一年两熟制。研究区烤烟大田生育期约140 d，4月中下旬移栽，7月中旬打顶，9月下旬烟叶采收结束。结合当地烤烟种植对生育期的划分经验，将烤烟大田生育期大致划分为4个阶段：4月中下旬～5月上旬为移栽还苗期，5月中旬～6月中旬为伸根期，6月中下旬～7月中旬为旺长期，7月中下旬～9月下旬为成熟采收期。

1.2 观测方法

涡度相关通量观测系统包括EC150开路式CO2/H2O气体分析仪、CSAT3型三维超声风速仪、HMP45C空气温湿度探头及CNR4净辐射传感器，通过数据采集器（CR5000）24 h连续自动采集，采样频率为10 Hz。涡度相关通量观测系统可以测量CO2绝对密度以及水汽密度、空气温度、大气压力。此外结合三维超声风速仪和净辐射传感器，还可以测量三维风速、超声空气温度以及太阳净辐射。系统在运行过程中，超声风速仪高频响应三维风速和虚温，CO2/H2O分析仪高频响应CO2和H2O，数据采集器实时采集这些变量数据，并对其做同步处理，之后在线计算得到感热通量、水汽通量、CO2通量、动量通量、摩擦风速，以及这些数据所需的协方差/均值等，并将计算结果保存在数据采集单元，同时各种高频变量的原始数据也会保存在数据采集单元中。最后将保存的数据进行分析和插补，从而获得理想通量数据。

1.3 数据的处理与计算

本研究基于研究区2015—2016年通量观测数据，根据有关方法计算获得大田生育期气象数据，包括降水、气温、饱和水汽压差（VPD）、净辐射，以及烟田生态系统水、碳通量数据，包括CO2通量、生态系统呼吸（Reco）、总初级生产力（GPP）、蒸散量（ET）及WUE[16,23-25]，文中日变化数据均采用生育期内平均值。目前关于30 min CO2通量观测数据的处理已基本形成一个完整的程序，即首先利用原始数据在3次坐标旋转、密度校正（WPL校正）及储存项计算的基础上计算出湍流通量，然后对通量数据进行质量控制、剔除异值，最后再对缺失数据进行插补。

生态系统GPP计算公式为[26-28]：

其中，NEP为净生态系统生产力，Reco为生态系统呼吸。Reco与温度相关，计算公式为[29-30]：

其中，R10为参考呼吸速率，Ta为气温（℃），k为系数。R10与k值通过拟合夜间CO2通量与气温的指数关系获得（图1）。此外，Reco包括夜间和日间两部分，日间Reco的计算通常利用夜间Reco与气温所建立的指数函数外推得到[26]。

生态系统WUE（gC·kg-1H2O）的计算公式为[31]：

图1 夜间CO2通量对气温的响应曲线Fig. 1 Curve of nighttime CO2 fl uxes against air temperature

2 结果与分析

2.1 气象因子状况

图2为不同季节气温、降水、净辐射、VPD的日均值动态，表1为不同生育期气温、降水、净辐射、VPD的统计结果。由图2a和表1可知，研究区大田生育期日均气温21.2℃，在4月中旬～5月下旬整体呈上升趋势，6～7月相对稳定，8月上旬～8月中旬气温再次升高，8月下旬以后迅速下降。其中旺长期和成熟采收期日均气温相对较高，分别为22.3℃和22.2℃，移栽还苗期日均气温最低，仅18.7℃。

由图2b和表1可知，研究区大田生育期降水总量为476 mm，但季节分配不均匀，其中旺长期降水量相对集中，占整个大田生育期的38.4%；成熟采收期和移栽还苗期降水量分别为142 mm和108 mm，分别占整个大田生育期的29.8%和22.7%；伸根期降水量最少，仅占整个大田生育期总降水量的9%。

由图2c和表1可知，研究区烟田净辐射整体呈现先升高后下降的趋势，在整个大田生育期日均净辐射为365 W·m-2，其中旺长期最高，为426 W·m-2，移栽还苗期和成熟采收期相对较低，分别为318 W·m-2和305 W·m-2。与烟田日均净辐射在不同生育期的变化规律相反，VPD在整个大田生育期呈现先降低后升高的变化趋势，整个大田生育期日均VPD为1602 Pa。其中，旺长期最低，为1369 Pa，移栽还苗期最高，为1716 Pa。

图2 气温、降水、净辐射及VPD日均值的季节变化特征Fig. 2 Seasonal trends of daily air temperature, precipitation, net radiation and VPD

表1 不同生育期气温、降水、净辐射和VPD统计结果Tab.1 Statistical results of air temperature, precipitation, net radiation and VPD at different growing stages

2.2 水、碳通量变化特征

2.2.1 水、碳通量日变化特征

由图3a可知，烟田生态系统CO2通量的日变化规律明显（负值表示碳吸收）。日出后，烟株CO2吸收能力逐渐增强，在12:00～14:00达到峰值；此后，CO2吸收量有一个小幅度下降的过程，该现象在旺长期表现最明显；15:00以后，CO2吸收量逐渐下降，并在日落前后（大约19:00）接近于0；CO2通量日动态在不同生育期间存在显著差异，伸根期CO2通量的峰值最低，单位时间内平均碳吸收量为0.38 gC·m-2·h-1，峰值出现的时间为16:00左右；成熟采收期CO2通量的峰值最高，单位时间内平均碳吸收量为1.1 gC·m-2·h-1，峰值出现的时间为12:00左右。

由图3b可知，烟田GPP的日变化特征呈“升高-稳定-下降”的规律，且峰值出现时间比CO2通量晚约1 h，在14:00～15:00。不同生育期GPP日变化幅度存在明显差异，旺长期和成熟采收期GPP日变化幅度较大，日间峰值相对较高，分别为1.37 gC·m-2·h-1和1.36 gC·m-2·h-1；伸根期日间GPP峰值最低，仅0.75 gC·m-2·h-1，且峰值出现的时间也最晚，大约在16:00左右。

由图3c可知，烟田ET日动态呈“单峰”曲线变化规律。日出后，烟田ET缓慢增加，大约在15:00左右达到峰值，随后烟田ET逐渐下降，在夜间23:00达到最小值并趋于稳定。不同生育期烟田ET的日变化幅度虽然也存在差异，但除了旺长期日间ET在12:00～15:00期间明显高于其它生育期外，移栽还苗期、伸根期和成熟采收期ET的日变化曲线基本保持一致。

图3 CO2通量、GPP和ET的日变化特征Fig.3 Daily trend of CO2 fl ux, GPP and ET

2.2.2 水、碳通量季节变化特征

结合图4a和表2可知，研究区日均CO2通量具有显著的季节变化特征，整个大田生育期平均值为4.13 gC·m-2·d-1。其中，移栽还苗期碳吸收量相对较低，日均值为4.02 gC·m-2·d-1；伸根期碳吸收能力最弱，烟田整体呈碳排放状态，日均碳排放量为0.46 g C·m-2·d-1；进入旺长期后，烟田碳吸收能力逐渐增强，并在7月下旬即旺长期后期和成熟采收期前期达到峰值；烟田旺长期和成熟采收期碳吸收量日均值分别为4.27 gC·m-2·d-1和 6.51 gC·m-2·d-1。此后，随烟叶的成熟采收，烟田碳吸收量开始降低。

结合图4b和表2可知，研究区日均GPP呈 “双峰”曲线变化规律，整个大田生育期内平均值为10.79 gC·m-2·d-1。烟株移栽后，烟田CO2吸收量增加，GPP逐渐上升；进入伸根期后，烟田整体呈碳排放状态，该时期烟田GPP日均值最低，仅为6.63 gC·m-2·d-1；6月中旬进入旺长期后，烟田碳吸收量迅速上升，日均GPP迅速升高；进入成熟采收期后，GPP的日均值最高，达13.51 gC·m-2·d-1；9月份后，随烟叶采收，烟田日均GPP迅速下降。

结合图4c和表2可知，研究区日均ET整体呈现逐渐增加的趋势，整个大田生育期ET总量为586.26 mm，日均值为3.58 mm·d-1。在生育期初期，烟田ET日均值相对较低，日均值为3.04 mm·d-1；进入旺长期后，烟田ET日均值小幅度下降，但下降程度未达到显著性水平；7月下旬后，烟田进入成熟采收阶段，ET日均值在整个大田生育期最高，达4.36 mm·d-1，该时期ET总量占整个大田生育期总量的45.4%。

图4 CO2通量、GPP和ET日均值的季节变化特征Fig.4 Seasonal trends of daily CO2 fl ux, GPP and ET

表2 不同生育期CO2通量、GPP和ET日均值统计结果Tab.2 Statistic results of daily CO2 fl ux, GPP and ET in different growing stages

2.3 WUE变化特征

2.3.1 WUE日变化特征

由图5和表3可知，研究区烟田WUE的日变化动态呈先升高，后降低，最后又小幅度增加的规律。日出后，烟田WUE迅速升高并在8:00～9:00左右达到峰值；随后WUE逐渐下降，到14:00左右趋于稳定，16:00以后又有小幅度上升。整个大田生育期WUE日间变化范围为0.33～3.68 gC·kg-1H2O，变异系数为35.62%。其中，旺长期WUE日变化范围最大，变异系数为53.17%，日间最大值达8.61 gC·kg-1H2O，最小值为1.72 gC·kg-1H2O；伸根期WUE日变化范围为0.27～3.16 gC·kg-1H2O，变化范围最小，变异系数为41.00%。

图5 烟田WUE的日变化特征Fig.5 Daily trend of WUE

表3 WUE在不同生育期的统计结果Tab. 3 Statistic results of WUE in different growing stages

2.3.2 WUE季节变化特征

由图6和表4可知，烟田WUE与GPP的季节变化特征相对一致。整个大田生育期WUE日均值为3.6 gC·kg-1H2O，其中成熟采收期WUE日均值最高，为4.47 gC·kg-1H2O，伸根期WUE日均值最低，为2.04 gC·kg-1H2O。4月中旬，烟株移栽后，WUE缓慢增加；5月中旬，烟株进入伸根期，烟田WUE呈下降趋势；6月中旬以后，烟田WUE迅速升高，并在7月中旬到达峰值；7月中旬～8月下旬，烟田WUE相对稳定；进入9月份以后，烟田WUE开始降低。

图6 WUE日均值季节变化特征Fig. 6 Seasonal trend of daily WUE

表4 不同生育期WUE日均值统计结果Tab.4 Statistic results of daily WUE in different growing stages

2.4 烟田生态系统WUE的影响因素

前人研究表明，气温、净辐射和VPD对生态系统WUE影响显著[7,19-20]。对于水分状况较好的生态系统来说，生态系统WUE会随气温的升高而逐渐增加[32]。研究区烟田WUE与气温呈二次曲线相关（图7a），气温大约21℃是烟田WUE随气温变化的转折点，当气温低于21℃时，WUE与气温呈正相关，当气温超过21℃后，两者呈负相关。

净辐射是控制植物冠层与大气之间水碳交换的重要因素，生态系统通过净辐射的变化可以间接影响WUE[29]。由图7b可知，研究区烟田WUE与净辐射成二次曲线相关。随净辐射的增强，烟田生态系统WUE呈现先迅速增加后逐渐降低的变化规律，净辐射强度为200 W·m-2时是研究区烟田WUE随净辐射变化的转折点。

VPD作为气孔与外界进行水、碳交换的重要驱动力，对生态系统WUE影响显著[33]。 由图7c可知，研究区烟田WUE与VPD呈负指数相关，随VPD的升高，烟田WUE迅速下降，当VPD升高到大约2500 Pa后，烟田WUE逐渐趋于稳定。

图7 气温、净辐射和VPD对烟田生态系统WUE的影响Fig. 7 Relationships between air temperature, net radiation, VPD and WUE respectively

3 讨论

3.1 气温、净辐射和VPD对烟田WUE的影响机理

植物光合与蒸腾作用在群体和生态系统尺度上与气温呈显著正相关[34-35]，在土壤水分条件适宜的情况下，生态系统GPP和ET均随气温的升高而增加[36]。由图8a可知，烟田GPP与气温呈二次相关，而ET与气温呈正指数相关。气温大约25℃后，GPP增加的速率开始下降，当气温超过29℃后，GPP逐渐趋于稳定并出现小幅度降低，说明高温下烟株光合作用会受到抑制。结合图7a可知，当气温相对较低时，烟田GPP和ET均随气温的升高而增加，但GPP对气温变化的响应更为敏感，此时GPP的增加成为WUE提高的主导因素；当气温相对较高时，烟株光合作用相关酶代谢活性降低、碳吸收受到抑制，同时烟田土壤蒸发逐渐加剧，WUE开始下降。

前人研究表明，净辐射的变化直接作用于GPP与ET过程，从而间接影响WUE动态[29,37]。由图8b可知，净辐射与GPP呈二次曲线相关，而净辐射与ET呈线性正相关。在净辐射相对较低时，烟田GPP随净辐射的增强而逐渐增加，当净辐射强度升至500 W·m-2左右时，GPP趋于稳定；此后，随净辐射的进一步增强，GPP开始降低。此外，净辐射的增强与气温升高密切相关，两者导致ET呈线性增加。结合烟田WUE随净辐射的变化规律（图7b），随净辐射的增强和气温的升高，烟田水碳交换能力逐渐增强，GPP上升，WUE缓慢增加。当净辐射超过一定的强度后，烟株光合作用超过光饱和点出现光抑制，GPP降低，而此时烟田ET尤其是地面蒸发的线性增加成为WUE降低的主要原因。

由图8c可知，研究区烟田VPD与GPP、ET均呈二次曲线相关。在VPD小于2000 Pa时，GPP随VPD的升高而迅速增加；当VPD上升至2000 Pa左右时，GPP达到峰值；随着VPD的继续升高，GPP开始降低。而ET随VPD上升始终保持增加的趋势，但当VPD上升至2500 Pa后，其增加的速率显著降低，这与气孔关闭对蒸发的抑制有关[38]。叶片对CO2吸收速率随VPD的升高而增加[33]，但当VPD超过一定的阈值后，气孔关闭对光合碳吸收的调控成为限制GPP增加的主要因素[39]。综合来看，VPD对GPP和ET的影响存在敏感度差异，从而导致WUE与VPD呈负指数相关（图7c）。在VPD较低时，GPP对VPD变化的响应更为敏感，而ET的响应则相对迟缓，因此烟田WUE迅速升高；随VPD的继续升高，烟田光合作用过程受到限制，GPP开始下降，而ET仍继续增加，导致烟田WUE逐渐降低。

图8 气温、净辐射和VPD与烟田生态系统GPP、ET的关系Fig. 8 Relationships between air temperature, net radiation, VPD and GPP, ET respectively

3.2 不同时间尺度上烟田生态系统WUE影响因素分析

烟田WUE的日变化主要受到冠层导度的影响，烟田冠层导度体现了烟株与大气之间的水碳交换能力，较高的冠层导度意味着烟株与外界拥有较快的水气交换速率[40,19]。日出后，气温、净辐射和VPD逐渐升高，冠层导度也随之增加，烟田光合作用增强，GPP上升驱动WUE逐渐升高。当气温和光照强度达到一定的程度以后，光合作用产物受“光合午休”影响生产速率下降，但烟田土壤蒸发却随气温的升高而逐渐强烈，WUE受到限制开始降低。

烟田WUE的季节变化主要受烟田LAI的影响[41]。LAI的增加不仅增加了有效光合面积从而提高了生产力水平，还能通过增强冠层阻抗抑制土壤无效水分蒸发，从而降低蒸发占ET的比例，提高了WUE。烟株生育期初期，烟田LAI较低，地表裸露程度较大，土壤蒸发成为限制烟田WUE的主要因素。在伸根期，烟田的降水量较少，而净辐射强度和气温相对较高，烟株地上部分生长缓慢导致碳吸收能力较弱，而根系呼吸、光呼吸及土壤呼吸作用强烈，不利于光合产物的积累，整个该阶段烟田生态系统碳交换量较低、收支不平衡，呈碳排放状态。在旺长期，烟田的降水集中，日间辐射充足，烟田LAI迅速增加，烟株叶片光合作用和根系吸水能力随之提高，水碳代谢活跃，夜间温度较低导致了生态系统呼吸和夜间土壤水分蒸发的降低，烟田WUE逐渐升高并达到峰值。在成熟采收前期，虽然烟叶的采收降低了烟田LAI，但烟株中上部叶光合作用的自我调控使生态系统生产能力维持在较高的水平，导致烟田WUE在一段时间内保持较高的水平。但在成熟采收期后期，烟株中上部叶采收结束，烟田LAI下降导致地表裸露程度逐渐增加，土壤蒸发再次成为烟田WUE降低的主导因素。

4 结论

本研究基于涡度相关通量观测数据，对研究区烟田生态系统水、碳通量及WUE在不同时间尺度上的变化特征进行了分析，并结合气象因子、LAI和冠层导度探讨了引起烟田生态系统WUE变化的原因及相关变异机理，主要结论如下：

（1）烟田WUE的日变化为先迅速升高达到日间峰值，后逐渐下降，在傍晚前有一个小幅度增加的过程；烟田WUE的季节变化为生育初期较低，进入旺长期逐渐升高，到旺长期后期和成熟采收期前期阶段达到峰值，随后又逐渐下降。

（2）烟田WUE与气温和净辐射均呈二次曲线相关，与VPD呈负指数相关。GPP和ET两者随气温、净辐射和VPD变化的转折点时间的差异，导致了烟田WUE的时间动态特征。随气温、净辐射和VPD升高，烟株光合作用增强，GPP的提高成为烟田WUE增加的主导因素；当净辐射、气温和VPD达到一定的阈值后，烟田光合、蒸腾速率开始下降，而土壤蒸发的增加成为WUE降低的主导因素。

（3）叶片光合碳吸收能力和土壤蒸发是限制烟田WUE提高的主要因素。在土壤水分状况较好的条件下，提高烟田WUE可通过协调烟田个体和群体关系、提高冠层光能截获、结合烟株不同生育期需水规律降低无效蒸发等途径。

[1]蔡寒玉，汪耀富. 土壤水分对烤烟形态和耗水特性的影响[J]. 灌溉排水学报, 2005, 24 (1): 38-41.CAI Hanyu, WANG Yaofu. E ff ects of soil moisture on fl uecured tobacco leaf morphology and water consumption characteristics[J]. Journal of Irrigation and Drainage, 2005,24 (1): 38-41.

[2]KUGLITSCH FG, REICHSTEIN M, BEER C, et al.Characterisation of ecosystem water-use efficiency of European forests from eddy covariance measurements[J].Biogeosciences Discussions, 2008, 82(2): 204-12.

[3]段淑辉, 杨亿军, 刘建利, 等. 烟草需水规律研究进展[J].中国烟草科学, 2012, 33(4): 99-105.DUAN Shuhui, YANG Yijun, LIU Jianli, et al. Advance in water requirement of tobacco[J]. Chinese Tobacco Science,2012, 33(4): 99-105.

[4]刘贞琦, 伍贤进, 刘振业. 土壤水分对烟草光合生理特性影响的研究[J]. 中国烟草学报, 1995, 2( 3): 44-49.LIU Zhenqi, WU Xianjin, LIU Zhenye. The effects of soil water content on photosynthetic characteristics of tobacco[J]. Acta Tobacaria Sinica, 1995, 2(3): 44-49.

[5]SINCLAIR T R, BENNETT J M. Water use efficiency in crop production[J]. Biology Science, 1984, 34(1): 36-40.

[6]NIU S L, XING X R, ZHANG Z, et al. Water use efficiency in response to climate change: from leaf to ecosystem in a temperate steppe[J]. Global Change Biology, 2011, 17(2):1073-82.

[7]TONG XJ, LI J, YU Q, et al. Ecosystem water use efficiency in an irrigated cropland in the North China Plain[J]. Journal of Hydrology, 2009, 374(3): 329-37.

[8]谈晶晶, 邵孝侯, 陈立华, 等. 不同供水条件下烟田的水分利用效率[J]. 河海大学学报: 自然科学版, 2012, 04:446-450.TAN Jingjing, SHAO Xiaohou, CHEN Lihua, et al.Water use efficiency of tobacco fi eld under di ff erent water supply[J]. Journal of HoHai University: Natural Science Edition, 2012, 04: 446-450.

[9]ITO A, INATOMI M. Water-use efficiency of the terrestrial biosphere: a model analysis focusing on interactions between the global carbon and water cycles [J]. Journal of Hydrometeorology, 2012, 13(2): 681-94.

[10]REICHSTEIN M, DINH N, RUNNING S et al.Computation of European carbon balance components through synergistic use of CARBOEUROPE eddy covariance, MODIS remote sensing data and advanced ecosystem and statistical modeling[J]. Child & Adolescent Social Work Journal, 2003, (8): 12-39.

[11]PAPAL D, VALENTINI R. A new assessment of European forest carbon exchanges by eddy fl uxes and arti fi cial neural network spatialization[J]. Global Change Biology, 2002, 9:525-535.

[12]HASON P J, AMTHOR J S, Wullschleger S D, et al.Carbon and water cycle simulations for an upland oak forest using 13 stand-level models: inter-model comparisons and evaluations against independent measurements[J].Ecological Monographs, 2004, 74: 443-489.

[13] 朱治林, 孙晓敏, 袁国富, 等. 非平坦下垫面涡度相关通量的校正方法及其在China FLUX中的应用[J]. 中国科学: 地球科学, 2004, A(02): 37-45.ZHU Zhilin, SUN Xiaomin, YUAN Guofu, et al. Correction method of eddy covariance flux on non-flat underlying surface and its application in China FLUX[J]. Chinese Journal of Earth Science, 2004, A (02): 37-45.

[14]WILLIAMS CA, REICHSTEIN M, BUCHMANN N, et al. Climate and vegetation controls on the surface water balance: Synthesis of evapotranspiration measured across a global network of fl ux towers[J]. Water Resources Research,2012, 48(6): 327-40.

[15]BALDOCCHI DD, RYU Y. A synthesis of forest evaporation fluxes-from days to years-as measured with eddy covariance[M]. Springer Netherlands, 2011.

[16]FALGE E, BALDOCCHI D, OLSON R, et al. Gap fi lling strategies for defensible annual sums of net ecosystem exchange[J]. Agricultural & Forest Meteorology, 2001a,107(1): 43-69.

[17]FALGE E, BALDOCCHI D, OLSON R, et al. Gap fi lling strategies for long term energy fl ux data sets[J]. Agricultural& Forest Meteorology, 2001b, 107(1): 71-7.

[18]战领, 杨汉波, 雷慧闽. 基于通量观测数据的玉米水碳交换量及水分利用效率分析[J]. 农业工程学报, 2016,32(Supp.1): 88-93.ZHAN Ling, YANG Hanbo, LEI Huimin. Analysis of corn water consumption, carbon assimilation and ecosystem water use efficiency based on flux observations[J].Transactions of the Chinese Society of Agricultural Engineering, 2016, 32(Supp.1): 88-93.

[19]TANG XG, DING Z, LI HP, et al. Characterizing ecosystem water-use efficiency of croplands with eddy covariance measurements and MODIS products[J]. Ecological Engineering, 2015, 85: 212-217.

[20]ZHANG FM, JU WM, SHEN SH, et al. How recent climate change influences water use efficiency in East Asia[J].Theoretical and Applied Climatology, 2014, 116(1/2): 359-370.

[21] 戴衍晨, 王瑞, 张继光, 等. 恩施烟区烟田生态系统CO2通量变化特征分析[J]. 中国烟草学报, 2015, 21(4): 63-69.DAI Yanchen, WANG Rui, ZHANG Jiguang, et al.Characteristic analysis of CO2flux variation in Enshi tobacco-growing fi eld ecosystem[J]. Acta Tobacaria Sinica,2015, 21(4): 63-69.

[22]马君红. 四川省烤烟品质区划与风格特色定位研究[D].河南农业大学, 2014.MA Junhong. Quality division and style characteristic location of flue-cured tobacco of Si Chuan Province[D].Henan Agriculture University, 2014.

[23] 李春, 何洪林, 刘敏, 等. China FLUX CO2通量数据处理系统与应用[J]. 地球信息科学学报, 2008, 10(5): 557-565.LI Chun, HE Honglin, LIU Min, et al.. The design and application of CO2flux data processing system at China FLUX[J]. Geo-informatics on Science, 2008, 10(5): 557-565.

[24]BLACK TA, ETHIER GJ, DEWITT GB. Annual and seasonal variability of sensible and latent heat fl uxes above a coastal Douglas- fi r forest, British Columbia, Canada.E.R.Humphreys[J]. Agriculture Forest Meteorol, 2003, 115:109-125.

[25]WILSON KB, BALDOCCHI DD. Seasonal and inter-annual variability of energy fl uxes over a broad leaved temperate deciduous forest in North America[J]. Agricultural & Forest Meteorology, 2000, 100(1): 1-18.

[26]JANSSENS IA, LANKREIJER H, MATTEUCCI G, et al. Productivity overshadows temperature in determining soil and ecosystem respiration across European forests[J].Global Change Biology, 2001, 7(3): 269-78.

[27]SAIGUSA N, YAMAMOTO S, MURAYAMA S, et al.Gross primary production and net ecosystem exchange of a cool-temperate deciduous forest estimated by the eddy covariance method[J]. Agricultural & Forest Meteorology,2002, 112(3-4): 203-15.

[28]HIRANO T, HIRATA R, FUJINUMA Y, et al. CO2and water vapor exchange of a larch forest in northern Japan[J].Tellus Series B-chemical & Physical Meteorology, 2003,55(2): 244-57.

[29]REICHSTEIN M, FALGE E, BALDOCCHI D, et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm[J]. Global Change Biology, 2005, 11(9): 1424-39.

[30]LLOYD J, TAYLOR JA. On the Temperature Dependence of Soil Respiration[J]. Functional Ecology, 1994, 8(3): 15-23.

[31] 胡中民, 于贵瑞, 王秋凤, 等. 生态系统水分利用效率研究进展[J]. 生态学报, 2009, 29(3): 498-507.HU Zhongmin, YU Guirui, WANG Qiufeng, et al.Ecosystem level water use efficiency: A review[J]. Acta Ecologica Sinica, 2009, 29(3): 498-507.

[32] 李辉东, 关德新, 袁凤辉, 等. 科尔沁草甸生态系统水分利用效率及影响因素[J]. 生态学报, 2015, 35(2): 478-88.LI Huidong, GUAN Dexin, YUAN Fenghui, et al. Water use efficiency and its in fl uential factor over Horqin Meadow[J].Acta Ecologica Sinica, 2015, 35(2): 478-488.

[33]BATISTA KD, ARAJOW L, ANTUNES WC, et al.Photosynthetic limitations in coffee plants are chiefly governed by diffusive factors[J]. Trees, 2012, 26(2): 459-68.

[34]周海燕, 黄子琛. 不同时期毛乌素沙区主要植物种光合作用和蒸腾作用的变化[J]. 植物生态学报, 1996, 20(2):120-131.ZHOU Haiyan, HUANG Zichen. Changes in Photosynthesis and Transpiration of Main Constructive Plants in Ms Us Desert. Chinese Journal of Plant Ecology, 1996, 20(2): 120-131.

[35]MORIANA A, VILLALOBOS FJ, FERERES E. Stomatal and photosynthetic responses of olive (Olea europaea L.)leaves to water deficits[J]. Plant Cell and Environment,2002, 25(3): 395-405.

[36] 李辉东, 关德新, 袁凤辉, 等. 科尔沁草甸生态系统水分利用效率及影响因素[J]. 生态学报, 2015, 35(2): 478-88.LI Huidong, GUAN Dexin, YUAN Fenghui, et al. Water use efficiency and its in fl uential factor over Horqin Meadow[J].Acta Ecologica Sinica, 2015, 35(2): 478-488.

[37]刘晨峰, 张志强，孙阁，等. 基于涡度相关法和树干液流法评价杨树人工林生态系统蒸发散及其环境响应[J].植物生态学报, 2009, 33(4): 706-718.LIU Chenfeng, ZAHNG Zhiqiang, SUN Ge, et al.Quantifying evapotranspiration and biophysical regulations of a poplar plantation assessed by eddy covariance and sapfl our methods[J]. Chinese Journal of Plant Ecology, 2009,33(4): 706-718.

[38]KATUL GG, ELLSWORTH DS, LAI CT. Modelling assimilation and intercellular CO2from measured conductance: a synthesis of approaches[J]. Plant Cell &Environment, 2000, 23(12): 1313-28.

[39]MAROCO JP, PEREIRA JS, CHAVES MM. Stomatal responses to leaf-to-air vapour pressure de fi cit in sahelian species[J]. Functional Plant Biology, 1997, 24(3): 381-7.

[40] 林同保, 王志强, 宋雪雷, 等. 冬小麦农田二氧化碳通量及其影响因素分析[J]. 中国生态农业学报, 2008, 16(6):1458-1463.LIN Tongbao, WANG Zhiqiang, SONG Xuelei, et al. CO2flux and impact factors in winter wheat field ecosystem[J]. Chinese Journal of Eco-Agriculture, 2008,16(6): 1458-1463.

[41]胡中民. 中国草地样带典型生态系统水分利用效率的时空变异特征研究[D]. 中国科学院地理科学与资源研究所, 2008.HU Zhongmin. Spatiotemporal dynamics of ecosystem water use efficiency across China Grassland Transect[D].Institute of Geographic Sciences and natural resources research, Chinese Academy of Science, 2008.

:XU Tongqing, XU Yimin, WANG Chengdong, et al. Dynamic characteristics of water use efficiency in Panxi tobacco planting area based on fl ux observtion data [J]. Acta Tabacaria Sinica, 2017,23(3)

*Corresponding author.Email：taojian_ecology@163.com

Dynamic characteristics of water use ef fi ciency in Panxi tobacco planting area based on fl ux observation data

XU Tongqing1,2, XU Yimin1, WANG Chengdong1, LIU Guangliang1, WANG Songfeng1, DONG Jianxin1, WANG Shusheng1, TAO Jian1*
1 Tobacco Research Institute, Chinese Academy of Agricultural Sciences, Qingdao 266101, China;2 Graduate School, Chinese Academy of Agricultural Sciences, Beijing 100081,China

Dynamics of water-carbon flux and WUE of tobacco field and effect of biological factors were investigated based on flux observation data from Panxi tobacco planting area in 2015-2016. Results showed that tobacco fi eld WUE showed obvious diurnal and seasonal variations derived from dynamics of water-carbon flux. In daytime, WUE increased rapidly to maximum at 8:00-9:00, then decreased to minimum before evening. For the seasonal variation, the highest WUE was found in the stage of long term and mature harvest period. Statistical analysis indicated that there were quadratic relationships between WUE and air temperature as well as between WUE and net radiation, but there was an exponential decline relationship between WUE and vapor pressure de fi cit (VPD). The driving e ff ects of air temperature, VPD and net radiation on the increase of WUE showed temporal di ff erences. During growing period, VPD was the main driving factor on increasing WUE at the beginning, followed by air temperature and net radiation. The temporal di ff erences were mainly related to di ff erent sensitivity of gross primary production (GPP) and evapotranspiration (ET) to air temperature, VPD and net radiation.Therefore, to identify characteristics of WUE and its driving factors may signi fi cantly improve water resources management in tobacco fi elds of the study area.

water use efficiency; water-carbon fl ux; fl ux observation; tobacco fi eld eco-system

徐同庆，徐宜民，王程栋，等. 基于通量观测的攀西烟区烟田水分利用效率特征研究[J]. 中国烟草学报，2017, 23（3）

国家自然科学基金青年科学基金项目（41501054）；中国农业科学院烟草研究所青年科学基金项目（2015A02）；中国农业科学院科技创新工程（ASTIP-TRIC03）

徐同庆（1993—），研究生，主要研究方向：生态系统生产力与水分利用效率，Tel：0532-66715598，Email：xutongqing2015@163.com

陶 健（1983—），Tel：0532-66715598，Email：taojian_ecology@163.com

2017-01-05；< class="emphasis_bold">网络出版日期：

日期：2017-05-16