張 偉，張 宏*，王亁蘊

(1．四川師范大學地理與資源科學學院，四川成都610101; 2．四川師范大學工學院，四川成都610101)

青藏高原東緣高寒草地土壤碳空間異質性

張 偉1，張 宏1*，王亁蘊2

(1．四川師范大學地理與資源科學學院，四川成都610101; 2．四川師范大學工學院，四川成都610101)

運用經典統(tǒng)計學和地統(tǒng)計學的方法，探討青藏高原東緣高寒地區(qū)典型的高寒草甸、高寒灌叢樣地中小尺度下土壤碳的空間變異特征．結果表明，研究區(qū)土壤有機碳變異系數隨著樣地尺度的增加而增大，變化范圍為13．4%～28．4%，均屬中等變異．在高寒草甸和高寒灌叢不同采樣尺度下土壤有機碳均具有良好的空間自相關性;半變異函數C0/(C0+C1)值介于57．93%～80．92%，具有中等程度空間相關性，說明區(qū)域土壤有機碳空間變異的主要原因為人為活動和微域內變異等隨機性因素．在相同尺度下高寒灌叢群落土壤有機碳的分形維數值要小于高寒草甸群落;在研究區(qū)，隨著采樣尺度的增大，植被分布、群落蓋度等隨機因素逐漸減弱，地形起伏等結構因素有逐漸加強趨勢．

土壤碳;地統(tǒng)計學;空間變異;尺度效應

土壤是陸地生態(tài)系統(tǒng)中最大的碳庫，其貯存的有機碳約為大氣碳庫的 2倍，植物碳庫的 3倍［1－3］，是全球碳循環(huán)的重要構成．在全球碳收支研究中精確估算陸地碳庫對二氧化碳的吸收有著至關重要的作用［4－5］．但土壤有機碳庫量的估算在不同學者間存在較大差別，全球土壤碳庫量的估算值為1 500～2 300 Pg［3］，中國土壤有機碳庫量的估算值為50～185 Pg［6－7］，表現出極大的不確定性．由于資料和實驗的限制，尺度選擇和轉換是碳儲量估算差異的重要原因之一．在同一研究區(qū)域、不同空間尺度上的土壤特征變異規(guī)律存在較大差異［8－9］，在認識尺度效應的基礎上，對同一研究區(qū)土壤碳進行多尺度研究，揭示土壤碳在連續(xù)空間上的變異規(guī)律，有助于土壤碳的精確估算．

青藏高原的土壤和植被因其所處高海拔、低溫的地理環(huán)境對氣候變化極其敏感，是全球變化研究的敏感區(qū)，青藏高原高寒草地生態(tài)系統(tǒng)有不同于其他生態(tài)系統(tǒng)獨特的生物地球化學過程［10］，在全球變化研究中占有特殊地位．該區(qū)域常年低溫，極大地減緩了土壤有機質的分解速率，土壤碳儲量巨大，占全球總碳庫的2．4% ～2．6%［11］，約占全國土壤碳儲量的55．6%，是我國土壤碳密度最高的地區(qū)之一［12］．近年來，國內外學者對該區(qū)域的土壤有機碳儲量、分布和轉化及其影響因素極為關注，已開展了大量的研究工作［13－15］，土壤碳空間變異的單一尺度和多尺度研究也均有報道［16］．本文在中小尺度下通過高密度網格嵌套采樣，從不同群落類型、樣地大小、采樣密度等角度探討土壤碳的空間異質性變化，揭示青藏高原東緣高寒地區(qū)不同草地類型土壤碳的空間分布及其變異特征，有助于認識該地區(qū)土壤碳在不同空間范圍上的變異特征，最終為大尺度的土壤碳分布和估算模型模擬和真實信息獲取提供實踐基礎，對評價青藏高原在全球環(huán)境變化下生物地球化學循環(huán)的響應有重要的意義．

1 材料與方法

