秦巴山地NPP及對氣候變化響應(yīng)的多維地帶性與暖溫帶-亞熱帶界線

趙 芳, 張久陽, 劉思遠(yuǎn), 王增艷,*, 王黎歡, 顧浩婷, 李萬隆

1 河南大學(xué)環(huán)境與規(guī)劃學(xué)院, 開封 475004

2 黃河中下游數(shù)字地理技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室, 開封 475004

秦嶺-大巴山(秦巴山地)地處我國中部,承接青藏高原與東部平原,連接長江流域與黃河流域,是我國的南北過渡帶。作為連接中國東西部的唯一一條大尺度生態(tài)廊道,秦巴山地影響我國生態(tài)地理格局、生物區(qū)系演化、自然資源分布[1]。受地理位置、山系分布等綜合作用,秦巴山地形成了復(fù)雜多樣且具有過渡性質(zhì)的植被體系和氣候類型[2-4],導(dǎo)致植被對氣候變化響應(yīng)的復(fù)雜性。

植被凈第一性生產(chǎn)力(Net Primary Productivity, NPP),是在單位面積單位時(shí)間上綠色植被固定的總能量與自身呼吸消耗掉的有機(jī)物之差,反映了自然環(huán)境條件下植被群落的生產(chǎn)能力[5-6]。研究發(fā)現(xiàn),秦巴山地地表NPP分布及氣候響應(yīng)隨經(jīng)度、緯度等的變化在不同地段呈現(xiàn)一定的規(guī)律性[7-8],可表征秦巴山地復(fù)雜的生態(tài)地理格局。秦嶺山地氣候和植被區(qū)系由南向北發(fā)生了明顯的變化,氣候由亞熱帶向暖溫帶過渡,地帶性植被從常綠闊葉林逐漸轉(zhuǎn)變?yōu)槁淙~闊葉林[9-13],而地表植被NPP呈現(xiàn)減少的趨勢,其中南段大巴山區(qū)雨水豐沛、氣候適宜,為NPP高值區(qū)[14]。秦巴山地植被NPP數(shù)值呈現(xiàn)出西高東低的分布,年際變化以上升趨勢為主,且大部分地區(qū)在波動(dòng)中表現(xiàn)出較為顯著的增長趨勢[15-17]。伏牛山NPP東部數(shù)值較小,中西部數(shù)值較大,隨時(shí)間推移數(shù)值增大[18]。秦巴山地NPP及氣候響應(yīng)隨緯度、經(jīng)度、高度、坡向等多種分異因素的變化呈現(xiàn)的規(guī)律性也叫做多維地帶性,體現(xiàn)了秦巴山地復(fù)雜的地域分異,是造成我國南北分界線位置存在爭議,比如秦嶺主脊[10, 19]、秦嶺南坡[9, 13]和秦嶺北坡[11, 20-21]的因素之一。

但是現(xiàn)有關(guān)于秦巴山地NPP及其氣候響應(yīng)的研究多集中在單一山體、單一流域,如蔣沖等對陜南漢江流域及大巴山植被NPP變化的分析[14]、李敏等對神農(nóng)架林區(qū)生態(tài)系統(tǒng)NPP的估算[22]及陳志超等對伏牛山兩側(cè)NPP的分析[18]等,冠名秦巴山地的文獻(xiàn)幾乎都局限于陜西境內(nèi)的秦巴山地[23-24],極少將其作為整體進(jìn)行分析,導(dǎo)致對秦巴山地植被NPP的變化認(rèn)識(shí)不夠全面。實(shí)際上秦巴山地是一個(gè)完整的、由南向北過渡的地域單元,其南部是亞熱帶,北部是暖溫帶,植被具有天然的南北過渡性質(zhì),研究秦巴山地NPP時(shí)空分異及其對氣候變化響應(yīng)的多維地帶性,對探索秦巴山地復(fù)雜的生態(tài)地理格局、評價(jià)中國陸地生態(tài)系統(tǒng)的復(fù)雜性和多樣性具有重要的意義[25-26],也可為解決亞熱帶與暖溫帶分界線具體位置在哪里的問題提供新的研究思路。

