王 利，徐翠玲，蔣鎰竹

(1.長(zhǎng)安大學(xué)地質(zhì)工程與測(cè)繪學(xué)院，陜西 西安 710054；2.長(zhǎng)安大學(xué)地球科學(xué)與資源學(xué)院，陜西 西安 710054)

0 引 言

植被是陸地生態(tài)系統(tǒng)與氣候系統(tǒng)生物圈的重要組成部分，也是氣候和生態(tài)環(huán)境變化的指示器[1-2]。植被與大氣、水等其他圈層之間相互作用，通過(guò)光合、呼吸等生物物理過(guò)程實(shí)現(xiàn)土壤-植被-大氣之間的物質(zhì)、能量和動(dòng)量交換[3- 4]，對(duì)植被的發(fā)展、物候等具有顯著影響。葉面積指數(shù)(Leaf Area Index，LAI)是指單位水平土地面積上植物葉片表面積的一半[5-7]。葉面積指數(shù)可以反映植被個(gè)體或群體的長(zhǎng)勢(shì)特征[8]，與植被的光合作用、蒸騰作用、降雨攔截和碳通量等生物物理過(guò)程密切相關(guān)，在植被動(dòng)態(tài)模擬觀測(cè)、生態(tài)系統(tǒng)循環(huán)研究、陸地表面蒸散、農(nóng)業(yè)遙感估產(chǎn)、碳蓄積量、水土保持等方面應(yīng)用廣泛?？傊?，植被葉面積指數(shù)的生態(tài)學(xué)意義明顯，利用葉面積指數(shù)研究區(qū)域尺度上的植被時(shí)空變化對(duì)區(qū)域生態(tài)環(huán)境時(shí)空演變尤為重要[9]。

近年來(lái)，國(guó)內(nèi)外學(xué)者對(duì)葉面積指數(shù)的時(shí)空動(dòng)態(tài)變化進(jìn)行了一系列研究。例如，Dusseux等基于SPOT和Quickbird遙感影像估算葉面積指數(shù)，利用葉面積指數(shù)的時(shí)間演變來(lái)識(shí)別草地的動(dòng)態(tài)變化[10]。任宏昌等采用MODIS LAI產(chǎn)品，分析中國(guó)地區(qū)LAI的時(shí)空變化特征及其對(duì)氣候因子響應(yīng)，結(jié)果表明，LAI時(shí)空變化及其對(duì)氣候條件響應(yīng)具有一致性[11]。丹利等發(fā)現(xiàn)，貴州地區(qū)植被的LAI對(duì)降水的響應(yīng)程度大于氣溫，尤其是冬季和春季，在夏季和秋季植被的分布主要受氣溫的影響[9]。段利民等研究發(fā)現(xiàn)，錫林河流域近30年LAI總體呈上升趨勢(shì)，高程對(duì)LAI的影響遠(yuǎn)大于坡度和坡向?qū)AI的影響[12]。張宇佳等發(fā)現(xiàn)，降水、氣溫、地形、植被類型和物候期等因子對(duì)LAI的變化均有顯著影響，且各因子對(duì)植被LAI的影響程度不同[13]。

黑河流域是我國(guó)第二大內(nèi)陸河流域，地處西北內(nèi)陸干旱和半干旱區(qū)，生態(tài)環(huán)境十分脆弱[14]。由于黑河流域植被種類較多，且不同植被類型具有明顯的季節(jié)性與地域性；同時(shí)，LAI的變化趨勢(shì)能夠較好的反映植被在數(shù)量與結(jié)構(gòu)方面的變化，因而LAI可以作為衡量植被變化的敏感性指標(biāo)，這對(duì)黑河流域的生態(tài)環(huán)境保護(hù)以及探究LAI在黑河流域中的應(yīng)用具有重要意義[15]。因此，為評(píng)估黑河流域近年來(lái)植被的變化情況，采用美國(guó)NASA發(fā)布的MODIS LAI數(shù)據(jù)和土地覆蓋類型數(shù)據(jù)，利用趨勢(shì)分析及相關(guān)分析等方法，揭示2001年～2016年黑河流域LAI的時(shí)空變化特征，同時(shí)結(jié)合同期氣象資料分析葉面積指數(shù)變化的主要原因，以期為黑河流域生態(tài)環(huán)境保護(hù)與治理提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)概況

