馬尾松光合作用光響應(yīng)曲線模型擬合

陳雪蓮,徐六一,郝焰平,章 健,姜春武,潘 婷

(安徽省林業(yè)科學研究院,安徽 合肥 230031)

馬尾松光合作用光響應(yīng)曲線模型擬合

陳雪蓮,徐六一*,郝焰平,章 健,姜春武,潘 婷

(安徽省林業(yè)科學研究院,安徽 合肥 230031)

為了探討不同的光響應(yīng)曲線模型對馬尾松無性系不同生長期光合特性的適用性,采用LI-6400便攜式光合儀器分別測定了7月、9月和11月馬尾松的光響應(yīng)曲線,并選用二次多項式回歸模型(QPM)、分段函數(shù)模型(SFM)、直角雙曲線模型(RHM)、非直角雙曲線模型(NRHM)、直角雙曲線修正模型(MRHM)和指數(shù)函數(shù)模型(EM)來進行最適光響應(yīng)曲線模型的選擇。結(jié)果表明：要準確獲得馬尾松不同生長期的光響應(yīng)參數(shù)需選擇不同的光響應(yīng)曲線模型,7月、9月和11月份馬尾松最適合的光響應(yīng)曲線模型分別是MRHM、EM和QPM。

馬尾松；光響應(yīng)曲線；模型擬合

Abstract： The photo-response curves ofMassonpinein July, September and November were respectively measured by using LI-6400 portable photosynthesis system. These photo-response curves were fitted by quadratic polynomial regression model (QPM), subsection function model (SFM), rectangular hyperbola model (RHM), non-rectangular hyperbola model (NRHM), modified rectangular hyperbola model (MRHM), and exponential function model (EM), respectively, and the optimum models were selected out based on their fitting effect. The results illustrated that different photo-response curve models should be selected to accurately obtain the photo-response parameters ofM.pinein different growth periods. The optimum models fitting the photo-response curves ofM.pinein July, September and November were MRHM, EM and QPM, respectively.

Keywords：Massonpine; Photo-response curve; Model fitting

光合作用是植物生長過程中的關(guān)鍵生理過程,光照強度是其主導(dǎo)因子。光響應(yīng)曲線是凈光合速率隨著光合有效輻射改變而變化的反應(yīng)曲線,通過對其進行測定,可以分析計算出最大光合速率(Pnmax)、光飽和點(LSP)、光補償點(LCP)、表觀量子效率(AQY)和暗呼吸速率(Rd)等重要參數(shù)來幫助了解植物光合機構(gòu)的運轉(zhuǎn)狀況、植物的光合效率及其對光的適應(yīng)性等生理機制[1-3]。為了獲得這些參數(shù),國內(nèi)外研究人員已建立了諸如二次多項式回歸模型、直角雙曲線模型、直角雙曲線修正模型、非直角雙曲線模型和指數(shù)模型等許多光響應(yīng)曲線模型[3]。這些模型都存在著對某一個或某幾個參數(shù)估計欠佳的缺點。由于不同的模型有不同的意義,而不同植物的生理特性也有差別,因此在研究對象植物的光合生理特性時要考慮模型的適用性,選擇最佳的擬合模型來估計光合參數(shù),提高估計參數(shù)的準確性以便更好地認識其光合生理特性[4]。

馬尾松(Pinusmassoniana)是我國特有的主要造林針葉樹種,分布廣、生長快、材質(zhì)好,是高產(chǎn)的木質(zhì)產(chǎn)品和高效的非木質(zhì)產(chǎn)品以及高耐逆境的優(yōu)良鄉(xiāng)土樹種,在我國的木材工業(yè)、造紙工業(yè)、林產(chǎn)化學工業(yè)、醫(yī)療保健品工業(yè)等生產(chǎn)原料的提供、荒山綠化等生態(tài)環(huán)境建設(shè)方面,都有其不可替代的作用,在我國社會經(jīng)濟建設(shè)中占居十分重要的地位[5]。目前在針對馬尾松光合生理方面的研究中,一些研究者在分析光響應(yīng)曲線時有的采用非直角雙曲線模型進行擬合[6],有的采用直角雙曲線修正模型進行擬合[7],雖然這些模型對馬尾松的擬合程度比較高,但都沒有對光響應(yīng)參數(shù)估計的合理性進行分析。研究光響應(yīng)模型對馬尾松的適用性,對提高馬尾松光響應(yīng)參數(shù)的合理性,提高同類型試驗結(jié)果的可比性有一定的科學意義和實用價值。本研究以馬尾松為研究對象,在不同的季節(jié)對其進行光響應(yīng)測定,檢驗6種典型的光響應(yīng)模型的擬合優(yōu)度,選擇最適的光合作用-光響應(yīng)模型,為馬尾松的栽培和撫育管理提供理論參考。

