黃明陽(yáng) 李媛 雍佳

摘要 利用位于南京北郊的南京信息工程大學(xué)觀測(cè)基地2010年各層（0、5、10、15、20、40、80、160、320 cm）地溫?cái)?shù)據(jù)，采用統(tǒng)計(jì)分析方法對(duì)南京北郊地溫的變化特征進(jìn)行了分析，揭示南京市浦口地區(qū)地溫與時(shí)間、深度的變化關(guān)系。結(jié)果表明，南京北郊地區(qū)表層地溫以及淺層地溫有明顯的季節(jié)變化和日變化，呈近似于正弦曲線(xiàn)的變化趨勢(shì)。從地表到20 cm深時(shí)，地溫的日變化逐漸減弱；當(dāng)深度在40～320 cm時(shí)，地溫日變化已不明顯。在垂直方向上，各層地溫日變化幅值隨著深度增加而減小，隨著深度的增加峰值出現(xiàn)的時(shí)間逐漸滯后。

關(guān)鍵詞 地溫；變化特征；差分進(jìn)化算法；垂直結(jié)構(gòu)

中圖分類(lèi)號(hào) P468.0+21 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2017）03-0222-04

Abstract By using ground temperature data of each layer（0， 5， 10， 15， 20， 40， 80， 160， 320 cm） in 2010 observed by Nanjing University of Information Science and Technology located in northern suburb， the variation characteristics of ground temperature were analyzed by statistical analysis method， the variation of ground temperature with time and depth in Pukou Region， Nanjing City was revealed. The results showed that， there were obvious seasonal and diurnal variations of surface temperature and shallow ground temperature in northern suburb of Nanjing， the change trend was similar to the sine curve. From land surface to 20 cm depth， the daily change of ground temperature decreased gradually； from 40-320 cm， the daily variation of ground temperature was not obvious. In vertical direction， the diurnal variation of ground temperature in each layer decreased with the increase of depth， the appearance of peaks gradually delayed with depth increasing.

Key words Ground temperature； Change characteristics； Differential evolution algorithm； Vertical structure

地溫即土壤溫度，是土壤熱量狀況的綜合表征指標(biāo)[1]。地溫受太陽(yáng)輻射、大氣循環(huán)、降水等因素作用，與大氣進(jìn)行著能量的交換，形成了地溫周期性的日變化和季節(jié)變化[2]。研究地溫的變化規(guī)律，有助于揭示土壤對(duì)農(nóng)業(yè)生產(chǎn)和生態(tài)環(huán)境的影響。同時(shí)，陸地下墊面過(guò)程也是引起氣候變化的重要因子之一。因此，研究地溫的變化規(guī)律也能為天氣、氣候的研究提供依據(jù)。

前人對(duì)地溫已進(jìn)行了大量研究，并取得一定的成果[3-6]，如Zhang等[3]研究20世紀(jì)加拿大土壤溫度發(fā)現(xiàn)，土壤溫度對(duì)氣溫和降水復(fù)雜的響應(yīng)過(guò)程會(huì)對(duì)氣候變化產(chǎn)生明顯的影響。史學(xué)正等[7]利用1951—1980年我國(guó)672個(gè)氣象站點(diǎn)40和80 cm年均土壤溫度數(shù)據(jù)，以土壤溫度內(nèi)插法估算了這些站點(diǎn)50 cm深度的年均土壤溫度狀況，并應(yīng)用于中國(guó)土壤系統(tǒng)分類(lèi)；其次，通過(guò)分析土壤溫度與環(huán)境因子（經(jīng)緯度、海拔）的相關(guān)關(guān)系，建立多元回歸方程對(duì)土壤溫度進(jìn)行預(yù)測(cè)。目前，關(guān)于地溫的研究越來(lái)越多，但針對(duì)南京地區(qū)的研究還不多。筆者利用中國(guó)氣象局綜合大氣探測(cè)基地（南京）觀測(cè)站的地溫資料，采用統(tǒng)計(jì)分析方法對(duì)南京北郊地溫的變化特征進(jìn)行了分析，揭示南京市浦口地區(qū)地溫與時(shí)間、深度的變化關(guān)系，為浦口地區(qū)天氣預(yù)報(bào)、氣候分析、農(nóng)業(yè)氣象、作物生產(chǎn)、新型環(huán)保能源開(kāi)發(fā)、政府決策等提供參考。

1 資料與方法

中國(guó)氣象局綜合探測(cè)培訓(xùn)實(shí)習(xí)基地（簡(jiǎn)稱(chēng)觀測(cè)基地）位于南京市浦口區(qū)盤(pán)城鎮(zhèn)南京信息工程大學(xué)西苑校園內(nèi)，占地1.3萬(wàn)m2，觀測(cè)場(chǎng)地及其周?chē)h(huán)境符合有關(guān)規(guī)定。該研究利用該測(cè)站2010年1月1日—12月31日深度為0、5、10、15、20、40、80、160、320 cm每小時(shí)測(cè)得的地溫值。數(shù)據(jù)中有147次記錄缺測(cè)（缺測(cè)記錄為散點(diǎn)隨機(jī)離散分布），可分析時(shí)次為8 613次。