近年來,隨著現(xiàn)代對地觀測、GIS和遙感技術(shù)的發(fā)展,秦巴山地的生物地理研究由過去的樣點(diǎn)和樣區(qū)的調(diào)查逐漸發(fā)展到使用遙感技術(shù)對整個(gè)區(qū)域生物生產(chǎn)力和生物量的調(diào)查?；诖?本文使用2000—2015年秦巴山地地表植被NPP遙感數(shù)據(jù)和氣象站點(diǎn)數(shù)據(jù),從經(jīng)度、緯度、高度和坡向等角度分析秦巴山地地表植被NPP及其對氣候變化的響應(yīng)的多維地帶性,深入討論我國南北過渡地帶的復(fù)雜生態(tài)地理格局,為暖溫帶-亞熱帶界線的界定提供新的佐證。

1 數(shù)據(jù)和方法

1.1 數(shù)據(jù)及處理

1.1.1DEM數(shù)據(jù)

本研究采用的數(shù)字高程模型(DEM)數(shù)據(jù)(來源于https://earthexplorer.usgs.gov/)空間分辨率為30m,分塊下載獲取秦巴山地不同區(qū)域數(shù)據(jù),將原始DEM拼接為一幅完整的秦巴山地高程數(shù)據(jù),如圖1所示。

圖1 秦巴山地DEM和氣象站點(diǎn)的分布

1.1.2NPP數(shù)據(jù)

研究采用的2000—2015年植被NPP數(shù)據(jù)源于美國NASA EOS/MODIS的MOD17A3產(chǎn)品(http://www.nasa.gov/)。該數(shù)據(jù)基于的BIOME-BGC生態(tài)模型,將衛(wèi)星獲得的土地覆蓋、部分光合有效輻射(FPAR)和葉面積指數(shù)(LAI)數(shù)據(jù)作為輸入的地表植被信息,模擬獲得的全球年NPP數(shù)據(jù)。該數(shù)據(jù)與同期MODIS NDVI數(shù)據(jù)的對比分析顯示,二者整體變化趨勢一致,且相關(guān)性較強(qiáng),證明該NPP數(shù)據(jù)和NDVI數(shù)據(jù)一樣具有反映區(qū)域內(nèi)植被覆蓋狀況和生長情況的能力。其時(shí)間分辨率為年,空間分辨率為1km×1km[27]。該數(shù)據(jù)已在不同區(qū)域NPP估算、植被長勢、環(huán)境評價(jià)等研究中得到驗(yàn)證和廣泛應(yīng)用[28-30]。在MRT軟件的支持下,將MOD17A3數(shù)據(jù)進(jìn)行格式轉(zhuǎn)換(HDF轉(zhuǎn)化為TIFF)、重投影、裁剪、比例換算等,得到單位為gC/m2的年均NPP數(shù)據(jù)。不過需要說明的是該NPP模擬數(shù)據(jù)只包括植物地上生產(chǎn)力部分而未包含地下生產(chǎn)力部分,未來仍需要進(jìn)一步修訂。

1.1.3氣象站點(diǎn)數(shù)據(jù)

氣象數(shù)據(jù)來自國家氣象信息中心(http://data.cma.cn/)的秦巴山地1951—2017年118個(gè)氣象站點(diǎn)每日氣溫和降水?dāng)?shù)據(jù),經(jīng)過月均計(jì)算得到月均溫、月均降水?dāng)?shù)據(jù),舍去數(shù)據(jù)缺測的部分站點(diǎn),實(shí)際使用93個(gè)氣象站點(diǎn)數(shù)據(jù),氣象站點(diǎn)主要分布在甘肅、陜西、河南、四川、湖北、重慶6個(gè)省市(如圖1)。

1.2 方法

1.2.1相關(guān)分析

相關(guān)性分析方法可以直觀反映要素間的相關(guān)關(guān)系。將NPP與氣溫和降水進(jìn)行相關(guān)分析,根據(jù)每個(gè)像元的NPP值與氣象因子的相關(guān)系數(shù)Rxy確定NPP與氣溫和降水之間的關(guān)系：

(1)

1.2.2趨勢分析

通過一定時(shí)間段內(nèi)(2000—2015年)NPP像元的數(shù)值變化,建立基于像元的一元線性回歸方程,NPP年際變化速率可以通過該方程的趨勢線斜率K來確定,其中：

(2)

K為趨勢線斜率,即NPP的年際變化速率,單位為gC m-2a-1；yi表示第i年的NPP值；i=1, 2, 3, …,16；n=16為監(jiān)測年數(shù)。當(dāng)K>0時(shí),表示秦巴山地NPP呈增加趨勢,反之則呈減少趨勢。

1.2.3地形因子分析

為了展示秦巴山地NPP與經(jīng)緯度和地形的關(guān)系,本研究選取33.95°N(經(jīng)過迭山、太白山、蟒嶺、伏牛山)、32.67°N(經(jīng)過雪寶頂、米倉山、十堰市等地區(qū)),107.81°E(經(jīng)過米倉山、漢江、太白山)和110.53°E(經(jīng)過神農(nóng)架、武當(dāng)山、漢江等地區(qū))4條剖面線(如圖1)對應(yīng)像元點(diǎn)的DEM數(shù)據(jù)和NPP數(shù)據(jù),制作出秦巴山地NPP典型剖面線分布圖,反映NPP隨地形和經(jīng)緯度的變化。

1.2.4空間統(tǒng)計(jì)分析

對氣象站點(diǎn)數(shù)據(jù)進(jìn)行Kriging空間插值,獲得與植被NPP柵格大小一致(1km×1km)、投影方式(D_WGS_1984)相同的氣象柵格數(shù)據(jù),與秦巴山地NPP進(jìn)行相關(guān)分析研究秦巴山地2000—2015年NPP變化與氣溫、降水的關(guān)系。

為了分析NPP與海拔關(guān)系,將秦巴山地DEM數(shù)據(jù)基于不同海拔進(jìn)行重分類,分為0—500、500—1000、1000—1500、1500—2000、2000—2500、2500—3000、3000—3500、>3500m共8個(gè)級(jí)別,統(tǒng)計(jì)各分段NPP均值,研究NPP沿海拔梯度的變化。

為了分析秦巴山地NPP及氣候響應(yīng)的地域分異,基于秦嶺、大巴山主山脊線和主要河流(嘉陵江、漢江)的位置,將秦巴山地劃分為秦嶺南坡、秦嶺北坡、大巴山南坡、大巴山北坡、西秦嶺,分區(qū)統(tǒng)計(jì)秦巴山地NPP及與氣溫和降水的相關(guān)在秦嶺南坡、秦嶺北坡,大巴山南坡和大巴山北坡的變化(由于西秦嶺復(fù)雜的地理結(jié)構(gòu),本文不做討論),為南北分界線的確定提供佐證。

2 結(jié)果與分析

2.1 秦巴山地地表植被NPP時(shí)空格局

根據(jù)2000—2015年多年年均NPP分布圖(圖2),秦巴山地NPP值存在明顯的區(qū)域差異,總體表現(xiàn)為南高北低,中部高,東西低的特征。其中在岷山以東104°E—105°E之間達(dá)到最高值(約700—900 gC/m2),在秦巴山地中部105°E—107°E之間NPP值出現(xiàn)高值集聚區(qū),而在秦嶺山脈以東的華山、崤山、伏牛山、熊耳山等地區(qū)NPP值僅為100—300gC/m2,以及西部的羊拱山、岷山、迭山等地區(qū)NPP值也較低(約200—400gC/m2)。時(shí)間尺度上,秦巴山地NPP值隨時(shí)間變化如表1所示,2000—2015年以輕微變動(dòng)、基本不變及中度增長為主,占比達(dá)92.5%。NPP值中度減少和嚴(yán)重減少的區(qū)域主要分布在東部的豫西山地、神農(nóng)架、武當(dāng)山地區(qū),而NPP值增加區(qū)域主要在秦嶺山脈的中段地區(qū)。

圖2 秦巴山地多年年均NPP的分布圖

表1 秦巴山地NPP年變化率分級(jí)統(tǒng)計(jì)

2.2 秦巴山地地表NPP及其對氣候響應(yīng)的緯向分異

為了更加突出表現(xiàn)NPP及其對氣候響應(yīng)的緯向分異,本文選取兩條穿越秦嶺大巴山主要山脊線的典型經(jīng)線107.81°E(西線)和110.53°E(東線),分析秦巴山地NPP的緯向變化,結(jié)果如圖3所示,NPP總的趨勢為：從北向南,隨著緯度的降低而逐漸增加。具體來說,西線從34.13°N降低到31.46°N,NPP值由442gC/m2增加到647gC/m2；東線從34.51°N降低到31.03°N,NPP值由220gC/m2增加到576gC/m2。

圖3 秦巴山地107.81°E、110.53°E剖面線NPP變化圖

沿東西兩條經(jīng)線分別提取NPP與氣溫、降水的相關(guān)系數(shù),其緯向變化如圖4、圖5所示,西線(107.81°E)從北向南,隨著緯度的降低,NPP與降水的相關(guān)性逐漸減弱,NPP與氣溫的關(guān)系由負(fù)相關(guān)逐漸轉(zhuǎn)變?yōu)檎嚓P(guān),分界點(diǎn)位于漢江地區(qū),這與漢江兩側(cè)山地垂直帶結(jié)構(gòu)、類型和分布模式的變化有關(guān)[31]。東線(110.53°E)NPP與氣溫和降水分別呈現(xiàn)弱負(fù)相關(guān)和正相關(guān),但是在神農(nóng)架南坡NPP與降水的關(guān)系突變?yōu)樨?fù)相關(guān),氣溫突變?yōu)檎嚓P(guān)關(guān)系。

