黄土丘陵区沙棘、油松和刺槐光合生理特性及其环境适应性

吴 旭, 唐亚坤, 陈 晨, 贾 畅, 陈云明1,,*

1 中国科学院水利部水土保持研究所黄土高原土壤侵蚀与旱地农业国家重点实验室, 杨凌 712100 2 中国科学院大学,北京 100049 3 西北农林科技大学水土保持研究所, 杨凌 712100 4 西北农林科技大学林学院, 杨凌 712100

光合作用是影响植物生长发育的重要生理过程[1]。光照和水分作为光合作用的能量来源和原料,是影响植物物质累积和代谢过程的重要生态因子,植物对光强和水分的响应机制一直是植物生理生态学研究热点和难点问题[2-4]。光响应曲线反映光合速率与光强的相关关系,通过光响应模型拟合的光合参数,可用来判断植物光合能力及其受环境变化的影响程度。因此,通过对光响应曲线的测定及模拟,是研究不同植物光合特性及其环境的适应性的重要手段。

全球气候变化导致极端事件多发,其中干旱和强光等环境变化对植物生长发育的胁迫和制约作用更为突出。植物受光胁迫时,光合作用的光抑制造成约10%的碳损失,对植物的生存具有显著影响[5]。此外,水分胁迫可以显著降低与植物光合作用有关代谢生理过程的活性,影响植物正常生长[6]。目前有关光合特性的研究主要从光响应模型的选择和水分胁迫等角度开展。

常用的光响应模型主要有直角双曲线(Rectangular hyperbola,RH)[7]、非直角双曲线(Non-rectangular hyperbola,NRH)[8]、指数(Exponential,EM)[9]、直角双曲线修正(Modified rectangular hyperbola,MRH)[10]和指数修正模型(Modified exponential,MEM)[11]等。由于植物生理特性差异,不同植物适宜模型不同。其中,刘强等[12]研究表明,MRH适宜长白落叶松(Larixolgensis)。肖丹丹等[13]认为EM为榆属植物最优拟合模型。此外,水分胁迫对植物光响应曲线影响较大,水分过低或过高均会降低植物光合能力,且不同光响应模型的适用性也随之发生改变[14]。郎莹等[15]研究发现,MRH能较好地拟合各种水分下山杏(Armeniacasibirica)光响应过程及其生理参数;RH和NRH仅适用于正常水分条件下。而周多多等[16]认为水分有利时胡杨(Populuseuphratica)应用MRH和NRH较好;水分亏缺条件下采用MRH和RH更为适合。因此,筛选不同水分条件下植物适宜的光响应模型,对于探究不同植物光合适应性及其环境适宜性具有重要意义。

然而,不同植物对环境的响应与适应特性存在较大差异[17-18]。为适应不同的环境,植物在光补偿点(LCP)、光饱和点(LSP)、暗呼吸(Rd)、光合能力(Pnmax)等光合生理特性上会发生相应变化。黄土丘陵区退耕还林以来营造了大面积的沙棘(Hippophaerhamnoides)、油松(Pinustabuliformis)和刺槐(Robiniapseudoacacia)等人工纯林,由于树种单一,加上集中连片分布,导致林地地力衰退,林分稳定性降低。相比之下,混交林可明显提高林地生产力和林分稳定性[19]。经过长期的造林实践,沙棘不但能够适应恶劣的自然环境,而且具有较强固氮能力,能为其他植物生长提供适宜的养分条件,故多与其他树种混交种植[20]。目前,已有关于沙棘的研究主要侧重于沙棘与其混交树种的水文效应及水土保持功能[21-22],而关于沙棘、油松和刺槐3个树种光合生理生态的研究也仅仅局限在纯林[23-25]。此外,相比纯林,混交林可以显著改变林内光照和水分等环境条件,最终可能导致林内各组成树种光合生理特性的差异。因此,有必要对不同种植模式条件下,植物光合生理特性对水分环境变化的响应进行进一步探究,旨在为科学评价黄土高原地区人工林光能利用稳定性和可持续提供依据。