1．1 研究區(qū)概況 研究區(qū)位于青藏高原東緣，四川省阿壩州紅原縣境內，區(qū)內地勢為東南向西北傾斜，平均海拔3 600 m．氣候屬大陸性高原氣候，氣溫日較差大，全年平均氣溫1．1℃，最熱月平均氣溫10．9℃，最冷月平均氣溫－10．3℃，全年霜凍期長;年均降雨量791．95 mm且主要集中在5～10月，年均蒸發(fā)量1 262．5 mm，降水小于蒸發(fā)，空氣干燥;日照時間長，年均日照時間2 158．7 h，太陽輻射強，太陽輻射年總量為6 194 MJ/m2．當地的主要植被類型為灌叢和草甸，主要的草甸類型有沼澤草甸、亞高山灌叢草甸和帶有高寒性質的亞高山草甸．區(qū)內土壤以沼澤土、高山草甸土和亞高山草甸土為主，其中亞高山草甸土占全縣土地總面積的55．4%［17］．

1．2 樣品采集與分析 樣地設置在紅原縣阿木柯河圍欄禁牧的瑪莫草場．草甸群落樣地(海拔: 3 521 m):群落優(yōu)勢種主要有:四川嵩草(Kobresiasetchwanesis)、絲穎針茅(StipaCapillacea)、西藏嵩草(Kobresiatibetica)、矮嵩草(Kobresiahumilis)等，群落總蓋度在 95%左右．灌叢群落樣地(海拔: 3 522 m):建群種為高山繡線菊(Spiraeaalpina)，伴生有西藏忍冬(Loniceratibetica)、金露梅(Dosiphorafruticosa)和藏沙棘(Hippophaethibetana)等，群落的總蓋度在98%左右．

野外取樣時間處于高寒草地生長旺盛期(7月下旬)．利用網格和大樣地套小樣地的方法在樣地內設點采樣，高寒草甸樣地分別設置2 m×1 m (M0．1)、10 m×5 m(M0．5)和20 m×10 m(M1)3種樣地尺度，在樣地內分別按0．1、0．5、1 m 3種采樣間距取樣，高寒灌叢樣地分別設置10 m×5 m (S0．5)、20 m×10 m(S1)和30 m×15 m(S1．5)3種樣地尺度，在樣地內分別按0．5、1、1．5 m 3種采樣間距進行采樣，取樣深度都為0～20 cm，采集土樣共計1 200個，將土樣編號裝入布袋帶回實驗室，經風干、研磨，過篩，制備裝瓶備用．土壤有機碳含量采用重鉻酸鉀氧化－外加熱法測定，每個土樣3次重復．

1．3 數據處理 本文采用地統(tǒng)計學方法進行分析．在地統(tǒng)計學中，區(qū)域化變量在一定尺度上的相關程度和空間變異可以通過半方差函數的塊金值、基臺值等參數來反映．相關統(tǒng)計學的方法和原理可參考文獻［18－19］．描述性統(tǒng)計利用SPSS13．0完成，自相關分析、分形維數、模型擬合以及空間插值圖利用GS+7．0和ArcGIS9．3完成．

2 結果與分析

2．1 土壤有機碳經典統(tǒng)計分析 從表1可以看出，研究區(qū)高寒草甸土壤有機碳含量為85．42～119．92 g/kg，高寒灌叢土壤有機碳含量為72．09～92．32g/kg，2種群落下土壤有機質含量差異顯著，在不同植被覆蓋，不同采樣尺度下土壤有機碳變異系數在13．4%～28．4%之間，均屬中等變異(0．1＜C．V．＜1)，且變異系數隨樣地尺度的增加而增大．通過對偏度、峰度的觀察和Kolomogorov－Semirnov (K－S)檢驗，在5%的檢驗水平下，高寒灌叢和高寒草甸不同樣地尺度6組土壤有機碳測量數據均服從正態(tài)分布，滿足地統(tǒng)計學分析要求．

表1 青藏高原東緣高寒地區(qū)土壤碳統(tǒng)計特征值Table 1 Statistical characteristics of soil organic carbon of Alpine Area in Eastern Qinghai－Tibetan Plateau