圖4 秦巴山地107.81°E降水、氣溫與年均NPP的相關(guān)系數(shù)變化圖

圖5 秦巴山地110.53°E降水、氣溫與NPP相關(guān)系數(shù)變化圖

2.3 秦巴山地NPP及對其氣候響應(yīng)的經(jīng)向分異

為了突出秦巴山地NPP與氣候響應(yīng)的經(jīng)向分異,本文對秦巴山地沿33.96°N和32.67°N兩條緯線的年均NPP與經(jīng)度的關(guān)系進(jìn)行了研究,由圖6可以看出,隨著經(jīng)度的升高,NPP值先增加后減少,在白水江處(105°E)達(dá)到峰值,最低值則出現(xiàn)在極高山區(qū)、秦嶺和大巴山東部。根據(jù)兩條剖面線上NPP與氣溫、降水的相關(guān)系數(shù)分布圖,在大巴山一線(圖7),除大巴山的西部以及東部的十堰市部分地區(qū),NPP與降水呈正相關(guān),與氣溫相關(guān)性較弱,其中85.0%地區(qū)的相關(guān)系數(shù)都在0.2以下。秦嶺一線(圖8)除伏牛山以東和西秦嶺的部分地區(qū)以外,NPP與降水在大部分地區(qū)(78.26%)呈正相關(guān),相關(guān)性先增加后減少,在蟒嶺附近達(dá)到峰值(0.49)。而NPP與氣溫的相關(guān)系數(shù)隨著經(jīng)度升高由正相關(guān)逐漸變?yōu)樨?fù)相關(guān),轉(zhuǎn)折點(diǎn)位于兩當(dāng)縣。

