Hydrologic and water isotope characterization of a regulated Canadian Shield river basin

Stable water isotope surveys have increasingly been integrated into river basins studies, but fewer have used them to evaluate impact of hydropower regulation. This study applies hydrologic and water isotope survey approaches to a Canadian Shield river basin with both regulated and natural flows. Historical streamflow records were used to evaluate the influence of three hydroelectric reservoirs and unregulated portions of the basin on downstream flows and changes in water level management implemented after an extreme flood year (1979). In 2013, water isotope surveys of surface and source waters (e.g., rainfall, groundwater, snowmelt) were conducted to examine spatial and temporal variation in contributions to river flow. Seasonal changes in relative groundwater contribution were assessed using a water‐isotope mass balance approach. Within the basin, two regulated reservoirs exhibited inverted hydrographs with augmented winter flows, whereas a third exhibited a hydrograph dominated by spring snowmelt. In 2013, spatial variation in rain‐on‐snow and air temperatures resulted in a critical lag in snowmelt initiation in the southern and northern portions of the basin resulting in a dispersed, double peak spring hydrograph, contrasting with 1979 when a combination of rain‐on‐snow and coincident snowmelt led to the highest flood on record. Although eastern basin reservoirs become seasonally enriched in δ¹⁸O and δ²H values, unregulated western basin flows remain less variable due to groundwater driven baseflow with increasing influence downstream. Combined analysis of historical streamflow (e.g., flood of 1979, drought of 2010) and the 2013 water isotope surveys illustrate extreme meteorological conditions that current management activities are unable to prevent. In this study, the influence of evaporative fractionation on large surface water reservoirs provides important evidence of streamflow partitioning, illustrating the value of stable water isotope tracers for study of larger catchments.     

A multidirectional semi‐analytical method for analysis of laterally loaded pile groups in multi‐layered elastic strata

A semi‐analytical method for calculating the response of single piles and pile groups subjected to lateral loading is developed in this paper. Displacements anywhere in the soil domain are tied to the displacements of the piles through decay functions. The principle of virtual work and the calculus of variations are used to derive the governing differential equations that describe the response of the piles and soil. The eigenvalue method and the finite difference technique are used to solve the system of coupled differential equations for the piles and soil, respectively. The proposed method takes into account the soil surface displacement along and perpendicular to the loading direction and produces displacement fields that are very close to those produced by the finite element method but at lower computational effort. Compared with the previous method that considered only the soil displacement along the loading direction, accounting for the multi‐directional soil displacement field produces responses for the piles and soil that are closer to those approximated by the finite element method. Copyright © 2016 John Wiley & Sons, Ltd.     

Syn‐metamorphic folding in the Tauern Window,Austria dated by Th–Pb ages from individual allanite porphyroblasts

High‐precision ²³²Th–²⁰⁸Pb dates have been obtained from allanite porphyroblasts that show unambiguous microstructural relationships to fabrics in a major syn‐metamorphic fold in the SE Tauern Window, Austria. Three porphyroblasts were analysed from a single garnet mica schist from the Peripheral Schieferhülle in the core of the Ankogel Synform, one of a series of folds which developed shortly before the thermal peak of Alpine epidote–amphibolite facies metamorphism: allanite grain 1 provided two analyses with a combined age of 27.7 ± 0.7 Ma; grain 2, which was slightly bent and fractured during crenulation, provided two analyses with a combined age of 27.7 ± 0.4 Ma; a single analysis from grain 3, which overgrew an already crenulated fabric, gave an age of 28.0 ± 1.4 Ma. The five ²³²Th–²⁰⁸Pb ages agree within error and define an isochron with an age of 27.71 ± 0.36 Ma (95% confidence level; MSWD = 0.46). The results imply that the crenulation event was in progress in a short interval (<1 Ma) c. 28 Ma, and that the Ankogel Synform was forming at this time. The thermal peak of regional metamorphism in the SE Tauern Window was probably attained shortly after 28 Ma, only c. 5 Ma after eclogite facies metamorphism in the central Tauern Window. Metasediment may contain allanite porphyroblasts with clear‐cut microstructural relationships to fabric development and metamorphic crystallization; for such rocks, ²³²Th–²⁰⁸Pb dating on microsamples offers a powerful geochronological tool.     

Metals discharged during different flow conditions from a mixed agricultural‐forest catchment (NW Spain)

Metal loads were determined from water samples collected under different streamflow conditions (baseflow and storm events) in a rural catchment (NW Spain) during 4 years. A study at annual, seasonal and storm‐event scales was carried out. In all analysed scales, the export order was Fe > Al > Mn > Zn > Cu. A high inter‐annual, seasonal and storm‐event scale variability of metal load was observed. The total metal loads in stream were higher during baseflow conditions than during storm events, which only represented 4% of the duration of the study period and 25% of streamflow. During storm events, both Al and Fe loads accounted 45% of the total load of the study period, whereas Mn, Cu and Zn loads represented 42%, 33% and 24%, respectively. This highlights the role of high flows on metal export. Only four big events exported around 30% of load of each metal transported in events. At all time scales, a prevalence of export of particulate metals over dissolved metals was observed, more pronounced for Al, Fe and Mn than for Cu and Zn. The export of metals in the Corbeira catchment is influenced by runoff and, to a lesser extent, by the rainfall amount. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical and analytical observations on long and infinite slopes

