華北地區(qū)水位下降是否會(huì)減緩氣溫上升
——淺部地溫影響的數(shù)值模擬分析

尹鳳玲， 張懷， 石耀霖

中國(guó)科學(xué)院計(jì)算地球動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室， 北京 100049

?

華北地區(qū)水位下降是否會(huì)減緩氣溫上升
——淺部地溫影響的數(shù)值模擬分析

尹鳳玲， 張懷， 石耀霖*

中國(guó)科學(xué)院計(jì)算地球動(dòng)力學(xué)重點(diǎn)實(shí)驗(yàn)室， 北京 100049

華北地區(qū)由于長(zhǎng)期持續(xù)的地下水過量開采，導(dǎo)致了大面積地下水位大幅下降，引發(fā)地面塌陷、地下水質(zhì)污染等一系列地質(zhì)環(huán)境問題，這些現(xiàn)象早已為人們所熟知和關(guān)注.然而地下水位下降還會(huì)造成百米量級(jí)淺部地溫及其梯度的變化，因此即使來自地球深部的大地?zé)崃髅芏葲]有變化，年度平均的從表淺部位通過地表實(shí)際傳導(dǎo)進(jìn)入大氣的熱流密度會(huì)減小，這是中外文獻(xiàn)中尚未見討論過的問題.我們通過數(shù)值模擬發(fā)現(xiàn)假定大地?zé)崃髅芏炔蛔兊臈l件下，華北數(shù)萬平方公里地下水位下降會(huì)造成百米尺度內(nèi)的地溫降低，從而傳入大氣的熱流密度降低40%以上，且會(huì)持續(xù)數(shù)百年以上的時(shí)間.這種長(zhǎng)時(shí)間大范圍的傳導(dǎo)入大氣的熱流密度變化對(duì)環(huán)境會(huì)造成什么影響是一個(gè)十分值得關(guān)注的問題.這一預(yù)測(cè)在一定程度上得到了氣象站地溫觀測(cè)數(shù)據(jù)的支持，但由于目前氣象觀測(cè)站只有3.2 m深度范圍內(nèi)的地溫資料，累計(jì)不超過5、60年，中間還有10余年的間斷，而且表淺深度地溫受地表多種因素的影響也較大，這些資料難以對(duì)我們關(guān)心的地下水位下降引起流入大氣的熱流密度變化這一問題提供直接確鑿的數(shù)據(jù)來進(jìn)行分析，因此今后有必要開展對(duì)地下數(shù)十乃至數(shù)百米地溫進(jìn)行持續(xù)精確的監(jiān)測(cè)工作.

地下水； 淺層地溫; 熱流密度; 數(shù)值模擬

1 引言

20世紀(jì)70年代以來，華北地區(qū)由于地下水超采嚴(yán)重和降水等補(bǔ)給不足，地下水位持續(xù)下降，已形成面積達(dá)7×104km2的漏斗群并在繼續(xù)擴(kuò)大(何慶成等，2006)，尤其90年代以后，該區(qū)地下水位以平均1 m·a-1(何慶成等，2006)左右的速率下降. 這一日益嚴(yán)重的狀況引起了國(guó)內(nèi)外專家學(xué)者的關(guān)注.劉中培等(2012)對(duì)石家莊平原區(qū)的地下水位年內(nèi)及年際動(dòng)態(tài)變化進(jìn)行研究，認(rèn)為人工開采和降雨是影響地下水位的兩個(gè)主要因素.劉伯華等(2008)分析了廊坊市地下水下降漏斗的成因，概括出地下水位下降漏斗的緩變性、積累性和不可逆性等特點(diǎn).徐勝?gòu)?qiáng)和劉鳳嶺(2010)對(duì)邯鄲市的淺層漏斗區(qū)地下水現(xiàn)狀及存在的問題進(jìn)行了闡述，并提出了可行的補(bǔ)救措施.王娟等(2009)對(duì)邢臺(tái)地區(qū)的地下水位進(jìn)行動(dòng)態(tài)變化趨勢(shì)分析，指出寧柏隆漏斗區(qū)和滏東平原漏斗區(qū)面積在進(jìn)一步擴(kuò)張, 中心埋深在進(jìn)一步增加.Liu等(2011)整體上論述了整個(gè)華北地區(qū)水資源狀況及存在的問題，并對(duì)其可持續(xù)利用與發(fā)展提出了可行的方法.Shimada等(2002)以日本關(guān)東平原地下水超采情況為參照，認(rèn)為目前華北地下水恢復(fù)工作更加困難并且需要更長(zhǎng)的時(shí)間.

地下水位下降并非單一存在的現(xiàn)象，它可能引發(fā)其他一系列連鎖性的問題，如地面沉降、海水入侵、水質(zhì)污染、生態(tài)環(huán)境的改變和氣候效應(yīng)等(Zektser et al.，2005；Liu et al.，2001；李志敏等，2010).地下水位埋深與土壤含水量呈顯著的負(fù)相關(guān)( 魏彬等, 2013)，土壤的含水量不同其熱學(xué)性質(zhì)亦有所不同(Nobel and Geller，1987；高青等，2005)，因?yàn)楸頊\土壤溫度關(guān)系到農(nóng)作物和植被的生長(zhǎng)，地下水位下降對(duì)近地表土壤溫度會(huì)造成什么樣的影響值得探討，近年來也有學(xué)者對(duì)此做了研究討論.孫緒金等(2003)認(rèn)為地表由于降雨滲入地下形成徑流對(duì)長(zhǎng)期受日照、太陽輻射而形成的高溫的近地表地層起冷卻降溫作用，當(dāng)?shù)叵滤陆禃r(shí)近地表的這一冷卻作用則消失，所以導(dǎo)致近地表土壤溫度會(huì)上升，但他沒有給出理論上的證明或?qū)嶋H觀測(cè)資料的支持.鮑艷等(2007)對(duì)西北干旱半干旱這一嚴(yán)重缺水地區(qū)進(jìn)行氣候分析，認(rèn)為當(dāng)?shù)叵滤陆禃r(shí)，當(dāng)?shù)氐牡乇矸凑章试黾?，地表接受的輻射能增加，致使地表溫度升?傅志敏等 (2006，2010a，2010b)根據(jù)近20年來漏斗地區(qū)表淺地溫有所升高的觀測(cè)事實(shí)，對(duì)地溫升高量與蒸發(fā)量、入滲補(bǔ)給量、降水量、日照、相對(duì)濕度和地下水埋深等多個(gè)物理量進(jìn)行了相關(guān)性分析，也認(rèn)為地下水下降造成表淺土壤溫度升高.兩個(gè)變量變化趨勢(shì)存在相關(guān)，并不能證明二者一定存在物理因果聯(lián)系.他忽略了地下水位不下降的地區(qū)同樣存在表淺地溫上升的現(xiàn)象，并且升溫幅度比水位下降區(qū)更大；忽略了所有地區(qū)均存在氣候變暖、因而地溫也相應(yīng)增高的基本事實(shí)，把表淺地溫升高認(rèn)為是地下水開采造成的，這一結(jié)論是值得商榷的.地溫變化與地下水開采兩者之間的物理關(guān)系需待進(jìn)一步研究.

