張家口市土地利用多功能性動態變化及影響因素*

孫丕苓，許月卿，劉慶果，劉 超

(1.中國農業大學資源與環境學院，北京 100193；2. 曲阜師范大學地理與旅游學院，山東日照 276800；3.黔江新華中學，重慶 409000)

0 引言

土地是人地系統交互耦合的載體，土地利用在人類社會經濟活動中表現出不同的功能特性。土地利用功能(land use function，LUFs)是指不同土地利用方式所提供的私人和公共的產品與服務[1]。土地利用功能的多樣性，即為土地利用多功能性(multi-functionality of land use)，是指一個區域土地利用功能及其經濟、社會和生態功能的狀態和表現，是評價土地利用變化對其功能影響的重要概念和方法體系[2-4]。近年來，隨著城市化與工業化進程的加快，土地資源的稀缺性日益顯現，土地利用方式和土地利用格局發生巨大變化，人地關系日益緊張。與此同時，土地利用的學術研究從關注土地利用格局變化向關注土地利用功能變化轉變[5]，土地利用多功能性研究逐漸成為土地科學研究的重要命題。開展土地利用多功能性研究對提升土地利用效率和實現土地資源可持續利用具有重要現實意義。

自20世紀90年代末以來，國內外學者圍繞土地利用的可持續性影響評估[6-8]、土地利用多功能性的概念與內涵[9-10]、分析框架[3， 11]、識別與分類[12]、度量與評價[5]、變化監測[13]、動態演變[14-15]等方面展開大量研究，也有部分研究成果關注土地利用多功能性的時空格局、空間差異及影響因素研究[16-18]。土地利用多功能性評價作為土地多功能利用研究的基礎和關鍵，旨在度量土地多元化利用所提供的產品、服務和功能，主要通過構建指標體系，運用不同數理分析方法計算土地利用多功能性指數來衡量土地利用功能的多樣化程度。目前，土地利用多功能性領域較為常見的評價方法主要有綜合指數法[18-19]、模糊綜合評價法、改進突變級數法[17]、灰色關聯投影法[15]、全排列多邊形綜合圖示法[20-21]等。其中，綜合指數法操作簡單，在案例研究得到較為廣泛應用[22]。已有的案例研究主要集中在國家、區域、省、市、縣等尺度[23]，評價的目的側重測度土地利用多功能性程度以及各功能之間的協調性[5]，較少涉及對土地利用多功能性影響因素的定量識別。

鑒于此，文章以張家口市為例，基于綜合指數模型對1985—2014年張家口市土地利用多功能性進行綜合評價，并運用灰色關聯模型分析土地利用多功能動態變化的影響因素，以期為該區域土地資源可持續利用提供科學依據。

1 研究區概況與數據來源

1.1 研究區概況

張家口市位于河北省西北部，地處內蒙古高原向華北平原過渡地帶，地理位置介于113°50′E～116°30′ E， 39°30′N～42°10′ N之間，土地總面積為3.68萬km2，是環京津貧困帶的典型區域(圖1)。全市海拔范圍320～2 841km，地勢由西北向東南傾斜，含壩上高原、壩下低山盆地兩個地貌單元，壩上高原區崗梁、湖泊、灘地和草坡相間分布，波狀高原景觀顯著； 壩下低山盆地區地形復雜，山巒起伏，山地、丘陵、河川、盆地相間分布[24]。氣候屬溫帶大陸性季風氣候，多年平均氣溫7.8°C，多年平均降水量350mm。1985—2014年張家口市社會、經濟及生態水平發生明顯變化。至2014年，張家口市總人口468.56萬人，實現生產總值1 358.51億元，全社會固定資產投資完成1 422.61億元，城鎮居民人均可支配收入21 651元，農村居民人均可支配收入7 462元。

圖1 研究區位置及高程圖

1.2 數據來源

該研究涉及的數據主要源于《張家口經濟年鑒(65年特別版)》(2014年)、《張家口經濟年鑒》(2015年)、《河北經濟年鑒》(1986—2015年)、《河北農村統計年鑒》(1995—2015年)的統計指標或根據相關數據整理換算得到。考慮到張家口市行政區劃的變遷，需要進行行政區域調整處理，根據2014年張家口市行政區域范圍對1993年前統計指標進行調整和整理計算。

2 研究方法

2.1 土地利用多功能性評價

2.1.1 指標體系構建

