熱帶北緣地區(qū)橡膠林水分利用效率的季節(jié)差異及歸因分析

郭沈沈，吳志祥，劉文杰，趙 豪，李運(yùn)帥，王譽(yù)博，王 鵬

（1.海南大學(xué) 生態(tài)與環(huán)境學(xué)院，?？?570228; 2.中國(guó)熱帶農(nóng)業(yè)科學(xué)院 橡膠研究所/ 海南儋州熱帶農(nóng)業(yè)生態(tài)系統(tǒng)國(guó)家野外科學(xué)觀測(cè)研究站，海南 儋州 571737）

水資源是地球-生物-化學(xué)-生態(tài)系統(tǒng)的核心載體，水資源空間差異是形成生態(tài)系統(tǒng)植被區(qū)域特征、影響總初級(jí)生產(chǎn)力分異的關(guān)鍵因素。水分利用效率（WUE）指的是植物光合作用中消耗單位質(zhì)量水分所固定的光合產(chǎn)物的總量，通常以總初級(jí)生產(chǎn)力（GPP）與對(duì)應(yīng)時(shí)期蒸散（ET）水分損失的比率[1-2]來(lái)表示，在生態(tài)系統(tǒng)水平上，WUE能夠反映植物的水分利用策略[3-4]，以及植物在光合碳同化過程中獲取的碳與消耗的水之間的權(quán)衡，即所謂的碳水耦合特征[5-7]。WUE的高低主要由植被的生產(chǎn)力和蒸散量決定，即表示碳封存和水資源消耗之間的關(guān)系，是生態(tài)系統(tǒng)碳和水循環(huán)耦合的重要生態(tài)指標(biāo)[8-9]。生態(tài)系統(tǒng)WUE的時(shí)間動(dòng)態(tài)主要受自然地理、植物功能類型、干擾歷史、氣候和環(huán)境因素以及植物控制氣孔的生理能力等因素的影響[10-11]。環(huán)境因素、生物因素及二者的交互作用以不同方式影響GPP和ET，進(jìn)而導(dǎo)致生態(tài)系統(tǒng)WUE發(fā)生變化[12]。鑒于持續(xù)的氣候變化，對(duì)生態(tài)系統(tǒng)WUE的深入了解將提高預(yù)測(cè)碳循環(huán)和水循環(huán)的能力，為區(qū)域水碳管理和生態(tài)環(huán)境的可持續(xù)發(fā)展提供參考[13-14]。

