張神寶，王永芳,，郭恩亮,，包慧娟

（1.內(nèi)蒙古師范大學(xué)地理科學(xué)學(xué)院, 內(nèi)蒙古 呼和浩特 010022；2.內(nèi)蒙古自治區(qū)蒙古高原災(zāi)害與生態(tài)安全重點實驗室,內(nèi)蒙古 呼和浩特 010022；3.蒙古高原氣候變化與區(qū)域響應(yīng)高校重點實驗室, 內(nèi)蒙古 呼和浩特 010022）

植被在陸地生態(tài)系統(tǒng)中具有極其重要的地位，是聯(lián)系土壤、大氣和水分等要素的紐帶，具有明顯的時間變化特征和空間異質(zhì)性[1]。植被覆蓋度(fractional vegetation cover, FVC)是指植被在地面的垂直投影面積占研究區(qū)總面積的百分比，是衡量地表植被覆蓋狀況的一個重要定量信息，也是氣候變化、土地利用變化等因素綜合作用結(jié)果的重要指征[2]。因此，植被覆蓋變化研究也成為區(qū)域生態(tài)環(huán)境效益評估、水土保持動態(tài)監(jiān)測、氣候變化影響研究等領(lǐng)域的重要手段之一[3-4]。

早期的植被覆蓋監(jiān)測方法一般包括調(diào)查采樣法、目測估計法與儀器觀測法[5]。這些方法主觀性強、耗時長，較適用于小區(qū)域尺度的FVC 監(jiān)測。遙感數(shù)據(jù)的出現(xiàn)使得連續(xù)時間序列、大區(qū)域尺度的植被覆蓋度監(jiān)測工作得到迅速發(fā)展[6-7]。但目前仍存在遙感數(shù)據(jù)使用門檻高、數(shù)據(jù)處理過程復(fù)雜、空間分辨率低等問題[8]。谷歌地球引擎平臺(google earth engine, GEE)擁有海量衛(wèi)星影像數(shù)據(jù)與地理數(shù)據(jù)集，是世界上最先進的PB 數(shù)量級地理信息處理分析及可視化平臺，能夠有效解決植被遙感監(jiān)測數(shù)據(jù)處理過程中的難題，在長時序的植被覆蓋變化研究中應(yīng)用廣泛[9]。例如岳奕帆等[10]基于GEE 平臺，以Landsat 遙感數(shù)據(jù)作為數(shù)據(jù)源，對舟曲縣22 年的植被覆蓋度變化進行了分析。Boothroyd 等[11]在GEE平臺支持下，對河岸植被變化、河流平面形態(tài)與瞬時水文的相關(guān)性進行了研究。李晶等[12]基于GEE平臺，利用Landsat 長時序的遙感影像對黃河流域34 年植被覆蓋時空變化進行了監(jiān)測。以上研究均借助GEE 平臺得到了時間連續(xù)、精度較高的FVC監(jiān)測結(jié)果。

奈曼旗地處科爾沁沙地南緣，是我國北方典型的農(nóng)牧交錯區(qū)之一。惡劣的氣候條件加之不合理的人類活動導(dǎo)致該區(qū)土地利用結(jié)構(gòu)發(fā)生巨大變化，致使區(qū)域沙化嚴(yán)重，生態(tài)環(huán)境尤為脆弱[13]。為綜合治理沙化土地，促進區(qū)域生態(tài)環(huán)境質(zhì)量提升，國家與地方政府先后在奈曼旗實施了“三北”防護林、退耕還林草、退牧還草等多個工程措施，并取得了顯著的生態(tài)環(huán)境恢復(fù)效果[14]。植被覆蓋動態(tài)變化及其驅(qū)動因子分析一直是該區(qū)熱點問題之一[15-16]，但由于遙感影像處理過程繁瑣、計算量大等問題，該區(qū)基于較高分辨率遙感影像的長時間序列FVC 動態(tài)變化研究尚少。因此，本研究借助GEE 平臺，利用1991－2020 年Landsat 長時序影像數(shù)據(jù)，結(jié)合同期氣象數(shù)據(jù)，研究奈曼旗FVC 時空演變規(guī)律及其對氣候因子和人類活動的響應(yīng)關(guān)系，以期為奈曼旗生態(tài)修復(fù)治理和高質(zhì)量發(fā)展提供切實有效的科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

