滩地杨树人工林皆伐后蒸发散与产流变化*

高升华，汤玉喜，唐 洁，唐学君，3，田 甜，易洪波，陈筑衢，张旭东**

(1.国家林业局林木培育重点实验室，林木遗传育种国家重点实验室，中国林业科学研究院林业研究所，北京 100091；2.湖南省林业科学院，湖南 长沙 410004； 3.国家林业局华东林业调查规划设计院，浙江 杭州 310019； 4.岳阳市君山区林业局，湖南 岳阳 414000)

滩地杨树人工林皆伐后蒸发散与产流变化*

高升华1，汤玉喜2，唐 洁2，唐学君1，3，田 甜1，易洪波4，陈筑衢4，张旭东1**

(1.国家林业局林木培育重点实验室，林木遗传育种国家重点实验室，中国林业科学研究院林业研究所，北京 100091；2.湖南省林业科学院，湖南 长沙 410004； 3.国家林业局华东林业调查规划设计院，浙江 杭州 310019； 4.岳阳市君山区林业局，湖南 岳阳 414000)

[目的]揭示杨树人工林皆伐对滩地蒸发散和产流的影响。[方法]基于涡度相关系统对长江滩地杨树人工林皆伐前后水汽通量连续3年(2010―2012年)的观测数据，通过对比皆伐前、后气候条件相似的2个时段(1整年)的蒸发散，揭示皆伐后研究区蒸发散的变化，并基于水量平衡反推研究区产流的变化。[结果]皆伐后土壤温度和水位上升，土壤表层含水量全年均值减小约0.03；研究区蒸发散皆伐前、后具有相似的日变化规律和季节动态特征，但皆伐后的全年蒸散量仅为皆伐前的66.3%；皆伐后研究区产流率(产流量/降雨量)从皆伐前的0.53上升至0.62；皆伐前、后7、8、12月的干旱指数(潜在蒸发散/降雨量)均大于1，其他时期均小于1。[结论]滩地杨树人工林皆伐后滩地蒸发散减少而产流率增加，加剧夏季干旱的可能性有所降低，但洪水爆发期间削减洪峰的能力也减弱。

杨树人工林；长江滩地；皆伐；蒸发散；涡度相关；产流

森林植被对径流的影响是森林水文学研究中争论的焦点[1]。尽管经历了上百年研究，森林植被变化对流域径流的影响至今仍未达成共识。最新研究表明，森林植被对径流的影响因流域水资源状况和流域特征(如面积、坡度、地表覆盖等)的不同而不同[2]；而人工林与天然林结构有很大差异，因而其水文效应与天然林也有巨大不同，因此，以天然林的研究结论来指导人工造林或恢复是不可取的[3]。目前，我国人工林面积已占全国森林总面积的1/3[4]，且人工林采伐与更新的频率远高于天然林；然而，目前关于人工林采伐与更新对流域水文过程影响的研究十分不足，从而限制了对我国林业生态工程水文效应的评估和人工林科学经营管理方案的制定。

蒸发散(ET)是森林生态系统水量平衡的重要组成部分，它占流域森林生态系统水分输出的40%～90%[3, 6]；然而，在传统的森林水文研究中，蒸发散往往仅作为降雨、径流和土壤储藏量的余项进行估算[6, 7]，而缺乏直接观测，这大大增加了森林水文效益研究的不确定性，也使蒸发散成为生态水文系统中最亟待研究的部分[8]。准确测算林地蒸发散的时空变化，对揭示森林影响径流的机理、开发流域水文模型和制定森林经营管理方案具有十分重要的意义[9]。近20年，涡度相关技术越来越多地被用于生态系统蒸发散的直接观测，而对人工林采伐前后的蒸发散观测却很少。