1 材料與方法

1.1試驗地概況

試驗地設(shè)在合肥市廬陽區(qū)大楊鎮(zhèn)國家林業(yè)高科技開發(fā)中心內(nèi)抗松材線蟲病馬尾松無性系保存圃,地處江淮之間,地理位置位于117°03′44″～117°25′33″ E, 31°41′54″～31°57′59″ N,為亞熱帶濕潤季風氣候,四季分明,氣候溫和,雨量適中,梅雨顯著,夏雨豐沛,年均溫度為15.7 ℃,極端最低氣溫-20.6 ℃,極端最高氣溫41.0 ℃,無霜期227 d,年平均降水量在988.4 mm。土壤為黃棕壤,質(zhì)地粘重,呈微酸性。

1.2試驗材料

材料來源于安徽省馬尾松抗松材線蟲病抗性育種項目選出的5年生馬尾松嫁接無性系。嫁接苗的接穗來自合肥市廬陽區(qū)大楊鎮(zhèn)國家林業(yè)高科技開發(fā)中心內(nèi)抗松材線蟲病馬尾松無性系保存圃中8年生同一種源的馬尾松實生苗,并經(jīng)過2次人工接種松材線蟲測定。

1.3光響應(yīng)曲線的測定

在7月份、9月份和11月份,選擇晴朗的天氣,在8：00～11：00采用LI-6400光合儀自帶的LED紅藍光源測定凈光合速率(Pn),形成Pn對光照強度的響應(yīng)曲線。在0～2000mol/(m2·s)范圍內(nèi),光照強度設(shè)定16個梯度,從低到高依次是0、50、80、100、120、150、200、400、600、800、1000、1200、1400、1600、1800、2000mol/(m2·s)。設(shè)定光照強度變化時最小穩(wěn)定時間120 s,測量結(jié)果變異率小于0.05時,儀器自動進行測定。測定時的CO2濃度為當前大氣CO2濃度(400±20)mol/mol2,葉溫跨度10～45 ℃,空氣濕度20%～60%。測定9個無性系,分為3組,每組3個重復(fù),每個重復(fù)測定樹木中部向陽部位5～6根松針,記錄3組穩(wěn)定數(shù)據(jù),最后的統(tǒng)計數(shù)據(jù)取平均值。

1.4光響應(yīng)曲線的擬合模型

選擇6個典型的光響應(yīng)曲線模型對Pn和I觀測值進行擬合。

(1)二次多項式回歸模型(QPM):

Pn=aI2+bI+c

(2)分段函數(shù)模型(SFM)[8]:

Pn=α′I+RdI≤150mol/(m2·s)

Pn=aI2+bI+cI>150mol/(m2·s)

(3)直角雙曲線模型(RHM)[9]：

(4)非直角雙曲線模型(NRHM)[10]：

(5)直角雙曲線修正模型(MRHM)[11]：

(6)指數(shù)函數(shù)模型(EM)[12]：

式中,Pn為凈光合速率；I為光合有效輻射；a、b、c均為系數(shù)；α為初始量子效率(即植物光合作用對光響應(yīng)曲線在I=0時的初始斜率)；α′為表觀量子效率；Pnmax為最大凈光合速率；Rd為暗呼吸速率；β為光抑制項,是修正系數(shù)；r為光飽和項,其值為α/Pnmax；θ(0<θ≤1)為光響應(yīng)曲線的曲角。