關(guān)于華北地區(qū)因采補(bǔ)失衡造成大面積地下水漏斗群及其可能引發(fā)的連鎖性問題的研究有很多，但是利用數(shù)值模擬的方法定量分析地下水位下降對(duì)近地表溫度和熱流影響的研究國(guó)內(nèi)外尚未見到.本文從數(shù)值模擬的角度出發(fā)，利用有限元數(shù)值方法，求解非穩(wěn)態(tài)熱傳導(dǎo)方程，在來自深部大地?zé)崃髅芏炔蛔兊臈l件下，分析華北地區(qū)地下水位持續(xù)下降對(duì)淺部地溫和經(jīng)地表流入大氣的熱流的影響，為探討地下水位下降對(duì)氣候和環(huán)境的長(zhǎng)期影響提供新的思考和研究課題.

2 數(shù)值模擬與結(jié)果分析

2.1 地?zé)崃骰咎卣?/p>

我們觀測(cè)到的地下溫度分布主要由以下幾個(gè)因素控制.

第一，來自地球深部的熱量，即大地?zé)崃髅芏?它在人類歷史的時(shí)間尺度內(nèi)一般可以視為恒定的，在地殼表層放射性元素含量很低的情況下，熱導(dǎo)率為常數(shù)時(shí)，它控制的地溫曲線為隨深度增加的直線(如果放射性元素產(chǎn)熱率不容忽視則地溫—深度曲線為拋物線)，如圖1a所示.如果存在兩種水平分層的均勻介質(zhì)，則在各層中地溫曲線為斜率不同的直線，但兩層中的地?zé)崃髅芏?地溫梯度與熱導(dǎo)率之乘積)相同，來自地殼深部的大地?zé)崃髅芏纫话阍?5～90 mW·m-2(Vitorello and Pollack，1980).

第二，地表溫度變化的影響.如果地表溫度按余弦函數(shù)變化T=ΔTcos(ωt)，則在地層熱導(dǎo)率、比熱等熱物理性質(zhì)不變時(shí)，地下深度y，時(shí)刻t的溫度變化T為

(1)

與地表溫度周期變化相聯(lián)系的熱流密度隨深度變化為

(2)

其中，λ為土壤熱導(dǎo)率.可以看出在接近地表處，地表溫度達(dá)到極大值之后，熱量從地表流向地下；在地表溫度達(dá)到極小值后，熱量從地下流向地表.對(duì)于日變化，熱流幅度最大可達(dá)近60 W·m-2；對(duì)于年變化，熱流幅度最大可達(dá)約3 W·m-2.但隨著深度增加這一熱流量也迅速呈指數(shù)減小.

以上討論中的定量表達(dá)式一般均假定熱導(dǎo)率、比熱和密度等為均勻常數(shù)，如果它們時(shí)空分布有變化，特別淺地表由于降水的滲透、含水量的變化和冰凍等，則會(huì)出現(xiàn)更加復(fù)雜的情況.

圖1 大地?zé)崃骱湍甓葰鉁刈兓C合作用下地溫的基本特征(a)來自地球深部的熱量形成的地溫梯度恒定的溫度分布； (b) 地表溫度年周期變化傳播到地下引起的在春夏秋冬四個(gè)相位的地溫分布； (c) 真實(shí)的地下溫度分布，它等于來自地球內(nèi)部的熱量和地表周期性變化傳播到地下形成的溫度的疊加.Fig.1 General characteristics of ground temperature under the combined action of geothermal heat flux and periodical variation of ground surface temperature(a) Temperature distribution produced by heat from the deep Earth with a constant geo-temperature gradient; (b) Temperature-depth curves at four different phases (corresponding to spring, summer, autumn, winter) due to the annual variation of ground surface temperature; (c) Real geo-temperature profile, which can be expressed as the sum of that produced by heat from the deep Earth and that produced by the effects of periodical temperature changes on the ground.

第三，地下水活動(dòng)的影響.地下水的活動(dòng)包括水平運(yùn)動(dòng)和垂直運(yùn)動(dòng)，它們對(duì)淺表地層的地溫變化有著很大的影響，但具體情況復(fù)雜，需要個(gè)別分析.本文將討論地下水位持續(xù)降低的影響.

2.2 數(shù)值模型

本文采用有限元方法，建立數(shù)值模型，計(jì)算地下水位下降時(shí)，由于由水飽和的“濕土”變?yōu)椴煌潭鹊摹案赏痢保瑹釋?dǎo)率和熱容量等熱學(xué)性質(zhì)發(fā)生變化，從而在地表淺層產(chǎn)生的溫度暫態(tài)變化.盡管土壤“干”“濕”的變化在實(shí)際情況下是一個(gè)漸變的過程，但為突出過程的物理特征，本模型把其簡(jiǎn)化為突變過渡.簡(jiǎn)化后的土壤模型便于說明地下水位下降對(duì)土壤溫度和地?zé)崃髅芏鹊挠绊戇@一物理問題的機(jī)理，且對(duì)數(shù)值模擬的結(jié)果的數(shù)量級(jí)也不會(huì)造成大的影響.待以后有了優(yōu)質(zhì)的不同深度地溫觀測(cè)的時(shí)間序列和地下熱學(xué)性質(zhì)的資料，可以用比較精確的模型進(jìn)行更加精確的計(jì)算.

數(shù)值模型取計(jì)算深度為500 m的一個(gè)柱體，初始地下水埋深為5 m，以深部熱流為0.06 mW·m-2(臧紹先等，2002)的穩(wěn)態(tài)溫度剖面作為初始溫度值，在左右兩側(cè)絕熱以及深部熱流和地表溫度不變的邊界條件下，求解瞬態(tài)熱傳導(dǎo)方程，計(jì)算華北地區(qū)水位以平均1 m·a-1左右的速率下降時(shí)“濕土”變?yōu)椤案赏痢边^程中地下的溫度和地?zé)崃鞯膭?dòng)態(tài)變化；然后假定在地下水位埋深降低到55 m后人們采取措施，地下水位保持在55 m埋深不再下降，繼續(xù)計(jì)算這種條件下后續(xù)溫度變化.網(wǎng)格尺寸為0.01 m，時(shí)間步長(zhǎng)為0.01 a.各向同性的二維齊次非穩(wěn)態(tài)熱傳導(dǎo)方程表達(dá)式為

(3)

通過有限元方法數(shù)值求解上述方程，得到地下各點(diǎn)溫度場(chǎng)和熱流場(chǎng)的分布.

