王霞 劉德璽 李永濤 王振猛 馬丙堯

摘要：為進(jìn)一步研究黃河三角洲鹽堿地區(qū)土壤-植被-大氣連續(xù)體的作用機(jī)理并為熱量平衡研究提供基礎(chǔ)數(shù)據(jù)，本試驗(yàn)以黃河三角洲鹽堿地區(qū)長(zhǎng)期生長(zhǎng)的刺槐人工林為對(duì)象，利用黃河三角洲森林生態(tài)系統(tǒng)定位研究站2016年度觀測(cè)數(shù)據(jù)，研究了該地區(qū)刺槐人工林土壤熱通量的變化特征及其與地表凈輻射的關(guān)系。結(jié)果表明：10、20 cm土層土壤熱通量日變化基本呈現(xiàn)“S”型變化，且隨土壤深度增加振幅減小;兩層土壤熱通量的月累積總量在3—8月為正值，其它月份為負(fù)值;兩層的年總土壤熱通量分別為-23.85、-12.97 MJ/m2，占地表凈輻射（720.07 MJ/m2）的比例分別為-3.31%、-1.80%，土壤均表現(xiàn)為熱源;在日尺度、月尺度水平上，土壤熱通量與地表凈輻射分別呈極顯著和顯著正相關(guān);土壤熱通量對(duì)地表凈輻射的反饋存在延時(shí)效應(yīng)，延時(shí)5 h的10 cm土層土壤熱通量和延時(shí)9 h的20 cm土層土壤熱通量與地表凈輻射相關(guān)性最高。本研究得出土壤熱通量與地表凈輻射關(guān)系密切，這對(duì)研究土壤-植被-大氣連續(xù)體中地氣能量交換和熱量平衡具有重要意義。

關(guān)鍵詞：黃河三角洲;刺槐;土壤熱通量;凈輻射

中圖分類號(hào)：S152.8 ?文獻(xiàn)標(biāo)識(shí)號(hào)：A ?文章編號(hào)：1001-4942（2020）01-0088-06

Abstract In order to further study the action mechanism of soil-plant-atmosphere continuum and provide basic data for studying the heat balance in the Yellow River Delta of saline-alkali area， based on the 2016 annual data of forest ecosystem positioning research station of the Yellow River Delta， the long-term plantation of Robinia pseudoacacia was taken as research objects to study the variation characteristics of soil heat flux and its relationship with surface net radiation. The results suggested that the diurnal variations of soil heat flux at 10- and 20-cm depth mostly showed S shapes， and the amplitude decreased with the depth of soil layers. The monthly accumulation of soil heat fluxes of both layers was positive from March to August and negative in the other months. The annual total soil heat fluxes at 10-cm and 20-cm depth were -23.85 and -12.97 MJ/m2 respectively with the proportion of total soil heat fluxes to the surface net radiation（720.07 MJ/m2） as -3.31% and -1.80% respectively， indicating that the soil of both layers was heat source at the annual temporal scale. There were extremely significant correlations between soil heat fluxes and the surface net radiation at daily temporal scale， while significant correlations at monthly temporal scale. There was time-lag effect in the feedback of soil heat flux to surface net radiation， the soil heat fluxes lagged 5 h at 10-cm depth and 9 h at 20-cm depth had the most correlation with the surface net radiation. In this study， we concluded that the soil heat flux was in close relationship with the surface net radiation， which had great significance to studying the energy exchange between ground and atmosphere among soil-plant-atmosphere continuum and heat balance.

Keywords Yellow River Delta;Robinia pseudoacacia;Soil heat flux;Net radiation

土壤作為植物賴以生存的重要環(huán)境要素，是土壤-植被-大氣連續(xù)體（soil-plant-atmosphere continuum，簡(jiǎn)稱SPAC）的基礎(chǔ)。土壤熱通量是表征地-氣能量交換的主要物理參量，是森林生態(tài)系統(tǒng)能量平衡方程中的重要組成部分，與系統(tǒng)的能量閉合程度息息相關(guān)[1]。研究土壤熱通量的變化特征，對(duì)于分析森林生態(tài)系統(tǒng)的能量平衡以及區(qū)域小氣候、土壤氣候和土壤生態(tài)環(huán)境的形成有重要意義。目前，國(guó)內(nèi)有關(guān)森林生態(tài)系統(tǒng)土壤熱通量的研究，主要集中于寒帶[2]、寒溫帶[3，4]和亞熱帶[1，5-10]等地區(qū)，研究對(duì)象有原始林、人工純林和混交林等，而針對(duì)黃河三角洲鹽堿地區(qū)森林生態(tài)系統(tǒng)土壤熱通量的研究還較少。