土地利用系統是一個具有特定結構和功能的自然—社會—經濟復合生態系統，土地利用功能是土地利用系統的綜合體現，土地利用功能也稱土地系統功能[5]。該文借鑒FAO《可持續土地利用評價綱要》的框架思路，基于系統論與可持續發展三維度理念，結合張家口市農牧交錯的特征，將土地利用功能劃分為社會、經濟、生態功能3類。參照已有研究成果，將土地利用功能納入以人為主體的邏輯框架，從社會、經濟、生態3個維度選取結構性與效益性相結合的指標構建土地利用多功能性評價指標體系(表1)。(1)社會功能是指土地利用系統提供人類居住空間、就業機會等基本物質保障功能，滿足人類發展需求的能力，主要體現在空間承載功能、社會保障功能和就業支持功能等3項子功能。該文選取人口密度、城鎮人均住房面積、農村人均住房面積表征空間承載功能，以城鄉居民收入平衡指數、農民人均純收入、萬人醫療床位數作為社會保障功能的表征指標，通過從業人口密度、第三產業從業人口比重表征就業支持功能。(2)經濟功能是指土地利用系統提供基本生產資料及經濟基礎設施等的能力，包括農業生產、牧業生產、經濟發展和交通服務等4項子功能。該文選取糧食單產和人均鮮活農副產品產量表征農業生產功能，以人均畜產品產量和地均大牲畜承載量表征牧業生產功能，通過經濟密度、第三產業比重表征經濟發展功能，以路網密度和人均城鎮道路面積表征交通服務功能。(3)生態功能是指土地利用系統維持人類生存的自然條件及其效用的能力，主要包括資源供給功能、環境凈化功能和生態維持功能等3項子功能。該文選取人均耕地面積、人均草地面積表征資源供給功能，以廢污水排放強度、固體廢物產生強度、化肥使用強度和農藥使用強度表征環境凈化功能，通過森林覆蓋率、市區人均公共綠地面積表征生態維持功能。

表1 張家口市土地利用功能指標體系

目標層準則層要素層指標層(單位)標準值標準值依據性質權重土地利用總功能社會功能(0.31)空間承載功能人口密度(人/km2)142全國平均值+0.05 城鎮人均住房面積(m2)32.9國家標準①+0.04 農村人均居住面積(m2)37.1國家標準②+0.04 社會保障功能城鄉居民收入平衡指數(—)0.8理想化標準[15]+0.04 農民人均純收入(元)8000國內公認值+0.03 萬人醫療床位數(張/萬人)48.5全國平均值+0.03 就業支持功能從業人口密度(人/km2)81.5全國平均值+0.04 第三產業從業人口比重(%)40.6全國平均值+0.04 經濟功能(0.33)農業生產功能糧食單產(kg/hm2)6450冀北地區較大值+0.05 人均鮮活農副產品產量③(kg/人)1495冀北地區平均值+0.03 牧業生產功能人均禽畜產品產量(kg/人)80全國平均值+0.04 地均大牲畜承載量(頭/hm2)3.31河北省平均值+0.03 經濟發展功能經濟密度(萬元/km2)713.96全國平均值+0.04 第三產業比重(%)45國內公認值+0.05 交通服務功能路網密度(m/km2)904冀北地區較大值+0.05 人均城鎮道路面積(m2)11.5最佳人居環境標準+0.04 生態功能(0.36)資源供給功能人均耕地面積(hm2)0.08國際公認值+0.05 人均草地面積(hm2)0.33全國平均值+0.04 環境凈化功能廢污水排放強度(萬m3/km2)2參考文獻[15, 25]—0.04 固體廢物產生強度(T/km2)120參考文獻[15, 26]—0.05 化肥使用強度(kg/hm2)225安全使用標準[14, 27]—0.05 農藥使用強度(kg/hm2)2.3安全使用標準[14, 27]—0.05 生態維持功能森林覆蓋率(%)40國際公認值+0.04 市區人均公共綠地面積(m2)18生態城市標準+0.04 注: ①②城鎮、農村人均住房面積標準來源于《中國民生發展報告2012》; ③該研究中鮮活農副產品產量主要是指鮮果、蔬菜產量

2.1.2 綜合測度方法

(1)數據標準化。為消除不同量綱數據對綜合評價的影響，對指標數據進行標準化處理：

(1)

(2)

式(2)中，Zij為各指標標準化值，Xij為指標的實際值，X0j為指標的標準化值。i表示年份，j表示指標個數，下同。標準值根據各指標的特性和水平，參照國際、國內公認的指標標準，國家、地方和行業規定的相關標準或全國與地方的平均水平等途徑進行確定(表1)。