本文以沙棘、油松和刺槐为对象,选取典型干湿季,测定其在不同水分条件下的光响应曲线,采用RH、NRH、EM、MRH及MEM模型对其进行拟合分析。主要目标：(1)明确不同水分条件下3个树种的光响应过程及其适宜拟合模型;(2)探讨沙棘、油松和刺槐在不同种植模式下的光合生理生态适应性。以期为黄土丘陵区植被恢复与生态重建提供理论依据。

1 材料与方法

1.1 研究区概况

研究区位于陕西安塞县中国科学院安塞水土保持综合试验站的山地试验场(109°19′E,36°51′N,海拔1068—1309 m),属典型暖温带大陆性季风气候,年均气温8.8 ℃,年均降水量500 mm左右,年际年内分布不均。潜在蒸发在600—800 mm之间变化。土壤以黄绵土为主;主要木本植物有沙棘、油松、刺槐、侧柏(Platycladusorientalis)、榆树(Ulmuspumila)等;草本植物有白羊草(Bothriochloaischaemum)、铁杆蒿(Artemisiagmelinii)、甘草(Glycyrrhizauralensis)等。

1.2 试验材料

以山地试验场野外长期固定监测林地为基础,选取沙棘、油松和刺槐纯林及沙棘*油松、沙棘*刺槐混交林为研究对象。各样地均为2000年左右栽植的人工林,其基本特征见表1。

表1 样地基本特征

混交林1和混交林2分别代表：沙棘*油松和沙棘*刺槐混交林;叶面积指数于生长季内每月月末测量(5—10月),样地调查(胸径/地径、树高等)于生长季前进行(4月);DBH: Diameter at breast height；DGH: Diameter at ground height

1.3 土壤体积含水量

在混交林1和混交林2及油松、沙棘纯林样地内利用ECH2O土壤含水量监测系统长期监测土壤水分动态变化,测定深度为150 cm,分层测量(0—10 cm、10—20 cm、20—50 cm、50—100 cm、100—150 cm)。对于刺槐纯林样地,在干湿季采用土钻法取样烘干法测量土壤水分变化,深度与ECH2O保持一致。

1.4 叶面积指数

于2017年生长季内(5—10月)每月月末,采用LAI-2200 C植物冠层分析仪进行测量。

1.5 光响应曲线

本研究中将同一区域相邻间隔区作为重复样地,即：每个样地均先选取20 m×20 m的样方,然后再划分为20 m×6 m的3个小样方作为样地重复。并在各小样方内按照平均胸径/地径-平均树高-平均冠幅选取标准木代表整个林分的光合特征(共21株)。

2017年7月和9月,选择晴朗天气9:00—11:30,从各标准木中选取中上部成熟向阳叶片,采用Li-Cor 6400便携式光合作用测定系统(Li-Cor, USA)进行光响应曲线测定(每株标准木测定1次光响应曲线)。使用红蓝光源,CO2浓度设定为环境CO2浓度(约400 μmol/mol),光强梯度设定为1500、1200、1000、800、600、400、200、150、100、50、20、0 μmol m-2s-1。

1.6 光响应曲线模型

本研究采用RH[7]、NRH[8]、EM[9]、MRH[10]及MEM[11]模型拟合各树种叶片光响应曲线,公式如下：

(1)

(2)

(3)

(4)

Pn(I)=λe-βI-γe-ξI

(5)

式中,Pn为净光合速率(μmol m-2s-1);I为光强即PAR(μmol m-2s-1);Ic为光补偿点即LCP(μmol m-2s-1);Pnmax为最大净光合速率(μmol m-2s-1);Rd为暗呼吸速率(μmol m-2s-1);α为初始量子效率即Φ(mol/mol);θ为曲角,反映光合曲线弯曲程度的凸度;β(m2s μmol-1)和γ(m2s μmol-1)是独立于I的系数。

1.7 模型评价指标

采用相关系数(R2)和平均相对误差绝对值(MAPE)比较5种模型的拟合优度。MAPE计算公式如下[26]：