2．2 不同采樣尺度下土壤有機碳空間自相關分析空間自相關是指同一個變量在不同空間位置上的相關性［20］．本文采用Moran’s I自相關指數分析土壤有機碳的空間自相關性隨采樣尺度的變化．研究結果顯示，土壤有機碳在不同取樣尺度下均呈現一定的結構性，高寒草甸和高寒灌叢不同采樣尺度下土壤有機碳的自相關函數變化趨勢基本一致，均隨著滯后距離增大，正相關方向逐漸向負方向變化，但自相關的狀態(tài)和變化程度不同．在高寒草甸土壤有機碳Moran’sI系數波動較大，自相關范圍和自相關系數Moran’sI最大值都隨著樣地尺度的增大而增大，拐點分別出現在0．93、1．4、3．6 m左右，Moran’sI最大值分別為0．158、0．295、0．308;高寒灌叢土壤有機碳Moran’sI系數變化較規(guī)則，自相關范圍隨著樣地尺度的增大而增大，拐點分別出現在2．5、5、11 m左右，但自相關系數Moran’sI最大值卻隨著樣地尺度的增大而減小，Moran’sI最大值分別為0．464、0．438、0．233．由于高寒灌叢群落樣地具有一定的坡度，而高寒草甸群落樣地較前者地勢平坦，在相同采樣尺度下，高寒灌叢土壤有機碳自相關指標值均大于高寒草甸，結構性因素(地形起伏)對研究區(qū)內的高寒灌叢土壤有機碳的空間異質影響明顯要高于高寒草甸群落．

2．3 不同采樣尺度下土壤有機碳的空間結構特征分析 半方差函數參數既能反映出分析對象空間相關的類型又能說明空間分布的結構，同時能夠量化空間相關的范圍，其中，塊金值(C0)表示隨機部分引起的空間異質性，基臺值(C0+C1)表示總的空間異質性，隨機因素引起的變異與總變異的關系通過塊基比(C0/(C0+C1))能夠得到反映，當塊基比較高時，表明引起空間變異的主要作用為隨機部分;當比值較低時，則表明引起空間變異的主導因素為空間自相關部分．在高寒草甸樣地，采樣尺度M0．1和M0．5有機質含量半方差函數最佳擬合模型符合球狀模型，M1符合指數模型．從表2可以看出，隨著采樣尺度的增大半方差函數的變程呈增大趨勢，而塊金值、基臺值和塊基比均隨采樣尺度增大先減小后增大，塊基比分別為64．58%、61．8%、80．92%，M0．1表現出各向同性，而M0．5和M1表現出各向異性，各向異性比分別為4．43和1．73．在高寒灌叢樣地，采樣尺度S0．5和S1．5有機質含量半方差函數最佳擬合模型符合球狀模型，S1符合指數模型．隨著采樣尺度的增大半方差函數的變程呈增大趨勢，塊金值、基臺值也隨著采樣尺度的增大而增大，塊基比則隨采樣尺度的增大而減小，塊基比分別為68．46%、66．33%、57．93%．灌叢樣地各采樣尺度(S0．5、S1和S1．5)均表現出各向異性，各向異性比分別為1．8、1．68和2．26．本研究區(qū)內的高寒灌叢和高寒草甸土壤有機碳在不同采樣尺度下塊基比值均較大，介于57．93%～80．92%，具有中等程度空間相關性，而且表明在該樣地區(qū)域內土壤有機碳的空間變異主要是由微域內變異和人為活動等隨機性因素引起的．

表2 不同樣地尺度的空間變異特征對比分析Table 2 Comparison of spatial varibility under different plot scales

2．4 不同采樣尺度下土壤有機碳的分形維數分析

分形維數是利用半方差函數的雙對數曲線的斜率求得的，是對空間異質性程度的表征參數之一［21］，若空間異質性程度越強，則分形維數的數值越小，反之則增大，當分形維數值趨于2時，表明區(qū)域化變量基本上是空間不相關的，若等于2，可以說整個樣本是空間獨立的［21］．從表2可以看出，高寒草甸群落下土壤有機碳的分形維數值較高，均大于1．9，其中高寒草甸群落采樣尺度M0．1的分形維數值最大，達到1．97．隨著采樣尺度的增大，土壤有機碳的分形維數值在高寒草甸群落下表現為減小趨勢，而在高寒灌叢群落下表現為先減小后增大的變化;在相同尺度下的土壤有機碳的分形維數值，高寒灌叢群落均小于高寒草甸群落，說明高寒灌叢群落土壤有機碳比高寒草甸群落具有更強的空間異質性．