作為熱帶地區(qū)的“水泵”，橡膠林的蒸散量高于附近天然林的蒸散量，甚至?xí)?dǎo)致種植區(qū)雨季的土壤儲(chǔ)水量不足以維持高蒸散量而引起橡膠園旱季缺水，有研究認(rèn)為，種植區(qū)域在旱季面臨的水資源短缺問題主要是橡膠樹的“水泵”功能引起的[15-18]。因此，橡膠種植園的碳、水循環(huán)之間的耦合過程越來(lái)越受到人們的關(guān)注。為增加對(duì)橡膠園WUE的年度格局變化、年際趨勢(shì)以及WUE調(diào)節(jié)因素的認(rèn)識(shí)，Lin等[19]對(duì)云南省橡膠園進(jìn)行研究后發(fā)現(xiàn)，年際尺度上WUE與光合有效輻射呈負(fù)相關(guān)，并且WUE更多地受溫度、凈輻射和胸徑增長(zhǎng)率的控制，其中胸徑增長(zhǎng)率更能解釋W(xué)UE的變化狀況，說(shuō)明橡膠園碳、水循環(huán)耦合過程的主要影響因素是碳封存而非蒸散耗水量；在Tan等[20]的研究中，溫度是橡膠園WUE季節(jié)性變化的主要控制因素，太陽(yáng)輻射次之，葉面積指數(shù)隨溫度升高而增加，也對(duì)WUE具有重要影響，飽和水汽壓差（VPD）在短時(shí)間內(nèi)對(duì)WUE的影響較為明顯。目前，我國(guó)關(guān)于橡膠林WUE的研究較少[19-20]，且研究多聚焦于WUE與某些影響因素的變化規(guī)律及相關(guān)性、多基于定性的方式分析WUE的驅(qū)動(dòng)機(jī)制。例如，Li等[21]分析不同尺度下WUE與環(huán)境因素的關(guān)系后發(fā)現(xiàn)，在生態(tài)系統(tǒng)層面，葉面積指數(shù)是WUE變化的主要驅(qū)動(dòng)因素，溫度通過影響葉片蒸騰和光合作用，對(duì)WUE有很強(qiáng)的控制作用，即溫度的升高導(dǎo)致蒸騰作用增強(qiáng)，葉片固有WUE下降，生態(tài)系統(tǒng)WUE增加；除此之外，土壤水分也是WUE的重要影響因素，Nie等[17]的研究表明，生態(tài)系統(tǒng)WUE與土壤含水量之間存在顯著負(fù)相關(guān)的關(guān)系，且VPD越高負(fù)相關(guān)越明顯。主要是由于在植被生長(zhǎng)高峰期，隨著土壤水分的增加，ET增長(zhǎng)速度快于GPP，從而導(dǎo)致WUE隨土壤含水量的增加而降低。上述研究結(jié)果表明，植被冠層活動(dòng)和VPD在調(diào)控WUE中具有重要作用，如Zhao等[12]的研究結(jié)果表明， VPD和冠層導(dǎo)度是調(diào)控生態(tài)系統(tǒng)WUE最重要的因素，并在對(duì)北半球不同植被類型進(jìn)行研究后發(fā)現(xiàn)，WUE對(duì)氣溫、降水和VPD具有負(fù)響應(yīng)，對(duì)冠層導(dǎo)度和生態(tài)系統(tǒng)呼吸則呈正響應(yīng) 。然而，各因素對(duì)WUE季節(jié)性變化的具體貢獻(xiàn)尚不清楚，因此筆者擬在參考TRM的分析框架和基于EC-LUE模型、Penman-Monteith模型的基礎(chǔ)上，量化關(guān)鍵的大氣和地表生物物理因素在WUE季節(jié)性變化中的作用，以期為熱帶地區(qū)橡膠林生態(tài)系統(tǒng)保護(hù)、水資源合理利用提供理論依據(jù)。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況及數(shù)據(jù)本研究數(shù)據(jù)來(lái)源于2010年1月至2019年12月的渦度協(xié)方差測(cè)量數(shù)據(jù)，數(shù)據(jù)采集點(diǎn)位于海南省的農(nóng)業(yè)農(nóng)村部儋州熱帶作物科學(xué)觀測(cè)實(shí)驗(yàn)站（19°32′16″N，109°28′06″E），平均海拔114 m。研究區(qū)域位于海南省西北部，屬于典型熱帶海島季風(fēng)性氣候，具有明顯的雨旱季[22]。根據(jù)前人研究，筆者將該地區(qū)每年5—9月定義為雨季，10月至翌年4月為旱季[23]；2010—2019年平均溫度21.5～28.5 ℃；年平均降水量1 476.3 mm，且主要集中在7—9月份，占全年降水的70%以上。研究區(qū)域橡膠林植物群落組成簡(jiǎn)單，主要為橡膠林喬木層和林下草本層，冠層平均高度為20 m[24]。

于2009年11月起對(duì)該通量站點(diǎn)進(jìn)行近地面通量觀測(cè)，從地下1 m至觀測(cè)塔頂端50 m，測(cè)量數(shù)據(jù)包含地上和地下兩部分，具體測(cè)量?jī)x器及安裝見表1。

表1 海南省儋州市的農(nóng)業(yè)農(nóng)村部儋州熱帶作物科學(xué)觀測(cè)實(shí)驗(yàn)站渦度通量系統(tǒng)數(shù)據(jù)及測(cè)量?jī)x器安裝

1.2 水分利用效率的計(jì)算生態(tài)系統(tǒng)WUE在一定程度上反映了生態(tài)系統(tǒng)碳循環(huán)和水循環(huán)的耦合關(guān)系，并將光合作用和蒸散緊密聯(lián)系在一起[25]，通常采用GPP與ETET的比值進(jìn)行計(jì)算[10,12,25-26]：

式中，WUE為水分利用效率，單位為 g·kg-1，GPP為生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力，單位為g·m-2，ET為生態(tài)系統(tǒng)蒸散（橡膠林生態(tài)系統(tǒng)整體蒸散，包括冠層蒸騰量、裸地蒸發(fā)量和冠層截留水量），單位為kg·m-2。

GPP、ET及其中變量計(jì)算參考相關(guān)文獻(xiàn)[26-34]，將WUE公式計(jì)算解析為由不同影響因素組成的函數(shù)：

