环境科学  2024, Vol. 45 Issue (4): 2107-2117   PDF    
丹江流域山区地表水-地下水水化学特征及其影响因素
张子燕1,2, 伏永朋1, 王宁涛1,2, 谭建民1,2, 刘亚磊1     
1. 中国地质调查局武汉地质调查中心(中南地质科技创新中心), 武汉 430205;
2. 中国地质大学地质调查研究院, 武汉 430074
摘要: 丹江流域是南水北调中线工程重要水源涵养区, 深入研究丹江流域水化学特征及其影响因素, 对区域水资源保护及合理利用具有重要意义. 通过对丹江流域水样采集与分析, 综合运用描述性统计、Piper三线图、Gibbs图、离子比例系数法及氢氧同位素的方法, 分析丹江流域水化学特征及其影响因素. 结果表明, 丹江流域水体整体呈弱碱性, 主要水化学类型为HCO3-Ca型, 中游离子质量浓度及水化学类型变化频繁. 流域水体均受大气降水补给, 地表水受蒸发作用影响明显, 水化学组分受岩石风化溶解、离子交换和人类活动共同影响. Ca2+和Mg2+源于方解石、白云石及石膏等碳酸盐和硫酸盐矿物溶解, 农业活动是水中硝酸盐的主要来源. 地表水与地下水频繁转化和支流汇入影响是丹江干流空间分布特征差异化的主要原因.
关键词: 丹江流域      水化学特征      同位素      岩石风化      离子交换     
Hydrochemical Characteristics and Influencing Factors of Surface Water and Groundwater in the Mountainous Area of Danjiang River Basin
ZHANG Zi-yan1,2 , FU Yong-peng1 , WANG Ning-tao1,2 , TAN Jian-min1,2 , LIU Ya-lei1     
1. Wuhan Center, China Geological Survey (Central South China Innovation Center for Geosciences), Wuhan 430205, China;
2. Institute of Geological Survey, China University of Geosciences, Wuhan 430074, China
Abstract: The Danjiang River Basin is an important water conservation area in the middle route of the South-to-North Water Transfer Project. Studying the hydrochemical characteristics and influencing factors of the Danjiang River Basin is of great significance for the protection and rational utilization of regional water resources. In this study, water samples in the Danjiang River Basin were collected and analyzed. Descriptive statistics, Piper trilinear diagram, Gibbs diagram, ion ratio coefficient, and hydrogen and oxygen isotopes were utilized to analyze the hydrochemical characteristics and influencing factors of the Danjiang River Basin. The results showed that the water body in the study area was weakly alkaline, and the main hydrochemical type was HCO3-Ca. The ion mass concentration and hydrochemical types changed frequently in the middle reaches. The water bodies in the basin were all recharged by atmospheric precipitation, and surface water was obviously affected by evaporation. The hydrochemical compositions were influenced by rock weathering and dissolution, ion exchange, and human activities. Ca2+ and Mg2+ originated from the dissolution of carbonate and sulfate minerals such as calcite, dolomite, and gypsum. Agricultural activities were the main source of nitrate in the water. The frequent exchange of surface water and groundwater and the influence of tributary inflow were the main reasons for the spatial distribution differentiation of the Danjiang main stream.
Key words: Danjiang River Basin      hydrochemical characteristics      isotope      rock weathering      ion exchange     

水资源的日益紧缺使得水的经济价值和社会效益更加突出[1, 2], 水在经济社会发展中的作用也日益凸显[3, 4], 同时人类活动也影响着水资源质量[5 ~ 7]. 相较于平原区水资源, 山区地表水和地下水呈现出较大时空差异[8 ~ 10]. 开展自然和人为活动共同作用下山区地表水-地下水水化学特征及影响因素研究, 是山区水资源保护及合理利用的基础.

在流域尺度上, 国内外学者通常采用同位素结合水化学方法研究水化学特征[11]. 具体体现在, 水化学类型法用于查明地下水水质状况, 探明地下水水化学类型时空分布规律, 认识区域的水化学特征[12 ~ 14];离子比例系数分析法用于探明地下水化学成分变化原因[15, 16]. 在水文地质学中, 同位素技术主要应用于定年、确定物质来源及过程反演等[17]. 氢氧同位素组成特征在不同来源水体中差异显著, 因此氢氧同位素主要用于确定水的来源及成因类型[18].

丹江流域是南水北调中线工程丹江口水库的重要水源涵养区, 因此其水质成为国家重点关注对象. 丹江口水库始建于1958年, 2013年丹江口水库大坝加高后, 于2014年开始蓄水和调水[19]. 自水库开始蓄水调水后, 水源地周边环境在自然和人类活动影响下发生显著改变. 首先表现在水域面积明显增加, 水域空间分布改变;其次是土地覆被变化, 耕地和建设用地变化最为显著, 导致植被类型亦发生变化;除此之外, 由消落带空间范围变化引起的农业面源污染残留问题, 也给丹江流域的水质安全带来隐患[20 ~ 24]. 此外, 丹江流域地层分布广泛, 主要为元古代武当岩群、南华系和震旦系, 早古生代寒武系、奥陶系、志留系、晚古生代泥盆系和石炭系, 中生代白垩系及新生代古近系、新近系和第四系也皆有分布.