黑河流域是我國(guó)西北地區(qū)第二大內(nèi)陸河流域，位于河西走廊中部[16]，地理位置介于37°～43°N，97°～102°E之間(見圖1)。其干流全長(zhǎng)821 km，流域面積約14.29萬(wàn)km2，橫跨青海、內(nèi)蒙古、甘肅3個(gè)省份，上游屬于青海省祁連縣，中游屬于甘肅省山丹、民樂(lè)、張掖、臨澤、高臺(tái)、肅南和酒泉等市縣，下游屬于甘肅省金塔市和內(nèi)蒙古自治區(qū)額濟(jì)納旗[17]。黑河流域位于歐亞大陸中部，遠(yuǎn)離海洋，周圍高山環(huán)繞，流域氣候主要受中高緯度的西風(fēng)帶環(huán)流控制和極地冷氣團(tuán)影響，氣候干燥，降水稀少且集中，多大風(fēng)，日照充足，太陽(yáng)輻射強(qiáng)烈，晝夜溫差較大，為典型的大陸性氣候[18]。流域內(nèi)垂直高度轉(zhuǎn)變較大，海拔分布呈現(xiàn)出由南到北順次下降的趨向。上游海拔一般在2.2 km以上，主要以高山、冰川、森林、草甸和草原為主，植被覆蓋度較高，中游海拔在1.0～2.0 km之間，下游海拔在0.98～1.2 km間，中、下游主要是海拔為1.3～2.3 km的走廊地帶，分布著大面積的戈壁和沙漠，中下游地帶性植被為溫帶小灌木和半灌木荒漠植被，綠洲植被大多以胡楊、檉柳、沙棗和梭梭為主，植被種類簡(jiǎn)單，覆蓋度小[19]。

圖1 黑河流域地理位置

1.2 數(shù)據(jù)來(lái)源

葉面積指數(shù)數(shù)據(jù)采用美國(guó)國(guó)家航空航天局提供的MODIS/Terra衛(wèi)星4級(jí)8d合成陸地產(chǎn)品MOD15A2H，數(shù)據(jù)格式為HDF，空間分辨率為500 m，時(shí)間跨度為2001年～2016年[20]。LAI值的有效范圍為0～100，比例因子為0.1。

土地覆蓋類型數(shù)據(jù)采用美國(guó)國(guó)家航空航天局提供的MODIS/Terra&Aqua衛(wèi)星3級(jí)土地覆蓋數(shù)據(jù)MCD12Q1，空間分辨率為500 m，時(shí)間分辨率為1a[21]，時(shí)間跨度為2001年～2016年。該數(shù)據(jù)采用國(guó)際地圈生物圈計(jì)劃(IGBP)分類方法，將土地覆蓋類型分為常綠闊葉林、常綠針葉林、落葉闊葉林、針葉闊葉林、混交林、郁閉灌叢、稀疏灌叢、多樹草原、稀樹草原、草地、永久濕地、農(nóng)田、城市和建筑用地、農(nóng)田/自然植被混合、冰雪、裸地或低植被覆蓋區(qū)域、水體等17種類型[22]。

氣象數(shù)據(jù)采用同期黑河流域內(nèi)山丹、祁連、張掖、野牛溝、托勒、高臺(tái)、鼎新、金塔、酒泉、馬鬃山、額濟(jì)納旗、剛察、門源、阿拉善右旗、永昌和玉門等16個(gè)氣象站點(diǎn)的逐日實(shí)測(cè)數(shù)據(jù)，包括平均氣溫(℃)和降水量(mm)。黑河流域邊界數(shù)據(jù)和數(shù)字高程模型數(shù)據(jù)來(lái)源于國(guó)家青藏高原科學(xué)數(shù)據(jù)中心(http:∥data.tpdc.ac.cn)[23-24]。

1.3 研究方法

首先，為反映植被LAI的時(shí)間變化特點(diǎn)，依據(jù)研究區(qū)植被的物候特征，采用3月～5月為春天、6月～8月為夏天、9月～11月為秋天、12月到次年2月為冬日的季節(jié)劃分原則，用最大值合成法(MVC)對(duì)葉面積指數(shù)進(jìn)行逐月、逐年合成，完成月、年葉面積指數(shù)提取并求取多年葉面積指數(shù)的月平均值及年平均值，以分析葉面積指數(shù)的年內(nèi)變化特征，同時(shí)利用變異系數(shù)(Coefficient of Variance，CV)表征葉面積指數(shù)的年際變化[25]。

