許 丹 伍維模，2* 王家強(qiáng) 李志軍，2 武建林

(1 塔里木大學(xué)植物科學(xué)學(xué)院，新疆 阿拉爾843300)(2 新疆生產(chǎn)建設(shè)兵團(tuán)塔里木盆地生物資源保護(hù)與利用重點(diǎn)實(shí)驗(yàn)室，新疆 阿拉爾 843300)

胡楊(Populus euphratica)是楊柳科(Salicaceae)楊屬落葉喬木，典型的潛水旱中生至中生植物，為我國首批確定的388 種珍稀瀕危植物中的漸危種，是新疆維吾爾自治區(qū)人民政府列為優(yōu)先保護(hù)的瀕危物種之一［1－2］。胡楊對荒漠生境具有高度的適應(yīng)性，是重要的荒漠河岸林樹種，是研究古植物學(xué)和植物對氣候變化響應(yīng)的極有價值的材料，也是改良楊柳科其他樹種的優(yōu)良抗逆基因資源。全世界60%的胡楊天然林集中分布在中國，而中國91.1%胡楊林集中分布在新疆塔里木河流域［1］。胡楊在維持荒漠區(qū)生態(tài)平衡和南疆區(qū)域經(jīng)濟(jì)發(fā)展中起著舉足輕重的作用［3］，它具有重要的生物學(xué)與生態(tài)學(xué)保護(hù)利用價值。

葉綠素是廣泛存在于綠色植物中的最主要色素，是光合作用的捕光物質(zhì)，在光合作用中發(fā)揮著重要的生理功能［4］。森林冠層葉綠素含量直接控制著森林能量和物質(zhì)循環(huán)，直接反映著森林的健康和脅迫情況［5］。植物在受到病蟲害威脅或外界環(huán)境脅迫時，葉片葉綠素含量、光合作用等會發(fā)生變化。葉綠素含量的變化提供了植物與環(huán)境關(guān)系的重要信息［6］。利用植物的光譜反射率與葉綠素含量間的相關(guān)性可以進(jìn)行葉綠素含量的估測［7，8］，這種非破壞性、可重復(fù)地在一定空間尺度上連續(xù)地監(jiān)測葉綠素的技術(shù)與傳統(tǒng)的化學(xué)分析方法相比較，具有快速、動態(tài)和實(shí)時的優(yōu)點(diǎn)，已經(jīng)在植物葉綠素的動態(tài)監(jiān)測上得到了廣泛地應(yīng)用［9－11］。如黑麥草冠層葉綠素a、葉綠素b、葉綠素a+b 與紅邊位置之間均存在顯著相關(guān)性［11］。石韌等(2008)測定了健康落葉松與遭受病害落葉松的冠層光譜反射率，根據(jù)光譜反射率與落葉松冠層光合色素含量的相關(guān)性，建立了診斷落葉松健康狀況的遙感方法［12］。

當(dāng)前，對塔里木河天然胡楊葉片葉綠素含量及其與光譜反射率關(guān)系的研究并不多。雖然，王海珍和韓路(2008)測定了阿瓦提縣天然胡楊葉片的葉綠素含量［13］，但其結(jié)果僅有7月份的，缺乏對胡楊整個生長季節(jié)葉綠素含量的變化的研究。袁月等(2009)也僅測定了艾比湖天然保護(hù)區(qū)胡楊5月份的葉綠素含量［14］。從2000年開始實(shí)施塔里木河下游應(yīng)急輸水工程以后，較多的研究人員關(guān)注塔里木河下游胡楊葉綠素與土壤水分和應(yīng)急輸水后的響應(yīng)。在塔里木河下游輸水河道，距離河道越遠(yuǎn)(從100 米至500 米)或者隨著地下水埋深的加深，水分脅迫程度逐漸加重，胡楊葉片葉綠素含量也越低［15－18］。