针对丹江流域水化学特征研究, 较多集中于丹江口库区范围[19 ~ 22], 部分以丹江流域内支流为研究对象展开讨论[23, 24]. 关注丹江流域全域的研究相对较少, 赵培等[25]研究也主要针对地表水化学特征进行讨论. 因此, 本文以丹江流域地表水、井水和泉水为研究对象, 在丹江流域全域范围内, 探讨地表水和地下水化学特征, 并进一步分析其影响因素和来源, 以期为认识丹江流域水资源保护及合理开发利用提供科学参考.

1 材料与方法 1.1 研究区概况

丹江是汉江最大的支流, 发源于秦岭南麓, 流经陕西、河南和湖北这3省, 最终汇入丹江口水库. 流域范围在109°30′~111°30′E, 32°30′~34°30′N, 海拔在100~2 164 m之间, 流域面积1.07万km2[23]. 流域内有老灌河、滔河、淇河、银花河、武关河等多条支流, 河流长度约430 km.

研究区气候类型为亚热带与温带的混合型气候, 多年平均气温在11~14℃[26]. 多年平均降雨量为791 mm, 降水时空分布不均, 主要集中在7~9月, 高山区降水量最多, 其次是低山丘陵区, 河谷川道地区降水最少[24, 26]. 主要土地利用类型为林地和草地, 土地利用方式在蓄水前后变化较大. 张田田等[27]研究表明, 2000~2018年丹江流域耕地面积占比不断下降, 林地面积先减少后增加, 草地则呈现先增加后减少再增加的趋势, 而水体面积在丹江口水库蓄水之后逐渐增加.

丹江流域地势整体西北高东南低, 地貌构成分别为:低山丘陵地貌、河谷川原地貌和中山地貌[26]. 丹江流域南北两侧为变质岩区, 中上游主要出露大面积砂砾岩, 下游地区主要出露大面积的碳酸盐岩.

1.2 样品采集和分析

本研究于2020年12月对丹江流域干流、主要一级支流及地下水开展取样工作, 完成取样82件, 其中包括51个地表水、21个井水和10个泉水, 具体分布如图 1所示. R1~R21为丹江干流上游至下游地表水样点, R22~R25、R26~R27、R28~R33、R34~R46和R47~R51分别为银花河、武关河、淇河、老灌河和滔河不同支流的样点, 井水和泉水样点编号分别为J1~J21和Q1~Q10. 根据河流位置, 本文将丹江干流分为上游商州-南召段(R1~R6)、中游南召-寺湾段(R7~R15)和下游寺湾-马镫段(R16~R21).

图 1 丹江流域水样点分布示意 Fig. 1 Distribution of sampling sites in Danjiang River Basin

本研究中, 水样分析由中国地质调查局武汉地质调查中心(中南地质科技创新中心)完成, K+、Na+、Ca2+、Mg2+、Cl-、SO42-和NO3-采用离子色谱仪(ICS-5000+DC, 美国)分析, HCO3-采用滴定法测定. δ2H和δ18O采用液态水同位素分析仪(Los Gatos Research IWA-35EP, 美国)进行测定, 采用VSMOW标准.

为探明不同水体水化学特征的影响因素, 采用Piper三线图分析水化学类型, 根据Gibbs图和离子比例系数等方法识别研究区水化学特征;利用氢氧同位素确定水的来源及成因类型.

2 结果与分析 2.1 水化学特征及水化学类型

丹江流域水化学组成特征见表 1. 整体来看, 丹江流域水体ρ(TDS)均小于1 000 mg·L-1, 其中井水ρ(TDS)平均值最大, 为532.9 mg·L-1, 泉水次之, 为351.1 mg·L-1, 地表水最小, 为296.9 mg·L-1. 根据TDS的变异系数发现, 井水(62.74%)与地表水(26.31%)和泉水(23.56%)的变异系数相差较大, 表明井水中溶解离子来源可能不同于地表水和泉水.

表 1 丹江流域地表水、井水和泉水水化学参数统计1) Table 1 Statistics of hydrochemical parameters of surface water, well water, and spring water in Danjiang River Basin

地表水、井水和泉水的pH变化范围分别为7.7~8.4、7.2~8.2和7.3~8.4, 平均值分别为8.1、7.6和7.8, 以中性至偏碱性环境为主, 地表水平均值更偏碱性. 地表水和地下水的主要阳离子为Ca2+, 主要阴离子为HCO3-.

阳离子中, 除了K+, 离子质量浓度平均值都表现出地下水高于地表水的规律;阴离子中, 井水中的离子质量浓度平均值始终高于其他水体类型;Na+、Ca2+、Mg2+和HCO3-的平均值变化规律同TDS平均值变化一致. 主要离子质量浓度变化幅度均呈现出地下水(井水、泉水)大于地表水的规律. 除Ca2+、Mg2+和HCO3-, 井水中Cl-、NO3-, 泉水中K+、Na+、SO42-和NO3-的变异系数均在100%以上, 进一步说明地下水(井水、泉水)中的主要离子来源丰富. 地表水中TN(223.04%)及地下水NO3-(井水:171.01%, 泉水:107.42%)的较高变异系数一定程度上也反映了人类活动对水质的影响.

