劉壽東 王璐瑤 胡凝 王詠薇 王亮 趙翔

摘要 利用2018年4月—2019年4月南京盤城大孔徑閃爍儀（Large Aperture Scintillometer，LAS）觀測數(shù)據(jù)，分析了城鎮(zhèn)感熱通量的時空變化特征及影響因素。結(jié)果表明：1）南京城鎮(zhèn)感熱通量呈單峰型日變化特征，白天明顯大于夜間，且白天晴天明顯大于陰天，夜間晴天略小于陰天，晴、陰天小時感熱通量年平均分別在2.25～200.53 W·m-2、13.10～132.52 W·m-2波動。2）城鎮(zhèn)感熱通量夏季明顯大于冬季，8月晝、夜分別為112.19、23.54 W·m-2，2月晝、夜分別為35.57、11.57 W·m-2。3）晴天白天條件下，不同風向（通量貢獻源區(qū)）城鎮(zhèn)感熱通量存在顯著差異，即隨著不透水層占比的增加，凈輻射分配到感熱通量的比例明顯提高，當占比大于60％時提高趨勢不明顯。4）以莫寧-奧布霍夫長度判斷大氣穩(wěn)定度為標準，C2n法在計算感熱通量的5種大氣穩(wěn)定度判斷方法中的誤判率較低且數(shù)據(jù)源于LAS，是比較適宜城鎮(zhèn)夜間大氣穩(wěn)定度的判斷方法。5）在影響城鎮(zhèn)感熱通量的地表參數(shù)中，有效高度變化的影響最大，風速變化的影響較大特別在秋冬季節(jié)更為明顯，波文比變化對城鎮(zhèn)感熱通量的影響較小，溫度、地表粗糙度和零平面位移變化的影響可忽略不計。

關(guān)鍵詞城鎮(zhèn);LAS;感熱通量;時空特征;環(huán)境因子;地表參數(shù)

地氣交換不僅影響區(qū)域熱量收支，還為不同尺度的天氣系統(tǒng)提供物質(zhì)、動力和熱力來源（馬紅云等，2018;張超等，2018;彭舒齡等，2019;蘇愛芳等，2019;孫永等，2019）。感熱通量是地氣交換的重要組成部分，準確測量地表感熱通量對于了解大氣運動以及局地氣候特征的形成機制十分重要（Weber and Kordowski，2010）。

遙感是獲取區(qū)域尺度感熱通量的主要手段（Li et al.，2009），但是由于受到遙感影像時間分辨率較低的制約，難以獲得高頻連續(xù)的感熱通量（Marx et al.，2008;劉紹民等，2010），且遙感手段需要借助地面觀測來優(yōu)化模型參數(shù)、驗證估算結(jié)果（楊凡，2012）。渦度相關(guān)儀（Eddy Covariance，EC）和大孔徑閃爍儀（Large aperture scintillometer，LAS）是兩種常用的感熱通量地面觀測儀器。隨著城市化的發(fā)展，城鎮(zhèn)下墊面異質(zhì)性越來越高，無法滿足EC方法相對均一下墊面的基本假設條件（Baldocchi et al.，2001;Massman and Lee，2002）。與EC相比，LAS可直接觀測5 km以內(nèi)的區(qū)域感熱通量，非常適用于非均勻地表感熱通量的長期觀測（Baldocchi et al.，2001;Hoedjes et al.，2002;Meijninger et al.，2002;盧俐等，2005，2009;宮麗娟等，2009;白潔等，2010a，2010b;劉紹民等，2010;楊凡，2012;Zieliński et al.，2014，2018;Ward et al.，2014;Zhang and Zhang，2015;孫根厚等，2016;Crawford et al.，2017;Ward，2017）。

LAS是基于莫寧-奧布霍夫相似性理論，結(jié)合氣象因素和地表參數(shù)，將空氣折射率結(jié)構(gòu)參數(shù)（C2n）轉(zhuǎn)換為感熱通量，因此計算過程中各參數(shù)的取值或參數(shù)化方案的選擇是感熱通量計算結(jié)果不確定性的重要來源（徐安倫等，2017）。已有的研究結(jié)果表明，有效高度和風速對LAS感熱通量觀測結(jié)果的影響最大，當有效高度±5 m、風速±1 m·s-1時，可以分別使感熱通量變化±6.3％～6.5％、±3.5％（Crawford et al.，2017）;不同的穩(wěn)定性函數(shù)造成的感熱通量誤差可達16％（Andreas，1988;Bruin et al.，1995）。但是，氣候條件不同，各種參數(shù)對于感熱通量的影響也有所區(qū)別。以波文比（β）為例，在地表干燥即β>0.6時，水汽的影響可以忽略;但在濕潤地表下即β較小時，水汽對溫度結(jié)構(gòu)參數(shù)C2T的準確性有很大影響，從而影響感熱通量（Moene，2003;盧俐等，2009）。

迄今為止，LAS在農(nóng)田、草地、森林、沙漠、黑河流域、海河流域等地區(qū)感熱通量的觀測研究較多（Kleissl et al.，2009;張勁松等，2010;朱治林等，2010;Liu et al.，2011，2013，2018;孫根厚等，2016;徐安倫等，2017），而在城鎮(zhèn)下墊面感熱通量的觀測研究尚少，主要集中在歐洲和亞洲的溫帶地區(qū)且觀測時間較短（Lagouarde et al.，2006;Wood et al.，2013;Ward et al.，2014;Zieliński et al.，2014，2018;Jacobs et al.，2015;Lee et al.，2015;Zhang and Zhang，2015;Crawford et al.，2017;Ward，2017）。因此，本項研究通過亞熱帶季風氣候區(qū)城鎮(zhèn)LAS的長期觀測，分析城鎮(zhèn)下墊面感熱通量時空變化特征及影響因子，為明晰城鎮(zhèn)地氣交換機制提供依據(jù)，為城市感熱通量的觀測研究提供可行方法和有效數(shù)據(jù)。