圖6 秦巴山地年均NPP沿32.67°N和33.96°N剖面線的變化

圖7 秦巴山地NPP與降水和氣溫的相關(guān)系數(shù)沿32.67°N的變化

圖8 秦巴山地NPP與降水和氣溫的相關(guān)系數(shù)沿33.96°N 的變化

2.4 秦巴山地NPP及其對氣候響應(yīng)的高程分異

對NPP值按不同海拔進(jìn)行重分類,發(fā)現(xiàn)各海拔分區(qū)的年均NPP值隨高程的增加,呈現(xiàn)先增加后降低的趨勢,NPP均值從453.96gC/m2增加到海拔3000—3500m處達(dá)到最大值,為539.87gC/m2,而后開始逐漸降低,至海拔3500—5528m處達(dá)到最小值,為435.32gC/m2。將秦巴山地各海拔NPP年均變化率(K)值進(jìn)行統(tǒng)計(jì)分析,發(fā)現(xiàn)NPP值在各海拔均呈升高趨勢,并且隨著海拔升高K值先升高后降低,在海拔0—1500m內(nèi)K值為1.87gC m-2a-1—3.01gC m-2a-1,在1500—2000m達(dá)到最大值(5.07gC m-2a-1),2500m以上K值降低至1.50—2.62gC m-2a-1。

