吳 哲，陳 歆，劉貝貝，初金鳳，彭黎旭

(1.中國熱帶農業(yè)科學院環(huán)境與植物保護研究所，海南?？?571101;2.中國熱帶農業(yè)科學院環(huán)境影響評價與風險分析研究中心，海南?？?571101;3.海南大學環(huán)境與植物保護學院，海南海口 570228;4.河北工程大學城市建設學院，河北邯鄲 056038)

不同土地利用/覆蓋類型下海南島產水量空間分布模擬

吳 哲1，2，3，陳 歆1，2，劉貝貝1，2，初金鳳4，彭黎旭1

(1.中國熱帶農業(yè)科學院環(huán)境與植物保護研究所，海南海口 571101;2.中國熱帶農業(yè)科學院環(huán)境影響評價與風險分析研究中心，海南海口 571101;3.海南大學環(huán)境與植物保護學院，海南?？?570228;4.河北工程大學城市建設學院，河北邯鄲 056038)

設定3種不同的情境，應用InVEST模型，利用水量平衡原理模擬土地利用/覆蓋類型變化下海南島產水量的空間分布。結果表明:海南島年均產水深可達979.04 mm;其中中部地區(qū)產水量最高，西部地區(qū)產水量最低;當林地面積逐漸減少、城市建成區(qū)面積增長1倍時，產水量的增加量最大，年均產水深可達1004.06 mm;當耕地面積逐漸減少、城市和建成區(qū)面積增長1倍時，年均產水深則為1002.22 mm。

產水量;InVEST模型;土地利用/覆蓋;海南島

水資源是人類生活生產中最關鍵的資源之一，水資源的分布直接決定了人口、耕地、企業(yè)等的區(qū)域布局，也會對地區(qū)經濟社會的發(fā)展產生極其重要的影響。因此，研究分析區(qū)域水資源的時空分布與產水量、合理開發(fā)和利用水資源是非常有必要的。然而，產水是一個非常復雜的過程，它與降水、蒸發(fā)、土地利用/覆蓋、土壤滲透和植被蒸騰等因素密切相關，目前應用于產水量分析的水文模型主要有MIKE SHE 模型［1］、TOPMODEL 模型［2］和 SWAT 模型［3］等，這些模型應用范圍廣，模擬精度較高;但是參數(shù)繁多，在實際應用中仍然面臨很多問題。2010年，由斯坦福大學、TNC和世界自然基金會合作聯(lián)合開發(fā)的模型——InVEST(the integrate valuation of ecosystem services and tradeoffs tool)模型［4］被引入中國，該模型對于模型參數(shù)、特性數(shù)據(jù)要求較低;并且輸入數(shù)據(jù)量少，輸出量大;基于GIS平臺，可以將量化的生態(tài)系統(tǒng)服務功能以地圖的形式表達，不再用煩瑣的計算公式、程序語言以及過多文字來表達［5］。InVEST模型擁有多個模塊，采用多級設計，可以進行多尺度、多情景的分析，已被成功應用于夏威夷群島［6］、坦桑尼亞［7］、北京山區(qū)［8-9］、西苕溪［10］、瀾滄江［11］等區(qū)域的生態(tài)系統(tǒng)服務功能及其價值的評估。

海南島水資源豐富，屬熱帶季風性氣候，受東北和西南季風影響，熱帶風暴和臺風頻繁，由于海南島城市排水系統(tǒng)存在諸多問題，內澇嚴重。隨著海南國際旅游島的建設，城區(qū)面積將不斷擴大，也將改變產水量的空間格局。筆者設定3種不同情境，應用InVEST模型產水量模塊分析土地利用/覆蓋類型變化下海南島產水量的空間分布，為城市的健康發(fā)展、水資源區(qū)劃保護和合理開發(fā)提供科學依據(jù)。

1 研究區(qū)概況

海南島是我國僅次于臺灣島的第二大島，地處熱帶，介于北緯 18°10'～ 20°10'、東經 108°37'～111°03'之間，總面積約3.4萬km2;地勢表現(xiàn)為中間高、四周低;森林覆蓋率為 60.5%［12］，排名全國前5;屬典型的熱帶季風海洋性氣候，年均降水量1500 mm，雨量充沛，但時空分布不均。全島有3大河流:南渡江、昌化江、萬泉河，3大河流流域面積約占全島面積的47%。海南島的土壤普查面積約為308.4萬hm2，其中磚紅壤占土地總面積的53.42%;其次是赤紅壤，占總面積的10.01%;其余的自然土壤還有黃壤、燥紅土、新積土、濱海沙土和火山灰土等。

