Hydrogeological investigation of karst system properties by common use of diverse methods: a case study of Lička Jesenica springs in Dinaric karst of Croatia

Li?ka Jesenica is a sinking river situated in the Dinaric karst environment of the Lika region of Croatia. The two main karst springs, Veliko Vrelo and Malo Vrelo, provide the majority of Li?ka Jesenica's water. Because of the quality and abundance of the water they contain, these springs are strategically important to the public water supply. Previous knowledge regarding the springs` karst system has been negligible. Therefore, a bulk hydrogeological research program was conducted with the purposes of establishing protection zones and assessing the springs' water‐supply potential. The research work comprised the following methods: hydrogeological mapping, setting up of monitoring stations for evaluating discharge and water quality, time series analyses (hydrograph–chemograph), hydrochemical analyses, and geophysical surveys. The main results and conclusions of the research include the following: (i) confirmation that both springs drain a common aquifer; (ii) determination of the highly vulnerable nature of the karst systems; (iii) identification of the Veliko Vrelo Spring's more direct connection to the aquifer‐drainage system in comparison to the Malo Vrelo; (iv) estimation of the positions and depths of the springs` main feeding conduits; (v) confirmation of the extent of the geologically presumed catchment area; and (vi) finally, establishment of sanitary protection zones. Copyright © 2012 John Wiley & Sons, Ltd.     

Assessment of annual high‐water events for the Mackenzie River basin,Canada

River ice break‐up is known to have important morphological, ecological and socio‐economic effects on cold‐regions river environments. One of the most persistent effects of the spring break‐up period is the occurrence of high‐water events. A return‐period assessment of maximum annual nominal water depths occurring during the spring break‐up and open‐water season at 28 Water Survey of Canada hydrometric sites over the 1913–2002 time period in the Mackenzie River basin is presented. For the return periods assessed, 13 (14) stations are dominated by peak events occurring during the spring break‐up (open‐water) season. One location is determined to have a mixed signal. A regime classification is proposed to separate ice‐ and open‐water dominated systems. As part of the regime classification procedure, specific characteristics of return‐period patterns including alignment, and difference between the 2 and 10‐year events are used to identify regime types. A dimensionless stage‐discharge plot allows for a contrast of the relative magnitudes of flows required to generate maximum nominal water‐depth events in the different regimes. At sites where discharge during the spring break‐up is approximately one‐quarter or greater than the magnitude of the peak annual discharge, nominal water depths can be expected to exceed those occurring during the peak annual discharge event. Several physical factors (location, basin area, stream order, gradient, river orientation, and climate) are considered to explain the differing regimes and discussed relative to the major sub‐regions of the MRB. Copyright © 2008 John Wiley & Sons, Ltd and Her Majesty the Queen in right of Canada.     

A high‐altitude temporary spring in a compartmentalized carbonate aquifer: the role of low‐permeability faults and karst conduits

The aim of this research was to refine the actual conceptual model related to the activation of high‐altitude temporary springs within the carbonate Apennines in southern Italy. The research was carried out through geophysical, hydrogeological, hydrochemical and isotopic investigations at the Acqua dei Faggi experimental site during five hydrologic years. The research demonstrated that, in carbonate aquifers where low‐permeability faults cause the aquifer system to be compartmentalized, high‐altitude temporary springs may be recharged by groundwater. In such settings, neither surface water infiltration in karst systems nor perched temporary aquifers play a role of utmost importance. The rare (once or a few time a year) activation of such springs is due to the fact that groundwater unusually reach the threshold head that allows the spring to flow. The activation of the studied high‐altitude temporary spring also depended on relationships between a low‐permeability fault core and a karst system that locally interrupts the low‐permeability barrier. In fact, when the hydraulic head did not reach the karst system, the concentrated head loss within the fault core did not allow the spring to flow, because the groundwater entirely flowed through the fault towards the downgradient compartment. Copyright © 2009 John Wiley & Sons, Ltd.     

THM analysis of a large‐scale heating test incorporating material fabric changes

Engineered barriers are basic elements in the design of repositories for the isolation of high‐level radioactive waste. This paper presents the thermo‐hydro‐mechanical (THM) analysis of a clay barrier subjected to heating and hydration. The study focuses on an ongoing large‐scale heating test, at almost full scale, which is being carried out at the CIEMAT laboratory under well‐controlled boundary conditions. The test is intensely instrumented and it has provided the opportunity to study in detail the evolution of the main THM variables over a long period of time. Comprehensive laboratory tests carried out in the context of the FEBEX and NF‐PRO projects have allowed the identification of the model parameters to describe the THM behaviour of the compacted expansive clay. A conventional THM approach that assumes the swelling clay as a single porosity medium has been initially adopted to analyse the evolution of the test. The model was able to predict correctly the global THM behaviour of the clay barrier in the short term (i.e. for times shorter than three years), but some model limitations were detected concerning the prediction of the long‐term hydration rate. An additional analysis of the test has been carried out using a double structure model to describe the actual behaviour of expansive clays. The double structure model explicitly considers the two dominant pore levels that actually exist in the FEBEX bentonite and it is able to account for the evolution of the material fabric. The simulation of the experiment using this enhanced model provides a more satisfactory reproduction of the long‐term experimental results. It also contributes to a better understanding of the observed test behaviour and it provides a physically based explanation for the very slow hydration of the barrier. Copyright © 2011 John Wiley & Sons, Ltd.     