(2)指標權重確定。運用主成分分析法確定指標權重：

(3)

式(3)中，Wj′為未進行歸一化處理的指標權重；gq為第q個主成分對總方差的貢獻率；aqj為指標在第q個主成分中的系數；m為主成分個數。該研究按累計貢獻率達到85%選取主成分個數。歸一化處理后得到各指標權重Wj(表1)：

(4)

(3)評價模型。運用綜合指數模型進行土地利用多功能性評價，計算公式為：

(5)

式(5)中，Fi表示土地利用多功能性指數。運用等間距法對土地利用多功能性水平進行劃分，將其劃分為4個等級，即低水平[0， 0.25]，較低水平(0.25， 0.50]，較高水平(0.50， 0.75]，高水平(0.75， 1]。

2.2 土地利用多功能性影響因素分析

2.2.1 影響因素選取

已有研究表明，土地利用多功能性受自然條件、社會經濟、區域政策等多種因素的共同影響[17-18]，因此，該文選擇年均降雨量、自然災害成災率、城鎮化水平、居民人均可支配收入、人均社會消費品零售額、外商直接投資、科技支出比重、廢污水排放達標率、固體廢物綜合利用率、造林面積等作為影響因子，分析其對土地利用多功能性的影響。

2.2.2 灰色關聯度模型

根據灰色關聯系統理論，將張家口市土地利用總功能與社會、經濟、生態系統視為灰色系統，運用灰色關聯模型衡量指標間的關聯性特征及差異性程度，通過計算灰色關聯度來判定主要影響因素[28]。具體模型為：

x0=x0(k)k=1， 2，…，n=x0(1)，x0(2)，…，x0(n)(k=1， 2，…，n)

(6)

xi=xi(k)k=1， 2，…，n=xi(1)，xi(2)，…，xi(n)(i=1， 2，…，m)

(7)

(8)

式(8)中，δi為比較數列xi對參考數列x0的關聯系數，μ為分辨系數，μ∈[0， 1]，一般取0.5。

3 結果與分析

3.1 土地利用功能變化分析

根據上述公式，計算得到1985—2014年張家口市土地利用總功能及單項功能指數，并繪制其動態變化圖(圖2)以及土地利用總功能中各單項功能的組成比重圖(圖3)。

3.1.1 土地利用總功能

由圖2可知， 1985—2014年張家口市土地利用總功能呈現波動增強趨勢，功能指數介于0.456～0.625之間，土地利用多功能性總體上處于較低水平和較高水平。近30年間張家口市土地利用總功能變化明顯分為兩個階段。

細胞處理結束后，于2000 g/min離心收集細胞，經PBS洗滌后使用TRIzol法提取總RNA，經M-MLV反轉錄酶將mRNA反轉錄為cDNA，然后采用SYBR法進行定量PCR檢測。擴增條件為預變性94℃30 s；循環反應94℃ 10 s，60℃ 30 s，35個循環；冷卻至4℃。PCR引物序列如下：ATG4B上游引物5'-GGTG TGGACAGATGATCTTTGC-3'，下游引物 5'-CCAAC TCCCATTTGCGCTATC-3'。檢測結果采用 2-ΔΔCt法分析后作圖。

第一階段為1985—2000年，張家口市土地利用多功能性主要處于較低水平，土地利用功能以生態功能為主導(圖3)。期間，張家口市土地利用總功能總體上有所減弱，總功能指數介于0.456～0.505。由于歷史方面的影響， 1995年以前張家口市處于相對閉塞狀態，社會經濟發展水平較低，區域發展資源環境依賴性較強，使得土地利用社會、經濟功能在波動中緩慢提升，生態功能則在波動中減弱。1995年張家口市實行對外開放以來，社會經濟得到較快發展，土地利用社會、經濟、生態功能均有所增強但具有波動性。

第二階段為2000—2014年，張家口市土地利用多功能性總體上處于較高水平，土地利用總功能在波動中增強，經歷了波動提升—快速提升的演變過程，總功能指數從0.501提升到0.625。2000—2008年張家口市土地利用多功能性在較低水平與較低水平之間波動，土地利用主導功能由環境功能向經濟功能轉變； 2008—2014年張家口市土地利用多功能性處于較高水平，土地利用功能以經濟功能和社會功能為主導(圖3)。2000年以來，張家口市土地利用社會、經濟功能水平不斷提高，其主導性逐漸增強。尤其是隨著京津冀協同發展和“村村通”工程的開展，張家口市城鎮化進程加快，基礎設施不斷完善，產業結構持續優化升級，土地利用向高效、集約化方向發展，土地利用社會經濟功能逐步提升。同時，退耕還林還草、京津風沙源治理等生態措施的實施，推動了生態環境功能的提升。