數(shù)值模擬中用到土壤的熱導(dǎo)率和體積熱容量?jī)蓚€(gè)熱性質(zhì)參數(shù).土壤的熱導(dǎo)率與其含水量、礦物質(zhì)組成及溫度等因素有關(guān)，其中含水量對(duì)土壤熱導(dǎo)率的影響最大(Becker et al.，1992).水位平均下降速率為1 m·a-1，相對(duì)瞬間抽水來說十分緩慢.這種情況下，水位長(zhǎng)期持續(xù)下降后，土壤釋放出大部分水，土壤的含水量在10%以下(王大純，1980).雖然不同土壤類型其熱導(dǎo)率具體值不同，但土壤熱導(dǎo)率隨含水量的變化趨勢(shì)相似.平均來講，飽和土壤的熱導(dǎo)率一般略大于1.00 W/m℃-1，含水量在10%以下的土壤熱導(dǎo)率為0.5 W/m℃-1左右(Nobel and Geller，1987；Becker et al.，1992；Sepaskhah and Boersma，1979；Bristow et al.，1998；Parikh，1979；Al Nakshabandi and Kohnke，1965；楊淑貞等，1986).土壤的體積熱容量與含水量呈線性相關(guān)，隨含水量線性增加，飽和土壤的體積熱容量為3.50×106J/m3℃-1左右，10%以下的體積熱容量為1.50×106J/m3℃-1左右(Bristow et al.，1998；Abu-Hamdeh，2003).因此本文取干土的體積熱容量為1.64×106J/m3℃-1，熱導(dǎo)率為0.45 W/m℃-1；濕土的體積熱容量3.50×106J/m3℃-1，熱導(dǎo)率為1.00 W/m℃-1.

2.3 數(shù)值模擬結(jié)果及分析

2.3.1 地下水位下降對(duì)地下淺部溫度的影響

初始時(shí)刻地下水埋深為5 m，由于干濕土壤的熱導(dǎo)率不同，地溫的空間分布在地下5 m處分為斜率不同的兩條直線.假設(shè)地下水位持續(xù)下降50年(即下降50 m)后水位不再變化，數(shù)值計(jì)算得到的不同時(shí)間的地溫深度曲線如圖2所示.圖3給出三種不同深度類型的溫度隨時(shí)間的變化：一種是始終位于地下水位之上的“干土”，以3 m深度處的溫度變化為例(圖3a)，地下水位持續(xù)下降會(huì)引起它的溫度逐漸降低，地下水停止下降后溫度會(huì)緩緩回升；另一種是原來位于地下水面之下的“濕土”，然后地下水位下降后變?yōu)椤案赏痢?，?0 m深度處的點(diǎn)為例(圖3b)，開始它的溫度上升，但是地下水持續(xù)下降該點(diǎn)變?yōu)樘幱谒幻嬷系摹案赏痢焙?，則溫度類似3a也要下降，在地下水位停止下降后，才會(huì)逐漸緩緩上升；最后一種是始終處于地下水位以下的“濕土”，以70 m深度的點(diǎn)為例(圖3c)，由于它遠(yuǎn)離原來的地下水面，所以當(dāng)水位開始下降時(shí)還沒有影響到它，溫度幾乎沒有變化，但在數(shù)百年的長(zhǎng)時(shí)間尺度內(nèi)會(huì)緩緩上升.這些層位的溫度降低幅度一般不超過0.3 ℃.這樣微小的降溫在3 m內(nèi)地表氣象地溫記錄中會(huì)被區(qū)域氣候變暖的大背景所掩蓋，在分析時(shí)需格外注意.

地溫下降的物理原因是，原來穩(wěn)態(tài)溫度場(chǎng)的情況下，各層接收到的從下面?zhèn)鬟f上來的熱量與向上傳遞出去的熱量處于動(dòng)態(tài)平衡狀態(tài)，各個(gè)地層保持它們一定的穩(wěn)態(tài)溫度場(chǎng).由于干土層熱導(dǎo)率低，因此穩(wěn)態(tài)時(shí)溫度梯度較高，濕土層熱導(dǎo)率高，穩(wěn)態(tài)時(shí)溫度梯度較低.但是當(dāng)?shù)叵滤唤档秃螅瓉淼臐裢磷優(yōu)楦赏梁?，這一轉(zhuǎn)變層的熱導(dǎo)率降低，阻滯了下面熱量傳入的速率，造成濕土層得到的熱量減少，溫度降低.而新的水位面以下的濕土，由于通過界面?zhèn)鬟f上去熱量減少，積聚的熱量使溫度有微弱升高.這種轉(zhuǎn)變層熱導(dǎo)率降低，接收到的熱量速率減小，無力繼續(xù)為其上面的層位按原有速率提供熱量，是造成地下水位面上層位的溫度降低的物理原因.

圖2 地溫動(dòng)態(tài)變化(a) 水位開始下降0～100 a的40 m深度內(nèi)的地溫動(dòng)態(tài)變化過程；(b) 0～10 a在16 m深度內(nèi)的地溫狀態(tài)變化更細(xì)致的描述.Fig.2 Dynamic subsurface temperatures Dynamic ground temperature as the groundwater level drops(a) over the course of 100 years within a depth of 40 m, and (b) from 0 to 10 years within a depth of 16 m.

圖3 不同層位的溫度動(dòng)態(tài)變化(a) 始終在水位面以上的點(diǎn)(depth=3 m)的溫度動(dòng)態(tài)變化； (b) 開始水位面以上，后來在水位面以下的點(diǎn)(depth=30 m)的溫度動(dòng)態(tài)變化； (c) 始終在水位面以下的點(diǎn)(depth=70 m)的溫度動(dòng)態(tài)變化.Fig.3 Dynamic temperatures at different depths(a) Temperature of the point (depth=3m) always above the water level varies with time; (b) Temperature of the point (depth=30 m) at the beginning above the water level and then below the water level varies with time; (c) Temperature of the point (depth=70 m) always below the water level varies with time.

如果采取措施使水位不再降低，則來自地下深部的熱量逐漸將轉(zhuǎn)變層以上的干土加熱，它們的溫度將緩緩上升，直到達(dá)到新的平衡態(tài)為止.也就是說，始終位于干土層的層位溫度將隨著地下水位下降而下降，直到水位下降停止后一段時(shí)間才會(huì)逐漸緩緩上升.始終位于濕土層的層位溫度將持續(xù)上升.而轉(zhuǎn)變層位的土層，在轉(zhuǎn)變前位于地下水面下的濕土層內(nèi)，溫度上升；轉(zhuǎn)變?yōu)楦赏梁鬁囟认陆?，水位下降停止后?jīng)過少許時(shí)間，再次轉(zhuǎn)為上升，直到上升到新的穩(wěn)態(tài)為止.這種升高重新達(dá)到新的穩(wěn)態(tài)需要上千年的時(shí)間，是一個(gè)持續(xù)很長(zhǎng)時(shí)間的緩慢過程.

2.3.2 地下水位下降對(duì)淺部和流入大氣的熱流密度的影響

初始大地?zé)崃髅芏?0 mW·m-2, 當(dāng)?shù)叵滤怀掷m(xù)下降50年即下降50 m后，從圖4a顯示從地表流入大氣的熱流密度減小了約25 mW·m-2，減小到了原來的60%左右，而且持續(xù)時(shí)間很長(zhǎng)，即使現(xiàn)在保持平均55 m的水位不再降低，恢復(fù)到穩(wěn)態(tài)的地溫梯度也需要上千年的時(shí)間，如圖4b所示.這種長(zhǎng)期持續(xù)的經(jīng)地表流入大氣的熱流密度降低對(duì)氣候和環(huán)境會(huì)有什么影響呢？