模型(1)在光響應(yīng)曲線研究中僅作為一個純數(shù)學模型被應(yīng)用,根據(jù)二次多項式的特性可以計算出最大凈光合速率[Pnmax=c-(b2/4a)]、暗呼吸速率(Rd)和光飽和點(LSP=-b/2a)；對曲線中光合有效輻射強度(PAR)小于150mol/(m2·s)的點進行線性擬合,與橫軸的交點即為光補償點(LCP)；模型(2)中,利用線性方程可得到α值和Rd值,并計算出光補償點(LCP=-Rd/α),再利用拋物線方程計算出LSP和Pnmax。模型(3)和(4)是沒有極值的函數(shù),無法根據(jù)方程求出LCP和LSP的解析解。因此只能利用直線方程擬合弱光下[I≤150mol/(m2·s)]的光響應(yīng)數(shù)據(jù),進行簡單線性回歸分析,求得的回歸曲線與直線y=Pnmax交點對應(yīng)的x坐標值即為LSP,然后再計算出LCP[13]。模型(5)是一個具有極值的函數(shù),上述其他幾個模型中包含的幾個參數(shù)都可通過該模型求解。若β>0,則光合參數(shù)LSP、LCP、Pnmax均可用該模型擬合所得參數(shù)求出；若β<0,則應(yīng)有(β+r)/β>0,否則亦無法直接求出LSP和Pnmax。模型(6)的LSP以光合速率達到最大值99%時的光照強度來計算,Pn=0時的光強為光補償點[14]。

1.5模型評價指標

為了比較所選用的6個光響應(yīng)曲線模型對馬尾松的擬合優(yōu)度,本研究除了選用決定系數(shù)(R2)來評價模型的擬合效果外,還選用了均方根誤差(RMSE)、平均相對誤差絕對值(MAE)和赤池信息量準則(AIC)3個參數(shù)[15-18]來進行評價:

1.6數(shù)據(jù)處理

本研究所得的實驗數(shù)據(jù)應(yīng)用Excel 2007進行初步的處理,光響應(yīng)曲線模型擬合通過SPSS 17.0軟件完成。

2 結(jié)果與分析

2.1馬尾松7月份的光響應(yīng)曲線及模型擬合比較

7月份馬尾松的實測凈光合速率對光合有效輻射強度的響應(yīng)曲線如圖1所示。凈光合速率隨光照強度的上升達到某個值后,凈光合速率則不再隨光照強度的上升而進一步上升,而是達到一種趨于飽和的狀態(tài),即為飽和趨近型[19]。利用二次多項式回歸(SFM)、分段函數(shù)(SFM)、直角雙曲線(RHM)、非直角雙曲線(NRHM)、修正的直角雙曲線(MRHM)和指數(shù)函數(shù)(EM)6種光響應(yīng)模型進行擬合,擬合曲線如圖1所示。從圖1中可以看出QMP模型的擬合值與實測值(MV)偏差比較大,在光照強度約600mol/(m2·s)以下至光補償點以上的階段,其Pn值都低于實測值,在光照強度達到約1800mol/(m2·s)時擬合曲線向下彎曲,雖然擬合值低于實測值,但較好地展現(xiàn)了在強光階段的光抑制現(xiàn)象。在光補償點以上的弱光階段以及光照強度達到約600mol/(m2·s)以上時其他幾種模型的擬合值與實測值都比較接近,只是在拐點的位置比實測值略低。

圖1 7月份不同模型的光響應(yīng)曲線

在表1中列出了6種模型擬合得到的光響應(yīng)參數(shù)值和實際測量值,并進行了擬合優(yōu)度的比較。由表1可知,6種模型中只有NRHM模型擬合的表觀量子效率α值不理想,達到了α的理論上限值0.125,這有悖于生理學意義。對于最大凈光合速率(Pnmax),QMP、SFM和MRHM模型的擬合值與實際值都比較接近,RHM、NRHM和EM的擬合值比實測值都要高。光飽和點的實測值約是1200mol/(m2·s),與之最接近的是QMP和MRHM模型的擬合值,SFM、NRHM和EM的擬合值都比實際值高,RHM的擬合值更遠高于實測值。光補償點的實測值約是55mol/(m2·s),QMP和MRHM的擬合值與其比較接近,其余的幾個模型擬合值要么遠大于實測值,要么遠小于實測值。

通過比較擬合優(yōu)度可以看出,QMP的決定系數(shù)R2最低,擬合程度最低；MRHM和EM模型的擬合程度最高,而且MRHM模型的其他幾個擬合優(yōu)度值也幾乎都是最低的。因此,通過綜合考察各模型的擬合值與實測值的接近程度以及擬合優(yōu)度的比較可以看出直角雙曲線修正模型(MRHM)是馬尾松7月份最適用的光響應(yīng)模型。

