毛烏素沙地油蒿光系統(tǒng)II多時(shí)間尺度的環(huán)境響應(yīng)特征

靳 川，蔣 燕，李鑫豪，徐銘澤，高圣杰，魏寧寧，賈 昕，田 赟，查天山

（1. 北京林業(yè)大學(xué)水土保持學(xué)院寧夏鹽池毛烏素沙地生態(tài)系統(tǒng)國家定位觀測研究站，北京 100083；2. 北京林業(yè)大學(xué)水土保持國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，北京 100083）

0 引言

作為維持生態(tài)系統(tǒng)物質(zhì)代謝和能量循環(huán)的基礎(chǔ)，植物光合作用非常容易受到環(huán)境波動(dòng)的影響，不同強(qiáng)度和持續(xù)時(shí)間的環(huán)境脅迫會(huì)協(xié)同抑制光合過程[1-2]。光系統(tǒng)II（Photosystem II，PSII）能量分配被認(rèn)為是光合過程中最為敏感的部分，PSII天線色素分子吸收光能有3個(gè)去向：1）進(jìn)行光化學(xué)反應(yīng)，通過碳同化合成有機(jī)物[3]；2）通過葉黃素循環(huán)以熱能形式將過剩光能耗散[4-5]；3）重新釋放一個(gè)光子，即產(chǎn)生葉綠素?zé)晒鈁6-7]。3種途徑相互競爭權(quán)衡，且受光合有效輻射（Photosynthetically Active Radiation，PAR）、空氣溫度（Ta）、飽和水汽壓差（Vapor Pressure Difference，VPD）和土壤含水率（Soil Water Content，SWC）等環(huán)境因素影響[5,8-9]。上述環(huán)境變量通過各自波動(dòng)周期、頻率和振幅調(diào)控植物PSII能量分配[10-11]。在分-秒時(shí)間尺度上，PAR增加激活葉黃素循環(huán)，植物通過改變調(diào)節(jié)性熱耗散（ΦNPQ）和非調(diào)節(jié)性熱耗散（ΦNO）的比例耗散過剩光能[12]；在小時(shí)時(shí)間尺度上，VPD升高可導(dǎo)致蒸騰加速，植物調(diào)節(jié)氣孔減少水分損失[13]，同時(shí)也不可避免地抑制了實(shí)際光化學(xué)效率（ΦPSII）；在日尺度上，由PAR、Ta和VPD的天循環(huán)驅(qū)動(dòng)PSII能量分配[14]；在幾天至數(shù)月時(shí)間周期上，高輻射、寒潮、熱浪和降水等天氣事件以及植物的物候期都會(huì)導(dǎo)致PSII能量分配變化[10,15]。盡管有上述認(rèn)識，但目前關(guān)于PSII能量分配的相關(guān)研究大多集中在室內(nèi)短期控制試驗(yàn)，難以模擬自然中真實(shí)的環(huán)境動(dòng)態(tài)[16-17]。具有高太陽輻射、極端溫度和水文梯度跨度大等特點(diǎn)的荒漠生態(tài)系統(tǒng)則更加復(fù)雜[18]?；哪募竟庹者h(yuǎn)超植物所需，導(dǎo)致光合機(jī)構(gòu)產(chǎn)生氧化損傷[19]；極端溫度抑制了相關(guān)酶的活性，造成光合效率降低[20-21]；空氣干濕循環(huán)引起植物氣孔調(diào)節(jié)[13,22]，導(dǎo)致熒光淬滅[23]。另外，荒漠植物也常因長期缺水造成PSII失活，產(chǎn)生可逆的光損傷[20-21,24]。因此，荒漠植物PSII能量分配與各環(huán)境因子在多時(shí)間尺度上很可能具有明顯關(guān)聯(lián)特征。