2 模型結構和驅動數(shù)據(jù)

2.1 基本原理

式中:Yxj為第j土地利用類型、柵格x的產水量;Exj為第j土地利用類型、柵格x的每年實際蒸發(fā)量;Px為柵格x的年降雨量。

Exj/Px計算水量平衡中的蒸散部分，計算公式為

式中:ωx為植物可利用水量與年均降雨量的比值;Rxj為第j土地利用類型、柵格x的無量綱的干燥指數(shù);Z為季節(jié)性因子，即表征季節(jié)性降雨特征的常數(shù);Wx為可利用水量;kxj為第j土地利用類型、柵格x的植被蒸散系數(shù);ET0x為柵格x可能蒸散量。

2.2 數(shù)據(jù)來源和處理

InVEST模型產水量模塊需要的驅動數(shù)據(jù)有氣象資料(降水量與可能蒸散量)、土地利用/土地覆蓋類型、土層厚度和土壤有效水量以及生物物理參數(shù)資料等。

a.氣象資料來自于海南島18個長期氣象觀測站1980—2010年30年的觀測數(shù)據(jù)。海南島年均降水量空間分布不均，西部少雨，最小的站點年均降水量約為961 mm;東部多雨，最大的站點年均降水量約為2363 mm。根據(jù)克里金差值法制作出年均降水量柵格圖。由于部分氣象要素的缺失，根據(jù)氣象干旱等級國家標準推薦的Thornthwaite方法進行可能蒸散量ET0的計算:

式中:Ti為i(i=1，2，…，12)月平均溫度，℃;H為年熱量指數(shù);A為常數(shù)。

b.土地利用/土地覆蓋圖來源于國際科學數(shù)據(jù)服務平臺海南島TM遙感影像，經過波段合成、鑲嵌、校正、裁剪和解譯獲得。根據(jù)數(shù)據(jù)分析要求將土地利用/覆蓋類型分為:常綠闊葉林、落葉闊葉林、常綠針葉林、落葉針葉林、混交林、草原、灌叢、濕地、城市建成區(qū)、裸地、水體、耕地等。

c.土壤厚度數(shù)據(jù)來源于中國西部環(huán)境與生態(tài)科學數(shù)據(jù)中心中國1∶100萬土壤圖，通過提取土層厚度的數(shù)據(jù)并用ArcGIS中Lookup模塊制作出土層厚度的柵格圖。海南島的土層厚度為1m左右。結合中國1∶100萬土壤圖與世界土壤數(shù)據(jù)庫中土壤質地的分類，將粉粒、砂粒、黏粒、有機質等質量分數(shù)導入土壤有效水量經驗公式式(10)［15］，提取得到土壤有效水量柵格圖。

式中:Wx為可利用水量;ωsa為砂粒質量分數(shù)，%;ωsi為粉粒質量分數(shù)，%;ωcl為黏粒質量分數(shù)，%;ωOM為有機質質量分數(shù)，%。

車牌識別技術可以準確地識別圖像中的阿拉伯數(shù)字、字母還有漢字，在識別了車輛之后，還能夠將車牌識別的結果傳送到云端的交通綜合系統(tǒng)，交通綜合系統(tǒng)就能夠確定車輛在校園內的停留時間，并提示對車輛進行收費。通過使用該機系統(tǒng)，能夠促進車輛監(jiān)控管理進一步自動化和智能化，提高效率的同時，也能減少人力的投入。

d.生物物理參數(shù)資料反映了土地利用/覆蓋類型的屬性，包括土地利用/覆蓋編碼、作物根系深度、蒸散系數(shù)。作物根系深度數(shù)據(jù)根據(jù)Canadell等［16］對全球范圍內植被最大根深的研究獲得。蒸散系數(shù)則根據(jù)聯(lián)合國糧農組織作物蒸散系數(shù)指南、InVEST用戶指南以及Allen等［17］估算獲得。海南島作物根系深度和蒸散系數(shù)見表1。

表1 海南島作物根系深度和蒸散系數(shù)

2.3 情境的設定