1.3 歸因方法及其優(yōu)化在本研究中，主要探討2個(gè)問題：（1）WUE雨季至旱季的變化狀況；（2）各大氣因素和生物物理因素在WUE變化中起到的作用及具體貢獻(xiàn)。在進(jìn)行分析時(shí)，將WUE雨季至旱季的變化視為由多個(gè)自變量變化所導(dǎo)致的結(jié)果，泰勒公式可以將復(fù)雜函數(shù)使用多項(xiàng)式來(lái)近似表達(dá)，因此利用多元函數(shù)的泰勒公式擬合出近似值（接近實(shí)測(cè)ET計(jì)算所得WUE變化值）。結(jié)合拉格朗日中值定理可知函數(shù)增加量、自變量增加量及導(dǎo)數(shù)之間的關(guān)系為：f(x+Δx)-f(x)=f'(x+θΔx)Δx，在此基礎(chǔ)上，參考TRM方法的分析框架[5,35]，量化關(guān)鍵大氣和地表生物物理因素在WUE變化中的作用，對(duì)其公式進(jìn)行一階求導(dǎo)可得：

式中，Δ表示雨旱季差值（例如，ΔG=G旱季-G雨季，表示G從雨季至旱季的變化量），偏微分項(xiàng)表示W(wǎng)UE對(duì)每個(gè)變量變化的敏感性（例如，表示W(wǎng)UE對(duì)G變化的敏感性）。

歸因方法是基于一階泰勒公式展開的，并未將高階項(xiàng)和交叉項(xiàng)考慮在內(nèi)，通常只有當(dāng)各歸因變量的變化很小時(shí)，這種方法引起的誤差才可以忽略不計(jì)，但本研究中部分影響因素的季節(jié)變化較大（表2）。為減少誤差，筆者參考加權(quán)平均方法來(lái)優(yōu)化結(jié)果[36]：

表2 雨季至旱季各影響因素的變化量

式中，X為模型中使用的最優(yōu)偏微分項(xiàng)，m為平均權(quán)重因子（取-200至200，以0.01為步長(zhǎng)，本研究m取值為0.04），x雨季和x旱季分別為僅使用雨季和僅使用旱季數(shù)據(jù)計(jì)算所得偏微分項(xiàng)。

2 結(jié)果與分析

2.1 ET、GPP和WUE的季節(jié)變化從圖1可知，同溫度（峰值位于雨季，且雨季平均溫度高于旱季6.67 ℃）和降水（多年雨季平均日降水高達(dá)7.57 mm，遠(yuǎn)高于旱季平均日降水1.50 mm）的年內(nèi)變化趨勢(shì)類似，ET和GPP具有明顯的季節(jié)性差異（P＜0.01），整體呈“倒U”型曲線，峰值通常位于雨季中期，雨季ET日均為7.93 kg·m-2，旱季日均3.31 kg·m-2，相差4.62 kg·m-2（圖1-a，圖1-f）；GPP雨季均值為20.89 g·m-2，旱季均值為10.83 g·m-2，相差10.06 g·m-2（圖1- b，圖1-g）；整體而言，WUE旱季的波動(dòng)幅度較大，平均值約為4.73 g·kg-1，雨季均值約為3.52 g·kg-1，雨旱季差異顯著（P＜0.01），旱季W(wǎng)UE平均高于雨季1.21 g·kg-1（圖1-c，圖1-h）。

圖1 橡膠林蒸散（ET）、生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力（GPP）、水分利用效率（WUE）、溫度（Ta）、降水（Pr）年內(nèi)變化及雨旱季對(duì)比

2.2 環(huán)境因素的季節(jié)差異由各環(huán)境因素雨旱季的對(duì)比及顯著性差異分析結(jié)果（圖2）可知，各因素均具有顯著性雨旱季差異（P＜0.01），其中，溫度、凈輻射、土壤熱通量、VPD、冠層光合有效輻射吸收比例、空氣動(dòng)力學(xué)阻力、太陽(yáng)輻射的雨季平均值高于旱季，地表阻抗的雨季平均值低于旱季；2010—2019 年雨季溫度平均值為29.11 ℃，旱季為22.44 ℃，相差6.67 ℃；凈輻射雨季均值為359.13 W·m-2，旱季為249.45 W·m-2，相差109.68 W·m-2；雨季平均土壤熱通量為9.48 W·m-2，旱季平均土壤熱通量為5.86 W·m-2，相差3.62 W·m-2；VPD的雨季平均值1.28 kPa，旱季平均值0.78 kPa， 差值為0.50 kPa；雨季植被冠層光合有效輻射吸收比例均值為0.67，旱季為0.63，雨季高于旱季0.04；太陽(yáng)輻射雨季均值為445.61 W·m-2，旱季為326.01 W·m-2，相差119.6 W·m-2；空氣動(dòng)力學(xué)阻力雨季均值32.14 s·m-1，旱季均值25.19 s·m-1，雨季高于旱季7.14 s·m-1；地表阻抗雨季均值為149.03 s·m-1，旱季均值260.46 s·m-1，雨季低于旱季111.43 s·m-1。