頻譜處理可檢驗(yàn)數(shù)據(jù)序列的時(shí)間動(dòng)態(tài)格局[25]，如小波分析。小波分析方法已經(jīng)廣泛應(yīng)用于多個(gè)生態(tài)相關(guān)領(lǐng)域，適用于探究多組時(shí)間序列數(shù)據(jù)之間關(guān)系[25-27]。然而，目前幾乎沒有研究將小波分析應(yīng)用于葉綠素?zé)晒獾拈L期監(jiān)測數(shù)據(jù)上，尤其缺乏植物對環(huán)境波動(dòng)在多時(shí)間尺度下響應(yīng)過程和調(diào)節(jié)機(jī)制的理解認(rèn)識。因此，本研究對寧夏鹽池毛烏素沙地優(yōu)勢荒漠物種油蒿（Artemisia ordosica）進(jìn)行葉綠素?zé)晒忾L期原位連續(xù)監(jiān)測，采用小波分析方法探究油蒿PSII能量分配在多時(shí)間尺度下的環(huán)境響應(yīng)特征，旨在為荒漠植物對不同環(huán)境脅迫的敏感程度、適應(yīng)能力和生理可塑性提供數(shù)據(jù)支撐，同時(shí)為荒漠地區(qū)植被恢復(fù)物種選擇提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于寧夏吳忠市鹽池荒漠生態(tài)系統(tǒng)定位研究站（37°53"08""N，107°25"46"" E，平均海拔1530 m)，是半干旱區(qū)和干旱區(qū)之間的交錯(cuò)帶，由于過度放牧等人為因素該區(qū)植被已嚴(yán)重退化，近20年的封育使該區(qū)植被得到了有效的恢復(fù)，這也被認(rèn)為是該區(qū)荒漠化逆轉(zhuǎn)的重要標(biāo)志。該區(qū)屬溫帶大陸性季風(fēng)氣候，平均氣溫為8.3 ℃，晝夜溫差大，全年晝夜平均溫差可達(dá)25 ℃以上。年平均潛在蒸散量為2024 mm，年平均降水量為287 mm，降雨稀少且不均勻，主要集中在6—9月，4—5月容易發(fā)生春季干旱。土壤以風(fēng)沙土和黑壚土為主，密度為1.54 g/cm3。盛行西北風(fēng)，年平均風(fēng)速2.6 m/s。研究地植被群落以旱生植物為主，優(yōu)勢灌木物種有油蒿（Artemisia ordosica）、沙柳（Salix psammophila）和楊柴（Hedysarum mongolicum），觀測期內(nèi)相對覆蓋度分別為53 %、17 %和12 %。

1.2 試驗(yàn)材料

油蒿廣泛分布于中國北部及西北部荒漠地區(qū)，因其根系發(fā)達(dá)耐沙埋、耐旱，在固定流沙和促進(jìn)群落演替中發(fā)揮著重要作用[17]。油蒿通常在4月萌發(fā)，8月積累生物量達(dá)到最大，9月末開始凋落。樣地內(nèi)油蒿平均密度為1.43株/m2，平均高度為34.7 cm，平均冠幅為0.64 m。選取3株長勢相同且具有代表性（高度為50 cm，冠幅為0.8 m，約50分枝左右）的油蒿作為樣本植株。避免監(jiān)測期內(nèi)人為因素的干擾，保證樣本植株處于完全自然生長的狀態(tài)，并進(jìn)行物候的同步觀測。

1.3 葉綠素?zé)晒獗O(jiān)測與計(jì)算

于2016年4月28日至10月3日，采用多通道連續(xù)監(jiān)測熒光儀（MONITORING-PAM2000，WALZ公司，德國）對油蒿樣本植株進(jìn)行24 h原位連續(xù)監(jiān)測，測量間隔為30 min。用熒光探頭的特制葉夾固定葉片，調(diào)整探頭方位使葉片充分受光，每隔7 d定期檢查葉片并及時(shí)更換。測定葉綠素?zé)晒鈪?shù)包括飽和脈沖下熒光（Fm"）和實(shí)時(shí)熒光（Fs）。將夜間的Fm"、Fs分別作為暗處理最大熒光（Fm）和暗處理最小熒光（Fo）。根據(jù)Kramer公式計(jì)算PSII能量分配[28]，實(shí)際光化學(xué)效率（ΦPSII，%）、調(diào)節(jié)性熱耗散（ΦNPQ，%）和非調(diào)節(jié)性熱耗散（ΦNO，%）的計(jì)算方法如下：