地下水位下降50 m時(shí)，地表流入大氣的熱流密度減小了25 mW·m-2，標(biāo)準(zhǔn)大氣壓下空氣比熱為1.0×103J/kg℃-1，密度為1.29 kg·m-3，地表向大氣傳遞的熱流密度減少量可以使與周圍絕熱的底面積為1 m2、高度為7752 m的等密度的大氣柱體(它在地面產(chǎn)生1個(gè)大氣壓力)溫度值一年降低0.08 ℃.而研究區(qū)域目前在全球變暖條件下的氣溫增長(zhǎng)速率每年也不過為0.03 ℃·a-1，因此持續(xù)數(shù)百年的經(jīng)地表流入大氣的熱流密度減小對(duì)氣候產(chǎn)生的影響恐怕不可輕易忽略.當(dāng)然實(shí)際大氣運(yùn)動(dòng)情況非常復(fù)雜，我們不可能單獨(dú)隔離一個(gè)大氣柱體觀測(cè)這種效應(yīng)，但這一計(jì)算結(jié)果分析表明，華北地區(qū)長(zhǎng)期持續(xù)的經(jīng)地表流入大氣的熱流密度降低提出了一個(gè)全新的需要研究的課題.

圖4 經(jīng)地表流入大氣的熱流密度的動(dòng)態(tài)變化Fig.4 Dynamic heat flux from ground surface to the atmosphere

3 實(shí)測(cè)氣象地溫資料分析

我們的理論計(jì)算是否能得到實(shí)際觀測(cè)的佐證呢？答案是比較困難的.我國(guó)缺少數(shù)十米、數(shù)百米地溫的長(zhǎng)期觀測(cè)資料，僅部分地震前兆臺(tái)站近十年開始有這方面的資料.華北一些縣的氣象觀測(cè)站有不深于3.2 m的地溫記錄，雖然記錄從1950或60年代開始，但在1966年到1979年(文化大革命及其后)有14年的中斷.而且這種表淺深度的地溫受地表溫度年變化影響十分顯著，一年以上周期的地表溫度變化可以波及到這一深度.除了氣候變化的影響外，數(shù)十年來中國(guó)城市化過程中，許多氣象臺(tái)的環(huán)境也發(fā)生了很大變化.理論計(jì)算的來自深部大地?zé)崃鞯挠绊懕谎谏w在地表溫度變化的影響之下，因此要精確提取來自深部的信號(hào)是非常困難的.然而，一些氣象地溫?cái)?shù)據(jù)仍然可以看出一些有意思的現(xiàn)象.

根據(jù)國(guó)家氣象信息中心提供的河北省一些氣象觀測(cè)站的1956—2011年共56年的逐年平均氣溫和自建站以來至2011年地下1.6 m、3.2 m深度的年平均地溫?cái)?shù)據(jù)，其測(cè)量精度均為0.1 ℃，我們對(duì)河北省主要漏斗區(qū)1956—2011年的年平均氣溫和相應(yīng)臺(tái)站自建站以來的地溫?cái)?shù)據(jù)進(jìn)行趨勢(shì)分析，如圖5所示，石家莊(a,b)、保定(c,d)、邢臺(tái)(e,f)、滄州(g,h)4個(gè)主要地下水漏斗區(qū)的氣溫分別以0.035 ℃·a-1，0.032 ℃·a-1，0.027 ℃·a-1，0.043 ℃·a-1的年平均速率變暖，4個(gè)觀測(cè)站氣溫年平均增長(zhǎng)速率為0.034 ℃·a-1.4個(gè)地區(qū)的1.6 m地溫的平均年變化速率分別為-0.022 ℃·a-1，0.0029 ℃·a-1，0.0075 ℃·a-1，-0.0006 ℃·a-1，4個(gè)觀測(cè)站平均為-0.0032 ℃·a-1.3.2 m地溫的平均年變化速率分別為-0.0048 ℃·a-1，0.0121 ℃·a-1，0.0163 ℃·a-1，0.0054 ℃·a-1，4個(gè)觀測(cè)站平均為0.0073 ℃·a-1.可見，1.6 m地溫和3.2 m地溫的年增溫速率遠(yuǎn)低于當(dāng)?shù)氐臍鉁卦鰷厮俾?，甚至石家莊地區(qū)地下1.6 m和3.2 m以及邢臺(tái)1.6 m深度出現(xiàn)年平均地溫降低的趨勢(shì).

非主要漏斗區(qū)張家口氣象觀測(cè)站則以0.048 ℃·a-1的年平均速率變暖，1.6 m地溫的平均年變化速率為0.034 ℃·a-1，3.2地溫的平均年變化速率為0.039 ℃·a-1.可見，非主要漏斗區(qū)近地表溫度上升速率較高，且與氣溫增溫速率比較接近(圖6).

在華北非主要地下水漏斗區(qū)近地表溫度上升速率接近氣溫增溫速率，而華北地下水主要漏斗區(qū)近地表溫度上升速率一般低于氣溫增溫速率，甚至在有的臺(tái)站反而出現(xiàn)地溫下降的現(xiàn)象，這是否與近半個(gè)世紀(jì)以來地下水位持續(xù)降低的現(xiàn)象有關(guān)，值得注意.

氣象1.6 m和3.2 m地溫的變化受地表溫度變化的影響較大，如果能有數(shù)十米乃至數(shù)百米深度地溫長(zhǎng)期變化監(jiān)測(cè)的資料將更加有利于理論計(jì)算的驗(yàn)證和分析，但是缺乏這類資料.現(xiàn)在一些地震前兆觀測(cè)臺(tái)站開始注意觀測(cè)井孔內(nèi)的溫度的記錄，我們期待在經(jīng)過數(shù)十年資料積累后，它們能在這類問題分析中發(fā)揮重要作用.

圖5 1956—2011年石家莊(a,b)、保定(c,d)、邢臺(tái)(e,f)和滄州(g,h)氣象觀測(cè)站氣溫和1.6 m、3.2 m的地溫?cái)?shù)據(jù)OAT表示實(shí)測(cè)氣溫值，LFAT表示氣溫線性擬合曲線，OGT表示實(shí)測(cè)地溫值，LFGT表示地溫線性擬合曲線.Fig.5 Air temperature and the ground temperature at depths of 1.6 and 3.2 m at Shijiazhuang (a,b), Baoding (c,d), Xingtai (e,f) and Cangzhou (g,h) Stations from 1956 to 2011OAT stands for observed air temperature, LFAT stands for linear fitting of air temperature, OGT stands for observed geo-temperature, and LFGT stands for linear fitting of geo-temperature.

圖6 1956—2011年張家口氣象觀測(cè)站氣溫和(a) 1.6 m、(b) 3.2 m的地溫?cái)?shù)據(jù)OAT表示實(shí)測(cè)氣溫值，LFAT表示氣溫線性擬合曲線，OGT表示實(shí)測(cè)地溫值，LFGT表示地溫線性擬合曲線.Fig.6 Air temperature and the ground temperature at depths of (a) 1.6 and (b) 3.2 m at Zhangjiakou Station from 1956 to 2011OAT stands for observed air temperature, LFAT stands for linear fitting of air temperature, OGT stands for observed geo-temperature, and LFGT stands for linear fitting of geo-temperature.

4 討論

