文章信息
- 王海珍, 韩路, 徐雅丽, 牛建龙, 于军
- Wang Haizhen, Han Lu, Xu Yali, Niu Jianlong, Yu Jun
- 灰胡杨叶片气孔导度特征及数值模拟
- Characteristics of Stomatal Conductance of Populus pruinosa and the Quantitative Simulation
- 林业科学, 2016, 52(1): 136-142
- Scientia Silvae Sinicae, 2016, 52(1): 136-142.
- DOI: 10.11707/j.1001-7488.20160116
-
文章历史
- 收稿日期:2014-11-04
- 修回日期:2015-10-21
-
作者相关文章
2. 新疆生产建设兵团塔里木盆地生物资源保护利用重点实验室 阿拉尔 843300
2. Key Laboratory of Protection and Utilization of Biological Resource in Tarim Basin, Xinjiang Production&Construction Groups Alar 843300
气孔是植物发育的具有复杂调节功能的器官(王玉辉等,2000),可调控植物体与外界环境之间的水气交换(郑云普等,2015),并随着所处的生境不同而发生变化,不仅是影响植物光合气体交换过程的关键要素,也是调控土壤-植被-大气连续体之间物质和能量交换的关键环节(唐凤德等,2008)。气孔开闭在某种程度上反映了植物的新陈代谢状况,其对环境响应的敏感度已成为评价植物抗旱性能的重要指标(司建华等,2008)。气孔调节在植物适应全球气候变化和逆境生存中扮演着关键角色,探讨植物气孔运动的外界环境因素和内在生理活动的影响及气孔开闭对生境变化的响应机制,有助于了解植物光合作用的限制因素、水分利用效率以及拟合作物生产力。同时,气孔导度模拟与模型构建可为揭示植物与大气间物质循环、能量流动过程和植被对环境变化的生理生态适应机制提供理论依据,从而进一步深入阐明植物气体交换特性与气孔调控机制(鱼腾飞等,2012),对构建全球不同类型生态系统的水碳循环耦合模型、促进植被恢复与生态环境改善具有重要意义。
灰胡杨(Populus pruinosa)是极端干旱荒漠区特有的珍稀渐危木本植物,集中分布于我国新疆塔里木河流域的上游区域,与胡杨(P.euphratica)、柽柳(Tamarix ramosissima)等植物构成独特的荒漠河岸林景观,在防风固沙、防治荒漠化和维系绿洲生态安全、保障社会经济持续发展中发挥着重要作用。近60年以来,受全球气候变暖和南疆大规模毁林开荒、大水漫灌等人类活动的影响,塔里木河河水径流量持续减少,地下水埋深不断降低,荒漠生态系统格局与过程被改变,植被稀疏,长势衰败,土地沙化,因此,挽救塔里木河流域极端干旱荒漠区的灰胡杨林迫在眉捷。目前,国内外对胡杨的生理生态研究较多(Gries et al., 2003; Chen et al., 2011; Overdieck et al., 2013),主要集中于光合作用动态过程、影响因素(于军等,2008)、光响应参数拟合(伍维模等,2007)、蒸腾耗水特性(王海珍等,2010)及生理生态适应机制等方面(Bogeat-Triboulot et al., 2007)。尽管已有研究对灰胡杨的光合气体交换特性和生理生态适应机制有所阐述,但对灰胡杨气孔导度的数值模拟尚未见报道。鉴于此,本文以塔里木盆地灰胡杨为研究对象,根据2年野外气体交换参数的实测数据,探讨其气孔导度对干旱荒漠环境的生态响应方式,并以国际上2类典型的气孔模型模拟和建立不同时段灰胡杨叶片的气孔导度模型,以期为构建适合极端干旱荒漠区植物水碳循环耦合模型及揭示荒漠植物与大气间物质、能量的交换机制提供理论依据。
1 材料与方法 1.1 试验材料在新疆叶尔羌河下游阿瓦提县丰收三场(80°25' E,39°40' N)大面积集中分布的天然灰胡杨林内选择地势平坦、长势良好的灰胡杨样方(20 m×20 m),测定每木胸径等指标。根据灰胡杨径级分布选出无病虫害、枝叶繁茂的5~8株作为标准木。
1.2 试验方法于2012年7—9月、2013年6—9月中旬晴朗天气下,利用光合测定仪(LI-6400,USA)测定标准木树冠外围小枝上位置相同成熟叶的气孔导度(Gs)、净光合速率(Pn)和环境因子日动态,每株标准木测定3个叶片,从8:00开始每隔2 h测定1次相关气体交换参数,于20:00结束。
1.3 数据处理采用Orign8.0软件作图,并利用Statistic统计软件进行参数估算。模型拟合效果依据决定系数(R2)和均方根误差(root mean sguare error,RESM)(钟楚等,2013)来判定。
2 结果与分析 2.1 灰胡杨叶片气孔导度动态植物叶片气孔导度(Gs)对诸如光合有效辐射、温度、水汽压亏缺及环境中CO2浓度等环境因子反映十分敏感(Kim et al., 1991)。由图 1可知,不同年份生长季各月灰胡杨叶片Gs的日变化均呈单峰曲线,随着光合有效辐射(PAR)增加,Gs逐渐上升,于10:00—12:00达到1天中的最高峰,随后一直持续下降。2012,2013年中除9月Gs峰值均出现在上午10:00外,其余月份Gs峰值出现时间不同。而且不同月份Gs也不同,2012年日平均Gs依次为9月>7月>8月,2013年则为8月>9月>7月>6月; 9月Gs日变幅最大(37.01%),6月日变幅最小。可见,晴朗天气中午由于PAR、大气温度(Tair)和叶内外水汽压差(VPD)都达到了1天内的较高值,此时气孔适时关闭可避免植物体内水分大量散失和维持木质部生理活动。这表明在极端干旱荒漠区,灰胡杨气孔的自我调控是对极限逆境的一种自我调整和适应。
1)光合有效辐射(PAR)H 由图 2A可知,灰胡杨Gs随PAR升高而增大,二者呈幂函数关系。PAR<100 μmol ·m-2s-1时Gs上升较快,PAR>100 μmol ·m-2s-1时上升速度减慢,表明高PAR对灰胡杨Gs有制约作用。
2)空气温度(Tair)H 由图 2B可知,灰胡杨Gs与Tair呈二次曲线关系。Gs随Tair升高逐渐增大,在25 ℃左右达到最大,此后随Tair升高Gs逐渐减小,表明灰胡杨可通过降低Gs开张程度来避免高温伤害与树体水分损失,以适应极端干旱荒漠生境。
3)相对湿度(RH)H 由图 2C可知,灰胡杨Gs与RH呈二次曲线关系。清晨RH较高,Gs随RH上升而增大,在RH为53%左右达到最大,此后随RH上升逐渐下降。维持较大气孔导度的RH在30%~55%之间,过高与过低均会降低Gs,这是由于RH受PAR和Tair的控制而引起Gs变化。
4)水汽压亏缺(VPD)H VPD是温度和空气湿度的综合反应(李炜等,2012)。由图 2D可知,灰胡杨随VPD的升高而增大,之后又随VPD升高逐渐降低。这是因为低VPD下灰胡杨蒸腾耗水相对较弱,水分供给基本能满足植物水分需求,未明显引起保卫细胞失水收缩; 而高VPD下因蒸腾失水过多,
水分供应不足引起保卫细胞大量失水,导致Gs逐渐关闭。过高VPD下,植物气孔的适时关闭是减少水分利用数量、提高水分利用效率和维持木质部生理活动的关键行为(李炜等,2012)。
5)大气CO2 浓度(Ca)H 由图 2E可知,灰胡杨Gs随Ca的升高而增大,至400 μmol ·mol-1左右开始逐渐降低,二者呈二次曲线关系。
2.3 环境影响因子分析根据2年野外测定的Gs,PAR,VPD,Tair,RH和大气CO2浓度(Ca)数据,采用逐步回归方法寻找影响灰胡杨叶片Gs的主要环境因子。由表 1可知,灰胡杨叶片气孔导度的回归方程表现出极显著的回归关系,2年中VPD,Tair,PAR出现频率较高,而RH,Ca仅出现3次。全天、上午和下午不同时段各环境因子对Gs的影响程度不同,全天Gs主要由VPD,Tair,PAR决定,上午主要由VPD,Tair决定,而下午则由VPD,Tair,RH决定。总体来看,整个生长季6—9月,PAR,VPD,Tair 3个环境因子对Gs的影响最为显著,而Ca,RH相对较弱。
影响植物Gs的环境因子众多,作用不同。目前,国际上有2类构建机制不同的代表性Gs模型(Jarvis,1976; Ball et al., 1987; Leuning,1995)。
1)Jarvis气孔导度模型H
影响荒漠河岸林灰胡杨叶片Gs的主要环境因子为PAR,Tair和VPD(表 1),因此,众多学者通过科学研究建立了Gs与环境因子(PAR,VPD,Tair)之间的数量关系(Kim et al., 1991; Turner et al., 1984; Hofstra et al., 1969),并将其应用到Jarvis气孔导度模型(Jarvis,1976)中,形成了Gs对环境因子的响应数学模型:
${G_s} = PAR{{{\alpha _3}T_{air}^2 + {\alpha _4}{T_{air}} + {\alpha _5}} \over {\left( {{\alpha _1} + PAR} \right)\left( {{\alpha _2} + VPD} \right)}}$
(1)
将野外测定的相关数据(PAR,VPD,Tair,Gs)代入方程(1),拟合参数(表 2)与构建模型如下:
全天:${G_s} = PAR{{0.001T_{air}^2 + 0.013{T_{air}} - 0.090} \over {\left( {260.443 + PAR} \right)\left( { - 0.219 + VPD} \right)}}$
上午:${G_s} = PAR{{0.001T_{air}^2 + 0.014{T_{air}} - 0.173} \over {\left( {246.369 + PAR} \right)\left( { + 0.256 + VPD} \right)}}$
下午:${G_s} = PAR{{0.006T_{air}^2 + 0.370{T_{air}} + 6.582} \over {\left( {134.766 + PAR} \right)\left( { - 0.949 + VPD} \right)}}$
该模型对全天、上午和下午时段灰胡杨叶片气孔导度变异程度的总体解释能力分别为69.1%,62.2%和63.3%,参数估算时模拟值与实测值的相关系数分别为0.829 6,0.799 8和0.550 7,达到极显著水平(P<0.01)。这说明非线性气孔导度模型能够较好说明荒漠河岸林灰胡杨叶片气孔导度与环境因子之间的动态关系。
2)Leuning-Ball气孔导度模型Ball等(1987)研究指出Gs与Pn呈线性相关,提出了代表性的线性Gs模型:
${G_s} = m{{{P_n} \times RH} \over {{C_a}}} + b$。
Leuning(1995)指出Ball方程存在的缺陷,认为Gs对VPD的响应灵敏性强于对RH,因此改进Ball模型用Ca-Γ (CO2补偿点)代替Ca,改进后的Leuning-Ball气孔导度模型为:
${G_s} = m{{{P_n}} \over {\left( {{C_a} - \Gamma } \right)\left( {1 + VP{D_s}/VP{D_0}} \right)}} + {G_{s0}}$
基于王海珍等(2015)估算灰胡杨叶片在25 ℃时的Γ 值和2年野外测定的PAR,VPD,Gs,Ca,Pn数据,代入Leuning-Ball
气孔导度模型,参考唐凤德等(2008)方法估算出相关参数(表2)。不同时段灰胡杨叶片Gs模型为:
全天:${G_s} = 161.125{{{P_n}} \over {\left( {{C_a} - \Gamma } \right)\left( {1 + VP{D_s}/0.066} \right)}} + 0.099$
上午:${G_s} = 246.355{{{P_n}} \over {\left( {{C_a} - \Gamma } \right)\left( {1 + VP{D_s}/0.035} \right)}} + 0.139$
下午:${G_s} = 85.694{{{P_n}} \over {\left( {{C_a} - \Gamma } \right)\left( {1 + VP{D_s}/0.073} \right)}} + 0.162$
该模型对全天、上午和下午时段灰胡杨叶片气孔导度变异程度的总体解释能力分别为53.5%,30.6%和44.5%,参数估算时模拟值与实测值的相关系数分别为0.730 7,0.797 7和0.672 1,达极显著水平(P<0.01)。从模型解释能力来看,Leuning-Ball气孔导度模型拟合下午时段灰胡杨叶片Gs的精准度明显高于上午时段,由于此模型中引入了RH和Ca,与逐步回归分析结果一致。
综合比较上述模型,不同时段Jarvis气孔导度模型的解释能力均优于Leuning-Ball气孔导度模型。H 2.5 2类代表性模型精度对比H W S 由图 3A可知,Jarvis气孔导度模型模拟未参与模型构建的灰胡杨Gs实测数据的模拟值与实测值的相关系数r(0.830 4)、直线斜率b(0.670 5)均大于Leuning-Ball气孔导度模型的r(0.730 7)和b(0.533 8)。同时,利用Jarvis和Leuning-Ball气孔导度模型对灰胡杨叶片Gs日动态进行模拟计算,Jarvis气孔导度模型的模拟值与实测值较接近且日变化规律更相似(图 3B),表明Jarvis非线性气孔导度模型优于Leuning-Ball线性气孔导度模型,这是由于气孔导度对环境因子的响应是一种非线性的(表 1)。从不同时段2类模型Gs拟合值的均方根误差(RMSE)(表 3)可看出,Jarvis气孔导度模型拟合全天、上午与下午时段的效果均优Leuning-Ball气孔导度模型,但下午时段Leuning-Ball气孔导度模型拟合效果优于上午时段,这与不同时段Gs拟合模型的R2表现一致。此外,Gs和Pn间呈线性关系这一试验条件是Ball模型构建的前提(Ball et al., 1987)。本研究中灰胡杨叶片Gs与Pn间呈极显著正相关,但二者呈二次曲线关系(图 4),这与Ball模型构建机制不一致,从而导致Ball模型拟合效果不如Jarvis模型。
1)不同年份生长季各月灰胡杨叶片Gs日变化均呈单峰曲线,Gs随着环境因子的变化而变化,10:00—12:00达峰值。不同月份Gs峰值出现时间与大小存在明显差异,日平均Gs依次为8月>9月>7月>6月(2013年),这与塔里木河洪水发生时间、低水多变的生境有关。
2)荒漠河岸林灰胡杨叶片Gs与环境因子关系较密切,逐步回归分析表明,PAR,VPD,Tair3个环境因子对Gs的影响最为显著,其次为Ca与RH。上午、下午不同时段影响灰胡杨Gs的主要环境因子不同,上午Gs主要由VPD,Tair,PAR决定,下午还受Ca,RH的影响。
3)利用Jarvis模型和Leuning-Ball模型拟合上午、下午不同时段灰胡杨Gs动态,结果表明Jarvis模型拟合上午、下午时段灰胡杨Gs的效果相似; 但Leuning-Ball模型拟合下午时段灰胡杨Gs的效果优于上午,因为上午、下午不同时段影响Gs的环境因子不同及此模型中引入了RH或Ca。
4)拟合效果表明,以Jarvis非线性模型为模板建立的不同时段灰胡杨气孔导度模型优于Leuning-Ball线性模型,这是因为灰胡杨Gs和Pn间呈二次曲线关系,与Ball模型构建的前提条件不一致导致的。塔里木极端干旱荒漠区灰胡杨叶片气孔导度的最优响应模型为:Gs=PAR(0.001Tair2+ 0.013 Tair -0.090)/[(260.443+PAR)(-0.219+VPD)]。
[1] |
李炜, 司建华, 冯起, 等. 2013. 胡杨(Populus euphratica)蒸腾耗水对水汽压差的响应. 中国沙漠, 33(5):1377-1384. (Li W, Si J H, Feng Q, et al. 2013. Response of transpiration to water vapour pressure differentials of Populus euphratica. Journal of Desert Research, 33(5):1377-1384.[in Chinese])(2) |
[2] |
司建华,常宗强,苏永红,等. 2008. 胡杨叶片气孔导度特征及其对环境因子的响应. 西北植物学报, 28(1):125-130. (Si J H, Chang Z Q, Su Y H, et al. 2008. Stomatal conductance characteristics of Populus euphratica leaves and response to environmental factors in the extreme arid region. Acta Boreali-Occidentalia Sinica, 28(1):125-130.[in Chinese])(1) |
[3] |
唐凤德, 武耀祥, 韩士杰, 等. 2008. 长白山阔叶红松林叶片气孔导度与环境因子的关系. 生态学报, 28(11):5649-5655. (Tang F D, Wu Y X, Han S J, et al. 2008. Relationship of stomatal conductance of leaf with environmental factors in broad-leaved Korean pine forest at Changbai Mountain. Acta Ecological Sinica, 28(11):5649-5655.[in Chinese])(2) |
[4] |
王玉辉, 周广胜. 2000. 羊草叶片气孔导度对环境因子的响应模拟. 植物生态学报, 24(6):739-743. (Wang Y H, Zhou G S. 2000. Analysis and quantitative simulation of stomatal conductance of Aneurolepidium chinense. Acta Phytoecologica Sinica, 24(6):739-743.[in Chinese])(1) |
[5] |
王海珍, 韩路, 徐雅丽, 等. 2010. 控水处理对胡杨、灰胡杨生长特性及水分利用效率的影响. 农业工程学报, 26(增2):195-201. (Wang H Z, Han L, Xu Y L, et al.2010. Effects of water stress on growth characteristics and water use efficiency of Populus euphratica and Populus pruinosa. Transactions of the CSAE, 26(Supp. 2):195-201.[in Chinese])(1) |
[6] |
王海珍, 韩路, 徐雅丽, 等. 2015. 不同温度下灰胡杨叶片气孔导度对光强响应的模型分析. 生态环境学报, 24(5):741-748. (Wang H Z, Han L, Xu Y L, et al.2015. Model analysis of the stomatal conductance response to light in Populus pruinosa at different temperatures in the Taklimakan desert. Ecology and Environmental Sciences, 24(5):741-748.[in Chinese])(1) |
[7] |
伍维模, 李志军, 罗青红, 等. 2007. 土壤水分胁迫对胡杨、灰叶胡杨光合作用-光响应特性的影响. 林业科学, 43(5):30-35. (Wu W M, Li Z J, Luo Q H, et al. 2007. Effects of soil water stress on light response curves of photosynthesis of Populus euphratica and Populus pruinosa. Scientia Silvae Sinicae, 43(5):30-35.[in Chinese])(1) |
[8] |
于军, 王海珍, 周正立, 等. 2008. 塔里木荒漠优势树种气体交换特性与环境因子的关系研究 西北植物学报, 28(10):2110-2117. (Yu J,Wang H Z, Zhou Z L, et al. 2008.Gas exchange characteristics of desert dominant plants and relationships between them and environmental factors in Tarim basin. Acta Boreali-Occidentalia Sinica, 28(10):2110-2117.[in Chinese]).(1) |
[9] |
鱼腾飞, 冯起, 司建华. 2012. 极端干旱区多枝柽柳叶片气孔导度的环境响应模拟. 植物生态学报, 36(6):483-490. (Yu T F, Feng Q, Si J H. 2012. Simulating responses of leaf stomatal conductance to environmental factors for Tamarix ramosissma in an extreme arid region of China. Chinese Journal of Plant Ecology, 36(6):483-490.[in Chinese])(1) |
[10] |
钟楚, 朱勇. 2013. 烟草气孔导度对光强的响应. 中国生态农业学报, 21(8):966-972. (Zhong C, Zhu Y. 2013.Response of stomatal conductance to light in tobacco plants. Chinese Journal of Eco-Agriculture, 21(8):966-972.[in Chinese])(1) |
[11] |
郑云普, 徐明, 王建书, 等. 2015. 玉米叶片气孔特征及气体交换过程对气候变暖的响应 作物学报, 41(4):601-612. (Zheng Y P, Xu M, Wang J S, et al. 2015. Responses of the stomatal traits and gas exchange of maize leaves to climate warming. Acta Agronomica Sinica, 41(4):601-612.[in Chinese]).(1) |
[12] | Ball J T, Woodrow I E, Berry J A. 1987. A model predicting stomatal conductance and its contribution to the control of photosynthesis under different environmental conditions//Biggins J. Progress in Photosynthesis Research. Netherland:Martinus Nijhoff Publishers, Dordrecht, 4:221-224.(3) |
[13] | Bogeat-Triboulot M B, Brosche M, Renaut J, et al. 2007. Gradual soil water depletion results in reversible changes of gene expression, protein profiles, ecophysiology, and growth perforrnance in Populus euphratica, a poplar growing in arid regions. Plant Physiology, 143(2):876-892.(1) |
[14] | Chen Y P, Chen Y N, Xu C C, et al. 2011. Photosynthesis and water use efficiency of Populus euphratica in response to changing groundwater depth and CO2 concentration. Environmental Earth Sciences, 62(1):119-125.(1) |
[15] | Gries D, Zeng F, Foetzki A, et al. 2003. Growth and water relations of Tamarix ramosissima and Populus euphratica on Taklamakan desert dunes in relation to depth to a permanent water table. Plant, Cell and Environment, 26(5):725-736.(1) |
[16] | Hofstra G, Hesketh J D. 1969. The effect of temperature on stomatal aperture in different species. Canadian Journal of Botany, 47(8):1307-1310.(1) |
[17] | Jarvis P G.1976. The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field. Philosophical Transactions of the Royal Society of London B:Biological Sciences, 273(2):593-610.(2) |
[18] | Kim J, Verma S B.1991. Modeling canopy stomatal conductance in a temperature grassland ecosystem. Agricultural and Forest Meteorology, 55(1/2):149-166.(2) |
[19] | Leuning R. 1995. Critical appraisal of a combined stomatal-photosynthesis model for C3 plants. Plant, Cell and Environment, 18(4):339-355.(2) |
[20] | Overdieck D, Ziche D, Yu R D. 2013. Gas exchange of Populus euphratica leaves in a riparian zone. Journal of Arid Land, 5(4):531-541.(1) |
[21] | Turner N C, Schulze E D, Gollan T. 1984.The response of stomata and leaf gas exchange to vapour pressure deficits and soil water content. Oecologia,63(3):338-342.(1) |