式中Fm"為飽和脈沖下熒光，F(xiàn)s為實(shí)時(shí)熒光，F(xiàn)m為暗處理最大熒光。

采用式（4）計(jì)算最大光化學(xué)效率（Fv/Fm）

式中Fo為暗處理最小熒光。

1.4 環(huán)境氣象因子監(jiān)測與計(jì)算

光合有效輻射由距熒光儀100 m處通量塔（高度為6 m）配有的輻射傳感器（PAR-LITE，Kipp & Zonen公司，荷蘭）測得?？諝鉁囟群拖鄬穸扔伤蠝貪穸葌鞲衅鳎℉MP155A，Vantaa公司，芬蘭）測得，采用式（5）和式（6）[16]計(jì)算飽和水汽壓差（VPD，kPa）

式中e為飽和空氣壓力，kPa；RH為相對濕度，%；Ta為空氣溫度，℃。

土壤含水率由樣地內(nèi)土壤溫濕度傳感器（ECH2O-5TE，Decagon Devices公司，美國）測得，深度為30 cm。降雨量由距樣地200 m處空地內(nèi)雨量筒（TE525W，Campbell公司，美國）測得。所有環(huán)境氣象因子測量間隔均為30 min。

1.5 數(shù)據(jù)處理與分析

本研究使用連續(xù)小波變換（Constant Wavelet Transform，CWT）和小波相干性（Wavelet Coherence，WTC）分析熒光參數(shù)的時(shí)間序列。CWT是將時(shí)間序列數(shù)據(jù)在固定間隔的變量進(jìn)行平移或伸縮等變換，用母小波來表示。CWT可以反映時(shí)間序列數(shù)據(jù)在不同時(shí)間尺度的波動(dòng)[29]，數(shù)量為N的離散變量xn的積分譜 ()xnWs定義為

式中δt為時(shí)間間隔，為復(fù)共軛尺度，s為小波變換尺度，n為數(shù)組1至N，n"為正整數(shù)取值。本研究選取Morlet小波作為母小波，采用式（8）計(jì)算xn的小波功率譜（Sn）[26,30]，得到其貢獻(xiàn)最大的時(shí)間周期

定義2個(gè)時(shí)間序列數(shù)據(jù)之間的交叉小波功率譜（Cn），相位角譜（An）和WTC功率譜（Rn2）[26,30]

式中x為環(huán)境因子，y為植物熒光參數(shù)，為變量x與y的交叉小波變換，S為時(shí)間尺度的平滑處理，和分別為的實(shí)部和虛部[31]（在時(shí)頻域圖中，An可以用箭頭來繪制，若x與y正相關(guān)，向右的箭頭表示x與y的相位是一致的，向左則相反。而向下的箭頭表示y滯后于x90°或270°，向上的箭頭表示x滯后于y90°或270°。每360°為1個(gè)時(shí)間周期）。

2 結(jié)果與分析

2.1 環(huán)境因子與油蒿PSII能量分配參數(shù)動(dòng)態(tài)

PAR求日最大值和平均值，Ta求日最大值、日最小值和日平均值，VPD和SWC求日平均值，降雨量求日總值。如圖1所示，觀測期內(nèi)PAR和Ta趨勢相近，在5 —7月呈上升趨勢，8—9月呈下降趨勢。PAR最高達(dá)1963μmol/(m2·s)（6月8日），7月和8月出現(xiàn)持續(xù)多天低值（PAR＜400μmol/(m2·s)）。Ta最高值達(dá)32.1 ℃（7月29日），最低值為2.9 ℃（5月6日），出現(xiàn)1次大幅降溫事件（8月22日）。VPD在5—7月呈周期性波動(dòng)，8—9月較低，日平均值最高達(dá)2.5 kPa（7月25日）。SWC整體呈持續(xù)下降趨勢，在8月20日降雨量達(dá)到52.6 mm的降水后有較大增加。如圖2所示，ΦPSII在7月和8月有2處明顯的下降，日平均值最低為0.45%（7月16日）。ΦNPQ和ΦNO趨勢大致相同，7月和8月有2處明顯增高。ΦNPQ日平均值最高達(dá)23 %（7月16日），ΦNO日平均值最高達(dá)36%（7月3日）。Fv/Fm呈周期性波動(dòng)，最高0.82（5月31日），最低0.68（5月3日）。觀測期內(nèi)，各葉綠素?zé)晒鈪?shù)和環(huán)境因子具有明顯季節(jié)波動(dòng)。