杨树是我国主要的造林树种之一[5, 10]，由于具有较短的轮伐期和较强的环境适应能力而被广泛用于长江中下游速生丰产林、滩地抑螺防病林和防护林的营建[11-13]。杨树人工林轮伐期一般10年左右，其采伐频率远高于杉木、马尾松等造林树种的人工林；然而，杨树人工林采伐更新对区域耗水和产流的影响尚未被揭示，使得难以准确评估长江中下游地区大面积营建杨树人工林所产生的区域水文效应。本研究基于涡度相关法对滩地杨树人工林皆伐前后生态系统蒸发散进行连续观测研究，旨在验证以下假设：(1)滩地杨树人工林皆伐后生态系统蒸发散减少；(2)杨树人工林皆伐后滩地产流率增加。

1 研究区概况

研究区位于湖南省岳阳市君山区长江外滩(29°31′35″ N，112°55′22″ E)，海拔31 m，系中亚热带向北亚热带过渡的气候区，具有典型的季风气候特征，春夏多雨，秋季多旱，冬季寒冷，全年盛行北风。长期年均气温16.8℃，年均降水量1 400 mm。土壤类型为江湖洲滩特有的潮土类型。

研究区面积60 hm2，呈南北向楔形，其中，东西向宽约800 m，南北长约1 500 m[14]，地势较平整。伐前林分营造于2000年，树种为美洲黑杨(PopulusdeltoidesMarsh)。该林分于2011年9月全部皆伐，皆伐前平均树高19.5 m，胸径23 cm，郁闭度达0.85，林下优势种为益母草(Leonurusarternisia(Laur.)S.Y.Hu)，其盖度达90%以上。新林于2012年2月初种植完成，树种仍为美洲黑杨，平均苗高3.5 m，胸径2.5 cm。2011年7月—2012年8月，林下植被主要有狗牙根(Cynodondactylon(L.)Pers.)、堇菜(ViolaverecundaA.Gray)、辣蓼(PolygonumflaccidumL.)、风轮菜(Clinopodiumchinense(Benth.) O. Ktze.)等。2010年淹水55 d(2010年6月28日—8月21日)，最大淹水深度2.8 m；2012年淹水40 d(2012年7月11日—8月20日)，最大淹水深度2.90 m；2011年未淹水。林地地表水随长江洪水退却而通过径流方式退出林地，2次退水后林地草本植物全部枯蔫并贴服地面。

2 研究方法

研究区建有30 m通量观测塔1座，以CSAT-3三维风速仪(Campbell Scientific, Inc., USA)和LI-COR 7500 CO2/H2O分析仪(LI-COR, Inc., USA)分别观测三维风速和CO2/水汽浓度。观测高度根据冠层高度进行调整，一般高于冠层高度2～3 m。采伐前观测高度为21.0 m；于2011年7月4日调整至3.5 m，此时通量塔以南的所有树木和通量塔的其它3个方向50 m以内的所有树木均皆伐完；2012年3月后观测高度为6.5 m。环境因子观测包括地面以上部分和地面以下部分，主要包括辐射、降雨、空气温湿度、土壤温度、土壤含水量、地下水位等[15]。其中，采用TE525雨量桶(Campbell Scientific, Inc., USA)观测降雨量，采用2个CS616土壤含水量传感器(Campbell Scientific, Inc., USA)观测土壤15 cm深处的土壤体积含水量(VWC，%)，用HOBO自记式水位温度记录仪(Onset Computer Corporation, USA)观测水位(WT，m)。