1 資料和方法

1.1 觀測實驗方法

LAS、EC的觀測于2018年4月—2019年4月在江蘇省南京市江北新區(qū)盤城鎮(zhèn)南京信息工程大學進行（118°42′ E，32°12′ N）。南京市屬于北亞熱帶濕潤氣候，四季分明，雨水充沛。觀測實驗區(qū)土地利用類型見圖1?？梢娬麄€試驗區(qū)以建筑物、街道等不透水層為主占比約42.85％，植被占比約35.82％，具有較明顯的城鎮(zhèn)下墊面特征。

LAS（荷蘭Kipp ＆ Zonen公司生產(chǎn)的MKII型）架設在校園內(nèi)建筑物樓頂，其中發(fā)射端（T）和接收端（R）分別距地面36 m和26 m，直線距離為1.1 km，采樣頻率為1 Hz;EC架設在LAS接收端以北臨近高塔上，安裝高度為60 m，采樣頻率為10 Hz。此外，高塔上還架設了三層溫度、風速和風向傳感器，安裝高度分別為10 m、30 m和60 m。

1.2 觀測數(shù)據(jù)處理及計算方法

LAS發(fā)射端在發(fā)出一定波長和直徑的光束后，由于大氣溫度、濕度和氣壓波動引起大氣折射系數(shù)波動，光束強度會發(fā)生改變，由此得到空氣折射率結(jié)構(gòu)參數(shù)（C2n）。大孔徑閃爍儀光束位于近紅外波段，對溫度變化更加敏感，因此C2n與溫度結(jié)構(gòu)參數(shù)C2T有關(guān)（Thiermann and Grassl，1992;Crawford et al.，2017），可以用來表征大氣中湍流運動特征。

C2T=C2n T2-0.78×10-6P21+0.03β。? （1）

其中：T是空氣溫度（單位：K），采用架設在EC高塔上30 m高度處的溫度值;P是氣壓（單位：Pa），來源于渦度相關(guān)系統(tǒng);β是波文比。因LAS發(fā)射波長為880 nm的電磁波，主要受氣溫的影響，水汽影響可忽略，這里β取觀測數(shù)據(jù)的平均值0.68。

根據(jù)MOST相似理論，由式（2）結(jié)合莫寧-奧布霍夫長度L=u2* TgkT*通過迭代計算可得到溫度尺度（T*）（Bruin et al.，1995）。將T*代入后，HLAS=-ρcpu*T*即可以得到感熱通量HLAS，最后采用足跡模型來確定源區(qū)范圍（Meijninger et al.，2002）。

T*=C2T （Zeff-d）2／3fTZeff-dL1／2。? （2）

其中：u*是摩擦速度（單位：m·s-1）（Zieliński et al.，2014）;Zeff是有效高度（單位：m），利用Evation中的數(shù)字地形模塊，測量光程路徑上的建筑物高度，結(jié)合空間權(quán)重函數(shù)，求得Zeff為16.70 m;d是零平面位移，利用風廓線計算得出d為4.56 m;ξ=Zeff-dL是層結(jié)參數(shù)，fT是穩(wěn)定度函數(shù)，當大氣不穩(wěn)定時，fT（ξ）=4.9（1-6.1ξ）-2／3;當大氣穩(wěn)定時，fT（ξ）=4.9（1+2.2ξ）2／3（Andreas，1988;Bruin et al.，1995）。本研究利用10 m和60 m高度上的溫度梯度判斷大氣穩(wěn)定與否。當T／z<0時大氣為不穩(wěn)定狀態(tài)，T／z>0則為穩(wěn)定狀態(tài)。

LAS數(shù)據(jù)采用Evation軟件處理，包括降水數(shù)據(jù)以及C2n超過飽和上限（3.88×10-13）及信號強度過?。?.02×10-17）數(shù)據(jù)的剔除。EC數(shù)據(jù)采用Eddy pro軟件處理，包括野點值去除、延遲時間的校正、坐標旋轉(zhuǎn)和WPL校正等。通過數(shù)據(jù)質(zhì)量控制后用于城鎮(zhèn)感熱通量分析的為2018年4月2日至2019年4月30日的感熱通量值;2018年4月19日至6月13日由于儀器維修數(shù)據(jù)缺失。

2 結(jié)果與分析

2.1 城鎮(zhèn)感熱通量時間變化特征

觀測實驗期間南京盤城鎮(zhèn)全年典型晴天和陰天條件下感熱通量HLAS和凈輻射的平均日變化情況見圖2。從圖中可以看出，城鎮(zhèn)感熱通量與凈輻射變化趨勢一致，白天明顯大于夜間，呈現(xiàn)出明顯的單峰型日變化特征。典型晴天條件下，白天日出后，隨著凈輻射的不斷增強，HLAS也在快速增大，在中午凈輻射達到峰值495.40 W·m-2后約1 h出現(xiàn)HLAS最大值200.53 W·m-2，感熱通量出現(xiàn)滯后現(xiàn)象主要是因為建筑群熱容量大，熱量儲存多，而當不透水層占比超過50％時，隨著不透水層比例的增加，感熱通量的滯后現(xiàn)象越明顯（Grimmond and Oke，1995;Grimmond et al.，2004;Ward et al.，2014;Crawford et al.，2017）;此后隨著凈輻射減弱HLAS快速減小。晴天夜間，由于建筑群、水泥路面熱容量更大，使得城鎮(zhèn)地區(qū)存儲了較多熱量，造成城鎮(zhèn)夜間大氣層結(jié)不穩(wěn)定，導致HLAS夜間為正但整體上小于11.00 W·m-2，最小僅2.25 W·m-2。典型陰天條件下，凈輻射和HLAS日變化特征與典型晴天基本一致，主要區(qū)別在于HLAS增大與減小的速率小于典型晴天，白天凈輻射、HLAS的峰值均明顯低于晴天分別為345.16 W·m-2、132.52 W·m-2。而在陰天夜間，由于云層存在保溫作用造成地氣間有明顯溫差，凈輻射和HLAS雖整體上小于20 W·m-2但略大于晴天，HLAS最小為13.10 W·m-2，這與在低丘紅壤區(qū)的觀測結(jié)果基本一致（李陽等，2017）。此外，由于晴天條件下地表吸收更多的熱量，同時城鎮(zhèn)下墊面不透水層占比大、蒸發(fā)較弱，導致晴天感熱通量與凈輻射的比例HLAS／Rn也略低于陰天。