奈曼旗(120°19′40″～121°35′40″ E，42°14′40″～43°32′20″ N)位于內(nèi)蒙古自治區(qū)通遼市西南部，區(qū)域 面積約為8 144.44 km2(圖1a)。該區(qū)海拔174～789 m，地勢由南向北逐漸降低，其南部為低山丘陵區(qū)，中部偏南為黃土臺地區(qū)，北部為沙區(qū)[17]。奈曼旗境內(nèi)有西遼河、教來河、孟克河等多條河流，屬于大陸性半干旱的季風(fēng)氣候，冬季寒冷干燥，夏季炎熱多雨，年均氣溫為7.1～7.6 ℃，由東南向西北溫度逐漸增高(圖1c)。年降水量介于332～456 mm，由東南向西北逐漸降低(圖1d)。奈曼旗常年多風(fēng)，以冬季和春季尤其明顯，冬季多西北風(fēng)，春季多西南風(fēng)。栗鈣土為奈曼旗地帶性土壤，但受沙漠化影響，風(fēng)沙土廣泛分布[17]。境內(nèi)植被以山杏(Armeniaca sibirica)、白茅(Imperata cylindrica)、大針茅(Stipa grandis)、羊草(Leymus chinensis)與黃柳(Salix gordejevii)為主[18]。奈曼旗共有14 個鄉(xiāng)鎮(zhèn)單位，2020 年人口為44.47 萬人，產(chǎn)業(yè)結(jié)構(gòu)為“三、一、二”，生產(chǎn)總值為124.53 億元。該區(qū)為典型的農(nóng)牧交錯帶，2020 年耕地和草地面積分別占全旗面積的45.27%和45.65%[19]。

圖1 研究區(qū)位置及氣候因子空間分布圖Figure 1 Location of the study area and spatial distribution of climatic factors

1.2 數(shù)據(jù)來源