研究期间，在常规样地调查的基础上，以中分辨率卫星遥感(MODISNDVI)对植被进行动态监测，所用数据集为MYD13Q1，其时间和空间分辨率分别为16 d和250 m。数据来源于美国地质调查网站(http://www.usgs.gov/)。

采用EdiRe软件(version 1.5.0.32；Robert Clement，University of Edinburgh，Edinburgh，UK；http://www.geos.ed.ac.uk/abs/research/micromet/EdiRe/)进行通量计算和数据质量控制。主要过程包括：去除了超出仪器测量范围和超过4倍标准差的异常值，运用坐标轴二次旋转法使垂直风速均值为0[16]，进行超声虚温校准(SND correction)[17]，频率损失校准，WPL校准[18]等。

通量计算后，对30 min的通量数据进行质量控制：(1)去除夜间摩擦风速过低时的通量值[19]，2010年1月—2011年6月和2011年7月—2012年9月摩擦风速的阈值分别为0.14、0.10 m·s-1。(2)通过通量贡献率、湍流积分特征、稳态测试3项指标对30 min的通量数据进行综合评价，其中,当湍流积分特征或稳态测试值>100%或皆伐后通量贡献率<80%时，该30 min数据标记为不合格。在2011年5—8月皆伐期间，每天均根据采伐进度单独设定每个方向的通量风浪区长度。当研究区的通量贡献率<80%时，该30 min的记录将弃而不用。因此，5—6月的通量主要来源于未采伐区域，而7—8月通量主要来源于皆伐后的采伐迹地,最终水汽通量的数据保存率为55.4%。

以潜热通量插补后的数据估算累计蒸发散。对<2.5 h的空缺以线性内插法插补，而对≥2.5 h的空缺则采用滑动窗平均法进行插补[20]，窗口大小分别为7、15、30 d。

单日潜在蒸发散以世界粮农组织推荐的公式FAO Penman-Monteith进行计算：

(1)

式中：PET为潜在蒸发散(mm)；Δ为饱和水汽压随温度变化的斜率(kPa·℃)；Rn为净辐射(MJ·m-2)；G为土壤热通量；γ为干湿表常数(kPa·℃)；Ta为空气温度(℃)；es为饱和蒸汽压(kPa)；ea为实际蒸汽压(kPa)；u2为2 m 高的风速(m·s-1)；900为转换系数。基于式(1)计算的潜在蒸发散表示0.12 m高的草坪在不受水分限制时的最大蒸发散，而杨树人工林和采伐迹地的潜在蒸发散较草地更大，因此，需乘以转化系数K。由于伐前期8月土壤水分充足，且植被指数恢复较好，因此，取该月实际蒸发散与草地潜在蒸发散的比值作为K值，即1.42。采伐迹地在2012年6月25日—7月11日由于地下水位的抬升蒸发散有明显的增大，因此，将此阶段的实际蒸发散与草地蒸发散的比值作为采伐迹地的K值，即1.19，而淹水期间由于植被较少而不进行转换。

由于受长江水位的顶托作用，其生态系统水量平衡公式为：

ΔS=P-ET-(O-I)

(2)

式中：ΔS为生态系统中储水量的变化量(mm)；P为降雨量(mm)；ET为蒸发散(mm)；I和O分别为从长江流入系统内的水量和系统内流入至长江的水量，(O-I)为生态系统与长江的净交换水量(R，mm)。当R>0时，表示从生态系统流出的水量，即生态系统的产流量(WY)。ET通过涡度相关系统直接观测的水汽通量(g·m-2·s-1)获得，P通过雨量桶(TE52，Campbell Scientific, Inc., USA)进行计量。由于伐前林分结构较其前一年没有明显变化，因此,认为其全年ΔS≈ 0[2]，而采伐迹地ΔS采用式(3)进行计算[6]。

ΔS=ΔVWC15×300+ΔWT×θd

(3)

式中：ΔVWC15为地表15 cm深处土壤体积含水量的变化量，代表地面以下300 mm深的土层平均体积含水量变化量；ΔWT为平均地下水位的变化量(mm)；θd为土壤非毛细管孔隙度，取0.4。

选取2010年1月—2012年8月蒸发散和环境因子的连续观测数据进行采伐前、后的比较研究，其中，2010年1月—2011年6月为伐前期，2011年7月—2012年8月为采伐迹地期。在进行全年比较时，伐前期和采伐迹地分别选择2010年1月—12月(P1)和2011年9月―2012年8月(P2)，分别包含了伐前期和采伐迹地期的淹水月份(7―8月)。选取2010—2013年全年降雨量数据计算4年年均降雨量和单月平均降雨量以计算干旱指数。

3 结果与分析

3.1 环境因子动态

总体而言，2010年和2012年的区域气象条件相似，即雨水充沛，且具有相似的淹水情形(图1)。P2阶段的平均气温同该区域长期(1981—2010年)年均气温相当(16.8℃)，但较P1阶段低0.16℃。两阶段的降雨量总体无明显差异(约2 200 mm)，但均比该区域长期年均降水量(约1 400 mm)高出约800 mm(表1)。

P1和P2分别为选定的伐前一年时段和采伐迹地一年时段；虚线为采伐前和采伐后的分界线；阴影部分分别为采伐前和采伐后的淹水时段。P1 and P2 indicate a whole year before and after clearcut, respectively; the dash line is the border between before and after clearcut; the shaded parts are indicate the inundation periods before and after clearcut.图1 环境因子和蒸发散的动态变化Fig.1 Dynamics of environmental factors and ET

皆伐后土壤温度升高，且土壤含水量下降(表1)。尽管P2阶段平均气温较P1阶段更低，但平均土壤温度(17.85℃)却较P1阶段高0.52℃，主要由于皆伐后，地表植被覆盖减少导致土壤直接接受太阳辐射的能量增加。在总降雨量相当的情况下，P2阶段土壤含水量全年均值却较P1阶段减小3.11%，一方面由于地表覆盖的减少使地表温度升高和地表风速增大，从而加速了地表水分的蒸发；另一方面由于采伐后以草本植物为主，其根系分布较浅(<30 cm)，因而，对土壤表层的水分消耗更快。P2和P1阶段土壤温度和土壤含水量的变化与皆伐期间采伐迹地和未采伐区同步对照观测结果一致[12]。

森林采伐会导致地下水位的抬升[21]，而本研究区地下水位受地表植被和长江水位的共同调节。在皆伐期间(2011年7月前后)地下水位急剧抬升，但由于此时恰逢长江水位快速上涨，地下水位的抬升难以完全归因于皆伐；然而，此后水位一直较高直至次年春季到来(图1b)。此外，P2阶段的平均长江水位较P1阶段下降0.36 m，但平均地下水位却上升0.36 m，因此推测，皆伐导致地下水位的抬升。这是由于蒸发散减少导致地下水消耗减少(图1e和表2)，以及采伐迹地活体根系所触及的地下水深度较伐前减小，因而使深层次的水受植物蒸腾影响较小。

表1 2个阶段环境因子比较

皆伐前后林下植被的优势草本植物发生变化，且草本物种数量有所增加[12]。从2010年和2011年6—8月的NDVI数据推测，病虫害降低了伐前人工林的植被指数(图1d)，加之淹水导致伐前林下草本植物的死亡，使2011年7月NDVI降至全年最低水平(0.16)；淹水期间采伐迹地NDVI降低至全年最低(约0)。2012年5月末的除草作业使采伐迹地的NDVI大幅下降，但很快恢复，表明该区域草本植物生长旺盛。

3.2 蒸发散的动态变化

皆伐前后的蒸发散均呈现明显的季节变化，主要表现为生长季强而非生长季弱，而生长旺季明显下降(图1e)。P1阶段全年蒸发散最大月份为8月(183.0 mm)，5月次之(157.5 mm)，最小月份为1月和2月(约20 mm)。P2阶段全年最大月份为5月(151.0 mm)，7月次之(135.8 mm)，最小月份为1月(20.3 mm)。P1阶段淹水月份(7—8月)总蒸发散(293.6 mm)约为P2阶段(147.3 mm)的2.0倍。病虫害和淹水分别是采伐前和采伐后生长旺季ET下降的最主要原因。

在生长季和非生长季皆伐前后蒸发散的日变化特征(图2)均显著。在非生长季(图2a)，采伐前的人工林在2010年和2011年的变化规律相似，在下午14：00—15：00出现较弱的“午休”现象；采伐迹地(2012年)的蒸发散在白天较伐前更大，且“午休”现象不明显(图2a)。病虫害爆发前，在2010年和2011年的生长季, 伐前林的蒸发散较采伐迹地同期的更大，最大值分别为2012年的1.37和1.67倍，月总蒸散量分别为2012年同期的1.5和1.7倍(图2b)。在病虫害爆发的月份(图2c)，伐前林分蒸发散日最大值与采伐迹地相当，且2011年当月的月通量降至采伐迹地同期的1.1倍，此时为当年受病虫害影响最严重的月份(图1d)，说明病虫害对杨树人工林蒸发散具有重要影响。淹水期间(7—8月)伐前林日最大蒸发散略高于未淹水月份的采伐迹地(图2d)，但显著高于淹水期间的采伐迹地；总蒸发散分别为未淹水伐前林和淹水期采伐迹地的1.3和2.0倍。

图2 不同阶段蒸发散的典型月份平均日变化Fig.2 Mean diurnal variation of ET in typical months

3.3 皆伐对水量平衡的影响

皆伐导致研究区水量输出的结构发生变化(表2)。P1和P2阶段的总降雨量无明显差异，然而，P2阶段的蒸发散下降明显，仅为P1阶段的65.9%，而产流量上升为P1阶段的1.17倍。此外，由于皆伐后地下水位的上升和土壤含水量的下降(表1)，土壤储水量较P1阶段增加了134.7 mm。因而，蒸发散的减少一方面增加了地下水的储藏，另一方面增加了产流。根据所观测的蒸发散推算，P2阶段产流量较P1阶段增加190.8 mm，产流率(R/P)上升至P1阶段的1.18倍(表2)。

伐前林地淹水期间ET>P，导致R<0，表明淹水期间洪水贡献了约60 mm的蒸散量；而采伐迹地在淹水期间(2012年7—8月)的蒸散率(0.578)远比未淹水年份同期(2011年7—8月，0.785)的小，同时淹水期采伐迹地的植被指数(约为0)也远比未淹水年份同期(0.4～0.6)的小，因而,推测采伐迹地期的植被蒸腾作用在蒸发散中占较高的比例。

表2 各阶段的水量平衡

3.4 干旱指数的变化

PET/P常被用作干旱程度的指示参数[2]。相对于2010—2013年这4年的平均年降雨量(P_m)(1 863.8 mm)，P1和P2阶段的全年干旱指数(PET/P_m)分别为0.62和0.48；降雨量在一年内分布不均匀，其中前6个月占全年的66.4%，导致P1和P2阶段的PET/P_m前半年小于0.5(图3b)。伐前林和采伐迹地在7—8月和12月的干旱指数均大于1，都可能加剧区域干旱，而7—8月恰逢长江洪水期，因此，杨树人工林由于具有较大的蒸散耗水和较大的地下储水潜力而恰好起到削弱洪峰的作用。从P2阶段的PET/P_m看，采伐迹地在7月也可能起到削弱洪峰作用(图3b)，而实际上在淹水期间，植被全部被淹没，蒸发散远小于潜在蒸发散，因此，难以起到削洪作用(图3a)。相反，冬季草本植物依然比较旺盛(图1d)，蒸发散较伐前林高(图2a)，因此采伐迹地更可能促使南方冬旱(图3a)。

ET1和ET2分别为P1和P2阶段的蒸发散；PET1和PET2分别为P1和P2阶段的潜在蒸发散；P_m为2010—2013四年平均年降雨量。ET1 and ET2 are the ET in period P1 and P2, respectively; PET1 and PET2 are the potential ET in period P1 and P2, respectively; P_m was the mean yearly precipitation during 2010—2013.图3 降水的季节分布和干旱指数的季节变化Fig. 3 The mean monthly precipitation in the latest four years and dryness index seasonal variation

4 讨论

国外一些观测研究显示，森林采伐后生态系统的蒸发散将先减小，之后随着再造林或采伐迹地的自然演替，蒸发散将逐渐上升。Sun等[6]基于涡度相关系统对比观测了采伐迹地前2 a和火炬松幼林的蒸发散动态变化，发现采伐迹地的蒸发散较火炬松幼林小23%。Liu等[22-23]对火烧作业后的3、15和80 a的研究区进行对比观测，发现3 a的研究区比80 a的研究区蒸发散小33%，但15 a的蒸发散和80 a的蒸发散无明显差异，并预测蒸发散在最初几年会有所下降，但随着植被的恢复，中林龄蒸发散将会超过老龄林。此外，Chen等[24]对20、40和450 a的红杉林夏季的水汽通量对比观测时也发现，450 a的红杉林蒸散量比20、40 a林分的蒸散量分别高70%和12.2%；但也有观测结果显示，皆伐并未对蒸发散有明显改变，如Gholz等[21]对湿地松的采伐迹地、中龄林和成熟林进行对比观测发现，尽管叶面积指数和冠层结构有很大不同，3种林分的年蒸发散几乎相同。本站点11 a杨树人工林皆伐后，采伐迹地的蒸发散较伐前下降了34.1%，与Sun等[6]的研究结果最相似。根据上述研究结果推测，随着幼林根系和叶片的发展，这一比例将可能下降，但需要更长时间尺度的观测验证。

本研究显示，11 a杨树人工林皆伐后，产流效率上升至伐前的1.18倍，和目前较普遍观点相一致，即森林覆盖的减少有增加地表径流的趋势[3, 9]。此外，本站点7 a杨树人工林和11 a林的年蒸发散分别为958[25]和1 033 mm，比南亚热带(鼎湖山)季风阔叶林的(1 324 mm)小[26]，与同气候区的会同杉木林(848～983 mm)[27-28]十分接近，而略比同气候区千烟洲站的人工针叶林(640～820 mm)[29]高，且明显比北京地区11 a杨树人工林(585 mm)[25]高。因此，可以推测，在水分充足的长江中下游，杨树的耗水与同区域其他人工林的差异不大。因而，营造杨树人工林是否会导致区域干旱的问题归根结底是森林植被覆盖与径流关系的问题。Farley研究了26个配对集水区的植被对水的影响，发现在草地或灌木上造林将可能多耗150～200 mm或15%的降水[30]，而本研究中11 a杨树人工林较采伐迹地多消耗水190.8 mm，但仅多耗8.7%的降雨，主要由于本研究区降雨量较大。根据Zhou等[2]的研究，在湿润地区，植被的变化对产流率的影响更小，因而，在南方水分较充足的区域，植树造林对产流率的影响比在干旱区更小。

值得注意的是，滩地杨树人工林在7—8月长江丰水季由于具有更大的地下蓄水潜力和更大的蒸散量而比采伐迹地更能起到削减洪峰的作用；然而，它同时也暗示，在夏季雨水偏少的情况下，非长江滩地区域的杨树人工林比草地更有可能加剧区域季节性干旱。此外，造林、采伐、再造林对流域水文的影响并不完全相同[31]，需要更多针对性研究。

5 结论

滩地杨树人工林皆伐后，土壤温度有所上升，土壤含水量下降，且地下水位上升；皆伐前、后的蒸发散具有相似的日变化规律和季节动态特征，但采伐后，全年蒸散量减少至采伐前的66.3%；同时，产流率增加至采伐前的1.18倍。此外，采伐前、后的7、8、12月的干旱指数(PET/P)均大于1，而其它月份均小于1。本研究结果表明，滩地杨树人工林皆伐会导致流域蒸发散减少和流域产流率增加；由于丰水季节具有更大的地下蓄水潜力和更大的蒸散量，滩地杨树人工林在长江洪水爆发期间可更好的减少长江径流量而起到削减洪峰的作用，相反，当丰水季长江径流量较小时，杨树人工林相比采伐迹地或草地更可能加剧区域夏季干旱。

[1] 周延辉. 森林对径流影响的综述[J]. 地理科学进展, 1990,9(3):31-37.

[2] Zhou Guoyi, Wei Xiaohua, Chen Xiuzhi,etal. Global pattern for the effect of climate and land cover on water yield[J]. Nature Communications, 2015,6:5918.

[3] 刘世荣, 孙鹏森, 王金锡, 等. 长江上游森林植被水文功能研究[J]. 自然资源学报, 2001,16(5):451-456.

[4] 国家林业局森林资源管理司. 第八次全国森林资源清查主要结果(2009-2013年)[EB/OL]. http://www.forestry.gov.cn/main/65/content-659670.html.

[5] Sun G, Noormets A, Gavazzi M J,etal. Energy and water balance of two contrasting loblolly pine plantations on the lower coastal plain of North Carolina, USA[J]. Forest Ecology and Management, 2010,259(7):1299-1310.

[6] Wilson K B, Hanson P J, Mulholland P J,etal. A comparison of methods for determining forest evapotranspiration and its components: sap-flow, soil water budget, eddy covariance and catchment water balance[J]. Agricultural and forest Meteorology, 2001,106(2):153-168.

[7] Kool D, Agam N, Lazarovitch N,etal. A review of approaches for evapotranspiration partitioning[J]. Agricultural and Forest Meteorology, 2014,184:56-70.

[8] 张志强, 余新晓, 赵玉涛, 等. 森林对水文过程影响研究进展[J]. 应用生态学报, 2003,14(1):113-116.

[9] 国家林业局. 中国森林资源简况——第八次全国森林资源清查[G].国家林业局, 2013.

[10] 盛炜彤. 中国人工林及其育林体系[M]. 北京: 中国林业出版社, 2014: 533.

[11] 张旭东, 彭镇华, 周金星. 抑螺防病林生态系统抑螺机理的研究进展[J]. 世界林业研究, 2006,19(3):38-43.

[12] 高升华, 张旭东, 汤玉喜, 等. 滩地美洲黑杨人工林皆伐对地表甲烷通量的短期影响[J]. 林业科学, 2013,49(1):7-13.

[13] 彭镇华, 江泽慧. 中国新林种：抑螺防病林研究[M]. 北京: 中国林业出版社, 1995.

[14] 高升华, 张旭东, 汤玉喜, 等. 滩地人工林幼林不同时间尺度CH4通量变化特征——基于涡度相关闭路系统的研究[J]. 生态学报, 2016,36(18):5912-5921.

[15] 魏 远, 张旭东, 江泽平, 等. 湖南岳阳地区杨树人工林生态系统净碳交换季节动态研究[J]. 林业科学研究, 2010,23(5):656-665.

[16] Rebmann C, Kolle O, Heinesch B,etal. Data acquisition and flux calculations[M]//Aubinet M, Vesala T, Papale D. Eddy covariance: a practical guide to measurement and data analysis. London: Springer Atmospheric Sciences, 2012:59-84.

[17] Schotanus P, Nieuwstadt FTMf, De Bruin HAR. Temperature measurement with a sonic anemometer and its application to heat and moisture fluxes[J]. Boundary-Layer Meteorology, 1983,26(1):81-93.

[18] Webb E K, Pearman G I, Leuning R. Correction of flux measurements for density effects due to heat and water vapour transfer[J]. Quarterly Journal of the Royal Meteorological Society, 1980,106(447):85-100.

[19] Schmid H P. Ecosystem-atmosphere exchange of carbon dioxide over a mixed hardwood forest in northern lower Michigan[J]. Journal of Geophysical Research, 2003,108(D14).

[20] Falge E, Baldocchi D, Olson R,etal. Gap filling strategies for defensible annual sums of net ecosystem exchange[J]. Agricultural and forest meteorology, 2001,107(1):43-69.

[21] Gholz H L, Clark K L. Energy exchange across a chronosequence of slash pine forests in Florida[J]. Agricultural and Forest Meteorology, 2002,112(2):87-102.

[22] Liu H, Randerson J T, Lindfors J,etal. Changes in the surface energy budget after fire in boreal ecosystems of interior Alaska: An annual perspective[J]. Journal of Geophysical Research: Atmospheres (1984-2012), 2005,110(D13).

[23] Liu H, Randerson J T. Interannual variability of surface energy exchange depends on stand age in a boreal forest fire chronosequence[J]. Journal of Geophysical Research: Biogeosciences (2005-2012), 2008,113(G1).

[24] Chen J, Ustin S L, Suchanek T H,etal. Net ecosystem exchanges of carbon, water, and energy in young and old-growth Douglas-fir forests[J]. Ecosystems, 2004,7(5):534-544.

[25] Xiao J f, Sun G, Chen J,etal. Carbon fluxes, evapotranspiration, and water use efficiency of terrestrial ecosystems in China[J]. Agricultural and Forest Meteorology, 2013,182:76-90.

[26] 黄忠良, 孔国辉, 余清发, 等. 南亚热带季风常绿阔叶林水文功能及其养分动态的研究[J]. 植物生态学报, 2000,24(2):157-161.

[27] 康文星. 杉木人工林蒸发散及能量消耗规律的研究[J]. 中南林学院学报, 1993,13(1):74-80.

[28] 康文星, 邓湘雯, 赵仲辉. 林冠截留对杉木人工林生态系统物质循环的影响[J]. 林业科学, 2006,42(12):1-5.

[29] 魏焕奇, 何洪林, 刘 敏, 等. 基于遥感的千烟洲人工林蒸散及其组分模拟研究[J]. 自然资源学报, 2012,27(5):778-789.

[30] Farley K A, Jobbágy E G, Jackson R B. Effects of afforestation on water yield: a global synthesis with implications for policy[J]. Global Change Biology, 2005,11(10):1565-1576.

[31]Brown A E, Zhang L, McMahon T A,etal. A review of paired catchment studies for determining changes in water yield resulting from alterations in vegetation[J]. Journal of hydrology, 2005,310(1):28-61.

(责任编辑：詹春梅)

Variations of Evapotranspiration and Water Yield in a Post-clearcutting Riparian Poplar Plantation

GAO Sheng-hua1, TANG Yu-xi2, TANG Jie2, TANG Xue-jun1,3, TIAN Tian1, YI Hong-bo4, CHEN Zhu-qu4, ZHANG Xu-dong1

(1.Key laboratory of Tree Breeding and Cultivation, State Forestry Administration; State Key Laboratory of Tree Genetics and Breeding; Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China; 2.Hu’nan Academy of Forestry, Changsha 410004, Hu’nan, China; 3.East China Forest Inventory and Planning Institute, State Forestry Administration, Hangzhou 310019, Zhejiang, China; 4.Junshan Forestry Bureau of Yueyang City, Yueyang 414000, Hu’nan, China)

[Objective]To reveal the effects of clearcutting a poplar plantation at the floodplain of Yangtze River on evapotranspiration (ET) and water yield (WY). [Method] Based on 3-years’consecutive observation data of water vapor flux from a poplar plantation on the floodplain of Yangtze River during the years before and after clearcutting by eddy covariance system, the variation ofETpost-clearcutting was investigated by comparing theETs in two periods (a whole year before and after clearcutting, respectively) with similar weather conditions. Then the variation ofWYpost-clearcutting was calculated by water balance. [Result] (1) After clearcutting, both soil temperature and ground water level increased obviously while the soil water content decreased by about 0.03. (2) TheETsbefore and post clearcutting had the similar seasonal and diurnal dynamics; Nevertheless, the annualETpost-clearcutting was only 66.3% of that before clearcutting. (3) TheWYratio (WY/precipitation) increased from 0.53 to 0.62 after clearcutting. (4) The dryness index (potentialET/precipitation) was larger than 1 in July, August, and December, while less than 1 in other months in both before and post clearcutting period. [Conclusion] Clearcutting a poplar plantation on the floodplain of Yangtze River leads to a decrease ofETand an increase ofWYratio. As a result, it is less likely to intensify regional drought in summer, but the capacity of the floodplain for flood peak clipping during flooding would also decline.

poplar plantation；floodplain of Yangtze River；clearcutting；evapotranspiration；eddy covariance；water yield

10.13275/j.cnki.lykxyj.2017.03.017.03.018

2015-06-29

国家科技支撑课题“长江防护林质量调控与高效经营技术研究与示范(2015BAD07B04)”和“生态经济型血防林构建技术研究与示范(2015BAD07B07)”

高升华，男，中国林业科学研究院林业研究所，助理研究员.研究方向：林业生态工程与技术.E-mail: shenghua1128@yeah.net

S715

A

1001-1498(2017)03-0486-08

* 感谢盛炜彤先生和徐梅卿先生对论文提出的宝贵意见!

** 通讯作者：张旭东，研究员，博士生导师.E-mail: zhxdcaf@163.com