表2 秦巴山地NPP均值和變化率隨海拔變化統(tǒng)計(jì)

2.5 秦巴山地NPP及其對氣候響應(yīng)的坡向分異

對秦嶺南北坡和大巴山南北坡年均NPP值隨高度分布的統(tǒng)計(jì)結(jié)果表明,在不同坡向NPP值隨高程的變化存在明顯的差異,如圖9所示。隨高度的增加,大巴山年均NPP值呈現(xiàn)先增加后降低的趨勢。比如,大巴山南坡年均NPP值先由500m以下的530.72gC/m2增加到500—1000m的568.70gC/m2,再降低至2000m以上的450—460gC m-2a-1；大巴山北坡先由470.20gC/m2增加至503.47gC/m2再降低至426.33gC/m2；秦嶺則呈現(xiàn)更復(fù)雜趨勢。秦嶺NPP先由0—500m的373.96—421.71gC/m2降低至1000—1500m的307.84—384.33gC/m2,再增加到437.44—478.01gC/m2,再降低至274.27—422.72gC/m2。

在海拔2000m以下,從大巴山南坡到秦嶺北坡,NPP呈現(xiàn)遞減趨勢,由494.59—568.70gC/m2降低至373.96—307.84gC/m2；2000m以上,秦嶺與大巴山之間的NPP值差異明顯減少,比如2000—2500m,大巴山南北坡分別為454.75、429.32gC/m2,秦嶺南坡和北坡分別為444.91、379.09gC/m2,各坡向之間NPP值較為相近。

為了突出秦巴山地NPP值的坡向分異,將各坡向NPP變化率均值進(jìn)行統(tǒng)計(jì),結(jié)果如圖9所示,2000—2015年秦巴山地各坡向NPP值均呈增加趨勢但存在明顯的差異,大巴山南北坡增長率分別為2.63、1.72gC m-2a-1,秦嶺南北坡增長率分別為3.45gC m-2a-1、2.03gC m-2a-1,秦嶺增長較大巴山更為明顯。

圖9 大巴山和秦嶺的NPP變化率分坡向統(tǒng)計(jì)

在海拔2000m以下,大巴山北坡增長率明顯高于南坡,其中北坡為1.29—3.34gC m-2a-1,南坡為0.12—2.29gC m-2a-1；秦嶺南坡增長率明顯高于北坡,其中南坡為2.58—6.59gC m-2a-1,北坡為1.34—6.32gC m-2a-1。在海拔2000—3000m,大巴山南北坡差異明顯,南坡為1.59—1.85gC m-2a-1,北坡為2.66—3.34gC m-2a-1；但秦嶺南北坡的差異較小,南坡為6.02—6.59gC m-2a-1,北坡為6.32—6.45gC m-2a-1。