2010年，隨著海南國際旅游島建設正式步入正軌，城市化進程將不斷加快。目前，海南島城市規(guī)模較小，城鎮(zhèn)化率偏低，因此擴大城市的空間框架，提高城市的規(guī)模等級是加快城市化進程的重要步驟。2010年，海南島城市建成區(qū)面積比2000年增長了62.63%。依據(jù)海南土地利用總體規(guī)劃、統(tǒng)計年鑒建成區(qū)的擴展趨勢，并結合近20年中國城市建成區(qū)的發(fā)展格局，預計至2020年，海南島城市建成區(qū)面積將比2010年增長1倍。然而城區(qū)面積的增長必然引起耕地、林地等土地利用類型面積的縮減，從而影響區(qū)域水循環(huán)的過程。本文設定3種不同情境，分析土地利用變化對產水量的影響。

情境1:2010年海南島土地利用/覆蓋狀況;

情境2:假設其他因子不變，林地向城市建成區(qū)轉化，林地面積減少，至2020年各個市(縣)城市建成區(qū)面積翻番;

情境3:假設其他因子不變，耕地向城市建成區(qū)轉化，耕地面積減少，至2020年各個市(縣)城市建成區(qū)面積翻番。

3 研究結果與分析

3.1 模型校驗

調節(jié)季節(jié)性因子，取值范圍為1～10，對產水量模塊的輸出結果進行校驗。根據(jù)歷年海南省水資源公報以及海南省統(tǒng)計年鑒，經過對產水量模塊的反復調試，發(fā)現(xiàn)當季節(jié)性因子取值為7.5時，較接近海南島多年平均產水量。年均產水量為335億m3，折算為多年平均產水深為979.04 mm。

圖1 3種情景下海南島年產水深空間分布

3.2 結果分析

根據(jù)3種情境的設定，海南島年產水深、年產水量空間分布模擬結果如圖1、圖2所示。

根據(jù)圖1、圖2和表2情境1的模擬結果可知，海南島中部地區(qū)產水量最大，可達12239.67m3/hm2，西部地區(qū)產水量最少，最低值僅為為5021.18 m3/hm2，可見海南島產水量的空間分布是極不平衡的。中部地區(qū)雨量充沛，年均降水量可達2 000 mm以上，并且該地區(qū)地勢較高，是多條河流的發(fā)源地。西部地區(qū)處于東南季風的背風坡，年均降水量在1000mm以下，而蒸散量較大，產水量較少。由情境2的模擬結果可見，隨著林地向城市建成區(qū)轉化，產水量將大幅增加。年均產水深可達1004.06mm，最大產水量為12744.99 m3/hm2，最小產水量為 5 392.16 m3/hm2。每個市(縣)的產水量都出現(xiàn)大幅上漲，其中?？谑械漠a水量由情境1的9 117.95 m3/hm2增長至9570.00m3/hm2，三亞市的產水量由8651.52 m3/hm2上升至8 873.95 m3/hm2。由情境3的模擬結果來看，耕地面積減少，城市建成區(qū)面積增長1倍，年均產水深將增長到 1 002.22 mm，最大產水量為12243.12 m3/hm2，最小產水量為5838.43 m3/hm2。

圖2 3種情景下的海南島年產水量空間分布

表2 海南島各市(縣)產水量

由以上分析可知城市建成區(qū)面積的擴展，林地和耕地面積的下降對海南島年均產水量有較大影響。熱帶森林具有很深的根系，能夠吸取深層土壤的水分;且林地土壤的滲透性強，對降雨有較強的攔截作用。耕地較于林地對降雨的攔截能力較弱，根系較淺，因此當耕地轉化為城市建設用地時，產水量的增加量小于林地轉化為城市建設用地的增加量。城市建成區(qū)的地表一般由混凝土、瀝青、水泥等構成，下滲幾乎為零。隨著區(qū)域城市化快速發(fā)展，城區(qū)面積不斷擴大，大面積的天然植被和耕地被住宅、街道、工廠等商業(yè)用地代替，蒸散減少，地下水位降低，不透水面積也將大量增加［18-20］。城市化過程破壞了土壤結構，降低了土壤含水量，因此城市建成區(qū)的產水量較大。

3.3 模擬結果驗證

根據(jù)中國水資源公報、海南省水資源公報、海南省統(tǒng)計年鑒以及海南島各個水文站點實測徑流量等相關數(shù)據(jù)，情景1的模擬結果與實測結果較為一致。各個市(縣)多年產水量模擬值的相對誤差在0.9%～36.1%之間。

4 結論