2.2馬尾松9月份的光響應(yīng)曲線及模型擬合比較

9月份馬尾松的實測凈光合速率對光合有效輻射強度的響應(yīng)曲線如圖2所示。凈光合速率隨光照強度的上升而升高,最后達到一種趨于飽和的狀態(tài),也屬于飽和趨近型。但是與7月份不同的是,在低光強階段Pn隨著PAR的升高迅速升高后并沒有迅速達到光飽和狀態(tài),而且經(jīng)過一個緩慢的過程后才趨于飽和。利用6種光響應(yīng)模型進行擬合,擬合曲線如圖2所示。從圖2中可以看出SFM和MRHM模型的擬合值與實測值(MV)偏差比較大,在給定的光照強度范圍內(nèi),SFM呈現(xiàn)一直升高的趨勢,沒有出現(xiàn)明顯的光飽和特征,與實際情況不符。RHM模型的擬合值在由弱光向飽和光強過渡的過程中與實測值有偏差,達到飽和光強后其值比實測值略低。其余幾個模型的擬合值與實測值都比較接近。

圖2 9月份不同模型的光響應(yīng)曲線

從表2中可以看出,MRHM模型在對9月份測定的光響應(yīng)數(shù)據(jù)進行擬合時,由于β=0,使得修正的直角雙曲線模型變成直角雙曲線模型,導(dǎo)致不能再利用公式直接求出Pnmax、LSP和LCP的值。排除MRHM模型,其他5種模型擬合的α值都比較理想。QMP和EM模型擬合的Pnmax值與實測值比較接近,SFM和RHM模型擬合的Pnmax值比實測值略大,NHM模型擬合的Pnmax值達到了50mol/(m2·s),遠大于實測值。馬尾松光飽和點的實測值是1800mol/(m2·s),只有EM模型的擬合值與其最接近,其他幾個值或遠大于實測值或遠小于實測值,都不理想。光補償點的實測值約是100mol/(m2·s),SFM模型的擬合值是負值,這顯然不符合生物學常識,RHM和NHM模型的擬合值都是0.47mol/(m2·s),遠低于實測值,EM模型的擬合值比QMP模型的擬合值更接近真實值。

表1 馬尾松7月份各光響應(yīng)參數(shù)的測量值及模型擬合值的擬合優(yōu)度比較

表2 馬尾松9月份各光響應(yīng)參數(shù)測量值及模型擬合值的擬合優(yōu)度比較

再分析比較擬合優(yōu)度發(fā)現(xiàn),雖然QMP模型的R2值不是最高的,但其RMSE、MAE和AIC值卻都是最低的,表明該模型具有一定的擬合優(yōu)勢。EM模型的擬合優(yōu)度值雖均略高于QMP模型,但其擬合值更接近實際值,所以選擇EM模型作為9月份最適用的光響應(yīng)模型。

2.3馬尾松11月份的光響應(yīng)曲線及模型擬合比較

11月份馬尾松的實測凈光合速率對光合有效輻射強度的響應(yīng)曲線如圖3所示。在給定的光照強度范圍內(nèi),凈光合速率隨光照強度的上升一直緩慢升高直至達到光飽和狀態(tài)。弱光階段不像7月份和9月份時Pn隨著光強的增大而迅速升高,沒有呈現(xiàn)出明顯的線性關(guān)系。利用6種光響應(yīng)模型進行擬合,擬合曲線如圖3所示。從圖3中可以看出NHM和MRHM模型的擬合值與實測值(MV)偏差非常大,RHM模型在由弱光向強光的上升階段的擬合值比實際值略有出入,EM、QMP和SFM模型的擬合值與實際值都比較接近。