對NPP與氣候因子的相關(guān)系數(shù)按海拔進(jìn)行分坡向統(tǒng)計(jì)得到圖10,結(jié)果表明,1000m以下,大巴山地表NPP和氣溫呈現(xiàn)正相關(guān)；在1000—2500m范圍內(nèi),秦嶺地表NPP值與氣溫呈現(xiàn)負(fù)相關(guān),大巴山地區(qū)NPP值和氣溫相關(guān)性較弱。秦嶺和大巴山地表NPP與降水均呈正相關(guān)的關(guān)系,且秦嶺地區(qū)NPP值與降水的相關(guān)性整體強(qiáng)于大巴山地區(qū)。

圖10 秦巴山地年均NPP與氣溫和降水相關(guān)系數(shù)按海拔分坡向統(tǒng)計(jì)圖

3 討論

3.1 秦巴山地地表NPP及對氣候變化響應(yīng)的多維地帶性

秦巴山地地表NPP時(shí)空分布及對氣候變化的響應(yīng)在緯向、經(jīng)向、海拔和坡向等方面具有多維地帶性。2000m以下,由南向北,從大巴山南坡到秦嶺北坡地表NPP呈現(xiàn)逐漸遞減規(guī)律,存在明顯的緯度地帶性。這主要是由于由南向北秦巴山地基帶由大巴山南坡的常綠闊葉林,逐漸過渡到大巴山北坡的常綠—落葉闊葉混交林、秦嶺北坡的落葉闊葉林[32-33],生物量與葉面積指數(shù)的關(guān)系在一定范圍內(nèi)成比例正相關(guān),常綠林NPP值一般大于落葉林NPP值[34]。自西向東隨經(jīng)度的增加,NPP值先增加后減少,在白水江自然保護(hù)區(qū)達(dá)到頂峰,呈現(xiàn)經(jīng)度地帶性規(guī)律。這與白水江自然保護(hù)區(qū)的位置有關(guān),這一區(qū)域位于青藏高原高寒區(qū)、中國東部季風(fēng)濕潤區(qū)和西北干旱區(qū)三大地理分區(qū)的過渡帶,既具有亞熱帶特征帶—亞熱帶常綠闊葉林,又具有青藏高原東緣成分如岷江冷杉為主的暗針葉林及高山灌叢草甸[31,35],加上自然歷史的變遷,該地區(qū)物種極為豐富,使得秦巴山地白水江自然保護(hù)區(qū)是山地垂直帶分帶最多的地區(qū)。

隨海拔增加,秦巴山地整體NPP值和變化率整體呈先增大后減少的規(guī)律,呈現(xiàn)垂直地帶性。這與植被垂直分異有關(guān)。秦巴山地1000m以下的地區(qū)大多被開墾為農(nóng)耕植被帶,植被NPP較少,1000—3500m為山地森林帶,植被NPP增多,到3500m以上以高山灌叢、矮曲林和草甸為主[31],植被稀少,植被NPP減少。2000—2015年,秦巴山地各坡向NPP均呈增加趨勢,但是秦嶺增長較大巴山更為明顯。2000m以下,秦嶺南坡增長率明顯高于北坡,但是大巴山北坡增長率明顯高于南坡,2000—3000m,秦嶺南北坡差異較小,但是大巴山南北坡差異明顯,這與趙芳等[31]提出的大巴山比秦嶺對秦巴山地植被的地帶性分異的影響更強(qiáng)的結(jié)論一致。

3.2 暖溫帶-亞熱帶界線問題