a.海南島水資源豐富，但時空分布不均，年均產水深可達979.04 mm。中部地區(qū)產水量最多，產水量可達12239.67 m3/hm2，西部地區(qū)產水量最少，產水量僅為5021.18 m3/hm2。

b.當林地面積逐漸減少、城市建成區(qū)面積增長1倍時，產水量的增加量最大，年均產水深可達1004.06 mm，最大產水量和最小產水量分別為12744.99 m3/hm2和5392.16 m3/hm2。

c.當耕地面積逐漸減少，城市建成區(qū)面積增長1倍時，年均產水深增加到1 002.22 mm，最大產水量和最小產水量分別為12 243.12 m3/hm2和5838.43 m3/hm2。

d.不同土地利用/覆蓋方式的轉換，其水文效應有較大的差別。其中林地向城市建成區(qū)轉化，其產水增加量最大，水文效應的敏感性也最高。海南島城市排水系統(tǒng)較為脆弱，內澇頻繁。暴雨期間，多地道路被水淹沒。因此，隨著海南島國際旅游島的建設，城市面積的擴展應將產水量等水文因子考慮在內，對水資源進行合理的規(guī)劃和管理，確定合適的城市規(guī)模和發(fā)展速度，盡量減少對林地的侵占，增加城市綠化面積，從而減少不透水面積的增加，增強城市的防洪抗災能力。

［1］王盛萍，張志強，GE Sun，等.基于MIKESHE模型的潮河流域土地利用與降水變化對水文的影響評價［J］.生態(tài)與農村環(huán)境學報，2012，28(3):320-325.(WANG Shengping，ZHANG Zhiqiang，GE Sun，et al.Assessing hydrologicalimpacts ofchanges in land use and precipitation in Chaohe Watershed using MIKESHE model［J］.Journal of Ecology and Rural Environment，2012，28(3):320-325.(in Chinese))

［2］王潤，劉洪斌，武偉.TOPMODEL模型研究進展與熱點［J］.水土保持研究，2005，12(1):47-48.(WANG Run，LIU Hongbin，WU Wei.Recent advances and research of TOPMODEL ［J］.Research of Soil and Water Conservation，2005，12(1):47-48.(in Chinese))

［3］丁晉利，鄭粉莉.SWAT模型及其應用［J］.水土保持研究，2004，11(4):128-130.(DING Jinli，ZHENG Fenli.SWAT model and its application［J］.Research of Soil and WaterConservation，2004，11(4):128-130.(in Chinese))

［4］KAREIVA P，TALLIS H，RICKETTS T H，et al.Natural capital:theory and practice of mapping ecosystem services［M］.New York:Oxford University Press，2011.

［5］海南省統(tǒng)計局.海南統(tǒng)計年鑒2011［M］.北京:中國統(tǒng)計出版社，2011.

［6］吳哲，陳歆，劉貝貝，等.InVEST模型及其應用的研究進展［J］.熱帶農業(yè)科學，2013，33(4):58-62.(WU Zhe，CHEN Xin，LIU Beibei，et al.Research progress and application of InVEST model［J］.Chinese Journal of Tropical Agriculture，2013，33(4):58-62.(in Chinese))［7］GOLDSTEINA J H，Integrating ecosystem-service tradeoffs into land-use decisions［J］.PNAS，2012，109(19):7565-7570.

［8］FISHER B，TURNER R K，BURGESS N D，et al.Measuring，modeling and mapping ecosystem services in the Eastern Arc Mountains of Tanzania［J］.Progress in Physical Geography，2011，35(5):595-611.

［9］楊芝歌，周彬，余新曉，等.北京山區(qū)生物多樣性分析與碳儲量評估［J］.水土保持通報，2012，32(3):42-46.(YANG Zhige，ZHOU Bin，YU Xinxiao，et al.Biodiversity analysisand carbon storage assessments in Beijing mountainous areas［J］.Bulletin ofSoiland Water Conservation，2012，32(3):42-46.(in Chinese))

［10］周彬，余新曉.基于InVEST模型的北京山區(qū)土壤侵蝕模擬［J］.水土保持研究，2010，17(6):9-13.(ZHOU Bin，YU Xinxiao.Soil erosion simulation in mountain areas of Beijing based on InVEST model［J］.Research of Soil and WaterConservation，2010，17(6):9-13.(in Chinese))

［11］ZHANG Canqiang.Water yield of Xitiaoxi River basin based on InVEST modeling［J］.Journal of Resources and Ecology，2012，3(1):50-54.