安徽滁河断裂带温泉的水化学和同位素特征及成因分析

笔者等以安徽滁河断裂带内的6个温泉为主要研究对象,分析了水样的水化学特征,利用氢氧同位素对温泉的补给高程进行估算,并提出温泉的成因模式。研究区温泉阳离子以Ca2+和Mg2+为主,根据SO2-4和HCO-3的相对含量的不同,可以将水样分为两组,A组水样（富HCO-3）的主要离子的质量浓度均低于B组（富SO2-4）水样,A组水样的水化学类型为HCO-3—Ca2+·Mg2+; B组水样的（除AH14为SO2-4—Ca2+外）水化学类型为SO2-4—Ca2+·Mg2+。A组水样的稀土元素含量高于B组,二者均在NASC标准化图解上表现出平坦型的配分模式,且都表现出轻稀土富集和Eu正异常的特征。水样的氢氧稳定同位素组成表明温泉的补给来源都是大气降水,补给区温度约为13~15 ℃。A组温泉的补给高程为120~160 m低于B组温泉的200~260 m,且A组温泉的热储温度为45~70℃,低于B组温泉的热储温度70~105 ℃。地下水经历深循环获得大地热流加热后沿断裂带上升出地表。     

黄河支流汾河流域水资源开发利用现状及生态环境问题

【研究目的】汾河是黄河第二大支流,也是山西省的第一大河,流域内水资源供需矛盾突出,分析水资源开发利用现状及其生态环境问题是进行流域生态修复的前提。【研究方法】本文在分析汾河流域水资源特征及其开发利用现状的基础上,系统总结了汾河径流量衰减、岩溶大泉断流和水质恶化等生态环境问题,并对其成因进行了分析。【研究结果】研究表明：汾河流域多年平均水资源量为33.59亿m³,其中地下水资源是水资源的主要组成部分,约占72%;2005年以后由于跨流域调水、地下水压采等汾河流域综合治理措施的实施,水资源的供水结构发生了较大的变化,地表水的供水比例由最初的30%提高到55%,地下水供水比例由原来的62%降低到目前的37%。整体上,汾河流域的水资源开发利用程度高达80%以上,水资源的过度开发已导致汾河干流断流、入黄径流量大幅衰减、岩溶大泉断流等严重的生态环境问题。其中,汾河流入黄河径流量从1955至2018年衰减程度达63.5%,衰减的原因主要是降水量的减少和岩溶大泉的流量衰减;汾河流域内8个岩溶大泉的总流量从1956至2018年的衰减程度达69%,50%的岩溶大泉已在不同时期断流,岩溶大泉的水质恶化问题也非常严重,如晋祠泉和龙子祠泉的TDS和SO₄^2-呈逐年升高的趋势,煤矿开采是造成岩溶泉水SO₄^2-含量快速升高的主要原因。【结论】汾河流域的水资源供需矛盾十分突出,虽然通过跨流域调水等生态修复措施实现了汾河干流全年不断流、地下水位止降回升和地表水环境质量初步改善,但生态环境恶化的趋势依然严峻。创新点：分析了山西省汾河流域近15年的水资源特征及其供水结构的变化规律;系统总结了汾河流域的生态环境问题,并对其成因进行了探讨。     

藏南隆子县模麓温泉群水文地球化学特征及成因机制研究

藏南地区地热资源丰富,是喜马拉雅地热带的重要组成部分,有望成为新的地热资源开发靶区。本文以藏南桑日-错那活动构造带内模麓温泉群为研究对象,以水化学和氢氧氚同位素为研究方法,分析模麓温泉群的水岩作用、热储温度、补给来源及径流时间,揭示了地热水的成因机制。模麓地热水pH在6.6～7.2之间,TDS为1 908mg/L～2 326 mg/L,水化学类型以HCO3·Cl-Na型和HCO3·Cl-Na·Ca型为主。地热水中主要阴阳离子来源于硅酸盐矿物的溶解和少量地球深部物质。利用硅-焓方程法和硅-焓图解法计算的初始热储温度为198℃～256℃,冷水混入比例为68%～85%。此外,对地热水中的Li、B、F等微量元素分析得出,研究区温泉水中微量组分除来自水-岩作用外,应该还与深部流体的混入有关,且该地区的氢氧同位素特征表明地下水补给主要来源于大气降水,补给高程为5 652m～5 664m,模麓地热水中的氚含量<0.5TU,表明其地热水为老水,有更长的径流时间,为水-岩作用提供了充足的时间,而宿麦郎曲河水为新水,补给径流时间短。研究区地热水与围岩遮拉组砂板岩发生水-岩作用,进行离子交换作用,在地...     

