段 菁，云海英，李 鑫，滕 飛，納 青

(內(nèi)蒙古自治區(qū)林業(yè)信息中心，內(nèi)蒙古 呼和浩特 010010)

內(nèi)蒙古地區(qū)是我國重要的牧業(yè)生產(chǎn)基地，該地區(qū)的主要覆蓋植被為草地和林地，其中可利用的草原面積約為6.8×107hm2，包含100余種適口性強(qiáng)、利用價(jià)值高的牧草[1]。隨著內(nèi)蒙古地區(qū)養(yǎng)殖業(yè)、畜牧業(yè)以及礦產(chǎn)業(yè)發(fā)展迅速，頻繁的人類活動(dòng)對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境造成極其嚴(yán)重的影響。自1978年始，退耕還林工程、“三北”防護(hù)林工程等多項(xiàng)國家重點(diǎn)工程的推進(jìn)使當(dāng)?shù)厣鷳B(tài)環(huán)境得到改善[2-3]。本研究構(gòu)建了內(nèi)蒙古草地初級(jí)凈生產(chǎn)力的估測模型，對(duì)2005—2019年間內(nèi)蒙古地區(qū)草地生產(chǎn)力進(jìn)行分析，預(yù)測不同區(qū)域、不同草地類型的生產(chǎn)力未來變化趨勢(shì)，為維護(hù)該地區(qū)生態(tài)平衡，指導(dǎo)當(dāng)?shù)匦竽翗I(yè)科學(xué)發(fā)展提供參考。

1 內(nèi)蒙古草地生產(chǎn)力估算

1.1 研究區(qū)域概況

內(nèi)蒙古地區(qū)所處緯度較高，主要以溫帶大陸性季風(fēng)氣候?yàn)橹?，降水空間分布極不均勻，年氣溫差距較大[4-5]。內(nèi)蒙古地區(qū)的主要植被覆蓋類型為草地和林地。從草地分布的空間特性來看，主要草地類型為溫性荒漠類、低地鹽化草甸類、溫性荒漠類、暖性灌草叢類、溫性草甸草原類、溫性草原類6種。6種主要草地類型占內(nèi)蒙古地區(qū)總草地面積的98.14%。東北部的草地地理?xiàng)l件較優(yōu)，牧草質(zhì)優(yōu)產(chǎn)高，以溫性草甸草原類、溫性草原類、低地鹽化草甸類和溫性山地草甸類為主。南部和中部地區(qū)牧草營養(yǎng)較為豐富，但是密度較低且產(chǎn)量不理想，只適合飼養(yǎng)牛羊等牲畜[6-7]。陰山北部和鄂爾多斯高原西部氣候干燥，牧草種類較少且產(chǎn)量較低。西部地區(qū)主要為荒漠化草原，適合飼養(yǎng)駱駝[8]。

1.2 內(nèi)蒙古草地NPP估算

草地NPP估算需要考慮光合有效輻射(APAR,Absorbed Photosynthetic Active Radiation)和實(shí)際光能利用率兩部分[9-10]。APAR的數(shù)值取決于太陽總輻射量、植被覆蓋指數(shù)(NDVI，Normalized Difference Vegetation Index)和光合有效輻射(PAR，Photosynthetically Active Radiation)。

(x,y,t)+FPARSR(x,y,t))

(1)

式(1)是光合有效輻射的計(jì)算公式。APAR(x,y,t)、SOL(x,y,t)分別表示在空間位置(x,y)處的像元在t月份吸收的光合有效輻射和太陽總輻射能量，單位均為MJ/m2。FPARNDVI、FPARSR均會(huì)對(duì)植被層對(duì)入射光合有效輻射的吸收比例造成影響，其計(jì)算公式如公式(2)所示。

(2)

式(2)中的NDVI(x,y,t)表示在空間位置(x,y)處像元的NDVI數(shù)值，NDVImax、NDVImin分別表示草地植被覆蓋率的最大值和最小值。根據(jù)內(nèi)蒙古地區(qū)現(xiàn)狀，選取NDVImax、NDVImin分別為0.634和0.023[11]。SRmax、SRmin分別表示為比值植被指數(shù)的最大值和最小值，內(nèi)蒙古地區(qū)比值植被指數(shù)的最大值和最小值取值為4.46和1.05[12]。

ε(x,y,t)=Tε1(x,y,t)×Tε2(x,y,t)×

Wε(x,y,t)×εmax

(3)