從表3中可以看出,同9月份一樣,MRHM模型在對11月份測定的光響應(yīng)數(shù)據(jù)進行擬合時,由于β=0,使得修正的直角雙曲線模型變成直角雙曲線模型,導(dǎo)致不能再利用公式直接求出Pnmax、LSP和LCP的值。6種模型中只有NRHM模型擬合的表觀量子效率α值不理想,同7月份擬合的值一樣,達到了理論上限值 0.125,也是不符合生理學意義的。通過擬合優(yōu)度比較也可以看出,雖然NRHM的決定系數(shù)R2值較高,但是RMSE、MAE和AIC值也很高,而且其Pnmax、LSP和LCP的擬合值與實測值又相差很大,明顯沒有擬合優(yōu)勢。排除MRHM和NRHM模型,其它幾個模型擬合的α值都比較理想,Pnmax擬合值與實測值也都比較接近。光飽和點的實測值約是1800mol/(m2·s),RHM模型擬合的光飽和點遠低于實測值,QMP和EM的擬合值與實測值都比較接近,但是SFM模型擬合的LCP值與實測值相差比較大,EM模型擬合的LCP值與實測值最接近。

通過分析比較各模型的模擬值與實測值的接近程度,可以發(fā)現(xiàn)QMP模型和EM模型的擬合值與實測值都比較接近,但是再比較這兩個模型的擬合優(yōu)度可以發(fā)現(xiàn),QMP模型的決定系數(shù)R2值不但略大于EM模型,而且其他幾個擬合優(yōu)度值都低于EM模型,具有更好的擬合優(yōu)勢,所以選擇QMP模型作為11月份最適合的光響應(yīng)曲線模型。

表3 馬尾松11月份各光響應(yīng)參數(shù)測量值及模型擬合值的擬合優(yōu)度比較

圖3 11月份不同模型擬合的光響應(yīng)曲線

3 討論

光響應(yīng)曲線模型是光合作用對環(huán)境因子響應(yīng)模擬系統(tǒng)的關(guān)鍵,通過構(gòu)建理想的光合光響應(yīng)模型,準確計算出光飽和點、光補償點、表觀量子效率、暗呼吸速率和最大凈光合速率等指標,可以更客觀準確地反應(yīng)光響應(yīng)過程。自然界中不同的植物對光照強度的適應(yīng)性不同,耐陰植物和喜光植物必然有不同的光響應(yīng)曲線[20]。本研究選用6種模型來擬合馬尾松不同生長期的光響應(yīng)曲線,但是決定系數(shù)R2的差異性并不明顯。雖然R2值的大小可以作為選擇光響應(yīng)模型的重要指標,但不能說明模型所擬合的光合參數(shù)是否符合實際情況[21]。因此本研究引入了均方根誤差(RMSE)、平均絕對誤差(MAE)和赤池量準則(AIC)值3項評價指標,與R2相結(jié)合,同時比較擬合值與實測值的接近程度,來確定馬尾松不同生長期擬合效果最佳的光響應(yīng)曲線模型。

二次多項式回歸模型是從統(tǒng)計學角度出發(fā),來分析凈光合速率(Pn)與光合有效輻射(PAR)之間的關(guān)系。雖然有些研究[22-23]認為二項式回歸模型擬合的光合速率在超過光飽和點后下降不符合生理學意義,存在著局限型,但本研究中該模型擬合的馬尾松11月份的光響應(yīng)參數(shù)卻與實際值最接近,尤其是擬合的光飽和點幾乎與實測值一致,擬合優(yōu)度最高；其他研究也有類似的結(jié)論[24-26]。

分段函數(shù)的特點是將弱光階段和強光階段分開,分別擬合弱光條件下凈光合速率與光強的線性關(guān)系和強光條件下光響應(yīng)曲線的拋物線特點。但是利用分段函數(shù)進行擬合的一個關(guān)鍵點是如何界定弱光條件,段愛國等的研究[8,19]發(fā)現(xiàn)弱光區(qū)域的界定對光合光響應(yīng)曲線的特征參數(shù)能產(chǎn)生顯著的影響。張利陽等的研究[27]認為,對光抑制顯現(xiàn)不是很明顯的毛竹,采用分段函數(shù)擬合可以準確地反映其光響應(yīng)實質(zhì),并將毛竹的弱光區(qū)域界定為50～100mol/(m2·s)。本研究根據(jù)凈光合速率隨光強上升的特點,將弱光強度定義為PAR≤150mol/(m2·s),擬合效果較好,在擬合11月份馬尾松的光響應(yīng)參數(shù)時,其決定系數(shù)高達0.99,但其擬合的光補償點卻低于實測值。