ρ(SiO2)在地表水、井水和泉水中的变化范围分别为:2.35~9.79、5.73~26.70和8.16~75.60 mg·L-1, 平均值分别为6.16、3.35和5.18 mg·L-1, 泉水中的SiO2变异系数(113.99%)远大于地表水(28.85%)和井水(25.92%). 如图 2所示, 研究区主要分布碳酸盐岩类裂隙含水层、碎屑岩类裂隙含水层、岩浆岩裂隙含水层和变质岩裂隙含水层等含水岩组, 且地下水与地表水转换频繁, 因此泉水中可能存在硅酸盐的溶解.

图 2 水文地质剖面示意 Fig. 2 Hydrogeological profile in study area

利用Piper三线图了解丹江流域的水文地球化学机制. 如图 3所示, 大部分阳离子分布在Ca2+区, 其次是Mg2+区, 只有一个井水样品分布在Na++K+区;从阴离子看, 大部分阴离子样品集中在HCO3-区;整体来看, 大部分样品分布在CaHCO3区, 因此丹江流域主要水化学类型为HCO3-Ca型, 个别井水是混合型. Piper图显示, 地下水中Ca2+占比整体而言较地表水高, 这可能与地下水与围盐相互作用, 导致碳酸盐岩溶解有关.

图 3 丹江流域地表水、井水和泉水Piper图 Fig. 3 Piper map of surface water, well water and spring water in Danjiang River Basin

为明晰丹江流域水化学的空间分布特征, 沿丹江干流展开讨论. 如图 4所示, 结合水化学类型可知, 丹江上游商州-南召段(R1~R6)及下游寺湾-马镫段(R16~R21)水化学类型单一, 均为Ca·Mg-HCO3. 中游南召-寺湾段(R7~R15)离子质量浓度之和呈波形变化, 水化学类型变化频繁.

图 4 丹江干流和支流地表水离子质量浓度 Fig. 4 Ion concentrations in surface water of main stream and tributaries of Danjiang River

结合支流离子质量浓度变化情况发现, 丹江中游在武关河支流汇入后, R10处八大离子质量浓度之和急剧降低至265 mg·L-1, 水化学类型转变为Ca·Mg-HCO3. 此后八大离子质量浓度之和沿河流方向逐渐增加, 在R12处八大离子质量浓度之和达到丹江干流最大值, 为630 mg·L-1, 其中Na+、Ca2+、Cl-及HCO3-贡献较大.

2.2 同位素特征

表 2显示了丹江流域水体同位素的统计特征. 在地表水、井水和泉水中, δ2H和δ18O稳定同位素的平均值呈现出泉水最高, 地表水次之, 井水最低的趋势. 而在δ2H和δ18O稳定同位素的变化范围方面, 则呈现出井水最大, 地表水居中, 泉水最小的规律. 具体而言, 地表水、井水和泉水的δ2H和δ18O的分布范围分别为-63.1‰ ~ -56.5‰和-9.1‰ ~ -8.0‰、-67.2‰ ~-51.4‰和-10.0‰ ~ -7.8‰、-61.0‰ ~ -54.5‰和-8.9‰ ~ -8.1‰.

表 2 丹江流域山区地表水、井水和泉水氢氧稳定同位素和氘盈余统计特征/‰ Table 2 Statistical characteristics of hydrogen and oxygen stable isotopes and deuterium surplus of surface water, well water, and spring water in Danjiang River Basin/‰

氘盈余是由于水体蒸发破坏2H与18O之间平衡分馏产生的[28]. 氘盈余计算基于全球大气降水线[29], 是反映水体二次蒸发的重要参数[30, 31]. 研究区内井水和泉水氘盈余平均值大于10‰, 地表水氘盈余平均值低于10‰.

为更好揭示地表水、井水和泉水δ2H和δ18O稳定同位素之间的联系, 除统计特征外, 还对水样点的分布特征进行了分析. 如图 5所示, 本文中当地大气降水线(LMWL)参考周毅等[32]在长江流域的大气降水线:δ2H = 7.41 δ18O + 6.04. 大部分地下水样品位于全球大气降水线(GMWL)[29]上方, 地表水样品位于当地大气降水线下方. 运用最小二乘法拟合得到丹江流域山区地表水蒸发线为:δ2H = 6.110 1 δ18O-7.045 2(R2 = 0.825 4);地下水蒸发线为:δ2H = 6.970 2 δ18O+2.251 1(R2 = 0.963 0).

图 5 丹江流域地表水、井水和泉水δ2H-δ18O同位素关系 Fig. 5 The δ2H-δ18O isotopic relationship of surface water, well water, and spring water in Danjiang River Basin

3 讨论 3.1 主控因素