刺槐（Robinia pseudoacacia Linn.），因其較強(qiáng)的耐鹽堿能力和適應(yīng)性，成為黃河三角洲鹽堿地區(qū)的主要代表樹種。近年來(lái)對(duì)該地區(qū)刺槐人工林土壤的研究多集中于土壤基本理化性質(zhì)、土壤酶活性和土壤微生物特征[11-14]，還有土壤水分生態(tài)特征[15]、土壤水鹽動(dòng)態(tài)變化[16，17]等方面，而關(guān)于熱通量的研究還未見報(bào)道。本試驗(yàn)以黃河三角洲鹽堿地區(qū)長(zhǎng)期生長(zhǎng)的人工刺槐林為研究對(duì)象，探討不同土層土壤熱通量的變化特征，分析土壤熱通量與地表凈輻射的關(guān)系，旨在為進(jìn)一步研究該區(qū)域人工刺槐林SPAC的作用機(jī)理和熱量平衡提供參考。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)區(qū)域位于山東省東營(yíng)市河口區(qū)，屬暖溫帶半濕潤(rùn)大陸性季風(fēng)氣候。年均氣溫12.9 ℃，年均無(wú)霜期234 d、凍土期44 d。全年降水分布不均，年均降水量500～600 mm，其中夏季降水約占全年降水量的70%，蒸降比為3.5∶ 1。地下水埋深2～3 m。土壤類型以鹽化潮土和濱海鹽土為主。土壤鹽分組成以氯化物為主。該區(qū)主要造林樹種有刺槐（Robinia pseudoacacia Linn.）、絨毛白蠟（Fraxinus velutina Torr.）、白榆（Ulmus pumila Linn.）、臭椿[Ailanthus altissima （Mill.）Swingle]、國(guó)槐（Sophora japonica Linn.）等，林下植被主要有翅堿蓬[Suaeda salsa （Linn.）Pall.]、白茅[Imperata cylindrica （Linn.）Beauv.]、小飛蓬[Conyza canadensis （Linn.）Cronq.]、狗牙根[Cynodon dactylon （Linn.）Pers.]、藜（Chenopodium album Linn.）、羅布麻（Apocynum venetum Linn.）等。

研究地點(diǎn)為黃河三角洲森林生態(tài)系統(tǒng)定位研究站孤島輔站，位于濟(jì)南軍區(qū)黃河三角洲綜合訓(xùn)練基地十二分場(chǎng)刺槐林內(nèi)。刺槐林為1985年春人工栽植，株行距為2 m×3 m，平均胸徑18.74 cm，平均樹高8.95 m。

1.2 研究方法

人工刺槐林內(nèi)設(shè)有25 m綜合梯度觀測(cè)塔（包含森林小氣候觀測(cè)系統(tǒng)）。其中，地下10 cm、20 cm處設(shè)有HFP01（SC）土壤熱通量傳感器（Hukseflux，Netherlands），用于測(cè)定土壤熱通量;凈輻射傳感器CNR4（Kipp&Zonen，Netherlands）設(shè)置在距地面1.5 m處，用于測(cè)定地表凈輻射能量。所有儀器測(cè)得數(shù)據(jù)由數(shù)據(jù)采集器CR1000（Campbell，USA）自動(dòng)記錄，數(shù)據(jù)采集頻率為10 min，24 h不間斷采集。

1.3 數(shù)據(jù)處理與分析

本研究選取2016年1月1日—2016年12月31日的觀測(cè)數(shù)據(jù)，以0∶ 00為次日劃分點(diǎn)。利用Microsoft Excel 2013軟件進(jìn)行數(shù)據(jù)統(tǒng)計(jì)、制圖和分析。

2 結(jié)果與分析

2.1 土壤熱通量的日變化特征