(6)

2 数据处理

使用SigmaPlot 10.0(Systat Software Inc., San Jose, CA, USA)进行曲线拟合,采用SPSS 16.0(SPSS Inc., Chicago, IL, USA)非线性回归模块估计各模型参数;在对光响应数据进行正态分布和方差齐性检验的基础上,采用One-Way ANOVA对土壤含水量进行方差分析,使用Paried-Sample T Test对同一树种不同水分条件下及同一水分不同造林模式条件下各树种叶片Pn光响应曲线进行显著性检验,使用一般线性模型(General Linear Model)对不同水分、造林模式条件下各树种光响应曲线拟合参数进行显著性检验。

根据实测值绘制Pn光响应曲线,并依据曲线走势估计Pnmax、Rd、LCP和LSP,作为各树种叶片Pn光响应特征的实测值。

3 结果与分析

3.1 水分条件

2017年生长季内(5月1—10月31日)研究区总降水量为442.8 mm,降水分布不均,其中70.6%的降水集中在7月22日—10月31日;而仅有3.16%的降水发生在6月6日—7月11日(图1)。第一次观测时期(7月12—14日),纯林沙棘、油松、刺槐及混交林1、混交林2样地0—150 cm土层土壤含水量均值分别为：0.073、0.087、0.058、0.087 m3/m3和0.088 m3/m3。与7月相比,第二次观测时期(9月17—19日)各样地土壤含水量均显著增加(P<0.05)分别为：0.112、0.163、0.139、0.119、0.122 m3/m3。

图1 生长季内降水量和土壤体积含水量。Fig.1 Changes in rainfall and soil volumetric moisture content in the studied plots during the growing season 箭头表示光响应曲线观测时段；D,干旱时期Dry;W,湿润时期Wet

3.2 光响应过程

不同种植模式条件下,各树种叶片Pn对光强(PAR)响应规律基本相似,随PAR的增大Pn先迅速增大(PAR < 200 μmol m-2s-1)后缓慢增大至光强临界点即LSP,其对应的Pn为Pnmax,此后Pn随PAR变化因土壤水分条件和树种不同具有较大差别(图2)。9月各树种叶片Pn均显著高于7月(P<0.05),且当PAR > 800 μmol m-2s-1左右时,7月纯林刺槐叶片Pn表现出明显下降趋势,当PAR > 1200 μmol m-2s-1左右时,7月纯林沙棘和混交林2刺槐叶片Pn也有减小趋势。

除9月纯林油松和混交林1油松无显著差异外,纯林和混交林其他各树种叶片Pn均有显著性差异(P<0.05)(表2)。

图2 不同水分条件下各树种叶片净光合速率的光响应曲线Fig.2 Net photosynthetic rate(Pn)-light response curves of each tree species under different water conditions

3.3 光响应过程拟合

利用RH、NRH、MRH、EM和MEM模型对不同水分条件下各树种叶片Pn光响应过程及其特征参数拟合结果(图3,4、表3)。其中,9月各模型MAPE均小于7月(图4)且决定系数(R2> 0.99以上)均大于7月(表3)。此外,通过实测值和模型拟合值对比分析,发现在拟合不同水分条件下各树种叶片光响应曲线时,RH、NRH和EM模型拟合的Pnmax大于实测值,LSP小于实测值,MEM拟合的Rd和LCP与实测值接近,MRH拟合的Pnmax和LSP与实测值接近,整体拟合效果最好(图3,4、表3)。

表2 不同种植模式条件下各树种净光合速率t检验结果

D,干旱时期Dry;W,湿润时期Wet

图3 不同光响应模型对各树种光合速率光响应曲线的模拟Fig.3 Simulation of photosynthetic rate-light response curve of each tree species by five models under different soil water conditions RH: Rectangular hyperbola；NRH: Non-rectangular hyperbola；EM: Exponential model；MRH: Modified rectangular hyperbola；MEM: Modified exponential model

表3 各树种光合生理指标的模型估计值与实测值

