趙一行,高潤藝,李海東,揣小偉,郭曉敏
(1.河北省地礦局國土資源勘查中心,河北 石家莊 050081;2.南京大學(xué)地理與海洋科學(xué)學(xué)院,江蘇 南京 210023;3.生態(tài)環(huán)境部南京環(huán)境科學(xué)研究所,江蘇 南京 210042;4.自然資源部碳中和與國土空間優(yōu)化重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210023;5.中國科學(xué)院地理科學(xué)與資源研究所,北京 100101)
全球變暖對生態(tài)環(huán)境產(chǎn)生一系列負(fù)面影響,低碳減排已成為國際科學(xué)界和各國政府關(guān)注的焦點(diǎn)[1]。由于經(jīng)濟(jì)快速發(fā)展和歷史遺留的能源結(jié)構(gòu)問題等原因,我國碳排放量位居世界首位[2]。2020年9月,國家主席習(xí)近平在第75屆聯(lián)合國大會一般性辯論上宣布中國二氧化碳排放量力爭于2030年前達(dá)到峰值、2060年前實(shí)現(xiàn)碳中和,此后又在多個重大國際場合反復(fù)重申,并強(qiáng)調(diào)要堅決落實(shí)?!吨泄仓醒腙P(guān)于制定國民經(jīng)濟(jì)和社會發(fā)展第十四個五年規(guī)劃和二○三五年遠(yuǎn)景目標(biāo)的建議》將“碳排放達(dá)峰后穩(wěn)中有降”列入中國2035年遠(yuǎn)景目標(biāo)。中央經(jīng)濟(jì)工作會議把做好碳達(dá)峰碳中和工作列為2021年重點(diǎn)任務(wù)之一,對應(yīng)對氣候變化工作做出明確部署[3]。
土地利用變化對碳排放具有重要影響[4]。植被和土壤在調(diào)節(jié)全球碳平衡和緩解氣候變化方面發(fā)揮了重要作用,土地利用變化會導(dǎo)致生態(tài)系統(tǒng)植被和土壤的變化,高生物量用地類型如林地、草地等生態(tài)用地向低生物量用地類型的轉(zhuǎn)變往往伴隨著生態(tài)系統(tǒng)碳儲量的釋放[5]。具體來看,土地利用變化、毀林等活動已經(jīng)導(dǎo)致了陸地生態(tài)系統(tǒng)碳儲量向大氣的釋放。目前,土地利用的碳排放約占全球碳排放總量的10%,是地球系統(tǒng)重要的碳排放來源,而國土空間格局優(yōu)化是推進(jìn)碳中和的重要手段之一[6-10]。20世紀(jì)90年代以來,土地利用變化及其對氣候變化的影響一直是國內(nèi)外的學(xué)術(shù)熱點(diǎn)[11-15]。目前已有多位學(xué)者對此展開研究,如楊慶媛[16]對我國1996—2005年土地利用變化與碳排放的研究表明,土地非農(nóng)化利用和開發(fā)會增加碳排放;赫曉慧等[17]采用灰色預(yù)測模型,動態(tài)評估了我國中原城市群在不同情景下的碳儲量變化趨勢,結(jié)果顯示合理的土地管理政策有助于區(qū)域碳蓄積;GAO等[18]發(fā)現(xiàn),造成長江三角洲城市群碳儲量損失的主要原因是耕地城市化,且在未來4種城市擴(kuò)張方式下碳儲量都將進(jìn)一步損失??傮w而言,土地利用變化制約著生態(tài)系統(tǒng)植被和土壤碳排放,土地利用造成的碳排放占全球碳排放總量的10%,而國土空間格局優(yōu)化是推進(jìn)碳減排的重要手段之一。
以土地利用為切入點(diǎn)開展碳源碳匯效應(yīng)研究,探索低碳用地發(fā)展路徑與策略,將對碳中和的實(shí)現(xiàn)起到積極的推進(jìn)作用,符合國家綠色發(fā)展和生態(tài)文明建設(shè)戰(zhàn)略需求。江陰市作為中國經(jīng)濟(jì)最發(fā)達(dá)的縣級市之一,近年來建設(shè)用地擴(kuò)張劇烈,給當(dāng)?shù)氐耐恋毓?yīng)和區(qū)域低碳發(fā)展帶來了巨大壓力。該研究以江陰市為例,對2000—2015年土地利用變化情況及其對應(yīng)的碳儲量釋放展開研究,并以2015年土地利用碳儲量為基準(zhǔn),為未來江陰市土地利用優(yōu)化方案提出相關(guān)政策建議。
江陰市為江蘇省省轄縣級市,由無錫市代管,位于無錫市北側(cè),介于北緯31°40′34″~31°57′36″,東經(jīng)119°59′~120°34′30″之間,區(qū)位條件優(yōu)越,長期位居全國百強(qiáng)縣前列。2020年江陰市實(shí)現(xiàn)地區(qū)生產(chǎn)總值4 113.75億元,人均地區(qū)生產(chǎn)總值達(dá)24.88萬元。2022年江陰成為我國GDP十強(qiáng)縣,排名第2位。
研究使用的土地利用數(shù)據(jù)來自于中科院地理科學(xué)與自然資源研究所,數(shù)據(jù)精度為30 m,一級土地利用分類的綜合準(zhǔn)確率高于93%,二級土地利用分類的綜合準(zhǔn)確率高于90%。筆者僅使用土地利用一級分類,將土地利用類型分為耕地、林地、草地、水域、建設(shè)用地、未利用地6類。植被碳密度數(shù)據(jù)主要基于筆者所在課題組近年來的800多個研究項目,幾乎涵蓋了所有的植被類型,最終形成了全國植被碳密度圖[19-20],筆者從中截取了江陰市的植被碳密度圖,在AcrGIS 10.5平臺將江陰市植被碳密度圖與土地利用類型圖進(jìn)行交叉分析后,對江陰市各用地類型的植被碳密度進(jìn)行了測算。江蘇省地質(zhì)調(diào)查院和中科院南京土壤研究所提供了江陰市140個不同深度土壤剖面的樣點(diǎn)數(shù)據(jù),包括每個樣點(diǎn)的土壤類型、土壤有機(jī)碳(TOC)含量、容重、經(jīng)緯度、土地利用類型、地表景觀等信息。
由于研究數(shù)據(jù)主要來源于2005年左右,為了使數(shù)據(jù)具有一定的代表性,提高土地結(jié)構(gòu)優(yōu)化的準(zhǔn)確性,故選取2000—2015年作為歷史研究期。
土壤有機(jī)碳密度是指單位面積(1 m2)一定厚度土層中的有機(jī)碳含量,由于排除了面積因素的影響而以土體體積為基礎(chǔ)來計量,土壤碳密度已成為評價土壤中有機(jī)碳的重要指標(biāo)。該研究分別計算表層(0~20 cm)土壤碳密度和碳儲量,具體計算方法如下:
(1)首先計算每個樣點(diǎn)的碳密度。對于具有一定剖面深度(D)的土壤樣點(diǎn)而言,其有機(jī)碳密度(CD)計算公式為
CD=TD×ρ。
(1)
式(1)中,CD為深度D處土壤樣點(diǎn)的碳密度,kg·m-3;TD為深度D處土壤有機(jī)質(zhì)含量,%;ρ為土壤容重,kg·m-3。
(2)利用ArcGIS 10.5軟件將140個土壤樣點(diǎn)信息導(dǎo)入數(shù)字化后的江陰市土壤圖,該圖共包含44個土壤亞類,每個土壤斑塊包含了數(shù)量不等的土壤樣點(diǎn),采取空間插值的方法生成不同時期的全境土壤有機(jī)碳密度圖。
(3)得到碳密度及各類型土壤面積數(shù)據(jù)后,即可求出江陰市總土壤碳儲量(Ctotal),其計算公式為
(2)
式(2)中,Ci為不同土壤亞類的碳密度,kg ·m-3;Si為不同土壤亞類的面積,m2。
由表1可知,2000—2015年,江陰市土地利用發(fā)生較大變化,耕地面積逐年增加,建設(shè)用地面積逐年減少;林地和水域面積變化較小,但總體而言有所增加;未利用地面積幾乎沒有變化;草地面積在2000—2010年期間變化不大,但在2015年有所減少。2000年江陰市耕地占比達(dá)67.25%,為江陰市最主要的土地利用類型,建設(shè)用地次之,但與耕地面積相比差距較大,僅占22.63%,林地、草地、水域等占比很小,總共只占10%左右。2000—2015年江陰市耕地面積占比持續(xù)下降,同時建設(shè)用地占比增加。2015年江陰市主要土地利用類型仍然是耕地和建設(shè)用地,但兩者之間的差距縮小了許多,耕地占比減少到51.71%,建設(shè)用地占比則增加到37.46%。2000—2015年林地、草地、水域以及未利用地占比幾乎沒有變化。
表1 2000—2015年江陰市土地利用類型分布
根據(jù)2000—2015年江陰市土地利用轉(zhuǎn)移矩陣(表2),建設(shè)用地的擴(kuò)張主要是通過占用耕地的方式實(shí)現(xiàn),建設(shè)用地轉(zhuǎn)入面積221.20 km2,其中有149.43 km2來自耕地。從空間分布來看,耕地向建設(shè)用地的轉(zhuǎn)移主要是在原有建設(shè)用地的基礎(chǔ)上向外擴(kuò)張。耕地由于建設(shè)用地的擴(kuò)張大量減少,還有部分耕地轉(zhuǎn)化為林地、水域以及未利用地,轉(zhuǎn)出面積為157.75 km2,轉(zhuǎn)入面積僅為5.84 km2。
表2 2000—2015年江陰市土地利用轉(zhuǎn)移矩陣
從空間分布(圖1)來看,2000—2015年面積增長最為明顯的是建設(shè)用地,主要是在江陰北部以及東部原本的建設(shè)用地基礎(chǔ)上大面積增加,南部和西部原有的建設(shè)用地也向周邊擴(kuò)張。林地和草地僅在江陰中部小面積分布,2000—2010年其空間分布和占地面積基本無變化,但到2015年江陰南部有較明顯的一塊區(qū)域從建設(shè)用地和耕地轉(zhuǎn)化為林地。未利用地在2000—2010年間僅在江陰西部有小面積分布,2015年則在北部和東部的建設(shè)用地中出現(xiàn)了新的未利用地。
對江陰市植被碳密度數(shù)據(jù)及土地利用數(shù)據(jù)進(jìn)行疊加分析,得出江陰市各土地利用類型的植被有機(jī)碳密度。耕地、林地、草地、水域、建設(shè)用地、未利用地的植被有機(jī)碳密度分別為0.55、1.94、0.21、0.35、0.01、0.01 kg·m-2,平均碳密度值為1.79 kg·m-2。各土地利用類型的植被有機(jī)碳密度大小依次為林地>耕地>水域>草地>未利用地>建設(shè)用地。
根據(jù)江陰市土地利用變化轉(zhuǎn)移矩陣以及不同土地利用類型的植被碳儲量特征,測算出江陰市土地利用變化對植被碳儲量的影響。需要說明的是,對每一類土地利用類型來說,這里計算的碳儲量變化值僅代表各類用地增量部分的植被碳儲量,并沒有考慮土地利用轉(zhuǎn)化前后的碳儲量增減(土壤碳儲量計算同理),2000—2015年土地利用變化帶來的植被碳儲量變化見表3。
如果不考慮土地利用轉(zhuǎn)化前后的植被碳儲量增減,由于土地利用變化,2000—2015年江陰市植被有機(jī)碳儲量減少75 009.49 t,土地利用變化造成了植被有機(jī)碳的損失。各土地利用類型轉(zhuǎn)出類型匯總分析結(jié)果顯示,耕地轉(zhuǎn)出類型合計使植被有機(jī)碳儲量減少最多,達(dá)76 963.41 t,其中耕地轉(zhuǎn)變?yōu)榻ㄔO(shè)用地導(dǎo)致植被碳儲量減少80 693.98 t,林地、草地的轉(zhuǎn)出使植被有機(jī)碳儲量分別減少1 206.01、69.30 t。建設(shè)用地、水域的轉(zhuǎn)出使植被有機(jī)碳儲量分別增加3 137.57、24.91 t。其中耕地向建設(shè)用地和水域、林地向建設(shè)用地和未利用地、草地向建設(shè)用地、水域向建設(shè)用地的轉(zhuǎn)化使植被有機(jī)碳減少,其中建設(shè)用地對其他土地利用類型的占用使植被有機(jī)碳損失最為嚴(yán)重,應(yīng)該引起重視。
將江陰市土壤碳密度數(shù)據(jù)在ArcGIS 10.5平臺進(jìn)行克里金插值,得到江陰市土壤碳密度分布圖(圖2)。再利用ArcGIS 10.5軟件將土壤碳密度圖與土地利用類型圖進(jìn)行疊加分析,分區(qū)統(tǒng)計不同土地利用類型表層(0~20 cm)土壤平均碳密度,得到江陰市土壤碳儲量(表4)。
表4 江陰市土壤碳儲量變化
根據(jù)江陰市土地利用轉(zhuǎn)移矩陣以及不同土地利用類型土壤碳密度差異,測算出江陰市土地利用變化對土壤碳儲量的影響(表5)。表5顯示,如果不考慮土地利用轉(zhuǎn)化前后的土壤碳儲量增減,由于土地利用變化,2000—2015年江陰市土壤有機(jī)碳儲量減少7 641.56 t。各土地利用轉(zhuǎn)出類型匯總分析結(jié)果顯示:耕地轉(zhuǎn)出類型合計使土壤有機(jī)碳儲量減少最多,達(dá)8 188.48 t,草地轉(zhuǎn)出類型使土壤有機(jī)碳儲量減少225.28 t,建設(shè)用地轉(zhuǎn)出類型使土壤有機(jī)碳儲量減少60.99 t,未利用地轉(zhuǎn)出類型使土壤有機(jī)碳儲量減少0.47 t。林地、水域的轉(zhuǎn)出合計使土壤有機(jī)碳儲量分別增加36.62、797.03 t。其中耕地向林地、水域、建設(shè)用地,林地向水域、建設(shè)用地,草地向建設(shè)用地,建設(shè)用地向水域,未利用地向建設(shè)用地的轉(zhuǎn)化使土壤有機(jī)碳減少,不利于碳匯的形成,其中建設(shè)用地對耕地的占用使土壤有機(jī)碳損失最為嚴(yán)重,占減少總量的90.50%,應(yīng)該引起重視。
表5 江陰市土地利用變化對土壤碳儲量的影響
基于上述研究結(jié)果,整理2000—2015年江陰市土地利用變化及其碳儲量變化,匯總得到江陰市整體的生態(tài)系統(tǒng)碳儲量變化(表6)。2000—2015年江陰市生態(tài)系統(tǒng)碳儲量整體減少8.3萬t,其中植被碳儲量減少7.5萬t,占生態(tài)系統(tǒng)碳儲量損失的90.75%。耕地是變化最為劇烈的土地利用類型,主要是向建設(shè)用地的轉(zhuǎn)出,這也導(dǎo)致了最大的碳儲量損失,約8.8萬t。雖然僅有0.63 km2的林地發(fā)生了變化,但其導(dǎo)致了0.12萬t的生態(tài)系統(tǒng)碳儲量損失。這表明在未來對土地利用管理的過程中,應(yīng)當(dāng)進(jìn)一步重視對林地資源的保護(hù),控制其轉(zhuǎn)出。同時,草地的轉(zhuǎn)出也導(dǎo)致了生態(tài)系統(tǒng)碳儲量的損失,而其他地類的轉(zhuǎn)出則有助于碳儲量的有效蓄積。
表6 江陰市土地利用變化及其對生態(tài)系統(tǒng)碳儲量的影響
土地利用結(jié)構(gòu)優(yōu)化是指通過土地利用類型、土地管理方式等人為因素改變地表植被類型,間接改變土壤物理結(jié)構(gòu)、土壤微生物、土壤團(tuán)聚體等性質(zhì),從而影響土壤有機(jī)碳的含量和穩(wěn)定性,以期提高土地生態(tài)系統(tǒng)有機(jī)碳儲量。該研究采用線性規(guī)劃方法,建立土地利用結(jié)構(gòu)優(yōu)化模型,并提出優(yōu)化措施。模型將土地利用的碳儲量或者碳排放量作為目標(biāo)函數(shù),將不同土地利用類型的面積作為決策變量,模型的一般表達(dá)公式為
(3)
式(3)中,Xj為決策變量,即土地利用類型面積,km2;aij為約束系數(shù);bj為約束常數(shù)。
土地利用結(jié)構(gòu)優(yōu)化的目的是使土地利用結(jié)構(gòu)趨于合理,該研究選取碳儲量為目標(biāo)效益函數(shù),選擇不同土地利用類型的面積作為決策變量,根據(jù)當(dāng)?shù)氐淖匀缓蜕鐣?jīng)濟(jì)狀況建立約束條件,進(jìn)行土地利用結(jié)構(gòu)優(yōu)化設(shè)計。依據(jù)前文中確定的主要土地利用類型,選取以下6個決策變量:耕地X1、林地X2、草地X3、水域X4、建設(shè)用地X5、未利用地X6。
考慮到生態(tài)環(huán)境建設(shè)的需要,對2030年土地利用結(jié)構(gòu)進(jìn)行重新調(diào)整。
(1)總面積:江陰市為濱港口花園城市,在不發(fā)生行政區(qū)劃變化的前提下,總面積基本保持不變。將江陰市2030年總面積確定為977.60 km2。
(2)耕地:在低碳情景下應(yīng)強(qiáng)化耕地資源保護(hù),實(shí)施嚴(yán)格的耕地保護(hù)政策。一方面建設(shè)用地的外延擴(kuò)展必然占用大量的耕地資源,耕地數(shù)量保持下降趨勢。隨著農(nóng)村人口的下降,農(nóng)村居民點(diǎn)用地需求從理論上是減少的,可以通過土地復(fù)墾將其整理為耕地,所以耕地的減少速度會得到遏制。以基本農(nóng)田保護(hù)為原則并結(jié)合2000年以后的耕地變化速度,將2030年江陰市的耕地保有量確定為351.55 km2,此保有量主要是根據(jù)糧食需求進(jìn)行的預(yù)測,筆者將其作為耕地的下限,將現(xiàn)有的耕地面積503.39 km2作為2030年耕地面積的上限,據(jù)此建立耕地的約束條件方程:351.55≤X1≤505.54。
(3)林地:考慮到生態(tài)環(huán)境建設(shè)要求,森林覆蓋率要不斷提高,2030年林地面積至少應(yīng)不低于現(xiàn)有水平(29.54 km2),將此作為林地面積的下限。根據(jù)2000年以后林地面積的平均變化速度,將增加3.27 km2后的林地面積作為2030年林地面積的上限,即32.82 km2, 據(jù)此建立林地的約束條件方程:29.54≤X2≤32.82。
(4)草地:江陰市的草地面積較少,且變化緩慢。鑒于草地面積已所剩有限以及政府對生態(tài)環(huán)境保護(hù)的重視,未來草地面積減少的趨勢會得到遏制。據(jù)此假設(shè)2030年草地面積會有所增加,面積不低于1.70 km2。但考慮到經(jīng)濟(jì)發(fā)展的需要,且依據(jù)2000年以后草地面積的平均增長速度,預(yù)計2030年草地面積可增加0.34 km2,將增加后的草地面積作為2030年面積的上限,即2.04 km2。據(jù)此建立草地的約束條件方程:1.70≤X3≤2.04。
(5)水域:2000—2015年水域面積發(fā)生少量的變化,期間有增有減,變化不顯著。水域的擴(kuò)張主要是受近年來水產(chǎn)養(yǎng)殖業(yè)發(fā)展的影響,許多草地被挖成坑塘;而水域轉(zhuǎn)為耕地是由于隨著江蘇濱海海岸線的東移,水產(chǎn)養(yǎng)殖用地(尤其是海水養(yǎng)殖用地)轉(zhuǎn)為水田用地。結(jié)合2015年實(shí)際水域面積,預(yù)計2030年水域面積不少于72.58 km2。且根據(jù)水域面積的歷年變化情況,預(yù)計2030年水域面積將增加2.00 km2,將增加后的面積作為上限值74.44 km2,建立水域的約束條件方程:72.58≤X4≤74.58。
(6)建設(shè)用地:未來經(jīng)濟(jì)增長對城鎮(zhèn)用地的需求強(qiáng)勁,城鎮(zhèn)用地將保持快速增長的趨勢,按照2000年以后建設(shè)用地的平均增長趨勢預(yù)測,2030年建設(shè)用地將達(dá)到511.34 km2。按照2005—2010年的增長速度,2030年建設(shè)用地為521.82 km2。據(jù)此建立建設(shè)用地的約束條件方程:511.34≤X5≤521.82。
(7)未利用地:江陰市未利用地面積較小,考慮土地利用的高效集約化趨勢,江陰市2030年未利用地面積會減少,至少會小于現(xiàn)有未利用地面積,據(jù)此建立未利用地的約束條件方程: 0≤X6≤2.05。
根據(jù)對各地類有機(jī)碳密度的計算結(jié)果,建立陸地生態(tài)系統(tǒng)碳儲量最大化的目標(biāo)函數(shù):max(Z)=4.87X1+6.23X2+5.13X3+4.15X4+4.28X5+4.54X6。根據(jù)以上目標(biāo)函數(shù)和各變量的約束條件,采用Lingo 14.0軟件對以上方程進(jìn)行求解,得到2030年江陰市碳儲量最大化的土地利用結(jié)構(gòu)優(yōu)化方案(表7)。
表7 2030年碳儲量最大化的土地利用優(yōu)化方案
由表7可知,在碳儲量最大化的目標(biāo)函數(shù)下,碳儲量最大為445.25萬t。2030年耕地面積減少,林地、草地、建設(shè)用地面積均增加,水域用地面積保持不變,未利用面積減少至0。其中,建設(shè)用地的碳儲量最大,為218.85萬t;其次為耕地,碳儲量為174.78萬t;林地和水域的碳儲量分別為20.45萬和30.12萬t;草地碳儲量最少,為1.05萬t。
將2030年江陰市低碳土地利用結(jié)構(gòu)優(yōu)化方案分別與2015年土地利用結(jié)構(gòu)、自然增長模式下確定的2030年土地利用結(jié)構(gòu)進(jìn)行比較,分析優(yōu)化方案的減排效果(表8)。
表8 基于碳儲量最大化的土地利用結(jié)構(gòu)優(yōu)化減排效果
表8顯示,基于碳儲量最大化的土地利用結(jié)構(gòu)方案能提高江陰市的碳儲量水平,到2030年其碳儲量比自然增長模式下增加0.29萬t,說明優(yōu)化模式下江陰市陸地生態(tài)系統(tǒng)碳儲量明顯增加。
對江陰市2000—2015年的土地利用變化及陸地生態(tài)系統(tǒng)土壤和植被的碳儲量進(jìn)行了估算,同時對未來江陰市土地利用提出了優(yōu)化方案。結(jié)果顯示:
(1) 2000—2015年江陰市16.98%的土地面積發(fā)生變化,其中耕地向建設(shè)用地的轉(zhuǎn)變最為顯著。
(2)江陰市土地利用變化共造成生態(tài)系統(tǒng)碳儲量減少82 651.10 t,主要來自于耕地、林地向外轉(zhuǎn)出,尤其是耕地向建設(shè)用地的轉(zhuǎn)變造成碳儲量損失87 611.00 t。
(3)基于歷史土地利用變化,碳儲量最大化目標(biāo)下的土地結(jié)構(gòu)優(yōu)化方案可使2030年江陰市生態(tài)系統(tǒng)碳儲量較自然增長模式增加0.29萬t。
為進(jìn)一步減少人類活動對碳排放的影響,未來需繼續(xù)重視土地管理從而減少生態(tài)系統(tǒng)碳儲量的釋放,具體可考慮如下對策:
(1)加強(qiáng)自然林地生態(tài)保護(hù),嚴(yán)格控制林地轉(zhuǎn)出,加強(qiáng)農(nóng)田林網(wǎng)建設(shè);同時,針對江陰市草地的下降趨勢,要減少對草地的人為破壞,保護(hù)原生草地。
(2)改善耕作制度,采取合理的農(nóng)作物輪作方式,提高土壤有機(jī)碳含量。
(3)加強(qiáng)對建設(shè)用地的總量控制,確定江陰市各城鎮(zhèn)的蔓延邊界,規(guī)范建設(shè)用地占用耕地的行為,防止建設(shè)用地擴(kuò)張破壞原有自然生態(tài)系統(tǒng)。
(4)江陰市旅游資源種類豐富,獨(dú)占性旅游資源有“徐霞客”“要塞”“華西”等,然而目前優(yōu)勢不明顯,未來應(yīng)進(jìn)一步優(yōu)化產(chǎn)業(yè)結(jié)構(gòu),大力發(fā)展旅游業(yè)。