雖然遙感技術(shù)已經(jīng)在植物葉綠素的動態(tài)監(jiān)測上得到了廣泛的應(yīng)用，但是，關(guān)于塔里木河流域上游胡楊葉片的光譜反射率與葉綠素含量之關(guān)系的基礎(chǔ)研究還十分地缺乏。因此，本研究的目的是通過實(shí)地測量塔里木河上游天然胡楊葉片的葉綠素含量和可見光－近紅外光譜反射率，分析葉綠素含量的季節(jié)變化和不同波長光譜反射率與葉綠素含量的相關(guān)性，從而為利用遙感技術(shù)來監(jiān)測胡楊葉綠素的變化，評價天然胡楊的健康狀況和環(huán)境脅迫提供基礎(chǔ)理論依據(jù)。

1 研究材料與方法

1.1 研究區(qū)概況

研究區(qū)位于塔里木河流域上游源流區(qū)葉爾羌河下游，全年干旱少雨，蒸發(fā)強(qiáng)烈，日照充足，熱量豐富，降水稀少，溫差大，適應(yīng)多種植物生長發(fā)育。年日照時數(shù)平均2 915.1 小時，年平均氣溫10.8℃，7月 份 平 均 氣 溫25.3℃，1月 份 平 均 氣 溫－8.3℃，年際溫差33.6℃。全年平均無霜期204 天，全年平均降水量42.7 毫米，年蒸發(fā)量2 110.5毫米，相差49.4 倍。土壤主要是河相沉積而成，以草甸為主，有較高的肥力。研究區(qū)樣地位于新疆阿瓦提縣胡楊林保護(hù)區(qū)艾買卻克墩鎮(zhèn)護(hù)林站天然胡楊林，北緯40°17'、東經(jīng)80°21'，是塔里木河三大源游之一葉爾羌河的下游。胡楊林下地表主要有三種類型分別是鹽結(jié)殼、風(fēng)沙土和自然溝。林下植被主要是甘草、檉柳等。先對艾買卻克墩鎮(zhèn)護(hù)林站天然胡楊林區(qū)進(jìn)行二天的野外踏察，結(jié)合衛(wèi)星遙感影像，建立了30 塊監(jiān)測樣地(圖1)，用GPS 接收機(jī)測定了樣地的經(jīng)度、緯度。樣地大小設(shè)為50 m ×50 m。在每個樣地中選擇一棵中心樹，用紅油漆進(jìn)行標(biāo)記。

1.2 胡楊光譜數(shù)據(jù)的采集

采用合肥儀思特光電技術(shù)有限公司生產(chǎn)的ISI921VF－512 野外地物光譜輻射計測量胡楊葉片的可見光－近紅外反射光譜。測量的波譜范圍是393～1 095 nm，共512 個波段，光譜采樣間隔為1.4 nm，光譜分辨率為3.5 nm，視場角為2.50。為減小太陽高度角變化對光譜測量結(jié)果的影響，測量時間選擇在12:00～16:00。儀器探頭垂直向下，先測量標(biāo)準(zhǔn)板，然后測定胡楊葉片。葉片放在有陽光直接照射的林間空地的黑布上。測量日期分別是2011年5月22日、6月17日、7月12日、8月3日、9月25日、和10月15。測量當(dāng)天是晴朗無云無風(fēng)的天氣。每個月份均測量了90 個樣本。每個樣地測量3 個樣本。每個樣本是從一棵胡楊中部冠層的枝條上隨機(jī)選擇的10 個葉片，以10 個葉片的光譜反射率的平均值作為樣本值。

圖1 天然胡楊林研究區(qū)樣地分布圖

1.3 葉綠素含量測定

按照Sims and Gamon(2002)的方法［9］，用丙酮:Tris 緩沖液(體積比為80:20)浸提胡楊葉片中的葉綠素。將新鮮葉片剪成長寬約2 mm的碎片。準(zhǔn)確稱取0.200 0 克，放入具塞刻度試管中，然后加入20 ml 丙酮:Tris 緩沖液，再放入冰柜中于4℃下避光、密封浸提48 小時以上。待葉片全部變成白色后進(jìn)行比色測定537 nm，647 nm 和663 nm 波長下浸提液的吸光度。按下列公式計算葉綠素a、葉綠素b的濃度。

