陳燕紅,蔡芫鑌,仝 川
1 福州大學(xué)至誠(chéng)學(xué)院, 福州 350002 2 福建師范大學(xué)濕潤(rùn)亞熱帶生態(tài)-地理過程教育部重點(diǎn)實(shí)驗(yàn)室, 福州 350007 3 福建師范大學(xué)地理科學(xué)學(xué)院, 福州 350007 4 福州大學(xué)環(huán)境與資源學(xué)院, 福州 350116
城市化通常表現(xiàn)為城市人口的大量聚集,廣泛的土地開發(fā)活動(dòng)和大幅減少的綠色空間[1- 3]。在中國(guó),2015年的城市化率為56.10%,預(yù)計(jì)城市人口在2030年將增長(zhǎng)70%[4]??焖俪鞘谢瘜?dǎo)致的城市熱島現(xiàn)象,是影響人類健康和生態(tài)系統(tǒng)功能最嚴(yán)重的環(huán)境問題之一[2,5]。位于城市中的綠色空間,包括林/草地、水體和濕地,在維護(hù)城市生態(tài)環(huán)境中扮演著重要的角色[6],其所帶來的生態(tài)效益受到廣泛關(guān)注。許多研究表明城市綠色空間可作為減輕熱島效應(yīng)的有效方法之一。研究提出綠地通過蒸騰作用和遮陰作用能夠有效降低周邊環(huán)境溫度[7- 9];地表水體由于本身的高熱容性和流動(dòng)性,以及較低的熱輻射率,具有降溫、增濕的作用,既是城市中的“冷島”,也對(duì)周邊環(huán)境具有顯著的降溫増濕作用[10- 12];濕地被稱為“地球之腎”、“天然水庫(kù)”和“天然物種庫(kù)”,在調(diào)節(jié)區(qū)域氣候特征方面亦具有重要作用[13]。因此,關(guān)注城市不同綠色空間類型的時(shí)空演化與合理配置,緩解城市熱島效應(yīng),改善人居環(huán)境已然成為當(dāng)前的重要任務(wù)[14- 17]。
相比傳統(tǒng)的大氣溫度測(cè)量技術(shù),遙感和地理信息系統(tǒng)技術(shù)正廣泛應(yīng)用于城市熱島效應(yīng)研究中[18]。遙感分析技術(shù)能直觀、動(dòng)態(tài)地獲取下墊面信息和模擬其未來演化方向,并以其他研究方法無(wú)法比擬的絕對(duì)優(yōu)勢(shì)成為了研究城市冷島強(qiáng)度的重要技術(shù)手段。城市景觀的形成和擴(kuò)張?jiān)黾恿藷釐u的威脅,由其引起的地表溫度時(shí)空變異的定量研究也是多樣化的[19- 21]。土地利用類型/覆蓋和地表溫度之間的關(guān)系在不同的研究條件中是不一致的。多數(shù)學(xué)者對(duì)城市的下墊面、土地利用類型、植被覆蓋度、綠色空間格局、不透水面等影響因子已有了大量研究,卻較少?gòu)木G色空間演化過程與地表溫差方面討論城市冷熱島與城市綠色空間的關(guān)系,對(duì)城市綠色空間不同演化類型引起的溫度效應(yīng)的定量研究也相對(duì)匱乏。因此,如何量化城市綠色空間類型演化對(duì)地表溫度變化的影響仍然是多個(gè)城市應(yīng)對(duì)高溫威脅的關(guān)注熱點(diǎn)問題。本研究選取快速城市化背景下的福州主城區(qū)為研究對(duì)象,利用遙感和地理信息系統(tǒng)軟件平臺(tái)及相關(guān)統(tǒng)計(jì)分析法,著眼于城市綠色空間與地表熱環(huán)境演化特征,開展城市綠色空間演化過程的溫度效應(yīng)定量化研究。研究不僅能有效指導(dǎo)城市綠色空間存量調(diào)控建設(shè),而且還可為相關(guān)城市規(guī)劃設(shè)計(jì)指標(biāo)控制提供理論指導(dǎo)與規(guī)范,具有重要的理論學(xué)術(shù)價(jià)值,為同類研究和規(guī)劃實(shí)踐提供案例借鑒。
福州市位于中國(guó)東南沿海(25°15′—26°39′N,118°08′—120°31′E),是省內(nèi)最大的工業(yè)城市,也是全省政治、經(jīng)濟(jì)、文化中心。近年來,隨著城鎮(zhèn)化進(jìn)程加速和福州城市發(fā)展“東擴(kuò)南進(jìn)、沿江向海”戰(zhàn)略的實(shí)施,福州主城區(qū)下墊面用地類型發(fā)生了巨大的變化。在1996—2005年間福州城區(qū)以年均6%—7%的速度擴(kuò)展[22- 23],主城區(qū)中建筑用地過多,綠地偏少且破碎嚴(yán)重[24],濕地面積不斷減少[25],城市建筑和綠地之間的平均溫差可達(dá)4℃[26],建筑不透水面與地表溫度之間表現(xiàn)出明顯的指數(shù)函數(shù)關(guān)系,高不透水面覆蓋地區(qū)的升溫要明顯快于低不透水面覆蓋地區(qū),城市熱島嚴(yán)重[27- 28]。自2007年起夏熱時(shí)段的平均地表溫度高達(dá)35—40℃,2013年更是成為“四大火爐”之首[29]。本研究選取福州主城區(qū)為研究對(duì)象,主要包含鼓樓區(qū)、臺(tái)江區(qū)、晉安區(qū)和倉(cāng)山區(qū)四區(qū)的建成區(qū)部分,總面積約為257.80km2。
本研究采用數(shù)據(jù)源分別是覆蓋福州主城區(qū)的1993年6月26日、2000年6月29日和2008年7月5日的三景LandsatTM5遙感影像數(shù)據(jù)與2013年8月4日和2016年7月27日的兩景Landsat8遙感影像數(shù)據(jù)。所有影像均來源于中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心,無(wú)云區(qū)能100%覆蓋研究區(qū),影像質(zhì)量較好。圖像預(yù)處理主要包括幾何校正、大氣校正和影像裁剪等。
1.3.1土地利用/覆蓋分類
根據(jù)中國(guó)土地資源分類系統(tǒng)的分類標(biāo)準(zhǔn),結(jié)合福州市實(shí)際土地利用情況,將研究區(qū)的土地利用類型劃分為林/草地、水體、農(nóng)田、濕地、建設(shè)用地和裸地6種類型。本次5個(gè)年份遙感影像總體分類精度分別為:84.08%、81.66%、82.70%、80.62%和80.28%;Kappa系數(shù)分別達(dá)到0.81、0.78、0.79、0.77和0.76;分類等級(jí)良好,分類結(jié)果可信度較高。再結(jié)合國(guó)內(nèi)外研究學(xué)者對(duì)綠色空間的界定[30- 33]與本研究區(qū)域的特點(diǎn),選取下墊面中對(duì)典型的林/草地、水體和濕地為綠色空間研究對(duì)象。
1.3.2綠色空間演化過程的界定
綠色空間的演化過程可分為以下4種情況:綠色空間擴(kuò)張,即非綠色空間類型向綠色空間的轉(zhuǎn)變;綠色空間損失,即綠色空間向非綠色空間的轉(zhuǎn)變;綠色空間交換,即綠色空間之間的相互轉(zhuǎn)換;綠色空間不變,即綠色空間類型未發(fā)生改變。通過ArcGIS 10.2軟件工具對(duì)5個(gè)年份的綠色空間專題圖進(jìn)行疊置分析,可得到1993—2000、2000—2008、2008—2013、2013—2016和1993—2016的5個(gè)時(shí)段的綠色空間動(dòng)態(tài)演化情況(圖1)。
圖1 研究區(qū)綠色空間演化過程變化圖Fig.1 Temporal variations of the urban green space evolution of Fuzhou city
1.3.3地表溫度反演與溫度標(biāo)準(zhǔn)化分級(jí)
在利用Landsat衛(wèi)星遙感數(shù)據(jù)反演地表溫度時(shí),國(guó)內(nèi)外普遍采用的方法主要有熱輻射傳輸方法、單窗算法、劈窗算法和單通道算法,大量研究證明了這些算法的可行性[34- 38]。其中,覃志豪等人提出的單窗算法(2001年)簡(jiǎn)便易行,所涉參數(shù)較易獲取,反演精度能夠滿足要求,故而本研究選用此方法進(jìn)行地表溫度的反演[38]。對(duì)Landsat TM 5和Landsat OLI 8的熱紅外波段反演5景遙感影像的地表溫度,具體步驟如下:
(1)熱紅外波段灰度值(DN值)轉(zhuǎn)化為相應(yīng)的輻射強(qiáng)度值:
Lλ=gain×DN+offset
(1)
式中,Lλ表示傳感器接收到的熱輻射強(qiáng)度,W m-2sr-1μm-1;DN表示熱紅外波段的像元灰度值,即Landsat 5的第6波段和Landsat 8 TIRS載荷的第10波段的亮度值;gain和offset分別為熱紅外波段對(duì)應(yīng)的增益和偏置。
(2)熱輻射強(qiáng)度值轉(zhuǎn)化為亮度溫度:
(2)
式中,Tb為像元的亮度溫度,K;K1和K2為常量,對(duì)于Landsat 5衛(wèi)星,K1=607.76W m-2sr-1μm-1,K2=1260.56K;對(duì)于Landsat 8衛(wèi)星,K1=774.89W m-2sr-1μm-1,K2=1321.08K。
(3)計(jì)算地表溫度(LST):
(3)
C=τ×ε
(4)
D=(1-τ)[1+τ(1-ε)]
(5)
式中,TS為像元的實(shí)際地表溫度,K;a和b為常量,a=-67.355351,b=0.458606;Ta為大氣平均作用溫度,K;τ為大氣透射率;ε為地表比輻射率。
在標(biāo)準(zhǔn)大氣狀態(tài)下(天空晴朗、沒有渦旋作用),大氣平均作用溫度是地面附近氣溫的線性函數(shù),用熱帶平均大氣公式(6)近似地推算Ta[39]:
Ta=17.9769+0.91715T0
(6)
式中,T0為地面附近氣溫,K。
采用NASA提供的大氣校正參數(shù)計(jì)算器和MODTRAN軟件獲取大氣透射率。
采用歸一化植被指數(shù)閾值法[40]估計(jì)地表比輻射率ε:
(7)
(8)
dε=(1-εu)(1-Pν)Fεν
(9)
式中,εv為植被比輻射率,取值0.99;εu為城市表面比輻射率,取值0.92;dε為包括自然表面的幾何分布和內(nèi)部反射的影響;F為形狀因子,取值0.55;NDVImax取0.5,NDVImin取0.2。
(4)溫度標(biāo)準(zhǔn)化分級(jí)
為了消除誤差增加不同時(shí)像溫度數(shù)據(jù)的可比較性,研究采用地表溫度標(biāo)準(zhǔn)化分級(jí)公式(10)對(duì)絕對(duì)溫度進(jìn)行處理。后采用Jenks提出的自然間斷點(diǎn)分級(jí)法[41]對(duì)TN進(jìn)行重分類可得到1993—2016年研究區(qū)地溫等級(jí)分布圖(圖2)。其中,將低溫區(qū)和次低溫區(qū)合并為冷島區(qū)域,次高溫區(qū)、高溫區(qū)和超高溫區(qū)合并為熱島區(qū)域。
(10)
式中,TN為地表溫度標(biāo)準(zhǔn)化值;LST為地表溫度,min及max分別代表其最小值和最大值。
圖2 研究區(qū)溫度等級(jí)分布圖Fig.2 Distribution of different temperature degrades
1.3.4統(tǒng)計(jì)分析
(1)疊置分析
首先分別進(jìn)行兩個(gè)年份土地利用/覆蓋類型的變化檢測(cè)。利用ArcGIS 10.2中的Map Calculator工具將兩個(gè)時(shí)相的土地利用類型分別編號(hào),相同地類編號(hào)相同,再分別將兩時(shí)相的分類結(jié)果圖作如下式疊加計(jì)算,由此可以獲得前后兩時(shí)相土地利用變化情況。其中:i1和i2分別是前后兩個(gè)年份土地利用類型圖中每個(gè)地類的編號(hào)。如數(shù)字1—6分別代表水體、林/草地、農(nóng)田、建設(shè)用地、濕地和裸地;而疊置后的兩位數(shù)如11表示水體不變,12則表示水體演化成林/草地,以此類推。
i1 =i1×10+i2
(11)
(2)重心分析
研究通過ArcGIS 10.2技術(shù)平臺(tái)的Mean Center工具實(shí)現(xiàn)冷/熱島重心、城市重心和綠色空間重心的提取[21,42],從而進(jìn)一步揭示城市不同用地類型的空間演變規(guī)律與地表熱環(huán)境演變趨勢(shì)的關(guān)系。
(3)剖面線分析
分別選取1993年、2000年、2008年和2016年4個(gè)年份的地表溫度等級(jí)和相應(yīng)年份的土地利用/覆蓋類型,以及1993—2016年的地表溫差和對(duì)應(yīng)的土地利用/覆蓋變化作為分析對(duì)象,提取剖面沿線的地類與地表熱環(huán)境特征變化分析,以更有效揭示城市沿線不同用地類型引起的溫度效應(yīng)。
對(duì)各年份土地利用類型進(jìn)行疊置分析,得到城市綠色空間在不同研究階段的演化情況(圖1、表1)。從面積變化大小來看,各研究階段中,綠色空間不變面積>綠色空間損失面積>綠色空間擴(kuò)張面積>綠色空間交換面積;1993—2016年,城市綠色空間不變、損失、擴(kuò)張和交換的面積演化速率分別為1.66、2.02、0.21、0.10km2/a。結(jié)合1993—2016年的土地利用轉(zhuǎn)移情況表(表2)發(fā)現(xiàn),綠色空間損失部分主要轉(zhuǎn)化為建設(shè)用地,其中尤以林/草地?fù)p失最多,表現(xiàn)為:林/草地(16.80km2)>水體(11.74km2)>濕地(1.76km2)。演化過程中,綠色空間擴(kuò)張與損失的面積差均為負(fù)值,且在2000—2008年這一階段相差最多,為22.21km2??臻g分布上看,面積差值最大的區(qū)域?yàn)閭}(cāng)山區(qū)(-33.23km2),其次是晉安區(qū)(-6.03km2)和鼓樓區(qū)(-1.96km2),最小的是臺(tái)江區(qū)(-0.55km2)。
表1 研究區(qū)1993—2016年不同綠色空間演化過程面積變化統(tǒng)計(jì)結(jié)果
表2 研究區(qū)1993—2016年土地利用轉(zhuǎn)移變化面積統(tǒng)計(jì)結(jié)果/km2
從面積大小上來看(表3),各研究階段中,不變的冷島區(qū)>不變的熱島區(qū)> 由冷島區(qū)轉(zhuǎn)為熱島區(qū)>由熱島區(qū)轉(zhuǎn)為冷島區(qū)。此外,不變的冷島區(qū)面積波動(dòng)較大,呈現(xiàn)出先下降后上升的趨勢(shì),而不變的熱島區(qū)的面積變化趨勢(shì)則相反。再者,由冷島區(qū)向熱島區(qū)轉(zhuǎn)化的面積波動(dòng)范圍為3.45—19.72km2,由熱島區(qū)向冷島區(qū)轉(zhuǎn)化的面積則較少,范圍在0.21—1.93km2之間。
1993—2016年,冷島不變區(qū)、熱島不變區(qū)、冷島損失區(qū)和熱島損失區(qū)的年變化量分別達(dá)到5.33、0.49、1.47、0.12km2/a。從溫度等級(jí)轉(zhuǎn)移情況上來看(表4),不同的溫度分布區(qū)有較大變化,冷島區(qū)面積大幅減少,熱島區(qū)面積顯著增加。低溫區(qū)與次低溫區(qū)面積損失分別為25.41km2和40.40km2,年變化率為1.10km2/a和1.76km2/a。其中,次高溫區(qū)面積增加最為顯著,其次是高溫區(qū)和超高溫區(qū)。熱島區(qū)從1993年開始逐漸上升,至2008年達(dá)到最大值,后又下降,至2016年又有所回升。熱島區(qū)面積增加最為明顯的是倉(cāng)山區(qū),其次是晉安區(qū)和鼓樓區(qū);而臺(tái)江區(qū)的熱島區(qū)面積則有所減少。這與綠色空間損失與擴(kuò)張的空間分布情況一致。因此,城市冷/熱島區(qū)分布變化與下墊面土地利用/覆蓋變化情況息息相關(guān)。
表3 研究區(qū)各階段冷/熱島斑塊轉(zhuǎn)移情況統(tǒng)計(jì)結(jié)果
UCI: 冷島區(qū)Urban cold island; UHI: 熱島區(qū)Urban heat island
表4 1993—2016年研究區(qū)溫度等級(jí)轉(zhuǎn)移情況統(tǒng)計(jì)結(jié)果/km2
LST: 地表溫度Land surface temperature
2.3.1基于空間重心法的城市冷/熱島演化
通過空間重心法分析城市冷/熱島重心與城市重心的遷移變化情況(圖3)發(fā)現(xiàn),2000年之前,城市重心和熱島重心轉(zhuǎn)移方向不完全一致;2000年以后,城市熱島重心與城市重心偏移方向基本吻合,且逐漸聚集;而城市冷島重心則越來越偏離熱島和城市重心。結(jié)合不同綠色空間類型的重心位置遷移情況(圖3)發(fā)現(xiàn),林/草地、水體和濕地3種綠色空間類型均發(fā)生了較大的位移,位移最大的是濕地,表明其受到的干擾最為強(qiáng)烈;其次是水體,最后是林/草地。這說明城市熱島重心與城市重心轉(zhuǎn)移方向有關(guān),而城市冷島的重心與城市綠色空間分布關(guān)系密切。城市下墊面土地利用/覆蓋類型的變化是造成城市地表溫度擾動(dòng)的重要因素。
圖3 1993—2016研究區(qū)綠色空間、城市重心和冷/熱重心位置變化圖Fig.3 Movement track of the gravity center of UGS, city and cold/heat island from 1993 to 2016圖中的X、Y軸分別為地理坐標(biāo)(WGS84)轉(zhuǎn)換為工程坐標(biāo)(Beijing54)后的橫、縱坐標(biāo)值
2.3.2城市下墊面土地利用變化的溫度效應(yīng)
從不同溫度等級(jí)與用地類型分布變化情況來看(圖4),1993年和2000年,研究區(qū)以低溫區(qū)和次低溫區(qū)為主,分布著大量的城市綠色空間和較少的建設(shè)用地類型;以建設(shè)用地為主的次高溫、高溫和超高溫區(qū)分布較少。2008年后,大量的綠色空間被建設(shè)用地占用,次高溫區(qū)和高溫區(qū)面積顯著擴(kuò)大,在浦下洲濕地范圍甚至出現(xiàn)了超高溫區(qū)。2016年,超高溫區(qū)消失,但仍有較多的高溫區(qū)分布。
從地表溫差與用地變化關(guān)系上看出(圖5),1993—2016年,福州主城區(qū)地表溫差呈峰谷交替態(tài)勢(shì)出現(xiàn)。峰值主要由綠色空間轉(zhuǎn)化為建設(shè)用地或裸地的綠色空間損失過程引起的,其溫差高于平均地表溫差;波谷位置多由綠色空間不變過程引起,溫差低于平均地表溫差;此外,由建設(shè)用地或裸地轉(zhuǎn)化為綠色空間的擴(kuò)張過程所引起的地表溫度降溫較小,這主要由于綠色空間演化方式并非是成片擴(kuò)張,而是在建設(shè)用地中零星式的擴(kuò)散,其引起的溫差必然受到周邊用地類型的影響。定量研究不同綠色空間演化過程引起的溫度效應(yīng)尤為重要。
圖4 1993—2016剖面線不同溫度等級(jí)與用地類型分布變化圖Fig.4 Profile analysis of the relationship between LST grade and LULC from 1993 to 2016溫度等級(jí)1—6分別代表低溫、次低溫、中溫、次高溫、高溫和超高溫;用地類型1—6分別代表水體、林/草地、農(nóng)田、建設(shè)用地、濕地和裸地;LST(land surface temperature),表示地表溫度;LULC(land use/land cover),表示土地利用/覆蓋
圖5 研究區(qū)1993—2016年地表溫差與用地變化關(guān)系剖面圖Fig.5 Profile analysis of the relationship between LST difference and evolution type from 1993 to 2016
2.3.3不同綠色空間演化過程的溫度效應(yīng)
圖6 1993—2016年不同綠色空間演化類型溫度效應(yīng)圖 Fig.6 Temperature effect of different types of UGS evolution from 1993 to 2016
采用隨機(jī)采樣的方式,對(duì)1993—2016年間,研究區(qū)100組不同綠色空間演化類型及相應(yīng)的地表溫度變化數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。從結(jié)果上看(圖6),不同綠色空間演化類型引起的降溫效應(yīng)依次排序?yàn)椋壕G色空間擴(kuò)張>綠色空間不變>綠色空間交換>綠色空間損失。值得注意的是,綠色空間損失引起的平均地表溫差約為+7.0℃,表現(xiàn)為升溫作用;而綠色空間擴(kuò)張帶來的平均地表溫差約為-5.0℃,表現(xiàn)為降溫作用,二者產(chǎn)生的平均升降溫作用并不相同。因此,由用地類型轉(zhuǎn)變帶來的溫度變化并非簡(jiǎn)單的“對(duì)等”關(guān)系。這是由于在快速的城市化進(jìn)程中,對(duì)綠色空間的擠占通常是成片且較為聚集的,這很大程度上削弱了綠色空間的冷島效果。相反,綠色空間的擴(kuò)張往往是通過邊緣式或飛地式的方式實(shí)現(xiàn),周邊環(huán)境的變化使其并不能完全抵消先前由于綠色空間的異地?fù)p失所帶來的升溫??傊?城市綠色空間損失和擴(kuò)張之間存在較大的不均衡性,由此引起了地表溫差的不均勻分布。而針對(duì)城市下墊面土地利用/覆蓋類型的優(yōu)化調(diào)控,尤其是對(duì)綠色空間存量的科學(xué)調(diào)配,值得學(xué)者們進(jìn)一步研究。
本研究以福州市主城區(qū)綠色空間為研究對(duì)象,利用疊置分析法、空間重心法和剖面分析等方法,定量分析福州城市綠色空間演化過程及其引起的地表溫度變化。結(jié)果表明:(1)1993—2016年,福州主城區(qū)下墊面土地利用/覆蓋變化顯著,城市建設(shè)用地面積大幅增加,綠色空間數(shù)量明顯減少,城市地表升溫明顯。(2)城市下墊面地表溫度等級(jí)空間分布變化較大,冷島區(qū)面積大幅減少,熱島區(qū)面積顯著增加。(3)城市熱島重心與城市重心的遷移方向有關(guān);城市冷島重心則與城市綠色空間的分布關(guān)系密切。(4)城市綠色空間的不同演化過程所引起的地表溫度變化不盡相同,綠色空間擴(kuò)張和綠色空間損失分別有著最大的降溫和升溫效果。
因此,采用多方法耦合分析不同綠色空間演化類型引起的地表溫差,可定量出不同綠色空間演化類型引起的地表溫度變化。本研究以獨(dú)特的視覺和方法,表征了城市下墊面不同溫度等級(jí)和土地利用/覆蓋類型的空間分布變化及關(guān)系,有效識(shí)別了不同綠色空間演化過程所引起的地表溫度變化差異,為城市綠色空間健康發(fā)展提供了嶄新思路和科學(xué)依據(jù)。
綜上所述,城市綠色空間作為城市的“綠脈”,在緩解城市熱島效應(yīng)方面具有重要作用??焖俚某鞘谢l(fā)展過程中,下墊面的土地利用/覆蓋變化是引起福州主城區(qū)地表溫度變化的重要因素。而城市綠色空間損失的日益加劇與不合理的綠色空間布局?jǐn)U張方式,一定程度上加劇了城市熱島效應(yīng)。在未來,應(yīng)進(jìn)一步將城市綠色空間統(tǒng)籌進(jìn)城市復(fù)合生態(tài)系統(tǒng)加以規(guī)劃和管理,著力提升城市下墊面用地類型在質(zhì)量和數(shù)量方面的配置水平,以有效緩解城市熱島效應(yīng),為綠色城市可持續(xù)、健康發(fā)展提供支持。
致謝:本研究是在福州大學(xué)至誠(chéng)學(xué)院、福建師范大學(xué)濕潤(rùn)亞熱帶生態(tài)-地理過程教育部重點(diǎn)實(shí)驗(yàn)室、福建師范大學(xué)地理科學(xué)學(xué)院和福州大學(xué)環(huán)境與資源學(xué)院的共同合作下完成。