9: Processes of tectonism, magmatism and mineralization: Lessons from Europe

Metallogenic provinces in Europe range in age from the Archaean to the Neogene. Deposit types include porphyry copper and epithermal Cu–Au, volcanic-hosted massive sulphide (VMS), orogenic gold, Fe-oxide–Cu–Au, anorthosite Fe–Ti-oxide and sediment-hosted base-metal deposits. Most of them formed during short-lived magmatic events in a wide range of tectonic settings; many can be related to specific tectonic processes such as subduction, hinge retreat, accretion of island arcs, continental collision, lithosphere delamination or slab tear. In contrast, most sediment-hosted deposits in Europe evolved in extensional, continental settings over significant periods of time. In Europe, as elsewhere, ore formation is an integral part of the geodynamic evolution of the Earth's crust and mantle. Many tectonic settings create conditions conducive to the generation of water-rich magma, but the generation of ore deposits appears to be restricted to locations and short periods of change in temperature and stress, imposed by transitory plate motions. Crustal influence is evident in the strong structural controls on the location and morphology of many ore deposits in Europe. Crustal-scale fault–fracture systems, many involving strike-slip elements, have provided the fabric for major plumbing systems. Rapid uplift, as in metamorphic core complexes, and hydraulic fracturing can generate or focus magmatic–hydrothermal fluid flow that may be active for time spans significantly less than a million years. Once a hydrologically stable flow is established, ore formation is strongly dependent on the steep temperature and pressure gradients experienced by the fluid, particularly within the upper crust. In Europe, significant fracture porosity deep in the crystalline basement (1%) is not only important for magmatic–hydrothermal systems, but allows brines to circulate down through sedimentary basins and then episodically upward, expelled seismically to produce sediment-hosted base-metal deposits and Kupferschiefer copper deposits. Emerging research, stimulated by GEODE, can improve the predicting power of numerical simulations of ore-forming processes and help discover the presence of orebodies beneath barren overburden.     

Test of a distributed modelling approach to predict flood flows in the karst Suoimuoi catchment in Vietnam

The major obstacles for modelling flood processes in karst areas are a lack of understanding and model representations of the distinctive features and processes associated with runoff generation and often a paucity of field data. In this study, a distributed flood-modelling approach, WetSpa, is modified and applied to simulate the hydrological features and processes in the karst Suoimuoi catchment in northwest Vietnam. With input of topography, land use and soil types in a GIS format, the model is calibrated based on 15 months of hourly meteorological and hydrological data, and is used to simulate both fast surface and conduit flows, and groundwater discharges from karst and non-karst aquifers. Considerable variability in the simulation accuracy is found among storm events and within the catchment. The simulation results show that the model is able to represent reasonably well the stormflows generated by rainfall events in the study catchment.     

Geochemical and isotopic investigation of the aquifer system in the Djerid-Nefzaoua basin, southern Tunisia

In the Djerid-Nefzaoua region, southern Tunisia, about 80% of agricultural and domestic water supply is provided by the complex terminal (CT) aquifer. However, 20% of this demand is provided by other hydraulically connected aquifers, namely the continental intercalaire (CI) and the Plio-Quaternary (PQ). Overexploitation of the CT aquifer for agricultural practices has contributed to the loss of the artesian condition and the decline of groundwater level which largely increased the downward leakage from the shallow PQ aquifer. Excess irrigation water concentrates at different rates in the irrigation channels and in the PQ aquifer itself. Then, it returns to the CT aquifer and mixes with water from the regional flow system, which contributes to the salinization of the CT groundwater. A geochemical and isotopic study had been undertaken over a 2-years period in order to investigate the origin of waters pumped from the CT aquifer with an emphasis on its hydraulic relationships with the underlying and the overlying CI and PQ aquifers. Geochemistry indicates that groundwater samples collected from different wells show an evolution of the water types from Na-Cl to Ca-SO₄-Cl. Dissolution of halite, gypsum and anhydrite-bearing rocks is the main mechanism that leads to the salinization of the groundwater. Isotopic data indicate the old origin of all groundwater in the aquifer system. Mixing and evaporation effects characterizing the CT and the PQ aquifers were identified using δ²H and δ¹⁸O relationship and confirmed by the conjunction of δ²H with chloride concentration.     