其次，基于像元尺度，采用Slope趨勢(shì)分析法，模擬研究區(qū)中每個(gè)柵格單元的變化趨勢(shì)，分析黑河流域2001年～2016年16年間LAI的空間變化趨勢(shì)[26-27]。

最后，對(duì)平均氣溫與降水量數(shù)據(jù)采用Kriging方法進(jìn)行空間插值[28]。在此基礎(chǔ)上，分析2001年～2016年LAI與平均氣溫和降水量的相關(guān)系數(shù)，以探索長(zhǎng)時(shí)序LAI與平均氣溫、降水量之間的相關(guān)程度及其空間分布規(guī)律。由于地理系統(tǒng)中各要素間相互影響，采用偏相關(guān)系數(shù)進(jìn)一步研究LAI與單個(gè)氣候因子之間的相關(guān)性，以探索導(dǎo)致黑河流域LAI變化的主要因素[29-31]。

2 結(jié)果與分析

2.1 黑河流域LAI時(shí)間變化分析

2.1.1LAI年內(nèi)變化

黑河流域LAI年內(nèi)變化如圖2所示。在整個(gè)流域內(nèi)，由于2月、3月、4月氣溫較低和降水量少，植被大多處于無(wú)葉或少葉狀態(tài)，LAI值較低，4月份LAI值達(dá)到全年最低；4月以后，隨著氣溫逐漸回升，降水量逐漸增加，植被開始生長(zhǎng)，葉面積指數(shù)呈逐漸增大的趨勢(shì)；到5月以后，葉面積指數(shù)大幅上升，進(jìn)入夏季以后，光照充足，降水豐沛，有利于植被生長(zhǎng)，葉面積指數(shù)達(dá)到全年最大值，即7月的3.49；8月以后氣溫開始下降，植被進(jìn)入落葉期，LAI逐漸降低；從9月開始葉面積指數(shù)降低速度較快，但由于下游地區(qū)存在大量稀疏植被，11月LAI又開始逐漸增加。分歧季候，葉面積指數(shù)的空間散布差異性分歧。結(jié)合LAI年內(nèi)變化情況(見圖3)，春季、秋季葉面積指數(shù)值較低，夏季、冬季葉面積指數(shù)值較高；同時(shí)中上游地區(qū)葉面積指數(shù)值較高，植被覆蓋度較高，下游葉面積指數(shù)較低，植被覆蓋度很低。

圖2 2001年～2016年黑河流域1月～12月LAI平均值

圖3 2001年～2016年黑河流域1月～12月LAI平均值空間分布

2.1.2LAI年際變化

2001年～2016年黑河流域LAI年際變化呈現(xiàn)出在波動(dòng)中上升的趨勢(shì)(見圖4)，年LAI值變化范圍為2.34～2.90，多年LAI均值為2.61。LAI值從2001年的2.34增加到2016年的2.82，升幅為20.51%，總體呈上升趨勢(shì)。其中，最小值和最大值分別出現(xiàn)在2001年和2015年，對(duì)應(yīng)的LAI值分別為2.34和2.90；2001年和2013年葉面積指數(shù)平均值明顯低于多年LAI均值，這是由于這2年發(fā)生了干旱。2001年1月至7月，流域范圍內(nèi)大部分地區(qū)出現(xiàn)了連旱天氣，2013年黑河流域降水量少且夏季氣溫偏高，阻礙了植被生長(zhǎng)，故葉面積指數(shù)較低[32-33]。

圖4 2001年～2016年黑河流域葉面積指數(shù)年平均值

2.2 黑河流域LAI空間分布變化趨勢(shì)分析