RH: Rectangular hyperbola；NRH: Non-rectangular hyperbola；EM: Exponential model；MRH: Modified rectangular hyperbola；MEM: Modified exponential model；D: Dry；W: Wet；Pnmax: Maximum net photosynthetic rate；Rd: Dark respiration rate；LCP: Light compensation point；LSP: Light saturation point

图4 5种模型对不同水分条件下各树种光响应曲线拟合的稳定性分析Fig.4 Stability analysis of the five models for the light responses curve of each tree species under different water conditions

3.4 水分条件和造林模式对光合参数的影响

依据5种模型对光合参数求解精度,对各树种叶片Rd和LCP分析采用MEM模型,对Pnmax和LSP分析采用MRH模型(表3),且因MRH整体拟合效果最好,故应用该模型求解的Φ、β和γ值进行表观量子效率特性及光饱和、光抑制现象分析(表5)。对不同水分、造林模式条件下各个光合参数方差分析结果(表6)。

7月各树种叶片Rd均较9月降低,且同一水分条件下,沙棘、刺槐叶片Rd整体低于油松(表3)。此外,沙棘与油松混交后,叶片Rd受水分条件、造林模式及二者交互作用影响(P<0.01),与刺槐混交后,受水分条件及二者交互作用影响(P<0.05);水分条件和造林模式对油松Rd有显著影响(P<0.05);刺槐仅受二者交互作用影响(P<0.05)(表6)。7月各树种叶片LSP均显著低于9月(P<0.05),而LCP均高于9月(表3)。油松LSP受水分条件及二者交互作用影响(P<0.01);沙棘和刺槐受造林模式、水分条件及二者交互作用影响(P<0.01)。油松LCP主要受水分条件影响(P<0.01);沙棘与油松、刺槐混交后均受造林模式和水分条件的影响(P<0.01);刺槐受造林模式、水分条件及二者交互作用影响(P<0.01)(表6)。7月各树种叶片Pnmax均显著低于9月(P<0.05),且油松、刺槐Pnmax受造林模式、水分条件及二者交互作用影响(P<0.01);沙棘与油松混交后,受造林模式和水分条件的影响(P<0.01),而与刺槐混交后,仅受水分条件影响(P<0.01)(表3,6)。

9月各树种叶片Φ均高于7月,且同一水分条件下,沙棘和刺槐叶片Φ整体大于油松(表5)。油松和沙棘Φ仅受土壤水分条件的影响(P<0.01);刺槐受造林模式和水分条件的影响(P<0.05)(表6)。此外,9月各树种叶片β和γ值均小于7月,且7月沙棘和刺槐叶片β和γ值均显著大于油松(P<0.05)。油松、刺槐β仅受水分条件影响(P<0.05);沙棘β受造林模式、水分条件及二者交互作用影响(P<0.01);沙棘γ受水分条件影响(P<0.01),刺槐受水分条件及二者交互作用影响(P<0.05)(表5,6)。

表5 MRH模型对各树种光响应曲线拟合的表观量子效率(Φ)、光抑制项(β)和光饱和项(γ)

Table 5 Estimation parameters(apparent quantum efficiency, the item of photo-inhibition (β) and light saturation (γ))of the light response curve of each tree species by the MRH model

样地Sample plot水分条件Water condition光合参数Photosynthetic parameterΦ(mol mol-1)β/(m2 s mol-1)γ/(m2 s mol-1)纯林沙棘D1.01×10-22.35×10-47.22×10-3Pure plantation of H. rhamnoidesW4.61×10-21.74×10-41.17×10-3混交林1 沙棘D1.43×10-28.64×10-46.17×10-3Mixed plantation 1 of H. rhamnoidesW5.33×10-22.29×10-42.52×10-3混交林2 沙棘D1.75×10-28.37×10-45.16×10-3Mixed plantation 2 of H. rhamnoidesW5.54×10-22.54×10-41.45×10-3纯林油松D4.19×10-34.18×10-53.66×10-3Pure plantation of P. tabuliformisW4.33×10-21.90×10-52.21×10-3混交林1 油松D6.14×10-33.54×10-53.84×10-3Mixed plantation 1 of P. tabuliformisW3.49×10-22.19×10-52.75×10-3纯林刺槐D4.23×10-23.40×10-46.58×10-3Pure plantation of R. pseudoacaciaW5.75×10-21.53×10-41.84×10-3混交林2 刺槐D3.63×10-24.67×10-48.55×10-3Mixed plantation 2 of R. pseudoacaciaW5.62×10-22.60×10-41.55×10-3