2．5 不同尺度下土壤有機碳的空間插值分析 利用ArcGIS軟件中的地統(tǒng)計模塊對不同采樣尺度下的土壤有機碳分布的Kriging插值結果如圖1，圖2在高寒草甸群落中，采樣尺度M0．1下土壤有機碳分布較均勻，沒有明顯方向趨勢，在M0．5中則表現出明顯的方向性，土壤有機碳呈南北向條帶狀分布，其中的低值區(qū)與樣地內的裸露斑塊的位置和延伸方向一致，在M1中表現出東西向的帶狀分布趨勢．在高寒灌叢群落樣地中，采樣尺度S0．5下，土壤碳含量表現出由灌叢下向外逐漸減少的趨勢，有明顯的“沃島”效應，與其它地區(qū)的研究一致［22］，隨著采樣尺度的增大，插值圖逐漸掩蓋了小尺度下的變異，由于地形的影響有機碳含量呈現出沿坡面向下逐漸增加，在坡底聚集［23］，地形起伏在景觀尺度上控制著土壤養(yǎng)分的分布狀況，是影響土壤空間異質性的一個重要因子［24］．

3 討論

3．1 不同植被群落對土壤有機碳的影響 在本研究區(qū)內高寒草甸表層土壤有機碳含量為85．42～119．92 g/kg，高寒灌叢土壤有機碳含量為72．09～92．32 g/kg，2種群落下土壤有機質含量差異顯著，高寒草甸土壤有機碳平均含量高于高寒灌叢土壤，略低于三江平原濕地土壤有機碳含量［25］，但遠高于亞熱帶濕潤區(qū)濕地和半干旱區(qū)主要草原土壤有機碳含量［26－27］，土壤有機碳的含量和分布與研究區(qū)的群落組成、蓋度、地形、水熱條件和土壤特性有一定關系［28］．高寒灌叢和高寒草甸帶所處海拔高，冷濕的氣候環(huán)境，有機質分解率低，利于有機碳的累積．在草甸群落中，植物地下根系生物量集中于表層，且每年死亡率高，周轉速率快，利于有機碳的積累;而灌叢根系較深，地下生物量周轉速率低，使其有機質和養(yǎng)分富集較慢．

3．2 高寒草甸群落土壤有機碳空間變異特征 高寒草甸群落土壤有機碳在不同取樣尺度下均呈現一定的結構性，自相關范圍和自相關系數Moran’s I最大值都隨著樣地尺度的增大而增大，表明高寒草甸土壤有機碳含量隨著采樣尺度的增大由氣候、地形起伏等結構性因素影響增強;隨著采樣尺度的增大半方差函數的變程呈增大趨勢，而塊金值、基臺值和塊基比均隨采樣尺度增大先減小后增大，M0．1表現出各向同性，而M0．5和M1表現出各向異性，在M0．5中土壤有機碳的方向變異與樣地內裸露斑塊的出現，斑塊大小和分布方向基本一致，在M1中有機碳的主要變異方向與所處的坡向基本一致．分形維數值隨采樣尺度的增大而緩慢減小，其中高寒草甸群落采樣尺度M0．1的分形維數值最大，為1．97說明其空間異質性較弱，趨于均質這與前面的半方差函數分析結果一致．

3．3 高寒灌叢群落土壤有機碳空間變異特征 高寒灌叢土壤有機碳Moran’sI系數變化較規(guī)則，自相關范圍隨著樣地尺度的增大而增大，但自相關系數Moran’sI最大值卻隨著樣地尺度的增大而減小，表明高寒灌叢土壤有機碳隨著采樣尺度的增大人為活動等隨機因素影響逐漸加強．由于高寒灌叢群落樣地具有一定的坡度，而高寒草甸群落樣地較前者地勢平坦，在相同采樣尺度下，高寒灌叢土壤有機碳自相關指標值均大于高寒草甸，結構性因素對高寒灌叢土壤有機碳的空間異質影響明顯要高于高寒草甸群落．在高寒灌叢樣地，隨著采樣尺度的增大半方差函數的變程、塊金值、基臺值均呈增大趨勢，塊基比則隨采樣尺度的增大而減小，空間相關性增強，這是由于隨著采樣尺度的增大灌叢叢間和叢下表現出來的有機碳含量分布差異［29］，逐漸被掩蓋，局部地形和地貌的影響更加突出．由于灌叢分布和局部地形的影響，灌叢樣地各采樣尺度均表現出各向異性，分形維數值隨采樣尺度的增大呈現出先減小后增大的變化趨勢，但在相同尺度下高寒灌叢群落小于高寒草甸群落，表明灌叢群落入侵草甸植被后會使土壤養(yǎng)分空間變異增大［30］．