四川稻城地区温泉流体地球化学特征及地震活动性分析

稻城地区位于青藏高原东南缘的川滇地块北部,为揭示该区域温泉流体地球化学特征以及其和地震活动性之间的关系,本次研究采集了稻城地区6个温泉的水样以及逸出气体样品。对温泉水中离子组分和浓度,温泉逸出气体组分及气体同位素进行了测试,得到以下认识。研究区温泉水化学类型主要为HCO₃-Na和HCO₃-Na·Ca型,通过阳离子温标估算热储温度在74℃～159℃之间,循环深度在2.2 km～5.0 km之间。温泉气体中CO₂主要是由储层中的碳酸盐岩受热分解或溶解产生的,氦来自幔源组分的比例较低,在0.4%～2.4%之间,研究区的温泉是由沿断裂带渗入的大气降水经地壳深部的热源加热形成的。在稻城地区,循环深度以及幔源气体贡献率等不同的温泉流体地球化学特征与地震活动性之间有很好的对应关系,并且研究区地震活动性要弱于周缘鲜水河断裂地区等深部流体上涌地区。同时在区域尺度上,未来位于断裂交汇部位的稻城仲堆温泉一带的地震活动性最值得关注。     

雅鲁藏布江色东普沟冰崩—堵江—溃决灾害数值模拟

为了揭示雅鲁藏布江色东普沟2018年10月17日冰崩—堵江—溃决灾害链的动力演化过程,基于Massflow数值模拟仿真平台,使用Fortran编程语言,根据研究区域地质条件特征对程序进行二次开发以优化Voellmy模型,模拟冰崩—泥石流动力过程;将模拟泥石流得到的堰塞坝体嵌入地形中,运用ArcGIS计算堰塞湖范围及体积,通过Manning模型模拟堰塞湖溃决洪水动力过程。采用分段模拟法再现冰崩—泥石流—堵江—堰塞湖—溃坝的完整动力过程,对泥石流运动过程中的流速、流深,坝体高度,溃决洪水的流深、流速等参数进行定量化研究,为色东普流域的防灾减灾工作提供有效支撑。为了揭示雅鲁藏布江色东普沟2018年10月17日冰崩-堵江-溃决灾害链的动力演化过程,采用Massflow数值模拟仿真平台,以Fortran语言为编程手段,根据研究区域地质条件特征对程序进行二次开发优化Voellmy模型,模拟冰崩-泥石流动力过程;将模拟泥石流所得到的堰塞坝体嵌入地形中,采用ArcGIS计算堰塞湖范围及体积,通过曼宁模型模拟堰塞湖溃决洪水动力过程。采用分段模拟法再现冰崩-泥石流-堵江-堰塞湖-溃坝的完整动力过程,对泥石流运动过程中的流速、流深,坝体高度,溃决洪水的流深、流速等参数进行定量化研究,为色东普流域的防灾减工作提供有效支撑。     

INTELLIGENT ALGORITHMS FOR SOLVING CNOP AND THEIR APPLICATIONS IN ENSO PREDICTABILITY AND TROPICAL CYCLONE ADAPTIVE OBSERVATIONS

Some intelligent algorithms (IAs) proposed by us, including swarm IAs and single individual IAs, have been applied to the Zebiak-Cane (ZC) model to solve conditional nonlinear optimal perturbation (CNOP) for studying El Ni?o – Southern Oscillation (ENSO) predictability. Compared to the adjoint-based method (the ADJ-method), which is referred to as a benchmark, these IAs can achieve approximate CNOP results in terms of magnitudes and patterns. Using IAs to solve CNOP can avoid the use of an adjoint model and widen the application of CNOP in numerical climate and weather modeling. Of the proposed swarm IAs, PCA-based particle swarm optimization (PPSO) obtains CNOPs with the best patterns and the best stability. Of the proposed single individual IAs, continuous tabu search algorithm with sine maps and staged strategy (CTS-SS) has the highest efficiency. In this paper, we compare the validity, stability and efficiency of parallel PPSO and CTS-SS using these two IAs to solve CNOP in the ZC model for studying ENSO predictability. The experimental results show that CTS-SS outperforms parallel PPSO except with respect to stability. At the same time, we are also concerned with whether these two IAs can effectively solve CNOP when applied to more complicated models. Taking the sensitive areas identification of tropical cyclone adaptive observations as an example and using the fifth-generation mesoscale model (MM5), we design some experiments. The experimental results demonstrate that each of these two IAs can effectively solve CNOP and that parallel PPSO has a higher efficiency than CTS-SS. We also provide some suggestions on how to choose a suitable IA to solve CNOP for different models.