熱量由表層或淺層土壤向深層土壤傳輸時(shí)，土壤熱通量為正值，即土壤表現(xiàn)為熱匯;熱量由深層土壤向淺層或表層土壤傳輸時(shí)，土壤熱通量為負(fù)值，即土壤表現(xiàn)為熱源。選擇1、4、7、10月作為代表月份，選取代表月份逐日每半小時(shí)的瞬時(shí)數(shù)據(jù)，每月相同時(shí)刻數(shù)據(jù)進(jìn)行平均處理，研究土壤熱通量的日變化特征（圖1）。

整體來(lái)說(shuō)，10、20 cm土層土壤熱通量日變化規(guī)律基本一致。其中，10 cm土層土壤熱通量呈典型的“S”型變化趨勢(shì)，波動(dòng)較大，振幅大于20 cm土層，最小值出現(xiàn)在7∶ 00—9∶ 30，最大值出現(xiàn)在15∶ 30—17∶ 30。20 cm土層土壤熱通量（1月除外）也呈“S”型變化，但波動(dòng)程度低于10 cm土層，振幅較小，說(shuō)明日尺度上深層土壤的地-氣能量交換強(qiáng)度較淺層土壤小;而1月份20 cm土層土壤熱通量變化較平緩，波動(dòng)較小，說(shuō)明該時(shí)段大氣與該層次土壤熱量交換強(qiáng)度很小;20 cm土層最小值出現(xiàn)在9∶ 30—12∶ 30，最大值出現(xiàn)在20∶ 00—21∶ 00，與10 cm土層相比較，極值出現(xiàn)時(shí)間有一定延遲，這是因?yàn)闊崃吭谕寥纼?nèi)傳遞需要一定的時(shí)間（圖1）。

由圖1可以看出，10、20 cm土層土壤熱通量1月均為負(fù)值，平均值分別為-7.29、-6.07 W/m2，兩層土壤表現(xiàn)為熱源，說(shuō)明兩層土壤分別向其上層土壤釋放熱量;4月，10、20 cm土層土壤熱通量平均值分別為5.07、4.99 W/m2，總體來(lái)說(shuō)該時(shí)段兩層土壤表現(xiàn)為熱匯，說(shuō)明兩層土壤分別向其下層土壤傳遞熱量;7月，10、20 cm土層土壤熱通量均為正值，平均值分別為4.77、5.09 W/m2，兩層土壤表現(xiàn)為熱匯，說(shuō)明兩層土壤分別向其下層土壤傳遞熱量;10月，10、20 cm土層土壤熱通量平均值分別為-4.64、-3.92 W/m2，總體來(lái)說(shuō)該時(shí)段兩層土壤表現(xiàn)為熱源，說(shuō)明兩層土壤分別向其上層土壤釋放熱量。

2.2 土壤熱通量的月變化特征

一年中，10、20 cm土層土壤熱通量月累計(jì)總量呈現(xiàn)一致的月變化特征。1、2月和9—12月10、20 cm土層土壤熱通量累積量均為負(fù)值，兩層土壤都表現(xiàn)為熱源，說(shuō)明兩層土壤分別向其上層土壤釋放熱量;3—8月為正值，兩層土壤都表現(xiàn)為熱匯，說(shuō)明兩層土壤分別向其下層土壤傳遞熱量;10 cm土層土壤熱通量累積量負(fù)向絕對(duì)值的最大值出現(xiàn)在11月（-23.18 MJ/m2），正向最大值出現(xiàn)在4月（13.57 MJ/m2），絕對(duì)值最小值出現(xiàn)在8月（2.80 MJ/m2），該層年總熱通量累積量為-23.85 MJ/m2，在年尺度上土壤表現(xiàn)為熱源;20 cm土層土壤熱通量累積量負(fù)向絕對(duì)值的最大值出現(xiàn)在11月（-20.48 MJ/m2），正向最大值出現(xiàn)在7月（14.07 MJ/m2），絕對(duì)值最小值出現(xiàn)在9月（-1.56 MJ/m2），該層年總熱通量累積量為-12.97 MJ/m2，其絕對(duì)值小于10 cm土層，在年尺度上土壤表現(xiàn)為熱源（圖2）。

參 考 文 獻(xiàn)：

[1] 胡兆永，陳書軍，龐軍柱，等. 秦嶺火地塘林區(qū)油松林土壤熱通量[J]. 東北林業(yè)大學(xué)學(xué)報(bào)，2014，42（6）：53-57.