直角雙曲線模型、非直角雙曲線模型和直角雙曲線修正模型是目前應(yīng)用比較廣泛的光響應(yīng)曲線模型。直接雙曲線和非直角雙曲線都是沒有極點的漸近線,所以擬合的光合速率隨光強呈現(xiàn)不斷上升的趨勢,而且無法直接求出最大凈光合速率和光飽和點的解析解,需要借助其他方法來計算[29-30]。非直角雙曲線的擬合考慮了曲線的凸度,拐點更為明顯,擬合得到的光合速率逐漸趨于平緩,使得光飽和的特征更為明顯,因此非直角雙曲線模型的擬合結(jié)果要好于直角雙曲線的擬合結(jié)果[1,19,28]。研究發(fā)現(xiàn)非直角雙曲線模型對落葉松[31]、深山含笑[32]、栓皮櫟[33]、油茶[34]等植物的光響應(yīng)曲線擬合效果較好。直角雙曲線修正模型是葉子飄[35]在直角雙曲線模型的基礎(chǔ)上建立起來的,其擬合的曲線存在極值,能夠直接求解最大凈光合速率和光飽和點。光合速率在達到飽和光強后呈現(xiàn)非線性下降的趨勢,較好地表現(xiàn)出植物光合作用的光抑制現(xiàn)象,具有較廣泛的適用性。在荷花[36]、甘薯[37]、黃枝油杉[38]、野生小蓬竹[39]、歐美楊[40]等植物光響應(yīng)曲線擬合的研究中得到了證實。本研究的結(jié)果也表明,馬尾松7月份的光響應(yīng)參數(shù)利用直角雙曲線修正模型擬合效果最佳。

指數(shù)函數(shù)模型雖然應(yīng)用得不是很廣泛,但有研究者認為[41]其在弱光環(huán)境下對光響應(yīng)過程的LCP有較好的擬合效果,肖丹丹等的研究[42]認為指數(shù)模型是擬合榆屬植物光響應(yīng)曲線的最優(yōu)模型。本研究的結(jié)果也表明,馬尾松9月份光響應(yīng)曲線的最佳擬合模型是指數(shù)模型。

4 結(jié)論

馬尾松在不同的生長期光合特性不同,需要選用不同的光響應(yīng)曲線模型來準確擬合光響應(yīng)參數(shù)。

7月份最適合的馬尾松光響應(yīng)曲線模型是直角雙曲線修正模型(MRHM),其擬合的初始量子效率(α)為0.066mol/(m2·s),最大凈光合速率(Pnmax)為5.13mol/(m2·s),暗呼吸速率(Rd)為3.40mol/(m2·s),光飽和點(LSP)為1148.24mol/(m2·s),光補償點(LCP)為51.47mol/(m2·s)；9月份最適合的馬尾松光響應(yīng)曲線模型為指數(shù)函數(shù)模型(EM),其擬合的初始量子效率(α)為0.003mol/(m2·s),最大凈光合速率(Pnmax)為2.58mol/(m2·s),暗呼吸速率(Rd)為0.25mol/(m2·s),光飽和點(LSP)為1868.67mol/(m2·s),光補償點(LCP)為83.33mol/(m2·s)；11月份最適合的馬尾松光響應(yīng)曲線模型是二次多項式回歸(SFM),其擬合的初始量子效率(α)為0.007mol/(m2·s),最大凈光合速率(Pnmax)為1.38mol/(m2·s),暗呼吸速率(Rd)為0.54mol/(m2·s),光飽和點(LSP)為1800mol/(m2·s),光補償點(LCP)為138.98mol/(m2·s)。導(dǎo)致這種差異性的原因可能是與不同生長期的水熱條件有關(guān),這種相關(guān)性還有待于進一步深入探討。

本研究結(jié)果是根據(jù)測定同一種源不同生長期的光響應(yīng)特性而來的,具有一定的局限型。至于不同種源在不同時期的光響應(yīng)特性是否有顯著差異性還有待于進一步研究。

[1] 葉子飄,于強.光合作用光響應(yīng)模型的比較[J].植物生態(tài)學報,2008,32(6):1356-1361.

[2] 王榮榮,夏江寶,楊吉華,等.貝殼砂生境干旱脅迫下杠柳葉片光合光響應(yīng)模型比較[J].植物生態(tài)學報,2013,37(2):111-121.