Interest in the mechanics of landslides has led to renewed evaluation of the infinite slope equations, and the need for a more general framework for estimating the factor of safety of long and infinite slopes involving non‐homogeneous soil profiles. The paper describes finite element methods that demonstrate the potential for predicting failure in long slope profiles where the critical mechanism is not necessarily at the base of the soil layer. The influence of slope angle is also examined in long slopes, leading to some counter‐intuitive conclusions about the impact of slope steepness on the factor of safety. Copyright © 2010 John Wiley & Sons, Ltd.     

Diatom-inferred depth models in 8 Canadian boreal lakes: inferred changes in the benthic:planktonic depth boundary and implications for assessment of past droughts

Assessment of past drought in boreal regions, a region predicted to be strongly affected by climate warming, can provide insights into future availability of water. However, limited instrumental data and paleoclimatic data are available for this assessment. To address this lack of data in the boreal region of northwest Ontario, a regional study of lakes in the Winnipeg River Drainage Basin was initiated. Diatom-inferred (D-I) depth models were developed based on surface samples collected along a depth gradient within 8 small boreal lakes. Weighted-averaging and modern analog approaches provided robust within-lake depth models for each of the study lakes, with bootstrapped r² values ranging from 0.90 to 0.98, and root-mean-squared-errors of prediction (RMSEP) between 1.1 and 2.5 m. Large differences in the estimated depth optima for three representative, but common diatom species across our 8 study lakes suggested that within-lake calibration datasets are more appropriate for inferring past drought based in depth models than a regional multi-lake calibration dataset, and that light and related variables are controlling factors governing the maximum depth of benthic taxa. A down-core application of the D-I depth models on a near-shore core from Meekin Lake, retrieved near the present-day ecotone between the benthic and planktonic diatom assemblages indicated highly similar trends in inferred depth (r = 0.96). The models have significant correlations with other metrics of changes in depth including diatom species richness (r = 0.74–0.78) and evenness (r = 0.76–0.8), thereby allowing a check on the strength and direction of the depth inferences down-core. Near-shore cores located near the benthic:planktonic transition is a sensitive region that can provide estimates of past droughts in lakes where such inferences have been difficult to estimate.     

Complex hydraulic interactions between compartmentalized carbonate aquifers and heterogeneous siliciclastic successions: a case study in southern Italy

The sedimentary sequences containing lithologic units with low permeability represent hydrogeologic systems, which, as of now, have been little studied despite their diffusion worldwide. A hydrogeologic study, aimed to assess the main factors controlling the groundwater flow dynamics in such systems and their hydraulic interactions with nearby carbonate aquifers, has been carried out in Longano (Isernia, Southern Italy). The analysis of the hydraulic heads, combined with the regimes of the springs and the electric conductivity of the groundwater, mainly reflect vertical and lateral heterogeneities of the media in terms of hydraulic properties. In particular, the flow system is controlled by lateral heterogeneities, which characterize a surficial horizon made up of clayey colluviums and talus deposits, separated from the deeper saturated, fissured bedrock. One‐to‐ten relationships in hydraulic heads, monitored in piezometers crossing the fissured media, further uphold the crucial role played by the lateral contrasts of permeability in controlling the flow dynamics. On the whole, significant interactions with the nearby carbonate hydrostructure take place. Nevertheless, the heterogeneities of the siliciclastic succession and surficial horizon, coupled with the compartmentalization of the carbonate system, lead to a complex hydrogeological scenario. In a wider perspective, this study gives information of utmost importance in order to improve the implementation of mathematical models and configuration of tapping works within these heterogeneous and complex settings. Copyright © 2015 John Wiley & Sons, Ltd.     

Scaling procedures for heterogeneous unconfined aquifers

Complex flows in heterogeneous confined and unconfined aquifers is a phenomenon that continues to present difficulties in flow mapping and modelling in the field, laboratory, and through numerical simulations. It is often the case with complicated phenomena that transformative scaling and reduction of the problem through symmetry is of great efficacy in the formation of predictive models in both the laboratory and computational settings. A detailed a study of the application of a broad class of Lie scaling transformations on a set of equations representing the groundwater flows in heterogeneous confined and unconfined aquifers has produced a set of scaling relationships between the spatial variables, hydrologic variables, and parameters. The set of scaling transformations preserve the structure of the equations in the sense that the scaling transformations leave the initial‐boundary value system representing the invariant groundwater flows. This theoretical approach elucidates not only the scaling relationships but also the properties that hydrologic variables and parameters must satisfy in order for calling to be possible. Validation of the theory developed is carried out through a series of four numerical simulations using the USGS modflow ‐2005 software package. The results of these experiments demonstrate that the derived scaling transformations can effectively form predictive models of large‐scale phenomena at small scales with negligible error in many cases. Comments on the limitations of the approach and directions for future research are made in the closing sections. Copyright © 2016 John Wiley & Sons, Ltd.     

Modelling sediment deposition and phosphorus retention in a river floodplain

Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitute an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development of a suitable prediction model is justified, and the present work is devoted to the validation and application of a predictive soluble reactive phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed two‐dimensional shallow‐water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated values of the three experiments were significant (r² ≥ 0.62; p ≤ 0.05), and a Nash–Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during two real flooding events within the same river reach. The comparison between calculated and measured sediment depositions showed a significant linear regression (p ≤ 0.05; r² = 0.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyse floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river and controls the river matter contribution to the floodplain. Copyright © 2014 John Wiley & Sons, Ltd.