一般而言, 矿物风化、蒸发结晶和大气降水是影响地下水化学的主要因素, 这是因为不同矿物具有不同的溶解速率、交换和大气来源, 而Gibbs图可以表征以上水化学影响因素[33]. 图 6显示, 丹江流域水体水化学组分主要受岩石风化作用影响. 然而Gibbs图中相较于Na+/(Na++Ca2+), Cl-/(Cl-+HCO3-)相对较小, 这主要与水中的Ca2+和Mg2+与岩石中的Na+发生了交换作用有关[34], 这对于研究区水体中污染物的迁移及追溯Na+来源具有一定的指示作用.

图 6 丹江流域地表水和井水和泉水Gibbs图 Fig. 6 Gibbs diagram of surface water, well water, and spring water in Danjiang River Basin

3.1.1 矿物溶解与沉淀

绘制研究区离子成分比例系数图, 如图 7(a)所示, 大部分水样点分布于1∶1线附近, 部分井水样点分布于1∶1线上方. 根据Abdelkader等[34]研究表明, 若N(Ca2+ + Mg2+)/N(HCO3- + SO42-)接近于1, 则Ca2+和Mg2+应是来自于方解石、白云石及石膏等碳酸盐和硫酸盐矿物溶解;若比值大于1, 还应有硅酸盐溶解. 因此, 丹江流域山区地表水、井水和泉水水化学特征的主要影响因素为碳酸盐和硫酸盐的风化溶解, 其中井水的水化学特征还受硅酸盐溶解影响, 如J7、J8、J11、J17和J18均位于冲积层上, 井水成分可能在流经硅酸盐岩地层时发生改变, J20位于以变石英角斑岩为主的硅酸盐岩地层上, 推测井水成分受该地层影响.

图 7 丹江流域地表水、井水和泉水不同离子之间的关系 Fig. 7 Relationship between different ions in surface water, well water, and spring water in Danjiang River Basin

图 7(b)所示, 地表水和泉水样点基本位于2∶1线上方及1∶1与2∶1线中间, 井水样点基本位于2∶1线上方. 根据Mayo等[35]研究表明, 当N(Ca2+)/N(Mg2+)值接近1, 表明溶解的碳酸盐矿物主要为白云石, 方解石的溶解会增大该比值;若比值大于2, 则还存在硅酸盐的溶解[25]. 因此丹江流域水体溶解的碳酸盐矿物为方解石、白云石, 井水中硅酸盐风化占主导地位, 个别井水水样点位于1∶1线下方, 表明其主要受到方解石等贫镁矿物影响. 如图 7(c)所示, 大部分水样点位于1∶1线下方, 根据Jalali[36]研究表明, 当Na+-Cl-关系偏离1∶1线, 说明Na+来源于其他过程.

根据Tesoriero等[37]研究发现, 硅酸盐的风化、钾长石和钠长石的溶解是HCO3-的主要来源. 从图 7(d)中可知, 水样点基本位于1∶1线下方, 结合前面的讨论, 丹江流域地表水和地下水中的HCO3-主要来源于硅酸盐的风化.

3.1.2 离子交换

图 8(a)所示, 所有样品沿着斜率为-1.59的直线分布, 并且具有中等强度的相关性R2 = 0.687 4, 表明离子交换是丹江流域水质的影响因素之一, 其中Na+被释放, Ca2+被截留. 此外图 8(b)中显示, 相较于Na+-Cl-关系, 在Ca2+-SO42-关系中, Ca2+占主导地位的趋势更加明显, 表明离子交换过程中, 陆源性SO42-与Ca2+结合沉淀, 导致石膏欠饱和. 根据图 8(c)中HCO3--Na+关系所示, 丹江流域Na+另一来源为碳酸盐矿物的溶解. 图 8(d)显示, 所有样品都位于1∶1线以下, 进一步表明Na++K+对研究区水化学的主要贡献.

图 8 丹江流域地表水、井水及泉水离子比例关系 Fig. 8 Ion proportion of surface water, well water, and spring water in Danjiang River Basin

氯碱指数(CAI1和CAI2)反映阳离子交替吸附的方向和强弱;两者均为正, 表明水中的Na+和K+与土壤中的Ca2+和Mg2+交换;两者均为负, 表明水中的Ca2+和Mg2+与土壤中的Na+和K+发生交换作用;指数绝对值大小代表交替强弱. 图 9显示84%样点的氯碱指数值为负值, 且CAI1绝对值高, CAI2接近0, 表明水中Ca2+-Mg2+与围岩中的Na+-K+发生交换.

图 9 丹江流域地表水、井水和泉水的氯碱指数 Fig. 9 Chlor-alkali index map of surface water, well water, and spring water in Danjiang River Basin

3.1.3 人类活动影响

硝酸盐是地下水中最普遍的污染物之一[38 ~ 40], 长期饮用含高浓度硝酸盐的水, 会增加人类罹患胃癌的可能性[41]. ρ(NO3-)在天然水体中一般低于10 mg·L-1, 但其具有高溶解性和易于迁移的特点, 随着人类活动影响会导致水体中NO3-质量浓度增加[42], 因此水中NO3-一般可表征人类活动对水化学的影响[20]. 农业上肥料及农药不当使用, 生活污水和工业废水的任意排放等均是地下水硝酸盐污染的重要原因[43 ~ 45]. 根据前面的讨论可知, 硝酸盐含量分布变化规律为:井水 > 地表水 > 泉水. 以往研究也表明, 丹江流域内土地利用类型大部分经由耕地转化为其他类型[19 ~ 23, 27], 因此丹江流域内水体中硝酸盐可能来源于耕地种植过程中使用的农业氮肥.