PAR、Ta和VPD都呈現(xiàn)先升高后降低的日變化趨勢，PAR峰值在13:00左右，Ta和VPD峰值約在16:00。PAR在5—7月較高，Ta在5月和9月較低，VPD在6—7月較高。ΦPSII呈現(xiàn)與ΦNPQ和ΦNO相反的日變化趨勢，峰值出現(xiàn)時(shí)間與PAR一致。ΦPSII在6—7月較低，ΦNPQ在5 —7月的峰值較高，ΦNO在7月的峰值較高（圖3）。

2.2 環(huán)境因子與油蒿PSII能量分配參數(shù)在各時(shí)間周期內(nèi)變化

除SWC以外，所有變量全局小波功率譜在1 d時(shí)間周期均出現(xiàn)強(qiáng)峰值，并且在4～40 d時(shí)間周期內(nèi)分別出現(xiàn)多次振蕩。葉綠素?zé)晒鈪?shù)（ΦPSII、ΦNPQ和ΦNO）峰值相對環(huán)境因子（PAR、Ta、VPD和SWC）較小。SWC既無日尺度譜峰，也無季節(jié)尺度譜峰。變量PAR、ΦPSII、ΦNPQ和ΦNO在小于1 d的時(shí)間周期內(nèi)出現(xiàn)多次峰值，但其功率遠(yuǎn)低于1 d時(shí)間周期處（圖4）。

小波功率譜可以清楚地顯示各變量在觀測期內(nèi)的時(shí)頻域局部特征。在0.5 d時(shí)間周期處，PAR出現(xiàn)顯著小波功率區(qū)域。在1 d時(shí)間周期處，PAR、Ta和VPD均出現(xiàn)顯著小波功率區(qū)域，區(qū)域非連續(xù)且強(qiáng)烈震蕩。Ta在4～32 d周期內(nèi)、VPD在8～32 d周期以及SWC在1～32 d周期內(nèi)均出現(xiàn)小波功率值的升高（圖5）。在1 d時(shí)間周期處，ΦPSII、ΦNPQ和ΦNO均出現(xiàn)強(qiáng)烈震蕩的顯著區(qū)域。ΦPSII在8～32 d周期內(nèi)、ΦNPQ在4～32 d周期以及ΦNO在16～32 d周期內(nèi)也呈現(xiàn)小波功率值的升高（圖6）。

2.3 油蒿PSII能量分配參數(shù)與環(huán)境因子的小波相關(guān)性

油蒿PSII能量分配參數(shù)與環(huán)境因子在多個(gè)時(shí)間尺度上存在顯著小波相關(guān)區(qū)域，分布不均且波動(dòng)明顯。在日尺度上，ΦPSII與PAR、Ta、VPD和SWC小波相關(guān)區(qū)域相位角分 別 為(12.32±3.56)°、(47.88±21.20)°、(73.25±13.01)°和(183.69±96.81) °，即ΦPSII滯后于PAR、Ta、VPD和SWC分別為(0.82±0.24)、(3.19±1.41)、(4.88±0.87)和(12.25±6.45) h（圖7；表1）。ΦNPQ滯后于PAR、Ta、VPD和SWC分別為(0.62±0.27)、(3.45±1.57)、(4.94±2.64)和(11.87±4.98) h（圖 8；表1）。ΦNO滯后于PAR、Ta、VPD和SWC分別為(0.72±0.32)、(3.27±1.43)、(3.83±1.69)和(8.07±6.97) h（圖 9；表 1）。油蒿PSII能量分配參數(shù)（ΦPSII、ΦNPQ和ΦNO）滯后于PAR、Ta、VPD和SWC分別平均為43 min、3.3、4.6和10.72 h。