該研究利用差分進(jìn)化算法（Differential Evolution）將2010年地溫?cái)?shù)據(jù)進(jìn)行擬合：

y=u+Jj=1ajsin（bjx+cj）（1）

式中，y為某一層地溫值（℃）；x為時(shí)間（h）；u為常數(shù)；系數(shù)aj、bj和cj分別代表其諧波成分的振幅、頻率與相位，其中由頻率bj可得到曲線(xiàn)變化周期Tj=1/bj。在擬合的基礎(chǔ)上，可計(jì)算擬合曲線(xiàn)與實(shí)際曲線(xiàn)的相關(guān)系數(shù)（擬合度）。J為實(shí)際計(jì)算中的諧波成分的個(gè)數(shù)。若取J=1，則僅進(jìn)行1個(gè)諧波成分的分析；取J=2，則同時(shí)進(jìn)行2個(gè)諧波成分的分析。

2 結(jié)果與分析

2.1 地溫的季節(jié)變化特征

總體來(lái)看（圖1），各深度層地溫均有很明顯的季節(jié)變化，夏季地溫明顯高于冬季，春秋季為單調(diào)上升（或下降）趨勢(shì)。這是由于我國(guó)處于北半球，太陽(yáng)輻射冬天弱、夏天強(qiáng)，導(dǎo)致地溫也隨之變化。隨著深度的增加，地溫曲線(xiàn)在逐漸平滑， 表明深層地溫受地表影響逐漸減小。0～80 cm地溫，2月份及其前后期間存在明顯的、連續(xù)幾天的低值區(qū)域，變化幅度較小，沒(méi)有明顯的日變化，這是由于持續(xù)的多云天氣造成的；6—7月存在成段的“V”型區(qū)，此時(shí)南京正處于江淮梅雨或多云連續(xù)陰雨天氣期間，氣溫偏低，日變化也較緩；9月初出現(xiàn)一個(gè)明顯的高峰區(qū)，此時(shí)南京處于“秋老虎”天氣，之后氣溫下降，地溫總體趨勢(shì)也隨之下降。

2.2 地溫的日變化特征

由圖2可知，地面溫度日變化幅度最大，越往深層，地溫日變化幅度越小。距地表較近（0、5、10、15 cm）的地溫有明顯的日變化，1天只有1個(gè)峰值和1個(gè)谷值；而20、40、80、160、320 cm地溫日變化已不明顯。地表的地溫日變化一般在12：00—13：00達(dá)到峰值，隨著土壤深度的增加，達(dá)到峰值的時(shí)間點(diǎn)逐漸推移，具有時(shí)間滯后性。距地表較近的地溫曲線(xiàn)達(dá)到溫度峰值和谷值的時(shí)間不同，這是由于一天中太陽(yáng)輻射變化導(dǎo)致地溫也隨之變化。秋季地表溫度日變化幅度最大，達(dá)16.85 ℃；冬季地表溫度日變化幅度最小，僅為10.13 ℃。

2.3 地溫的定量分析 從表1可看出，J=1時(shí)，a1sin（b1x+c1）對(duì)應(yīng)季節(jié)變化規(guī)律，季節(jié)變化周期為年周期，不是365 d，而是360～375 d，這可能是由于該研究?jī)H取2010年數(shù)據(jù)的隨機(jī)性造成的。J=2時(shí)，在a1sin（b1x+c1）的基礎(chǔ)上，加入了a2sin（b2x+c2），對(duì)應(yīng)日變化。日變化周期為24 h，這與常識(shí)一致。

從地表到深層相關(guān)系數(shù)r1變大，擬合效果變好，擬合曲線(xiàn)與實(shí)際曲線(xiàn)的相似度更高，可見(jiàn)各層地溫受季節(jié)變化影響，并且由于深層受到日變化或其他因素影響程度較小，則深層地溫變化更符合正弦規(guī)律。加上了日變化過(guò)程，0 cm地溫?cái)M合曲線(xiàn)與實(shí)際曲線(xiàn)的相關(guān)系數(shù)有所提高，說(shuō)明加上日變化后0 cm地溫的擬合效果變好，表明地表溫度既有季節(jié)變化，又有一定的日變化規(guī)律。J=2時(shí)的相關(guān)系數(shù)較J=1時(shí)有明顯提高。而在進(jìn)一步擬合過(guò)程中，當(dāng)取J=3時(shí)，即在a1sin（b1x+c1）+a2sin（b2x+c2）的基礎(chǔ)上，加入了a3sin（b3x+c3），分析是否存在第3種重要的變化規(guī)律影響地溫變化時(shí)，不能得到一個(gè)較穩(wěn)定的b3，同時(shí)r3也提高得很小，所以該研究未列出J=3時(shí)的擬合結(jié)果。由此推斷，南京北郊地溫主要是2種變化方式，由周期判斷為季節(jié)變化和日變化。

