Groundwater vulnerability mapping: A GIS and fuzzy rule based integrated tool

Contamination of groundwater has become a major concern in recent years. Since testing of water quality of all domestic and irrigation wells within large watersheds is not economically feasible, one frequently used monitoring strategy is to develop contamination potential maps of groundwater, and then prioritize those wells located in the potentially highly contaminated areas for testing of contaminants. However, generation of contamination potential maps based on groundwater sensitivity and vulnerability is not an easy task due inherent uncertainty. Therefore, the overall goal of this research is to improve the methodology for the generation of contamination potential maps by using detailed landuse/pesticide and soil structure information in conjunction with selected parameters from the DRASTIC model. The specific objectives of this study are (i) to incorporate GIS, GPS, remote sensing and the fuzzy rule-based model to generate groundwater sensitivity maps, and (ii) compare the results of our new methodologies with the modified DRASTIC Index (DI) and field water quality data. In this study, three different models were developed (viz. DI_fuzz, VI_fuzz and VI_{fuzz_ped}) and were compared to the DI. Once the preliminary fuzzy logic-based (DI_fuzz) was generated using selected parameters from DI, the methodology was further refined through VI_fuzz and VI_{fuzz_ped} models that incorporated landuse/pesticide application and soil structure information, respectively. This study was conducted in Woodruff County of the Mississippi Delta region of Arkansas. Water quality data for 55 wells were used to evaluate the contamination potential maps. The sensitivity map generated by VI_{fuzz_ped} with soil structure showed significantly better coincidence results when compared with the field data.     

Climatic response of thick leaf spruce (Picea crassifolia) tree-ring width at different elevations over Qilian Mountains, northwestern China

Tree-ring cores of thick leaf spruce (Picea crassifolia) taken from four sites at different elevations, in the middle of the Qilian Mountains, in the arid and semi-arid region of northwestern China, were used to develop four tree-ring width chronologies using standard dendrochronological methods. Results indicate that with increasing altitude the chronologies’ year-to-year variations decreased. Hence, the sensitivity of the tree-ring chronologies to climate decreases with altitude. Further analysis showed that the significant limiting factor on tree growth is spring precipitation. Measurements of stomata density and leaf dry weight suggest the species’ ecological adaptation strategy changes with elevation. At high elevation the metabolic rate of thick leaf spruce decreases, thus showing the effect of the climate.     

GIS-based detachment susceptibility analyses of a cut slope in limestone, Ankara—Turkey

Due to the rapidly growing population of the city of Ankara (Turkey) and increased traffic congestion, it has become necessary to widen the Ankara-Eskişehir (E-90) highway connecting the newly built areas west of the city to the city center. During widening, several cut slopes were formed along the highway route. As a result, some instability problems (small-sized rock falls/sliding, sloughing, raveling) produced detachment zones along a cut slope in highly jointed, folded and sheared limestone, causing local degradation of the cut slope. Identification of the areas that are likely to detach from the cut slope in the future is considered to be very important for the application of remedial measures. For this purpose, the relationships between the existing detachment zones and various parameters (e.g., point load strength index, weathering, block size, daylighting, shear zone) were investigated using GIS-based statistical detachment susceptibility analyses in order to predict the further aerial extension of the detachment zones with time. During the overlay analyses, statistical index and weighting factor methods were used. The outcomes of the analyses were compared and evaluated with the field observations to check the reliability of the methods and to assess the detachment zones that may develop in the future. The detachment susceptibility map without the block-size layer gives the best result and indicates some risky zones where detachments are likely to occur in the future. Recommendations on remedial measures of the cut slope should consider these risky zones.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Sensitivity of a hydraulic model to channel erosion uncertainty during extreme flooding

Recent research into flood modelling has primarily concentrated on the simulation of inundation flow without considering the influences of channel morphology. River channels are often represented by a simplified geometry that is implicitly assumed to remain unchanged during flood simulations. However, field evidence demonstrates that significant morphological changes can occur during floods to mobilize the boundary sediments. Despite this, the effect of channel morphology on model results has been largely unexplored. To address this issue, the impact of channel cross‐section geometry and channel long‐profile variability on flood dynamics is examined using an ensemble of a 1D–2D hydraulic model (LISFLOOD‐FP) of the ~1 : 2000 year recurrence interval floods in Cockermouth, UK, within an uncertainty framework. A series of simulated scenarios of channel erosional changes were constructed on the basis of a simple velocity‐based model of critical entrainment. A Monte‐Carlo simulation framework was used to quantify the effects of this channel morphology together with variations in the channel and floodplain roughness coefficients, grain size characteristics and critical shear stress on measures of flood inundation. The results showed that the bed elevation modifications generated by the simplistic equations reflected an approximation of the observed patterns of spatial erosion that enveloped observed erosion depths. The effect of uncertainty on channel long‐profile variability only affected the local flood dynamics and did not significantly affect the friction sensitivity and flood inundation mapping. The results imply that hydraulic models generally do not need to account for within event morphodynamic changes of the type and magnitude of event modelled, as these have a negligible impact that is smaller than other uncertainties, e.g. boundary conditions. Instead, morphodynamic change needs to happen over a series of events to become large enough to change the hydrodynamics of floods in supply limited gravel‐bed rivers such as the one used in this research. Copyright © 2014 John Wiley & Sons, Ltd.     

Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

华北陆块晋冀Al- Fe- Au- Pb- Zn- Ag- Cu- 煤成矿带主要地质成矿特征及潜力分析

华北陆块晋冀成矿带是我国重要的金属和能源矿产资源基地,近年来新的找矿进展和预测成果不断涌现,使之成为我国重点找矿地区之一,需要对其开展系统的成矿区划和部署研究工作。本文以新的研究成果和找矿突破为基础,通过综合分析,在该地区划分新的重要成矿部署区带,厘定边界并将其命名为"华北陆块晋冀Al-Fe-Au-Pb-Zn-Ag-Cu-煤成矿带"。结合区域成矿地质背景的综合研究,在区内划分了10个成矿亚带,初步建立了Al-Fe-Au-Pb-Zn-Ag-Cu-煤成矿带成矿谱系,其成矿作用具有多期次发育和多系统控制的特点,在时间演化上表现为多旋回性、继承性、新生性和叠加性特点,主要成矿期有前寒武纪,寒武纪-晚二叠世,晚三叠世-晚白垩世和新生代,其中晚侏罗-早白垩世是内生金属矿床成矿大爆发期。该带产出的矿床类型复杂多样,与岩浆有关成矿作用有岩浆岩型、斑岩型、矽卡岩型、海(陆)相火山岩型、岩浆热液型和石英脉型,变质成矿作用有沉积变质型、绿岩建造型,层控成因矿床有碳酸盐岩中热液型,沉积成矿作用有海(陆)相沉积型、砂矿型、风化壳型矿床等。结合区域找矿进展和潜力评价最新成果,认为本区铝土矿、铁、金、锰、钼、铜、铅锌、煤矿潜力巨大,为下一步勘查部署的主攻矿种,同时在本区划分了26个重点远景区,对研究区下一步矿产勘查部署工作有一定指导意义。