[3] Fang L D, Zhang S Y, Zhang G C, et al. Application of five light-response model in the photosynthesis ofPopulusEuramericanacv. ‘Zhongline46’leaves[J]. Applied Biochemistry Biotechnology, 2015, 176(1): 86-100.

[4] 朱世忠.油松光合光響應(yīng)曲線的模型擬合比較研究[J].山西大學學報:自然科學版,2016,39(4):679-685.

[5] 秦國峰,周至春.中國馬尾松優(yōu)良種質(zhì)資源[M].北京:中國林業(yè)出版社,2012.

[6] 魏永成,劉青華,周志春,等.馬尾松不同產(chǎn)脂力無性系光合特性差異分析[J].東北林業(yè)大學學報,2015,43(2):36-41.

[7] 張向峰,王玉杰,李云霞,等.模擬酸雨對盆栽馬尾松葉片光合特性的影響[J].水土保持學報,2013,27(1):197-202.

[8] 段愛國,張建國,何彩云,等.干熱河谷主要植被恢復(fù)樹種干季光合光響應(yīng)生理參數(shù)[J].林業(yè)科學,2010,46(3):68-73.

[9] Lewis J D, Olszyk D, Tingeyd T. Seasonal patterns of photosynthetic light response inDoughlas-firseedlings subjected to elevated atmospheric CO2and temperature[J]. Tree Physiol, 1999, 19(4/5): 243-252.

[10] Thornley J H M．Dynamic model of leaf photosynthesis with acclimation to light and nitrogen[J]. Annals of Botany, 1998, 81: 421-430.

[11] Ye Z P. A new model for relationship between light intensity and the rate of photosynthesis inOryzasativa[J]. Photosynthetica, 2007, 45(4): 637-640.

[12] Bassman J, Zwier J C. Gas exchange characteristics ofPopulustrichocarpa,PopulusdeltoidesandPopulustrichocarpa×P.deltiodesclone[J]. Tree Physiology, 1991, 8: 145-159.

[13] 葉子飄.光合作用對光合CO2響應(yīng)模型的研究進展[J].植物生態(tài)學報,2014,34(6):727-740.

[14] 黃紅英,竇新永,孫蓓育,等.兩種不同生態(tài)型麻瘋樹夏季光合特性的比較[J].生態(tài)學報,2009,29(6):2861-2867.

[15] Chen Z Y, Peng Z S, Yang J, et al. A mathematical model for describing light-response curves inNicotianatabacumL[J]. Photosynthetica, 2011, 49: 467-471.

[16] Akaike H. A new look at the statistical model identification[J]. IEEE Transactions on Automatic Control, 1974, 19(6): 716-723.

[17] 高利霞,畢潤成,閆明.山西霍山油松林的物種多度分布格局[J].植物生態(tài)學報,2011,35(12):1256-1270.

[18] 劉澤彬,程瑞梅,肖文發(fā),等.不同淹水時間下中華蚊母樹光響應(yīng)特征及其模型比較[J].應(yīng)用生態(tài)學報,2015,26(4):1083-1090.

[19] 段愛國,張建國.光合作用光響應(yīng)曲線模型選擇及低光強屬性界定[J].林業(yè)科學研究,2009,22(6):765-771.

[20] 王雁,馬武昌.扶芳藤、紫藤等7種藤本植物光能利用特性及耐蔭性比較研究[J].林業(yè)科學研究,2007,17(3):305-309.

[21] 錢一凡,廖詠梅,權(quán)秋梅,等.4種光響應(yīng)曲線模型對3種十大功勞屬植物的實用性[J].植物研究,2014(5):716-720.

[22] 劉宇鋒,蕭浪濤,童建華,等.非直角雙曲線模型在光合光響應(yīng)曲線數(shù)據(jù)分析中的應(yīng)用[J].中國農(nóng)學通報,2005,21(8):76-79.

[23] 梁文斌,聶東伶,吳思政,等.短梗大參光合作用光響應(yīng)曲線及模型擬合[J].經(jīng)濟林研究,2014,32(4):38-44.

[24] 劉國華,王福升,丁雨龍,等.4種地被竹光合作用日變化及光響應(yīng)曲線[J].福建林學院學報,2009,29(3):258-263.