4 結論

1)研究區(qū)內，土壤有機碳含量較高，平均含量均高于70 g/kg，高寒草甸群落有機碳含量明顯高于高寒灌叢群落．由于采樣尺度的差異，土壤有機碳變異系數變化范圍為13．4%～28．4%，屬中等變異，空間變異隨樣地尺度的增加而增大．

2)高寒草甸群落土壤有機碳含量空間分布的自相關系數和自相關距離隨著采樣尺度的增大而增大，分形維數隨著采樣尺度的增大而減小，塊基比則隨著采樣尺度的增大先減小后增大．高寒灌叢群落土壤有機碳含量空間分布的自相關距離隨著采樣尺度的增大而增大，自相關系數和塊基比則隨著采樣尺度的增大而減小，分形維數的變化趨勢是先減小后增大．

3)在不同群落覆蓋，不同采樣尺度下，除了M0．1，其余各尺度下土壤有機碳空間分布均存在各向異性，隨著采樣尺度的增大，植被分布，群落蓋度等隨機因素逐漸減弱，地形起伏等結構因素逐漸加強，在本研究的最大采樣尺度上土壤有機碳的含量基本沿著坡面呈水平帶狀分布，且有由上到下逐漸增強的趨勢．

致謝 四川師范大學大精設備開放項目對本文給予了資助，謹致謝意．

［1］王琳，歐陽華，周才平，等．貢嘎山東坡土壤有機質及氮素分布特征［J］．地理學報，2004，59(6):1012－1019．

［2］GLEIXNER G．Soil organic matter dynamics:a biological perspective derived from the use of compound－specific isotopes studies［J］．Ecol Res，2013，28:683－695．

［3］DON A，RODENBEDK C，GLEIXNER G．Unexpected control of soil carbon turnover by soil carbon concentration［J］．Environ Chem Lett，2013，11:407－413．

［4］HOUGHTON R A，HACKLER J L，LAWRENCE K T．The U．S carbon budget:contributions from land use change［J］．Science，1999，285:574－578．

［5］FANG J Y，CHEN A P，PENG C H，et al．Changes in forest biomass carbon storage in China between 1949 and 1998［J］．Science，2001，292:2320－2322．

［6］王紹強，周成虎，李克讓，等．中國土壤有機碳庫及空間分布特征分析［J］．地理學報，2000，55(5):533－544．

［7］NI J．Carbon storage in terrestrial ecosystems of China:estimates at different spatial resolutions and their responses to climate change［J］．Climatic Change，2001，49(3):339－358．

［8］SIMON B B，JEFFREY E H，MICHELLE M W，et al．Multi－scale variability in soil aggregate stability:implications for understanding and predicting semi－arid grassland degradation［J］．Geoderma，2007，140(1/2):106－118．

［9］MATRIN Y，DAVID G R，ROSALINE N．Multi－scale characterization of soil variability within an agricultural landscape mosaic system in southern Cameroon［J］．Geoderma，2005，125(1/2):117－143．

［10］李凌浩，陳佐忠．草地生態(tài)系統(tǒng)碳循環(huán)及其對全球變化的響應I，碳循環(huán)的分室模型、碳輸入與貯存［J］．植物學通報，1998，15(2):14－22．

［11］王根緒，程國棟，沈永平．青藏高原草地土壤有機碳庫及其全球意義［J］．冰川凍土，2002，24(6):693－700．

［12］田玉強，歐陽華，徐興良，等．青藏高原土壤有機碳儲量與密度分布［J］．土壤學報，2008，45(5):933－942．

［13］陶貞，沈承德，高全洲，等．高寒草甸土壤有機碳儲量及其垂直分布特征［J］．地理學報，2006，61(7):720－728．

［14］徐紅靈，張宏，張偉．川西北高寒草甸土壤理化性質對土壤呼吸速率影響研究［J］．四川師范大學學報(自然科學版)，2012，35(6):835－841．

［15］ZHANG W，HANGH．Distribution characteristics of soil organic carbon of alpine meadow in the Eastern Qinghai－Tibet Plateau［J］．Wuhan University J Natural Sciences，2009，14(3):274－280．

［16］李元壽，張人禾，王根緒，等．青藏高原典型高寒草甸區(qū)土壤有機碳氮的變異特性［J］．環(huán)境科學，2009，30(6):1826－1831．