4 讨论

4.1 光响应模型适用性评价

光响应曲线模拟对了解光合作用光反应过程非常重要,被广泛应用于植物光合生理的研究中[9,15,26]。近年来,诸多学者已提出了不同的光响应模型,从经验模型到机理模型,从简单模型到复杂模型[12, 16,27]。不同模型在具体地应用中往往存在较大差异,本研究表明,不同水分条件下,RH和NRH模型拟合的各树种叶片Pnmax均大于实测值(表3),这与朱永宁等[28]、郎莹等[15]、陈志成等[29]以及周多多等[16]对不同水分条件下玉米(ZeamaysL.)、山杏、栾树(Koelreuteriapaniculata)及胡杨的研究结果类似,这可能主要是因为这两个模型曲线均为渐近线,没有极值,需要利用非线性最小二乘法对Pnmax进行估算,如果光响应曲线在饱和光强之后下降,估算值必然大于实测值[30]。此外,RH、NRH、EM和MEM模型因不能直接求解LSP,估算结果通常小于实测值,本研究这4个模型对不同水分条件下各树种的拟合效果也证实了这一结论(表3)。MRH模型函数因存在极值,能直接求解Pnmax和LSP,并能较好地拟合光抑制阶段的光响应过程,故逐渐被应用[15, 31],本研究也表明干旱时期MRH对纯林沙棘、刺槐以及混交林2刺槐叶片的光抑制现象拟合效果较好。

表6 各树种光响应曲线的拟合参数方差分析结果

此外,5种模型对不同水分条件下各树种叶片Pn光响应曲线的拟合优度均较高(R2> 0.71),其中9月R2均在0.99以上,说明对湿润时期拟合效果整体较好。然而,尽管R2的大小可以作为判断光响应模型的重要指标,却不能说明所得到的光合参数与实际情况相吻合[32]。因此,与R2相结合,本研究采用平均相对误差绝对值(MAPE)比较5种模型拟合效果,确定出不同水分条件下各树种叶片Pnmax和LSP优先选择MRH模型,Rd和LCP优先选择MEM模型。且5种模型拟合效果优劣顺序：MRH整体拟合效果最佳,其次是MEM、NRH和EM,而RH模型拟合效果最差(表3、图4),这与藏川杨(Populusszechuanica)[33]及长白落叶松[12]对不同模型适应性研究结果一致。

4.2 水分条件、造林模式与光合参数关系

光合作用是制约人工林生产力最重要的生理过程,影响植物光合作用的环境因子主要有水分、CO2、光照、温度等。有研究表明,混交造林后,可能由于混交林和纯林在自身物种组成上的差异,造成林分内部生态环境的不同,从而可能导致林分各组成树种光合生理特性的差异[34]。

