劉景明, 丁建麗, 包青嶺, 張子鵬, 姜磊鵬, 曲 藝
(1.新疆大學(xué)地理與遙感科學(xué)學(xué)院,新疆 烏魯木齊 830046;2.新疆綠洲生態(tài)重點(diǎn)實(shí)驗(yàn)室,新疆 烏魯木齊830046;3.智慧城市與環(huán)境建模自治區(qū)普通高校重點(diǎn)實(shí)驗(yàn)室,新疆 烏魯木齊 830046)
大規(guī)模的地下水開(kāi)采造成地面沉降、干旱期泉水和河流基流減少、咸水入侵、水質(zhì)退化、甚至全球海平面上升等問(wèn)題,并對(duì)世界各地的生態(tài)系統(tǒng)造成損害[1]。特別是內(nèi)陸干旱地區(qū),地下水對(duì)于調(diào)節(jié)水文循環(huán)和維系生態(tài)系統(tǒng)健康至關(guān)重要[2]。其中地下水的補(bǔ)給、交換機(jī)制和組成成分是水文學(xué)和水文地質(zhì)學(xué)的基本問(wèn)題之一。在傳統(tǒng)的常規(guī)手段較難獲取這些信息的情況下,同位素水文學(xué)技術(shù)已被證明是解決許多關(guān)鍵水文問(wèn)題和過(guò)程的有效工具,同時(shí)開(kāi)展水化學(xué)多方法聯(lián)合示蹤并相互驗(yàn)證可以提高評(píng)價(jià)精度[3-4]。
利用水化學(xué)和同位素方法對(duì)干旱區(qū)流域地下水補(bǔ)給和水循環(huán)的研究已有開(kāi)展。通過(guò)對(duì)比地下水與潛在水源的同位素組成可以示蹤干旱區(qū)地下水補(bǔ)給來(lái)源[5],Yapiyev等[6]通過(guò)水的氫氧同位素組成確定中亞地區(qū)的蒸發(fā)損失和地下水輸入;Joshi等[7-8]利用水同位素追蹤印度恒河流域地下水補(bǔ)給源;Jesiya 等[9-10]根據(jù)地下水的穩(wěn)定同位素時(shí)間序列數(shù)據(jù)研究地下水補(bǔ)給機(jī)制。同時(shí),結(jié)合水化學(xué)和氫氧穩(wěn)定同位素指標(biāo)成為定量評(píng)價(jià)地下水與河水交換的有效手段[11],文廣超等[12-15]學(xué)者分別研究了鄱陽(yáng)湖、孔雀河、伊犁河谷和巴音河流域地下水和地表水轉(zhuǎn)換關(guān)系。在艾比湖流域地下水的相關(guān)研究中,郝帥等[16-18]分析了流域地下水氫氧穩(wěn)定同位素的時(shí)空變化特征,揭示了水體穩(wěn)定同位素與環(huán)境要素關(guān)系;朱世丹等[19-20]結(jié)合水化學(xué)和氫氧穩(wěn)定同位素技術(shù),分析艾比湖主要入湖河流氫氧同位素及水化學(xué)的組成特征,但地下水涉及較少。因此,本次工作在前人研究的基礎(chǔ)上,綜合利用水化學(xué)和氫氧穩(wěn)定同位素開(kāi)展艾比湖流域地下水研究。
艾比湖流域作為“一帶一路”沿線尤為重要的核心區(qū)域,其水土安全問(wèn)題關(guān)系到國(guó)家戰(zhàn)略實(shí)施[21]。艾比湖流域具有典型的干旱區(qū)生態(tài)環(huán)境特征,近年來(lái)由干旱引發(fā)的湖泊面積萎縮,使湖濱周邊干涸湖底大面積裸露且嚴(yán)重鹽堿化[22]。而地下水資源對(duì)流域生態(tài)環(huán)境和區(qū)域水資源管理至關(guān)重要,提供地下水的穩(wěn)定同位素和水化學(xué)組成證據(jù),可以為跨流域調(diào)水和流域治理提供理論指導(dǎo)。研究的具體目標(biāo)是(1)艾比湖流域地下水同位素水化學(xué)特征;(2)確定補(bǔ)給源和補(bǔ)給特征;(3)了解艾比湖流域含水層系統(tǒng)的流動(dòng)系統(tǒng)。
艾比湖流域位于新疆維吾爾自治區(qū)博爾塔拉蒙古自治州精河縣北部,湖泊呈西北—東南走向,地理位置位于44°34′~45°08′N、82°35′~83°16′E 之間(圖1)。流域地勢(shì)較低,是一個(gè)半閉合的洼地,東邊開(kāi)口,其他三面被高山環(huán)繞,湖盆是盆地內(nèi)地表水和地下水的匯集中心。根據(jù)博爾塔拉政府網(wǎng)(http://www.xjboz.gov.cn/zjbz.htm)公布的2018 年博州水資源公報(bào)可知,全流域地表水開(kāi)采量從1980年的1.35×108m3增加到2018 年的9.98×108m3,地下水呈嚴(yán)重超采狀態(tài)。根據(jù)全州19 眼地下水觀測(cè)井統(tǒng)計(jì),全州地下水埋深平均下降了0.12 m(較上年增加0.26 m)。
圖1 采樣點(diǎn)區(qū)域示意圖Fig.1 Location of sampling sites in the study area
本研究于2020 年和2021 年7 月對(duì)流域地表水和地下水進(jìn)行采集,采樣點(diǎn)主要取自博河、精河斷面和艾比湖周邊管護(hù)站,共采集地下水樣本30個(gè)和地表水樣本15 個(gè)。所有水樣都收集在150 mL 的聚乙烯瓶中,并在現(xiàn)場(chǎng)用樣品水沖洗瓶子2 次。為避免樣品的擴(kuò)散和蒸發(fā)損失,密封并冷凍樣品瓶,野外任務(wù)結(jié)束后立即帶到實(shí)驗(yàn)室進(jìn)行同位素分析。采樣期間,使用手持全球定位系統(tǒng)接收器記錄了樣本坐標(biāo)和海拔。
水化學(xué)分析在室內(nèi)25 ℃進(jìn)行,電導(dǎo)率(Electrical conductance,EC)、溶解性總固體(Total dissolved solids,TDS)和鹽度(Salinity,SAL)采用電導(dǎo)率儀(Multi 3420 Set B,WTW GmbH,Germany)測(cè)定,酸堿度(pH)采用離子選擇性電極法測(cè)定,取3 次平均值作為最終的結(jié)果。水樣的氫氧同位素測(cè)試使用全自動(dòng)真空冷凝抽提系統(tǒng)(LI-2100,Automatic water extraction system)和液態(tài)水穩(wěn)定同位素分析(TIWA-45-EP,Los Gatos Research)完成樣本的提取和測(cè)試。每個(gè)樣本進(jìn)行6次分析,考慮記憶效應(yīng),采用后4 次的分析結(jié)果。測(cè)得的氫氧穩(wěn)定同位素含量為維也納標(biāo)準(zhǔn)平均大洋水(V-SMOW)的千分偏差,計(jì)算公式為:
式中:Rsample為樣本;Rstandard為標(biāo)樣。測(cè)試結(jié)果以相對(duì)V-SMOW 千分偏差表示(δ,‰),測(cè)試精度:δ2O為±0.33‰,δ18H為±0.5‰。
2.2.1 線性回歸方程的擬合最小二乘法(又稱最小平方法)是一種數(shù)學(xué)優(yōu)化技術(shù)。它通過(guò)最小化誤差(真實(shí)目標(biāo)對(duì)象與擬合目標(biāo)對(duì)象的差)的平方和尋找數(shù)據(jù)的最佳函數(shù)匹配[23]。最小二乘法可用于曲線擬合,利用最小二乘法建立不同區(qū)域地下水δ2H與δ18O之間的線性回歸方程,設(shè)擬合直線的公式為:
式中:y為δ2H 值;x為δ18O 值;b為方程截距;擬合直線的斜率為;計(jì)算出斜率后,根據(jù)和已經(jīng)確定的斜率k,利用待定系數(shù)法求出截距b。
2.2.2 地下水補(bǔ)給率的計(jì)算根據(jù)線性回歸方程的擬合結(jié)果,利用同位素守恒規(guī)律,假設(shè)地下水樣本的同位素組成是通過(guò)混合2 個(gè)端元組分得出的,即降水和地表水,使用雙端元混合模型估算流域地下水 的 補(bǔ) 給 率[7,24]。因 為δ18O 水 轉(zhuǎn) 化 率 精 度 高 于δ2H[25],因此公式可表示為:
式中:Qp、Qg、Qs分別為降水量、地下水量、河水量;δ18Op、δ18Og、δ18Os分別為降水量、地下水量、河水量的δ18O值。
2.2.3 瑞利分餾方程瑞利分餾模型本源于描述相反應(yīng)過(guò)程中混合在一起的兩種不同溶液的蒸餾過(guò)程,后被用于描述開(kāi)放系統(tǒng)中穩(wěn)定同位素的非平衡分餾過(guò)程。根據(jù)瑞利分餾公式,δ18O存在如下關(guān)系[26]:
式中:δ18O 為蒸發(fā)后樣品的同位素值;δ18O0為初始同位素值;ε為富集系數(shù);f為剩余水比例。
3.1.1 不同區(qū)域地下水氫氧同位素特征地下水采樣點(diǎn)主要分布在博河、精河和艾比湖周邊區(qū)域,由于精河上游取樣難度大,只取得中下游地下水。博河、精河區(qū)域井深15~100 m,艾比湖周邊管護(hù)站井深30~260 m。由表1均值可知,整個(gè)流域δ2H和δ18O最大值位于博河和精河中下游區(qū)域,艾比湖周邊區(qū)域次之,博河上游區(qū)域最小;氘盈余(d-excess)最大值位于艾比湖周邊區(qū)域,博河上游區(qū)域次之,博河和精河中下游區(qū)域最小。地下含水層分布在地表以下,因此地下水氫氧穩(wěn)定同位素的區(qū)域異質(zhì)性表明流域地下水受到補(bǔ)給源和水巖作用的影響。
表1 不同區(qū)域氫氧同位素值Tab.1 Isotopic values of hydroxide in different regions /‰
3.1.2 地下水氘盈余空間特征氘盈余(d-excess)最早由Dansgaard 定義并量化[27],d=δ2H-8δ18O,反應(yīng)了包括空氣濕度、海洋表面溫度和風(fēng)速的變化,也可反映該地區(qū)水-巖氧同位素交換程度[28]。整個(gè)流域d-excess 范 圍 從6.47‰ 到17.45‰,平 均 值 為10.67‰。為進(jìn)行流域d-excess 參數(shù)空間特征分析,將全球大氣降水線及其d=5、10、15 時(shí)的特征線繪出。從圖2中可以看出,博河、精河中下游區(qū)域主要位于d=5 和d=10 特征線之間,而艾比湖周邊地下水主要位于d=10 和d=15 特征線之間,表明流域不同區(qū)域地下水存在不同的循環(huán)過(guò)程。
圖2 艾比湖流域地下水d-excess特征Fig.2 Groundwater d-excess characteristics in Ebinur Lake Basin
圖3a顯示出整個(gè)流域內(nèi)地下水的d-excess在空間上的變化。水平分布上,博河、精河流域地下水采樣點(diǎn)分布在距河流0~8 km(圖3b)。不同距離下d-excess 表現(xiàn)出聚集特征,一定程度表明其地下水橫向連通性較好。精河流域J1 和J4 與艾比湖距離不同導(dǎo)致的同位素差異特征,與精河地表水同位素特征一致[20]。艾比湖周邊區(qū)域地下水位高(平均深度95 m),距河湖遠(yuǎn),且受外界環(huán)境因素影響較小,氫氧同位素整體顯現(xiàn)出貧化特征。其中ABH6距艾比湖較近(2.97 km),而d-excess相對(duì)較低,可能受到巖石含氧化合物的化學(xué)組分和含水層的溫度影響。距湖3~6 km內(nèi)d-excess差異性可能反映了區(qū)域含水層系統(tǒng)的非均質(zhì)和各向異性導(dǎo)致的有限橫向連通性[7]。垂直分布上(圖3c),除ABH8 和ABH12 極端樣本以外,回歸分析表明擬合關(guān)系不存在顯著性意義,相關(guān)系數(shù)極小,地下水d-excess 沒(méi)有表現(xiàn)出深度效應(yīng),未來(lái)需要更多的同位素指標(biāo)去分析驗(yàn)證。
圖3 艾比湖流域地下水d-excess隨河湖距離和深度的變化Fig.3 Changes of groundwater d-excess with distance of rivers and the lake and depth in Ebinur Lake Basin
3.2.1 地下水δ2H與δ18O相互關(guān)系圖4a為艾比湖流域地下水δ2H 與δ18O 相互關(guān)系,其δ2H 與δ18O 有顯著線性關(guān)系(R2=0.83,P<0.001),關(guān)系式為δ2H=4.79δ18O-24.31。艾比湖地區(qū)處于西風(fēng)帶,受季風(fēng)影響很弱[20],這使區(qū)域地下水同位素值偏離降水線和地表水線。根據(jù)Wang 等[29]在天山北坡、Hao 等[18]在艾比湖流域的降水穩(wěn)定同位素統(tǒng)計(jì)數(shù)據(jù)可知,天山北坡的大氣降水線為δ2H=7.51δ18O+1.95,流域大氣降水線為δ2H=6.69δ18O-6.53(圖4a)。干旱區(qū)氣候干燥,空氣濕度小,降水的云下蒸發(fā)強(qiáng)烈,再加上局地循環(huán)水汽的混合導(dǎo)致降水中同位素的動(dòng)力分餾現(xiàn)象,地下水δ2H 與δ18O 關(guān)系線整體位于天山北坡大氣降水線和流域大氣降水線的上方,斜率和截距均小于二者,且具有負(fù)截距,表明地下水接受降水補(bǔ)給的同時(shí)也經(jīng)歷了一定程度的蒸發(fā)富集,這與干旱區(qū)石羊河流域的研究一致[30]。其中博河地下水基本位于流域大氣降水線上,表明地下水主要來(lái)源于降水,精河地下水位于流域大氣降水線上方,表明該區(qū)域降水補(bǔ)給期間的蒸發(fā)分餾作用與博河不同。博河和精河地下水位于流域地表水線附近(圖4b),表明受到地表水補(bǔ)給的影響,博河和精河地下水主要補(bǔ)給源是地表水和降水;艾比湖周邊地下水位于流域大氣降水線上方,δ18O 較低,表現(xiàn)出貧化特征,且部分地下水處在流域地表水線的延伸線上,表明地下水主要來(lái)源于冰雪融水和降水,部分地下水受地表水補(bǔ)給。結(jié)合Tan 等[31]學(xué)者的研究發(fā)現(xiàn),蒸發(fā)對(duì)同位素的主要影響集中在上層土壤層,對(duì)地下水的影響較弱。因此,流域地下水的同位素變化主要受補(bǔ)給源的影響。
圖4 艾比湖流域地下水δ2H與δ18O相關(guān)關(guān)系Fig.4 Correlation between groundwater δ2H and δ18O in Ebinur Lake Basin
3.2.2 地下水影響因素分析基于地表水交換區(qū)δ18O 的平均值及式4 對(duì)地下水補(bǔ)給率進(jìn)行計(jì)算(博河、精河中下游地下水和地表水平均值為-10.31‰、-11.84‰,河湖交匯區(qū)地下水和地表水平均值為-11.39‰、-11.84‰),參考艾比湖流域降水δ18O 平均值-4.38‰[17],得出博河和精河中下游地表水對(duì)地下水的補(bǔ)給率為79.7%,降水補(bǔ)給率為20.3%。河湖交匯區(qū)地表水對(duì)地下水的補(bǔ)給率為93.9%,降水補(bǔ)給率為6.1%。該結(jié)果相較郝帥等[16]計(jì)算的中下游河水和地下水之間交換比例偏高,主要是季節(jié)灌溉回流補(bǔ)給作用大,也可能是參考降水δ18O 與實(shí)際采樣時(shí)間差異導(dǎo)致的。從地下水TDS與δ18O 關(guān)系圖可知,地下水同位素與TDS 的相關(guān)性較低,表明地表水和降水的補(bǔ)給對(duì)TDS 含量的影響較?。▓D5)。已知地下水最富集樣本δ18O 為-9.81‰,最貧化樣本δ18O 為-11.01‰,參考艾比湖流域富集系數(shù)為-83.33‰[18],根據(jù)瑞利分餾公式,計(jì)算得到剩余水比例為98%。再根據(jù)質(zhì)量平衡原理,假設(shè)蒸發(fā)是控制TDS 的主要因素,地下水TDS 含量應(yīng)為(C0對(duì)應(yīng)地下水最小TDS 含量)。這與地下水實(shí)際TDS 含量變化相差較大(圖5),因此必然存在其他因素影響地下水的TDS 變化。水中TDS 主要由水中的HCO3-、Ca2+的濃度決定,受采樣區(qū)范圍內(nèi)的巖石風(fēng)化作用的控制[32],因此需要結(jié)合水位和水質(zhì)數(shù)據(jù)進(jìn)一步評(píng)估。
圖5 艾比湖流域地下水TDS與δ18O關(guān)系Fig.5 Relationship between groundwater TDS and δ18O in Ebinur Lake Basin
3.2.3 地下水氘盈余沿程變化圖6a~b 顯示的是博河和精河流域d-excess 從上游到下游的空間變化趨勢(shì)。博河上游地下水d-excess 較高,中下游區(qū)域較為穩(wěn)定,下游區(qū)域存在波動(dòng)。精河中下游區(qū)域地下水d-excess 沿程逐漸降低,但整體比博河中下游高。博河上游區(qū)域B1樣點(diǎn)d-excess較高,且δ18O與郝帥等[16]在博河上游山谷冬季雪水氧同位素相近(均值-18.63‰),反映上游靠近山區(qū),主要補(bǔ)給源是冰川融水,氣溫低,蒸發(fā)弱。中下游區(qū)域d-excess 總體較低,表明補(bǔ)給地下水存在一定蒸發(fā)富集。根據(jù)地下水補(bǔ)給率的計(jì)算,可知存在蒸發(fā)富集的地表水對(duì)地下水補(bǔ)給頻繁。B5 到B7 位于中段區(qū)域,該區(qū)域的特點(diǎn)是城市化程度高,工農(nóng)業(yè)發(fā)達(dá),水利灌溉設(shè)施密集,這與埃及尼羅河流域密集農(nóng)業(yè)措施影響地下水穩(wěn)定同位素比率;水面積和灌溉面積的增加導(dǎo)致石羊河流域年蒸發(fā)量的增加研究類似[33-34],灌溉回流對(duì)地下水的補(bǔ)給可能是地下水同位素特征蒸發(fā)富集的一種解釋。下游B9和B10地下水d-excess向上波動(dòng),可能由于井較深(85 m),與地表水交換相對(duì)較少,蒸發(fā)分餾作用弱。B12非常接近艾比湖,湖水的交換作用可能導(dǎo)致同位素值偏高。J3 與J4 位于精河下游區(qū)域,與地表水轉(zhuǎn)換頻繁,且靠近艾比湖鹽區(qū),d-excess富集特征明顯。
圖6 博河和精河流域地下水d-excess沿程變化Fig.6 Groundwater d-excess variation in Bortala and Jing River Basin along the flow path
3.3.1 不同區(qū)域地下水水化學(xué)特征及沿程變化除博河B9(pH=6.39)和B10(pH=6.78)樣點(diǎn)外,整個(gè)流域地下水pH 范圍為7.21~9.23,屬于弱堿性;除鄰近艾比湖的ABH12樣點(diǎn)外,流域地下水含鹽量接近于1 dS·m-1。由均值可知,艾比湖周邊區(qū)域地下水各水化學(xué)指標(biāo)都大于博河和精河流域;由最小和最大值可知,艾比湖周邊區(qū)域地下水水化學(xué)參數(shù)變化范圍大于博河和精河流域地下水(表2)。
表2 艾比湖流域水化學(xué)參數(shù)特征Tab.2 Characteristics of water chemistry parameters in Ebinur Lake Basin
從水化學(xué)參數(shù)的沿程變化(圖7a~b)可知,博河、精河中下游和博河上游地下水水化學(xué)參數(shù)差異明顯,上游區(qū)域EC 均值小于200 μS·cm-1。中下游EC 均值大于200 μS·cm-1。艾比湖周邊地下水EC變化范圍最強(qiáng),博河次之,精河最弱。結(jié)合氫氧同位素特征,中下游區(qū)域和河湖交匯區(qū)地表水與地下水交換密切,同時(shí)考慮周邊農(nóng)田灌溉等人類活動(dòng)影響,這些區(qū)域是地下水防控和治理的重點(diǎn)區(qū)域。
圖7 艾比湖流域地下水水化學(xué)特征沿程變化Fig.7 Variation of hydrochemical characteristics of groundwater in Ebinur Lake Basin
3.3.2 不同區(qū)域地下水流動(dòng)特征結(jié)合朱世丹和雷米等[20,35]的研究發(fā)現(xiàn),除個(gè)別樣點(diǎn),精河、博河地表水的δ2H 與δ18O 從上游到下游逐漸偏正,表明地下水與地表水的流向基本一致。根據(jù)δ2H、δ18O 和d-excess,結(jié)合水化學(xué)指標(biāo)的變化范圍可將流域分為2 個(gè)不同的流動(dòng)系統(tǒng)(表3)。從EC 變化特征來(lái)看,流動(dòng)系統(tǒng)Ⅰ的離子含量可能要低于流動(dòng)系統(tǒng)Ⅱ,說(shuō)明2 個(gè)系統(tǒng)的地球化學(xué)環(huán)境不同。從同位素特征看,流動(dòng)系統(tǒng)Ⅰ的同位素d-excess 范圍比流動(dòng)系統(tǒng)Ⅱ更小。從井深來(lái)看,流動(dòng)系統(tǒng)Ⅰ地下水深度均在15~100 m 之間,而流動(dòng)系統(tǒng)Ⅱ地下水深度在30~260 m 之間,所以這2 個(gè)系統(tǒng)可能分屬于水力聯(lián)系不同的流動(dòng)系統(tǒng)。
表3 地下水流動(dòng)系統(tǒng)的氫氧同位素特征Tab.3 Isotopic characteristics of groundwater flow systems
本文結(jié)合線性回歸、雙端元混合模型和GIS 空間分析等方法,分析了艾比湖流域地下水氫氧同位素和水化學(xué)的空間異質(zhì)性,探討了艾比湖流域地下水補(bǔ)給特征以及影響補(bǔ)給的因素,揭示了艾比湖流域含水層的流動(dòng)特征。
(1)艾比湖流域不同區(qū)域地下水氫氧同位素存在空間異質(zhì)性??傮w來(lái)看,氫氧同位素值最大位于博河和精河中下游區(qū)域,艾比湖周邊區(qū)域次之,博河上游區(qū)域最小,流域地下水主要受到補(bǔ)給源和水巖作用的影響。水平分布上,博河、精河流域地下水d-excess 表現(xiàn)出聚集特征,而艾比湖周邊地下水表現(xiàn)出分異特征。垂直分布上,d-excess 沒(méi)有表現(xiàn)出深度效應(yīng)。
(2)艾比湖流域不同區(qū)域地下水存在不同的補(bǔ)給機(jī)制和補(bǔ)給特征。地下水δ2H 與δ18O 關(guān)系線的斜率和截距均小于天山北坡大氣降水線和流域大氣降水線,表明地下水受到一定程度的蒸發(fā)分餾。結(jié)合不同區(qū)域d-excess 特征可知,博河上游區(qū)域地下水主要受冰川積雪融水補(bǔ)給;博河和精河中下游地下水主要來(lái)源為地表水和降水,同時(shí)受巖層性質(zhì)、農(nóng)田開(kāi)發(fā)和灌溉措施影響較大;艾比湖周邊地下水主要來(lái)源于冰雪融水和降水,部分地下水受地表水補(bǔ)給。結(jié)合EC 變化特征可知,精河、博河中下游區(qū)域和河湖交匯區(qū)地下水是防控和治理的重點(diǎn)區(qū)域。
(3)不同區(qū)域地下水水化學(xué)特征沿程變化明顯。結(jié)合氫氧穩(wěn)定同位素差異性,判斷流域地下水存在2 個(gè)不同的流動(dòng)系統(tǒng)。流動(dòng)系統(tǒng)ⅠEC 在210.00~2500.00 μS·cm-1之 間,d-excess 在6.47‰~9.70‰之間,而流動(dòng)系統(tǒng)ⅡEC 在141.60~5260.00 μS·cm-1之間,d-excess 在9.61‰~17.45‰之間,表明含水層存在不同的水力聯(lián)系。
本文根據(jù)已有數(shù)據(jù)很好的揭示了艾比湖流域地下水特征,但是為了進(jìn)一步理清地下水多源補(bǔ)給區(qū)和補(bǔ)給方式、垂直分異規(guī)律,在如此大規(guī)模的含水層系統(tǒng)中,需要水力計(jì)算進(jìn)一步驗(yàn)證我們的結(jié)果。未來(lái)還要重復(fù)訪問(wèn)這些樣點(diǎn),增加樣本點(diǎn)和利用多種同位素示蹤,以理清地下水同位素水化學(xué)年內(nèi)、年際分異特征,以期為干旱半干旱地區(qū)跨流域調(diào)水、水權(quán)、入湖水量分配和流域治理等提供理論支撐。