基于TM資料對上海區(qū)域熱環(huán)境特征的分析

王桂玲，談建國，崔林麗，王成林，金賽花

(1.解放軍理工大學氣象學院，江蘇 南京 211101;2.上海市氣象科學研究所，上海 200030;3.上海市衛(wèi)星遙感與測量應用中心，上海 201100)

0 引言

城市化導致自然植被、土壤和水體被改變，必然造成地表物理特征和地氣之間能量交換的改變。地表熱通量與地表溫度、土壤濕度及植被狀況之間關系最為直接也較為復雜［1］。由于城市地表呈現尺度小、類型多的特點，各類地表又有其各自特性，它們的輻射、熱力及水汽分布各不相同，這就使得城市地區(qū)的觀測和模式研究都很困難，而目前衛(wèi)星遙感技術不僅能有效探測城市下墊面的溫度特征，而且還能獲取地表反照率及植被指數等地表特征，因此得到了越來越廣泛的應用［2－4］。目前，利用遙感資料從地表通量的角度定量地研究城市熱島形成機制，尤其城市中綠地和水體的生態(tài)效應問題還比較少，因此利用遙感資料進行城市地表熱環(huán)境的研究對減緩城市熱島效應也具有重要意義。

利用衛(wèi)星資料對地表熱通量的估算通常依據對地表輻射溫度處理不同分為單源模式和雙源模式［5－7］。盡管雙源模式物理意義比較合理，但公式中經驗系數較多且難以確定，導致誤差較大。地表能量平衡參數化(SEBAL)，是基于衛(wèi)星遙感測量的一種較新的地表熱通量參數化方法，由BRUTSAERT等提出的單源模式基礎上發(fā)展而來，該方法相對簡單，計算中除了風、氣溫外不需要其他資料，一直被人們應用并獲得了較好的結果［8－10］。

為了揭示城市熱島(UHI)形成的內在物理機制及其發(fā)生原因，文中以TM衛(wèi)星遙感資料為主要數據源，結合地面自動氣象站(AWS)實測氣象資料，利用SEBAL的方法估算了上海地區(qū)的地表凈輻射通量和土壤熱通量，分析了上海區(qū)域的地表溫度及地表凈輻射通量和土壤熱通量的空間分布及變化特征。

1 研究范圍及資料來源

1.1 研究區(qū)域概況

研究范圍是以上海(31.25°N，121.47°E)為中心，面積為157 km×150 km，包括了整個上海城區(qū)及其相鄰區(qū)域，其地理范圍為 30.57～31.93°N，120.76～122.17°E。

1.2 資料來源

衛(wèi)星遙感資料來源于2007年7月28日的TM衛(wèi)星地表特征資料，包括地表溫度、地表反照率、植被指數及葉面積指數和土地覆蓋類型資料。由于氣象資料分辨率較粗，在計算通量過程中將TM衛(wèi)星資料由30 m分辨率合成為500 m分辨率。在研究區(qū)域內共有83個地面自動氣象站地面觀測資料可以利用(表1)。

表1 地面自動站及其地理位置

續(xù)表1

1.3 資料分析

TM地表覆蓋類型資料共有10種地表類型，圖1與表2為研究區(qū)域的地表覆蓋類型及說明。

圖1 上海地區(qū)的地表覆蓋類型

表2 TM地表覆蓋類型及其代碼

從圖中可以看出，該地區(qū)地表狀況復雜，整個地區(qū)的下墊面是極不均勻的，有水面、公園、林地、綠地、交通、商業(yè)區(qū)及居民區(qū)等多種類型，在城市中心區(qū)域，商業(yè)區(qū)和居民區(qū)以及交通面積所占比例較大。文中以商業(yè)區(qū)和居民區(qū)建筑以及交通面積的平均值代表中心城區(qū)特征值，而綠地面積平均值代表綠地特征值。