3.1.2 土地利用社會功能

1985—2014年張家口市土地利用社會功能緩慢增強，社會功能指數介于0.139～0.247之間，其變化過程經歷了波動提升和持續提升兩個階段(圖2)。1985—2000年張家口市社會功能指數從0.144波動提升到0.176。期間，張家口市空間承載功能與就業支撐功能穩步提升(圖4)。社會保障功能在波動中增強，這一時期城鄉居民收入水平均有所提高，但城鄉差距較大，城鄉居民收入平衡指數較低且處于波動下降狀態。與此同時，張家口市社會經濟發展較為落后，醫療衛生條件的改善落后于人口數量增長速度，萬人醫療床位數從1985年的27.21張/萬人降低到2000年的25.82張/萬人。

2001—2014年張家口市土地利用社會功能持續增強，社會功能指數從0.174增至0.247。期間，空間承載、社會保障、就業支撐功能均不斷增強(圖4)。2000年以來，張家口市人口數量與人口密度不斷增大，城鎮化建設和舊房改造同步進行，城鎮、農村人均住房面積分別從2001年的22.52m2、18.56m2增加到2014年的31.92m2、26.78m2。城鄉居民收入和醫療衛生水平逐漸提高。此外，京津冀協同發展和京張聯合申辦冬奧會為張家口市承接北京和天津的產業轉移提供新機遇，產業結構不斷優化升級。

3.1.3 土地利用經濟功能

1985—2014年張家口市土地利用經濟功能在波動中增強，經濟功能指數從0.067提升到0.232。由于歷史因素、自然條件、生態環境等方面的原因，張家口市土地利用經濟功能較弱。1985—2000年間，張家口市土地利用經濟功能緩慢增強，但處于明顯的波動狀態。2000年以來，張家口市土地利用經濟功能步入快速提升階段(圖2)。圖5顯示，研究期內，張家口市農業生產、牧業生產、經濟發展、交通服務功能均有所增強。(1)1985—2014年，農業生產和經濟發展功能處于明顯的波動狀態，張家口市自然災害頻繁，對其農業生產影響頗大，尤其是干旱、冰雹、大風、霜凍等氣象災害多發且嚴重，近30年來自然災害成災率介于40.03%～98.52%，人均糧食產量在103～363kg之間波動，低于全國400kg/人的平均水平。2002年以來，隨著菜籃子工程的實施，張家口市大力發展錯季蔬菜，致使農業生產功能明顯增強。產業結構的變化對經濟發展影響顯著，第三產業比重的波動變化導致經濟發展功能在波動中增強。(2)1985—2014年，牧業生產功能先在波動中增強后保持平穩， 1985—1995年間，張家口市草地沙化、退化嚴重，草地資源潛能難以發揮，導致其牧業生產功能提升緩慢。1995年張家口市對外開放以來，畜牧業發展迅速。2007年張家口市遭受嚴重的干旱、洪澇、風雹災害，造成牲畜飼草短缺，畜產品產量降低，牧業生產功能出現階段性極小值(0.042)。之后，張家口市以京張聯合申辦冬奧會為契機進行生態管護，通過政策引導推動舍飼養殖，牧業生產功能呈現相對穩定狀態。(3)1985—2014年，交通服務功能經過平穩變化后迅速提升。1985—2000年受地形條件和經濟發展水平限制，張家口市道路建設處于低水平狀態，路網密度在123.34～145.58m/km2之間徘徊，交通服務功能變化平緩； 2000—2014年間，張家口市實施了村村通道路工程，區域交通設施不斷完善，路網密度增大了2.81倍，交通服務功能顯著增強。

圖4 1985—2014年張家口市社會功能子功能變化趨勢 圖5 1985—2014年張家口市經濟功能子功能變化趨勢

圖6 1985—2014年張家口市生態功能子功能變化趨勢

3.1.4 土地利用生態功能