式(3)為實(shí)際光能利用率ε的計(jì)算公式。Tε1(x,y,t)、Tε1(x,y,t)為溫度脅迫系數(shù)，Wε(x,y,t)為水分脅迫系數(shù)。Tε1(x,y,t)可以反映溫度過高或過低情況下，植被光合作用效果降低而導(dǎo)致的植被初級(jí)生產(chǎn)力的降低值。Tε1(x,y,t)可以體現(xiàn)溫度從最適溫度Topt(x,y)向高溫或低溫變化時(shí)，植被光能利用率逐漸減小的趨勢(shì)。εmax為最大光能利用率，其具體數(shù)值受到草地類型的影響[13]。

(4)

式(4)是溫度脅迫系數(shù)Tε1(x,y,t)和Tε1(x,y,t)的計(jì)算公式，T(x,y)表示某個(gè)月份的平均氣溫。植被光合作用的最適溫度為研究區(qū)域內(nèi)一年中NDVI數(shù)值達(dá)到最大的月份的平均氣溫。

(5)

式(5)是水分脅迫系數(shù)Wε(x,y,t)的計(jì)算公式。EET(x,y,t)、PET(x,y,t)分別表示空間位置(x,y)處像元在t月份的實(shí)際蒸散量和潛在蒸散量，單位為mm。

為驗(yàn)證研究選用的CASA模型的有效性，選取研究區(qū)內(nèi)的1 000個(gè)隨機(jī)點(diǎn)，對(duì)比兩種方式估算的年度草地凈初級(jí)生產(chǎn)力(SNPP)結(jié)果，得到的草地NPP結(jié)果和MOD17A3-NPP結(jié)果擬合情況如圖1所示。草地NPP結(jié)果和MOD17A3-NPP結(jié)果的相關(guān)系數(shù)為0.935，具有較高的擬合度，調(diào)整之后的擬合優(yōu)化系數(shù)為0.874。隨機(jī)生成點(diǎn)中，溫性草甸草原類、溫性草原類、溫性荒漠草原類、溫性草原化荒漠類、溫性荒漠類、低地鹽化草甸類的有效隨機(jī)點(diǎn)數(shù)分別為150、397、132、33、3、181，6種類型草地的調(diào)整后擬合優(yōu)化系數(shù)分別為0.77、0.80、0.89、0.91、0.86、0.84。說明研究選用的CASA模型具有較高的精度，使用該模型模擬得到的內(nèi)蒙古草地生產(chǎn)力結(jié)果具有一定的科學(xué)性。

圖1 MOD17A3H和CASA模型估算的 NPP值結(jié)果相關(guān)性分析Fig.1 Correlation analysis of NPP value estimated by MOD17A3H and CASA model

2 內(nèi)蒙古草地NPP空間格局

2.1 草地SNPP均值空間分布特征

內(nèi)蒙古地區(qū)2005—2019年的SNPP均值整體上呈現(xiàn)出自西向東、自南向北逐漸增加的趨勢(shì)。阿拉善盟、巴彥淖爾市的SNPP值在0～100 g·C/(m2·a)，處于較低水平；呼倫貝爾市、興安盟和赤峰市部分地區(qū)的SNPP值在300 g·C/(m2·a)以上。

研究區(qū)NPP經(jīng)向變化規(guī)律見圖2所示。

圖2 內(nèi)蒙古草地平均NPP經(jīng)向變化規(guī)律Fig.2 Longitude variation of grassland average NPP in Inner Mongolia

從圖2所示的NPP經(jīng)向變化規(guī)律來看，地區(qū)NPP變化具有明顯的地帶性，上升趨勢(shì)較為明顯，變化速率約為12.69 g·C/(m2·a)。溫性草甸草原類草地的平均SNPP值最高，為309.98 g·C/(m2·a)。溫性荒漠類草地的平均SNPP值最低，為40.06 g·C/(m2·a)。結(jié)合各區(qū)域草地類型及環(huán)境特點(diǎn)可知，草甸草原和溫性草原類區(qū)域由于氣溫較低、降水較多，草地的初級(jí)凈生產(chǎn)力更高?；哪菰瓍^(qū)域的氣溫較高、降水較少，生態(tài)環(huán)境脆弱，草地的初級(jí)凈生產(chǎn)力較低。從經(jīng)緯度來看，內(nèi)蒙古地區(qū)的草地SNPP空間分布，SNPP值從西向東的變化速率大約為12.65 g·C/(m2·a·10°)。

2.2 草地SNPP時(shí)空演變特征

研究內(nèi)蒙古地區(qū)草地SNPP的時(shí)空演變特征，可以在像元尺度上計(jì)算每個(gè)柵格的變化趨勢(shì)。對(duì)15年間的SNPP值進(jìn)行擬合回歸獲取一元線性函數(shù)，函數(shù)的斜率(k值)即為年際變化率。研究使用F檢驗(yàn)方法驗(yàn)證SNPP值變化的可置信程度，其值為p。若某一區(qū)域內(nèi)k值小于0，p小于0.01，則該區(qū)域草地SNPP極顯著減少；p在0.01～0.05，則草地SNPP顯著減少；p大于0.05，則草地SNPP減小不顯著。若某一區(qū)域內(nèi)k值大于0，p大于0.05，則該區(qū)域草地SNPP增加不顯著；p在0.01～0.05，則草地SNPP顯著增加；p大于0.05，則草地SNPP極顯著增加。根據(jù)上述標(biāo)準(zhǔn)，統(tǒng)計(jì)分析了內(nèi)蒙古草地SNPP變化趨勢(shì)，如圖3所示。