Φ反映植物在弱光下吸收、转换和利用光能的能力[35]。通常,在适宜生长条件下,植物Φ在0.03—0.05之间[36]。本研究表明,9月各树种叶片Φ均高于7月,表明湿润时期,3个树种对弱光利用能力增强(表5)。在湿润条件下,各树种Φ在0.035—0.058之间,与王荣荣等[37]对杠柳(Periplocasepium)研究结果相近(0.042—0.066),位于一般植物Φ范围。而檧木(Araliaelata)Φ在0.004—0.029之间,山杏Φ在0.015—0.083之间[6,38],可见本研究3个树种在适宜水分条件下利用弱光的能力高于一般灌木或小乔木。其大小顺序为：混交林1油松(0.035)<纯林油松(0.043)<纯林沙棘(0.046)<混交林1沙棘(0.053)<混交林2沙棘(0.055)<混交林2刺槐(0.056)<纯林刺槐(0.0575)(表5),与纯林油松、刺槐相比,混交林1油松和混交林2刺槐叶片Φ分别下降19.40%和20.87%,而与纯林沙棘相比,混交林1沙棘和混交林2沙棘分别增加15.62%和20.17%,说明混交造林有助于沙棘提高对弱光的利用效率;刺槐和沙棘Φ整体大于油松,表明刺槐和沙棘在把光能转化为净能量的能力强于油松,即阔叶树种Φ大于针叶树种,这与其他研究结果相类似：在植物生长旺盛期测定杉木(Cunninghamialanceolata)的Φ小于栓皮栎(Quercusvariabilis)和枫香树(Liquidambarformosana)[39-40],也小于华东椴(Tiliajaponica)[41]。这可能主要和植物自身生物学特性及其适应环境的能力的差异有关[42]。

LSP和LCP反映植物对光照的需求,是判断植物适应高光强能力和耐荫性强弱的重要指标[43]。一般认为阴性植物LCP低于20 μmol m-2s-1,LSP低于500 μmol m-2s-1;而阳性植物LCP在50 μmol m-2s-1以上,LSP在800 μmol m-2s-1以上[43-44]。在本研究中,干旱时期各树种LCP在52.86—251.03 μmol m-2s-1之间,LSP在419.33—969.59 μmol m-2s-1之间;湿润时期各树种LCP在11.85—90.96 μmol m-2s-1之间,LSP在1030.02—1939.07 μmol m-2s-1之间;可见随水分环境变化,各树种对光照强度表现出一定的适应性和可塑性,即在适宜水分条件下对弱光与强光的利用能力都较高,因此具有较宽的光照生态幅,耐荫性和喜光性都增强,较适合在干旱、半干旱地区生长。其生态幅宽窄顺序为：混交林2刺槐(1895.81 μmol m-2s-1)>混交林2沙棘(1673.23 μmol m-2s-1)>纯林油松(1496.24 μmol m-2s-1)>混交林1沙棘(1466.6 μmol m-2s-1)>纯林刺槐(1302.27 μmol m-2s-1)>混交林1油松(1172.67 μmol m-2s-1)>纯林沙棘(977.23 μmol m-2s-1),与对应纯林相比,刺槐与沙棘混交后,刺槐叶片生态幅增加45.58%,沙棘增加71.22%;沙棘与油松混交后,油松叶片生态幅下降21.62%,而沙棘增加50.08%,说明水分适宜条件下,混交造林有利于刺槐和沙棘增加其对光强的适应范围。

此外,MRH模型中参数光抑制项(β)和光饱和项(γ)的生物学意义是光系统Ⅱ天线色素分子光量子吸收截面与其处于激发平均寿命的乘积,β和γ越大,表示植物越容易受到光抑制和光饱和[45]。本研究表明,与7月相比,9月各树种叶片β和γ均减小,表明湿润条件下各树种均不易出现光抑制和光饱和现象(图3、表5)。7月沙棘和刺槐β和γ均显著大于油松(P<0.05),表明干旱时期,相对于油松,沙棘和刺槐可能更易出现光抑制和光饱和现象(图2)。这与王荣荣等[37]对杠柳叶片β和γ随干旱胁迫变化研究结果类似,可能主要和针阔叶树种叶片结构特征有关。有研究表明,在针叶内部可能有光通量密度梯度的存在,这样就阻止了所有的光合单元达到饱和。针叶的厚度是600—1200 μm,仅仅传输1%—5%的可见光。而200—600 μm厚的杨树叶传输约5%—20%的可见光,比针叶大4—5倍[46-47]。

