Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

Groundwater mixing between different aquifer types in a complex structural setting discerned by elemental and stable isotope geochemistry

Karst aquifers are well known for their intricate stratigraphy and geologic structures, which make groundwater characterization challenging because flowpaths and recharge sources are complex and difficult to evaluate. Geochemical data, collected from ten closely spaced production wells constructed in two karst aquifers (Bangor Limestone (Mb) and Tuscumbia Limestone/Fort Payne Chert (Mftp)) in Trussville, north‐central Alabama, illustrate two distinctive groundwater end‐members: (1) higher major ion, dissolved inorganic carbon, conductivity, alkalinity concentrations, heavier δ¹³C ratios (max: −10.2 ± 0.2‰ Vienna Pee Dee Belemnite (PDB)) and lower residence times (mean: 19.5 ± 2 years, n = 2) of groundwater in the Mb aquifer and (2) lower constituent concentrations, lighter δ¹³C ratios (min: −13.4 ± 0.2‰ PDB) and longer residence times of groundwater (mean: 23.6 ± 2 years, n = 4) in the Mftp aquifer. Summer and fall data and the binary mixing model show aquifer inter‐flow mixing along solution fractures and confirms the distinctive groundwater geochemistry of the two aquifers. Lowering of static water levels over the summer (drawdown from 2 to 5.2 m) leads to more reducing groundwater conditions (lower Eh values) and slightly enriched δ¹⁸O and δD ratios during the fall [δ¹⁸O: −4.8 ± 0.1 to −5.4 ± 0.1‰ Vienna Standard Mean Oceanic Water (VSMOW), n = 9; δD: −25.4 ± 1 to −27.4 ± 1‰ VSMOW, n = 9] when compared with summer season samples (δ¹⁸O: −5.1 ± 0.1 to −5.7 ± 0.1‰ VSMOW, n = 11; δD: −25.0 ± 1 to −30.6 ± 1‰ VSMOW, n = 11). GIS analyses confirm the localized origin of recharge to the investigated aquifers. The combination of GIS, field parameters and geochemistry analyses can be successfully used to identify recharge sources, evaluate groundwater flow and transport pathways and to improve understanding of how groundwater withdrawals impact the sustainability and susceptibility to contamination of karst aquifers. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydraulic modeling of flood pulses in the Middle Bussento Karst System (MBSKS), UNESCO Cilento Global Geopark,southern Italy

Karst systems provide water for domestic and industrial uses and for generating hydropower, but they can also create fluvial hazards, such as upstream back‐flooding and downstream karst flash‐flood events. However, these hazards are difficult to foresee due to the complex recharge‐discharge processes as well as the lack of information on the inside of the system, which has often not been completely surveyed by speleologists or explored by boreholes. To overcome these difficulties, hydro‐chemical data from the monitoring system in the Middle Bussento Karst System (MBSKS), one of the first Experimental Karst Systems in southern Italy, were recorded and previously discussed. Based on shared background in flood karst hydraulic modeling, this paper describes the conceptual premises and rationale of a general‐purpose hydraulic model that is suitable both for the MBSKS and for other Mediterranean, multi‐recharge, mature, conduit‐dominated karst systems. To test the reliability of the model, simulations of time–space behavior and response are performed using natural and artificial flood pulses “as tracers”, considering a “pulse” as a significant variation in water quantity and/or quality. The results of the model explain the interactions between allogenic, autogenic, and anthropogenic recharges from differentiated sources and phreatic conduit systems. These results also clarify the overall response of karst springs at typical time scales of flood pulses. Table acronym name     

内蒙古自治区重要矿种成矿规律综述

2007~2013年开展的内蒙古矿产资源潜力评价项目,对铁、铜、铅锌、金、银、钨、钼、铬、镍等重要矿种进行了区域成矿规律总结,为矿产资源潜力评价提供了基础资料。文章即是对该项工作部分成果的概括总结。主要进展包括:在全国Ⅲ级成矿区带划分的基础上,首次进行了覆盖全自治区的Ⅳ级成矿亚带的划分,共划分出34个Ⅳ级成矿亚带;对内蒙古铁、铜等11个重要矿种的主要矿床类型及成矿特征进行了概述,对其时空分布规律做了归纳,认为全区70%以上的矿床数量和资源储量均集中在Ⅲ_5、Ⅲ_6、Ⅲ_8、Ⅲ_10和Ⅲ_11五个成矿区带。此外,不同矿种甚至同一矿种,由于成矿地质背景的差异,在不同的三级区带中的分布也不一样。主要成矿期为元古宙和中生代,次为太古宙和晚古生代,不同矿种的重要成矿期也不完全相同。从区域演化的角度探讨了构造与成矿的关系,认为不同的构造演化阶段形成不同的矿床类型和不同矿种的矿床,其中,古大陆边缘裂谷带以白云鄂博式铁_稀土元素矿床为代表,而大兴安岭岩浆岩带则以产出与燕山期中酸性火山侵入杂岩有关的多金属矿产为特点。     