在季節(jié)尺度上，ΦPSII滯后于PAR、Ta、VPD和SWC分 別 為(13.64±5.54)、(9.87±3.47)、(7.68±1.99)和(23.47±12.46) d（圖7；表2）。ΦNPQ滯后于PAR、Ta、VPD和SWC分別為(16.05±7.65)、(3.52±2.73)、(9.90±5.64)和(17.89±14.80) d（圖8；表2）。ΦNO滯后于PAR、Ta、VPD和SWC分別為(14.44±6.12)、(8.23±4.47)、(8.73±5.05)和(23.73±10.97) d（圖9；表2）。油蒿PSII能量分配參數(shù)分別滯后于Ta、VPD、PAR和SWC 7.2、8.8、14.7和21.6 d。

表1 PSII能量分配參數(shù)在日尺度上滯后環(huán)境因子的時(shí)間 Table 1 Lag time of PSII energy partitioning parameters relative to environmental factors on daily scale h

表2 PSII能量分配參數(shù)在季節(jié)尺度上滯后環(huán)境因子的時(shí)間 Table 2 Lag time of PSII energy partitioning parameters relative to environmental factors on seasonal scale d

3 討論

植物通過熱耗散機(jī)制消耗光合過程中的過剩光能，減少活性氧以避免產(chǎn)生光損傷[23]。將熱耗散分為ΦNPQ和ΦNO可獲得更精準(zhǔn)的定量信息，ΦNPQ表征熱耗散等光保護(hù)調(diào)節(jié)機(jī)制的強(qiáng)弱[28]，ΦNO則指示光抑制或光損傷的程度[10,28]。日尺度上，油蒿PSII能量分配參數(shù)與各環(huán)境因子均有很強(qiáng)相關(guān)性（圖7～圖9），對PAR變化最敏感且受其主導(dǎo)。PAR增加導(dǎo)致更多過剩光能需要耗散，因此ΦNPQ和ΦNO與PAR正相關(guān)，同時(shí)不可避免抑制了ΦPSII，這一系列過程是瞬間發(fā)生的[5,16]。油蒿ΦPSII在12:00左右發(fā)生明顯降低，在夜晚恢復(fù)（圖3）。表明在高輻射下，雖然油蒿光化學(xué)反應(yīng)過程受到了限制，但其PSII反應(yīng)中心實(shí)際并沒有產(chǎn)生損害[32-33]。這與同地區(qū)不同物種相關(guān)研究結(jié)果不一致，研究發(fā)現(xiàn)沙柳和楊柴因高輻射表現(xiàn)出嚴(yán)重光抑制，并造成永久性光損傷。這可能與植物形態(tài)學(xué)有關(guān)，沙柳葉呈葉線形，楊柴為互生扁平小葉，而具有針形葉的油蒿則有利于避免過量太陽輻射。油蒿PSII能量分配分別滯后Ta和VPD 3.3和4.6 h（表1）。研究表明Ta降低抑制葉綠體捕光蛋白酶活性，從而影響植物光合生理狀態(tài)[20-21]。清晨ΦPSII下降（圖3），此時(shí)Ta處于全天最低，油蒿PSII光反應(yīng)相關(guān)酶活性受到抑制，導(dǎo)致ΦPSII早于PAR開始降低[21,34]。VPD增加導(dǎo)致氣孔調(diào)節(jié)，蒸騰速率降低[13,22]，抑制了水的裂解[24]，也提高了跨類囊體質(zhì)子梯度[35]，最終造成ΦPSII降低，ΦNPQ和ΦNO增加。油蒿PSII能量分配在1 d時(shí)間周期內(nèi)滯后SWC約10.72 h，遠(yuǎn)高于其他環(huán)境因子（表1）。這說明短時(shí)間內(nèi)油蒿對于土壤水的敏感程度較低[10,34]。

