Sediment and particulate phosphorus characteristics in grassed waterways from row crop corn and alfalfa fields collected by manual University of Exeter samplers and automatic sampling

Phosphorus (P) export from agricultural lands above known threshold levels can result in adverse impacts to receiving water quality. Phosphorus loss occurs in dissolved and sediment‐bound, or particulate phosphorous (PP), forms, with the latter often dominating losses from row‐cropped systems. To target practices, land managers need good computer models and model developers need good monitoring data. Sediment monitoring data (e.g. radiometric finger printing and sediment P sorption capacity) can help identify sediment source areas and improve models, but require more sediment mass than is typically obtained by automatic sampling. This study compares a simple suspended sediment sampler developed at the University of Exeter (UE) with automatic sampling in intermittent channels draining corn and alfalfa fields. The corn field had a greater runoff coefficient (27%) than alfalfa (11%). No differences were found in enrichment ratios (sediment constituent/soil constituent) in PP (PP_ER) or percent loss on ignition (LOI_ER) between paired UE samplers on corn. The median LOI_ER for the UE samplers (1·9%) did not differ significantly (p > 0·13) from the automatic sampler (2·0%). The PP_ER from the UE samplers was on average 20% lower than the automatic samplers. A correlation (r² = 0·75) was found between sediment PP and % LOI from automatic samplers and UE samplers for particles < 50 µm, while for > 50 µm PP concentration did not change with changes in % LOI. Sediment ammonium‐oxalate extractable metals were similarly related to LOI, with the strongest correlation for iron (r² = 0·71) and magnesium (r² = 0·70). Copyright © 2011 John Wiley & Sons, Ltd.     

Stage‐discharge uncertainty derived with a non‐stationary rating curve in the Choluteca River,Honduras

Uncertainty in discharge data must be critically assessed before data can be used in, e.g. water resources estimation or hydrological modelling. In the alluvial Choluteca River in Honduras, the river‐bed characteristics change over time as fill, scour and other processes occur in the channel, leading to a non‐stationary stage‐discharge relationship and difficulties in deriving consistent rating curves. Few studies have investigated the uncertainties related to non‐stationarity in the stage‐discharge relationship. We calculated discharge and the associated uncertainty with a weighted fuzzy regression of rating curves applied within a moving time window, based on estimated uncertainties in the observed rating data. An 18‐year‐long dataset with unusually frequent ratings (1268 in total) was the basis of this study. A large temporal variability in the stage‐discharge relationship was found especially for low flows. The time‐variable rating curve resulted in discharge estimate differences of ? 60 to + 90% for low flows and ± 20% for medium to high flows when compared to a constant rating curve. The final estimated uncertainty in discharge was substantial and the uncertainty limits varied between ? 43 to + 73% of the best discharge estimate. Copyright © 2010 John Wiley & Sons, Ltd.     

Quantifying the effects of climate variability and human activities on runoff from the Laohahe basin in northern China using three different methods

Much attention has recently been focused on the effects that climate variability and human activities have had on runoff. In this study, these effects are quantified using three methods, namely, multi‐regression, hydrologic sensitivity analysis, and hydrologic model simulation. A conceptual framework is defined to separate the effects. As an example, the change in annual runoff from the semiarid Laohahe basin (18 112 km²) in northern China was investigated. Non‐parametric Mann‐Kendall test, Pettitt test, and precipitation‐runoff double cumulative curve method were adopted to identify the trends and change‐points in the annual runoff from 1964 to 2008 by first dividing the long‐term runoff series into a natural period (1964–1979) and a human‐induced period (1980–2008). Then the three quantifying methods were calibrated and calculated, and they provided consistent estimates of the percentage change in mean annual runoff for the human‐induced period. In 1980–2008, human activities were the main factors that reduced runoff with contributions of 89–93%, while the reduction percentages due to changes in precipitation and potential evapotranspiration only ranged from 7 to 11%. For the various effects at different durations, human activities were the main reasons runoff decreased during the two drier periods of 1980–1989 and 2000–2008. Increased runoff during the wetter period of 1990–1999 is mainly attributed to climate variability. This study quantitatively separates the effects of climate variability and human activities on runoff, which can serve as a reference for regional water resources assessment and management. Copyright © 2011 John Wiley & Sons, Ltd.     

Spatio‐temporal distribution of near‐surface and root zone soil moisture at the catchment scale