［12］CHEN Long，XIE Gaodi，ZHANG Changshun，et al.Modelling ecosystem water supply services across the Lancang River Basin［J］.Journal of Resources and Ecology，2011，2(4):322-327.

［13］DONOHUE R J，RODERICK M L，MC VICAR T R.On the importance ofincluding vegetation dynamics in Budyko’s hydrological model［J］.Hydrology and Earth System Sciences，2007，11(2):983-995.

［14］李斌，李麗娟，覃馭楚，等.基于Budyko假設評估洮兒河流域中上游氣候變化的徑流影響［J］.資源科學，2011，33(1):70-76.(LI Bin，LI Lijuan，QIN Yuchu，et al.Impacts of climate variability on streamflow in the upper and middle reaches of the Taoer River based on the Budyko hypothesis［J］.Resources Science，2011，33(1):70-76.(in Chinese))

［15］ ZHOU W Z，LIU G H，PAN J J，et al.Distribution of available soil water capacity in China［J］.Journal of Geographical Sciences，2005，15(1):3-12.

［16］CANADELL J，JACKSON R B，EHLERINGER J B.Maximum rooting depth of vegetation types at the global scale［J］.Oecologia，1996，108(4):583-595.

［17］ALLEN R G，PEREIRA L S，RAES D，et al.Crop evapotranspiration-guidelines for computing crop water requirements［C］//FAO irrigation and drainage.Rome:FAO-Food and Agriculture Organization of the United Nations，1998.

［18］李麗娟，姜德娟，李九一，等.土地利用/覆被變化的水文效應研究進展［J］.自然資源學報，2007，22(2):211-224.(LI Lijuan，JIANG Dejuan，LI Jiuyi，et al.Advances in hydrological response to land use/land cover change［J］.Journal of Natural Resource，2007，22(2):211-224.(in Chinese))

［19］李曉兵.國際土地利用/土地覆蓋變化的環(huán)境影響研究［J］.地球科學 進展，1999，14(4):395-400.(LI Xiaobing. International research on environmental consequence of land use/cover change［J］.Advance in Earth Sciences，1999，14(4):395-400.(in Chinese))

［20］王艷君，呂宏軍，施雅風，等.城市化流域的土地利用變化對水文過程的影響:以秦淮河流域為例［J］.自然資源學 報，2009，24(1):30-36.(WANG Yanjun，LU Hongjun，SHI Yafeng，et al.Impact of land use changes on hydrological processes in urbanized basin:a case study in Qinhuai River Basin［J］.Journal of Natural Resource，2009，24(1):30-36.(in Chinese))

Simulation of spatial distribution of water yield of Hainan Island with different types of land use/land cover

WU Zhe1，2，3，CHEN Xin1，2，LIU Beibei1，2，CHU Jinfeng4，PENG Lixu1
(1.Environment and Plant Protection Institute，Chinese Academy of Tropical Agricultural Sciences，Haikou 571101，China;2.Environmental Assessment and Risk Analysis Center，Chinese Academy of Tropical Agricultural Sciences，Haikou 571101，China;3.College of Environment and Plant Protection，Hainan University，Haikou 570228，China;4.School of Urban Construction，Hebei University of Engineering，Handan 056038，China)

In this study，the InVEST model was used to simulate the spatial distribution of water yield of Hainan Island under conditions of changing land use/land cover types in three scenarios according to the principle of water balance.The results show that the annual average water yield depth of Hainan Island was 979.04mm.The central area of the island had the highest water yield and the western area had the lowest water yield.When the forest area was gradually reduced and the urban built-up area was doubled，the water yield increased most significantly，with an annual average water yield depth of 1 004.06 mm.When the cultivated area was gradually reduced and the urban built-up area was doubled，the annual average water yield depth reached 1002.22mm.

water yield;InVEST model;land use/land cover;Hainan Island

TV213.4

A

1004-6933(2014)03-0009-05

10.3969/j.issn.1004-6933.2014.03.003

中央級公益性科研院基本科研業(yè)務費專項(1630042013018);中國熱帶農業(yè)科學院本級基本科研業(yè)務費專項(1630042014004)

吳哲(1988—)，男，碩士研究生，研究方向為環(huán)境影響評價與風險評價。E-mail:zhoushanshi@163.com

彭黎旭，研究員。E-mail:Penglixu@hotmail.com

book=13,ebook=95

(收稿日期:2013-10-10 編輯:徐 娟)