以上計(jì)算中沒有考慮地下水在向上抽取的過程中，深部溫度較地表平均氣溫高的水?dāng)y帶的一部分熱量.一個(gè)問題是帶上來的這部分熱量是否可以使地表增溫以抵消地?zé)崃髅芏冉档偷睦鋮s效應(yīng)呢？抽水一年攜帶出的熱量為q=cρΔTu≈1.09 W·m-2，其中，水的比熱c=4.2×103J/(kg·℃)，水的密度ρ=103kg·m-3，地表與地下水的平均溫差ΔT=0.825 ℃，抽水速度u=1 m·a-1=3.1536×10-7m·s-1.抽出的水溫度夏季低于氣溫，冬季高于氣溫，如果以均勻速度抽水，年平均來說可能會(huì)向上帶出一些熱量補(bǔ)償因地下水位下降造成的地?zé)崃髅芏葴p少的影響.但由于地下水人工開采量特別是農(nóng)業(yè)灌溉用水各個(gè)季節(jié)不同(冬小麥春灌，以及夏季的玉米、蔬菜等)(Zhang and Pei，2003；Xu et al.，2005)，主要集中于春、夏季的4—7月(毛任釗等，2002)，這段時(shí)間地下水溫低于地表平均溫度，加劇了地表熱流密度的減少，對(duì)空氣起到了冷卻作用，而且抽出的部分水會(huì)流走，因此并不會(huì)抵消地?zé)崃髅芏葴p少的效應(yīng).

目前的計(jì)算表明，地表溫度的日變化會(huì)在1 m深度內(nèi)產(chǎn)生很大的熱擾動(dòng)，白天地表高溫出現(xiàn)后熱量會(huì)從地表向地下傳導(dǎo)，夜間地表低溫出現(xiàn)后熱量會(huì)從地下向地表傳達(dá)，雖然其熱流密度峰值很大，可達(dá)60 W·m-2，但影響深度不超過1 m；對(duì)于年變化，夏季高溫峰值后熱量會(huì)從地表向下傳導(dǎo)，冬季低溫峰值后熱量會(huì)從地下向地表傳導(dǎo)，熱流密度峰值熱流幅度也不小，最大可達(dá)約3 W·m-2，但影響深度不超過數(shù)米；這些周期性的熱流變化雖然很大，但長(zhǎng)時(shí)間平均起來，傳下去的和傳上來的熱量恰恰抵消.相反大地地?zé)崃髁棵芏群苄?，僅有60 mW·m-2，抽水造成地下水位下降造成的經(jīng)地表傳入大氣的熱流密度的減小量?jī)H有25 mW·m-2，但是它是長(zhǎng)年累月大面積持續(xù)作用的.

地?zé)崃髅芏鹊淖兓欠駮?huì)影響氣候變化，這是一個(gè)困難的問題，因?yàn)榇蟮責(zé)崃髅芏攘考?jí)太小，在地表完全被淹沒在地表溫度周期變化造成的表淺熱流的振蕩之中，雖然在數(shù)學(xué)分析中我們可以區(qū)分出地?zé)崃鞯姆至?，但是在?shí)際現(xiàn)象的觀測(cè)中很難區(qū)分出各個(gè)組分的作用.然而我們前面提到過，地下水位下降造成經(jīng)地表傳入大氣的熱流密度減小25 mW·m-2，在與周圍絕熱條件下，足以使單位面積上造成1個(gè)大氣壓力的空氣柱體一年溫度降低0.08 ℃，而漏斗地區(qū)在全球變暖條件下觀測(cè)到的氣溫增長(zhǎng)速率每年也不過為0.03 ℃·a-1，這種作用應(yīng)該是不能忽略的.

值得注意的是，湯懋蒼等(湯懋蒼等，1981；吳士杰等，1982；胡澤勇等，1989；周亞軍等，1989；陳文海等，2001；賈效祿等，2011)從20世紀(jì)70年代以來致力于地溫、地?zé)岙惓?duì)氣候的影響的研究，經(jīng)過多年的思考和實(shí)踐，他們根據(jù)地?zé)岙惓５淖兓?guī)律，用地下如3.2 m深度的地溫距平場(chǎng)分布來進(jìn)行汛期旱澇預(yù)報(bào)，多年來汛期降水預(yù)報(bào)結(jié)果和實(shí)際旱澇分布吻合得較好(鄧國(guó)衛(wèi)等，2011；湯懋蒼等，1986，1988，1994).得到了他人的很高評(píng)價(jià)，趙紅巖等(2003)提到“在短期氣候預(yù)報(bào)中，湯懋蒼等利用前冬(12—2月)地溫距平場(chǎng)，預(yù)測(cè)我國(guó)汛期(6—8月)降水，已有近30年的歷史，預(yù)測(cè)水平在各種方法中名列前茅”.盡管其預(yù)報(bào)的機(jī)理仍然缺乏定量的理論，但地溫和氣候變化之間可能存在物理聯(lián)系，應(yīng)該是客觀事實(shí).還有一些學(xué)者探索過大地?zé)崃髟趨^(qū)域生態(tài)環(huán)境和氣候的演化過程中扮演著重要角色.匡耀求等(2002a)通過對(duì)比世界各地主要地理單元的大地?zé)崃魈卣髋c生態(tài)環(huán)境特征，詳細(xì)論述了大地?zé)崃髟诘乇憝h(huán)境和生態(tài)演變的重要影響，認(rèn)為大地?zé)崃鞯母叩蛯?duì)生態(tài)環(huán)境的優(yōu)劣有密切的關(guān)系，在一定程度上影響了區(qū)域氣候的演變.對(duì)我國(guó)西部地區(qū)(匡耀求等，2002b，2003)、鄂爾多斯盆地(匡耀求等，2004)以及滇池地區(qū)(匡耀求，2002)進(jìn)行區(qū)域性分析，進(jìn)一步論證了區(qū)域大地?zé)崃髋c生態(tài)環(huán)境之間的密切關(guān)系，大地?zé)崃魇怯绊憛^(qū)域生態(tài)環(huán)境和氣候的重要因素.

因此，華北地區(qū)由地下水位下降引起的地溫及地溫梯度和經(jīng)地表傳入大氣的熱流密度變化對(duì)區(qū)域性氣候和環(huán)境會(huì)造成什么樣的潛在影響，是一個(gè)值得注意的新的研究課題.這類研究需要對(duì)地溫、土壤含水情況、不同條件下土壤的熱物理性質(zhì)等有長(zhǎng)期的觀測(cè)和積累，這些都是長(zhǎng)期艱巨的任務(wù).

5 結(jié)論

1) 利用有限元數(shù)值模擬方法定量分析地下水下降對(duì)淺地表溫度的影響表明，雖然來自地球深部的大地?zé)崃髅芏葲]有變化，但是地下水位下降，會(huì)引起淺部地下溫度的降低以及接近地表熱流密度的減少.其物理原因是，大地?zé)崃鳜F(xiàn)在有相當(dāng)部分用于加熱淺部被擾動(dòng)的(低于穩(wěn)態(tài)應(yīng)該有的溫度的)土層，因此真正能夠經(jīng)過地表流入大氣的熱量小于深部大地?zé)崃髅芏?