當(dāng)深度≤10 cm時(shí)，地溫有明顯或較明顯的季節(jié)變化和日變化。而隨著深度的增加，相關(guān)系數(shù)差（r2-r1）在減小，a2sin（b2x+c2）項(xiàng)的加入使得地溫的擬合曲線(xiàn)與實(shí)際曲線(xiàn)的相似度有所提高，表明隨著深度的增加，a2sin（b2x+c2）項(xiàng)的絕對(duì)值在減小，地溫的日變化幅度在減小。在深度≥15 cm后，深層次的土壤并不能分析出一個(gè)固定的頻率，所以在表1中不再列出。從圖2可以看出，15和20 cm的地溫有日變化趨勢(shì)，但取J=2難以得到較穩(wěn)定的b2，r2也增加不明顯。這可能是由于15和20 cm地溫的日變化程度對(duì)于季節(jié)變化程度來(lái)說(shuō)過(guò)小，擬合時(shí)日變化波動(dòng)難以從季節(jié)變化中區(qū)分出來(lái)。為了證明這一猜想，該研究進(jìn)行了各層地溫日變化和季節(jié)變化程度的對(duì)比，結(jié)果發(fā)現(xiàn)（表2），各層地溫的日極差與年極差之比隨著深度的增加逐漸減小，即15和20 cm地溫日變化確實(shí)對(duì)于整年數(shù)據(jù)而言影響程度在變小，導(dǎo)致擬合過(guò)程中區(qū)分出日變化的難度提高，所以該研究擬合過(guò)程中沒(méi)有擬合出15 cm以下地溫的日變化。

2.4 地溫的垂直結(jié)構(gòu)

從圖2可以看出，隨著深度的增加，不同深度土壤的溫度日變化峰值出現(xiàn)的時(shí)間在后移。將2010年1月2日的地溫?cái)?shù)據(jù)按式（1）進(jìn)行擬合，結(jié)果發(fā)現(xiàn)（表3），不論實(shí)際測(cè)量還是擬合結(jié)果，地溫峰值出現(xiàn)的時(shí)間隨著深度的增加逐漸后移；初相位在變小，對(duì)應(yīng)的曲線(xiàn)峰值也將后移。同時(shí)，實(shí)測(cè)值的峰值與谷值之差（即極值差）隨深度增加也在減小，表明地溫受太陽(yáng)輻射產(chǎn)生的日變化幅度在減小。從擬合結(jié)果也可以看出，相關(guān)系數(shù)在逐漸增加，由于地溫受除太陽(yáng)輻射以外的因素影響在減小，正弦趨勢(shì)逐漸提高。

3 小結(jié)

（1）南京北郊地區(qū)表層地溫以及淺層地溫有明顯的季節(jié)變化和日變化，呈近似于正弦曲線(xiàn)的變化趨勢(shì)。

（2）從地表到20 cm深時(shí)，地溫的日變化逐漸減弱；當(dāng)深度在40～320 cm時(shí)，地溫日變化已不明顯。

（3）在垂直方向上，各層地溫日變化幅值隨著深度增加而減小，隨著深度的增加峰值出現(xiàn)的時(shí)間逐漸滯后。

參考文獻(xiàn)

[1] 李興榮，胡非，孫向明，等.深圳冬季多層地溫日變化特征[J].中國(guó)科學(xué)院研究生院學(xué)報(bào)，2009，26（3）：323-329.

[2] BOND-LAMBERTY B， WANG C K， GOWER S T. Spatiotemporal measurement and modeling of standlevel boreal forest soil temperatures[J].Agricultural and forest meteorology， 2005， 131：27-40.

[3] ZHANG Y，CHEN W J， SMITH S L， et al.Soil temperature in Canada during the twentieth century：Complex responses to atmospheric climate change[J/OL]. Journal of geophysical research，2005，110（3）[2016-09-07].http：//www.researchgate.net/publication/228761679_Soil_tempertature_in_Canade_during_the_twentieth_century_Complex_responses_to_atmospheric_climate_change.DOI.10.1029/2004JD004910.

[4] 湯懋蒼，尹建華，蔡潔萍.冬季地溫分布與春、夏降水相關(guān)的統(tǒng)計(jì)分析[J].高原氣象， 1986， 5（1）：40-52.

[5] 湯懋蒼，張建，王敬香，等.我國(guó)季平均的0.8 m地溫距平場(chǎng)與后一季降水的相關(guān)分析[J].氣象學(xué)報(bào)，1988， 46（4）：481-485.

[6] 湯懋蒼，張建.季平均3.2 m地溫距平場(chǎng)在汛期預(yù)報(bào)中的應(yīng)用[J].高原氣象，1994，13（2）：178-187.

[7] 史學(xué)正，鄧西海.我國(guó)土壤溫度狀況[M]//龔子同.中國(guó)土壤系統(tǒng)分類(lèi)進(jìn)展.北京：科學(xué)出版社，1993：353-360.