熱帶季節(jié)性濕潤林苦楝(Melia azedarach)徑向生長季節(jié)動態(tài)及其對環(huán)境因子的響應(yīng)

吳麗杰,Arisa Kaewmano,付培立,王文禮,范澤鑫,*

1 云南大學(xué)生態(tài)與環(huán)境學(xué)院, 昆明 650091 2 中國科學(xué)院西雙版納熱帶植物園熱帶森林生態(tài)學(xué)重點(diǎn)實驗室, 勐臘 666303 3 中國科學(xué)院大學(xué), 北京 100049

作為陸地生態(tài)系統(tǒng)的重要組成部分,樹木生長對環(huán)境變化十分敏感,能迅速對環(huán)境變化做出響應(yīng)[1-4]。溫度的升高和降水的變化可能會以正面或負(fù)面的方式影響樹木季節(jié)性生長模式。樹木徑向生長對氣候因子的響應(yīng)是樹輪氣候?qū)W研究的基礎(chǔ),然而傳統(tǒng)的樹木年輪研究僅能在年際時間尺度上評估環(huán)境因子對樹木生長的研究,缺乏樹木生長的生理基礎(chǔ)[5]。

高精度生長儀(High resolution dendrometer)有助于在沒有破壞性采樣的情況下,以高分辨率(微米級)連續(xù)自動監(jiān)測樹木徑向生長的年內(nèi)變化過程,有助于從生理生態(tài)學(xué)角度深入認(rèn)識不同時間尺度上樹木生長季節(jié)動態(tài)規(guī)律及其對環(huán)境因子的響應(yīng),近年來得到廣泛應(yīng)用[6-11]。在溫帶和高海拔地區(qū)的研究發(fā)現(xiàn),溫度是限制生長開始的最重要因素[12-13]。對青海云杉(Piceacrassifolia)徑向生長季節(jié)動態(tài)的研究發(fā)現(xiàn),土壤溫度是限制樹木生長開始的主要因素[14]。我國北方半干旱地區(qū)的常綠樹種側(cè)柏(Platycladusorientalis)的徑向生長主要受到水分的限制[15]。對祁連山脈祁連圓柏(Sabinaprzewalskii)的研究發(fā)現(xiàn),不同海拔梯度的祁連圓柏的徑向生長主要受到降水的限制,但在生長起止時間、生長速率及生長期長度上有較大的差異[16]。而對熱帶地區(qū)的研究也發(fā)現(xiàn)降水和土壤含水量對樹木的生長起著重要的作用[17-18]。對熱帶喀斯特森林樹木生長的研究表明,樹木生長主要受到降水和相對濕度的限制,且生長季初期的降水對樹木的生長尤為重要[19]。

喀斯特(巖溶)地貌在我國西南廣泛分布,云南西雙版納地區(qū)的喀斯特地貌面積高達(dá)19%。喀斯特地貌由于土層淺薄,土壤保水能力差,許多植物生長在裸露的巖石上,即使在多雨的濕潤地區(qū),植物也遭受頻繁受到干旱脅迫。受熱帶季風(fēng)氣候的影響,西雙版納地區(qū)一年中有近6個月的旱季[20]。熱帶季節(jié)性濕潤林主要分布與石灰?guī)r山中部,是滇南石灰?guī)r森林的代表植被類型,落葉樹種在喬木種類中約占30%[21]。研究喀斯特生境樹木年內(nèi)徑向生長動態(tài)及其對環(huán)境因子的響應(yīng),可為了解熱帶季節(jié)性濕潤林植物的生長策略及其對季節(jié)性干旱的適應(yīng)提供參考。本研究選取云南西雙版納熱帶季節(jié)性濕潤林的落葉樹種苦楝(Meliaazedarach)為研究對象,楝科植物多具有清晰的年輪,且在東南亞熱帶地區(qū)廣泛分布,在熱帶樹木年輪研究中具有較好的價值。研究選擇的苦楝,是熱帶喀斯特生境中的常見樹種,多分布于土層淺薄的石灰?guī)r山,旱季落葉。利用高精度樹木徑向生長測量儀,在2018和2019兩個年度測定了苦楝的徑向生長季節(jié)動態(tài),并與同步監(jiān)測的環(huán)境因子進(jìn)行相關(guān)分析,旨在 (1)確定苦楝(Meliaazedarach)徑向生長的年內(nèi)動態(tài)規(guī)律及其年際間差異,(2)分析苦楝樹木徑向生長與氣候因子的響應(yīng)關(guān)系,明確苦楝徑向生長物候特征和主要限制氣候因子。