Soil moisture is highly variable both spatially and temporally. It is widely recognized that improving the knowledge and understanding of soil moisture and the processes underpinning its spatial and temporal distribution is critical. This paper addresses the relationship between near‐surface and root zone soil moisture, the way in which they vary spatially and temporally, and the effect of sampling design for determining catchment scale soil moisture dynamics. In this study, catchment scale near‐surface (0–50 mm) and root zone (0–300 mm) soil moisture were monitored over a four‐week period. Measurements of near‐surface soil moisture were recorded at various resolutions, and near‐surface and root zone soil moisture data were also monitored continuously within a network of recording sensors. Catchment average near‐surface soil moisture derived from detailed spatial measurements and continuous observations at fixed points were found to be significantly correlated (r² = 0·96; P = 0·0063; n = 4). Root zone soil moisture was also found to be highly correlated with catchment average near‐surface, continuously monitored (r² = 0·81; P < 0·0001; n = 26) and with detailed spatial measurements of near‐surface soil moisture (r² = 0·84). The weaker relationship observed between near‐surface and root zone soil moisture is considered to be caused by the different responses to rainfall and the different factors controlling soil moisture for the soil depths of 0–50 mm and 0–300 mm. Aspect is considered to be the main factor influencing the spatial and temporal distribution of near‐surface soil moisture, while topography and soil type are considered important for root zone soil moisture. The ability of a limited number of monitoring stations to provide accurate estimates of catchment scale average soil moisture for both near‐surface and root zone is thus demonstrated, as opposed to high resolution spatial measurements. Similarly, the use of near‐surface soil moisture measurements to obtain a reliable estimate of deeper soil moisture levels at the small catchment scale was demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.     

A high strain‐rate constitutive model for sand and its application in finite‐element analysis of tunnels subjected to blast

The paper presents a constitutive model for simulating the high strain‐rate behavior of sands. Based on the concepts of critical‐state soil mechanics, the bounding surface plasticity theory and the overstress theory of viscoplasticity, the constitutive model simulates the high strain‐rate behavior of sands under uniaxial, triaxial and multi‐axial loading conditions. The model parameters are determined for Ottawa and Fontainebleau sands, and the performance of the model under extreme transient loading conditions is demonstrated through simulations of split Hopkinson pressure bar tests up to a strain rate of 2000/s. The constitutive model is implemented in a finite‐element analysis software Abaqus to analyze underground tunnels in sandy soil subjected to internal blast loads. Parametric studies are conducted to examine the effect of relative density and type of sand and of the depth of tunnel on the variation of stresses and deformations in the soil adjacent to the tunnels. Copyright © 2012 John Wiley & Sons, Ltd.     

Extreme value modelling of daily areal rainfall over Mediterranean catchments in a changing climate

Heavy rainfall events during the fall season are causing extended damages in Mediterranean catchments. A peaks‐over‐threshold model is developed for the extreme daily areal rainfall occurrence and magnitude in fall over six catchments in Southern France. The main driver of the heavy rainfall events observed in this region is the humidity flux (FHUM) from the Mediterranean Sea. Reanalysis data are used to compute the daily FHUM during the period 1958–2008, to be included as a covariate in the model parameters. Results indicate that the introduction of FHUM as a covariate can improve the modelling of extreme areal precipitation. The seasonal average of FHUM can improve the modelling of the seasonal occurrences of heavy rainfall events, whereas daily FHUM values can improve the modelling of the events magnitudes. In addition, an ensemble of simulations produced by five different general circulation models are considered to compute FHUM in future climate with the emission scenario A1B and hence to evaluate the effect of climate change on the heavy rainfall distribution in the selected catchments. This ensemble of climate models allows the evaluation of the uncertainties in climate projections. By comparison to the reference period 1960–1990, all models project an amplification of the mean seasonal FHUM from the Mediterranean Sea for the projection period 2070–2099, on average by +22%. This increase in FHUM leads to an increase in the number of heavy rainfall events, from an average of 2.55 events during the fall season in present climate to 3.57 events projected for the period 2070–2099. However, the projected changes have limited effects on the magnitude of extreme events, with only a 5% increase in the median of the 100‐year quantiles. Copyright © 2011 John Wiley & Sons, Ltd.     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

The impact of highly permeable layer on hydraulic system in a coastal aquifer