葉綠素a(chla)=0.01373A663－0.000897A537－0.003046A647

葉綠素b(chlb)=0.02405A647－0.004305A537－0.005507A663

公式中的A663表示在波長663 nm 時用1 cm的比色杯測量的浸提液的吸光度值。chla，chlb 分別表示浸提液中葉綠素a 和葉綠素b的摩爾濃度，單位是ū molml－1。根據(jù)葉綠素a、b的分子量(chla=893.5gmol－1，chlb=907.5gmol－1)和胡楊葉片的含水率、比葉重計算單位葉片鮮重的葉綠素含量(mg/g F w)、單位葉片干重的葉綠素含量(mg/g DW)，和單位葉片葉面積的葉綠素含量(mg/m2)。

1.4 葉片比葉重(cm2/gFw)和含水量的測定

將葉片用萬分之一天平準(zhǔn)確稱取重量(0.0001 g)后，在白紙上用鉛筆把葉片的形狀描繪下來。然后用哈爾濱光學(xué)儀器廠生產(chǎn)的QCJ—2A型數(shù)學(xué)式求積儀測量葉片圖形的面積。葉面積除以該葉片的鮮重即為比葉重。

先用萬分之一天平準(zhǔn)確稱取10 個胡楊新鮮葉片的鮮重，計為W0;然后放入鋁盒中，于80℃下烘干24 小時后稱葉片的干重，計為W1;葉片含水率φ以干基為準(zhǔn)，用下列公式計算。

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

用方差分析來檢驗(yàn)月份間胡楊葉片的葉綠素a、b 和葉綠素a +b 是否存在顯著差異。按照完全隨機(jī)設(shè)計試驗(yàn)方案，以6 個月份(5月、6月、7月、8月、9月、10月)為處理，30 個樣地作為重復(fù)。分別以葉綠素a、b 和葉綠素a+b 含量為自變量，各波長反射率為因變量，對葉綠素與反射率進(jìn)行相關(guān)分析。方差分析與相關(guān)分析使用的是DPS 11.5 軟件［21］。

2 結(jié)果與分析

2.1 胡楊葉片的葉綠素含量

圖2 2011年5～10月胡楊葉片葉綠素含量的比較

2.2 葉綠素色素含量與反射率的相關(guān)性

從圖2 及方差分析可知，在胡楊葉片的生長季節(jié)(5月～10月)內(nèi)，10月葉綠素a、葉綠素b，和葉綠素a+b的含量極顯著地低于其他月份。10月葉綠素a、葉綠素b，和葉綠素a+b 分別為47、18 和65 mg/m2，而8月的葉綠素a、葉綠素b，和葉綠素a +b 分別為277、86 和365 mg/m2。5月與8月相比，胡楊葉片的葉綠素a、葉綠素b，和葉綠素a+b的含量無顯著差異，8月略高，但未達(dá)到顯著水平。6月、7月和9月胡楊葉片的葉綠素a、葉綠素b，和葉綠素a+b的含量處于中間水平。從5月～10月，胡楊葉片的葉綠素a 明顯地高于葉綠素b的含量。用三種單位表示的葉綠素含量(mg/m2，mg/g F w，和mg/g DW)在不同月份之間的變化規(guī)律是相同的，均表現(xiàn)出10月份最低，8月份最高。5月份胡楊葉片的葉綠素含量要高于6月份，從而出現(xiàn)了隨月份變化，葉綠素呈現(xiàn)出“兩峰兩谷”型的現(xiàn)象，即:5月份葉綠素含量是一個“峰”，8月份是另外一個“峰”，6月出現(xiàn)了第一個“谷”，10月出現(xiàn)另外一個“谷”。塔里木河上游胡楊葉片葉綠素含量隨季節(jié)變化出現(xiàn)的“兩峰兩谷”現(xiàn)象還未見報道。在塔里木河上游阿瓦提縣境內(nèi)，據(jù)氣象資料的統(tǒng)計，10月下旬開始出現(xiàn)“初霜”，胡楊葉片普遍開始褪綠，變?yōu)榈S色，至金黃色，出現(xiàn)塔里木特有的金色胡楊的美麗景色。從10月開始，氣溫下降，日照時間縮短，胡楊葉片內(nèi)葉綠素降解，含量迅速下降，最后葉片也脫落，表明10月份葉綠素含量最低是符合胡楊的生長規(guī)律的。