荒漠植物需要應(yīng)對干旱、高輻射和極端溫度等多種環(huán)境脅迫[5,18,36]，PSII長期適應(yīng)機(jī)制的形成需要幾天甚至幾周時(shí)間[10,17]，主要通過調(diào)整光合色素濃度、葉綠素a/b比例或PSII特定蛋白質(zhì)含量等途徑實(shí)現(xiàn)[10,37]。Fv/Fm是表征植物PSII受脅迫程度的關(guān)鍵指標(biāo)，通常情況下在0.75～0.85之間，發(fā)生下降則表示PSII 受損失活[34]。油蒿Fv/Fm和PSII能量分配季節(jié)波動(dòng)明顯，但始終維持在一個(gè)整體水平上（Fv/Fm約為0.78）（圖2）。研究認(rèn)為油蒿可以通過自身修復(fù)調(diào)節(jié)，使其PSII恢復(fù)到正常狀態(tài)[16,28]。研究發(fā)現(xiàn)，油蒿PSII能量分配受物候影響，5月初處于萌芽期，葉綠素含量及活性低，對外界抵抗力較差，易產(chǎn)生光抑制[10,23]。7—8月受高溫高輻射影響，產(chǎn)生光抑制造成PSII損傷，此時(shí)油蒿處于長勢最旺盛的花果期，其生物量和葉片功能性狀達(dá)到最佳，PSII處于不斷受損和修復(fù)的過程[32-34]。9月末油蒿葉片老化進(jìn)入枯萎期，葉綠素含量、各種酶活性伴隨光合單元機(jī)構(gòu)衰退降低，低溫脅迫導(dǎo)致Fv/Fm低值出現(xiàn)[10,38]。

在季節(jié)尺度上，油蒿PSII能量分配滯后Ta和VPD相對PAR和SWC較短。說明在長時(shí)間尺度，油蒿PSII對Ta和VPD更為敏感。植物在高VPD下會(huì)調(diào)節(jié)氣孔關(guān)閉[13,22]，增強(qiáng)葉黃素循環(huán)等持續(xù)性熱耗散[23]，同時(shí)對PSII反應(yīng)中心造成可逆損害[32-34]。極端溫度或高太陽輻射則會(huì)打破PSII能量輸入和利用的平衡，造成PSII捕光蛋白損傷[39]。與其他研究比較，油蒿對極端溫度的適應(yīng)能力要高于荒漠物種沙柳和紅砂（Reaumuria soongorica）[8,36]?；哪貐^(qū)晝夜溫差大，油蒿可通過提高細(xì)胞內(nèi)滲透壓，增加原生質(zhì)抗凝結(jié)性，減緩代謝速率來適應(yīng)極端溫度[39]。油蒿PSII能量分配滯后PAR約14.7 d（表 2），油蒿對于高輻射具有一定的耐受能力[16,34]。干旱是影響植物生長生存的重要環(huán)境脅迫，水分限制著細(xì)胞內(nèi)一系列的生理生化反應(yīng)[13,24]。油蒿具有粗長發(fā)達(dá)的主根和茂密的側(cè)根[17,40]，PSII能量分配滯后SWC時(shí)間最長，約21.6 d（表2），可見其對干旱有著極強(qiáng)的耐受性，因此也經(jīng)常用作干旱區(qū)植被恢復(fù)的先鋒物種。

4 結(jié)論

1）日尺度上，油蒿可能通過葉黃素循環(huán)等熱耗散機(jī)制和氣孔導(dǎo)度調(diào)節(jié)PSII能量分配，主要受光合有效輻射（Photosynthetically Active Radiation，PAR）控制，同時(shí)受空氣溫度（Ta）和飽和水汽壓差（Vapor Pressure Difference，VPD）的調(diào)控，對土壤含水率（Soil Water Content ，SWC）敏感性較差。PSII能量分配參數(shù)滯后于PAR 43 min，且分別滯后于Ta、VPD和SWC 3.3、4.6和10.72 h。

2）季節(jié)尺度上，油蒿可能通過提高PSII修復(fù)能力、增加特定蛋白數(shù)量和色素濃度等方式調(diào)節(jié)PSII能量分配，對Ta和VPD變化較敏感，并且受物候期的影響。PSII能量分配參數(shù)分別滯后于Ta、VPD、PAR和SWC 7.2、8.8、14.7和21.6 d。

3）在荒漠惡劣環(huán)境下，油蒿可以通過自身調(diào)節(jié)不斷修復(fù)受損PSII（Fv/Fm維持在0.78附近），對干旱和高輻射表現(xiàn)出耐受性，具有一定的適應(yīng)能力和生理可塑性。