2001年～2016年黑河流域葉面積指數(shù)存在顯著的空間分布差異(見圖5)，葉面積指數(shù)上升的區(qū)域主要集中在流域中上游，下游地區(qū)由于土地沙化和水資源緊缺等原因，16年間葉面積指數(shù)呈下降趨勢(shì)。同時(shí)，黑河流域葉面積指數(shù)變異系數(shù)高值區(qū)域主要集中在黑河流域中上游，CV>1.0，部分區(qū)域達(dá)到了2.0以上，年際間波動(dòng)幅度較大，但分布范圍很小，這表明中上游地區(qū)LAI年際波動(dòng)范圍較大的區(qū)域很少；流域中上游地區(qū)大部分區(qū)域LAI年際間波動(dòng)程度較小，CV值介于0～1之間；連系其線性轉(zhuǎn)變趨向得悉，均為葉面積指數(shù)上升的原因。黑河流域下游大部分地區(qū)幾乎沒有較大的年際波動(dòng)，表明葉面積指數(shù)16年間表現(xiàn)相對(duì)較為平穩(wěn)。

圖5 2001年～2016年黑河流域葉面積指數(shù)空間變化趨勢(shì)

2.3 不同土地覆蓋類型LAI變化分析

2.3.1 不同土地覆蓋類型LAI年內(nèi)變化

2001年～2016年不同土地覆蓋類型1月～12月LAI平均值如圖6所示，農(nóng)用地、濕地、木本熱帶稀樹草原、熱帶稀樹草原、常綠針葉林、稀疏植被、稀疏灌叢年內(nèi)變化呈單峰變化趨勢(shì)，草地葉面積指數(shù)呈雙峰變化趨勢(shì)。農(nóng)用地、濕地、木本熱帶稀樹草原、熱帶稀樹草原、常綠針葉林葉面積指數(shù)在4月、5月迅速升高，7月達(dá)到最大值，8月葉面積指數(shù)開始逐漸降低，11月下降至1.0以下，之后降落速率變緩，冬日其葉面積指數(shù)無(wú)明顯變化。而草地和稀疏植被的葉面積指數(shù)年內(nèi)轉(zhuǎn)變較大，稀疏植被的葉面積指數(shù)從9月起迅速上升，12月到達(dá)最大值，到次年1月，葉面積指數(shù)才開始逐步下降；草地的葉面積指數(shù)在4月和5月敏捷升高，7月到達(dá)夏日最大值3.15，8月葉面積指數(shù)逐漸降低，10月以后又開始迅速上升，12月達(dá)到全年最大值3.33，然后到次年2月迅速下降至1.0以下。稀疏灌叢的葉面積指數(shù)在4月、5月逐步升高，6月達(dá)到最大值0.14，9月開始逐漸降低，下降速度緩慢，冬季葉面積指數(shù)無(wú)明顯變化。

圖6 2001年～2016年不同土地覆蓋類型1月～12月LAI平均值

2.3.2 不同土地覆蓋類型LAI年際變化

根據(jù)2001年～2016年MCD12Q1分類數(shù)據(jù)成果統(tǒng)計(jì)，黑河流域內(nèi)包含熱帶稀樹草原、常綠針葉林、草地、濕地、稀疏灌叢、農(nóng)用地、稀疏植被、木本熱帶稀樹草原、水體、城市及建筑區(qū)和冰雪11種土地覆蓋類型。2001年～2016年不同土地覆蓋類型LAI年平均值的變化如圖7所示，城市及建筑用地、冰雪和水體葉面積指數(shù)為0，常綠針葉林、稀疏植被、草地、濕地、木本熱帶稀樹草原、農(nóng)用地、熱帶稀樹草原7種土地覆蓋類型的葉面積指數(shù)總體均呈上升趨勢(shì)。其中，稀疏植被植被葉面積指數(shù)上升速度最快，為0.049a-1；木本熱帶稀樹草原次之，為0.039 a-1；草地葉面積指數(shù)上升速率最慢，為0.013 a-1?？傮w而言，黑河流域不同土地覆蓋類型的葉面積指數(shù)差異不大，除了熱帶稀樹草原和稀疏灌叢的葉面積指數(shù)多年平均值小于1.0，其他土地覆蓋類型的葉面積指數(shù)波動(dòng)范圍均介于[1，2]之間。

圖7 2001年～2016年不同土地覆蓋類型LAI年平均值

2.4 氣象因素對(duì)LAI變化的影響