圖3 2011年5月至10月胡楊葉片可見光－近紅外波段反射光譜

圖3 所示，胡楊葉片5月～10月的可見光－近紅外反射光譜符合綠色植物反射光譜曲線的特征。綠色植物在450～500 nm 為中心的藍(lán)波段及以650～680 nm 為中心的紅波段，由于光合色素強(qiáng)烈吸收太陽輻射能而形成兩個通常被稱為“藍(lán)谷”和“紅谷”的低反射區(qū)。在這兩個低反射區(qū)(“谷”)之間約550 nm 附近，植物葉片對輻射能吸收很少，因而形成了綠色反射峰即通常所說的“綠峰”。典型綠色植被光譜曲線還有一明顯的“紅邊”特征。植物光譜曲線在紅波段的葉綠素吸收谷與近紅外波段的高反射肩之間的陡升部分，約在700～750 nm 處，被稱為植物的“紅邊”區(qū)。胡楊葉片的可見光－近紅外反射光譜也具有“藍(lán)谷”，“紅谷”，“綠峰”和“紅邊”特征，其最高反射率在近紅外波段，約35～50%，最低反射率在藍(lán)波段，約為5～12%。在420 nm 和660 nm 附近各有一個反射率極小值(“藍(lán)谷”和“紅谷”);在550 nm 附近處有一個反射峰(“綠峰”);在650 nm～800 nm 之間反射率急劇增大(“紅邊”)，800 nm～950 nm 有高的反射率。10月份在吸收谷670 nm 左右的反射率高于其它5 個月，在750 nm 后明顯低于其它5 個月;6月、7月的光譜反射率在950 nm 波長后低于5月、8月和9月。不同月份間的光譜反射率在近紅外波段差異最大，藍(lán)、綠波段差異最小。6月至9月的反射光譜曲線具有較大的相似性，但是，10月份有反射光譜曲線與其它5 個月份的有比較明顯的差異，表現(xiàn)在近紅外波段和"紅邊"反射率顯著下降，低于43%，“紅谷”反射率上升，“綠峰”反射率下降。

不同波長位置處的不同月份葉綠素與反射率之間的相關(guān)系數(shù)表現(xiàn)出顯著的差異(圖4)。在紅光區(qū)，胡楊葉片葉綠素與反射率呈負(fù)相關(guān)關(guān)系，這可能是由于葉綠素吸收紅色光所致。而在紅邊區(qū)及近紅外區(qū)和部分藍(lán)光區(qū)，葉綠素含量與反射率呈正相關(guān)關(guān)系。經(jīng)統(tǒng)計檢驗(yàn)，最大相關(guān)系數(shù)和最小相關(guān)系數(shù)(負(fù)值)均達(dá)到了極顯著的水平(p＜0.01)。由圖4可知，在769 nm 處葉綠素a(mg/m2)、葉綠素b(mg/m2)、葉綠素a+b(mg/m2)含量與反射率呈極顯著正相關(guān)關(guān)系，相關(guān)系數(shù)最大(r 在0.34～0.41之間)。在645 nm 處葉綠素a (mg/g Dw)、葉綠素b (mg/g Dw)、葉綠素a+b (mg/g Dw)含量與反射率呈極顯著的負(fù)相關(guān)關(guān)系，相關(guān)系數(shù)最小(r 在－0.45～－0.35 之間)。

圖4 不同波長位置處胡楊光譜反射率與葉綠素的相關(guān)系數(shù)注:樣本容量為507，顯著水平為0.01 時的r 臨界值為0.12。

3 討論與結(jié)論