植物Rd越大说明其呼吸速率越大。本研究中,7月各树种叶片Rd值均小于9月,表明在干旱条件下,3个树种可能通过减弱呼吸作用以减少对光合产物的消耗,以适应土壤干旱条件;干旱时期各树种Rd大小顺序为：纯林油松(3.038 μmol m-2s-1)>混交林1油松(2.608 μmol m-2s-1)>纯林沙棘(1.318 μmol m-2s-1)>混交林1沙棘(1.294 μmol m-2s-1)>混交林2沙棘(1.127 μmol m-2s-1)>混交林2刺槐(1.082 μmol m-2s-1)>纯林刺槐(0.462 μmol m-2s-1),与对应纯林相比,油松与沙棘混交后,油松叶片Rd下降14.15%,沙棘下降1.82%;刺槐与沙棘混交后,刺槐叶片Rd增加134.2%,沙棘下降14.49%,说明在干旱时期,油松和沙棘混交有利于光合产物的积累;不同水分条件下,各树种Rd在0.46—3.92 μmol m-2s-1之间(表3),与檧木Rd(0.24—0.86 μmol m-2s-1)[6],山杏Rd(0.67—1.92 μmol m-2s-1)[38]相比,各树种叶片Rd偏高,这可能主要和树种在不同水分条件下对弱光环境的光合生理适应差异性有关[37]。

Pnmax是衡量植物叶片光合潜力的指标,在最佳的环境条件下,可表示其最大光合能力[48]。本研究表明,9月各树种叶片Pnmax均显著增高(P<0.05),说明湿润时期各树种叶片光合潜力增强。其大小顺序为：混交林1沙棘(18.721 μmol m-2s-1)>混交林2沙棘(18.1 μmol m-2s-1)>纯林油松(17.886 μmol m-2s-1)>纯林沙棘(16.931 μmol m-2s-1)>纯林刺槐(14.851 μmol m-2s-1)>混交林2刺槐(10.778 μmol m-2s-1)>混交林1油松(10.38 μmol m-2s-1),与对应各纯林相比,发现油松与沙棘混交后,油松Pnmax下降41.97%,沙棘增加10.57%;刺槐与沙棘混交后,刺槐Pnmax下降27.43%,沙棘增加6.46%,说明在适宜水分条件下,混交造林有助于沙棘提高其光合潜力,混交林1油松叶片光合潜力最低,这可能主要和不同树种叶片排列结构有关。有研究表明,针叶树种叶片密集排列导致的相互遮荫可能会限制针叶叶束的光合作用。如一个云杉叶束的投影面积仅是单个针叶投影面积之和的一半左右[46]。

黄土丘陵区属于季节性干旱事件多发的地区,研究不同种植模式下主要造林树种在不同水分条件下光合生理特性差异,可为该区域不同尺度植被生理生态过程模拟提供依据。本试验着重分析了生长季内典型时期各树种光响应过程及其模拟,因此只选取了生长季相对干旱和湿润两个有代表性的时段进行分析。而后续的分析还需基于更多水分梯度的观测或者结合降水事件,才能深入理解不同树种对区域水分环境变化的响应规律。

5 结论

(1)纯林和混交林沙棘、油松和刺槐在湿润时期叶片Pn均显著高于干旱时期(P<0.05)。

(2) 5种模型对湿润时期纯林和混交林各树种叶片光响应过程整体拟合效果较好(R2> 0.99),不同模型拟合效果的优劣顺序为：MRH模型对各树种叶片光响应过程整体拟合效果最佳,其次是MEM、NRH和EM模型,而RH模型拟合效果最差。

(3)与干旱时期相比,湿润时期各树种叶片均具有较高的光合生理参数值,表明适宜水分条件下,其光合潜力、呼吸活性以及对光能的利用能力较强。

(4)与纯林相比,混交林沙棘光合生理参数值整体呈上升趋势,表明混交造林有助于提高其光合潜力、对光能的利用能力及光照生态幅宽度;油松整体呈下降趋势,表明混交种植可能会降低其光合作用能力;刺槐变化无明显规律,其叶片Φ和Pnmax值下降,Rd增加,光照生态幅最宽,表明混交种植可能会增强其对强光的适应性,但会降低其对弱光的利用能力及光合潜力。