通過GEE 平臺(https://code.earthengine.google.com/)獲取了來源于美國地質(zhì)調(diào)查局網(wǎng)站(https://eartxplorer.usgs.gov/) 1991－2020 年 行 列 號 為(30，120)、(30，121)、(31，120)與(31，121)的Landsat-5 TM、Landsat-7 ETM + 、Landsat -8 OLI 植被生長季時間序列數(shù)據(jù)，空間分辨率為30 m、時間分辨率為16 d。借助GEE 平臺篩選出云量低于10%的遙感影像，利用奈曼旗行政界限進行裁剪處理。利用最大值合成算法[20]得到歸一化植被指數(shù)(normalized difference vegetation index，NDVI)年值數(shù)據(jù)集，采用像元二分模型[13]計算得到奈曼旗1991－2020 年FVC 數(shù)據(jù)集；氣象數(shù)據(jù)來源于內(nèi)蒙古自治區(qū)氣象局(http://nm.cma.gov.cn/)，收集了奈曼旗及其周邊共5 個氣象站點的氣象數(shù)據(jù)，采用反距離權(quán)重法(inverse distance weighted，IDW)，對氣象數(shù)據(jù)進行空間插值，獲取分辨率與FVC 相同為30 m 的氣候因子?xùn)鸥駭?shù)據(jù)集。

1.3 研究方法

1.3.1 像元二分模型

本研究采用像元二分模型計算了地表植被覆蓋度[21-22]，公式如下：

式中：FVC為植被覆蓋度，NDVIsoil為裸地或無植被覆蓋地區(qū)的NDVI 值，NDVIveg為純植被覆蓋地區(qū)的NDVI 值。根據(jù)研究區(qū)NDVI 灰度分布情況，以5%和95%置信水平截取的NDVI 上下限閾值作為NDVIsoil和NDVIveg。根據(jù)奈曼旗植被區(qū)實際情況，采用等間距分級法[23]，將FVC 劃分為5 個等級，包括極低覆蓋度(FVC < 0.2)、低覆蓋度(0.2 ≤ FVC <0.4)、中覆蓋度(0.4 ≤ FVC < 0.6)、高覆蓋度(0.6 ≤FVC < 0.8)、極高覆蓋度(FVC ≥ 0.8)。

1.3.2 Sen 氏斜率與Mann-Kendall 趨勢檢驗法

Sen 氏斜率法是Sen 在1968 年提出的計算時間序列數(shù)據(jù)趨勢分析的方法[24]，本研究利用該方法計算了奈曼旗FVC 變化斜率，表征其變化趨勢，計算公式為：

式中：xj和xi為 某一像元在i年 和j年的FVC 值。在本研究中 2020 ≥j≥i≥1991；Median代表所求序列的中位數(shù)，β>0表 示FVC 呈上升趨勢；β<0表示FVC 呈下降趨勢，β=0表示FVC 穩(wěn)定不變。

本研究采用Mann-Kendall 方法在0.05 置信水平上對FVC 的空間變化趨勢進行顯著性檢驗[25]，顯著 性 由Z值 表 示，若|Z| ≥1.96 為 顯 著，|Z| < 1.96 為非顯著。將Sen 氏斜率和Z值計算結(jié)果進行疊加[26]，將FVC 變化趨勢劃分為5 種變化類型(表1)。

表1 植被覆蓋度趨勢變化類型統(tǒng)計Table 1 Statistics of fractional vegetation cover (FVC) trend change types

1.3.3 偏相關(guān)分析

為了剔除FVC 與降水量、氣溫相關(guān)性計算過程中變量間的相互影響[27]，本研究采用偏相關(guān)分析法計算了奈曼旗FVC 與年平均氣溫、年降水量的偏相關(guān)系數(shù)。

式中：rxy、rxz、ryz分別代表x與y、x與z、y與z之間的線性相關(guān)系數(shù)；rxy,z表示剔除z變量后，x與y之間的相關(guān)系數(shù)。本研究x代表FVC，y代表降水量，z代表氣溫。當(dāng)rxy,z大于0，表示兩變量正相關(guān)，rxy,z小于0，表示兩變量負(fù)相關(guān)。本研究同樣采用Mann-Kendall 方法在0.05 置信水平上對相關(guān)性進行了顯著性檢驗，并將輸出結(jié)果劃分為非顯著負(fù)相關(guān)(r< 0，P> 0.05)、顯著負(fù)相關(guān)(r< 0，P≤ 0.05)、顯著正相關(guān)(r> 0，P≤ 0.05)與非顯著正相關(guān)(r> 0，P>0.05)共4 類。

1.3.4 殘差分析

本研究利用殘差分析方法量化了人類活動對奈曼旗FVC 變化的影響[28]。對逐像元FVC 與氣候因子進行回歸分析[29]，得到FVC 的預(yù)測值FVCf，并通過實測的FVCp與預(yù)測的FVCf之間的殘差σ表征人類活動對FVC 的影響。

式中：TEM為年均溫，PRE為年降水量，σ為殘差值。若σ為0，表明人類活動對FVC 無影響；若σ為正，表明人類活動促進FVC 增加；反之則認(rèn)為人類活動造成FVC 降低。

2 結(jié)果與分析

2.1 植被覆蓋時空演變特征分析

2.1.1 植被覆蓋度年際變化特征

研究時段內(nèi)，可將研究區(qū)平均FVC 年際變化分為3 個階段：1991－1998 年，區(qū)域FVC 呈現(xiàn)顯著上升趨勢(P< 0.05)，均值由0.4 上升至0.5，斜率為每年0.007 7；1999－2010 年，F(xiàn)VC 呈非顯著下降趨勢(P> 0.05)，均值介于0.4～0.5，斜率為每年-0.000 6；2011－2020 年，研究區(qū)FVC 均值增至0.5 以上，并呈非顯著上升趨勢(P> 0.05)，斜率為每年0.002 7。整體而言，1991－2020 年奈曼旗FVC 呈非顯著增加趨勢，斜率為每年0.003 1，說明近30 年研究區(qū)植被恢復(fù)與生長狀況良好(圖2)。

圖2 1991－2020 年奈曼旗平均FVC 年際變化Figure 2 Annual variation of mean FVC in Naiman Banner from 1991 to 2020

2.1.2 植被覆蓋度等級變化特征

由奈曼旗4 期FVC 各等級面積及其占比可知，1991 年奈曼旗植被以中覆蓋度及以下等級為主(表2)，其中低和極低覆蓋度等級共占研究區(qū)總面積的46.12%，高覆蓋與極高覆蓋度等級面積占比29.87%，說明奈曼旗將近一半的土地FVC 低于0.4。但從2000 年開始，極低覆蓋度等級面積持續(xù)下降，30 年來共減少317.59 km2，極高覆蓋度等級面積呈連續(xù)增加趨勢，共增加1 218.26 km2。2020 年低、極低覆蓋度等級共占研究區(qū)總面積的40.50%，高覆蓋與極高覆蓋度等級面積占比之和升至44.15%，說明研究區(qū)FVC 明顯提升。

表2 不同等級FVC 面積統(tǒng)計Table 2 FVC area statistics for different grades

根據(jù)1991－2020 年不同F(xiàn)VC 等級的轉(zhuǎn)移矩陣，1991－2020 年奈曼旗植被覆蓋各個等級相互轉(zhuǎn)化(表3)，植被覆蓋“低轉(zhuǎn)高”的面積遠(yuǎn)大于“高轉(zhuǎn)低”，表明近30 年來奈曼旗植被向好的方向轉(zhuǎn)化。FVC 等級降低以相鄰等級的轉(zhuǎn)化為主，中覆蓋轉(zhuǎn)化為低覆蓋面積最多，為398.15 km2，其次為低覆蓋轉(zhuǎn)化為極低覆蓋，面積為361.22 km2，說明FVC 低值區(qū)持續(xù)退化。FVC 等級上升的類型中不僅有相鄰等級間的轉(zhuǎn)化，還包含了FVC 呈兩個及兩個以上等級的提升。例如極低、低和中覆蓋轉(zhuǎn)化為極高覆蓋的FVC 面積分別為130.30、283.29 和523.60 km2，進一步表明30 年間研究區(qū)植被覆蓋度明顯提升。

表3 1991－2020 年FVC 等級的轉(zhuǎn)移面積統(tǒng)計Table 3 Statistics of the transferred area for FVC grades from 1991 to 2020

2.1.3 植被覆蓋度空間變化特征

通過對奈曼旗1991、2000、2010 及2020 年4 期FVC 分級圖進行分析，發(fā)現(xiàn)空間上FVC 分布具有明顯的區(qū)域差異，F(xiàn)VC 總體呈現(xiàn)南北高、東西低的空間分布格局(圖3)。奈曼旗北部與南部FVC 以中覆蓋度為主，東部與西部FVC 以極低、低覆蓋度為主；1991－2020 年，奈曼旗極高覆蓋區(qū)域在北部、南部逐漸擴張，以中、高覆蓋向極高覆蓋轉(zhuǎn)移為主，東部與西部的極低覆蓋與低覆蓋區(qū)域減少。2020 年奈曼旗北部FVC 改善明顯，西部FVC 基本穩(wěn)定不變，南部FVC 出現(xiàn)局部退化?？傮w上，奈曼旗FVC 低等級覆蓋區(qū)域減少，高等級覆蓋區(qū)域增多。

圖3 奈曼旗1991－2020 年FVC 等級變化空間分布Figure 3 Spatial distribution map of FVC grade changes in Naiman Banner from 1991 to 2020

根據(jù)奈曼旗FVC 趨勢變化空間分布圖(圖4)與趨勢變化類型統(tǒng)計表(表4)，F(xiàn)VC 退化區(qū)域占研究區(qū)總面積的29.29%，其中16.95%的區(qū)域顯著退化，在奈曼旗北部和東南部破碎狀分布。穩(wěn)定不變的區(qū)域面積占比12.70%，主要分布在葦蓮蘇鄉(xiāng)。研究時間段內(nèi)，研究區(qū)58.01%的FVC 有所改善，其中38.37%顯著改善，主要分布在奈曼旗中部及其偏南地區(qū)。總體而言，30 年間奈曼旗FVC 空間變化趨勢以改善為主，但北部的FVC 低值區(qū)呈持續(xù)退化的特征。

表4 FVC 趨勢變化類型及其顯著性檢驗結(jié)果統(tǒng)計Table 4 Statistics of FVC trend change types and their significance test results

圖4 1991－2020 年FVC 趨勢變化空間分布Figure 4 Spatial distribution map of FVC trend changes from 1991 to 2020

2.2 植被覆蓋驅(qū)動因子分析

2.2.1 植被覆蓋度對氣候變化的響應(yīng)

30 年年均氣溫與年降水量的時間變化趨勢(圖5)分析表明，奈曼旗年均氣溫呈非顯著上升趨勢(P>0.05)，增加速率為0.014 4 ℃·a-1；年降水量則呈非顯著下降趨勢(P> 0.05)，變化速率為0.543 2 mm·a-1，表明1991－2020 年奈曼旗氣候呈現(xiàn)暖干化趨勢，不利于植被生長與恢復(fù)。然而，在研究時段內(nèi)奈曼旗FVC 有所改善，因此可以認(rèn)為氣候因子對區(qū)域FVC 變化的影響并不大。但考慮到FVC 和氣候因子均存在空間異質(zhì)性，本研究利用偏相關(guān)法進一步分析了降水量、氣溫與FVC 在空間上的相關(guān)性。

FVC 與降水量的偏相關(guān)系數(shù)集中在-0.60～0.60 (圖6)，其中正相關(guān)區(qū)域58.03% (表5)，顯著正相關(guān)面積6.89% (P< 0.05)，在八仙筒、白音他拉、固日班花和青龍山等地分布。FVC 與年降水量負(fù)相關(guān)區(qū)域41.97%，其中顯著負(fù)相關(guān)面積僅3.31% (P<0.05)，零星分布于大沁他拉鎮(zhèn)與義隆永鎮(zhèn)；FVC 與氣溫偏相關(guān)系數(shù)也集中在-0.60～0.60。FVC 與年均溫正相關(guān)區(qū)域48.09%，其中顯著正相關(guān)面積僅2.64% (P< 0.05)，點狀分布于明仁蘇木與六號農(nóng)場。FVC 與氣溫負(fù)相關(guān)區(qū)域51.91%，其中顯著負(fù)相關(guān)面積2.87% (P< 0.05)，條狀分布在葦蓮蘇鄉(xiāng)。整體而言，相關(guān)系數(shù)計算結(jié)果表明氣候因子對FVC 有一定的影響，但結(jié)合顯著性統(tǒng)計結(jié)果與空間耦合度，本研究認(rèn)為氣候因子的影響不是FVC 變化的必然結(jié)果。

表5 奈曼旗FVC 與年降水量、年均氣溫的相關(guān)性及其顯著性結(jié)果分類統(tǒng)計Table 5 Classification statistics of correlation coefficients and significance results of FVC with annual precipitation and annual average temperature in Naiman Banner

圖6 奈曼旗FVC 與年降水量、年均氣溫偏相關(guān)性系數(shù)及其顯著性空間分布Figure 6 Correlation coefficients and significance of spatial distribution of FVC with annual precipitation and average annual temperature in Naiman Banner

2.2.2 植被覆蓋度對人類活動的響應(yīng)

奈曼旗是典型的農(nóng)牧交錯區(qū)，F(xiàn)VC 與人類活動密切相關(guān)。本研究利用殘差分析法研究了FVC 對該區(qū)人類活動的響應(yīng)。由奈曼旗FVC 殘差值空間分布圖(圖7)可以看出，殘差值為正值的區(qū)域占全旗面積的62.59%，與圖4 中FVC 呈改善的區(qū)域分布基本相同，主要分布在奈曼旗中部及偏南地區(qū)；殘差值為負(fù)的區(qū)域占37.41%，分布與圖4 中FVC 呈退化趨勢的區(qū)域的分布較一致?？傮w而言，在研究時段內(nèi)正向的人類活動多于負(fù)向的，且FVC 殘差值空間分布與FVC 空間趨勢變化基本一致。因此，人類活動為奈曼旗植被覆蓋變化的主導(dǎo)因素。

圖7 1991－2020 年奈曼旗FVC 殘差值空間分布Figure 7 Spatial distribution map of FVC residual values in Naiman Banner from 1991 to 2020

3 討論與結(jié)論

3.1 討論

本研究借助GEE 平臺，以Landsat-5 TM、Landsat-7ETM + 和Landsat -8 OLI 為數(shù)據(jù)源，利用像元二分模型計算了奈曼旗30 年長時序高精度的FVC 值，并分析了其時空演變特征。結(jié)果表明，1991－2020 年奈曼旗植被覆蓋經(jīng)歷了先增加、后減少、再增加的演變過程，整體上呈非顯著線性增加趨勢，這與杜佳夢[30]的研究結(jié)果較為一致?？臻g上，奈曼旗FVC表現(xiàn)為南北高、東西低，改善與退化面積占比分別為58.01%和29.29%，這與王旭洋等[16]的研究結(jié)果較為一致，但和崔珍珍等[31]的研究結(jié)果有所差異，可能與選取的遙感影像空間分辨率不一致有關(guān)。

本研究利用偏相關(guān)和殘差分析得人類活動是奈曼旗FVC 變化的主要驅(qū)動力，降水和氣溫則驅(qū)動了部分區(qū)域FVC 的變化，與呂家欣等[32]和Duan 等[33]的研究結(jié)果較為一致。此外，本研究發(fā)現(xiàn)FVC 區(qū)域差異明顯，F(xiàn)VC 退化區(qū)域主要集中在北部地區(qū)。奈曼旗北部為沙區(qū)，分布著廣泛的沙漠化土地，沙漠化敏感性[34]、生態(tài)脆弱性[35]及其發(fā)生災(zāi)害的風(fēng)險等級[36]均較高，且在該區(qū)人類活動也在不斷增強[37]，例如大沁他拉鎮(zhèn)、八仙筒鎮(zhèn)及白音他拉蘇木工業(yè)園區(qū)的設(shè)立使得城鎮(zhèn)擴張建設(shè)用地面積不斷增加[38]。而奈曼旗大面積FVC 的改善與該區(qū)實施的“三北”防護林、尤其是21 世紀(jì)初實施的退耕還林草等生態(tài)環(huán)境工程措施密不可分[33]。本研究發(fā)現(xiàn)奈曼旗FVC 改善區(qū)域與耕地的分布高度重疊[39]，這與Chen等[40]的研究結(jié)果較為一致，中國變綠的過程中，農(nóng)用地貢獻(xiàn)了32%。此外，奈曼旗青龍山鎮(zhèn)地勢高、坡度大，起沙風(fēng)日數(shù)較多[34]，且奈曼旗較其他區(qū)域水資源相對匱乏[41]，不適合大面積的植樹造林。另一方面，多數(shù)人工植被都處于成長初期，管理稍有不當(dāng)便會適得其反。這也導(dǎo)致了青龍山及其周邊地區(qū)出現(xiàn)了FVC 退化現(xiàn)象。

相較于傳統(tǒng)的遙感分析手段，本研究充分發(fā)揮了GEE 平臺在遙感監(jiān)測研究中數(shù)據(jù)存儲量大、運算速度快的特色優(yōu)勢，實現(xiàn)了區(qū)域尺度長時間序列、較高精度的被覆蓋狀況監(jiān)測，可為生態(tài)環(huán)境管理工作提供動態(tài)、客觀和科學(xué)性的基礎(chǔ)數(shù)據(jù)，并豐富區(qū)域植被狀況監(jiān)測技術(shù)手段。本研究也存在不足之處，自然驅(qū)動力僅考慮了氣溫、降水量兩個氣象因素。另一方面，從時空變化趨勢來看，奈曼旗在實施退耕還林(草)等一系列生態(tài)工程實施期間，F(xiàn)VC 呈現(xiàn)明顯的改善趨勢，但由于缺乏驗證數(shù)據(jù)，難以明確不同生態(tài)工程對FVC 影響。因此，今后有必要建立相關(guān)評價指標(biāo)體系定量分析不同生態(tài)工程對FVC的影響，即進一步量化人類活動對FVC 變化的影響。

3.2 結(jié)論

1) 1991－2020 年奈曼旗平均FVC 以每年0.003 1的速率呈非顯著增加趨勢，其中極高FVC 持續(xù)增加。并且，F(xiàn)VC 向高等級轉(zhuǎn)移面積遠(yuǎn)大于轉(zhuǎn)化為低等級面積，說明奈曼旗植被向好的方向轉(zhuǎn)化。

2) 30 年間奈曼旗FVC 呈現(xiàn)南北高、東西低的空間分布格局，區(qū)域差異明顯，改善與退化面積占比分別為58.01%和29.29%，其中北部FVC 低值區(qū)呈現(xiàn)退化趨勢，南部高值區(qū)則呈持續(xù)改善趨勢。

3) 1991－2020 年奈曼旗氣候具有暖干化趨勢，而FVC 呈現(xiàn)非顯著增加趨勢。兩者偏相關(guān)系數(shù)集中于-0.60～0.60，但通過顯著性檢驗的區(qū)域極少。而FVC 殘差值與趨勢變化高度一致，說明人類活動為奈曼旗FVC 變化的主導(dǎo)因素。