城鎮(zhèn)感熱通量存在明顯的季節(jié)性變化（圖3），這主要是由不同季節(jié)的凈輻射差異造成的，即感熱通量隨著凈輻射增強而增大。圖3a顯示，觀測期間城鎮(zhèn)春、夏、秋季白天的HLAS明顯高于冬季，8月白天月平均感熱通量最高為112.19 W·m-2，2月最低僅為35.57 W·m-2。但不同季節(jié)凈輻射中用于感熱支出的比例HLAS／Rn差異較大，觀測期間各月HLAS／Rn在39％～64％，表現(xiàn)為冬季高而夏季低。主要是因為冬季氣溫低、蒸散較弱，有更多的凈輻射分配給感熱支出，且城鎮(zhèn)居民取暖產(chǎn)生大量人為熱，使得感熱通量占凈輻射的比例更高;夏季則正好相反，氣溫高，降水多，蒸散強烈，凈輻射分配給感熱支出的比例低，7月白天的HLAS／Rn僅為39％;也正是這種能量收支和分配的差異造成了不同下墊面和不同季節(jié)的感熱通量有所不同（Ward et al.，2014;Crawford et al.，2017;Ward，2017;Zieliński et al.，2018）。

夜間，各月HLAS的平均值均為正值，同樣表現(xiàn)為夏季大而冬季小。7、8月夜間氣溫較高，地氣之間湍流交換作用非常強烈，月平均HLAS最高分別為25.29 W·m-2、23.54 W·m-2。冬季，由于城鎮(zhèn)存在大量人為熱，導致HLAS夜間為正但整體較小，2月夜間平均HLAS降至最小為11.57 W·m-2。

2.2 城鎮(zhèn)感熱通量的空間代表性

因為通量貢獻區(qū)分布在盛行風向一側(cè)，所以在不均勻下墊面上HLAS會隨著風向和大氣穩(wěn)定度的不同而發(fā)生變化。為了分析LAS觀測結(jié)果的空間代表性，選取向下短波輻射在200～400 W·m-2、降雨12 h之后的HLAS數(shù)據(jù)進行統(tǒng)計，以確保各個風向上HLAS的差異均來自干燥的地表、相近的太陽輻射（Kotthaus and Grimmond，2014;Crawford et al.，2017），統(tǒng)計結(jié)果見圖4。從圖中可以看出，HLAS大多在50～100 W·m-2，大于150 W·m-2的HLAS多出現(xiàn)在風向為ENE、E、ESE、SE時，而小于50 W·m-2的HLAS多分布在SSE、SSW、SW、WSW，這主要與城鎮(zhèn)的下墊面類型有關(guān)（王亮等，2020）。

為了進一步分析下墊面對城鎮(zhèn)感熱通量的影響，繪制不同風向下通量源區(qū)內(nèi)的不透水層（建筑物和街道）占比、植被占比情況（圖5）。

從圖5中可以看出，實驗區(qū)的北部、東部（0°～100°、240°～360°）以建筑物和街道為主，表面由水泥等不透水材料組成，不透水層占比最高平均達到70％，密集建筑使得HLAS較高、變化幅度大，常出現(xiàn)大于250 W·m-2的HLAS 值。而南部、西部（160°～230°）的植被占比較高，不透水層占比低平均僅為47％，HLAS也較低。主要是因為不透水層不僅自身的保水性差，且阻礙大氣與土壤水分的直接交換，在降水時水分因無法下滲而很快流走，沒有降水時土壤水分難以通過植物蒸騰和土壤蒸發(fā)穿過不透水層進入大氣而消耗更多能量。這一觀測結(jié)果與在倫敦（Crawford et al.，2017）、羅茲（Zieliński et al.，2018）等城區(qū)的觀測結(jié)果一致。由此可見，晴天白天條件下，不同風向（通量貢獻源區(qū)）的城鎮(zhèn)感熱通量存在顯著的差異，較大的感熱通量值多分布在建筑物和人為活動的高密度風向區(qū)，而較小的感熱通量值多分布在植被占比更大的風向區(qū)。王亮等（2020）指出，在城鎮(zhèn)地區(qū)不透水層占比對于能量分配有很大影響。圖6顯示，隨著不透水層占比增加，凈輻射中通過感熱交換方式散失能量的比例HLAS／Rn逐漸提高;而當不透水層占比超過60％時，HLAS／Rn提高趨勢不明顯，主要原因是不透水層占比超過60％時通量貢獻源區(qū)以不透水層為主，感熱通量中的城市存儲熱比例增加，凈輻射中的感熱支出部分減少（Lin et al.，2015）。

2.3 大氣穩(wěn)定度判斷方法的影響

大氣穩(wěn)定度的判斷結(jié)果對感熱通量觀測結(jié)果的準確性至關(guān)重要。為了比較不同大氣穩(wěn)定度判斷方法的影響，引入另外4種大氣穩(wěn)定度判斷方法（表1）進行計算分析。