1 材料與方法

1.1 研究樹種及樣地

本研究地點(diǎn)位于西雙版納綠石林森林公園(21°54′N,101°46′N,海拔680 m),是云南省西雙版納國家級自然保護(hù)區(qū)的一部分。公園內(nèi)石灰?guī)r叢生,土壤為石灰?guī)r基質(zhì)上發(fā)育的磚紅壤,呈斑塊狀散布于裸露巖石之間,地表崎嶇呈喀斯特地貌,植被屬于喀斯特?zé)釒Ъ竟?jié)性濕潤林,含有一定比例旱季落葉成分。西雙版納地區(qū)年均溫21.7 ℃,年平均降雨量1560 mm,存在明顯的干、濕之分,每年的5—10月為雨季,11月至翌年4月為干季,其中11月到次年2月為霧涼季,早晚有濃霧,彌補(bǔ)了降水的不足[22]。其中雨季降雨量占全年的83%,干季降雨量不足20%。本研究選取熱帶季節(jié)性濕潤林中自然分布的落葉樹種苦楝(Meliaazedarach)作為研究樹種。

1.2 樹木徑向生長動態(tài)監(jiān)測

本研究選取了在胸徑、樹高和生長狀況等方面特征相似的4株苦楝(Meliaazedarach)植株(平均胸徑為 (32±5) cm),在選取植株的胸高(<1.3 m)處安裝高精度生長儀(精度：±2 μm,DR型,德國Ecomatik),連續(xù)兩年(2018—2019年)持續(xù)監(jiān)測了苦楝的樹干徑向變化,以10 min為間隔自動記錄并保存數(shù)據(jù)。為了減少樹皮膨脹和收縮過程對樹干徑向監(jiān)測的影響,在安裝生長儀時剝?nèi)ゲ糠謽淦ぁ?/p>

1.3 氣象數(shù)據(jù)的獲取

氣象數(shù)據(jù)使用中國科學(xué)院西雙版納熱帶雨林生態(tài)系統(tǒng)研究站(21°55′N,101°16′E)監(jiān)測數(shù)據(jù),該氣象站距離森林公園監(jiān)測樣點(diǎn)3.5 km。氣象因子包括空氣溫度(Air Temperature,T)、相對濕度(Relative humidity,RH)、降水量(Precipitation,Prec)、光合有效輻射(Photosynthetic active radiation,PAR)、土壤溫度(Soil Temperature,ST),土壤含水量(Soil water content,SWC)、風(fēng)速(Wind speed,WS)等。氣象數(shù)據(jù)均為日平均值。根據(jù)空氣溫度(平均)和空氣相對濕度計算得出飽和水汽壓虧缺(VPD)：

式中,Ta為空氣溫度,RH為相對濕度。a,b,c是常數(shù),分別為0.611 kPa,17.052 (無單位) 和240.97 ℃。

1.4 數(shù)據(jù)分析

使用Gompertz函數(shù)來擬合樹木的徑向生長變化,確定兩個研究樹種的生長季[23]。Gompertz函數(shù)公式如下：

Y=Y0+Aexp(-eβ-kt)

式中,Y表示日徑向累積增長量;Y0是下漸近線;A為上漸近線;β為x軸的截距參數(shù);κ為變化速率參數(shù),其大小會影響生長季的長短;t為時間,以儒略日(DOY)表示,儒略日(Day of year, DOY)為一年中的一天。因此,(A-Y0)表示季節(jié)性生長總量,而(Yt-Yt-1)對應(yīng)于增長速率(mm/d),當(dāng)達(dá)到總生長量的10%和90%時,確定為生長開始和生長停止的時間[24]。