［17］四川省紅原縣志編撰委員會．紅原縣志［M］．成都:四川人民出版社，1996:25－29．

［18］王政權．地統(tǒng)計學及在生態(tài)學中的應用［M］．北京:科學出版社，1999:1－142．

［19］GOOVAERTS P．Geostatistical tools for characterizing the spatial variability of microbiological and physico－chemical soil properties［J］．Biol Fertil Soils，1998，27:315－334．

［20］劉慶，夏江寶，謝文軍．半方差函數與Moran’s I在土壤微量元素空間分布研究中的應用:以壽光市為例［J］．武漢大學學報(信息科學版)，2011，36(9):1129－1133．

［21］呂貽忠，李保國，崔燕．不同植被群落下土壤有機質和速效磷的小尺度空間變異［J］．中國農業(yè)科學，2006，39(8): 1581－1588．

［22］張宏，史培軍，鄭秋紅．半干旱地區(qū)天然草地灌叢化與土壤異質性關系研究進展［J］．植物生態(tài)學報，2001，25(3): 366 －370．

［23］李從娟，李彥，馬健．古爾班通古特沙漠土壤化學性質空間異質性的尺度特征［J］．土壤學報，2011，48(2):302－310．

［24］SEIBERT J，STENDAH I J，SORENSEN R．Topographical influences on soil properties in boreal forests［J］．Geoderma，2007，141:139－148．

［25］石福臣，李瑞利，王紹強，等．三江平原典型濕地土壤剖面有機碳及全氮分布與積累特征［J］．應用生態(tài)學報，2007，18(7):1425－1431．

［26］王艷芬，陳佐忠，LARRY T．人類活動對錫林郭勒地區(qū)主要草原土壤有機碳分布的影響［J］．植物生態(tài)學報，1998，22(6):545－551．

［27］吳琴，堯波，幸瑞新，等．鄱陽湖典型濕地土壤有機碳分布及影響因子［J］．生態(tài)學雜志，2012，31(2):313－318．

［28］傅華，陳亞明，王彥榮，等．阿拉善主要草地類型土壤有機碳特征及其影響因素［J］．生態(tài)學報，2004，24(3):653－659．

［29］劉長秀，張宏，澤柏．灌叢對川西北高寒草甸土壤資源的影響［J］．山地學報，2006，24(3):357－365．

［30］SCHLESINGER W H，PILMANIS A M．Plant－soil interactions in deserts［J］．Biogeochemistry，1998，42:169－187．

Spatial Heterogeneity of Soil Organic Carbon of Alpine Meadow in the Eastern Qinghai-Tibet Plateau

ZHANG Wei1，ZHANG Hong1，WANG Qianyun2

(1．College of Geography and Resource Science，Sichuan Normal University，Chengdu 610101，Sichuan; 2．Technical Faculty，Sichuan Normal University，Chengdu 610101，Sichuan)

Classical statistics and geostatistics were used for the analysis of the spatial variability of soil organic carbon of alpine shrub and meadow land in the eastern Qinghai-Tibet plateau with three sampling scales．The results showed that the coefficient of variation of soil organic carbon with different sampling scales were in the range of 13．4% ～28．4%，which was a moderate variation and the spatial variability of soil organic increased with the increasing of the sample scales．There was a fine spatial autocorrelation of the soil organic carbon in the alpine meadow and shrub land．In different scales，the value of C0/(C0+C1)of the soil organic carbon was in range of 57．93% ～80．92%，which indicated a medium spatial correlation，and it was shown that the spatial variability of soil organic carbon in the sample area was mainly caused by random factors，such as micro domain variation and human activities．In the same scale，fractal dimension of soil organic carbon of alpine shrub communities was less than alpine meadow community．In the sample areas，with the increase of sampling scale，such as vegetation distribution and community coverage，the random factors gradually weakened．However，the structure factors，such as the terrain，gradually strengthened．

soil organic carbon;geostatistics;spatial variability;scale effects

S159

A

1001－8395(2016)04－0602－06

10．3969/j．issn．1001－8395．2016．04．027

(編輯 鄭月蓉)

2015－03－31

國家自然科學基金(41371125)

*通信作者簡介:張 宏(1961—)，男，教授，主要從事生態(tài)資源開發(fā)與保護等方面的研究，E－mail:zh_zhhong@sicnu．edu．cn