2) 在華北漏斗區(qū)，經(jīng)地表傳入大氣的熱流密度減少雖然絕對(duì)量級(jí)不大，僅僅25 mW·m-2，但是其相對(duì)變化量卻不小，達(dá)到原來地?zé)崃髅芏鹊?0%.而且即使現(xiàn)在保持平均55 m的水位不再降低，其影響也將持續(xù)數(shù)百年乃至上千年的時(shí)間，只有在千年尺度才有可能緩緩恢復(fù)到原來的穩(wěn)態(tài)地溫梯度狀態(tài).數(shù)萬平方公里大面積，數(shù)百年長(zhǎng)以上時(shí)間的流入大氣的熱流密度降低，會(huì)對(duì)氣候變化產(chǎn)生什么影響是一個(gè)重要的問題.

3) 影響氣候變化的因素很多，氣溫變化反之又影響氣象地溫觀測(cè)結(jié)果，因此這個(gè)研究是一個(gè)困難的問題.加強(qiáng)對(duì)數(shù)十米到數(shù)百米的深度的地溫長(zhǎng)期監(jiān)測(cè)，對(duì)于了解地溫變化特征具有重要意義，今后應(yīng)該在地震前兆臺(tái)站觀測(cè)中注意積累這方面資料.并對(duì)地下巖石、土層的熱學(xué)性值做進(jìn)一步研究，以便進(jìn)一步完善數(shù)值模型研究.

致謝 感謝國(guó)家自然基金(NSFC41474067)和中國(guó)科學(xué)院、國(guó)家外國(guó)專家局創(chuàng)新團(tuán)隊(duì)國(guó)際合作伙伴計(jì)劃(KZZD-EW-TZ-19)資助，感謝國(guó)家氣象信息中心提供的數(shù)據(jù)，感謝外審專家提出的寶貴意見.

Abu-Hamdeh N H. 2003. Thermal properties of soils as affected by density and water content.BiosystemsEngineering, 86(1): 97-102.

Al Nakshabandi G, Kohnke H. 1965. Thermal conductivity and diffusivity of soils as related to moisture tension and other physical properties.AgriculturalMeteorology, 2(4): 271-279.Bao Y, Li Y H, Chen R S, et al. 2007. A suppositional experiment about the impact of groundwater drawdown on regional climate change in Northwest China.AridZoneResearch(in Chinese), 24(4): 434-440. Becker B R, Misra A, Fricke B A. 1992. Development of correlations for soil thermal conductivity.InternationalCommunicationsinHeatandMassTransfer, 19(1): 59-68.

Bristow K L. 1998. Measurement of thermal properties and water content of unsaturated sandy soil using dual-probe heat-pulse probes.AgriculturalandForestMeteorology, 89(2): 75-84.

Chen W H, Liu Y X, Tang M C. 2001. A “Diqi” map method of qualitative prediction in annual precipitation with geothermal anomaly in China.PlateauMeteorology(in Chinese), 20(4): 429-434.

Deng G W, Tang M C, Guo W D. 2011. Analysis on evolution of earth-atmosphere system in China from January 2009 to April 2010.PlateauMeteorology(in Chinese), 30(5): 1274-1278.

Fu Z M, Zhou Z F, Li M W, et al. 2006. Catastrophe criterion for ground temperature based on effect of groundwater table.JournalofJiangsuUniversity(NaturalScienceEdition) (in Chinese), 27(5): 441-445.

Fu Z M, Xiang Y. 2010. Brief analysis of heat island effect of shallow groundwater depression cone and relevant concepts.WaterResourcesandPower(in Chinese), 28(12): 30-32.

Fu Z M, Xiang Y, Zhou Z F. 2010. Critical groundwater depth influencing the surface air temperature.YellowRiver(in Chinese), 32(12): 7-9. Gao Q, Li M, Yu M, et al. 2005. Heat transfer influenced by containing water in wet soil.JournalofThermalScienceandTechnology(in Chinese), 4(2): 136-140.

He Q C, Liu W B, Li Z M. 2006. Land subsidence survey and monitoring in the North China Plain.GeologicalJournalofChinaUniversities(in Chinese), 12(2): 195-209.

Hu Z Y, Xiang W G, Tang M C. 1989. Geothermal anomaly as an important factor to climate anomaly.∥ Comprehensive Study of Cosmic, Earth and Biosphere (in Chinese). Beijing: Science and Technology of China Press, 257-260.

Jia X L, Yang Q, Zhao Y, et al. 2011. Characteristics of soil temperature at deep layers in Urumqi in latest 50 years.JournalofAridMeteorology(in Chinese), 29(2): 201-204.

Kuang Y Q. 2002. The effects of terrestrial heat flow on the environment of Lake Dianchi and new approach to control its pollution.BulletinofMineralogy,PetrologyandGeochemistry(in Chinese), 21(2): 86-91.Kuang Y Q, Huang N S, Zhu Z Y, et al. 2002a. Influence of the terrestrial heat flow on the evolution of the surface ecosystem and environment.GeologyinChina(in Chinese), 29(1): 86-95.

Kuang Y Q, Huang N S, Zhu Z Y, et al. 2002b. Terrestrial heat flow——an important Natural factor influencing evolution of ecological system and environment in the Western China.BulletinofMineralogy,PetrologyandGeochemistry(in Chinese), 21(1): 30-34.

Kuang Y Q, Huang N S, Wu Z F, et al. 2003. Effects of terrestrial heat flow on the environment and ecologic evolution of Western China and its significance.AdvanceinEarthSciences(in Chinese), 18(1): 22-29.Kuang Y Q, Huang N S, Zhu Z Y, et al. 2004. Change of terrestrial heat flow in Ordos Basin since late Mesozoic and its effects on the pattern and evolution of the ecologic environment.BulletinofMineralogy,PetrologyandGeochemistry(in Chinese), 23(4): 318-325.Li Z M, Sun B H. 2010. Analysis and discussion of the relationship between groundwater exceeding extraction and land subsidence in Cangzhou.GroundWater(in Chinese), 32(3): 35-36.

Liu B H, Chen B J, Gong J Z, et al. 2008. Cause analysis of groundwater funnel in Langfang of Hebei province.AridMeteorology(in Chinese), 26(2), 53-56.

Liu C M, Yu J J, Kendy E. 2001. Groundwater exploitation and its impact on the environment in the North China Plain.WaterInternational, 26(2): 265-272.

Liu J, Cao G L, Zheng C M. 2011. Sustainability of groundwater resources in the North China Plain. ∥Sustaining Groundwater Resources. Netherlands: Springer, 69-87. Liu Z P, Wang F Q, Yu F R. 2012. Variation of shallow groundwater level in Shijiazhuang Plain.South-to-NorthWaterTransfersandWaterScience&Technology(in Chinese), 10(5): 124-127.

Mao R Z, Liu X J, Lou H J. 2002. A study on the dynamic variation of shallow groundwater table in the saline and water deficit region in North China Plain.ProgressinGeography(in Chinese), 21(6): 564-572.

Nobel P S, Geller G N. 1987. Temperature modelling of wet and dry desert soils.TheJournalofEcology, 75(1): 247-258.

Parikh R J, Havens J A, Scott H D. 1979. Thermal diffusivity and conductivity of moist porous media.SoilScienceSocietyofAmericaJournal, 43(5): 1050-1052.