Urban stormwater thermal gene expression models for protection of sensitive receiving streams

Thermal impact of typical high‐density residential, industrial, and commercial land uses is a major concern for the health of aquatic life in urban watersheds, especially in smaller, cold, and cool‐water streams. This is the first study of its kind that provides simple easy‐to‐use equations, developed using gene expression programming (GEP) that can guide the assessment and the design of urban stormwater management systems to protect thermally sensitive receiving streams. We developed 3 GEP models using data collected during 3 years (2009–2011) from 4 urban catchments; the first GEP model predicts event mean temperature at the inlet of the pond; the second model predicts the stormwater temperature at the outlet of the pond; and the third model predicts the temperature of the stormwater after flowing through a cooling trench and before discharging to the receiving stream. The new models have high correlation coefficients of 0.90–0.94 and low prediction uncertainty of less than 4% of the median value of the predicted runoff temperatures. Sensitivity analysis shows that climatic factors have the highest influence on the thermal enrichment followed by the catchment characteristics and the key design variables of the stormwater pond and the cooling trench. The general method presented here is easily transferable to other regions of the world (but not necessarily the exact equations developed here); also through sensitivity and parametric analysis, we gained insight on the key factors and their relative importance in modelling thermal enrichment of urban stromwater runoff.     

Glacial flood pulse effects on benthic fauna in equatorial high‐Andean streams

Equatorial glacier‐fed streams present unique hydraulic patterns when compared to glacier‐fed observed in temperate regions as the main variability in discharge occurs on a daily basis. To assess how benthic fauna respond to these specific hydraulic conditions, we investigated the relationships between flow regime, hydraulic conditions (boundary Reynolds number, Re*), and macroinvertebrate communities (taxon richness and abundance) in a tropical glacier‐fed stream located in the high Ecuadorian Andes (> 4000 m). Both physical and biotic variables were measured under four discharge conditions (base‐flow and glacial flood pulses of various intensities), at 30 random points, in two sites whose hydraulic conditions were representative to those found in other streams of the study catchment. While daily glacial flood pulses significantly increased hydraulic stress in the benthic habitats (appearance of Re* > 2000), low stress areas still persisted even during extreme flood events (Re* < 500). In contrast to previous research in temperate glacier‐fed streams, taxon richness and abundance were not significantly affected by changes in hydraulic conditions induced by daily glacial flood pulses. However, we found that a few rare taxa, in particular rare ones, preferentially occurred in highly stressed hydraulic habitats. Monte‐Carlo simulations of benthic communities under glacial flood reduction scenarios predicted that taxon richness would be significantly reduced by the loss of high hydraulic stress habitats following glacier shrinking. This pioneer study on the relationship between hydraulic conditions and benthic diversity in an equatorial glacial stream evidenced unknown effects of climate change on singular yet endangered aquatic systems. Copyright © 2013 John Wiley & Sons, Ltd.     

What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Quantifying the impact of landscape changes on hydrological variables in the alpine and cold region using hydrological model and remote sensing data