采用5種大氣穩(wěn)定度判斷方法計算的HLAS日變化見圖7。由圖可見，四季5種方法計算的HLAS日變化趨勢一致，且白天HLAS數(shù)值接近，各種方法之間幾乎不存在差異。各種方法計算的HLAS日變化之間的差異主要有兩個，一是大氣穩(wěn)定度的轉(zhuǎn)換時間，這與在北京密云（盧俐等，2009）和比利時德河地區(qū)（Samain et al.，2012）的研究結(jié)果類似;二是夜間5種方法計算的HLAS數(shù)值差異明顯。從大氣穩(wěn)定度轉(zhuǎn)換時間看，M1方法判斷的轉(zhuǎn)換時間早晨最晚、傍晚最早，因此M1方法判斷的大氣不穩(wěn)定層結(jié)持續(xù)時間是最短的;而M5方法判斷夜間大氣始終為不穩(wěn)定層結(jié)，故M5方法判斷的大氣不穩(wěn)定層結(jié)持續(xù)時間是最長的，比M1方法平均長14.5 h。而M2、M3和M4方法的判斷結(jié)果較為接近，大氣穩(wěn)定度的轉(zhuǎn)換時間相差在0.5～2.5 h。

為了明確不同大氣穩(wěn)定度判斷方法的適用性，將5種判斷方法計算的HLAS與HEC在小時尺度上進行對比。圖8顯示，白天5種判斷方法計算的HLAS結(jié)果與HEC的一致性均很好。

5種大氣穩(wěn)定度判斷方法計算的HLAS與HEC之間差異主要在夜間（圖8）。M5方法是常用的大氣穩(wěn)定度判斷方法，但是用M5方法計算的觀測期間夜間HLAS值有約43％為HLAS>0而HEC<0，這可能是因為夜間有大量城市存儲熱存在，在這種情況下用M5方法判斷城鎮(zhèn)夜間大氣多為不穩(wěn)定層結(jié);若以莫寧-奧布霍夫長度L判斷的大氣穩(wěn)定度作為標準，M5方法夜間大氣穩(wěn)定度的誤判率高達60％，可見M5方法并不適用于城鎮(zhèn)夜間大氣穩(wěn)定度的判斷。而M2方法不采用觀測數(shù)據(jù)來判斷大氣穩(wěn)定度，與M1方法的判斷結(jié)果相似，夜間大氣層結(jié)均判斷為穩(wěn)定，但是城鎮(zhèn)白天存儲了較多的熱量，夜間也可能出現(xiàn)大氣層結(jié)不穩(wěn)定的現(xiàn)象，大氣穩(wěn)定度的誤判率達到34％，因此M1和M2方法也不適用于城鎮(zhèn)地區(qū)。用M4方法計算的夜間HLAS值有約28％為HLAS<0而HEC>0，主要是因為夜間湍流較弱，采用C2n最小值法判斷大氣穩(wěn)定度可能會存在誤判，夜間大氣穩(wěn)定度的誤判率為35％。而用M3方法計算的夜間HLAS與HEC的一致性較好，其夜間大氣穩(wěn)定度的誤判率最低為21％。綜合來看，雖然M4方法夜間大氣穩(wěn)定度誤判率高于M3方法，但明顯低于M5方法，且與其他方法相比C2n不需要額外的觀測數(shù)據(jù)就可以通過LAS觀測即可得到，因此在計算HLAS時采用M4方法判斷夜間大氣穩(wěn)定度是比較適宜的（Samain et al.，2012）。但從長遠考慮，選擇適用的大氣穩(wěn)定度判斷方法、顯著降低夜間誤判率尚需進一步探索。

2.4 環(huán)境因子和地表參數(shù)的影響

u*和T*是決定感熱通量的兩個主要因素，影響它們的環(huán)境因子有T、P、u和β，地表參數(shù)有d、z0（地表粗糙度）和Zeff。選取大氣不穩(wěn)定條件下的12 457個樣本，參考已有的研究，以±0.3、±1 ℃、±0.1 m·s-1、±0.1 m、±0.5 m、±1.5 m作為β、T、u、z0、d、Zeff的改變量（盧俐等，2009;Crawford et al.，2017），用來分析各參數(shù)變化對感熱通量的敏感性，結(jié)果如圖9所示。

從圖9中可以看出，在環(huán)境因子中，氣溫變化對HLAS的影響最小。當氣溫改變±1 ℃時，HLAS的變化小于（0.57±0.06）％，可見氣溫變化的影響可忽略不計。比較而言，風速和波文比的變化對HLAS的影響更為顯著。當u變化±0.1 m·s-1時，HLAS變化±（2.19±1.88）％;但在不同月份u的變化對HLAS的影響是有區(qū)別的，當年9月至次年1月，u改變±0.1 m·s-1時，HLAS平均變化±2.31％，明顯高于其他月份，主要是因為，秋、冬季的熱力作用較春、夏季弱，風速切變引起的機械湍流貢獻更大（張強等，2011）。另外，波文比β變化對HLAS的影響并不對稱，β取值越小其變化對感熱通量的影響越明顯，但隨著β的增大其變化對HLAS的影響逐漸減弱。已有的研究結(jié)果表明，在β>1時，其變化對HLAS的影響不大（Moene，2003;朱治林等，2010;Ward et al.，2014）。本次試驗結(jié)果顯示，實測的β有48％的數(shù)據(jù)大于0.68，用固定值0.68代替實測的β對HLAS的影響小于1.62％。因此，在城鎮(zhèn)或者其他干旱地表，β變化對HLAS的影響可以忽略不計，β采用固定值是適宜的。