[2] 李明財(cái)，羅天祥，郭軍，等. 藏東南高山林線冷杉原始林土壤熱通量[J]. 山地學(xué)報(bào)，2008，26（4）：490-495.

[3] 王美蓮，崔學(xué)明，韓鵬，等. 大興安嶺原始林區(qū)土壤熱通量變化特征的初探[J]. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào)，2010，31（4）：139-142.

[4] 韋昌雷，王立中，劉學(xué)爽，等. 黑龍江大興安嶺興安落葉松林土壤熱通量研究[J]. 防護(hù)林科技，2014（5）：5-7.

[5] 張利平，趙仲輝. 會(huì)同杉木人工林土壤熱通量特征[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào)，2010，30（5）：12-17.

[6] 孫成，江洪，陳健，等. 亞熱帶毛竹林土壤熱通量變異特征[J]. 土壤學(xué)報(bào)，2013，50（5）：966-973.

[7] 陳云飛，江洪，周國(guó)模，等. 覆蓋經(jīng)營(yíng)雷竹林的土壤熱通量季節(jié)變化特征[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào)，2013，30（2）：165-171.

[8] 孫成，江洪，劉玉莉，等. 安吉毛竹林生長(zhǎng)季土壤熱通量的變化特征[J]. 土壤通報(bào)，2014，45（3）：590-594.

[9] 王旭，周國(guó)逸，張德強(qiáng)，等. 南亞熱帶針闊混交林土壤熱通量研究[J]. 生態(tài)環(huán)境，2005，14（2）：260-265.

[10] 尹光彩，王旭，周國(guó)逸，等. 鼎湖山針闊混交林土壤熱狀況研究[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2006，27（3）：16-20.

[11] 馬風(fēng)云，白世紅，侯本棟，等. 黃河三角洲退化人工刺槐林地土壤特征[J]. 中國(guó)水土保持科學(xué)，2010，8（2）：74-79.

[12] 白世紅，馬風(fēng)云，李樹生，等. 黃河三角洲不同退化程度人工刺槐林土壤酶活性、養(yǎng)分和微生物相關(guān)性研究[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2012，20（11）：1478-1483.

[13] 夏江寶，許景偉，李傳榮，等. 黃河三角洲退化刺槐林不同改造方式對(duì)土壤酶活性及理化性質(zhì)的影響[J]. 水土保持通報(bào)，2012，32（5）：171-175，181.

[14] 韓躍. 黃河三角洲混交林生長(zhǎng)規(guī)律和土壤特性研究[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2015.

[15] 夏江寶，許景偉，李傳榮，等. 黃河三角洲退化刺槐林地的土壤水分生態(tài)特征[J]. 水土保持通報(bào)，2010，30（6）：75-80.

[16] 董海鳳，杜振宇，馬丙堯，等. 黃河三角洲人工林地土壤的水鹽動(dòng)態(tài)變化[J]. 水土保持學(xué)報(bào)，2013，27（5）：48-53.

[17] 董海鳳. 黃河三角洲長(zhǎng)期人工林地土壤特征與改良[D]. 泰安：山東農(nóng)業(yè)大學(xué)，2014.

[18] 張文斌. 塔克拉瑪干沙漠土壤熱通量變化特征分析[D]. 烏魯木齊：新疆師范大學(xué)，2016.

[19] 張宏，胡波，劉廣仁，等. 中國(guó)土壤熱通量的時(shí)空分布特征研究[J]. 氣候與環(huán)境研究，2012，17（5）：515-522.

[20] 牛曉棟，劉曉靜，劉世榮，等.亞熱帶-暖溫帶過(guò)渡區(qū)天然櫟林的能量平衡特征[J].生態(tài)學(xué)報(bào)，2018，38（18）：6701-6711.

[21] 陽(yáng)伏林，周廣勝. 內(nèi)蒙古溫帶荒漠草原能量平衡特征及其驅(qū)動(dòng)因子[J]. 生態(tài)學(xué)報(bào)，2010，30（21）：5769-5780.

[22] 王文杰，崔崧，劉瑋，等. 落葉松人工林土壤熱傳導(dǎo)季節(jié)變化及其與環(huán)境因子的關(guān)系[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2008，19（10）：2125-2131.