圖2 溫度（Ta）、凈輻射（Rn）、土壤熱通量（G）、飽和水汽壓差（VPD）、冠層光合有效輻射吸收比例（fPAR）、太陽(yáng)輻射（Sin）、空氣動(dòng)力學(xué)阻力（ra）、地表阻抗（rs）的雨旱季對(duì)比條形分別代表雨旱季多年日均值，誤差棒分別代表變量雨旱季的標(biāo)準(zhǔn)偏差。

2.3 歸因分析各影響因素對(duì)WUE變化的貢獻(xiàn)可見圖3。將從雨季至旱季的WUE變化分解為5個(gè)大氣因素（溫度、大氣壓、飽和水汽壓差、凈輻射、太陽(yáng)輻射）和4個(gè)生物物理因素（植被冠層光合有效輻射吸收比例、土壤熱通量、空氣動(dòng)力學(xué)阻力、地表阻抗）的貢獻(xiàn)，如圖3所示，第Ⅰ部分顯示模型估算GPP和模型估算ET計(jì)算所得ΔWUE與模型估算GPP和實(shí)測(cè)ET計(jì)算所得ΔWUE是相似的，WUE在旱季增加。第Ⅱ部分可知，大氣因素總貢獻(xiàn)和生物物理因素總貢獻(xiàn)皆為正向，這說(shuō)明，進(jìn)入旱季后，大氣因素和生物物理因素都傾向于增加WUE，并且生物物理因素的貢獻(xiàn)高于大氣因素的貢獻(xiàn)。

圖3 各影響因素對(duì)WUE變化量的貢獻(xiàn)量化圖

在大氣因素中（圖3中第Ⅲ部分，橙色條），溫度、VPD、凈輻射產(chǎn)生正增益（即三者的變化傾向于增加WUE），而太陽(yáng)輻射的變化傾向于減少旱季W(wǎng)UE（表3），并且被溫度、VPD和凈輻射的正增益所抵消，大氣壓的影響可忽略不計(jì)。

表3 WUE對(duì)大氣和生物物理因素變化的敏感性

在生物物理因素中（圖3中的第Ⅲ部分，藍(lán)色條），植被冠層光合有效輻射吸收比例的雨旱季差別不大（圖2-e, 表2），但WUE對(duì)其變化的敏感性較大，故產(chǎn)生負(fù)增益，傾向于減少旱季W(wǎng)UE（表3）；空氣動(dòng)力學(xué)阻力旱季的減少導(dǎo)致WUE（表3）的增加；地表阻抗的變化傾向于增加旱季W(wǎng)UE（表3），并抵消了植被冠層光合有效輻射吸收比例的負(fù)增益。土壤熱通量對(duì)WUE的貢獻(xiàn)則小到可以忽略不計(jì)。

第I部分中ΔWUEm代表GPP與模型ET計(jì)算所得WUE雨旱季差值；ΔWUE代表GPP與實(shí)測(cè)ET計(jì)算所得WUE雨旱季差值）。第II部分中，fa為大氣因素的總貢獻(xiàn)量，fs為生物物理因素的總貢獻(xiàn)量；第III部分中，橙色代表各大氣因素的貢獻(xiàn)，藍(lán)色代表各生物物理因素的貢獻(xiàn)。

3 討 論