以漢江為界,秦嶺和大巴山無論是NPP均值分布、南北坡對比,還是對氣候因子的響應(yīng)都呈現(xiàn)明顯的差異。中山地區(qū)(1000—2500m),秦嶺NPP與氣溫呈負(fù)相關(guān),而大巴山則呈現(xiàn)正相關(guān)或弱相關(guān)；秦嶺地區(qū)NPP與降水的相關(guān)性整體強(qiáng)于大巴山地區(qū)。這與先前蔣沖[14]、李敏[22]、陳志超[18]等作者在陜南、神農(nóng)架、伏牛山等地的研究結(jié)果對應(yīng)。漢江作為中山地區(qū)NPP與氣溫相關(guān)性由正相關(guān)性到負(fù)相關(guān)的轉(zhuǎn)折點(diǎn),降水由弱相關(guān)到正相關(guān)的轉(zhuǎn)折點(diǎn),是秦巴山地極為重要的界線,即暖溫帶-亞熱帶界線。這就意味著一方面全球變暖將對中山地區(qū)秦嶺植被產(chǎn)生不利影響,而對大巴山有利,另一方面降水的增多則更有助于秦嶺植被的生長。這也暗示了降水是秦嶺植被的限制性因子,降水對秦嶺的影響強(qiáng)于溫度,這與方精云等[36]提出的秦嶺是因?yàn)榻邓牟蛔愣纬傻募俣ǖ哪媳狈纸缇€的結(jié)果一致。這也使得漢江成為氣溫影響的重要界限,漢江以南氣溫變暖有利于植被生長,漢江以北則相反,而降水的增多則對漢江以北植被生長起到較大的促進(jìn)作用。

之前的研究者多使用氣候和植被來界定這一界線問題,而本文則為其提供了一種新的方法、新的佐證。由于地帶性是大尺度問題,先前的研究常使用大尺度指標(biāo),比如水平地帶的植被和氣候的分異來界定,而本文將地表植被NPP值引入秦巴山地氣候響應(yīng)研究,從多維角度分析NPP的時(shí)空格局與氣候變化的關(guān)系,揭示秦巴山地地理和生態(tài)的過渡性和復(fù)雜性,為暖溫帶-亞熱帶分界線的確定提供了一種新的方法和新的論據(jù)。未來還需要結(jié)合南北過渡地帶不同尺度植被、氣候和土壤等分異進(jìn)行系統(tǒng)的研究。

4 結(jié)論

本研究基于GIS與遙感技術(shù),利用MODIS遙感影像及氣象站點(diǎn)數(shù)據(jù),對秦巴山地2000—2015年NPP時(shí)空變化進(jìn)行動(dòng)態(tài)監(jiān)測,利用氣溫、降水?dāng)?shù)據(jù)對秦巴山區(qū)NPP的時(shí)空特征進(jìn)行分析,以此研究秦巴山區(qū)NPP對氣候變化的多維響應(yīng),主要研究結(jié)論如下：

(1)秦巴山地地表NPP存在明顯區(qū)域差異,整體表現(xiàn)為南高北低,中部高東西低。2000—2015年秦巴山地NPP在秦嶺山脈中段的大部分地區(qū)增加明顯,而在豫西山地、神農(nóng)架、武當(dāng)山地區(qū)NPP存在一定程度的減少。

(2)自南向北隨著緯度的增加,秦巴山地地表NPP值呈遞減趨勢,其中2000m以下的中低山地區(qū)表現(xiàn)更為顯著,呈現(xiàn)緯度地帶性。地表NPP與溫度的關(guān)系由負(fù)相關(guān)變?yōu)檎嚓P(guān),轉(zhuǎn)折點(diǎn)出現(xiàn)在漢江。自西向東隨經(jīng)度的增加,NPP值先增加后減少,秦嶺一線地表植被NPP與溫度由正相關(guān)變?yōu)樨?fù)相關(guān),與降水主要呈正相關(guān),相關(guān)性先增加后減少。

(3)隨海拔增加,秦巴山地NPP值和變化率均呈先增大后減少的趨勢,但秦嶺和大巴山不同坡向地表NPP時(shí)空分布及對氣候變化的響應(yīng)呈現(xiàn)明顯的差異。2000—2015年,秦巴山地各坡向NPP均呈增加趨勢,但是秦嶺增長較大巴山更為明顯。2000m以下,秦嶺南坡增長率明顯高于北坡,大巴山北坡增長率明顯高于南坡；2000—3000m,秦嶺南北坡差異較小,但是大巴山差異明顯。中山地區(qū)(1000—2500m),秦嶺NPP與氣溫呈負(fù)相關(guān),而大巴山則呈現(xiàn)正相關(guān)或弱相關(guān)；秦嶺地區(qū)NPP與降水的相關(guān)性整體強(qiáng)于大巴山地區(qū)。