高精度樹木生長儀記錄的數(shù)據(jù)不僅包括了由生長引起的樹干的不可逆膨脹,還包括了由于水分變化引起的樹木的可逆收縮。根據(jù)Downes等的樹干日循環(huán)法[25]可將樹木的徑向生長劃分為3個不同的階段：(1)收縮階段,即當(dāng)日最大值到當(dāng)日最小值之間的時間段,(2)恢復(fù)階段,樹干直徑從當(dāng)日最小值到次日最大值之間的時間段,(3)增長階段,樹干直徑超過早晨的最大值到隨后的最大值之間的部分?jǐn)U張階段。從中提取樹干的最大日收縮和日凈增長量,樹干日凈增長量(Daily radius increment,DR)為連續(xù)兩日的最大值之差,而樹干最大日收縮(Maximum daily shrinkage,MDS)為當(dāng)日樹干直徑最大值和最小值之差。

將2018—2019年生長季 (3—10月)的樹干徑向增長量與氣候因子做相關(guān)分析,利用R軟件的nlme軟件包建立線性混合效應(yīng)模型 (Linear mixed-effects modeling,LMM),將不同樣木個體作為隨機(jī)因子,樹木徑向日生長量為因變量,氣候因子(降水、空氣溫濕度、土壤溫度、20cm深度土壤含水量、光合有效輻射)作為自變量,通過標(biāo)準(zhǔn)化回歸系數(shù)評估分析不同氣候因子對樹木徑向生長的影響,以確定主要限制因子。統(tǒng)計分析利用R統(tǒng)計軟件完成[26]。

2 結(jié)果與分析

2.1 環(huán)境因子的變化

研究地點(diǎn)的氣候具有明顯的季節(jié)性特征,旱季和雨季氣候變化差異明顯。圖1展示了2018—2019年兩年的氣候因子季節(jié)動態(tài)。2018年的日平均氣溫為 22.3 ℃,極端最高溫為36 ℃,2018年年降水總量為1899.4 mm,略高于該地區(qū)的多年平均降水量。2019年年降水總量為832.6 mm,2019年年均氣溫為23 ℃略高于2018年,2019年極端最高溫達(dá)42.4 ℃(5月20日),而極端最低氣溫為5.1 ℃(12月10日)。2019年5月溫度為1959年以來最高,大于0.1 mm降水量為歷史以來最低。由于2019年高溫少雨,導(dǎo)致相對濕度和水汽壓虧缺的波動幅度遠(yuǎn)遠(yuǎn)大于2018年,并在2019年5月20日出現(xiàn)極值(相對濕度：56%,水汽壓虧缺：1.75 kPa)。

圖1 西雙版納勐侖地區(qū)2018—2019年溫度、降水、相對濕度、水汽壓虧缺和光合有效輻射因子變化(數(shù)據(jù)由西雙版納生態(tài)站提供)Fig.1 Seasonal variations of temperature, rainfall, relative humidity, vapour pressure deficit, photoactive radiation over the 2018—2019 in Menglun, Xishuangbanna (The data provided by Xishuangbanna Ecological Station)

2.2 樹干徑向季節(jié)變化

利用Gompertz函數(shù)對苦楝的徑向生長進(jìn)行擬合(圖2),結(jié)果表明：2018年苦楝的最大生長速率集中在5月中下旬(DOY134-144),而2019年苦楝的最大生長速率集中在6月初到7月中旬(DOY158-200)。,苦楝在2019年的年總生長量 (0.90 mm)略低于2018年的年總生長量(0.95 mm),但2019年苦楝的最大生長速率要高于2018年(圖2和圖3)。由于兩個生長季氣候差異,苦楝開始生長和結(jié)束生長的時間也存在不同,2018年苦楝徑向生長開始的時間為3月初(DOY64),2019年開始生長的時間為6月中旬(DOY163),2019年苦楝開始生長的時間要比2018年晚,可能與2019年4—5月的高溫干旱有關(guān)。2018年苦楝結(jié)束生長的時間為DOY155要早于2019年結(jié)束生長的時間(DOY242)。