從圖中還可以看出，不同月份地表參數(shù)的變化對HLAS的影響基本類似。其中，當Zeff改變±1.5 m時，HLAS變化為±（7.50±2.30）％，這一結(jié)果與在北京昌平小湯山（盧俐等，2009）和英國倫敦（Crawford et al.，2017）的敏感性試驗結(jié)果一致，即Zeff變化對HLAS的影響最大。在英國斯溫頓郊區(qū)，Zeff改變±2.25 m時HLAS變化±3.0％（Ward et al.，2014）;而在倫敦密集的城市中心，Zeff改變±5.0 m時，HLAS變化±6.4％（Crawford et al.，2017），可見在本項觀測實驗條件下有效高度Zeff變化對HLAS的影響更大。另外，z0改變±0.1 m和d改變±0.5 m對HLAS的影響均小于±2.80％，因此z0和d變化對HLAS的影響最小。

3 結(jié)論與討論

1）觀測實驗期間南京城鎮(zhèn)感熱通量與凈輻射變化趨勢一致，白天明顯大于夜間，呈現(xiàn)明顯的單峰型日變化特征;城鎮(zhèn)感熱通量白天晴天明顯大于陰天，夜間晴天略小于陰天，晴、陰天小時感熱通量年平均分別在2.25～200.53 W·m-2、13.10～132.52 W·m-2波動。城鎮(zhèn)感熱通量夏季明顯大于冬季，8月晝、夜平均感熱通量分別為112.19、23.54 W·m-2，2月晝、夜平均感熱通量分別為35.57、11.57 W·m-2。

2）晴天白天條件下，不同風向（通量貢獻源區(qū)）城鎮(zhèn)感熱通量存在顯著差異，即隨通量貢獻源區(qū)下墊面不透水層占比的增加凈輻射分配到感熱通量的比例明顯提高，當占比大于60％時提高趨勢不明顯。

3）采用5種大氣穩(wěn)定度判斷方法計算的四季城鎮(zhèn)感熱通量的日變化趨勢一致，且白天數(shù)值接近，各種方法之間的差異主要為大氣穩(wěn)定度的轉(zhuǎn)換時間不一致和夜間城鎮(zhèn)感熱通量不同。以莫寧-奧布霍夫長度L判斷大氣穩(wěn)定度為標準，溫度梯度法夜間大氣穩(wěn)定度誤判率最高達60％，Rn法、日出日落時間法和C2n法誤判率相近約35％，EC法的誤判率最低僅21％。而C2n法誤判率較低且數(shù)據(jù)源于LAS，是當前比較適宜判斷夜間大氣穩(wěn)定度的方法。但是從長遠考慮，選擇適用的大氣穩(wěn)定度判斷方法、顯著降低夜間誤判率尚需進一步探索。

4）影響城鎮(zhèn)感熱通量的地表參數(shù)中，有效高度變化的影響最大，風速變化的影響較大尤其在秋冬季節(jié)更為明顯，溫度、地表粗糙度和零平面位移變化的影響可忽略不計。波文比取值越小其變化對感熱通量的影響越明顯，但是隨著波文比的增大其變化對感熱通量的影響逐漸減弱，在城鎮(zhèn)等干旱地表波文比的變化對感熱通量的影響可以忽略不計，波文比采用固定值在城鎮(zhèn)是適宜的。

本項研究利用一年時間的LAS觀測資料進行了城鎮(zhèn)感熱通量時空變化及影響因子的分析研究，但國內(nèi)仍缺乏不同氣候區(qū)域、不同城市（城鎮(zhèn)）下墊面類型LAS感熱通量的長期觀測實驗，需要進一步開展同類觀測實驗，以得出更有說服力的結(jié)論。

參考文獻（References）

Andreas E L，1988.Estimating C2n over snow and sea ice from meteorological data［J］.J Opt Soc Am A，5（4）：481-495.doi：10.1364/JOSAA.5.000481.

白潔，劉紹民，丁曉萍，等，2010a.大孔徑閃爍儀觀測數(shù)據(jù)的處理方法研究［J］.地球科學進展，25（11）：1148-1165. Bai J，Liu S M，Ding X P，et al.，2010a.A study of the processing method of large aperture scintillometer observation data［J］.Adv Earth Sci，25（11）：1148-1165.doi：10.3969／j.issn.1001-7313.2009.02.006.（in Chinese）.

白潔，劉紹民，丁曉萍，2010b.海河流域不同下墊面上大孔徑閃爍儀觀測顯熱通量的時空特征分析［J］.地球科學進展，25（11）：1187-1198. Bai J，Liu S M，Ding X P，2010b.Temporal-spatial characteristics analysis of area-averaged sensible heat fluxes by large aperture scintillometer over Hai River Basin［J］.Adv Earth Sci，25（11）：1187-1198.doi：10.3969／j.issn.1001-7313.2009.02.006.（in Chinese）.

Baldocchi D，F(xiàn)alge E，Gu L H，et al.，2001.FLUXNET：a new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide，water vapor，and energy flux densities［J］.Bull Amer Meteor Soc，82（11）：2415-2434.doi：10.1175／1520-0477（2001）082<2415：fantts>2.3.co;2.

Bruin H A R，Hurk B J J M，Kohsiek W，1995.The scintillation method tested over a dry vineyard area［J］.Bound Layer Meteor，76（1／2）：25-40.doi：10.1007／BF00710889.

Crawford B，Grimmond C S B，Ward H C，et al.，2017.Spatial and temporal patterns of surface-atmosphere energy exchange in a dense urban environment using scintillometry［J］.Qurat J Roy Meteor Soc，143（703）：817-833.doi：10.1002／qj.2967.

宮麗娟，劉紹民，雙喜，等，2009.渦動相關(guān)儀和大孔徑閃爍儀觀測通量的空間代表性［J］.高原氣象，28（2）：246-257. Gong L J，Liu S M，Shuang X，et al.，2009.Investigation of spatial representativeness for surface flux measurements with eddy covariance system and large aperture scintillometer［J］.Plateau Meteor，28（2）：246-257.（in Chinese）.doi：CNKI：SUN：GYQX.0.2009-02-002.