黑河流域LAI及各氣象要素月平均值見表1。由表1可知，LAI從4月開始逐步上升，7月達(dá)到整年最大值3.49；8月～10月LAI開始逐漸降低，10～12月又開始逐漸增加，12月到次年4月逐漸降低并在4月達(dá)到全年最小值1.28。結(jié)合植被物候特點(diǎn)，黑河流域LAI表現(xiàn)出明顯的季節(jié)變化。同時(shí)，LAI的年內(nèi)轉(zhuǎn)變可以直接反映出植被的物候，從4月起頭，各種植被開始漸漸生長(zhǎng)；到7月，植被處于生長(zhǎng)發(fā)展旺盛期，此時(shí)LAI到達(dá)最大值；9月，植被進(jìn)入落葉期，植被覆蓋度逐漸減小，LAI逐漸降低；但由于下游植被覆蓋類型的不同，10月LAI又逐漸增加，12月LAI達(dá)到峰值，稀疏植被處于生長(zhǎng)旺盛期；12月以后，LAI逐漸降低，稀疏植被逐漸停止生長(zhǎng)，次年4月LAI達(dá)到最小值。此外，根據(jù)表1可以得知，黑河流域的月平均氣溫和月降水量的年內(nèi)變化趨勢(shì)與LAI的年內(nèi)變化趨勢(shì)并不一致，氣象要素年內(nèi)呈單峰趨勢(shì)，冬季平均氣溫較低，降水量較少，夏季平均氣溫較高，降水量較多。月平均氣溫最大值出現(xiàn)在7月，為19.93℃，最低值出現(xiàn)在1月，為-10.81 ℃；月降水量最大值也出現(xiàn)在7月，為47.89 mm，最小值出現(xiàn)在2月，為1.51 mm。

表1 黑河流域LAI及各氣象要素月平均值

黑河流域LAI及各氣象要素年際變化趨勢(shì)見圖8。由圖8可得，黑河流域2001年～2016年年平均氣溫和年降水量的年際變化趨勢(shì)與年平均LAI幾乎一致，總體都呈緩慢上升的趨勢(shì)，表明16年間黑河流域內(nèi)植被呈緩慢增加的趨勢(shì)，且氣候總體上呈暖濕化趨勢(shì)。其中，年平均氣溫和年降水量的升幅分別為10.83%和40.04%。

圖8 黑河流域LAI及各氣象要素年際變化趨勢(shì)

對(duì)黑河流域2001年～2016年年平均LAI與年平均氣溫和年降水量開展逐像元相干闡發(fā)、偏相干闡發(fā)，成果如圖9所示。

圖9 LAI與平均氣溫、降水量相關(guān)性空間分布

圖9a顯示了黑河流域年平均LAI與年降水量的相關(guān)系數(shù)的空間分布情況，其值介于-0.59～0.81之間。圖9c顯示了黑河流域年平均LAI與年降水量的偏相關(guān)系數(shù)的空間分布情況，其值介于-0.71～0.85之間。由圖可知，二者的相關(guān)系數(shù)和偏相關(guān)系數(shù)的極化現(xiàn)象嚴(yán)重且空間差異明顯。年平均LAI與年降水量呈正相關(guān)關(guān)系的區(qū)域主要分布在黑河流域中上游，包括祁連、野牛溝、托勒、山丹、張掖、高臺(tái)、酒泉、金塔，這些地區(qū)植被以森林、草甸草原為主，植被覆蓋度較高，且年平均氣溫較低，同期氣溫升高，降水量增多，有利于植被生長(zhǎng)。年平均LAI與年降水量呈負(fù)相關(guān)關(guān)系的地區(qū)分布在黑河流域下游大部分地區(qū)和中游少部分地區(qū)，包括鼎新、玉門、阿拉善右旗、額濟(jì)納旗、馬鬃山。結(jié)合氣象站點(diǎn)的統(tǒng)計(jì)數(shù)據(jù)，這些地區(qū)年降水量較少；同時(shí)，根據(jù)相關(guān)分析原理得知，相關(guān)系數(shù)的絕對(duì)值越大，變量之間的相關(guān)性越強(qiáng)；故表明降水對(duì)植被生長(zhǎng)影響較大。圖9b顯示了黑河流域年平均LAI與年平均氣溫的相關(guān)系數(shù)的空間分布情況，相關(guān)系數(shù)的值介于-0.68～0.86之間。圖9d顯示了黑河流域年平均LAI與年平均氣溫的偏相關(guān)系數(shù)的空間分布情況，偏相關(guān)系數(shù)介于-0.73～0.85之間。從空間分布上可以看出，年平均LAI與年平均氣溫呈正相關(guān)關(guān)系的區(qū)域主要分布在中上游地區(qū)，尤其是祁連、野牛溝、托勒、山丹、張掖、高臺(tái)地區(qū)，年平均LAI與年平均氣溫的正相關(guān)性較明顯。年平均LAI與年平均氣溫呈負(fù)相關(guān)關(guān)系的區(qū)域主要分布在中下游地區(qū)?？傮w而言，黑河流域的氣溫和降水對(duì)植被的生長(zhǎng)影響顯著，且二者對(duì)植被LAI的響應(yīng)程度幾乎相同。