3.2 不同水体补给来源及蒸发强度

δ18O与Cl-的关系可以指示地下水的水文运移作用[46, 47]. 水体δ18O与Cl-关系存在以下趋势:①Cl-质量浓度急剧增加, δ18O值变化不大, 是溶滤与蒸发混合作用的结果[48];②随着Cl-质量浓度的适度变化, δ18O值急剧增加. 这反映了蒸发作用对水体的影响, 因为蒸发作用可以在适度改变Cl-的情况下改变水体的δ18O值[49];③δ18O值不伴随Cl-质量浓度变化而变化.

Cl-δ18O的关系如图 10所示, 研究区地表水主要呈现出Cl-质量浓度适度变化, δ18O富集的特征, 因此地表水受蒸发作用较大. 根据δ18O的富集情况, 井水的Cl-δ18O关系可分为2组, 第一组δ18O富集, Cl-质量浓度接近, 该井水样点均位于南秦岭地层分区上, 含水层类型均为碎屑岩类裂隙孔隙水;第二组整体上Cl-质量浓度随δ18O值增加而增大, 4个井水样点分属4种含水层类型, 表明该变化可能是溶滤作用引起的;J4样点处Cl-质量浓度显著高于其他井水, 推测周边人类活动可能是该水样点Cl-质量浓度偏高的原因. 研究区泉水δ18O值的变化不伴随Cl-质量浓度变化, 表明泉水可能是由于侧向补给水与深层地下水混合所致.

图 10 丹江流域地表水、井水和泉水δ18O与Cl-关系 Fig. 10 Relationship between δ18O and Cl- of surface water, well water, and spring water in Danjiang River Basin

地表水、井水和泉水的δ2H和δ18O同位素关系显示, 这3种水体同位素组成接近大气降水线且同位素变化范围相似, 表明均接受大气降水补给且水力联系密切. 降水过程中强烈的蒸发会导致降水氘盈余变小[50]. 如图 11所示, 丹江流域地表水的蒸发作用较为强烈, 与Cl-δ18O关系得出的结果一致.

图 11 丹江流域地表水、井水和泉水d-excess与δ18O关系 Fig. 11 Relationship between d-excess and δ18O surface water, well water, and spring water in Danjiang River Basin

3.3 地表水与地下水转化关系

通过分析研究区水体水化学及同位素特征, 发现丹江中游和下游水化学特征空间差异显著. 以TDS和δ18O为示踪剂, 分析丹江中下游地表水和地下水转化关系.

图 12所示, 丹江中游地下水中TDS与δ18O呈波动型变化, 推测丹江中游地下水补给来源及循环过程存在差异. J12水样点处TDS值升高但δ18O迅速贫化, 推断其可能接受侧向径流补给;Q1水样点处TDS降低, δ18O富集且值高于地表水, 推断为接受地表水和大气降水补给;J4水样点处TDS显著升高而δ18O迅速贫化, 推断其可能接受侧向径流补给, 且受到人类活动干扰较大. Q8、Q9及J15水样点处TDS及δ18O呈增加趋势, 且表现出地下水值高于地表水, 两者之间存在一定的转化, 但补给关系还需借助水力梯度进一步明确.

图 12 研究区地下水与地表水δ18O与TDS沿程变化 Fig. 12 Changes in δ18O and TDS of groundwater and surface water in study area

丹江下游主要由老灌河支流汇入, 以老灌河为例研究丹江下游地表水与地下水转化关系. 该区域整体上TDS及δ18O呈下降趋势, 水样点J21处TDS及δ18O值与地表水接近, 推测该处地表水补给地下水;水样点Q4处δ18O值降低, 推测可能接受δ18O更为贫化的侧向径流补给.

4 结论

(1)丹江流域水体均受大气降水补给, 地表水受蒸发作用影响明显. 研究区大部分地表水和地下水呈弱碱性, HCO3-和Ca2+是不同水体的主要阴阳离子, 水化学组分受矿物溶解、离子交换和人类活动共同制约.

(2)研究区地下水水化学组成主要受岩石风化影响. Ca2+和Mg2+来自于方解石、白云石和石膏等碳酸盐和硫酸盐矿物的溶解;农业活动是水中硝酸盐的主要来源.

(3)研究区干流水化学特征空间差异显著, 表现为水化学类型复杂、离子质量浓度变化剧烈, 主要与地下水与地表水的频繁转化以及支流汇入有关.