2000—2014年張家口市城鎮化與工業化進程加快，國家采取了天然林保護、退耕還林還草、京津風沙源治理、濕地資源保護等生態措施，導致大量耕地向建設用地、林地轉變，資源供給功能持續減弱，生態維持功能不斷增強。期間，張家口市環境凈化功能先減弱后增強，整體呈減弱趨勢且具有較大波動性。一方面是由于工業化進程的加快和一系列環境保護措施的實施，廢污水排放強度和固體廢物產生強度先增大后降低，其對生態環境產生的壓力先升后降。另一方面張家口市依托氣候優勢進行錯季蔬菜生產，設施農業發展迅速，化肥、農藥使用加劇，土地生態環境不斷惡化。

3.2 土地利用多功能性變化的影響因素分析

基于1985—2014年張家口市土地利用多功能性變化的階段性特征，該研究運用灰色關聯模型測算不同階段張家口市土地利用總功能與各影響因子進行關聯分析，旨在探究不同因素對土地利用多功能性水平的影響程度。鑒于土地利用多功能性影響因子較多，且各因子對土地利用多功能性水平影響程度不同，根據相關研究[27]將關聯度大于0.75的因子確定為關鍵因子，對其進行具體分析。

表2 張家口市土地利用總功能與各影響因子的灰色關聯度

影響因素影響因子1985—2000年2000—2014年1985—2014年關聯度排序關聯度排序關聯度排序自然條件年均降水量0.78030.743100.7467自然災害成災率0.87020.76550.8082社會經濟因素城鎮化水平0.87610.79430.8221居民人均可支配收入0.57290.79240.7506人均社會消費品零售額0.559100.75960.7525外商直接投資0.67750.75080.7139科技支出比重0.60880.75370.70510環境因素廢污水排放達標率0.64970.81910.7923固體廢物綜合利用率0.71940.74990.7448政策因素造林面積0.67660.80720.7764

3.2.1 自然因素

由表2可知， 1985—2000年年均降雨量、自然災害成災率與土地利用總功能的關聯系數分別為0.780、0.870，分列第3、2位，這一時期張家口市社會經濟發展較為落后，其土地利用對自然環境依賴性較強，大部分地區農業生產靠天吃飯，年均降雨量、自然災害成災率對農牧交錯區的土地利用多功能性水平影響顯著。2000年以來，隨著生態建設工程的開展以及溫室大棚、噴灌等現代技術的發展，張家口市土地利用多功能性水平對自然環境的依賴性降低。2000—2014年年均降雨量、自然災害成災率對土地利用多功能性的影響減弱，其關聯度分別降為0.743、0.765，分列第10、5位，表明自然因素對張家口市土地利用多功能性水平具有重要影響，但影響程度有所減弱。綜觀1985—2014年自然災害成災率與土地利用總功能的關聯度均在0.75以上，是影響張家口市土地利用多功能性水平的關鍵影響因子。特殊的地理位置和復雜的地形特點，導致張家口市氣象災害頻發而嚴重，影響其土地多功能利用。

3.2.2 社會經濟因素

表2顯示， 1985—2000年城鎮化發展對土地多功能利用具有重要影響，城鎮化水平與土地利用總功能的關聯系數為0.876，列第1位，這一時期張家口市社會經濟發展水平相對落后，城鄉差距較大，城鎮化水平較高地區土地利用方式較為多樣,土地利用多功能性水平較高。2000—2014年影響張家口市土地多功能利用的社會經濟因素增多，其中城鎮化發展、居民收入與消費水平、區域投資、科技發展等社會經濟因素均對土地利用多功能性水平產生重要影響，城鎮化水平、居民人均可支配收入、人均社會消費品零售額、外商直接投資、科技支出比重與土地利用多功能性水平的關聯系數均達到0.75以上。2000年以來，張家口市城鎮化進程加快，基礎設施的不斷完善，城市經濟輻射帶動作用增強。隨著京津冀協同發展，張家口市承接京津產業梯度轉移，產業結構不斷優化升級，土地利用向高效化、集約化方向發展，土地利用效益持續提高。同時，張家口市鄉鎮企業、鄉村旅游業得到快速發展，吸納了大量農村剩余勞動力，居民收入與消費水平不斷提高，人們對高質量的產品與服務需求增加，土地利用功能趨于多元化方向發展。隨著科技投入的增加，張家口市農業現代化、機械化水平提高，促進了農業生產與經濟發展功能的提升； 廢污水治理和固體廢物綜合利用技術升級，環境凈化功能得到不同程度的提升。此外，這一時期張家口市對外開放程度進一步擴大，外商直接投資對經濟增長的拉動作用增強，并成為影響土地利用多功能性水平的關鍵因子。縱觀整個研究時段，城鎮化水平、居民收入與消費水平均對土地利用多功能性水平具有重要影響。