3 結(jié) 論

本文以黑河流域?yàn)檠芯繀^(qū)，分析LAI的時(shí)空變化及其與氣候因子的相關(guān)性。結(jié)果表明，2001年～2016年16年間黑河流域的氣候呈暖濕化趨勢(shì)，總體而言流域內(nèi)植被狀況也處于改善趨勢(shì)。16年間黑河流域植被LAI年內(nèi)呈雙峰分布，年際間總體呈上升趨勢(shì)，但空間差異明顯。在空間分布上，流域上游、中游植被活動(dòng)的強(qiáng)度明顯大于下游地區(qū)；因此，LAI增大的趨勢(shì)表現(xiàn)為上游強(qiáng)于中游，中游強(qiáng)于下游，且植被顯著改善區(qū)域主要集中分布在上游寒區(qū)試驗(yàn)區(qū)和中游人工綠洲試驗(yàn)區(qū)。期間，黑河流域不同土地覆蓋類型植被狀況表現(xiàn)出良好的改善趨勢(shì)，常綠針葉林、木本熱帶稀樹草原、稀疏植被、濕地、農(nóng)用地、草地、熱帶稀樹草原總體恢復(fù)情況較好。在LAI與氣候因子的相關(guān)性分析中發(fā)現(xiàn)，LAI對(duì)氣溫和降水的響應(yīng)強(qiáng)度大致相同。但黑河流域LAI與氣候因子的相關(guān)性空間差異明顯，中上游主要呈正相關(guān)關(guān)系，下游主要呈負(fù)相關(guān)關(guān)系。

其中，LAI低值年份均對(duì)應(yīng)于極端氣候條件，反映出氣候變化對(duì)黑河流域植被動(dòng)態(tài)變化的影響是明顯的，且黑河流域植被顯著恢復(fù)區(qū)也在一定程度上反映出生態(tài)環(huán)境治理工程對(duì)植被恢復(fù)起到了顯著的促進(jìn)作用，這表明氣候暖濕化和生態(tài)環(huán)境治理工程的有效實(shí)施對(duì)黑河流域的植被恢復(fù)起到了至關(guān)重要的作用，未來(lái)應(yīng)加大力度實(shí)施生態(tài)環(huán)境治理工程，例如濕地工程等，以保護(hù)流域生態(tài)環(huán)境。此外，本文只針對(duì)黑河流域植被總體的LAI進(jìn)行分析，沒有具體研究不同土地覆蓋類型在不同區(qū)域植被LAI的動(dòng)態(tài)變化及其原因，這有待進(jìn)一步嘗試，以便為相關(guān)部門制定對(duì)應(yīng)的生態(tài)環(huán)境保護(hù)與治理策略。同時(shí)，除了氣候因素外，葉面積指數(shù)還受到植被類型、植被搭配比重、植被生長(zhǎng)年齡、DEM等因素的影響，未來(lái)的研究應(yīng)全面考慮各個(gè)因子對(duì)葉面積指數(shù)的影響，進(jìn)行深入全面的研究。

致謝：本文所用的MODIS數(shù)據(jù)來(lái)源于美國(guó)國(guó)家航空航天局，黑河流域氣象站觀測(cè)數(shù)據(jù)由“溫室數(shù)據(jù)共享平臺(tái)”提供，黑河流域邊界數(shù)據(jù)和數(shù)字高程模型數(shù)據(jù)來(lái)源于 “國(guó)家青藏高原科學(xué)數(shù)據(jù)中心”(原“中國(guó)西部生態(tài)與環(huán)境科學(xué)數(shù)據(jù)中心”)，從“中國(guó)西部生態(tài)與環(huán)境科學(xué)數(shù)據(jù)中心”下載，在此一并表示感謝!