Grimmond C S B，Oke T R，1995.Comparison of heat fluxes from summertime observations in the suburbs of four North American cities［J］.J Appl Meteor，34（4）：873-889.doi：10.1175／1520-0450（1995）034<0873：cohffs>2.0.co;2.

Grimmond C S B，Salmond J A，Oke T R，et al.，2004.Flux and turbulence measurements at a densely built-up site in Marseille：heat，mass （water and carbon dioxide），and momentum［J］.J Geophys Res Atmos，109（D24）：D24101.doi：10.1029／2004JD004936.

Hoedjes J B，Zuurbier R M，Watts C J，2002.Large aperture scintillometer used over a homogeneous irrigated area，partly affected by regional advection［J］.Bound Layer Meteor，105（1）：99-117.doi：10.1023／A：1019644420081.

Jacobs C，Elbers J，Brolsma R，et al.，2015.Assessment of evaporative water loss from Dutch Cities［J］.Build Environ，83：27-38.doi：10.1016／j.buildenv.2014.07.005.

Kleissl J，Gomez J，Hong S H，et al.，2008.Large aperture scintillometer intercomparison study［J］.Bound Layer Meteor，128（1）：133-150.doi：10.1007／s10546-008-9274-1.

Kleissl J，Hong S H，Hendrickx J M H，2009.New Mexico scintillometer network：supporting remote sensing and hydrologic and meteorological models［J］.Bull Amer Meteor Soc，90（2）：207-218.doi：10.1175／2008bams2480.1.

Kohsiek W，Meijninger W L，Debruin H R，et al.，2006.Saturation of the large aperture scintillometer［J］.Bound Layer Meteor，121（1）：111-126.doi：10.1007／s10546-005-9031-7.

Kotthaus S，Grimmond C S B，2014.Energy exchange in a dense urban environment-Part I：temporal variability of long-term observations in central London［J］.Urban Climate，10：261-280.doi：10.1016／j.uclim.2013.10.002.

Lagouarde J P，Irvine M，Bonnefond J M，et al.，2006.Monitoring the sensible heat flux over urban areas using large aperture scintillometry：case study of Marseille city during the escompte experiment［J］.Bound Layer Meteor，118（3）：449-476.doi：10.1007／s10546-005-9001-0.

Lee S H，Lee J H，Kim B Y，2015.Estimation of turbulent sensible heat and momentum fluxes over a heterogeneous urban area using a large aperture scintillometer［J］.Adv Atmos Sci，32（8）：1092-1105.doi：10.1007／s00376-015-4236-2.

李陽，景元書，秦奔奔，2017.低丘紅壤區(qū)農(nóng)田水熱通量變化特征及氣候?qū)W足跡［J］.應用生態(tài)學報，28（1）：180-190. Li Y，Jing Y S，Qin B B，2017.Characteristics of water and heat fluxes and its footprint climatology on farmland in low hilly region of red soil［J］.Chin J Appl Ecol，28（1）：180-190.doi：10.13287／j.1001-9332.201701.032.（in Chinese）.

Li Z L，Tang R L，Wan Z M，et al.，2009.A review of current methodologies for regional evapotranspiration estimation from remotely sensed data［J］.Sensors，9（5）：3801-3853.doi：10.3390／s90503801.

Lin C Y，Su C J，Kusaka H，et al.，2015.Impact of an improved WRF-urban canopy model on diurnal air temperature simulation over northern Taiwan［J］.Atmos Chem Phys，15（20）：28483-28516.doi：10.5194／acpd-15-28483-2015.

劉紹民，李小文，施生錦，等，2010.大尺度地表水熱通量的觀測、分析與應用［J］.地球科學進展，25（11）：1113-1127. Liu S M，Li X W，Shi S J，et al.，2010.Measurement，analysis and application of surface energy and water vapor fluxes at large scale［J］.Adv Earth Sci，25（11）：1113-1127.doi：10.3969／j.issn.1001-7313.2009.02.006.（in Chinese）.

Liu S M，Xu Z W，Wang W Z，et al.，2011.A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem［J］.Hydrol Earth Syst Sci，15（4）：1291-1306.doi：10.5194／hess-15-1291-2011.

Liu S M，Xu Z W，Zhu Z L，et al.，2013.Measurements of evapotranspiration from eddy-covariance systems and large aperture scintillometers in the Hai River Basin，China［J］.J Hydrol，487：24-38.doi：10.1016／j.jhydrol.2013.02.025.

Liu S M，Li X，Xu Z W，et al.，2018.The Heihe integrated observatory network：a basin-scale land surface processes observatory in China［J］.Vadose Zone J，17（1）：1-21.doi：10.2136／vzj2018.04.0072.

盧俐，劉紹民，孫敏章，等，2005.大孔徑閃爍儀研究區(qū)域地表通量的進展［J］.地球科學進展，20（9）：932-938. Lu L，Liu S M，Sun M Z，et al.，2005.Advances in the study of areal surface fluxes with large aperture scintillometer［J］.Adv Earth Sci，20（9）：932-938.doi：10.3321／j.issn：1001-8166.2005.09.002.（in Chinese）.

盧俐，劉紹民，徐自為，等，2009.不同下墊面大孔徑閃爍儀觀測數(shù)據(jù)處理與分析［J］.應用氣象學報，20（2）：171-178. Lu L，Liu S M，Xu Z W，et al.，2009.Results from measurements of large aperture scintillometer over different surfaces［J］.J Appl Meteorol Sci，20（2）：171-178.doi：10.3969／j.issn.1001-7313.2009.02.006.（in Chinese）.

馬紅云，董軒，孫岑霄，2018.南京地區(qū)高溫天氣下城市化影響的模擬研究［J］.大氣科學學報，41（1）：67-76. Ma H Y，Dong X，Sun C X，2018.Numerical simulations of urbanization impacts under hot weather conditions in Nanjing［J］.Trans Atmos Sci，41（1）：67-76.doi：10.13878／j.cnki.dqkxxb.20170607001.（in Chinese）.

