Does upward seepage of river water and storm water runoff determine water quality of urban drainage systems in lowland areas? A case study for the Rhine–Meuse delta

The water quality of urban drainage ditches in lowlands in the Rhine‐Meuse delta was analysed with principal component analysis (PCA) during a dry period and a rain storm, and related to the seepage of polluted river water and effective impervious area (EIA). This was done in order to test the hypothesis that seepage of river water and storm water runoff from impervious areas strongly determine the water quality of urban drainage systems along large lowland rivers. Our analysis revealed that upward seepage of groundwater originating from rivers Rhine and Meuse was positively correlated with nitrate, potassium, sodium and chloride and negatively correlated with alkalinity, calcium, magnesium and iron. EIA was correlated with very few environmental variables (i.e. phosphate, pH and iron in the dry period and iron during the rain storm). Nickel and zinc concentrations generally exceeded the maximum allowable concentrations (MAC), while lead and phosphorus concentrations were just above the nutrient standards and MAC in a few locations during the rain storm. To optimize water quality in urban water systems, attention should be paid to all sources of pollution and not only to EIA. The impact of local groundwater seepage originating from large rivers in lowlands on the chemistry of urban water systems is often underestimated and should be taken into account when assessing water quality and improving water quality status. Copyright © 2009 John Wiley & Sons, Ltd.     

A continuum damage concrete model for earthquake analysis of concrete gravity dam–reservoir systems

In this study, the earthquake damage response of the concrete gravity dams is investigated with considering the effects of dam–reservoir interaction. A continuum damage model which is a second-order tensor and includes the strain softening behavior is selected for the concrete material. The mesh-dependent hardening technique is adopted such that the fracture energy dissipated is not affected by the finite element mesh size. The dynamic equilibrium equations of motion are solved by using the improved form of the HHT-α time integration algorithm. Two dimensional seismic analysis of Koyna gravity dam is performed by using the 1967 Koyna earthquake records. The effects of damage on the earthquake response of concrete gravity dams are discussed. Comparison of the Westergaard and Lagrangian dam–reservoir interaction solutions is made. The effects of viscous damping ratio on the damage response of the dam are also studied.     

Determination of groundwater discharge rates and water residence time of groundwater‐fed lakes by stable isotopes of water (18O, 2H) and radon (222Rn) mass balances

Lacustrine groundwater discharge (LGD) and the related water residence time are crucial parameters for quantifying lake matter budgets and assessing its vulnerability to contaminant input. Our approach utilizes the stable isotopes of water (δ¹⁸O, δ²H) and the radioisotope radon (²²²Rn) for determining long‐term average and short‐term snapshots in LGD. We conducted isotope balances for the 0.5‐km² Lake Ammelshainer See (Germany) based on measurements of lake isotope inventories and groundwater composition accompanied by good quality and comprehensive long‐term meteorological and isotopic data (precipitation) from nearby monitoring stations. The results from the steady‐state annual isotope balances that rely on only two sampling campaigns are consistent for both δ¹⁸O and δ²H and suggested an overall long‐term average LGD rate that was used to infer the water residence time of the lake. These findings were supported by the good agreement of the simulated LGD‐driven annual cycles of δ¹⁸O and δ²H lake inventories with the observed lake isotope inventories. However, radon mass balances revealed lower values that might be the result of seasonal LGD variability. For obtaining further insights into possible seasonal variability of groundwater–lake interaction, stable water isotope and radon mass balances could be conducted more frequently (e.g., monthly) in order to use the derived groundwater discharge rates as input for time‐variant isotope balances.     

Earthquake response of concrete arch dams: a plastic–damage approach

There are several alternatives to evaluate seismic damage‐cracking behavior of concrete arch dams, among which damage theory is the most popular. A more recent option introduced for this purpose is plastic–damage (PD) approach. In this study, a special finite element program coded in 3‐D space is developed on the basis of a well‐established PD model successfully applied to gravity dams in 2‐D plane stress state. The model originally proposed by Lee and Fenves in 1998 relies on isotropic damaged elasticity in combination with isotropic tensile and compressive plasticity to capture inelastic behaviors of concrete in cyclic or dynamic loadings. The present implementation is based on the rate‐dependent version of the model, including large crack opening/closing possibilities. Moreover, with utilizing the Hilber–Hughes–Taylor time integration scheme, an incremental–iterative solution strategy is detailed for the coupled dam–reservoir equations while the damage–dependent damping stress is included. The program is initially validated, and then, it is employed for the main analyses of the Koyna gravity dam in a 3‐D modeling as well as a typical concrete arch dam. The former is a major verification for the further examination on the arch dam. The application of the PD model to an arch dam is more challenging because the governing stress condition is multiaxial, causing shear damage to become more important than uniaxial states dominated in gravity dams. In fact, the softening and strength loss in compression for the damaged regions under multiaxial cyclic loadings affect its seismic safety. Copyright © 2013 John Wiley & Sons, Ltd.     

Improved spatial delineation of streambed properties and water fluxes using distributed temperature sensing