Pollack H N, Huang S, Shen P Y. 1998. Climate change record in subsurface temperatures: a global perspective.Science, 282(5387): 279-281.

Pollack H N, Huang S P. 2000. Climate reconstruction from subsurface temperatures.AnnualReviewofEarthandPlanetarySciences, 28(1): 339-365.Sass J. 1991. Climate plumbs the depths.Nature, 349(6309): 458.

Sepaskhah A R, Boersma L. 1979. Thermal conductivity of soils as a function of temperature and water content.SoilScienceSocietyofAmericaJournal, 43(3): 439-444. Shimada J, Tang C, Tanaka T, et al. 2002. Irrigation caused groundwater drawdown beneath the North China Plain.∥Proceedings of Int. Groundwater Conf., Darwin, Australia, 1-7.Sun X J, Zhou G, Sun X P, et al. 2003. The study of ground temperature change after groundwater exceeding extraction.JournalofNorthChinaInstituteofWaterConservancyandHydroelectricPower(in Chinese). 24(1): 37-39.

Tang M C, Wu S J. 1981. The prediction of precipitation trend in flood season and with soil temperature at deep layers.Meteorology(in Chinese), (8): 11-12.

Tang M C, Yin J H, Cai J P. 1986. A statistical relationship between the field of soil temperature in winter and of the precipitation in spring and flood season.PlateauMeteorology(in Chinese), 5(1): 40-52.

Tang M C, Zhang J, Wang J X, et al. 1988. The similarity between the seasonal anomalous maps of soil temperature and the precipitation of the subsequent season.ActaMeteorologicaSinica(in Chinese), 46(4): 481-485.

Tang M C, Zhang J. 1994. Seasonal mean soil temperature anomaly field at depth 3. 2m and its application in prediction for flood season.PlateauMeteorology(in Chinese), 13(2): 178-187.

Turcotte D L, Schubert G. 2002. Geodynamics. 2nd ed. Cambridge: Cambridge University Press, 150-152.

Vitorello I, Pollack H N. 1980. On the variation of continental heat flow with age and the thermal evolution of continents.JournalofGeophysicalResearch:SolidEarth(1978—2012), 85(B2): 983-995.

Wang D C, Zhang R Q, Shi Y H, et al. 1995. Fundamentals of Hydrogeology (designed for Hydrogeological and Engineering Geology). Beijing:Geological Publishing House, 8-13.

Wang J, Yan J P, Du J W. 2009. The tendency of groundwater dynamic variation and regional response in typical area of North China Plain.JournalofNaturalScienceofHunanNormalUniversity(in Chinese), 32(4): 102-107.

Wang J Y. 1992. Inversing climate change from subsurface temperatures——a frontier topic of contemporary theoretical geothermics.QuaternarySciences(in Chinese), (1): 36-39.

Wei B, Yimiti H, Wang Q F, et al. 2013. Correlation of burial depth of groundwater and soil water content in the Keriya Oasis, Xinjiang, China.JournalofDesertResearch(in Chinese), 33(4): 1110-1116.

Wu S J, Tang M C, Zhong Q. 1982. A thermodynamic-statistical model for the long-range rainfall forecasting at the flood season in China.PlateauMeteorology(in Chinese), 1(2): 18-24.

Xu S Q, Liu F L. 2010. The present state, existing problems and prevention measures of groundwater water funnel in Handan City.GroundWater(in Chinese), 32(4): 76-77.

Xu Y Q, Mo X G, Cai Y L, et al. 2005. Analysis on groundwater table drawdown by land use and the quest for sustainable water use in the Hebei Plain in China.AgriculturalWaterManagement, 75(1): 38-53.

Yang S Z, Zhang W R, Shen X J. 1986. Experimental conductivity of research on the thermal water-saturated porous rocks.ActaPetrologicaSinica(in Chinese), 2(4): 83-91.

Zang S X, Liu Y G, Ning J Y. 2002. Thermal structure of the lithosphere in North China.ChineseJ.Geophys. (in Chinese), 45(1): 56-66.

Zektser S, Loáiciga H A, Wolf J T. 2005. Environmental impacts of groundwater overdraft: selected case studies in the southwestern United States.EnvironmentalGeology, 47(3): 396-404.

Zhang X Y, Pei D, Hu C S. 2003. Conserving groundwater for irrigation in the North China Plain.IrrigationScience, 21(4): 159-166.

Zhao H Y, Xue W X, Zhu X N. 2003. The study on dryness and wetness prediction by earth-atmosphere chart in Northwest China.AridMeteorology(in Chinese), 21(4): 37-40.

Zhou Y J, Tang M C. 1989. The basic method of forecasting yearly precipitation with soil temperature.PlateauMeteorology(in Chinese), 8(3): 228-240.

附中文參考文獻(xiàn)

鮑艷, 李耀輝, 陳仁升等. 2007. 地下水位下降對(duì)區(qū)域氣候影響的虛擬試驗(yàn). 干旱區(qū)研究, 24(4): 434-440.

陳文海, 柳艷香, 湯懋蒼. 2001. 利用地?zé)岙惓６ㄐ灶A(yù)測(cè)年度降水之方法. 高原氣象, 20(4): 429-434.

鄧國(guó)衛(wèi), 湯懋蒼, 郭維棟. 2011. 2009年1月—2010年4月我國(guó)地氣系統(tǒng)的演變分析. 高原氣象, 30(5): 1274-1278.

傅志敏, 周志芳, 李明武等. 2006. 基于地下水埋深效應(yīng)的地溫突變判據(jù). 江蘇大學(xué)學(xué)報(bào) (自然科學(xué)版), 27(5): 441-445.

傅志敏, 向衍. 2010. 簡(jiǎn)析淺層地下水降落漏斗區(qū)熱島效應(yīng)及有關(guān)概念. 水電能源科學(xué), 28(12): 30-32.

傅志敏, 向衍, 周志芳. 2010. 影響近地表氣溫的地下水臨界埋深. 人民黃河, 32(12): 7-9.

高青, 李明, 于鳴等. 2005. 濕土壤含濕特性對(duì)傳熱影響研究. 熱科學(xué)與技術(shù), 4(2): 136-140.

何慶成, 劉文波, 李志明. 2006. 華北平原地面沉降調(diào)查與監(jiān)測(cè). 高校地質(zhì)學(xué)報(bào), 12(2): 195-209.

胡澤勇, 向衛(wèi)國(guó), 湯懋蒼. 1989. 地?zé)岙惓Ｊ菤夂虍惓５闹匾蜃?∥天地生綜合研究. 北京: 中國(guó)科學(xué)技術(shù)出版社, 257-260.

賈效祿, 楊青, 趙勇等. 2011. 近50a烏魯木齊市的深層地溫特征. 干旱氣象, 29(2): 201-204.

匡耀求. 2002. 大地?zé)崃鲗?duì)滇池環(huán)境的影響與滇池污染治理的新思路. 礦物巖石地球化學(xué)通報(bào), 21(2): 86-91.

匡耀求, 黃寧生, 朱照宇等. 2002a. 試論大地?zé)崃鲗?duì)地表環(huán)境與生態(tài)演變的影響. 中國(guó)地質(zhì), 29(1): 86-95.