圖3 內(nèi)蒙古草地SNPP值變化顯著性Fig.3 Significance of SNPP value change in Inner Mongolia

內(nèi)蒙古93.56%的地區(qū)15年間的SNPP值呈現(xiàn)增加趨勢(shì)，“顯著增加”和“極顯著增加”的比例分別為22.26%和29.76%?！帮@著減小”和“極顯著減小”的區(qū)域僅占內(nèi)蒙古地區(qū)的0.35%，主要集中在阿拉善盟的西部。從圖3中可以看出，按照草地類型劃分，“極顯著增加”區(qū)域占比較高的為低地鹽化草甸類和溫性草甸草原類，分別為48%和39.98%。占比最低的為溫性荒漠類，僅有8.44%?！皹O顯著減小”和“顯著減小”的比例在6種草地類型中占比均低于2%。除去溫性荒漠類的其他草地類型中，這2種變化趨勢(shì)不超過0.3%。

3 內(nèi)蒙古草地NPP時(shí)間變化特征分析

3.1 年變化特性

對(duì)2005—2019年內(nèi)蒙古地區(qū)草地SNPP均值的年際變化情況進(jìn)行統(tǒng)計(jì)，繪制NPP年際變化圖(圖4)。

圖4 內(nèi)蒙古地區(qū)15年間的草地NPP年際變化Fig.4 Interannual variation of grassland NPP in Inner Mongolia in 15 years

從圖4顯示的NPP年際變化情況可知，15年間的SNPP均值呈現(xiàn)出臺(tái)階式上升趨勢(shì)。2007—2008年、2012—2013年、2016—2018年為3個(gè)增長階梯。2006年的平均SNPP值最小，為27.48 g·C/(m2·a)；2018年的平均SNPP值最大，為43.54 g·C/(m2·a)。從圖中可以看出，15年間，300 g·C/(m2·a)以下的區(qū)域減少了12.84%，300 g·C/(m2·a)以上的區(qū)域增加了15.08%。當(dāng)?shù)厣鷳B(tài)環(huán)境有明顯改善，草地生產(chǎn)力穩(wěn)步提高。

6種草地類型的SNPP年際變化趨勢(shì)擬合函數(shù)見表1所示。

表1 內(nèi)蒙古地區(qū)6種草地類型的SNPP年際變化趨勢(shì)Tab.1 Interannual variation trend of SNPP of six grassland types in Inner Mongolia

6種主要草地類型的變化趨勢(shì)基本一致，均呈現(xiàn)出波動(dòng)上升趨勢(shì)。由擬合曲線的函數(shù)表達(dá)式可知，溫性草甸草原類的增速最高，其次是低地鹽化草甸類和溫性草原類。溫性荒漠類草地的增速最低，僅為0.709 6 g·C/(m2·a)。

3.2 月、季變化特性

受氣溫和濕度影響，內(nèi)蒙古地區(qū)草地NPP在每個(gè)月份的數(shù)值不同。1—7月份，氣溫逐漸升高且草地生長所需水分逐漸增加，草地NPP也逐漸提高。7月份之后，草地NPP整體處于下降趨勢(shì)。冬季氣溫低，降水降雪不斷增加，低溫高濕的環(huán)境不利于植被生長。

草地NPP的主要積累期在水熱搭配較好的4—10月，但受草地類型的影響，不同草地類型的NPP積累期峰值所處位置略有不同。6種草地類型的NPP在一年中各個(gè)月份的積累情況如圖5所示。溫性草甸草原類、溫性草原類和低地鹽化草甸類的NPP均值峰值出現(xiàn)在7月，其他3種草地類型的NPP均值峰值均出現(xiàn)在8月。該現(xiàn)象可為該地區(qū)草地管理提供一定參考，提高草地NPP的積累。

圖5 6種草地類型的NPP均值月變動(dòng)情況Fig.5 Monthly variation of average NPP value of six grassland types

3.3 未來變化趨勢(shì)

使用Hurst指數(shù)法預(yù)測內(nèi)蒙古草地NPP的未來變化趨勢(shì)，揭示草地SNPP變化規(guī)律。當(dāng)Hurst指數(shù)的數(shù)值越接近0，反持續(xù)性越強(qiáng)；H數(shù)值越接近1，持續(xù)性越強(qiáng)。

