Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

A modelling framework for evaluation of the hydrological impacts of nature‐based approaches to flood risk management,with application to in‐channel interventions across a 29‐km2 scale catchment in the United Kingdom

Nature‐based approaches to flood risk management are increasing in popularity. Evidence for the effectiveness at the catchment scale of such spatially distributed upstream measures is inconclusive. However, it also remains an open question whether, under certain conditions, the individual impacts of a collection of flood mitigation interventions could combine to produce a detrimental effect on runoff response. A modelling framework is presented for evaluation of the impacts of hillslope and in‐channel natural flood management interventions. It couples an existing semidistributed hydrological model with a new, spatially explicit, hydraulic channel network routing model. The model is applied to assess a potential flood mitigation scheme in an agricultural catchment in North Yorkshire, United Kingdom, comprising various configurations of a single variety of in‐channel feature. The hydrological model is used to generate subsurface and surface fluxes for a flood event in 2012. The network routing model is then applied to evaluate the response to the addition of up to 59 features. Additional channel and floodplain storage of approximately 70,000 m³ is seen with a reduction of around 11% in peak discharge. Although this might be sufficient to reduce flooding in moderate events, it is inadequate to prevent flooding in the double‐peaked storm of the magnitude that caused damage within the catchment in 2012. Some strategies using features specific to this catchment are suggested in order to improve the attenuation that could be achieved by applying a nature‐based approach.     

Remote estimation of terrestrial evapotranspiration by Landsat 5 TM and the SEBAL model in cold and high‐altitude regions: a case study of the upper reach of the Shule River Basin,China

Evapotranspiration (ET) is an important expenditure in water and energy balances, especially on cold and high‐altitude land surfaces. Daily ET of the upper reach of the Shule River Basin was estimated using Landsat 5 TM data and the Surface Energy Balance Algorithm for Land (SEBAL) model. Based on observations made at the Suli station, the algorithms of land surface temperature and soil heat flux in SEBAL were modified. Land surface temperature was retrieved and compared with ground truth via three methods: the radiative transfer equation method, the mono‐window algorithm, and the single‐channel method. We selected the best of these methods, mono‐window algorithm, for estimating ET. The average error of daily ET estimated by the modified SEBAL model and measured by the eddy covariance system was 16.4%, with a root‐mean‐square error of 0.52 mm d^?1. The estimated ET means were 3.09, 2.48, and 1.48 mm d^?1 on June 9 (DOY 160), June 25 (DOY 176), and July 27 (DOY 208) of the year 2010, respectively. The average estimated ET on the glacier surface of all days was more than 3 mm d^?1, a measurement that is difficult to capture in‐situ and has rarely been reported. This study will improve the understanding of water balance in cold, high‐altitude regions. Copyright © 2016 John Wiley & Sons, Ltd.     

Long‐term monitoring of a mercury contaminated estuary (Ria de Aveiro,Portugal): the effect of weather events and management in mercury transport

The main aim of this research was to assess the mercury transport from an estuarine basin with a background of anthropogenic contamination during a spring tidal cycle (year 2009) and compare it with two previous tidal cycles (years 1994 and 1999), as part of a long‐term monitoring program. Results showed that effective mercury transport occurs both in the dissolved and particulate fractions (0.18 and 0.20 kg per tidal cycle, respectively), and despite an overall decrease in environmental contamination, results more than double previous findings on particulate transport in the system. These findings result essentially from changes in the tidal prism (net export of 2 million m³ of water), given that both dissolved and particulate concentrations did not increase over time. Hydrodynamic simulations were performed to evaluate the effect of physical disturbance (dredging) and weather events (increased freshwater flow) in these processes, and results suggest the increased freshwater flow into the system as the main forcing function for the mercury transport increment. These results highlight the importance of long‐term monitoring programs, since despite an overall improvement in local contamination levels, the enhancement of transport processes through hydrological changes increases environmental pressure away from the contamination source. Copyright © 2012 John Wiley & Sons, Ltd.     

The groundwater recharge regime of some slightly metamorphosed neoproterozoic sedimentary rocks: an application of natural environmental tracers

Isotope data of precipitation and groundwater in parts of the Voltaian Basin in Northern Ghana were used to explain the groundwater recharge regime in the area. Groundwater recharge is an important parameter in the development of a decision support system for the management and efficient utilization of groundwater resources in the area. It is therefore important to establish the processes and sources of groundwater recharge. δ¹⁸O and δ²H data for local precipitation suggest enrichment relative to the Global Meteoric Water Line (GMWL) and indicate that precipitation takes place at a relative humidity less than 100%. The groundwater data plot on an evaporation line with a slope of 5, suggesting a high degree of evaporative enrichment of the precipitation in the process of vertical infiltration and percolation through the unsaturated zone into the saturated zone. This finding is consistent with the observation of high evapotranspiration rates in the area and ties in with the fact that significant clay fraction in the unsaturated zone limits vertical percolation and thus exposes the percolating rainwater to the effects of high temperatures and low humidities resulting in high evapotranspiration rates. Groundwater recharge estimates from the chloride mass balance, CMB, method suggest recharge in the range of 1.8–32% of the annual average precipitation in the form of rainfall. The highest rates are associated with areas where open wells encourage significant amount of groundwater recharge from precipitation in the area. In the northern parts of the study area, groundwater recharge is lower than 12%. The recharge so computed through the application of the CMB methodology takes on a spatial distribution akin to the converse of the spatial pattern of both δ¹⁸O and δ²H in the area. As such, the locations of the highest recharge are associated with the most depleted values of the two isotopes. This observation is consistent with the assertion that low vertical hydraulic conductivities slow down vertical percolation of precipitation down to the groundwater water. The percolating precipitation water thus gets enriched in the heavier isotopes through high evapotranspiration rates. At the same time, the amount of water that finally reaches the water table is considerably reduced. Copyright © 2013 John Wiley & Sons, Ltd.     

Torsional dynamic response of a large‐diameter pipe pile in viscoelastic saturated soil

An analytical solution is developed in this paper to investigate the dynamic response of a large‐diameter end‐bearing pipe pile subjected to torsional loading in viscoelastic saturated soil. The wave propagation in saturated soil and pile are simulated by Biot's two‐phased linear theory and one‐dimensional elastic theory, respectively. The dynamic equilibrium equations of the outer soil, inner soil, and pile are established. The solutions for the outer and inner soils in frequency domain are obtained by Laplace transform technique and the separation of variables method. Then, the dynamic response of the pile is obtained on the basis of the perfect contacts between the pile and the outer soil as well as the inner soil. The results in this paper are compared with that of a solid pile in elastic saturated soil to verify the validity of the solution. Furthermore, the solution in this paper is compared with the classic plane strain solution to verify the solution further and check the accuracy of the plane strain solution. Numerical results are presented to analyze the vibration characteristics and illustrate the effect of the soil parameters and the geometry size of the pile on the complex impedance and velocity admittance of the pile head. Finally, the displacement of the soil at different depth and frequency is analyzed. Copyright © 2014 John Wiley & Sons, Ltd.