A new method was developed for analysing and delineating streambed water fluxes, flow conditions and hydraulic properties using coiled fibre‐optic distributed temperature sensing or closely spaced discrete temperature sensors. This method allows for a thorough treatment of the spatial information embedded in temperature data by creating a matrix visualization of all possible sensor pairs. Application of the method to a 5‐day field dataset reveals the complexity of shallow streambed thermal regimes. To understand how velocity estimates are affected by violations of assumptions of one‐dimensional, saturated, homogeneous flow and to aid in the interpretation of field observations, the method was also applied to temperature data generated by numerical models of common field conditions: horizontal layering, presence of lateral flow and variable streambed saturation. The results show that each condition creates a distinct signature visible in the triangular matrices. The matrices are used to perform a comparison of the behaviour of one‐dimensional analytical heat‐tracing models. The results show that the amplitude ratio‐based method of velocity calculation leads to the most reliable estimates. The minimum sensor spacing required to obtain reliable velocity estimates with discrete sensors is also investigated using field data. The developed method will aid future heat‐tracing studies by providing a technique for visualizing and comparing results from fibre‐optic distributed temperature sensing installations and testing the robustness of analytical heat‐tracing models. Copyright © 2016 John Wiley & Sons, Ltd.     

Temporal variation of stable isotopes in a precipitation–groundwater system: implications for determining the mechanism of groundwater recharge in high mountain–hills of the Loess Plateau,China

The groundwater in shallow loess aquifers in high mountain–hills in the western Loess Plateau in China is almost the sole water resource for local residents. However, the question about how the loess groundwater naturally circulates in these high mountain–hills, characterized by low precipitation and high potential evaporation, remains unclear. The objectives of this study are to evaluate the application of hydrogen and oxygen isotopes to (1) examine temporal variations of the isotopic composition of precipitation and shallow groundwater and (2) uncover the mechanism of groundwater recharge in high mountain–hills. Results from 2 years of monitoring data show a difference in the stable isotopes for groundwater and local precipitation between the winter and summer periods. Similar to precipitation, stable isotopes in groundwater are observed to be depleted in winter and enriched in summer, particularly in oxygen isotope. A prominent characteristic is that H and O isotopes of groundwater show a very clear response to strong precipitation in the rainy season in 2013. The results highlight that local precipitation is the likely recharge source for groundwater in shallow loess aquifers. Annual recharge from local precipitation maintains the groundwater resource in the shallower loess aquifer. The mechanisms governing shallow loess groundwater recharge in high mountain–hills were evaluated. In addition to possible vertical slow percolation of soil water through the unsaturated zone, rapid groundwater recharge mechanisms have been identified as temporal preferential infiltration through sinkholes, slip surface or landslide surface and through the interface of loess layer and palaeo‐soils. Most groundwater can be recharged after a heavy rainy season. Copyright © 2015 John Wiley & Sons, Ltd.     

Non‐linear behavior of mass concrete in three‐dimensional problems using a smeared crack approach

A smeared crack approach has been proposed to model the static and dynamic behavior of mass concrete in three‐dimensional space. The proposed model simulates the tensile fracture on the mass concrete and contains pre‐softening behavior, softening initiation, fracture energy conservation and strain rate effects under dynamic loads. The validity of the proposed model has been checked using the available experimental results under static and dynamic loads. The direct and indirect displacement control algorithms have been employed under incremental increasing static loads. It was found that the proposed model gives excellent results and crack profiles when compared with the available data under static loads. The Koyna Dam in India has been used to verify the dynamic behavior of the proposed model. It was found that the resulting crack profiles were in good agreement with the available experimental results. Finally, the Morrow Point Dam was analyzed, including the dam–reservoir interaction effects, to consider its non‐linear seismic behavior. It was found that the resulting crack profiles were in good agreement with the contour of maximum principal stresses and no numerical instability occurred during the analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

Investigation of the thermal regime and subsurface properties of a tidally affected,variably saturated streambed

Temperature and moisture content in the variably saturated subsurface are two of the most important physical parameters that govern a wide variety of geochemical and ecological processes. An understanding of thermal and hydraulic processes and properties of transient vadose zones is therefore fundamental in the evaluation of such processes. Here, an investigation of the thermal regime and subsurface properties of a tidally affected, variably saturated streambed is presented. Field and laboratory measurements, as well as a forward numerical model, are jointly employed in the investigation. Temperature, soil moisture, surface level, and water level data were recorded in a transect perpendicular to a tidally driven stream. Frequency‐domain analysis of the subsurface temperature measurements revealed the rapid decay of the tidal temperature driver within the top ～30 cm of sediment. Several techniques were used to evaluate subsurface thermal and hydraulic properties, including thermal conductivity and the soil water retention curve. These properties were used to constrain a forward numerical model that included coupled treatment of relevant variable saturation thermal and hydraulic physics. Even though the investigated vadose zone is intermittent and relatively shallow ( 20 cm), the results illustrate how error can be introduced into heat‐transport calculations if unsaturated conditions are not taken into account.