Marx A，Kunstmann H，Schüttemeyer D，et al.，2008.Uncertainty analysis for satellite derived sensible heat fluxes and scintillometer measurements over Savannah environment and comparison to mesoscale meteorological simulation results［J］.Agric For Meteorol，148（4）：656-667.doi：10.1016/j.agrformet.2007.11.009.

Massman W J，Lee X，2002.Eddy covariance flux corrections and uncertainties in long-term studies of carbon and energy exchanges［J］.Agr Forest Meteor，113（1／2／3／4）：121-144.doi：10.1016／S0168-1923（02）00105-3.

Meijninger W M L，Hartogensis O K，Kohsiek W，et al.，2002.Determination of area-averaged sensible heat fluxes with a large aperture scintillometer over a heterogeneous surface-Flevoland field experiment［J］.Bound Layer Meteor，105（1）：37-62.doi：10.1023／A：1019647732027.

Moene A F，2003.Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation［J］.Bound Layer Meteor，107（3）：635-653.doi：10.1023／A：1022807617073.

彭舒齡，周樹道，衛(wèi)克晶，等，2019.京津冀地區(qū)一次強沙塵天氣過程的成因及特征［J］.大氣科學學報，42（6）：926-935. Peng S L，Zhou S D，Wei K J，et al.，2019.Causes and characteristics of a dust weather process in the Beijing-Tianjin-Hebei region［J］.Trans Atmos Sci，42（6）：926-935.doi：10.13878／j.cnki.dqkxxb.20180428001.（in Chinese）.

Samain B，F(xiàn)erket B V A，Defloor W，et al.，2011.Estimation of catchment averaged sensible heat fluxes using a large aperture scintillometer［J］.Water Resour Res，47（5）：33-44.W05536.doi：10.1029／2009WR009032.

Samain B，Defloor W，Pauwels V R N，2012.Continuous time series of catchment-averaged sensible heat flux from a large aperture scintillometer：efficient estimation of stability conditions and importance of fluxes under stable conditions［J］.J Hydrometeorol，13（2）：423-442.doi：10.1175／jhm-d-11-030.1.

蘇愛芳，施東雷，葛旭陽，2019.下墊面對鄭州城市強降水的影響：城市化及地形影響的數(shù)值模擬研究［J］.大氣科學學報，42（3）：434-446. Su A F，Shi D L，Ge X Y，2019.Numerical simulation of the influence from urbanization and orography on a severe rainfall event in Zhengzhou［J］.Trans Atmos Sci，42（3）：434-446.doi：10.13878／j.cnki.dqkxxb.20181116003.（in Chinese）.

孫根厚，胡澤勇，王介民，等，2016.那曲地區(qū)兩種空間尺度感熱通量的對比分析［J］.高原氣象，35（2）：285-296. Sun G H，Hu Z Y，Wang J M，et al.，2016.Comparison analysis of sensible heat fluxes at two spatial scales in Naqu area［J］.Plateau Meteor，35（2）：285-296.doi：10.7522／j.issn.1000-0534.2015.00088.（in Chinese）.

孫永，王詠薇，高陽華，等，2019.復雜地形條件下城市熱島及局地環(huán)流特征的數(shù)值模擬［J］.大氣科學學報，42（2）：280-292. Sun Y，Wang Y W，Gao Y H，et al.，2019.Numerical simulation of urban heat island and local circulation characteristics under complex terrain conditions［J］.Trans Atmos Sci，42（2）：280-292.doi：10.13878／j.cnki.dqkxxb.20180204001.（in Chinese）.

Thiermann V，Grassl H，1992.The measurement of turbulent surface-layer fluxes by use of bichromatic scintillation［J］.Bound Layer Meteor，58（4）：367-389.doi：10.1007／BF00120238.

王亮，胡凝，王詠薇，等，2020.南京地區(qū)大口徑閃爍儀與渦動相關(guān)儀感熱通量觀測對比［J］.生態(tài)學雜志，39（1）：315-325. Wang L，Hu N，Wang Y W，et al.，2020.Comparative analysis of sensible heat flux observed by large aperture scintillometer and eddy-covariance system in Nanjing［J］.Chin J Ecol，39（1）：315-325.doi：10.13292／j.1000-4890.202001.005.（in Chinese）.

Ward H C，2017.Scintillometry in urban and complex environments：a review［J］.Meas Sci Technol，28：（6）：064005.doi：10.1088／1361-6501／aa5e85.

Ward H C，Evans J G，Grimmond C S B，2014.Multi-scale sensible heat fluxes in the suburban environment from large-aperture scintillometry and eddy covariance［J］.Bound Layer Meteor，152（1）：65-89.doi：10.1007／s10546-014-9916-4.

Weber S，Kordowski K，2010.Comparison of atmospheric turbulence characteristics and turbulent fluxes from two urban sites in Essen，Germany［J］.Theor Appl Climatol，102（1／2）：61-74.doi：10.1007／s00704-009-0240-8.

Wesely M L，Alcaraz E C，1973.Diurnal cycles of the refractive index structure function coefficient［J］.J Geophys Res，78（27）：6224-6232.doi：10.1029／JC078i027p06224.

Wood C R，Kouznetsov R D，Gierens R，et al.，2013.On the temperature structure parameter and sensible heat flux over Helsinki from sonic anemometry and scintillometry［J］.J Atmos Ocean Technol，30（8）：1604-1615.doi：10.1175/JTECH-D-12-00209.1.

徐安倫，李建，彭浩，等，2017.洱海湖濱農(nóng)田下墊面大口徑閃爍儀與渦動相關(guān)儀測量的湍流熱通量對比分析［J］.高原氣象，36（1）：98-106. Xu A L，Li J，Peng H，et al.，2017.Comparison of sensible and latent heat fluxes over farmland underlying surface of Erhai lakeside region measured from large aperture scintillometer and eddy covariance system［J］.Plateau Meteor，36（1）：98-106.doi：10.7522／j.issn.1000-0534.2016.00009.（in Chinese）.