匡耀求, 黃寧生, 朱照宇等. 2002b. 大地?zé)崃鳌绊懳鞑凯h(huán)境與生態(tài)演變的重要自然因素. 礦物巖石地球化學(xué)通報(bào), 21(1): 30-34.

匡耀求, 黃寧生, 吳志峰等. 2003. 大地?zé)崃鲗?duì)中國(guó)西部環(huán)境與生態(tài)演變的影響及其研究意義. 地球科學(xué)進(jìn)展, 18(1): 22-29.

匡耀求, 黃寧生, 胡振宇等. 2004. 鄂爾多斯盆地中生代晚期以來大地?zé)崃鞯淖兓捌鋵?duì)生態(tài)環(huán)境格局和演變的影響. 礦物巖石地球化學(xué)通報(bào), 23(4): 318-325.

李志敏, 孫炳華. 2010. 滄州地下水超采與地面沉降關(guān)系的分析與探討. 地下水, 32(3): 35-36.

劉伯華, 陳寶江, 宮進(jìn)忠等. 2008. 廊坊市地下水位下降漏斗的成因. 干旱氣象, 26(2): 53-56.

劉中培, 王富強(qiáng), 于福榮. 2012. 石家莊平原區(qū)淺層地下水位變化研究. 南水北調(diào)與水利科技, 10(5): 124-127.

毛任釗, 劉小京, 婁華君. 2002. 華北平原缺水鹽漬區(qū)淺層地下水位動(dòng)態(tài)分析. 地理科學(xué)進(jìn)展, 21(6): 564-572.

孫緒金, 周革, 孫曉平等. 2003. 地下水超采引起的近地表地溫變化規(guī)律研究. 華北水利水電學(xué)院學(xué)報(bào), 24(1): 37-39.

湯懋蒼, 吳士杰. 1981. 用深層地溫預(yù)報(bào)汛期降水趨勢(shì). 氣象, (8): 11-12.

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

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

湯懋蒼, 張建. 1994. 季平均3. 2m地溫距平場(chǎng)在汛期預(yù)報(bào)中的應(yīng)用. 高原氣象, 13(2): 178-187.

王大純, 張人權(quán), 史毅虹等. 1995. 水文地質(zhì)學(xué)基礎(chǔ) (水文地質(zhì)及工程地質(zhì)專用). 北京: 地質(zhì)出版社, 8-13.

王娟, 延軍平, 杜繼穩(wěn). 2009. 華北平原典型區(qū)地下水動(dòng)態(tài)變化趨勢(shì)及其區(qū)域響應(yīng). 湖南師范大學(xué)自然科學(xué)學(xué)報(bào), 32(4): 102-107.

汪集旸. 1992. 根據(jù)地溫資料推斷氣候變化——當(dāng)代理論地?zé)嵫芯康囊粋€(gè)前沿課題. 第四紀(jì)研究, (1): 36-39.

魏彬, 海米提·依米提, 王慶峰等. 2013. 克里雅綠洲地下水埋深與土壤含水量的相關(guān)性. 中國(guó)沙漠, 33(4): 1110-1116.

吳士杰, 湯懋蒼, 鐘強(qiáng). 1982. 汛期降水預(yù)報(bào)的統(tǒng)計(jì)—熱力學(xué)模式. 高原氣象, 1(2): 18-24.

徐勝?gòu)?qiáng), 劉鳳嶺. 2010. 邯鄲市淺層漏斗區(qū)地下水水質(zhì)現(xiàn)狀、存在問題及防治措施. 地下水, 32(4): 76-77.

楊淑貞, 張文仁, 沈顯杰. 1986. 孔隙巖石熱導(dǎo)率的飽水試驗(yàn)研究. 巖石學(xué)報(bào), 2(4): 83-91.

臧紹先, 劉永剛, 寧杰遠(yuǎn). 2002. 華北地區(qū)巖石圈熱結(jié)構(gòu)的研究. 地球物理學(xué)報(bào), 45(1): 56-66.

趙紅巖, 薛萬孝, 祝小妮. 2003. 應(yīng)用 “地氣圖” 預(yù)測(cè)西北旱澇的初步研究. 干旱氣象, 21(4): 37-40.

周亞軍, 湯懋蒼. 1989. 用地溫做年度降水預(yù)報(bào)的一種改進(jìn). 高原氣象, 8(3): 228-240.

(本文編輯 汪海英)

Persistent drawdown of groundwater table in North China may reduce local climate warming rate：Numerical simulation and analysis of the impacts on shallow ground temperature

YIN Feng-Ling, ZHANG Huai, SHI Yao-Lin*

UniversityofChineseAcademyofSciences,KeyLaboratoryofComputationalGeodynamicsofChineseAcademyofSciences,Beijing100049,China

Due to long-term and persistent groundwater over-exploitation, the water table has been found significantly declined across an area of approximately 70000 km2in North China. Consequently, a series of geo-environmental problems have emerged and drawn a large amount of public attention. These include land subsidence, groundwater contamination, etc. We found that decline of groundwater level can result in significant variation of temperature and its gradient within depths of hundreds of meters. Therefore, the annual average heat flux conducted from shallow ground to atmosphere can be reduced, even though there is no change in heat flux from the deep Earth. This effect has yet discussed neither in domestic nor in foreign studies. Our numerical simulation shows that assuming the heat flux from deep Earth is constant, the depression of the water table leads to reduction of ground temperature within depths of hundreds of meters, and therefore the reduction of heat flux transferred to atmosphere as much as 40% in groundwater drawdown zone in North China after 50 years over-exploitation. It is small the amount of reduction of heat flux transferred from surface to atmosphere. However, it can last for a long period of hundreds of years, and cover a large area of 70000 km2. Its effect on climate may not be negligible. Meteorological temperature (up to 3.2 m depth) records show similar trend. Although air temperatures are influenced by many factors, the decline of heat flux transferred to atmosphere from shallow ground should be one of the factors to be studied. We suggest that more attention be paid to this effect, and systematic monitoring of geotherm at depths of several dozen to several hundred meters should be carried out.

Groundwater; Shallow ground temperature; Heat flow density; Numerical simulation

10.6038/cjg20151018.

Yin F L, Zhang H, Shi Y L. 2015. Persistent drawdown of groundwater table in North China may reduce local climate warming rate: Numerical simulation and analysis of the impacts on shallow ground temperature.ChineseJ.Geophys. (in Chinese),58(10):3649-3659,doi:10.6038/cjg20151018.

國(guó)家自然科學(xué)基金(NSFC41474067)和中國(guó)科學(xué)院，國(guó)家外國(guó)專家局創(chuàng)新團(tuán)隊(duì)國(guó)際合作伙伴計(jì)劃(KZZD-EW-TZ-19)資助.

尹鳳玲，女，1984年生，博士研究生，主要從事地球動(dòng)力學(xué)研究.E-mail:yin_fengling@126.com

*通訊作者 石耀霖，E-mail：shiyl@ucas.ac.cn

10.6038/cjg20151018

P314

2014-02-17，2015-09-17收修定稿

尹鳳玲， 張懷， 石耀霖. 2015. 華北地區(qū)水位下降是否會(huì)減緩氣溫上升——淺部地溫影響的數(shù)值模擬分析.地球物理學(xué)報(bào)，58(10):3649-3659,