圖2 苦楝(Melia azedarach)2018和2019年樹干徑向變化和Gompertz函數(shù)模擬的日生長速率Fig.2 Stem radius variation and daily growth rates modeled by Gompertz function of Melia azedarach for year 2018 and 2019LS-M6, LS-M7, LS-M8, LS-M9為研究樣樹樹號

圖3 苦楝(Melia azedarach)2018和2019年徑向生長物候參數(shù)箱型圖Fig.3 Boxplots of growth phenology parameters of Melia azedarach for year 2018 and 2019

2.3 樹干徑向季節(jié)對氣候因子的響應(yīng)

圖4 苦楝(Melia azedarach)2018年和2019年樹干徑向日凈增長量與氣候因子的Spearman相關(guān)分析Fig.4 Spearman correlation coefficients between daily growth rates of Melia azedarach and climatic factors for year 2018 and 2019圖中灰色短虛線為P<0.05顯著水平線;Prec：降水量PreciPitation;RH：相對濕度Relative humidity;SWC20：20cm土壤含水量Soil water content at 20cm;Ta_max：日最高氣溫 Daily maximum air temPerature;Ta_mean：日平均氣溫Daily mean air temPerature;Ta_min：日最低氣溫Daily minimum air temPerature;Ts20：20cm土壤溫度Soil temPerature at 20cm;PAR：光合有效輻射Photosynthetic active radiation;VPD：水汽壓虧缺VaPor Pressure deficit;WS：風(fēng)速Wind SPeed

將2018—2019年生長季期間苦楝的樹干直徑日凈增長量(DR)與氣候因子做相關(guān)分析(圖4),結(jié)果表明：2018年徑向日凈增長量與降水和相對濕度呈正相關(guān)關(guān)系,而與土壤含水量、空氣溫度(最高/最低/平均溫)、土壤溫度、光合有效輻射、水汽壓虧缺和風(fēng)速呈負(fù)相關(guān)關(guān)系;而2019年徑向日凈增長量與降水、相對濕度、土壤含水量、最低/平均氣溫和土壤溫度呈正相關(guān)關(guān)系,與最高氣溫、光合有效輻射、水汽壓虧缺和風(fēng)速呈負(fù)相關(guān)關(guān)系。不管是降水較多的濕潤年份(2018年),還是降水較少的干旱年份(2019年),苦楝的徑向日凈增長量均與降水和相對濕度呈正相關(guān)關(guān)系,但這種相關(guān)性在干旱年份(2019年)更強(qiáng)(圖5)。在7天的時間尺度上,2018—2019年生長季期間苦楝的樹干直徑日凈增長量與降水和相對濕度響應(yīng)趨勢與日尺度上一致,徑向生長量與降水和相對濕度仍然是正相關(guān)關(guān)系(圖6)。

圖5 苦楝(Melia azedarach)2018和2019年日生長量與降水和相對濕度的散點(diǎn)圖Fig.5 Scatter Plots between daily growth rates of Melia azedarach and preciPitation and relative humidity for year 2018 and 2019

圖6 苦楝(Melia azedarach)2018年和2019年樹干徑向凈增長量(7d)與氣候因子的Spearman相關(guān)分析Fig.6 Spearman correlation coefficients between weekly radial growth rates of Melia azedarach and climatic factors for year 2018 and 2019

利用線性混合效應(yīng)模型分析各個氣候因子對苦楝樹木徑向生長的貢獻(xiàn),通過模型選擇剔除部分不限制的氣候因子后,分別選擇出2018年和2019年顯著的氣候因子(圖7),結(jié)果表明：2018年徑向日凈增長量與降水和日平均氣溫呈正相關(guān)關(guān)系,與光合有效輻射和土壤含水量呈負(fù)相關(guān)關(guān)系;2019年徑向日凈增長量與降水、相對濕度和土壤溫度呈正相關(guān)關(guān)系,與土壤含水量呈負(fù)相關(guān)關(guān)系。不管是降水多的年份還是降水較少的年份,熱帶季節(jié)性濕潤林中苦楝的徑向生長都與降水有穩(wěn)定的正相關(guān)關(guān)系,說明苦楝的徑向生長主要受到了水分的影響。