楊凡，2012.農(nóng)田蒸散量的多尺度監(jiān)測與驗證研究［D］.石家莊：河北師范大學. Yang F，2012.Multi-scale research on monitoring and verification of farmland evapotranspiration［D］.Shijiazhuang：Hebei Normal University.（in Chinese）.

Zhang H，Zhang H S，2015.Comparison of turbulent sensible heat flux determined by large-aperture scintillometer and eddy covariance over urban and suburban areas［J］.Bound Layer Meteor，154（1）：119-136.doi：10.1007／s10546-014-9965-8.

張超，田榮湘，茆慧玲，等，2018.青藏高原4月感熱通量異常對長江以南夏季降水的影響［J］.大氣科學學報，41（6）：775-785. Zhang C，Tian R X，Mao H L，et al.，2018.Impact of the sensible heat flux anomaly over the Tibetan Plateau in April on summer precipitation in the south of the Yangtze River Region［J］.Trans Atmos Sci，41（6）：775-785.doi：10.13878／j.cnki.dqkxxb.20170124003.（in Chinese）.

張勁松，孟平，鄭寧，等，2010.大孔徑閃爍儀法測算低丘山地人工混交林顯熱通量的可行性分析［J］.地球科學進展，25（11）：1283-1290. Zhang J S，Meng P，Zheng N，et al.，2010.The feasibility of using LAS measurements of the sensible heat flux from a mixed plantation in the hilly zone of the North China［J］.Adv Earth Sci，25（11）：1283-1290.doi：10.3321／j.issn：0479-8023.2006.03.006.（in Chinese）.

張強，張杰，喬娟，等，2011.我國干旱區(qū)深厚大氣邊界層與陸面熱力過程的關(guān)系研究［J］.中國科學（地球科學），41（9）：1365-1374. Zhang Q，Zhang J，Qiao J，et al.，2011.Study on the relationship between deep atmospheric boundary layer and land surface thermal process in arid areas of China［J］.Sci Sin Terrae，41（9）：1365-1374.doi：10.3321／j.issn：1001-8166.2009.04.002.（in Chinese）.

朱治林，孫曉敏，賈媛媛，等，2010.基于大孔徑閃爍儀（LAS）測定農(nóng)田顯熱通量的不確定性分析［J］.地球科學進展，25（11）：1199-1207. Zhu Z L，Sun X M，Jia Y Y，et al.，2010.The uncertainty of sensible heat flux measurement over cropland using large aperture scintillometer （LAS）［J］.Adv Earth Sci，25（11）：1199-1207.doi：10.3969／j.issn.1001-7313.2009.02.006.（in Chinese）.

Zieliński M，F(xiàn)ortuniak K，Pawlak W，et al.，2014.Turbulent sensible heat flux in ódz＇，central Poland，obtained from scintillometer and eddy covariance measurements［J］.Meteor Z，22（5）：603-613.doi：10.1127／0941-2948／2013／0448.

Zieliński M，F(xiàn)ortuniak K，Pawlak W，et al.，2018.Long-term turbulent sensible-heat-flux measurements with a large-aperture scintillometer in the centre of ódz＇，central Poland［J］.Bound Layer Meteor，167（3）：469-492.doi：10.1007／s10546-017-0331-5.

·ARTICLE·

Spatial-temporal characteristics and influencing factors of sensible heat flux based on LAS in a town

LIU Shoudong1，2，WANG Luyao1，2，HU Ning1，2，WANG Yongwei1，2，WANG Liang1，2，ZHAO Xiang1，2

1Yale-NUIST Center on Atmospheric Environment，International Joint Laboratory on Climate and Environment Change （ILCEC），Nanjing University of Information Science ＆ Technology，Nanjing 210044，China;

2Key Laboratory of Meteorological Disaster，Ministry of Education （KLME）／Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters （CIC-FEMD），Nanjing University of Information Science ＆ Technology，Nanjing 210044，China

Abstract This study analyzes the spatial-temporal characteristics and influencing factors of sensible heat flux based on data observed by the Large Aperture Scintillometer （LAS） from April 2018 to April 2019 in Pancheng town.The results indicate the following：1） In the town，the sensible heat flux exhibits a unimodal diurnal variation pattern，significantly larger during the daytime compared to the nighttime.The sensible heat flux is significantly higher on clear days than on cloudy days and slightly lower at night compared to cloudy days.The average hourly sensible heat flux ranges between 2.25—200.53 W·m-2 and 13.10—132.52 W·m-2 for clear and cloudy days，respectively.2） The sensible heat flux is significantly higher in summer than in winter.In August，the values during the daytime and at night are 112.19 and 23.54 W·m-2，respectively.In February，the values during the daytime and at night are 35.57 and 11.57 W·m-2，respectively.3） On clear days with inflows from different directions，notable differences in sensible heat fluxes are observed in the town.As the proportion of impervious layers increases，the proportion of net radiation allocated to sensible heat fluxes also increases significantly.When the proportion exceeds 60％，the upward trend becomes less evident.4） Using the Monin-Obukhov length as a criterion for judging atmospheric stability，the C2n method exhibits a low misjudgment rate among the five methods when using LAS data，making it a more suitable method for determining atmospheric stability at night in the town.5） Among the surface parameters affecting sensible heat flux，changes in effective height have the greatest impact.Changes in wind speed have a larger effect，especially in autumn and winter.The change in Bowens ratio has a minor effect，while the effects of changes in temperature，roughness length，and displacement height are negligible.

Keywords town;LAS;sensible heat flux;spatial-temporal characteristics;environmental factors;surface parameters

doi：10.13878／j.cnki.dqkxxb.20200317001

（責任編輯：劉菲）