草地SNPP的平均Hurst指數(shù)為0.71。研究區(qū)域15.8%的區(qū)域呈現(xiàn)出反持續(xù)特征，84.2%的區(qū)域呈現(xiàn)出持續(xù)性特征。錫林郭勒盟、赤峰市和通遼市的大部分區(qū)域均呈現(xiàn)出“強(qiáng)持續(xù)性”和“極強(qiáng)持續(xù)性”；阿拉善盟西部和呼倫貝爾市東北部呈現(xiàn)出“極強(qiáng)反持續(xù)性”和“強(qiáng)反持續(xù)性”。研究區(qū)域草地SNPP的Hurst指數(shù)統(tǒng)計(jì)結(jié)果如圖6所示。從草地類型看，Hurst指數(shù)比重最高的為溫性荒漠草原類，為92.44%；Hurst指數(shù)比重最低的為溫性荒漠類，為73.55%。反持續(xù)性特征最為明顯的為溫性荒漠類，持續(xù)性特征最為明顯的為溫性荒漠草原類。

圖6 內(nèi)蒙古地區(qū)草地SNPP的Hurst指數(shù)空間分布Fig.6 Spatial distribution of Hurst index of grassland SNPP in Inner Mongolia

根據(jù)研究區(qū)域SNPP的Hurst指數(shù)和年際變化率，可將草地SNPP的未來變化趨勢(shì)劃分為如圖7所示的4種(圖7)。

圖7 內(nèi)蒙古地區(qū)草地SNPP未來變化趨勢(shì)Fig.7 Future trend of grassland SNPP in Inner Mongolia

年際變化率大于0，Hurst指數(shù)大于0.5時(shí)，草地SNPP值“持續(xù)增加”；Hurst指數(shù)小于0.5時(shí)，SNPP值“由增加變?yōu)闇p小”；年際變化率小于0，Hurst指數(shù)大于0.5時(shí)，草地SNPP值“持續(xù)減小”；Hurst指數(shù)小于0.5時(shí)，SNPP值“由減小變?yōu)樵黾印?。從圖中可以看出，內(nèi)蒙古大部分地區(qū)的草地SNPP值將“持續(xù)增加”。大約4.07%的區(qū)域會(huì)出現(xiàn)“持續(xù)減小”情況，主要集中在阿拉善盟西部地區(qū)。大約13.45%的區(qū)域呈現(xiàn)“由增加變?yōu)闇p小”趨勢(shì)，主要集中在呼倫貝爾市東北部。從草地類型看，SNPP“持續(xù)增加”占比最高的草地類型為溫性荒漠草原類，為91.68%；SNPP“持續(xù)減小”比重最高的草地類型為溫性荒漠類，為14.98%；SNPP“由增加變?yōu)闇p小”占比最高的類型為低地鹽化草甸類，其次為溫性草甸草原類，分別為19.43%和17.59%。需要采取一定人為措施提高溫性荒漠類草地的生產(chǎn)力，同時(shí)應(yīng)特別注意低地鹽化草甸類和溫性草甸類草地，防止其生產(chǎn)力由提升轉(zhuǎn)為降低。

4 結(jié)論

草原生態(tài)系統(tǒng)植被生產(chǎn)力的研究是全球氣候變化研究的熱點(diǎn)之一，也是區(qū)域政府制定可持續(xù)發(fā)展戰(zhàn)略的必要參考內(nèi)容。本次研究使用CASA模型對(duì)內(nèi)蒙古地區(qū)草地NPP在不同時(shí)間尺度下的時(shí)空變化特征進(jìn)行分析。使用Hurst指數(shù)法預(yù)測研究區(qū)域內(nèi)草地未來NPP變化趨勢(shì)。結(jié)果表明，內(nèi)蒙古93.56%的地區(qū)在15年間的SNPP值呈現(xiàn)增加趨勢(shì)，溫性草甸草原類草地的平均SNPP值最高。溫性草甸草原類、溫性草原類和低地鹽化草甸類的NPP均值峰值出現(xiàn)在7月，其他3種草地類型的NPP均值峰值在8月。Hurst指數(shù)比重最高的是溫性荒漠草原類，為92.44%。SNPP“由增加變?yōu)闇p小”占比最高的類型是低地鹽化草甸類，其次是溫性草甸草原類，分別為19.43%和17.59%。上述研究成果將為內(nèi)蒙古地區(qū)生態(tài)治理和經(jīng)濟(jì)發(fā)展規(guī)劃提供一定參考依據(jù)，但是有關(guān)草地生產(chǎn)力未來變化趨勢(shì)的預(yù)測并未經(jīng)過實(shí)際驗(yàn)證，后續(xù)還需要繼續(xù)跟蹤研究。