將自動站白天11時地表溫度與TM地表溫度進行對照，根據各個站點的地理位置確定其在研究區(qū)域中對應的位置，然后對TM地表溫度進行訂正(圖2)。由圖可見，在中心區(qū)域存在明顯的UHI現象，中心城區(qū)地表溫度比周圍地區(qū)明顯增高。圖3表示不同地表類型下的TM地表溫度分布狀況，可見不同類型地表，其地表溫度差異明顯，交通道路地表溫度最高，居民區(qū)和商業(yè)區(qū)地表溫度其次，水面最低。交通道路與水面地表溫差可達10℃，與綠地和林地相差約為4～5℃，由此可見城市道路是城市重要熱量來源。

2 地表熱通量的計算

地氣系統(tǒng)主要驅動力是輻射強迫，地表熱通量的確定對于理解地氣能量交換過程有著重要的作用。地表能量來源是凈輻射，其平衡方程如下:

Rn=H+LE+G

式中:Rn——地表凈輻射通量，W/m2;H——感熱通量;LE——潛熱通量;G——土壤熱通量。

2.1 凈輻射通量Rn

凈輻射通量為地表吸收的太陽短波輻射和地面有效輻射之和，EYMARD等評估了由衛(wèi)星資料估算輻射通量的方法［13］，其誤差約為25～50 W/m2，導致誤差的主要原因是不確定1 km尺度的地表反照率和比輻射率。城市區(qū)域內實測的輻射通量和反照率資料非常稀少，而衛(wèi)星技術能有效地提供城市區(qū)域輻射通量和反照率。利用全波段的地表反照率和地表溫度及植被指數資料計算得出凈輻射通量，方法如下:

Rn=K↓－K↑+L↓－L↑

=(1－α)τKsun↓－σε0Ts4+L↓

式中:α——地表反照率;τ——大氣透射率;Ksun——大氣上界處的太陽輻射，W/m2;σ——Stefan-Boltzmann 常數;ε0——地表比輻射率;Ts——地表溫度，K;L↓——向下的長波輻射通量，W/m2;參照文獻［14］計算L↓，地表比輻射率ε0則根據Van的經驗公式計算得出［15］。

2.2 土壤熱通量G

土壤熱通量是用于加熱或冷卻地表下層土壤的熱量，主要由土壤熱傳導率和土壤溫度梯度決定，它并不能由遙感測得。對于單源模式，土壤熱通量G可認為與到達地表凈輻射通量成正比，其比例系數Γ可由遙感獲得的植被特征估算得出［14］:

式中:Ts，α，Rn意義同上;NDVI——標準化差值植被指數。

3 結果分析

圖4為上海地區(qū)2007年7月28日凈輻射通量分布圖。計算結果表明，中心城區(qū)的平均凈輻射通量小于綠地，差值約為綠地的4%。這種差異更多地是由于地表反照率差異造成的，由于綠地植被較多，植被吸收太陽輻射用于光合作用，因此綠地的地表反照率低于中心城區(qū)的地表反照率。由圖可見中心城區(qū)的凈輻射通量與其周邊地區(qū)之間差異較小，根據WHITE和OKE等的實測資料結果也得出相似的結論［15－17］。

圖4 2007年7月28日上海地區(qū)凈輻射通量(W/m2)

計算結果反映了地表熱通量空間分布，由于中心城區(qū)的地表溫度高，向深層土壤傳送的熱量最多，因此中心城區(qū)的土壤熱通量最大，約高于綠地的10%，而果園、林地與綠地的土壤熱通量比較接近(圖5)。地表吸收的太陽凈輻射通量除了向下傳遞給深層土壤外，還以感熱和潛熱形式影響大氣，對感熱通量和潛熱通量的計算將在另文詳細闡述。

圖5 2007年7月28日上海地區(qū)土壤熱通量(W/m2)

4 結論

文中以TM衛(wèi)星遙感資料，結合AWS實測氣象資料，估算地表通量，從能量分配機制上分析了城市熱島效應產生的原因，同時研究表明了植被和水體具有減緩城市熱島效應。分析結果得到如下結論:

(1)交通道路地表溫度最高，其次為居民區(qū)和商業(yè)區(qū)，水面最低;