3.1 塔里木河上游胡楊葉片的葉綠素含量在生長季節(jié)內(nèi)不同月份之間存在顯著的差異，呈現(xiàn)出“兩峰兩谷”的現(xiàn)象，其中8月份葉綠素含量最高，10月份最低。

本研究測定了塔里木河上游天然胡楊葉片生長季節(jié)內(nèi)(5月至10月)的葉綠素含量，發(fā)現(xiàn)了葉綠素含量隨月份變化呈現(xiàn)出"兩峰兩谷"的現(xiàn)象。據(jù)王海珍和韓路測定胡楊葉片7月葉綠素含量為葉綠素a，葉綠素b 及葉綠素a +b 分別為1.39，0.34，和1.74 mg/g Fw。而本研究測得的葉綠素a，葉綠素b 及葉綠素a+b 分別為1.03，0.34 和1.37 mg/g Fw，比王海珍和韓路(2008)的葉綠素a 測定結(jié)果略低，而葉綠素b的結(jié)果相同。袁月等(2009)對新疆艾比湖自然保護(hù)區(qū)5月份胡楊葉片的葉綠素測定結(jié)果為葉綠素a 0.6～1.4 mg/g Dw，葉綠素b 0.2～0.4 mg/g Dw;而本研究結(jié)果為葉綠素a 2.26 mg/g Dw，葉綠素b 0.71 mg/g Dw，明顯地要高于袁月等(2009)的研究結(jié)果。

3.2 胡楊葉片的可見光－近紅外反射光譜符合綠色植物的反射光譜曲線特征。塔里木河上游胡楊葉片的反射光譜表現(xiàn)出綠色植物特有的反射光譜曲線特征，如“紅邊”，“紅谷”，“綠峰”等，且本研究的結(jié)果與其它研究人員測定的塔里河中下游的胡楊葉片反射光譜曲線結(jié)果是一致的［19，20］。進(jìn)入10月份之后，由于氣溫下降，胡楊葉片開始衰老，葉綠素含量明顯下降，反射光譜中的“紅邊”與近紅外波段反射率顯著下降，“紅谷”反射率抬升。

3.3 胡楊葉片葉綠素含量與可見光－近紅外光譜反射率之間的相關(guān)性因波長不同而不同。隨著波長的變化，不同波長位置處的相關(guān)系數(shù)呈現(xiàn)出“V”型變化。769 nm 和645 nm 是兩個具有特殊意義的波長位置。

植物葉片葉綠素含量的估測可以利用高光譜遙感數(shù)據(jù)與葉綠素的相關(guān)性來實(shí)現(xiàn)。如楊曦光將Hyperion 數(shù)據(jù)的像元反射率轉(zhuǎn)化為葉片反射率，從影像上反演出葉片水平的葉綠素含量［5］。通過本研究，作者認(rèn)為基于單個波長位置的反射率與葉綠素相關(guān)性來估測胡楊葉片葉綠素含量將會存在一定的局限性.因?yàn)?，單個波長處反射率與葉綠素的相關(guān)系數(shù)絕對值最大不超過0.5，同時，本研究也發(fā)現(xiàn)胡楊葉片在769 nm 處的葉綠素含量與反射率呈極顯著正相關(guān)關(guān)系，而在645 nm 處呈極顯著的負(fù)相關(guān)關(guān)系。因此，結(jié)合多個波長位置的反射率數(shù)據(jù)或者用高光譜指數(shù)(基于兩個或多個波長位置的反射率的組合指標(biāo))將會比用單一波長位置的反射率來估測葉綠素含量會更有效，這有待于下一步深入地分析研究。

［1］王世績.全球胡楊現(xiàn)狀、保護(hù)與恢復(fù)［J］.世界林業(yè)研究，1996(6):37－43.

［2］蘇培璽，張立新，杜明武，等.胡楊不同葉形光合特性、水分利用效率及其對加富CO2的響應(yīng)［J］.植物生態(tài)學(xué)報，2003，27(1):34－40.

