本文摘要：摘要：多年凍土約占地球陸地面積的1/4，凍結(jié)層上水作為一種多年凍土區(qū)的特殊水文地質(zhì)現(xiàn)象，其分布受局地因素的控制，且隨活動層的季節(jié)性凍融而變化，影響著地表水和地下水循環(huán)以及多年凍土環(huán)境水熱平衡。多年凍土將凍結(jié)層上水限制在一個狹窄的空間內(nèi)，在暖季，凍結(jié)層上
摘要：多年凍土約占地球陸地面積的1/4，凍結(jié)層上水作為一種多年凍土區(qū)的特殊水文地質(zhì)現(xiàn)象，其分布受局地因素的控制，且隨活動層的季節(jié)性凍融而變化，影響著地表水和地下水循環(huán)以及多年凍土環(huán)境水熱平衡。多年凍土將凍結(jié)層上水限制在一個狹窄的空間內(nèi)，在暖季，凍結(jié)層上水橫向和豎向的滲流傳熱將加劇多年凍土的退化，也會對上覆工程構(gòu)筑物的穩(wěn)定運營造成極大威脅。目前關于凍結(jié)層上水的研究主要集中在分布特征、變化規(guī)律、流量計算、滲流模擬、水熱耦合等方面，以區(qū)域分布特征為出發(fā)點，對凍結(jié)層上水的研究現(xiàn)狀進行了歸納和總結(jié)，并對其工程影響有關的滲流攜熱理論研究成果進行了梳理，最終對今后需要進行深入研究的方向進行了展望，這有助于全面理解凍結(jié)層上水在凍土區(qū)水文過程中的功能，為掌握凍結(jié)層上水的研究現(xiàn)狀提供了進一步的理論參考。
關鍵詞：凍結(jié)層上水;水熱耦合;水文地質(zhì);水量平衡;活動層
0引言
凍結(jié)層上水是指分布于多年凍土層上、活動層中的潛水，存在于平坦的分水嶺、平原和河流階地，隨著氣溫的變化而產(chǎn)生季節(jié)性凍融[1]。多年凍土層作為隔水層或弱透水層，很大程度上阻礙了凍結(jié)層上水的下滲和排泄，同時限制了地表水與凍土層下部融土層之間的水分交換和水力聯(lián)系，決定了凍結(jié)層上水的分布格局和循環(huán)方式[2]。
凍結(jié)層上水的形成、演化、運移和動力過程在水溫、滲流量等方面都具有不同于常規(guī)地區(qū)潛水的特殊性，不僅是寒區(qū)能水循環(huán)中的一個關鍵組成部分，而且與寒區(qū)生態(tài)環(huán)境變化密切相關，在多年凍土區(qū)的陸地水循環(huán)中有著特殊作用[3]。凍結(jié)層上水在含水層厚度、水量、水溫、相態(tài)、水化學特征及水動力特征等方面均具有季節(jié)性變化特性，且受局地因素如地形、地貌、植被、降雨、巖性等的控制而產(chǎn)生區(qū)域性變化。
凍結(jié)層上水和多年凍土相互制約、相互作用，使寒區(qū)凍土水文地質(zhì)條件具有特殊性：隨著活動層的凍結(jié)，凍結(jié)層上水的排泄占整個徑流消退的60%~80%[4-5];隨著活動層的融化，凍結(jié)層上水是多年凍土與大氣系統(tǒng)、植被系統(tǒng)、地表水進行水熱交換的重要紐帶，其空間分布控制著包氣帶含水量的變化特征和地表植被的分布類型[3]。全球變暖導致的凍土退化將增大凍結(jié)層上含水層厚度，勢必影響到凍結(jié)層上水的分布特征和水量變化，但因其動態(tài)過程的復雜性和觀測研究的諸多困難，尚缺乏對其運動規(guī)律、驅(qū)動因素和機制的系統(tǒng)認知[6-7]。
本文對凍結(jié)層上水的研究現(xiàn)狀進行了綜述，有助于系統(tǒng)地認知和掌握凍結(jié)層上水在多年凍土環(huán)境中的功能。從凍結(jié)層上水的賦存特征切入，分析了各類凍結(jié)層上水的分布、運移、補給以及排泄特征。以此為基礎，根據(jù)全球多年凍土的分布范圍，分別歸納了低緯度多年凍土區(qū)(青藏高原、祁連山)和高緯度多年凍土區(qū)(大興安嶺、俄羅斯以及北美)典型區(qū)域凍結(jié)層上水的分布特征。此外，還概括了凍結(jié)層上水對于多年凍土區(qū)既有工程構(gòu)筑物穩(wěn)定性的影響研究，并對凍結(jié)層上水在水循環(huán)系統(tǒng)中水量平衡和運移過程中滲流攜熱效應的理論研究和數(shù)值計算進行了總結(jié)和整理。
1凍結(jié)層上水的賦存特征
凍結(jié)層上水的定義是根據(jù)含水層與多年凍土層之間的空間關系而確定的[8-9]。1932年，俄羅斯學者將多年凍土區(qū)地下水分為凍結(jié)層上水、凍結(jié)層間水和凍結(jié)層下水，將凍結(jié)層上水又細分為凍結(jié)層上層滯水、季節(jié)性融化層水和凍結(jié)層上地下水3種[10]。第一種凍結(jié)層上水即為活動層中的水;第二種中，部分凍結(jié)層上水被季節(jié)性低溫層覆蓋;第三種凍結(jié)層上水為凍結(jié)層上多年融區(qū)的水。
隨著對其認知的發(fā)展，寒區(qū)地下水的分類問題仍然沒有完全解決，但對凍結(jié)層上水的形成和發(fā)展過程已有明確解釋。凍結(jié)層上在冷季發(fā)生凍結(jié)，當?shù)乇黹_始融化后，隨著融化厚度的增加，從土體內(nèi)排出的水分將漸增并聚集到上限之上，便形成了一定厚度的潛流層(即為凍結(jié)層上水)[11]。凍結(jié)層上水的凍結(jié)溫度在-0.5~0℃之間，其補給主要來源于大氣降水和融雪融水，其排泄一部分以下降泉的方式溢出地表并匯流成溪，一部分通過沼澤濕地垂直蒸發(fā)排泄，另一部分則通過融區(qū)補給凍結(jié)層下水。
由于凍結(jié)層上水的產(chǎn)流和補給都發(fā)生在活動層，因此其補給面積和分布面積基本重合[1]。多年凍土對于凍結(jié)層上水具有阻隔、蓄水調(diào)節(jié)和蒸發(fā)抑制的作用[12-13]，目前主要通過定性研究來分析二者的相互關系。多年凍土層產(chǎn)生的水文效應影響到了水循環(huán)動態(tài)規(guī)律和地下水的補給排泄關系，使得降水、地表水、地下水失去了循環(huán)通道，對其三者之間的轉(zhuǎn)換關系造成了制約[14]。
多年凍土的阻隔作用導致了土體導熱導水性能的改變，直接影響著凍結(jié)層上水的下滲和蒸發(fā)[15]，多年凍土的退化將增加地表水和凍結(jié)層上水流量，其中包括凍結(jié)層上水水量、地下水和地表水的交換量、基流流量的增大。相反地，凍結(jié)層上水隨著地形的變化在多年凍土層上部產(chǎn)生側(cè)向流動，流動過程所引起的傳熱效應將造成下伏多年凍土的退化[16-17]。總而言之，凍結(jié)層上水與多年凍土在熱交換、地表水-地下水動態(tài)循環(huán)及區(qū)域水文效應等方面有著密切聯(lián)系。
2低緯度多年凍土區(qū)
青藏高原被稱為“世界屋脊”、“亞洲水塔”、“世界第三極”，是我國重要的生態(tài)安全屏障，也是全球氣候變暖的“先兆區(qū)”和“放大器”[18-19]。青藏高原地區(qū)是世界上中低緯度海拔最高、面積最大的多年凍土區(qū)，在目前氣候背景下，多年凍土面積約為120×104km2，隨著高原平均氣溫的升高，青藏高原多年凍土將發(fā)生退化，在氣溫年增加0.02℃和0.052℃兩種不同氣候背景下50a后多年凍土面積分別減少8.8%和13.4%[20-21]。
此外，受地質(zhì)構(gòu)造、地形條件、水文特征和工程影響，青藏高原地區(qū)自然災害分布呈高強度、高頻率、突發(fā)性、季節(jié)性、準周期性特征[22]。而凍結(jié)層上水作為地下水的一種，成為了多年凍土區(qū)典型的水文地質(zhì)現(xiàn)象，在青藏高原凍土退化及災害發(fā)生過程中有著不可忽視的影響。青藏高原凍結(jié)層上水分布范圍廣，含水層以第四紀松散層、山區(qū)基巖風化破碎帶為主[23]。地下水的補給、流動和排泄均受到多年凍土的控制，同時，也受局地因素如地形、地貌、構(gòu)造、地表水文的影響，因此，每個水文地質(zhì)單元中的凍結(jié)層上水具有差異性。此外，由于多年凍土層厚度隨高程的增大而加厚，這進一步改變了地下水動態(tài)變化規(guī)律，使得凍結(jié)層上水在含水層厚度、流速、水溫、水文化學、動力學等方面都具有季節(jié)變化[24-26]。
凍結(jié)層上水含水層厚度受與凍土上限有關，在青藏高原地區(qū)一般為1~3m，尤其是河流和小湖塘附近，含水層厚度較大，儲水量較高。含水層的空間分布特征明顯受到局部地形的控制，地勢低洼地帶儲水量較大，在地勢較高地區(qū)受徑流及蒸發(fā)的影響，因儲水量低而無法形成統(tǒng)一的含水層。青藏高原地區(qū)含水層的形成從每年4月土層中的冰體融化開始，一直到9月份厚度達到最大[27]，凍土層上水的循環(huán)周期一般為1~4a[3]。
目前對于青藏高原地區(qū)凍結(jié)層上水的研究主要集中在三個區(qū)域：青藏工程走廊，黃河源區(qū)和祁連山地區(qū)。青藏工程走廊內(nèi)由于集中了5條重大線性工程而備受關注，凍結(jié)層上水的存在不但對工程穩(wěn)定性有著不可忽視的影響，還對走廊沿線凍土環(huán)境特征的演化過程產(chǎn)生著不同程度的影響[28-29];黃河源區(qū)作為黃河流域的水源地，分布有許多小盆地、徑流和湖泊，豐富的地表水和地下水使其在多年凍土區(qū)具有特殊的水文地質(zhì)特征[30-31]，作為流域的重要補給來源，凍結(jié)層上水在黃河源區(qū)的水循環(huán)過程中起著調(diào)控作用;祁連山地區(qū)分布有大量內(nèi)陸水系，下墊面特征復雜、冰川眾多、儲水量較大，凍結(jié)層上水廣泛發(fā)育，在水分疏導、地下水交互方面具有不可替代的作用。
2.1青藏工程走廊
與其他區(qū)域類似，青藏工程走廊內(nèi)的凍結(jié)層上水埋藏和分布狀況也受地形、地貌、地層、凍土上限深度等因素的制約，補給來源為大氣降水、融雪、冰川融水[32]。廊內(nèi)山區(qū)地帶如昆侖山區(qū)、可可西里山區(qū)、風火山區(qū)、開心嶺山區(qū)、唐古拉山區(qū)、頭二九山區(qū)基本無凍結(jié)層上水，相反，在一些河流盆地地區(qū)如楚瑪爾河斷陷盆地、通天河盆地，由于表層為粉砂層、卵礫石層，地形平緩，凍土上限深，凍結(jié)層上水較為發(fā)育[33]。以沱沱河為例，凍結(jié)層上水在其流域內(nèi)主要賦存于基巖類山區(qū)的活動層中，一般存在于每年5-9月，11月至次年4月全部凍結(jié)。
含水層厚度也受多年凍土層的控制，一般小于3 m。這類凍結(jié)層上水通過風化裂隙及構(gòu)造裂隙接受大氣降水的補給后，產(chǎn)生的徑流沿著凍結(jié)面進行流動，在山腳和山坡轉(zhuǎn)折處以線狀、片狀的方式泄于地表或順構(gòu)造裂隙帶直接補給凍結(jié)層下水[34-35]。在沱沱河北部及南部的高山與高海拔丘陵地區(qū)，由于降水較少、寒凍風化較弱、節(jié)理裂隙不發(fā)育等原因，導致補給水源不足，加之所處地區(qū)儲水空間有限，因此活動層內(nèi)土體的富水性較差[36]。
在沱沱河東部的極高山地區(qū)，由于降水充沛、寒凍風化強烈、節(jié)理裂隙發(fā)育，加之山體寬厚、巖屑坡廣泛分布，因此整個地區(qū)有利于大氣降水入滲補給，富水性較好。除凍結(jié)層上水的空間分布特征之外，還有對其時間變化特征的研究，通過研究發(fā)現(xiàn)，風火山地區(qū)典型高寒草甸坡面凍結(jié)層上水位的季節(jié)動態(tài)變化與活動層土壤溫度和水分變化類似，活動層土壤溫度控制了凍結(jié)層上水季節(jié)性動態(tài)變化格局。研究表明，距地表60cm深度以下土壤含水率和地下水補給來源決定了凍結(jié)層上水水位動態(tài)變化的位相和幅度[7]。
在青藏工程走廊內(nèi)，熱融湖塘作為富冰凍土區(qū)的一種常見地貌形態(tài)，其演化過程對于凍土退化、地表水和地下水動態(tài)平衡均具有重要影響，湖塘的存在及積水范圍的擴大成為了工程安全運營的潛在隱患，因此湖塘的演化特征在多年凍土區(qū)備受關注[37-38]，而凍結(jié)層上水在湖塘周邊廣泛發(fā)育，其排泄和補給過程與湖塘的變化密切相關。青藏工程走廊內(nèi)青藏公路兩側(cè)從西大灘到可可西里山之間分布著數(shù)百個熱融湖塘，且隨著多年凍土的退化，湖塘數(shù)量呈增加趨勢[39-40]。
凍結(jié)層上水的流入和流出影響著熱融湖塘的水文過程，為了評估凍結(jié)層上水的熱效應對于熱融湖塘的發(fā)展的影響，可以通過水位、水溫、凍結(jié)層上水流量的觀測，結(jié)合探地雷達、電阻率層析成像、鉆孔溫度監(jiān)測等方式對活動層厚度的時空變化以及湖塘周圍凍土分布進行研究。研究發(fā)現(xiàn)，廊內(nèi)典型湖塘周圍的凍結(jié)層上水每年融化時長持續(xù)5個月，湖岸邊溫度較低的凍結(jié)層上水將補給湖水，在這一過程中伴隨著與下伏凍土熱交換的產(chǎn)生，將引起活動層厚度和湖塘周圍凍土特征的定向差異[41]。
2.2黃河源區(qū)
黃河源區(qū)位于青藏高原東北部、青海省瑪多縣境內(nèi)，承擔著35%以上的黃河流域徑流量，是整個黃河流域重要的水源地，極易受到人為和自然因素的影響而發(fā)生凍土環(huán)境退化[42-43]。
黃河源區(qū)凍融和風化作用強烈，地表巖石發(fā)生破碎，產(chǎn)生大量的裂隙，這有利于凍結(jié)層上水的儲存。源區(qū)內(nèi)地下水水位近幾十年來呈現(xiàn)明顯下降的趨勢，甚至有些區(qū)域出現(xiàn)了凍結(jié)層上水消失的現(xiàn)象。多年凍土的退化使得冷生隔水層下移，導致了包氣帶土壤層的含水量減少，該區(qū)域出現(xiàn)了植被草場退化、植物多樣性減少、沼澤濕地萎縮、荒漠化加劇及黃河斷流等一系列生態(tài)環(huán)境問題[30,44]。整個源區(qū)內(nèi)凍結(jié)層上水的發(fā)育與凍土類型密切相關，主要發(fā)育在含冰量較高的凍土層上部的季節(jié)融化帶內(nèi)，埋藏深度受巖性、地層結(jié)構(gòu)、坡向、植被覆蓋率等諸多因素綜合作用控制，埋藏深度一般為0.5~2m[45]。
凍結(jié)層上水占水資源總量的88.3%，受凍結(jié)層上水供給的洼地和沼澤面積在逐年縮小。夏季和秋季入流量大，而春季和冬季由于活動層的凍結(jié)作用而入流量明顯降低[24]。黃河源區(qū)的凍結(jié)層上水主要分布在地形相對低洼的地表水和地下水匯集區(qū)，如巴顏喀拉山主峰-尕拉拉錯周圍的大型古冰川雪粒盆、皺瑪曲源頭的冰蝕洼地、以及沖積扇前緣及溝谷中[46]。
黃河源區(qū)的水文地質(zhì)結(jié)構(gòu)、地下水位的變化以及過度放牧和開采沙金等人為作用均會對凍結(jié)層上水的空間分布產(chǎn)生影響，這些作用不僅影響著凍土類型的退化，還會引發(fā)一系列生態(tài)問題[46]。源區(qū)內(nèi)凍結(jié)層上水的主要補給來源是大氣降水和地表水。在降水期內(nèi)，大氣降水和地表水給予地下水充足的水量補給，同時也伴隨著水熱交換，從而引起凍結(jié)層上界的消融和退化。凍結(jié)層上水主要通過河流、泉和蒸發(fā)排泄，有時通過融區(qū)補給凍結(jié)層下水[45]。從沼澤草甸排泄出的凍結(jié)層上水匯聚成溪流，成為源區(qū)主要河流的發(fā)源地，是河流的主要補給來源之一。
2.3祁連山
祁連山地區(qū)多年凍土廣泛發(fā)育，且溫度較高，是我國第二大內(nèi)陸河黑河的上游地區(qū)[47]。此地區(qū)凍地下水的補給來源主要是降水，在海拔3800m以上的多年凍土區(qū)，地下水以凍結(jié)層上水和凍結(jié)層下水形式存在;在海拔3800m以下的非多年凍土區(qū)，地下水主要以潛水形式存在[48]。基巖山區(qū)凍結(jié)層上水儲存于風化裂隙或孔隙中，含水層厚度隨巖性的差異而不同，厚度一般為1~2m，山區(qū)降水充沛，補給條件良好，凍結(jié)層上水較發(fā)育且水量豐富。
各盆地內(nèi)凍土平原凍結(jié)層上水則儲存于第四系冰磧-冰水相含泥砂礫石層中，含水層厚度0.5~2m[49]。凍結(jié)層上水與地表水、其他類型的地下水之間具有水力聯(lián)系。凍結(jié)層上水成為了冰川和積雪融水從高山往平原地區(qū)遷移的通道，作為蓄水層，凍結(jié)層上水對于水分的疏導能力產(chǎn)生的季節(jié)性變化主要體現(xiàn)在兩方面：在暖季沿著地勢從高高低流向產(chǎn)生溪流的溝谷地帶;在冷季流向因河冰和季節(jié)性凍結(jié)造成的臨時承壓含水層，這也說明冷季平原地區(qū)儲存的地下水和河槽中的基流進行了緩慢釋放[50]。
為了研究凍土和地下水的交互影響，通過對祁連山葫蘆溝流域三維水文地質(zhì)模型的計算結(jié)果表明，此流域地下水在基流的形成中占43%，地表溫度升高會使多年凍土的分布范圍減小，增加地下水向河流的排泄。整個流域內(nèi)凍結(jié)層上水距有明顯的流動性特征，在滲流過程中有少部分向下穿透了多年凍土層。凍結(jié)層下水水頭較高，而凍結(jié)層上水水頭較低，向上的水力梯度驅(qū)動著溫度更高的地下水流向地表，促進暗渠和融區(qū)的形成，增加地下水與地表水的交互作用[51]。
3高緯度多年凍土區(qū)
3.1大興安嶺
大興安嶺位于歐亞大陸多年凍土區(qū)的南部邊緣，海拔1100~1500m，凍結(jié)層上水在山間溝谷和河谷谷地最為發(fā)育[52]。區(qū)域內(nèi)植被茂密、森林眾多、降雨量大，水系有黑龍江和嫩江及其支流，如額木爾河、盤古河、呼瑪河、訥謨爾河、烏裕爾河等，濕地沼澤較為發(fā)育，地下水的賦存特征受地質(zhì)構(gòu)造、巖性、風化作用和地形地貌等因素的影響。凍結(jié)層上水主要分布于大興安嶺北部較大的河谷平原區(qū)，含水層厚度、富水性隨季節(jié)融化深度的增加而增大，此地區(qū)凍結(jié)層上水主要接受大氣降水入滲補給，豐水期接受地表水補給，枯水期排泄于地表水[53-54]。此區(qū)域凍結(jié)層上含水層的分布、厚度、演變程度既與地形地貌有關，又受植被覆蓋度的影響，還與地表粗糙度等地形因子有關。
植物根系及堆積的枯枝敗葉對降水起到良好的輸導入滲和涵養(yǎng)功能，完整基巖就構(gòu)成了穩(wěn)定的隔水底板，地下水在基巖風化帶里形成了凍結(jié)層上水。凍結(jié)層上水分布空間狹窄，富水條件普遍不佳，再加上接受降水補給時間短，水量貧乏，動態(tài)波動特征明顯[55]�；鶐r風化帶凍結(jié)層上水按照流量大小分為水量中等和水量貧乏兩類：中等水量的凍結(jié)層上水分布于伊勒呼里山分水嶺北部大片中低山區(qū)，含水層巖性較為復雜，主要有侏羅系火山角礫巖;水量貧乏的凍結(jié)層上水主要分布于伊勒呼里山分水嶺南部中山區(qū)以及分水嶺北部零星山區(qū)，含水層巖性主要為侏羅系火山巖[56]。
3.2俄羅斯
俄羅斯多年凍土面積占國土面積的一半以上，是北半球多年凍土區(qū)的重要組成部分[57]。俄羅斯典型地區(qū)的凍結(jié)層上水動態(tài)特征與氣溫變化密切相關，且含水層具有凍結(jié)層上水和凍結(jié)層間水的特點。在東西伯利亞勒拿河盆地，凍結(jié)層上水分布廣泛，占多年凍結(jié)層面積的40%~60%，凍結(jié)層上水和凍結(jié)層間水混合體的滯留時間為1~55a，多年凍土退化將增大地下水的補給[58]。
在季曼-伯朝拉地區(qū)，地下水水位動態(tài)變化與氣溫之間聯(lián)系密切。該地區(qū)凍結(jié)層上水水位一年最小值出現(xiàn)在冬季末期，一年之中，近地面大氣平均晝夜溫度超過0℃后，包氣帶水分交換過程和凍結(jié)層上水水位開始出現(xiàn)季節(jié)性變化[59]。俄羅斯科學院西伯利亞分院凍土研究所和俄羅斯地質(zhì)部的東北水文地質(zhì)分部共同研究了不同地貌和巖性條件下凍結(jié)層上水的季節(jié)性變化，總結(jié)了三種典型地區(qū)的水位變化特征。
當凍結(jié)層上水水位距地表較淺時，依靠積雪融化和大氣降水進行補給，水位在衰減時期也能上升;當水位距離地表較深時，基本不受大氣降水的影響，因此水位在夏季呈總體下降的趨勢;當水位與季節(jié)性融化層水資源枯竭無關，與含水層自上而下和自下而上的凍結(jié)有關[60]。除此之外，一些特殊地區(qū)如河床下融區(qū)產(chǎn)生的凍結(jié)層上水不受季節(jié)性凍融過程的直接影響，影響程度取決于河谷的地貌、凍土特征、水文地質(zhì)條件等[61]。高緯度多年凍土區(qū)人類活動較多，除了道路等線性工程之外，還分布有市政建筑等眾多基礎設施，城市的發(fā)展和擴張也對凍結(jié)層上水產(chǎn)生影響。
通過已有文獻可知，俄羅斯高緯度地區(qū)人類活動對凍結(jié)層上水補給的影響主要包括耕作區(qū)部分灌溉入滲和來自城市管道、蓄水系統(tǒng)水的滲漏兩方面，通過對俄羅斯北部城市雅庫茨克市區(qū)開展的綜合性水文監(jiān)測工作的研究表明，雖然人類活動因素產(chǎn)生的影響呈局部性特征，但卻使凍結(jié)層上水的水文情勢發(fā)生了本質(zhì)性變化[60]。
3.3北美
北美地區(qū)的多年凍土主要分布于阿拉斯加和加拿大北部，對于凍結(jié)層上水的研究主要集中在幾條主要河流的流域盆地如育空河、馬更些河以及環(huán)北極地區(qū)的康沃利斯島等區(qū)域[62]。
育空河全長約3300km，穿越了加拿大及阿拉斯加地區(qū)，流域內(nèi)多年凍土的分布和厚度差異性較大。隨著多年凍土的融化，近30a以來地下水占河流總流量的比例正以每年0.7%~0.9%的速度持續(xù)增大[63]。馬更些河是加拿大第一長河，馬更些河盆地占整個加拿大國土面積的近四分之一，盆地中有75%的區(qū)域分布有連續(xù)或不連續(xù)多年凍土[64]，在春季融化過程中，凍結(jié)層上水占整個排泄量的90%，這表明在此盆地中，凍結(jié)層上水是極地水文的重要組成部分[65]。
對于育空河盆地地區(qū)凍結(jié)層上水的研究主要集中在活動層凍融過程、氣候-湖泊-地下水之間關系、側(cè)向流動特征三個方面，研究方法以實地觀測和模擬計算為主，將二者結(jié)合，能夠透過現(xiàn)象揭示機理，為查明凍結(jié)層上水的變化特征奠定基礎。育空河盆地多年凍土區(qū)地下水向河流的排泄量逐年增大，地下水總量將有增大的趨勢，這首先表現(xiàn)為地下水在低洼地帶的排泄范圍增大，除此之外還發(fā)現(xiàn)，此地區(qū)凍結(jié)層上水對于基流的貢獻率有所降低[66]。
在氣候變化背景下，湖泊與地下水之間在融區(qū)形成過程中有著密切關系，氣溫變暖會加劇凍土融化、增加凍結(jié)層上水流量。在寒區(qū)，除了氣候是凍土環(huán)境變化的主導因素外，湖面面積和地下水流動在多年凍土融化和融區(qū)形成過程中有著重要影響。湖邊緣季節(jié)性的凍融循環(huán)影響著凍結(jié)層上水的運移和凍結(jié)過程，在湖面產(chǎn)流區(qū)域，多年凍土的分布影響著湖中水體在凍結(jié)層上下之間的交換，與此同時，凍土分布特征也受湖水及凍結(jié)層的影響[67]。
凍結(jié)層上水在坡面地區(qū)將產(chǎn)生側(cè)向流動，與垂向流動一并形成了其特殊的水文地質(zhì)現(xiàn)象，狼溪(WolfCreek)盆地坡面垂向的滲透和蒸發(fā)過程在南坡較為明顯，而北坡的下伏多年凍土阻止了深層滲漏，使得凍結(jié)層上水得以存儲[68]。加拿大的多年凍土主要分布在北部的極地地區(qū)，對于凍結(jié)層上水的研究主要集中在坡面、河流流域以及極地島嶼等典型區(qū)域。在極地地區(qū)，產(chǎn)生坡面地表徑流的原因在于凍結(jié)層上水水位的上升和活動層土壤的飽和，坡面徑流取決于凍結(jié)層上水水位，并與地表形態(tài)有關。
從徑流源頭來看，融雪所占比例最大，降雨次之，排泄方式依次為蒸發(fā)、側(cè)向徑流。側(cè)向徑流包括地表流和地下流，地下流則跟導水系數(shù)、活動層厚度、水位等因素有關[69]。在加拿大北部的礫石地貌中，有濕地、積水和徑流，所有的這些地貌都與凍結(jié)層上水的局地滲流有關。
礫石具有很高的單位產(chǎn)水量和導水系數(shù)，能夠促進水分流動，但是，淺層凍結(jié)鋒面影響著深層的滲透，因此，當有融雪和降雨的充足補給時，在地表與水位相交的地方便有滲流，且具有季節(jié)性變化特征。通過水量平衡計算可知，上坡礫石區(qū)域的滲流是濕地和積水的主要來源，連續(xù)性的滲流改變了植被特征和土壤含水率[70]。此外，凍結(jié)層上水對河流流量有一定程度的貢獻，在河道下方的活動層厚度通常更大，甚至出現(xiàn)融區(qū)，原因在于凍結(jié)層上水的攜熱效應[71]。
4凍結(jié)層上水對工程穩(wěn)定性的影響
由上述可知，多年凍土對于凍結(jié)層上水的影響特征主要有以下四個方面：①減緩水循環(huán)過程;②阻隔水分運動的垂向通道，影響地表水和地下水的轉(zhuǎn)換;③使含水層的變化具有季節(jié)性特點和動態(tài)規(guī)律;④凍脹作用改變了土壤結(jié)構(gòu)并增強了入滲能力[14]。隨著多年凍土區(qū)人類活動的加劇和各類工程的修筑，凍結(jié)層上水的存在及其滲流攜熱效應嚴重影響著各類工程基礎的熱收支特征，直接對工程穩(wěn)定性產(chǎn)生威脅，凍結(jié)層上水滲流、凍土退化及工程穩(wěn)定性之間的相互影響，但三者的定量關系有待研究。從已有文獻來看，凍結(jié)層上水對于工程穩(wěn)定性的影響研究主要集中在公路、鐵路以及機場。對于多年凍土區(qū)的道路工程而言，路基修筑后凍土上限上升，將有一部分凍結(jié)層上水在上限抬升處進入并橫穿路基。
隨著運營時間的增加，活動層和路基中凍結(jié)層上水(2~4℃)的滲流將造成多年凍土溫度升高，造成多年凍土退化，同時還會加速融區(qū)的出現(xiàn)，加劇路基不均勻沉降，導致路基病害的出現(xiàn)[76]。通過室內(nèi)模型試驗，結(jié)合現(xiàn)場監(jiān)測資料，能夠從凍結(jié)層上水夏季滲流導致路基下伏多年凍土融化的角度，進一步證明路面橫向裂縫、縱向裂縫以及路基不均勻沉降的成因[77]。針對凍結(jié)層上水對公路路基的影響，可采用防水材料，阻隔水分遷移和滲透通道，對路基下部過水進行隔水治理可有效降低路基沉降，同時降低因為活動層集水而造成的凍脹[78]。
5理論研究
對于凍結(jié)層上水的理論研究主要可以歸納為三個方面：一是水文模型的建立與改進;二是水量平衡計算，三是滲流攜熱效應。含水層土體的導水系數(shù)、液態(tài)水儲存能力、凍結(jié)層上水水位隨著土體的凍融與溫度密切相關，而且活動層的雙向凍結(jié)控制和改變了凍結(jié)層上水的滲流過程，向下凍結(jié)控制著地表徑流的產(chǎn)生和凍結(jié)層上水的補充，向上凍結(jié)影響著地下水尤其是凍結(jié)層上水的排出[87]。由此可見，凍結(jié)層上水的滲流過程和水位變化特征受局地因素的控制，極為復雜，因此截至目前，仍然沒有有效的數(shù)學模型和理論來預測多年凍土區(qū)地下水的變化特征[88]。
5.1水文模型
國際上第一代水文模型很好地描繪了流域產(chǎn)流的基本規(guī)律，如斯坦福模型、水箱模型、新安江模型，但這類模型忽略了地形、植被等的空間變化對流域產(chǎn)流過程的影響。隨后發(fā)展出了一批分布式水文模型，如TOPMODEL、SHE、VIC、DHSVM、GBHM等[89]。近十多年來，由于氣象學、生態(tài)學、冰凍圈等學科的興起和發(fā)展，分布式水文模型與多學科的交叉研究越來越多，這使得建模難度增加，分布式水文模型開始描述不同圈層之間的關系，即“多圈層水文模型”[90-91]。目前中國寒區(qū)水文模型的主要限制因素是缺乏全面、同步的觀測數(shù)據(jù)，可通觀測、試驗、模擬的方法，獲取寒區(qū)水文過程的基本規(guī)律、經(jīng)驗公式和關鍵參數(shù)[92]�；谖锢韺W基礎的水文數(shù)值模型能夠定性和定量地分析含水層凍融過程中的關鍵影響因素，且能夠反映出水文要素的變化。
在已有水文模型的基礎上，可以采用熱傳導方程模擬凍土上限在水文模型中的演變過程[93]，結(jié)合能量平衡計算地表熱狀況的空間分布特征，以此為基礎可以得到土壤剖面的融化深度[94]，也可以通過水分遷移模型和熱傳導模型的耦合，可以研究多年凍土形成和演化過程中的水力特性，同時研究地下水和地表水的季節(jié)性變化特征[95]，還可以采用達西方程、流體連續(xù)性方程、溶質(zhì)連續(xù)性方程、能量遷移方程來計算冰川作用下地下水流動[96]。此外，也可以通過在已有的簡易地下水模型(SIMGM)、陸面模式(GLM)中增加未凍水含量、相變、土層性質(zhì)等參數(shù)等方式進行改進[97]，建立面向?qū)ο蟮哪Ｐ拖到y(tǒng)，來模擬和預測中小尺度的區(qū)域中水循環(huán)特征[98]及對氣候的反饋[99-100]。
6小結(jié)與展望
凍結(jié)層上水的研究正從定性的特征分析向定量的機理研究發(fā)展，已有的觀測和模擬研究已經(jīng)揭示了凍結(jié)層上水的影響因素、賦存特征、補給排泄以及動態(tài)變化。凍結(jié)層上水的存在使多年凍土區(qū)的水文地質(zhì)條件變得更為復雜，凍結(jié)層上水與多年凍土層的相互影響、相互制約的交互作用將會對整個凍土區(qū)的生態(tài)環(huán)境產(chǎn)生深遠影響。在今后有必要在以下幾個方面進行深入研究：
(1)對局地因素條件下凍結(jié)層上水的水文要素時空變化特征進行研究。如在坡面地區(qū)，需要繼續(xù)深入研究凍結(jié)層上水對于坡面徑流的調(diào)控機理，可針對典型地區(qū)不同的水文地理和水文生態(tài)單元開展比較研究，用以獲得春季融化和秋季凍結(jié)期間凍結(jié)層上水對于坡面徑流的補給占比以及年度變化特征。
(2)對多年凍土區(qū)淺表層地下水與地表水的相互關系及轉(zhuǎn)化特征進行研究。在地勢緩和且凍結(jié)層上水發(fā)育地區(qū)，往往伴隨著一定的地表積水，應該對地表積水入滲至活動層以及凍結(jié)層上水在凍土退化過程中遷移至地表的機理、變化特征進行深入研究，揭示地表積水與凍結(jié)層上水之間的轉(zhuǎn)換規(guī)律。
(3)對多年凍土區(qū)路基穩(wěn)定性受凍結(jié)層上水的影響機理進行研究。多年凍土區(qū)路基/地基下部凍結(jié)層上水的存在成為其穩(wěn)定性的潛在隱患，有必要對路基沉降變化特征與凍結(jié)層上水之間的關系進行研究，并對凍結(jié)層上水的滲流攜熱效應對路基變形的影響程度進行定量研究，以此來修正有凍結(jié)層上水存在條件下的路基變形計算模型。
(4)對氣候變暖背景下凍結(jié)層上水與多年凍土交互作用及變化特征進行預測。需要對多年凍土與凍結(jié)層上水的相互作用進行研究，揭示活動層的凍融過程對多年凍土退化的影響，以及多年凍土對凍結(jié)層上水流動特征的限制，并結(jié)合氣候變化模式預測二者的動態(tài)發(fā)展規(guī)律，獲得更具有普適意義的凍結(jié)層上水動態(tài)規(guī)律和驅(qū)動機制。
(5)對凍結(jié)層上水水位、流速、構(gòu)成等基本參數(shù)在測試過程中的新方法和新手段進行探索。目前對于凍結(jié)層上水的監(jiān)測主要集中在水位和水溫方面，受客觀條件制約，凍結(jié)層上水的流量和流速很難被精確測量，針對這一問題，應采用新的監(jiān)測方法和傳感器，提高測量精度，從整體上掌握凍結(jié)層上水的存在范圍和流動特征，為理論研究提供更為精確的參數(shù)。
參考文獻(References):
[1]DingYongjian,ZhangShiqiang,ChenRensheng,etal.Introductiontohydrologyincoldregions[M].Beijing:SciencePress,2017.[丁永建,張世強,陳仁升,等.寒區(qū)水文導論[M].北京:科學出版社,2017.]
[2]MetcalfeRA,ButtleJM.SoilpartitioningandsurfacestorecontrolsonspringrunofffromaborealforestpeatlandbasininnorthcentralManitoba,Canada[J].HydrologicalProcesses,2001,15(12):2305-2324.
[3]ChengGuodong,JinHuijun.GroundwaterinthepermafrostregionsontheQinghai-TibetPlateauanditchanges[J].HydrologyandEngineeringGeology,2013,40(1):1-11.[程國棟,金會軍.青藏高原多年凍土區(qū)地下水及其變化[J].水文地質(zhì)工程地質(zhì),2013,40(1):1-11.]
[4]LiuGuangsheng,WangGenxu,SunXiangyang,etal.VariationcharacteristicsofstableisotopesinprecipitationandriverwaterinFenghuoshanpermafrostwatershed[J].AdvancesinWaterScience,2012,23(5):621-627.[劉光生,王根緒,孫向陽,等.多年凍土區(qū)風火山流域降水河水穩(wěn)定同位素特征分析[J].水科學進展,2012,23(5):621-627.]
[5]WangGX,HuHC,LiTB.Theinfluenceoffreeze-thawcyclesofactivesoillayeronsurfacerunoffinapermafrostwatershed[J].JournalofHydrology,2009,375(3/4):438-449.
作者：柴明堂，馬巍，穆彥虎
轉(zhuǎn)載請注明來自發(fā)表學術論文網(wǎng)：http:///nylw/28701.html