(2)中心城區(qū)地表的平均凈輻射通量與周邊相差不大，比綠地約低4%;

(3)中心城區(qū)的土壤熱通量比周邊要高，比綠地高約10%。

用TM衛(wèi)星遙感資料估算城市區(qū)域地表熱通量分布具有方便快捷的優(yōu)勢，故具有廣闊的應用前景，但是由于遙感技術易受到大氣狀況的影響，地表通量計算仍帶有經驗性，特別是因缺乏同期近地層通量觀測資料，沒有對地表通量結果進行驗證，這將有待于今后進一步研究工作的開展。

［1］CARLSON T N，GILLIES R.A method to make use of thermal infrared temperature and NDVI measurements to infer surface soil water content and fractional vegetation cover［J］.Remote Sensing Reviews，1994(9):161- 173.

［2］KATO S，YAMAGUCHI Y.Analysis of urban heat-island effect using ASTER and ETM+ data:Separation of anthropogenic heat discharge and natural heat radiation from sensible heat flux［J］.Remote Sensing of Environment，2005，99:44- 54.

［3］田武文，黃祖英，胡春娟.西安市氣候變暖與城市熱島效應問題研究［J］.應用氣象學報，2006，17(4):438- 443.

［4］謝苗苗，周偉，王仰麟，等.城市土地利用的熱環(huán)境效應研究——以寧波城區(qū)為例［J］.北京大學學報:自然科學版，2008，44(5):815- 821.

［5］BASTIAANSSEN W G M，MENENTI M，FEDDES R A，et al.A remote sensing surface energy balance algorithm for land(SEBAL)part 1.formulation［J］.Journal of Hydrology，1998(212- 213):198- 212.

［6］KUSTAS W P，HUMES K S，NORMAN J M.Single and dualsource modeling of surface energy fluxes with radiometric surface temperature［J］.Journal of Applied Meteorology，1996(35):110- 121.

［7］NORMAN J M，KUSTAS W P，PRUEGER J H，et al.Surface flux estimation using radiometric temperature:A dual-temperature-difference method to minimize measurements errors［J］.Water Resources Research，2000(36):2263- 2274.

［8］SU Z.The Surface Energy Balance System(SEBS)for estimation of turbulent heat fluxes［J］.Hydrology and Earth System Sciences，2002(6):85- 99.

［9］馬耀明，戴有學，馬偉強，等.干旱半干旱區(qū)非均勻地表區(qū)域能量通量的衛(wèi)星遙感參數化［J］.高原氣象，2004，23(2):139- 146.

［10］BASTIAANSSEN W G M.SEBAL-based sensible and latent heat fluxes in the irrigated Gediz Basin，Turkey［J］.Journal of Hydrology，2000(229):87- 100.

［11］WANG G L，JIANG W M.Analysis of urban/rural albedo from MODIS over East China［J］，Proceedings of SPIE，2006(6199):123- 134.

［12］王桂玲，蔣維楣，魏鳴.城市熱島效應的衛(wèi)星遙感分析［J］.南京氣象學院學報，2007，30(3):201- 207.

［13］EYMARD L，TACONET O.The methods for inferring surface fluxes from satellite data and their use for atmosphere model validation［J］.International Journal of Remote Sensing，1995(16):1907- 1930.

［14］PRATA A J.A new long － wave formula for estimating downward clear sky radiation at the surface［J］.Quarterly Journal of the Royal Meteorological Society，1996(122):1127- 1151.

［15］BOEGH E，SOEGAARD H，THOMSEN A.Evaluating evapotranspiration rates and surface conditions using Landsat TM to estimate atmospheric resistance and surface［J］.Remote Sensing of Environment，2002(79):329- 343.

［16］WHITE J M，EATON F D，AUER A H.The net radiation budget of the St.Louis Metropolitan area［J］.Journal of Applied Meteorology，1978，17:593- 599.

［17］OKE T R.The energetic basis of the urban heat island［J］.Quarterly Journal Royal Meteorological Society，1982，108(455):1- 24.