[25] 高瓊,王維有,梁冬,等.8個種源油松生長和光合生理指標的比較研究[J].北京林業(yè)大學學報,2014,36(2):97-93.

[26] 郭水良,方芳,黃華,等.外來入侵植物北美車前繁殖及光合生理生態(tài)學研究[J].植物生態(tài)學報,2004,28(6):787-793.

[27] 張利陽,溫國勝,王圣杰,等.毛竹光響應(yīng)模型適用性分析[J].浙江農(nóng)林大學學報,2011,28(2):187-193.

[28] 葉子飄,李進省.光合作用對光響應(yīng)的直角雙曲線修正模型和非直角雙曲線模型的對比研究[J].井岡山大學學報:自然科學版,2010,31(3):38-44.

[29] Walker D A, Jarvis P G, Farquhar, et al. Automated measurement of leaf photosynthetic O2evolution as a function of photon flux density[J]. Philosophical Transactions of the Royal Society London B, 1989, 323(1216): 313-326.

[30] Richardson A D, Berlyn G P. Changes in foliar spectral reflectance and chlorophyll fluorescence of four temperate species following branch cutting[J]. Tree Physiology, 2002, 22(7): 499-506.

[31] 許晨璐,孫曉梅,張守攻.本落葉松與長白落葉松及其雜種光合特性比較[J].北京林業(yè)大學學報,2012,34(4):62-66.

[32] 盤李均,冼桿標,陳偉光,等.深山含笑等18個景觀樹種光合生理特性研究[J].廣東林業(yè)科技,2013,29(4):33-37.

[33] 楊自立,馬履一,賈忠奎,等.不同供氮水平對栓皮櫟播種苗光響應(yīng)曲線的影響[J].北京林業(yè)大學學報,2011,33(5):56-60.

[34] 孔文娟,劉學錄,姚小華,等.4個油茶物種的光合特性研究[J].西南大學學報：自然科學版,2013,35(1):16-22.

[35] 葉子飄,于強.一個光合作用光響應(yīng)新模型與傳統(tǒng)模型的比較[J].沈陽大學學報,2007,38(6):771-775.

[36] 冷寒冰,秦俊,葉康,等.不同光照環(huán)境下荷花光合光響應(yīng)曲線模型比較[J].應(yīng)用生態(tài)學報,2014,25(10):2855-2800.

[37] 張磊,劉維正,辛國勝,等.3種專用型甘薯光合光響應(yīng)曲線及其模型擬合研究[J].中國農(nóng)學通報,2015,31(15):71-77.

[38] 柴勝豐,唐健民,楊雪,等.4種模型對黃枝由杉光合光響應(yīng)曲線的擬合分析[J].廣西科學院學報,2015,31(4):1-4.

[39] 廖小鋒,劉濟明,張東凱,等.野生小蓬竹的光合光響應(yīng)曲線及其模型擬合[J].中南林業(yè)科技大學學報,2012,32(3):124-128.

[40] 馬永春,方升佐.歐美楊107不同冠層光合特性的研究[J].南京林業(yè)大學學報:自然科學版,2011,35(4):39-42.

[41] 劉強,李鳳日,謝龍飛.人工長白落葉松冠層光合作用-光響應(yīng)曲線最優(yōu)模型[J].應(yīng)用生態(tài)學報,2016,27(8):2420-2428.

[42] 肖丹丹,左力輝,王茂進,等.5種榆屬植物光響應(yīng)曲線模型對比分析[J].中國農(nóng)業(yè)科技導(dǎo)報,2016,18(4):124-131.

(責任編輯：許晶晶)

ModelFittingofPhotosyntheticPhoto-responseCurveofMassonpine

CHEN Xue-lian, XU Liu-yi*, HAO Yan-ping, ZHANG Jian, JIANG Chun-wu, PAN Ting

(Forestry Research Institute of Anhui Province, Hefei 230031, China)

S791.248

A

1001-8581(2017)10-0024-07

2017-07-14

林業(yè)公益性行業(yè)科研專項項目“抗松材線蟲病馬尾松種源及抗性育種技術(shù)研究”(201004072)。

陳雪蓮(1984─),女,安徽鳳臺人,助理研究員,碩士,研究方向：林木遺傳育種。*通訊作者：徐六一。