塔里木盆地西南凹陷古新统巨厚海相石膏岩沉积特征及环境意义

塔里木盆地(简称塔,下同)西南凹陷古新世阿尔塔什组发育巨厚层海相石膏岩,夹薄层泥岩、粉砂岩及灰岩,是塔西南凹陷断续海侵环境下多期次蒸发沉积的产物。野外调查显示,该层海相石膏岩出现于皮拉里、阿尔塔什、麻扎塔格及大山口地区的阿尔塔什组露头剖面。石膏岩在凹陷内分布广泛,在西昆仑山前、南天山山前及麦盖提斜坡带均有发育且沉积厚度比较稳定。石膏岩中主要盐类矿物为石膏、硬石膏。扫描电镜分析发现,石膏岩中尚含石盐、钙芒硝及含钾镁的硫酸盐等;石膏岩样品中石盐、石膏、硬石膏等多呈细晶或自形微晶,推测阿尔塔什组沉积期,古盐湖曾出现过富钾卤水;通过对皮拉里剖面石膏岩样品进行地球化学分析,揭示该地区古新世古盐湖演化过程中出现2个相对富钾峰值。在古盐湖演化过程中,由于多期次特提斯海水的侵入,凹陷内部阿尔塔时期发育了碎屑岩—巨厚层石膏岩—碎屑岩的沉积韵律,古盐湖卤水表现为淡—咸—淡的变化规律。伴随着阿尔塔时期4次大规模的海侵,石膏岩沉积从西昆仑山前扩展到麦盖提斜坡地带,海侵范围也逐渐扩大。在阿尔塔什组顶部发育中厚层灰岩,显示在阿尔塔什组沉积晚期,塔西南凹陷沉积环境从●湖相逐渐向浅海相环境演变。     

前陆盆地钾盐矿床成因及模式——以西班牙北部埃布罗盆地为例

前陆盆地蕴藏有巨量的钾盐资源,而对于其内的钾盐矿床成因和模式还欠缺系统的总结。西班牙埃布罗盆地是由于伊比利亚和欧亚块体碰撞而形成的前陆盆地。始新世晚期(约36 Ma)海水完全从盆地退出后,在极端干旱气候作用下,由于碰撞造山导致盆地的封闭作用,在南比利牛斯前陆盆地系统的前渊带(即埃布罗盆地北部)形成了典型的厚层含钾石盐_光卤石的正常海相蒸发岩序列。后期受到构造挤压作用,钾盐地层以盐底劈的形式出露在背斜核部。埃布罗盆地钾盐成因是构造、气候和物源三者耦合作用的结果,与中国库车前陆盆地有很大的相似性。据此,作者建议可重点关注盆地南北盐丘地带苏维依组蒸发岩以及卤水的迁移方向。     

新疆东天山地区土屋和延东铜矿床斑岩叠加改造成矿作用

土屋和延东铜矿床位于东天山大南湖-头苏泉岛弧带南部,是中亚成矿带的重要组成部分。文章根据脉次穿插关系、蚀变矿物组合及矿物共生关系,将土屋和延东铜矿床均划分为斑岩成矿期、叠加改造期和表生期3个期次。土屋铜矿床的铜矿化形成于斑岩成矿期和叠加改造期,而延东铜矿床的铜矿化主要形成于叠加改造期;土屋和延东铜矿床伴生的钼矿化主要形成于叠加改造期。因此,笔者认为前人获得的辉钼矿Re-Os年龄(326.2~322.7 Ma)代表叠加改造期的成矿年龄,该期矿化与石英钠长斑岩((323.6±2.5)Ma)的侵入相关,而斑岩成矿期的矿化与斜长花岗斑岩(339~332 Ma)相关,成矿年龄为341.2~333.9 Ma。叠加改造期的存在,使得斑岩成矿期的蚀变分带可能受到了叠加和破坏。