圖7 苦楝(Melia azedarach)2018和2019年日生長量與氣候因子的混合效應(yīng)模型的標(biāo)準(zhǔn)化系數(shù)Fig.7 Standardized coefficients for the linear mixed-effects models to exPlain the variations of daily growth rates of Melia azedarach and climatic factors for year 2018 and 2019

3 討論

研究結(jié)果表明,苦楝的生長季較短,2018年苦楝徑向生長主要集中在3月初到6月初;而2019年苦楝徑向生長主要集中在6月中旬到8月底。本研究中苦楝的生長季較短可能與石灰?guī)r生境較為嚴(yán)重的干旱脅迫有關(guān)。由于我們計算徑向生長開始的時間是參照的是生長達(dá)到最大生長量的10%作為標(biāo)準(zhǔn),但實際上樹木形成層活動的時間可能會早于根據(jù)此標(biāo)準(zhǔn)計算的結(jié)果[27]。另外一方面,由于樹木形成層活動過程中還會經(jīng)歷細(xì)胞壁加厚這一過程,但這一階段較難通過生長儀監(jiān)測出來。由于年際間氣候的差異,樹木開始生長,結(jié)束生長以及持續(xù)生長的時間存在差異。2019年樹木開始生長和結(jié)束生長的時間都要晚于2018年,且2018年的生長持續(xù)時間要長于2019年,年總生長量也高于2019年(圖3)。熱帶季節(jié)性濕潤林苦楝的徑向生長在不同年份存在較大差異,研究地區(qū)在2019年生長季早期(4—6月)經(jīng)歷了嚴(yán)重的高溫干旱,降水量少,相對濕度低,水汽壓虧缺高,加上喀斯特土壤持水能力差,樹木可能經(jīng)歷嚴(yán)重的干旱脅迫,較低的土壤含水量不能滿足苦楝形成層細(xì)胞啟動生長所需的條件,因此徑向生長開始較2018年晚,年總生長量也較少。對熱帶北緣喀斯特森林的四個主要樹種的徑向生長研究發(fā)現(xiàn),徑向生長開始于旱季末期,在雨季加速,10月以后生長減緩甚至停止[19]。但在對厄瓜多爾南部不同熱帶森林中的Tabebuiachrysantha研究發(fā)現(xiàn),形成層活動期的長度不是確定絕對生長量的唯一因素[18]。