Quantifying the impact of landscape on hydrological variables is essential for the sustainable development of water resources. Understanding how landscape changes influence hydrological variables will greatly enhance the understanding of hydrological processes. Important vegetation parameters are considered in this study by using remote sensing data and VIC-CAS model to analyse the impact of landscape changes on hydrology in upper reaches of the Shule River Basin (URSLB). The results show there are differences in the runoff generation of landscape both in space and time. With increasing altitude, the runoff yields increased, with approximately 79.9% of the total runoff generated in the high mountains (4200–5900 m), and mainly consumed in the mid-low mountain region. Glacier landscape produced the largest runoff yields (24.9% of the total runoff), followed by low-coverage grassland (LG; 22.5%), alpine cold desert (AL; 19.6%), mid-coverage grassland (MG; 15.6%), bare land (12.5%), high-coverage grassland (HG; 4.5%) and shrubbery (0.4%). The relative capacity of runoff generation by landscapes, from high to low, was the glaciers, AL, LG, HG, MG, shrubbery and bare land. Furthermore, changes in landscapes cause hydrological variables changes, including evapotranspiration, runoff and baseflow. The study revealed that HG, MG, and bare land have a positive impact on evapotranspiration and a negative impact on runoff and baseflow, whereas AL and LG have a positive impact on runoff and baseflow and a negative impact on evapotranspiration. In contrast, glaciers have a positive impact on runoff. After the simulation in four vegetation scenarios, we concluded that the runoff regulation ability of grassland is greater than that of bare land. The grassland landscape is essential since it reduced the flood peak and conserved the soil and water.     

Landscape controls on long‐term runoff in subhumid heterogeneous Boreal Plains catchments

We compared median runoff (R) and precipitation (P) relationships over 25 years from 20 mesoscale (50 to 5,000 km²) catchments on the Boreal Plains, Alberta, Canada, to understand controls on water sink and source dynamics in water‐limited, low‐relief northern environments. Long‐term catchment R and runoff efficiency (RP^?1) were low and varied spatially by over an order of magnitude (3 to 119 mm/year, 1 to 27%). Intercatchment differences were not associated with small variations in climate. The partitioning of P into evapotranspiration (ET) and R instead reflected the interplay between underlying glacial deposit texture, overlying soil‐vegetation land cover, and regional slope. Correlation and principal component analyses results show that peatland‐swamp wetlands were the major source areas of water. The lowest estimates of median annual catchment ET (321 to 395 mm) and greatest R (60 to 119 mm, 13 to 27% of P) were observed in low‐relief, peatland‐swamp dominated catchments, within both fine‐textured clay‐plain and coarse‐textured glacial deposits. In contrast, open‐water wetlands and deciduous‐mixedwood forest land covers acted as water sinks, and less catchment R was observed with increases in proportional coverage of these land covers. In catchments dominated by hummocky moraines, long‐term runoff was restricted to 10 mm/year, or 2% of P. This reflects the poor surface‐drainage networks and slightly greater regional slope of the fine‐textured glacial deposit, coupled with the large soil‐water and depression storage and higher actual ET of associated shallow open‐water marsh wetland and deciduous‐forest land covers. This intercatchment study enhances current conceptual frameworks for predicting water yield in the Boreal Plains based on the sink and source functions of glacial landforms and soil‐vegetation land covers. It offers the capability within this hydro‐geoclimatic region to design reclaimed catchments with desired hydrological functionality and associated tolerances to climate or land‐use changes and inform land management decisions based on effective catchment‐scale conceptual understanding.     

Ressi experimental catchment: Ecohydrological research in the Italian pre-Alps

Ressi is a small (2.4 ha) forested catchment located in the Italian pre-Alps. The site became an experimental catchment to investigate the water fluxes in the soil–plant–atmosphere continuum and the impact of vegetation on runoff generation in 2012. The elevation of the catchment ranges from 598 to 721 m a.s.l. and the climate is humid temperate. The bedrock consists of rhyolites and dacites; the soil is a Cambisol. The catchment is covered by a dense forest, dominated by beech, chestnut, maple, and hazel trees. The field set up includes measurements of the rainfall in an open area, streamflow at the outlet, soil moisture at various depths and locations, and depth to water table in six piezometers at a 5- or 10-min interval. Samples of precipitation, stream water, shallow groundwater and soil water are collected monthly for tracer analysis (stable isotopes (²H and ¹⁸O), electrical conductivity and major ions), and during selected rainfall–runoff events to determine the contribution of the various sources to runoff. Since 2017, soil and plant water samples have been collected to determine the sources of tree transpiration. Data collected in the period 2012–2016 are publicly available. Data collection is ongoing, and the data set is expected to be updated on an annual basis to include the most recent measurements.     

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.