3.1 橡膠林WUE的變化特征本研究結(jié)果顯示，橡膠林WUE的年內(nèi)變化與水分充足的中國(guó)東部森林保持一致，即在GPP和ET較高的季節(jié)，WUE較低[37]。在本研究中，橡膠林WUE的年內(nèi)變化幅度較小，多年旱季W(wǎng)UE平均高于雨季1.21 g·kg-1，與同處于熱帶地區(qū)的鼎湖山針闊葉混交林雨旱季W(wǎng)UE差值相近[38]，平均日均值為4.16 g·kg-1，略高于云南地區(qū)熱帶森林WUE3.16 g·kg-1，旱季水分利用效率的增加則是ET比GPP更強(qiáng)烈的下降引起的[20]， 結(jié)合本研究結(jié)果，從雨季過渡至旱季，ET降低幅度58.26%，GPP減少48.16%。橡膠樹在雨季從上層土壤中提取水分，在旱季從深層土壤中提取水分。依次來(lái)響應(yīng)土壤水分的變化[25]。對(duì)深層土壤水分的利用使其在旱季仍能提供一定水分以供光合作用固碳活動(dòng)的正常進(jìn)行，因此，旱季橡膠林生態(tài)系統(tǒng)的GPP下降趨勢(shì)較緩。

另一方面，根據(jù)降水的年內(nèi)分布可知，研究區(qū)域的降水集中分布在7—9月份，且雨熱同期，年降水量達(dá)1 476.3 mm，雨季降水的增加、氣溫的上升導(dǎo)致ET明顯升高，但GPP增加有限。在同處于熱帶的亞馬遜雨林生態(tài)系統(tǒng)的研究中，已有結(jié)果表明ET 是森林WUE變化的主導(dǎo)因子[18]，即相同增幅的降水量對(duì)GPP的促進(jìn)作用較小，但大幅度增加蒸散用水量，使得該地區(qū)雨季的WUE低于旱季[39]。在同處于熱帶地區(qū)的研究中，鼎湖山針闊葉混交林WUE變化出現(xiàn)相似情況，即降雨量集中的濕季W(wǎng)UE低于干季[38]。

3.2 橡膠林WUE季節(jié)變化的主要影響因素歸因分析結(jié)果表明，WUE的變化是大氣因素和生物物理因素共同作用的結(jié)果。隨著旱季的到來(lái)，溫度的降低對(duì)WUE的增加產(chǎn)生正增益，符合已有結(jié)論，在對(duì)西雙版納橡膠林的研究中，蒸散對(duì)溫度降低更為敏感，且溫度 和WUE之間存在負(fù)相關(guān)的關(guān)系[20,40-41]；旱季VPD的減少引起WUE增加，主要原因是VPD和ET的正相關(guān)性，這一點(diǎn)不僅在Zhao等[12]對(duì)北半球不同森林類型的研究中得以體現(xiàn)，而且在亞馬遜熱帶雨林的研究中也有類似結(jié)論[29]；凈輻射量的降低更傾向于減少蒸散，因而旱季凈輻射的減少導(dǎo)致WUE的增加。旱季入射短波輻射的降低會(huì)導(dǎo)致植物光合有效輻射減少，光合作用活動(dòng)減弱，從而引起WUE的減少（產(chǎn)生負(fù)增益）。

旱季植被冠層光合有效輻射吸收比例的減少意味著植被冠層固定能量減少，植被GPP降低，導(dǎo)致WUE的減少；但是，植被冠層光合有效輻射吸收比例帶來(lái)的負(fù)增益完全被地表阻抗抵消。從歸因分析中的生物物理要素來(lái)看，地表阻抗的貢獻(xiàn)量最高。以往的研究結(jié)果表明，冠層導(dǎo)度和生態(tài)系統(tǒng)呼吸是影響森林生態(tài)系統(tǒng)WUE的2個(gè)關(guān)鍵生物因素[12]。氣孔是冠層與大氣進(jìn)行氣體交換的主要通道，并且可以有效控制水分流失與植物碳的固定及有機(jī)物生產(chǎn)之間的平衡[42-43]。在同樣位于熱帶地區(qū)研究中，亞馬遜熱帶雨林旱季氣孔關(guān)閉通常會(huì)導(dǎo)致比光合作用更強(qiáng)烈的蒸騰作用下降[44]。本研究結(jié)果也與之相符，進(jìn)入旱季后，氣孔導(dǎo)度下降，地表阻抗增加，從而導(dǎo)致WUE增加。結(jié)合現(xiàn)有研究認(rèn)為，氣孔是控制植物蒸騰作用的最直接的因素[48]，其開閉活動(dòng)對(duì)植物蒸騰作用的直接影響遠(yuǎn)大于其對(duì)光合作用生產(chǎn)有機(jī)物的影響，橡膠林主要通過氣孔調(diào)節(jié)來(lái)應(yīng)對(duì)進(jìn)入旱季后供水量的減少。