3.2.3 環境因素

2000年以前張家口市社會經濟發展較為落后，科技水平較低，廢污水與固體廢物處理能力較弱，廢污水排放達標率及固體廢物綜合利用率較低對土地利用多功能性水平影響較小。2000年后，隨著工業化進程的加快，廢污水排放與固體廢物產生強度提高，土地生態環境壓力增大。與此同時，張家口市科技進步為廢污水治理和固體廢物綜合利用提供了技術保障，有助于緩解土地生態系統面臨的壓力，影響土地環境凈化功能。表2顯示， 2000—2014年廢污水排放達標率與土地利用多功能性水平的關聯系數為0.819，固體廢物綜合利用率與土地利用總功能的關聯度為0.749，接近0.75。1985—2014年廢污水排放達標率與固體廢物綜合利用率與土地利用多功能性水平的關聯系數分別為0.792、0.744，表明環境治理與調控對土地多功能利用影響較大。

3.2.4 政策因素

區域政策反映了政府對土地資源配置的宏觀調控，是影響土地利用多功能性水平變化的重要因素。張家口市由昔日的軍事重鎮演變為首都綠色生態屏障，區域生態政策對土地利用多功能性水平產生較大影響。2000—2014年間，國家實施了天然防護林保護、退耕還林還草、京津風沙源治理、荒山造林等生態工程，在一定程度上改善了張家口市的生態環境，促進了生態維持功能的增強。此外，京津冀協同發展和京張聯合申辦冬奧會帶來的重大機遇，帶動了張家口市社會經濟發展，增強了其生態環境功能，進一步提升了土地利用多功能性水平。

綜合上述分析可知， 1985—2014年影響張家口市土地利用多功能性水平的關鍵因子類別、數目及影響程度均處于變化中，總體上表現為由自然、社會經濟因素向環境、政策因素擴展。

4 結論與討論

(1)1985—2014年張家口市土地利用多功能性總體上處于較低水平和較高水平，土地利用總功能變化具有明顯的階段性， 1985—2000年土地利用總功能在波動中減弱，土地利用功能以生態功能為主導； 2000—2014年土地利用總功能在波動中增強，土地利用主導功能由生態功能向經濟、社會功能轉變。

(2)1985—2014年張家口市土地利用單項功能變化趨勢各異，土地利用社會功能緩慢增強，經濟功能在波動中增強，生態功能顯著減弱并處于明顯的波動狀態。1985—2014年張家口市土地利用子功能變化趨勢差異顯著，其中空間承載、就業支撐、生態維持功能不斷增強，社會保障、農業生產、牧業生產、經濟發展功能在波動中增強，交通服務功能先緩慢提升后快速提升，資源供給功能逐漸減弱，環境凈化功能顯著減弱并處于明顯的波動狀態。

(3)1985—2014年影響張家口市土地利用多功能性變化的關鍵因子由自然、社會經濟因素向環境、政策因素擴展。1985—2000年影響土地利用多功能性變化的關鍵因子包括年均降雨量、自然災害成災率、城鎮化水平； 2000—2014年影響土地利用多功能性變化的關鍵因子主要包括自然災害成災率、居民人均可支配收入、人均社會消費品零售額、外商直接投資、科技支出比重、廢污水排放達標率和造林面積。

該文結合張家口市實際情況，基于可持續發展三維度理念構建的評價指標體系，具有可操作性，基本適用于張家口市土地利用多功能性評價。參照國際、國內公認指標標準以及國家、地方和行業規定的相關標準或全國與地方的平均水平等途徑確定指標標準值，采用目標標準化方法進行指標標準化有效避免了分析結果出現單邊上漲的假象。灰色關聯方法被廣泛應用于衡量指標間的關聯性特征及差異性程度，在該研究中應用灰色關聯模型能夠較好地識別土地利用多功能性變化的關鍵因子。同時，該文仍存在不足之處，該研究側重市域尺度的土地利用多功能時序變化分析，土地利用多功能性的空間差異需進一步研究[29]； 鑒于數據的可獲取性與數據統計口徑的一致性，所建立的指標體系基本反映了土地利用多功能性水平，社會保障、交通服務、資源供給功能等方面的指標尚顯不足，有待繼續完善深化。