This study is aimed to understand the hydraulic mechanism of coastal aquifer systems that include highly permeable layers (HPLs). These hydrologic conditions can be found in many volcanic islands that are composed of a series of lava flows discharged into sea or other standing body of water. In the first part, we developed a numerical model based on the geologic and hydrologic data obtained from the eastern Jeju Island, Korea, of which the aquifer contains clinker and hyaloclastite layers. The simulation results reproduced spatial location of fresh‐saline water interface, especially the abrupt decline of interface at the inland part and the thickness variation of transition zone along the cross‐section observed at the eastern Jeju coastal aquifer. We were able to find out that these phenomena are strongly related to the presence of the HPL. In the second part, quantitative analyses were conducted with the use of hypothetical models in order to understand the dynamic characteristics of coastal system that includes HPLs. A series of sensitivity studies were conducted to assess the effect of the horizontal length and vertical depth of HPL on the spatial location of the interface toe and the configuration of transition zone. Various case studies have shown that the seawater intruded into the inland more as the horizontal length of HPL was increased and its vertical depth was decreased. In other simulations including two HPLs, the vertical distance between these two HPLs primarily controlled the flow regime, flux variations, and the configuration of the transition zone. Finally, we performed simulations to evaluate the effect of a rising sea‐level. This study provides more understanding of how the presence of HPL controls the seawater intrusion processes, and the spatial configurations of fresh‐saline water interface at coastal aquifers. Copyright © 2012 John Wiley & Sons, Ltd.     

On seasonal approach to flood frequency modelling. Part II: flood frequency analysis of Polish rivers

The annual peak flow series of Polish rivers are mixtures of summer and winter flows. As Part II of a sequence of two papers, practical aspects of applicability of seasonal approach to flood frequency analysis (FFA) of Polish rivers are discussed. Taking A Two‐Component Extreme Value (TCEV1) model as an example it was shown in the first part that regardless of estimation method, the seasonal approach can give profit in terms of upper quantile estimation accuracy that rises with the return period of the quantile and is the greatest for no seasonal variation. In this part, an assessment of annual maxima (AM) versus seasonal maxima (SM) approach to FFA was carried out with respect to seasonal and annual peak flow series of 38 Polish gauging stations. First, the assumption of mutual independence of the seasonal maxima has been tested. The smoothness of SM and AM empirical probability distribution functions was analysed and compared. The TCEV1 model with seasonally estimated parameters was found to be not appropriate for most Polish data as it considerably underrates the skewness of AM distributions and upper quantile values as well. Consequently, the discrepancies between the SM and AM estimates of TCEV1 are observed. Taking SM and TCEV1 distribution, the dominating season in AM series was confronted with predominant season for extreme floods. The key argument for presumptive superiority of SM approach that SM samples are more statistically homogeneous than AM samples has not been confirmed by the data. An analysis of fitness to SM and AM of Polish datasets made for seven distributions pointed to Pearson (3) distribution as the best for AM and Summer Maxima, whereas it was impossible to select a single best model for winter samples. In the multi‐model approach to FFA, the tree functions, i.e., Pe(3), CD3 and LN3, should be involved for both SM and AM. As the case study, Warsaw gauge on the Vistula River was selected. While most of AM elements are here from winter season, the prevailing majority of extreme annual floods are the summer maxima. The upper quantile estimates got by means of classical annual and two‐season methods happen to be fairly close; what's more they are nearly equal to the quantiles calculated just for the season of dominating extreme floods. Copyright © 2011 John Wiley & Sons, Ltd.     

Response of a tunnel deeply embedded in a viscoelastic medium

This paper presents a three‐dimensional energy‐based solution for the time‐dependent response of a deeply embedded and unsupported semi‐infinite tunnel of circular cross‐section. The tunnel is taken to be excavated quasi‐instantaneously from an infinite rock body that initially exhibits an isotropic stress state and that is made up of a homogeneous, isotropic and viscoelastic material. The viscoelastic behaviour is modelled by means of Burger's model, and the rock is taken to behave volumetrically linear elastic and to exhibit exclusively deviatoric creep. This viscoelastic problem is transformed into the Laplace domain, where it represents a quasi‐elastic problem. The displacement fields in the new solution are taken to be the products of independent functions that vary in the radial and longitudinal directions. The differential equations governing the displacements of the system and appropriate boundary conditions are obtained using the principle of minimum potential energy. The solutions for these governing equations in the Laplace domain are then obtained analytically and numerically using a one‐dimensional finite difference technique. The results are then transformed back into the time domain using an efficient numerical scheme. The accuracy of the new solution is comparable with that of a finite element analysis but requires much less computation effort. Copyright © 2011 John Wiley & Sons, Ltd.