我們的研究發(fā)現(xiàn)不管是降水多的年份還是降水較少的年份,苦楝的徑向生長與降水和相對濕度呈正相關(guān)關(guān)系,與光合有效輻射(PAR)、水汽壓虧缺(VPD)和風(fēng)速(WS)呈負(fù)相關(guān)關(guān)系,這些環(huán)境因子都是直接(Prec、RH)或間接(PAR、VPD、WS)的方式影響樹木的水分供應(yīng),從而影響樹木的徑向生長,這表明了在熱帶喀斯特地區(qū)樹木的徑向生長受到了水分的限制。García-Cervigón等人對厄爾瓜多南部的熱帶干旱森林中5種樹種的樹干徑向變化研究發(fā)現(xiàn),木質(zhì)部的生長主要受降水量的限制[28]。Hu和Fan對熱帶喀斯特森林的4個樹種研究同樣也得到了相似的結(jié)論,降水和相對濕度是限制樹木徑向生長的主要因子,其中生長初期的降水量尤為重要[19]。Muller等人發(fā)現(xiàn),水分通過影響細(xì)胞的膨壓進(jìn)而影響細(xì)胞分裂和伸長生長[29]。由于細(xì)胞的分裂和擴(kuò)大取決于細(xì)胞的膨大,水分虧缺可能直接影響到木質(zhì)部形成層的分化和細(xì)胞生長[30]。厄瓜多爾南部不同熱帶森林中樹木的生長與降雨量以及水汽壓虧缺之間存在正相關(guān)關(guān)系,表明降雨對樹木的生長起著重要的作用[18]。對亞熱帶卡西亞松樹木徑向生長與環(huán)境因子關(guān)系的研究,也證明降水和相對濕度是樹木徑向增長的最重要限制因素,而水汽壓虧缺會影響樹木的水分平衡[31]。對我國祁連山東部青杄(Piceawilsorii)年內(nèi)徑向生長對氣候因子的響應(yīng)研究發(fā)現(xiàn)了相似的結(jié)果,在不同的時間尺度上青杄的徑向生長與降水呈正相關(guān),與溫度和水汽壓虧缺呈負(fù)相關(guān),青杄的徑向生長受到了水分的限制[32]。細(xì)胞的擴(kuò)大促進(jìn)了樹木的徑向增長[33],而細(xì)胞的擴(kuò)大主要是因為細(xì)胞內(nèi)水分含量導(dǎo)致的。在本研究中,苦楝的徑向生長與降水和相對濕度呈正相關(guān)關(guān)系,土壤含水量增高有利于樹干夜間補(bǔ)充水分,使得形成層細(xì)胞處于正壓的環(huán)境,導(dǎo)致形成層細(xì)胞增大[34]。苦楝的徑向生長與日最高氣溫、光合有效輻射、水汽壓虧缺呈負(fù)相關(guān)關(guān)系,這些環(huán)境因子的升高會導(dǎo)致蒸騰速率升高,增大了水分運(yùn)輸系統(tǒng)的負(fù)壓,樹木形成層區(qū)域細(xì)胞膨壓降低[35],導(dǎo)致氣孔關(guān)閉[36],進(jìn)而影響了植物碳水化合物的形成和徑向生長[37]。在旱季長度和干旱嚴(yán)重程度不同的兩種熱帶干旱森林的研究也發(fā)現(xiàn),較高的氣溫對徑向增量速率產(chǎn)生了負(fù)面影響,這可能是由于水汽壓虧缺的增加和蒸發(fā)散增加所致[38]。

4 結(jié)論

本研究利用高精度生長儀連續(xù)兩年監(jiān)測了西雙版納石灰?guī)r山的熱帶季節(jié)性濕潤林苦楝的徑向生長動態(tài)。由于年際間氣候存在差異,苦楝生長開始、結(jié)束以及持續(xù)生長的時間在不同年份之間存在較大差異?？嚅膹较蛏L與降水和相對濕度呈正相關(guān)關(guān)系,與光合有效輻射、水汽壓虧缺和風(fēng)速呈負(fù)相關(guān)關(guān)系,表明了在熱帶石灰?guī)r生境樹木的徑向生長受到了水分的限制。在極端干旱的年份(2019年),熱帶季節(jié)性濕潤林苦楝的徑向生長開始時間較晚,可能因為水分供應(yīng)不足導(dǎo)致的木質(zhì)部細(xì)胞分化延遲所致。在氣候變暖和區(qū)域干旱化的背景下,熱帶石灰?guī)r生境植物可能面對更為嚴(yán)重的干旱脅迫,結(jié)合徑向生長監(jiān)測、形成層活動和樹木光合水分生理測定,有望進(jìn)一步評估熱帶季節(jié)性濕潤林植物的生長策略和適應(yīng)性。楝科植物多具有清晰的年輪,且在東南亞熱帶地區(qū)廣泛分布,在熱帶樹木年輪研究中具有較好的價值,本文的研究結(jié)果可為深入認(rèn)識熱帶地區(qū)樹木生長的季節(jié)動態(tài)和氣候敏感性提供參考。

致謝：感謝馬宏和譚美毅在野外試驗中提供的幫助。感謝中國科學(xué)院西雙版納熱帶雨林生態(tài)系統(tǒng)研究站提供的氣象數(shù)據(jù)。