［3］程小玲，劉淑清.淺議塔里木胡楊保護(hù)區(qū)胡楊林的保護(hù)［J］.中南林業(yè)調(diào)查規(guī)劃，2004，23(1):34－36.

［4］盧東昱，崔新圖，黃鏡榮，等.葉綠素吸收光譜的觀測［J］.大學(xué)物理，2006，(1):53－56.

［5］楊曦光，范文義，于穎.基于PROSPECT +SAIL 模型的森林冠層葉綠素含量反演［J］.光譜學(xué)與光譜分析，2010，11，30(11):144－148.

［6］張友勝，張?zhí)K峻，李鎮(zhèn)魁.植物葉綠素特征及其在森林生態(tài)學(xué)研究中的應(yīng)用［J］.安徽農(nóng)業(yè)科學(xué)，2008，36 (3):1014－1017.

［7］童慶禧，張兵，鄭蘭芬.高光譜遙感的多學(xué)科應(yīng)用［M］.北京:電子工業(yè)出版社，2006:196.

［8］楊曦光，范文義，于穎.基于Hyperion 數(shù)據(jù)的森林葉綠素含量反演［J］.東北林業(yè)大學(xué)學(xué)報，2010，6，38(6):127－128.

［9］Daniel A.Sims and John A.Gamon.Relationship between leaf pigment content and spectral reflectance across a wide range of species，leaf structure and development stages.Remote Sensing，2002，81:337－354.

［10］伍維模，牛建龍，溫善菊，等.植物色素高光譜遙感研究進(jìn)展［J］.塔里木大學(xué)學(xué)報，2009，12，(4):67－74.

［11］楊紅麗，陳功，吳建付.氮肥水平對多花黑麥草葉綠素及高光譜反射率特征的影響［J］.云南農(nóng)業(yè)大學(xué)學(xué)報，2009，7，24(4):534－538.

［12］石韌，劉禮，高娜.用高光譜數(shù)據(jù)反演健康與病害落葉松冠層光合色素含量的模型研究－基于2005年吉林省敦化、和龍兩市落葉松冠層采樣測量數(shù)據(jù)［J］.遙感技術(shù)與應(yīng)用，2008，23(3):264－271.

［13］王海珍，韓 路.塔里木荒漠優(yōu)勢樹種的光合特性研究［J］.安徽農(nóng)業(yè)科學(xué)，2008，36(26):11205－11207.

［14］袁月，呂光輝，徐敏，等.干旱脅迫下不同胸徑胡楊生理特點(diǎn)分析［J］.新疆農(nóng)業(yè)科學(xué)，2009，46(2):299－305.

［15］孫衛(wèi).塔里木河下游應(yīng)急輸水后植物響應(yīng)研究［D］.新疆農(nóng)業(yè)大學(xué)碩士學(xué)位論文.新疆農(nóng)業(yè)大學(xué)，2004.

［16］莊麗，陳亞寧，陳明，等.模糊隸屬法在塔里木河荒漠植物抗旱性評論中的應(yīng)用［J］.干旱區(qū)地理，2005，28(3):367－372.

［17］都偉新，張超.塔里木河中下游胡楊葉綠素特征［J］.科技創(chuàng)新導(dǎo)報，2009，23:94.

［18］張青青，徐海量，葉茂，等.不同地下水埋深下胡楊葉片生理指標(biāo)變化特點(diǎn)［J］.植物研究，2009，29(4):453－459.

［19］肖正清，王樹東，古麗·加帕爾，等.塔里木河流域中下游胡楊典型葉片光譜響應(yīng)特征［J］.干旱區(qū)地理，2007，35(5):680－687.

［20］趙金，陳曦，古麗·加帕爾，等.塔里木河荒漠植被光譜可分性模擬［J］.中國沙漠，2009，29(2):270－278.

［21］唐啟義.DPS 數(shù)據(jù)處理系統(tǒng)，實(shí)驗(yàn)設(shè)計、統(tǒng)計分析及數(shù)據(jù)挖掘［M］.北京:科學(xué)出版社，2007.