参考文献
[1] 岳思羽. 水资源短缺地区河道生态基流的价值与时空变化研究——以渭河为例[D]. 西安: 西安理工大学, 2020.
[2] 黄硕俏. 黄河流域水资源经济价值空间分布特征及其影响因素分析[D]. 郑州: 郑州大学, 2020.
[3] 傅懿兵. 浅谈水资源科学管理对促进经济可持续发展的重要性[J]. 给水排水, 2022, 58(8): 212-213.
[4] 黄怡, 效存德, 苏勃, 等. 冰冻圈水资源对中国社会经济的重要性: 基于多级流域尺度分析[J]. 科学通报, 2020, 65(24): 2636-2650.
Huang Y, Xiao C D, Su B, et al. Importance of cryosphere water resources to China's socio-economy: based on multi-level basins analysis[J]. Chinese Science Bulletin, 2020, 65(24): 2636-2650.
[5] 黄昌硕, 耿雷华, 颜冰, 等. 水资源承载力动态预测与调控——以黄河流域为例[J]. 水科学进展, 2021, 32(1): 59-67.
Huang C S, Geng L H, Yan B, et al. Dynamic prediction and regulation of water resource carrying capacity: a case study on the Yellow River Basin[J]. Advances in Water Science, 2021, 32(1): 59-67.
[6] 陈吟, 赵莹, 杜鹏飞, 等. 近60年无定河流域降雨及人类活动对水沙变化的影响[A]. 见: 2022中国水利学术大会论文集(第一分册)[C]. 北京: 中国水利学会, 2022.
[7] 刘强, 尉飞鸿, 常康飞, 等. 黄河中游龙潼区间水沙变化特征及其对降水和人类活动的响应[J]. 应用基础与工程科学学报, 2023, 31(3): 555-568.
Liu Q, Yu F H, Chang K F, et al. Variation character of water and sediment and its response to rainfall and human activities in Longtong Region of the Middle Yellow River[J]. Journal of Basic Science and Engineering, 2023, 31(3): 555-568.
[8] Sun C J, Wang S Y, Chen W. Hydrochemical characteristics and the relationship between surface and groundwater in a typical 'mountain-oasis' ecosystem in central Asia[J]. Sustainability, 2022, 14(12). DOI:10.3390/su14127453
[9] 王继华, 于福荣, 田鹏州, 等. 南太行山前冲积平原潜水水化学特征及控制因素分析[J]. 干旱区资源与环境, 2023, 37(5): 130-138.
Wang J H, Yu F R, Tian P Z, et al. Analysis on hydrochemical characteristics of groundwater in alluvial plain of southern Taihang Mountains[J]. Journal of Arid Land Resources and Environment, 2023, 37(5): 130-138.
[10] 马惠莹. 石羊河流域山区-绿洲-荒漠不同水体水化学时空变化及控制因素[D]. 兰州: 西北师范大学, 2022.
[11] 雷米, 周金龙, 张杰, 等. 新疆博尔塔拉河流域平原区地表水与地下水水化学特征及转化关系[J]. 环境科学, 2022, 43(4): 1873-1884.
Lei M, Zhou J L, Zhang J, et al. Hydrochemical characteristics and transformation relationship of surface water and groundwater in the plain area of Bortala River Basin, Xinjiang[J]. Environmental Science, 2022, 43(4): 1873-1884.
[12] Lu S S, Zhou N Q, Jiang S M, et al. Combining hydrochemistry and environmental isotopes to study hydrogeochemical evolution of karst groundwater in the Jinci spring area, North China[J]. Carbonates and Evaporites, 2023, 38(2). DOI:10.1007/s13146-023-00859-9
[13] 张保祥, 张超. 水文地球化学方法在地下水研究中的应用综述[J]. 人民黄河, 2019, 41(10): 135-142.
Zhang B X, Zhang C. Progress on hydrogeochemical method applied in groundwater study[J]. Yellow River, 2019, 41(10): 135-142.
[14] 张广禄, 刘海燕, 郭华明, 等. 华北平原典型山前冲洪积扇高硝态氮地下水分布特征及健康风险评价[J]. 地学前缘, 2023, 30(4): 485-503.
Zhang G L, Liu H Y, Guo H M, et al. Occurrences and health risks of high-nitrate groundwater in typical piedmont areas of the North China Plain[J]. Earth Science Frontiers, 2023, 30(4): 485-503.
[15] Sheng D R, Meng X H, Wen X H, et al. Hydrochemical characteristics, quality and health risk assessment of nitrate enriched coastal groundwater in northern China[J]. Journal of Cleaner Production, 2023, 403. DOI:10.1016/j.jclepro.2023.136872
[16] Hu M Y, Zhou P J, Chen C Q. Hydro-geochemical evolution of groundwater in the central Yangtze River Basin, China[J]. Carbonates and Evaporites, 2023, 38(2). DOI:10.1007/s13146-023-00852-2
[17] 谢先军, 甘义群, 刘运德, 等. 环境同位素原理与应用[M]. 北京: 科学出版社, 2019.
[18] Du G L, Zhai B, Su S J, et al. Hydrogeochemical characteristics and the genesis of the No. Lu 32 well in the Jiaodong Peninsula, China[J]. Frontiers in Earth Science, 2023, 11. DOI:10.3389/feart.2023.1131539
[19] 李冰, 王亚, 郑钊, 等. 丹江口水库调水前后表层沉积物营养盐和重金属时空变化[J]. 环境科学, 2018, 39(8): 3591-3600.
Li B, Wang Y, Zheng Z, et al. Temporal and spatial changes in sediment nutrients and heavy metals of the Danjiangkou Reservoir before and after water division of the Mid-route Project[J]. Environmental Science, 2018, 39(8): 3591-3600.
[20] 杨潇, 方子樊, 何文武. 丹江口库周典型支流污染离子的迁移模式及机理研究[J]. 资源环境与工程, 2022, 36(2): 204-209, 243.
Yang X, Fang Z F, He W W. Study on migrational model and mechanism of pollution Ions in typical tributaries around Danjiangkou Reservoir[J]. Resources Environment & Engineering, 2022, 36(2): 204-209, 243.
[21] 徐慧. 丹江口水库调蓄水效应遥感监测分析——以土地覆被变化为例[D]. 武汉: 湖北大学, 2017.
[22] 王磊, 章昱, 徐帅, 等. 丹江口水库总磷浓度遥感反演及其时空特征研究[J]. 华南地质与矿产, 2019, 35(4): 449-456.
Wang L, Zhang Y, Xu S, et al. Remote sensing retrieval of total phosphorus concentration in Danjiangkou Reservoir and its spatial and temporal characteristics[J]. South China Geology, 2019, 35(4): 449-456.
[23] 王一然, 丁文峰, 张冠华. 基于SWAT模型的丹江流域土地利用变化对径流影响研究[J]. 水土保持研究, 2022, 29(4): 62-67, 74.
Wang Y R, Ding W F, Zhang G H. Study on the impact of land use change on runoff in Danjiang watershed based on SWAT model[J]. Research of Soil and Water Conservation, 2022, 29(4): 62-67, 74.
[24] 李明月. 秦岭丹江流域水源涵养与土壤保持功能评价[D]. 西安: 西北大学, 2021.
[25] 赵培, 杨维鸽, 李晓刚. 丹江河流枯水期水化学特征及离子来源分析[J]. 甘肃科学学报, 2023, 35(2): 26-30.
Zhao P, Yang W G, Li X G. Analysis on hydrochemical characteristics and iron sources in dry season of Dan River[J]. Journal of Gansu Sciences, 2023, 35(2): 26-30.
[26] 李秀清. 基于VIC模型的丹江流域水文模拟及水资源管理对策[D]. 西安: 西北大学, 2021.
[27] 张田田, 陈有超, 李潜, 等. 土地利用变化对丹江流域径流和泥沙时空格局的影响[J]. 长江流域资源与环境, 2022, 31(8): 1797-1811.
Zhang T T, Chen Y C, Li Q, et al. Effects of land-use change on the spatio-temporal patterns of runoff and sediment in the Danjiang River Basin[J]. Resources and Environment in the Yangtze Basin, 2022, 31(8): 1797-1811.
[28] Dansgaard W. Stable isotopes in precipitation[J]. Tellus, 1964, 16(4): 436-468. DOI:10.1111/j.2153-3490.1964.tb00181.x
[29] Craig H. Isotopic variations in meteoric waters[J]. Science, 1961, 133(3465): 1702-1703. DOI:10.1126/science.133.3465.1702
[30] Aemisegger F, Pfahl S, Sodemann H, et al. Deuterium excess as a proxy for continental moisture recycling and plant transpiration[J]. Atmospheric Chemistry and Physics, 2014, 14(8): 4029-4054. DOI:10.5194/acp-14-4029-2014
[31] Natali S, Doveri M, Giannecchini R, et al. Is the deuterium excess in precipitation a reliable tracer of moisture sources and water resources fate in the western Mediterranean? new insights from Apuan Alps (Italy)[J]. Journal of Hydrology, 2022, 614. DOI:10.1016/j.jhydrol.2022.128497
[32] 周毅, 吴华武, 贺斌, 等. 长江水δ18O和δD时空变化特征及其影响因素分析[J]. 长江流域资源与环境, 2017, 26(5): 678-686.
Zhou Y, Wu H W, He B, et al. Study on spatial and temporal variations of δ18O and δD in Yangtze river water and its factors[J]. Resources and Environment in the Yangtze Basin, 2017, 26(5): 678-686.
[33] Feth J H. Mechanisms controlling world water chemistry: evaporation-crystallization process[J]. Science, 1971, 172(3985): 870-871. DOI:10.1126/science.172.3985.870
[34] Abdelkader R, Larbi D, Rihab H, et al. Geochemical characterization of groundwater from shallow aquifer surrounding Fetzara Lake N. E. Algeria[J]. Arabian Journal of Geosciences, 2012, 5(1): 1-13. DOI:10.1007/s12517-010-0202-6
[35] Mayo A L, Loucks M D. Solute and isotopic geochemistry and ground water flow in the central Wasatch Range, Utah[J]. Journal of Hydrology, 1995, 172(1-4): 31-59. DOI:10.1016/0022-1694(95)02748-E
[36] Jalali M. Chemical characteristics of groundwater in parts of mountainous region, Alvand, Hamadan, Iran[J]. Environmental Geology, 2006, 51(3): 433-446. DOI:10.1007/s00254-006-0338-6
[37] Tesoriero A J, Spruill T B, Eimers J L. Erratum to "geochemistry of shallow ground water in coastal plain environments in the southeastern United States: implications for aquifer susceptibility" [Applied Geochemistry 19 (2004) 1471–1482][J]. Applied Geochemistry, 2005, 20(3). DOI:10.1016/j.apgeochem.2004.10.002
[38] Assa B G, Bhowmick A, Cholo B E. Modeling nitrogen balance for pre-assessment of surface and groundwater nitrate (NO3--N) contamination from N-fertilizer application loss: a case of the Bilate Downstream Watershed Cropland[J]. Water, 2023, 234(2). DOI:10.1007/s11270-023-06114-0
[39] Wakida F T, Lerner D N. Non-agricultural sources of groundwater nitrate: a review and case study[J]. Water Research, 2005, 39(1): 3-16. DOI:10.1016/j.watres.2004.07.026
[40] Yang P H, Wang Y Y, Wu X Y, et al. Nitrate sources and biogeochemical processes in karst underground rivers impacted by different anthropogenic input characteristics[J]. Environmental Pollution, 2020, 265. DOI:10.1016/j.envpol.2020.114835
[41] Fan A M, Steinberg V E. Health implications of nitrate and nitrite in drinking water: an update on methemoglobinemia occurrence and reproductive and developmental toxicity[J]. Regulatory Toxicology and Pharmacology, 1996, 23(1): 35-43. DOI:10.1006/rtph.1996.0006
[42] Whitehead P, Bussi G, Hossain M A, et al. Restoring water quality in the polluted Turag-Tongi-Balu River system, Dhaka: modelling nutrient and total coliform intervention strategies[J]. Science of the Total Environment, 2018, 631-632: 223-232. DOI:10.1016/j.scitotenv.2018.03.038
[43] 任坤, 潘晓东, 梁嘉鹏, 等. 碳氮氧同位素解析典型岩溶流域地下水中硝酸盐来源与归趋[J]. 环境科学, 2021, 42(5): 2268-2275.
Ren K, Pan X D, Liang J P, et al. Sources and fate of nitrate in groundwater in a typical karst basin: insights from carbon, nitrogen, and oxygen isotopes[J]. Environmental Science, 2021, 42(5): 2268-2275.
[44] 李霄, 王晓光, 柴璐, 等. 沉积盆地地下水无机氮来源示踪及其演化模式[J]. 中国环境科学, 2021, 41(4): 1856-1867.
Li X, Wang X G, Chai L, et al. Sources tracing and evolution model of inorganic nitrogen of groundwater in sedimentary basin[J]. China Environmental Science, 2021, 41(4): 1856-1867.
[45] 殷超, 杨海全, 陈敬安, 等. 基于水化学和氮氧同位素的贵州草海丰水期水体硝酸盐来源辨析[J]. 湖泊科学, 2020, 32(4): 989-998.
Yin C, Yang H Q, Chen J A, et al. Tracing nitrate sources with dual isotopes and hydrochemical characteristics during wet season in Lake Caohai, Guizhou Province[J]. Journal of Lake Sciences, 2020, 32(4): 989-998.
[46] Xie X J, Wang Y X. Su C L.et al. Influence of irrigation practices on arsenic mobilization: evidence from isotope composition and Cl/Br ratios in groundwater from Datong Basin, northern China[J]. Journal of Hydrology, 2012, 424-425: 37-47. DOI:10.1016/j.jhydrol.2011.12.017
[47] 曾邯斌, 苏春利, 谢先军, 等. 河套灌区西部浅层地下水咸化机制[J]. 地球科学, 2021, 46(6): 2267-2277.
Zeng H B, Su C L, Xie X J, et al. Mechanism of salinization of shallow groundwater in western Hetao irrigation area[J]. Earth Science, 2021, 46(6): 2267-2277.
[48] 徐蓉桢, 魏世博, 李成业, 等. 基于水化学与环境同位素的额济纳平原区域地下水循环规律解析[J]. 地学前缘, 2023, 30(4): 440-450.
Xu R Z, Wei S B, Li C Y, et al. Groundwater circulation in the Ejina Plain: insights from hydrochemical and environmental isotope studies[J]. Earth Science Frontiers, 2023, 30(4): 440-450.
[49] 李静, 朱佳俊, 梁杏, 等. 海州湾弱透水层孔隙水的化学特征与盐分演变[J]. 地球科学, 2022.
Li J, Zhu J J, Liang X, et al. Chemical signatures and salinity evolution of aquitard porewater in the Haizhou Bay[J]. Earth Science, 2022. DOI:10.3799/dqkx.2022.232
[50] Li Z X, Feng Q, Wang Q J, et al. Contributions of local terrestrial evaporation and transpiration to precipitation using δ18O and D-excess as a proxy in Shiyang inland river basin in China[J]. Global and Planetary Change, 2016, 146: 140-151. DOI:10.1016/j.gloplacha.2016.10.003