ALBERT

Description: We examined the relationships between large wood (LW) export and precipitation patterns and intensity by analyzing the data on the annual volume of LW removed from 42 reservoirs and the daily precipitation at or near the reservoir sites. We also calculated the effective precipitation by considering the antecedent precipitation. Both daily and effective precipitation data were used as explanatory variables to explain LW export. The model selection revealed that the precipitation pattern and intensity controlling LW export varied with latitude in the Japanese archipelago. In small watersheds with narrow channel widths and low discharges, mass movements, such as landslides and debris flows, are major factors in the production and transport of LW. In this case, the effective precipitation required to initiate mass movements regulated the LW export and did not vary with the latitude. In intermediate and large watersheds with wide channel widths and high stream discharges, heavy rainfall and subsequent floods regulated buoyant depth, influencing the initiation of LW movement. In southern and central Japan, intense rainfall accompanied by typhoons or localized torrential downpours causes geomorphic disturbances, which introduce abundant pieces of LW into the channels. However, these pieces continue to be removed by repeated rainfall events. Therefore, LW export is supply-limited and potentially produces less LW accumulation. Conversely, in northern Japan, where typhoons and torrential downpours are rare, LW export is transport-limited because LW pieces recruited by bank erosion, tree mortality, and windthrow accumulate and persist on valley floors. These pieces may be easily exported by infrequent flooding.

Description: We developed a method to measure in situ the isotopic composition of liquid water with minimal supervision and, most important, with a temporal resolution of less than a minute. For this purpose a microporous hydrophobic membrane contactor (Membrana) was combined with an isotope laser spectrometer (Picarro). The contactor, originally designed for degassing liquids, was used with N2 as a carrier gas in order to transform a small fraction of liquid water to water vapor. The generated water vapor was then analyzed continuously by the Picarro analyzer. To prove the membrane's applicability, we determined the specific isotope fractionation factor for the phase change through the contactor's membrane across an extended temperature range (8°C–21°C) and with different waters of known isotopic compositions. This fractionation factor is needed to subsequently derive the liquid water isotope ratio from the measured water vapor isotope ratios. The system was tested with a soil column experiment, where the isotope values derived with the new method corresponded well (R2 = 0.998 for δ18O and R2 = 0.997 for δ2H) with those of liquid water samples taken simultaneously and analyzed with a conventional method (cavity ring-down spectroscopy). The new method supersedes taking liquid samples and employs only relatively cheap and readily available components. This makes it a relatively inexpensive, fast, user-friendly, and easily reproducible method. It can be applied in both the field and laboratory wherever a water vapor isotope analyzer can be run and whenever real-time isotope data of liquid water are required at high temporal resolution.

Description: Transit time of discharge is a hydrological characteristic used in water resource management. Previous studies have demonstrated large spatial variation in the mean transit time (MTT) of stream base flow in meso-scale catchments. Various relationships between topography and MTT have been reported. Although it is generally assumed that base flow MTT is controlled by the depth of the hydrologically active layer that recharges a stream, this hypothesis has not been tested in field studies. This study confirmed that the depth of hydrologically active soil and bedrock controls spatial variation in MTT. The study used isotopic and geochemical tracer data gathered in the 4.27 km2 Fudoji catchment, central Japan. The results, together with previously documented relationships between topography and MTT, indicate that the depth of the hydrologically active layer is sometimes, but not always, related to topography. A comprehensive understanding of the factors that control base flow production in mountainous catchments will require further study of the water flow path depths that recharge streams.

Description: This paper examines the long-term historical changes in frequency and amplitude of hydroclimatic extremes in the Blue Nile basin using data from the second half of 20th century. The temporal variability of basin-wide rainfall extremes and river flow extremes from four gauging stations was investigated under the hypothesis of no trend and no persistence in time. On the basis of a quantile anomaly analysis method, decadal variations in extreme daily, monthly, and annual quantiles were studied, and the periods of statistical significance were identified. The analysis showed that high and low river flows and rainfall depths do not vary in time in a fully random way but show a particular variation pattern. Their extremes show significant decadal variations. The 1980s had statistically significant negative anomalies in extremes in comparison with the long-term reference period of 1964–2009, while the 1960s–1970s and the 1990s–2000s had positive anomalies, although less significant. There is neither consistent increasing nor decreasing trend in rainfall and flow extremes of recent years. Therefore, anticipated trends due to global warming could not be identified. Conversely, low-flow extremes show an increasing trend during the last decade, which could be related to the effect of water regulation works at the outlet of Lake Tana. Moreover, similar patterns and statistically significant correlations were found between climatic indices representing the Pacific and Atlantic Oceans and the Blue Nile rainfall and flow extremes. Changes that occur on the Pacific Ocean appear to be a main driver for the decadal oscillations in climate and related high and low Blue Nile water availability for Ethiopia, Sudan, and Egypt.

Description: Few studies exist on infiltration processes in badlands, although infiltration and subsurface lateral flows are known to contribute to soil erosion and to control slope instability. Our investigation was carried out in a 100 m² plot located in a 0.5 ha landslide in black marls (South-East France). An artificial sprinkling was performed with an intensity of 10 mm.h -1 during 66.4 h interrupted with 8.4 h. breaks. KBr and KCl were used as tracers. A pseudo-steady state was reached after 25-35 hours and 250-350 mm of rainfall. The runoff coefficient was 40% (ratio total runoff volume/total sprinkling water amount). Pre-event water (PE) contributed to the groundwater recharge at the very beginning of the experiment but PE contribution dropped steadily while the soil was saturating. After around 200 mm cumulative rainfall, PE contribution started to rise steeply before reaching a nearly constant value. This original mechanism implies an efficient transfer process of PE. It was assumed from the description of the material structure and from hydrological evidences that PE was mainly drained from a structure porosity made of the marl's flaked nature. Total pre-event water contributions ranged from 25 to 79 % (PE contribution was over 50 % in 2/3 of the observations wells). Over the recession phase, release of pre-event water occurred from the drainage of a texture porosity. The study showed that at the plot scale, infiltration processes proved well organised despite the high heterogeneity and anisotropy of the material. It was possible to propose a general conceptual model explaining the hydrological processes over time and area. The peculiar structure of regolith originating from black marl is preserved over a large part of the weathering time, so that the material structure (type, orientation of grains, small/large pores) remains a first order control of water flow generation in Black marl soils. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In this work, the multifractal properties of hourly rainfall data recorded at a location in Southern Spain have been related to the scale properties of the corresponding Intensity-Duration-Frequency (IDF) curves. Four parametric models for the IDF curves have been fitted to the quantiles of rainfall obtained by using the Generalized Pareto frequency distribution function with the extreme data series obtained for the same place. The scaling of the rainfall intensity moments has been analysed and the empirical moments scaling exponent function has been obtained. The corresponding values of q 1 and γ 1 have been empirical and theoretically calculated and compared to some characteristics of the different IDF models. Thus, the scaling behaviour of IDF curves has been analysed and the best model has been selected. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Major-ion compositions of groundwater are employed in this study of the water-rock interactions and hydrogeochemical evolution within a carbonate aquifer system. The groundwater samples were collected from boreholes or underground tunnels in the Ordovician limestone of Yanzhou Coalfield where catastrophic groundwater inflows can be hazardous to mining and impact use of the groundwater as a water supply. The concentration of TDS ranged from 961 to 3,555 mg/l and indicates moderately to highly mineralized water. The main water-type of the middle Ordovician limestone groundwater is Ca-Mg-SO 4 , with SO 4 2- ranging from 537 mg/l to 2,297 mg/l, and average values of Ca 2+ and Mg 2+ of 455.7 mg/l and 116.6 mg/l, respectively. The water samples were supersaturated with respect to calcite and dolomite and undersaturated or saturated with respect to gypsum. Along the general flow direction, deduced from increases of TDS and Cl - , the main water-rock interactions that caused hydrogeochemical evolution of the groundwater within the aquifer were the dissolution of gypsum, the precipitation of calcite, the dissolution or precipitation of dolomite, and ion exchange. Ion exchange is the major cause for the lower mole concentration of Ca 2+ than that of SO 4 2- . The groundwater level of Ordovician aquifer is much higher than that of C-P coal-bearing aquifers, so the potential flow direction is upward and the pyrite in coal is not a possible source of sulfate, additional data on the stable sulfur and oxygen isotopic composition of the sulfate may be helpful to identify its origin. Although ion exchange probably accounts for the higher mole concentration of Na + than that of Cl - , the dissolution of aluminosilicate can not be ruled out. The data evaluation methods and results of this study could be useful in other areas to understand flow paths in aquifers and to provide information needed to identify the origin of groundwater. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Quantitative evaluation of the effect of climate variability and human activities on runoff is of great importance for water resources planning and management in term of maintaining the ecosystem integrity and sustaining the society development. In this paper, hydro-climatic data from 4 catchments (i.e., Luanhe River Catchment, Chaohe River Catchment, Hutuo River Catchment and Zhanghe River Catchment) in the Haihe River Basin from 1957–2000 were used to quantitatively attribute the hydrological response (i.e., runoff) to climate change and human activities separately. To separate the attributes, the temporal trends of annual precipitation, potential evapotranspiration (PET) and runoff during 1957–2000 were first explored by the Man-Kendall test. Despite that only Hutuo River Catchment was dominated by significant negative trend in annual precipitation, all four catchments presented significant negative trend in annual runoff varying from −0.859 mm a -1 (Chaohe River) to −1.996 mm a -1 (Zhanghe River). Change points in 1977 and 1979 are detected by precipitation-runoff double cumulative curves method and Pettitt's test for Zhanghe River and other three rivers, respectively, and are adopted to divide data set into two study periods as the pre-change period and post-change period. Three methods including hydrological model method, hydrological sensitivity analysis method and climate elasticity method were calibrated with the hydro-climatic data during the pre-change period. And then, hydrological runoff response to climate variability and human activities were quantitatively evaluated with the help of the three methods and based on the assumption that climate and human activities are the only drivers for streamflow and are independent to each other. Similar estimates of anthropogenic and climatic effects on runoff for catchments considered can be obtained from the three methods. We found that human activities were the main driving factors for the decline in annual runoff in Luanhe River Catchment, Chaohe River Catchment and Zhanghe River Catchment, accounting for over 50% of runoff reduction. However, Climate variability should be responsible for the decrease in annual runoff in the Hutuo River Catchment. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Measurements of river water quality at Haridwar, India, taken during 2005 – 2006 show that the river water does not meet the WHO and Bureau of Indian Standards criteria of drinking water quality, especially with respect to total coliform and fecal-coliform. This study investigates the removal of pathogens at a river bank filtration (RBF) site in Haridwar. Using the quality of river water and the quality of abstracted water from a nearby production well, semi-empirical models based on the concept of filtration coefficient are developed and tested for their effectiveness in removing pathogens under varying bacteriological quality of source water. A two-tier model which includes the effect of clogged layer is developed to obtain an equivalent filtration coefficient. This coefficient is found to be linearly related with natural logarithm of the concentration of pathogens in the source water. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Nonpoint source pollution and hydromodification are the leading causes of impairment to our nation's rivers and streams. Roadside ditch networks, ubiquitous in both rural and urban landscapes, intercept and shunt substantial quantities of overland runoff and shallow groundwater to stream systems. By altering natural flowpaths, road ditches contribute not only to hydromodification but also potentially to nonpoint-source (NPS) pollution by acting as hydrologic links between agricultural fields and natural streams. Unfortunately, the impacts of these alterations on watershed hydrology and water quality are not well understood. Through a series of field measurements, including field surveys and discharge monitoring, this study examined the effect of road ditch networks on basin morphometry, field- and watershed-scale hydrology, and pollutant transport in a 38 km 2 agricultural watershed in south-central NY. Salient findings include: (i) 94% of road ditches discharged directly to natural streams, effectively doubling the drainage density, (ii) on average, road ditches increased peak and total event flows in their receiving streams by 78% and 57%, respectively, but displayed significant variation across ditches, (iii) ditches intercepted large quantities of surface and subsurface runoff from agricultural fields and therefore represent efficient conduits for the transport of agricultural NPS pollutants to sensitive receiving waterbodies. Our results provide useful information for hydrologists who wish to further understand how artificial drainage may be affecting watershed hydrology and for managers and engineers tasked with designing appropriate flood and NPS pollution control measures. Copyright © 2012 John Wiley & Sons, Ltd.

Description: We coupled the process-based NIES Integrated Catchment-based Eco-hydrology (NICE) model to an urban canopy model (UCM) and the Regional Atmospheric Modeling System (RAMS) in order to simulate the effect of urban geometry and anthropogenic exhaustion on the hydrothermal changes in the atmospheric/land and the interfacial areas of the Japanese megalopolis. The simulation was conducted with multi-scale in horizontally regional–urban–point levels, and in vertically atmosphere–surface–unsaturated–saturated layers. The model reproduced reasonably the observed hydrothermal values by using ground-truth data in various types of natural/artificial land covers. The simulated results also suggested that the latent heat flux in new water-holding pavement (consisting of porous asphalt and water-holding filler made of steel by-products based on silica compound) has a strong impact on hydrologic cycle and cooling temperature in comparison with the observed heat budget by newly incorporating the effect of water amount on the heat conductivity in the pavement. Furthermore, the model predicted the hydrothermal changes under two types of land cover scenarios to promote evaporation and to reduce air temperature against heat island phenomenon. Finally, we evaluated the relationship between the effect of groundwater use to ameliorate the heat island and the effect of infiltration on the water cycle in the catchment. These procedures to integrate the multi-scaled model simulation with political scenario based on the effective management of water resources as heat sink/source would be very powerful approaches to recovering a sound hydrologic cycle and create thermally-pleasing environments in the megalopolis. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Soil and vadose zone profiles are used as an archive of changes in groundwater recharge and water quality following changes in land use in an area of the Loess Plateau of China. A typical rain-fed loess-terrace agriculture region in Hequan, Guyuan is taken as an example and multiple tracers (chloride mass balance, stable isotopes, tritium, and water chemistry) are used to examine groundwater recharge mechanisms and to evaluate soil water chloride as an archive for recharge rate and water quality. Results show that groundwater recharge beneath natural uncultivated grassland, used as a baseline, is about 94–100 mm yr -1 and the time it takes for annual precipitation to reach water table through the thick unsaturated zone is from decades to hundreds of years (tritium free). This recharge rate is 2–3 orders of magnitude more than in the other semiarid areas with similar annual rainfall and with deep rooted vegetation and relatively high temperature. Most of the water that eventually becomes recharge originally infiltrated in the summer months. The conversion from native grassland to winter wheat has reduced groundwater recharge by 42–50% (50–55 mm yr -1 for recharge), and from winter wheat to alfalfa resulted in a significant chloride accumulation in the upper soil zone which terminated deep drainage. The paper also evaluates the time lag between potential recharge and actual recharge to aquifer and between increase in solute concentration in soil moisture and that in the aquifer following land-use change due to the deep unsaturated zone. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Shallow upland drains, grips, have been hypothesized as responsible for increased downstream flow magnitudes. Observations provide counterfactual evidence, often relating to the difficulty of inferring conclusions from statistical correlation and paired catchment comparisons; and the complexity of designing field experiments to test grip impacts at the catchment-scale. Drainage should provide drier antecedent moisture conditions, providing more storage at the start of an event; but, grips have higher flow velocities than overland flow so potentially delivering flow more rapidly to the drainage network. We develop and apply a model for assessing the impacts of grips upon flow hydrographs. The model was calibrated on the gripped case; then the gripped case was compared with the intact case by removing all grips. This comparison showed that even given parameter uncertainty, the intact case had significantly higher flood peaks and lower baseflows, mirroring field observations of the hydrological response of intact peat. The simulations suggest that this is because delivery effects may not to translate into catchment-scale impacts for three reasons. First, in our case, the proportions of flow path lengths that were hillslope were not changed significantly by gripping. Second, the structure of the grip network as compared with the structure of the drainage basin mitigated against grip-related increases in the concentration of runoff in the drainage network, although it did marginally reduce the mean timing of that concentration at the catchment outlet. Third, the effect of the latter upon downstream flow magnitudes can only be assessed by reference to the peak timing of other tributary basins, emphasizing that drain effects are both relative and scale dependent. However, given the importance of hillslope flow paths, we show that if upland drainage causes significant changes in surface roughness on hillslopes, then critical and important feedbacks may impact upon the speed of hydrological response. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Natural soil pipes are found in peatlands but little is known about their hydrological role. This paper presents the most complete set of pipe discharge data to date from a deep blanket peatland in northern England. In a 17.4-ha catchment, we identified 24 perennially-flowing and 60 ephemerally-flowing pipe outlets. Eight pipe outlets along with the catchment outlet were continuously gauged over an 18-month period. The pipes in the catchment were estimated to produce around 13.7 % of annual streamflow with individual pipes often producing large peak flows (maximum peak of 3.8 L s -1 ). Almost all pipes, whether ephemeral, perennially-flowing, shallow or deep (outlets 〉 1 m below the peat surface), showed increased discharge within a mean of 3 hours after rainfall commencement and were dominated by stormflow, indicating good connectivity between the peatland surface and the pipes. However, almost all pipes had a longer time period between hydrograph peak and return to baseflow compared to the stream (mean of 23.9 hours for pipes, 19.7 hours for stream). As a result, the proportion of streamflow produced by the pipes at any given time increased at low flows and formed the most important component of stream discharge for the lowest 10 % of flows. Thus, a small number of perennially-flowing pipes became more important to the stream system under low flow conditions and probably received water via matrix flow during periods between storms. Given the importance of pipes to streamflow in blanket peatlands, further research is required into their wider role in influencing stream-water chemistry, water temperature and fluvial carbon fluxes, as well as their role in altering local hydrochemical cycling within the peat mass itself. Enhanced piping within peatlands caused by environmental change may lead to changes in streamflow regime with larger low flows and more prolonged drainage of the peat. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The hydrology of Quebec (Canada) boreal fens is poorly documented. Many peatlands are located in watersheds with impounded rivers. In such cases, their presence influences reservoir inflows. In recent years, some fens have been subjected to an increase of their wet area, a sign that they may be evolving towards an aquatic ecosystem. This dynamic process is referred to as aqualysis. This paper presents the seasonal and monthly hydrological budgets of a small watershed including a highly aqualysed fen (James Bay region). Monitoring of precipitation ( P ), runoff ( Q ) and groundwater levels ( WL ) was conducted during the ice-free season. Three semi-empirical equations (Thornthwaite, Priestley-Taylor and Penman-Monteith) were used and compared to calculate potential evapotranspiration ( PET ). The first two equations, having fewer parameters, estimate higher PET values than the third equation. The use of pressure level gauges installed in wells, for the calculation of peatland water storage, is inconclusive. Swelling of peat, peat decomposition and plant composition could be responsible for non-negligible amounts of absorbed water, which are not entirely accounted for by well levels. The estimation of peat matrix water storage is potentially the largest source of error and the limiting factor to calculate water balances in this environment. The results show that the groundwater level and the water storage vary depending on the season and especially after a heavy rainfall. Finally, the results illustrate the complexity of water routing through the site and, thus, raise several questions to be resolved in the future. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Potential climate change effects on aspects of conjunctive management of water resources can be evaluated by linking climate models with fully integrated groundwater–surface water models. The objective of this study is to develop a modeling system that links global climate models with regional hydrologic models, using the California Central Valley as a case study. The new method is a supply and demand modeling framework that can be used to simulate and analyze potential climate change and conjunctive use. Supply-constrained and demand-driven linkages in the water system in the Central Valley are represented with the linked climate models, precipitation-runoff models, agricultural and native vegetation water use, and hydrologic flow models to demonstrate the feasibility of this method. Simulated precipitation and temperature were used from the GFDL-A2 climate change scenario through the 21st century to drive a regional water balance mountain hydrologic watershed model (MHWM) for the surrounding watersheds in combination with a regional integrated hydrologic model of the Central Valley (CVHM). Application of this method demonstrates the potential transition from predominantly surface water to groundwater supply for agriculture with secondary effects that may limit this transition of conjunctive use. The particular scenario considered includes intermittent climatic droughts in the first half of the 21st century followed by severe persistent droughts in the second half of the 21st century. These climatic droughts do not yield a valley-wide operational drought but do cause reduced surface water deliveries and increased groundwater abstractions that may cause additional land subsidence, reduced water for riparian habitat, or changes in flows at the Sacramento–San Joaquin River Delta. The method developed here can be used to explore conjunctive use adaptation options and hydrologic risk assessments in regional hydrologic systems throughout the world.

Description: An efficient approach is developed to analytically evaluate solute transport in a horizontal, divergent radial flow field with a multistep injection flow rate and an arbitrary input concentration history. By assuming a piecewise steady state flow and transforming the time domain to the cumulative injected flow domain, the concentration distribution is found to be completely determined by the total volume of injected flow and independent of specific flow rates. Thus, on the cumulative flow domain, the transport problem with a temporally varying velocity field can be transformed into a steady state flow problem. Linear convolution can then be applied on the cumulative injected flow domain to evaluate the solution for an arbitrarily time-dependent input concentration. Solutions on the regular time domain can be conveniently obtained by mapping the solution on the cumulative injected flow domain to the time domain. Furthermore, we theoretically examine the conditions for the assumption of piecewise steady state flow to be valid. On the basis of the critical time scale of the “pseudosteady state condition,” defined as when velocity changes accomplish 99% of their steady state differences, and the relative error in the mean travel time of plume front, we obtain conditions for neglecting the transitional period between two pumping steps. Such conditions include the following: (1) the duration of a pumping step, tp, must be longer than the critical time scale, tc, i.e., tp ≥ tc = 25r2S/T, where r is the radial distance, S is the storage coefficient, and T is the transmissivity, or similarly, a maximum problem domain needs to be defined for a given pumping strategy. (2) the maximum well pumping rate, qmax, should satisfy qmax ≤ πθT/25S, where θ is the effective porosity. When both conditions are satisfied, transitional periods may be neglected.

Description: Meander bends in alluvial rivers morphologically evolve toward meander cutoff with narrowing intra-meander necks, and this should steepen hydraulic gradients and intensify intra-meander hyporheic flux. This research used dye tracking and head loss measurements in a 1:500 planimetrically scaled laboratory river table to quantify the spatial and temporal intensification of intra-meander flux rates at two evolution ages. The younger meander bend, M1, had a sinuosity of 2.3, a river neck width of 0.39 cm, and 0.6% river slope, and the older meander bend, M3, had a sinuosity of 5.2, a river neck width of 0.12 cm, and 0.5% river slope. Flux into and out of the meander bend was estimated along the normalized curvilinear distance s *, with the meander neck at s * = 0.1 and s * = 0.9, the meander centroid at s * = 0.37 and s * = 0.63, and the apex at s * = 0.5. Between the meander centroid and neck we documented a 60% spatial intensification for M1 and a 90% spatial intensification for M3. Between M1 and M3 we documented a 135% temporal intensification at the neck and a 100% intensification at the centroid. Our empirical spatial and temporal intensification rates involving the M1 and the M3 scenario were 1 to 3 times lower than theoretical rates derived from a river evolution model with equivalent M1 and M3 planimetry. Over estimation by the theoretical model was attributed to exaggerated head loss caused by the model neglecting groundwater contributions to river stage. Hyporheic exchange provides critical ecosystem services and its spatial and temporal variation with meander evolution should be considered in river management. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Watershed models which combine hydrology and water quality are being widely utilized in integrated watershed management for the determination of best water management practices. In this study, the hydrology of the Lower Porsuk Stream Watershed in Turkey has been modeled with SWAT to determine optimal water management strategies. The calibration and validation process have been accomplished using data from two monitoring stations. The model has been run for the 1978–2009 period and while the 1998–2004 period has been used for calibration, the validation has spanned the whole period. The SWATCup calibration and uncertainity program has been utilized for this purpose. No significant differences have been detected among different iteration numbers in the calibration period. The monthly Nash-Sutcliffe and R 2 performance indicatiors for the upstream Esenkara station have been 0.74 and 0.88, respectively for the calibration period, and 0.87 and 0.87, respectively for the validation period. The Kıranharmanı station which is located close to the watershed outlet has shown values of 0.59 and 0.72, respectively for the calibration period, and 0.44 and 0.56, respectively for the validation period. There are uncertanities in the abstracted irrigation and groundwater quantities which have reflected in the results in the Kıranharmanı station which is more affected as it lies downstream of the irrigation areas. The effects of different irrigation practices on the flow regime have been also investigated. A scenario has been implemented in which drip irrigation wholly replaces conventional furrow and sprinkler irrigation. The scenario has shown increases in stream flows by 87 % for the whole year. The adoption of more efficient irrigation practices thus results in reducing the water stress induced by irrigation demands. With this study a modeling framework has been founded to aid water management applications in the Lower Porsuk Stream Watershed by generating scenarios for best management practices. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Many observations and studies have shown that water resources amount in the Hai River Basin decreased significantly over the last half of the twentieth century. This study attempts to attribute the observed changes in the water resources amount in the basin over a 40 year period (1961–2000) to different factors, including natural climate variability, climate change induced by anthropogenic forcing of greenhouse gas emissions (referred to as anthropogenic forcing hereafter), and local human activity. First, the temporal variation of the annual water resources amount in the basin during the past 40 years is analyzed by employing the moving-average method, the linear regression method, and the Mann-Kendall method. Second, through setting different scenarios, the effects on the water resources amount due to different factors, including natural climate variability, anthropogenic forcing, and local human activity, are obtained using the parallel climate model, the distributed hydrological model water and energy transfer processes in large river basins, and the statistical downscaling model. Third, the fingerprint-based attribution method is used to obtain the signal strengths of observed changes in water resources amount during 1961–2000 and changes in the water resources amount under different scenarios. Finally, by comparing the signal strengths, the observed changes in water resources amount in the basin can be attributed to different factors. The results indicate that natural climate variability and local human activity may be two factors responsible for the observed changes in the water resources amount during the past 40 years in the basin, with local human activity being the main factor and accounting for about 60% of the changes.

Description: Flood events can induce temporal changes in streambed elevation and particle-size composition, which may influence the bed's hydraulic properties and stream-aquifer fluxes during and after an event. This study combines a set of previously developed modeling approaches to create a synthetic flood event during which bed sediment is entrained and deposited as a function of hydraulic conditions and particle size. One simulated river reach in a state of approximate dynamic equilibrium is chosen to investigate the impacts of size-selective sediment transport on stream-aquifer interaction. Along this reach, the preferential entrainment of fine sediment during the flood's rising limb leads to overall bed coarsening, and increases in vertical hydraulic conductivity (Kbv) and downward fluxes of floodwater into the streambed. Progressively finer sediment layers are deposited during the event's falling limb, causing the redevelopment of a colmation (clogging) layer on the bed surface and a decline in overall Kbv by the event's conclusion. This reduction in Kbv leads to prolonged retention of event water in the streambed (after the reach reverts from losing to gaining river conditions) when compared with what is expected if pre-event Kbv values are used to estimate river-aquifer exchanges. This process of sequential bed coarsening and fining during a flood event provides a mechanistic explanation for the event size-and-duration threshold, inferred in some systems, that must be exceeded for significant amounts of flood recharge to occur. The major consequences of these processes—enhanced infiltration and prolonged floodwater retention—have potentially major implications for groundwater-surface water interactions, water quality, contaminant transport, and riparian biogeochemistry.

Description: Nonstationary oscillation (NSO) processes are observed in a number of hydroclimatic data series. Stochastic simulation models are useful to study the impacts of the climatic variations induced by NSO processes into hydroclimatic regimes. Reproducing NSO processes in a stochastic time series model is, however, a difficult task because of the complexity of the nonstationary behaviors. In the current study, a novel stochastic simulation technique that reproduces the NSO processes embedded in hydroclimatic data series is presented. The proposed model reproduces NSO processes by utilizing empirical mode decomposition (EMD) and nonparametric simulation techniques (i.e., k-nearest-neighbor resampling and block bootstrapping). The model was first tested with synthetic data sets from trigonometric functions and the Rössler system. The North Atlantic Oscillation (NAO) index was then examined as a real case study. This NAO index was then employed as an exogenous variable for the stochastic simulation of streamflows at the Romaine River in the province of Quebec, Canada. The results of the application to the synthetic data sets and the real-world case studies indicate that the proposed model preserves well the NSO processes along with the key statistical characteristics of the observations. It was concluded that the proposed model possesses a reasonable simulation capacity and a high potential as a stochastic model, especially for hydroclimatic data sets that embed NSO processes.

Description: Vegetation zonation and tidal hydrology are basic attributes of intertidal salt marshes, but specific links among vegetation zonation, plant water use, and spatiotemporally dynamic hydrology have eluded thorough characterization. We developed a quantitative model of an intensively studied salt marsh field site, integrating coupled 2-D surface water and 3-D groundwater flow and zonal plant water use. Comparison of model scenarios with and without heterogeneity in (1) evapotranspiration rates and rooting depths, according to mapped vegetation zonation, and (2) sediment hydraulic properties from inferred geological heterogeneity revealed the coupled importance of both sources of ecohydrological variability at the site. Complex spatial variations in root zone pressure heads, saturations, and vertical groundwater velocities emerged in the model but only when both sources of ecohydrological variability were represented together and with tidal dynamics. These regions of distinctive root zone hydraulic conditions, caused by the intersection of vegetation and sediment spatial patterns, were termed “ecohydrological zones” (EHZ). Five EHZ emerged from different combinations of sediment hydraulic properties and evapotranspiration rates, and two EHZ emerged from local topography. Simulated pressure heads and groundwater dynamics among the EHZ were validated with field data. The model and data showed that hydraulic differences between EHZ were masked shortly after a flooding tide but again became prominent during prolonged marsh exposure. We suggest that ecohydrological zones, which reflect the combined influences of topographic, sediment, and vegetation heterogeneity and do not emphasize one influence over the others, are the fundamental spatial habitat units comprising the salt marsh ecosystem.

Description: Quantification of biophysical parameters of urban trees is important for urban planning, and for assessing carbon sequestration and ecosystem services. Airborne lidar has been used extensively in recent years to estimate biophysical parameters of trees in forested ecosystems. However, similar studies are largely lacking for individual trees in urban landscapes. Prediction models to estimate biophysical parameters such as height, crown area, diameter at breast height, and biomass for over two thousand individual trees were developed using best subsets multiple linear regression for a study area in central Oklahoma, USA using point cloud distributional metrics from an Optech ALTM 2050 lidar system. A high level of accuracy was attained for estimating individual tree height (R2 = 0.89), dbh (R2 = 0.82), crown diameter (R2 = 0.90), and biomass (R2 = 0.67) using lidar-based metrics for pooled data of all tree species. More variance was explained in species-specific estimates of biomass (R2 = 0.68 for Juniperus virginiana to 0.84 for Ulmus parviflora) than in estimates from broadleaf deciduous (R2 = 0.63) and coniferous (R2 = 0.45) taxonomic groups—or the data set analysed as a whole (R2 = 0.67). The metric crown area performed particularly well for most of the species-specific biomass equations, which suggests that tree crowns should be delineated accurately, whether manually or using automatic individual tree detection algorithms, to obtain a good estimation of biomass using lidar-based metrics.

Description: In recent decades, copula functions have been applied in bivariate drought duration and severity frequency analysis. Among a number of potential copulas, Clayton has been mostly employed in drought analysis. In this research, we study the influence of the tail shape of various copula functions (i.e. Gumbel, Frank, Clayton, and Gaussian) on drought bivariate frequency analysis. The appropriateness of Clayton copula for the characterization of drought characteristics is also investigated. Drought data are extracted from standardized precipitation index (SPI) time series for four stations in Canada (La Tuque and Grande Prairie) and Iran (Anzali and Zahedan). Both duration and severity datasets are positively skewed. Different marginal distributions were first fitted to drought duration and severity data. The gamma and exponential distributions were respectively selected for drought duration and severity according to the positive skewness and Kolmogorov- Smirnov test. The results of copula modeling show that the Clayton copula function is not an appropriate choice for the employed datasets in the current study, and does not give more drought risk information than an independent model for which the duration and severity dependence is not significant. The reason is that the dependence of two variables in the upper tail of Clayton copula is very weak and similar to the independent case, while the observed data in the transformed domain of cumulative density function shows high association in the upper tail. Instead, the Frank and Gumbel copula functions show better performance than Clayton function for drought bivariate frequency analysis. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Soil hydraulic parameters were upscaled from a 30 m resolution to a 1 km resolution using a new aggregation scheme (described in the companion paper) where the scale parameter was based on the topography. When soil hydraulic parameter aggregation or upscaling schemes ignore the effect of topography, their application becomes limited at hillslope scales and beyond, where topography plays a dominant role in soil deposition and formation. Hence the new upscaling algorithm was tested at the hillslope scale (1 km) across two locations: (1) the Little Washita watershed in Oklahoma, and (2) the Walnut Creek watershed in Iowa. The watersheds were divided into pixels of 1 km resolution and the effective soil hydraulic parameters obtained for each pixel. Each pixel/domain was then simulated using the physically based HYDRUS-3-D modeling platform. In order to account for the surface (runoff/on) and subsurface fluxes between pixels, an algorithm to route infiltration-excess runoff onto downstream pixels at daily time steps and to update the soil moisture states of the downstream pixels was applied. Simulated soil moisture states were compared across scales, and the coarse scale values compared against the airborne soil moisture data products obtained during the hydrology experiment field campaign periods (SGP97 and SMEX02) for selected pixels with different topographic complexities, soil distributions, and land cover. Results from these comparisons show good correlations between simulated and observed soil moisture states across time, topographic variations, location, elevation, and land cover. Stream discharge comparisons made at two gauging stations in the Little Washita watershed also provide reasonably good results as to the suitability of the upscaling algorithm used. Based only on the topography of the domain, the new upscaling algorithm was able to provide coarse resolution values for soil hydraulic parameters which effectively captured the variations in soil moisture across the watershed domains.

Description: Hydropower accounts for about 20% of the worldwide electrical power production. In mountainous regions this ratio is significantly higher. In this study we present how future projected climatic forcing, as described in regional climate models (RCMs), will affect water resources and subsequently hydropower production in downstream hydropower plants in a glacierized alpine valley (Vispa valley, Switzerland, 778 km2). In order to estimate future runoff generation and hydropower production, we used error-corrected and downscaled climate scenarios from regional climate models (RCMs) as well as glacier retreat projections from a dynamic glacier model and coupled them to a physically based hydrological model. Furthermore, we implemented all relevant hydropower operational rules in the hydrological model to estimate future hydropower production based on the runoff projections. The uncertainty of each modeling component (climate projections, glacier retreat, and hydrological projection) and the resulting propagation of uncertainty to the projected future water availability for energy production were assessed using an analysis of variance. While the uncertainty of the projections is considerable, the consistent trends observed in all projections indicate significant changes to the current situation. The model results indicate that future melt- and rainfall-runoff will increase during spring but decline during summer. The study concludes by outlining the most relevant expected changes for hydropower operations.

Description: The effect of spatial concentration fluctuations on the reaction of two solutes, A + B ⇀ C, is considered. In the absence of fluctuations, the concentration of solutes decays as Adet = Bdet ∼ t−1. Contrary to this, experimental and numerical studies suggest that concentrations decay significantly slower. Existing theory suggests a t−d/4 scaling in the asymptotic regime (d is the dimensionality of the problem). Here we study the effect of fluctuations using the classical diffusion-reaction equation with random initial conditions. Initial concentrations of the reactants are treated as correlated random fields. We use the method of moment equations to solve the resulting stochastic diffusion-reaction equation and obtain a solution for the average concentrations that deviates from ∼t−1 to ∼t−d/4 behavior at characteristic transition time t*. We also derive analytical expressions for t* as a function of Damköhler number and the coefficient of variation of the initial concentration.

Description: Based on existing techniques in nonlinear physics that work in the Fourier domain, we develop a multivariate, wavelet-based method for the generation of synthetic discharge time series. This approach not only retains the cross-correlative structure of the original data (which makes it preferable to principal component methods that merely preserve the correlations) but also replicates the nonlinear properties of the original data. We argue that the temporal asymmetry of the typical hydrograph is the most important form of nonlinearity to preserve in the synthetic data. Using the derivative skewness as a measure of asymmetry and an example data set of 35 years of daily discharge data from 107 gauging stations in the United States, we compare two approaches that preserve the asymmetry of the original records. We generate synthetic data and then study the properties of fitting a generalized extreme value distribution to the annual maxima for a total flux time series. The synthetic series provides error bands for the fitted distribution that give a different way of assessing credible return periods. It is found that the best approach for studying extremes is to match the asymmetry of each series individually, rather than to formulate a global threshold criterion.

Description: This article focuses on household water use in Spain by analyzing the influence of a detailed set of factors. We find that, although the presence of both water-saving equipment and water-conservation habits leads to water savings, the factors that influence each are not the same. In particular, our results show that those individuals most committed to the adoption of water-saving equipment and, at the same time, less committed to water-conservation habits tend to have higher incomes.

Description: Wetlands are valuable ecosystems that provide many valuable services, yet many of these important ecosystems are at risk because of current trends in climate change. The Prairie Pothole Region (PPR) in the upper-midwest of the U.S. and south-central Canada, characterized by glacially-sculpted landscapes and abundant wetlands, is one such vulnerable region. According to regional/global climate model predictions, drought occurrence will increase in the PPR region through the 21st century and thus will probably cause the amount of water in wetlands to decline. Water surface area (WSA) of Kidder County, ND from 1984-2011 was measured by classifying TM/ETM+ images through the Modified Normalized Difference Water Index (MNDWI).We then developed a linear model based on the water surface area (WSA) of these wetlands and historical climate data, and used this to determine the wetland sensitivity to climate change and predict future wetlands WSA in the PPR. Our model based on Palmer Drought Severity Index (PDSI) of the current year (PDSI t-0 ) and two years previous (PDSI t-2 ) can explain 79% of the annual wetland WSA variance, suggesting a high sensitivity of wetlands to drought/climate change. We also predicted the PPR wetlands WSA in the 21st century under A1B scenario (a mid-carbon emission scenario) using simulated PDSI based on IPCC AR4 22-model ensemble climate. According to our prediction, the WSA of the PPR wetlands will decrease to less than half of the baseline WSA (defined as the mean wetlands WSA of the 2000s) by the mid of the 21st century, and to less than one-third by the 2080s, and will then slightly increase in the 2090s. This considerable future wetland loss caused only by climate change provides important implication to future wetland management and climate adaptation policy. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In hydrology, the storage-discharge relationship is a fundamental catchment property. Understanding what controls this relationship is at the core of catchment science. To date, there are no direct methods to measure water storage at catchment scales (10 1 -10 3  km 2 ). In this study, we use direct measurements of terrestrial water storage dynamics by means of superconducting gravimetry in a small headwater catchment of the Regen River, Germany, to derive empirical storage-discharge relationships in nested catchments of increasing scale. Our results show that the local storage measurements are strongly related to streamflow dynamics at larger scales (〉 100 km 2 ; correlation coefficient = 0.78-0.81), but at small scale no such relationship exists (~ 1 km 2 ; correlation coefficients = -0.11). The geologic setting in the region can explain both the disconnection between local water storage and headwater runoff, and the connectivity between headwater storage and streams draining larger catchment areas. More research is required to understand what controls the form of the observed storage-discharge relationships at the catchment scale. This study demonstrates that high-precision gravimetry can provide new insights into the complex relationship between state and response of hydrological systems. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In semiarid regions, the rooting strategies employed by vegetation can be critical to its survival. Arid regions are characterized by high variability in the arrival of rainfall, and species found in these areas have adapted mechanisms to ensure the capture of this scarce resource. Vegetation roots have strong control over this partitioning, and assuming a static root profile, predetermine the manner in which this partitioning is undertaken.A coupled, dynamic vegetation and hydrologic model, tRIBS + VEGGIE, was used to explore the role of vertical root distribution on hydrologic fluxes. Point-scale simulations were carried out using two spatially and temporally invariant rooting schemes: uniform: a one-parameter model and logistic: a two-parameter model. The simulations were forced with a stochastic climate generator calibrated to weather stations and rain gauges in the semiarid Walnut Gulch Experimental Watershed (WGEW) in Arizona. A series of simulations were undertaken exploring the parameter space of both rooting schemes and the optimal root distribution for the simulation, which was defined as the root distribution with the maximum mean transpiration over a 100-yr period, and this was identified. This optimal root profile was determined for five generic soil textures and two plant-functional types (PFTs) to illustrate the role of soil texture on the partitioning of moisture at the land surface. The simulation results illustrate the strong control soil texture has on the partitioning of rainfall and consequently the depth of the optimal rooting profile. High-conductivity soils resulted in the deepest optimal rooting profile with land surface moisture fluxes dominated by transpiration. As we move toward the lower conductivity end of the soil spectrum, a shallowing of the optimal rooting profile is observed and evaporation gradually becomes the dominate flux from the land surface. This study offers a methodology through which local plant, soil, and climate can be accounted for in the parameterization of rooting profiles in semiarid regions.

Description: There has been a recent debate in the hydrological community about the relative merits of the informal generalized likelihood uncertainty estimation (GLUE) approach to uncertainty assessment in hydrological modeling versus formal probabilistic approaches. Some recent literature has suggested that the methods can give similar results in practice when properly applied. In this note, we show that the connection between formal Bayes and GLUE is not merely operational but goes deeper, with GLUE corresponding to a certain approximate Bayesian procedure even when the “generalized likelihood” is not a true likelihood. The connection we describe relates to recent approximate Bayesian computation (ABC) methods originating in genetics. ABC algorithms involve the use of a kernel function, and the generalized likelihood in GLUE can be thought of as relating to this kernel function rather than to the model likelihood. Two interpretations of GLUE emerge, one as a computational approximation to a Bayes procedure for a certain “error-free” model and the second as an exact Bayes procedure for a perturbation of that model in which the truncation of the generalized likelihood in GLUE plays a role. The intent of this study is to encourage cross-fertilization of ideas regarding GLUE and ABC in hydrologic applications. The connection we outline suggests the possibility of combining a formal likelihood with a kernel based on a generalized likelihood within the ABC framework and also allows advanced ABC computational methods to be used in GLUE applications. The model-based interpretation of GLUE may also be helpful in partially illuminating the implicit assumptions in different choices of generalized likelihood.

Description: Satellite-passive microwave remote sensing has been extensively used to estimate snow water equivalent (SWE) in northern regions. Although passive microwave sensors operate independent of solar illumination and the lower frequencies are independent of atmospheric conditions, the coarse spatial resolution introduces uncertainties to SWE retrievals due to the surface heterogeneity within individual pixels. In this article, we investigate the coupling of a thermodynamic multilayered snow model with a passive microwave emission model. Results show that the snow model itself provides poor SWE simulations when compared to field measurements from two major field campaigns. Coupling the snow and microwave emission models with successive iterations to correct the influence of snow grain size and density significantly improves SWE simulations. This method was further validated using an additional independent data set, which also showed significant improvement using the two-step iteration method compared to standalone simulations with the snow model.

Description: In this paper we present a visual analytics approach for deriving spatio-temporal patterns of collective human mobility from a vast mobile network traffic data set. More than 88 million movements between pairs of radio cells—so-called handovers—served as a proxy for more than two months of mobility within four urban test areas in Northern Italy. In contrast to previous work, our approach relies entirely on visualization and mapping techniques, implemented in several software applications. We purposefully avoid statistical or probabilistic modeling and, nonetheless, reveal characteristic and exceptional mobility patterns. The results show, for example, surprising similarities and symmetries amongst the total mobility and people flows between the test areas. Moreover, the exceptional patterns detected can be associated to real-world events such as soccer matches. We conclude that the visual analytics approach presented can shed new light on large-scale collective urban mobility behavior and thus helps to better understand the “pulse” of dynamic urban systems.

Description: Existing analytical solutions to determine aquifer response to a change in stream stage are inappropriate where an unsaturated zone exists beneath the stream, as in the case of disconnected stream-aquifer systems. A better understanding of the relationship between aquifer response and transient stream stage in disconnected systems is therefore required, as this would also aid in the field determination of the status of connection between the stream and aquifer. We use a numerical model to examine transient stream stage and the corresponding water table response. Beneath disconnected streams, the magnitude of head change in the water table level is a balance between the cumulative infiltration during a flow event and the rate at which the water can disperse laterally. Increases in wave duration, stream width, and streambed permeability result in greater infiltrated water volume and therefore a higher peak response at the water table. Conversely, higher aquifer transmissivity and aquifer hydraulic conductivity allow the water to move laterally away from the stream faster, resulting in a smaller head change below the stream. Lower unsaturated storage results in a greater and faster aquifer response because the unsaturated zone can fill more quickly. Under some combinations of parameters, the magnitude of the disconnected head response is more than seven times greater than the change in stream stage driving streambed infiltration; an effect which can never occur beneath a connected stream. The results of this sensitivity analysis are compared to field data from a river in eastern Australia to determine periods of disconnection. Where the change in aquifer head is greater than the change in stream stage, disconnection between the stream and aquifer can be determined.

Description: The ability to understand and predict the flux and fate of sediment delivered to the sea by rivers remains an outstanding scientific challenge. Approaches to this challenge are necessarily synthetic, spanning wide ranges in spatial and temporal scales. Here a conventional sediment transport theory used by engineers and sedimentologists at reach and channel scales is applied at the basin scale. Specifically, a straightforward expression proposed by Bagnold and modified accordingly predicts the observed importance of combined wetness and steepness of a source basin as a control of sediment supply to the sea. The reasonable, key assumption underlying the application of sediment transport theory in this context is that the river-mouth sites for which suspended-sediment loads are reported are alluviated, and thus characterized by transport-limited flux of sediment. This analysis also indicates the potential significance of additional, as yet poorly documented factors constraining sediment supply to the sea. These factors, some of which appear to covary systematically with climate, include river-profile concavity, river-mouth channel width and friction, and the characteristic size of sediment in transport.

Description: Key documents such as the European Water Framework Directive and the U.S. Clean Water Act state that public and stakeholder participation in water resource management is required. Participation aims to enhance resource management and involve individuals and groups in a democratic way. Evaluation of participatory programs and projects is necessary to assess whether these objectives are being achieved and to identify how participatory programs and projects can be improved. The different methods of evaluation can be classified into three groups: (i) process evaluation assesses the quality of participation process, for example, whether it is legitimate and promotes equal power between participants, (ii) intermediary outcome evaluation assesses the achievement of mainly nontangible outcomes, such as trust and communication, as well as short- to medium-term tangible outcomes, such as agreements and institutional change, and (iii) resource management outcome evaluation assesses the achievement of changes in resource management, such as water quality improvements. Process evaluation forms a major component of the literature but can rarely indicate whether a participation program improves water resource management. Resource management outcome evaluation is challenging because resource changes often emerge beyond the typical period covered by the evaluation and because changes cannot always be clearly related to participation activities. Intermediary outcome evaluation has been given less attention than process evaluation but can identify some real achievements and side benefits that emerge through participation. This review suggests that intermediary outcome evaluation should play a more important role in evaluating participation in water resource management.

Description: A planning support system for land consolidation has been developed that has, at its heart, an expert system called LandSpaCES (Land Spatial Consolidation Expert System) which contains a “design module” that generates alternative land redistributions under different scenarios and an “evaluation module” which integrates GIS with multi-criteria decision making for assessing these alternatives. This paper introduces the structural framework of the latter module which has been applied using a case study in Cyprus. Two new indices are introduced: the “parcel concentration coefficient” for measuring the dispersion of parcels; and the “landowner satisfaction rate” for predicting the acceptance of the land redistribution plan by the landowners in terms of the location of their new parcels. These two indices are used as criteria for the evaluation of the land redistribution alternatives and are transferable to any land consolidation project. Moreover, a modified version of the ratio estimation procedure, referred to as the “qualitative rating method” for assigning weights to the evaluation criteria, is presented, along with a set of non-linear value functions for standardizing the performance scores of the alternatives and incorporating expert knowledge for five evaluation criteria. The application of the module showed that it is a powerful new tool for the evaluation of alternative land redistribution plans that could be implemented in other countries after appropriate adjustments. A broader contribution has also been made to spatial planning processes, which might follow the methodology and innovations presented in this paper.

Description: The key question that is asked in this study is “how are the three independent bias components of satellite rainfall estimation, comprising hit bias, missed, and false precipitation, physically related to the estimation uncertainty of soil moisture and runoff for a physically based hydrologic model?” The study also investigated the performance of different satellite rainfall products as a function of land use and land cover (LULC) type. Using the entire Mississippi river basin as the study region and the variable infiltration capacity (VIC)-3L as the distributed hydrologic model, the study of the satellite products (CMORPH, 3B42RT, and PERSIANN-CCS) yielded two key findings. First, during the winter season, more than 40% of the rainfall total bias is dominated by missed precipitation in forest and woodland regions (southeast of Mississippi). During the summer season, 51% of the total bias is governed by the hit bias, and about 42% by the false precipitation in grassland-savanna region (western part of Mississippi basin). Second, a strong dependence is observed between hit bias and runoff error, and missed precipitation and soil moisture error. High correlation with runoff error is observed with hit bias (∼0.85), indicating the need for improving the satellite rainfall product's ability to detect rainfall more consistently for flood prediction. For soil moisture error, it is the total bias that correlated significantly (∼0.78), indicating that a satellite product needed to be minimized of total bias for long-term monitoring of watershed conditions for drought through continuous simulation.

Description: Bayesian inference is used to study the effect of precipitation and model structural uncertainty on estimates of model parameters and confidence limits of predictive variables in a conceptual rainfall-runoff model in the snow-fed Rudbäck catchment (142 ha) in southern Finland. The IHACRES model is coupled with a simple degree day model to account for snow accumulation and melt. The posterior probability distribution of the model parameters is sampled by using the Differential Evolution Adaptive Metropolis (DREAM(ZS)) algorithm and the generalized likelihood function. Precipitation uncertainty is taken into account by introducing additional latent variables that were used as multipliers for individual storm events. Results suggest that occasional snow water equivalent (SWE) observations together with daily streamflow observations do not contain enough information to simultaneously identify model parameters, precipitation uncertainty and model structural uncertainty in the Rudbäck catchment. The addition of an autoregressive component to account for model structure error and latent variables having uniform priors to account for input uncertainty lead to dubious posterior distributions of model parameters. Thus our hypothesis that informative priors for latent variables could be replaced by additional SWE data could not be confirmed. The model was found to work adequately in 1-day-ahead simulation mode, but the results were poor in the simulation batch mode. This was caused by the interaction of parameters that were used to describe different sources of uncertainty. The findings may have lessons for other cases where parameterizations are similarly high in relation to available prior information.

Description: Operational land surface models (LSMs) compute hydrologic states such as soil moisture that are needed for a range of important applications (e.g., drought, flood, and weather prediction). The uncertainty in LSM parameters is sufficiently great that several researchers have proposed conducting parameter estimation using globally available remote sensing data to identify best fit local parameter sets. However, even with in situ data at fine modeling scales, there can be significant remaining uncertainty in LSM parameters and outputs. Here, using a new uncertainty estimation subsystem of the NASA Land Information System (LIS) (described herein), a Markov chain Monte Carlo (MCMC) technique is applied to conduct Bayesian analysis for the accounting of parameter uncertainties. The Differential Evolution Markov Chain (DE-MC) MCMC algorithm was applied, for which a new parallel implementation was developed. A case study is examined that builds on previous work in which the Noah LSM was calibrated to passive (L-band) microwave remote sensing estimates of soil moisture for the Walnut Gulch Experimental Watershed. In keeping with prior related studies, the parameters subjected to the analysis were restricted to the soil hydraulic properties (SHPs). The main goal is to estimate SHPs and soil moisture simulation uncertainty before and after consideration of the remote sensing data. The prior SHP uncertainty is based on the original source of the standard SHP lookup tables for the Noah LSM. Conclusions are drawn regarding the value and viability of Bayesian analysis over alternative approaches (e.g., parameter estimation, lookup tables) and further research needs are identified.

Description: It has long been known that surface gravity waves induce significant seepage through the porous layer found at the lake bottom. Away from coastal regions, however, the pressure signature of surface waves at the lake bottom is weak. We consider fully nonlinear internal gravity waves, whose long wavelength and slow motion implies a sustained and strong pressure perturbation even in the deep regions of the lake. We argue that internal waves can induce significant seepage through the sediment layer, in regions where surface gravity waves have negligible impact. The pressure profile at the fluid–porous layer interface is computed from the “exact” Dubreil-Jacotin-Long theory, giving a reliable profile even for large waves. This profile is used in conjunction with Darcy's law to compute the seepage within the porous region. We find that the geometric distribution of seepage is strongly controlled by both the ratio of porous media thickness to the horizontal length of the pressure perturbation, and the bottom topography, when it is present. Based on work on the interaction of internal solitary waves with the bottom boundary layer, we develop a model to account for the changes in permeability due to wave-induced instabilities in the bottom boundary layer and enhanced benthic turbulence. This turbulence acts to unplug the pores near the surface by lifting the detritus that clogs them. The resulting changes in permeability significantly enhance exchange between the free fluid and the porous medium on the downstream side of the wave.

Description: Investigations of solute transport in fractured rock aquifers often rely on tracer test data acquired at a limited number of observation points. Such data do not, by themselves, allow detailed assessments of the spreading of the injected tracer plume. To better understand the transport behavior in a granitic aquifer, we combine tracer test data with single-hole ground-penetrating radar (GPR) reflection monitoring data. Five successful tracer tests were performed under various experimental conditions between two boreholes 6 m apart. For each experiment, saline tracer was injected into a previously identified packed-off transmissive fracture while repeatedly acquiring single-hole GPR reflection profiles together with electrical conductivity logs in the pumping borehole. By analyzing depth-migrated GPR difference images together with tracer breakthrough curves and associated simplified flow and transport modeling, we estimate (1) the number, the connectivity, and the geometry of fractures that contribute to tracer transport, (2) the velocity and the mass of tracer that was carried along each flow path, and (3) the effective transport parameters of the identified flow paths. We find a qualitative agreement when comparing the time evolution of GPR reflectivity strengths at strategic locations in the formation with those arising from simulated transport. The discrepancies are on the same order as those between observed and simulated breakthrough curves at the outflow locations. The rather subtle and repeatable GPR signals provide useful and complementary information to tracer test data acquired at the outflow locations and may help us to characterize transport phenomena in fractured rock aquifers.

Description: : Exempt domestic water well provisions exist in sixteen states and allow landowners to use water wells producing relatively small amounts of water or water wells used for domestic or other specific purposes to avoid certain regulatory requirements. Exempt wells provisions apply mainly in the west, where most states use the prior appropriation rule (“first in time, first in right”) for ground water rights. Critics voice concerns over implications ranging from land use planning to water supply. This article explores existing state regulation and finds that the term “exempt well” is a misnomer. The vast majority of states impose numerous restrictions on exempt wells. Existing regulations include limitations on the quantity, location, and amount of irrigation use. Other provisions address construction standards, required filing of information and other practices. The review of regulations suggests that many steps to mitigate the impact of exempt wells are already in place.

Description: : National attention has turned to a case currently pending before the New Mexico Supreme Court. Bounds v. State of New Mexico involves a facial constitutional challenge to New Mexico's domestic well statute. The Plaintiff, Horace Bounds, claims that the statute, which requires the State Engineer to issue permits to any applicants for domestic purposes, violates the constitutional provision known as the Prior Appropriation Doctrine. The decision of the Court will likely have impacts across the country and beyond the legal realm. This decision could impact the well drilling industry, economic development, and administrative agencies. Numerous other western states have similar exempt well provisions, and the decision by the New Mexico Supreme Court could well lead to litigation or pushes for changes in legislation on other states as well.

Description: : Over the last decade, Montana's fastest growing basins have experienced a surge of reliance on a certain “permit-exempt well loophole,” which enables small ground water diversions to avoid having to go through the regulatory process that determines whether the proposed use would adversely impact senior water users and ground water supply. This loophole has generated significant conflict with senior right owners in areas where exempt wells have provided water for large-scale residential developments. As a partial solution, this article examines ground water mitigation exchanges designed to address similar issues in Washington and presents the 63rd Montana Legislature with recommendations for a Montana pilot project.

Description: : Most western states allow landowners to withdraw certain amounts of groundwater for specified purposes without obtaining a water right permit. As water demands increase in the West, some developers are using exempt wells to supply water to residential developments without acquiring the permits and water rights needed to build public water supply systems. This has led to concerns that the cumulative impact of these wells could impair senior water rights, create environmental problems, and threaten water quality in some areas. This article explores the considerations associated with mitigating exempt well impacts and argues that surgical approaches focused on specific issues and geographic areas (scalpels) are more politically and administratively feasible than broad efforts to repeal or reduce exemptions on a statewide basis (hammers). The views expressed in this article are those of the author alone and do not represent those of the Western States Water Council or its member states.

Description: : The saga over exempt wells in the western United States and Canada epitomizes a new type of water conflict – a spaghetti-western water war. The political melodrama stars local governments to serve as sheriff of water-supply planning duties. Exempt wells number in the millions, and herding the growing numbers is testing the mettle of the states and provinces responsible for the management, allocation, and protection of natural resources. The separation of laws governing ground water and surface water, coupled with changes in geography and geology within a jurisdiction, compound the administrative riddle and give rise to a broad spectrum of conflicts, from differing interpretations of hydrogeologic data, economic impacts associated with increasing the herd, to differing identities associated with the use of ground water from the exempt wells. Despite the political melodrama of exempt wells, there is room and willingness for other trails and paths to keep the herd intact. This paper describes the different breeds of conflicts associated with exempt wells and gives examples of how the mysterious stranger of collaborative decision making processes and water governance systems can ride into town and lead to successful water management and conflict resolution.

Description: : The water management practice of storing water underground for use during dry periods is a very effective technique that has been practiced for a long time. Another water management practice that has been around for a long time is collection of rainwater and storage of rainwater aboce ground. There have been significant advances more recently in our understanding of aquifer geology and technologies used for rainwater harvesting. This article combines our advance in knowledge of storing and recovery water from wells and new rainwater harvesting practices to recharge a low yielding exempt well with treated rainwater. The studies presented in this article demonstrate that exempt wells provide a valuable underground water storage option and may provide a new effective urban water management tool.

Description: Predicting long-term consequences of climate change on hydrologic processes has been limited due to the needs to accommodate the uncertainties in hydrological measurements for calibration, and to account for the uncertainties in the models that would ingest those calibrations and uncertainties in climate predictions as basis for hydrological predictions. We implemented a hierarchical Bayesian (HB) analysis to coherently admit multiple data sources and uncertainties including data inputs, parameters, and model structures to identify the potential consequences of climate change on soil moisture and streamflow at the head watersheds ranging from low to high elevations in the southern Appalachian region of the United States. We have considered climate change scenarios based on three greenhouse gas emission scenarios of the Interovernmental Panel on Climate Change: A2, A1B and B1 emission scenarios. Full predictive distributions based on hierarchical Bayesian models are capable of providing rich information and facilitating the summarization of prediction uncertainties. With predictive uncertainties taken into account, the most pronounced change in soil moisture and streamflow would occur under the A2 scenario at both low and high elevations, followed by the A1B scenario and then by the B1 scenario. Uncertainty in the change of soil moisture is less than that of streamflow for each season, especially at high elevations. A reduction of soil moisture in summer and fall, a reduction or slight increase of streamflow in summer, and an increase of streamflow in winter are predicted for all three scenarios at both low and high elevations. The hydrological predictions with quantified uncertainties from a HB model could aid more-informed water resource management in developing mitigation plans, and dealing with water security under climate change. Copyright © 2012 John Wiley & Sons, Ltd.

Description: This research proposes a combination of SWAT and MODFLOW, MD-SWAT-MODFLOW, to address the multi-aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi-aquifers so that the daily streamflow can be simulated correctly. A multi-aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in-between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river-aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD-SWAT-MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the LH-OAT sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi-aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio-temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Large areas of Amazonian evergreen forest experience seasonal droughts extending for three or more months, yet show maximum rates of photosynthesis and evapotranspiration during dry intervals. This apparent resilience is belied by disproportionate mortality of the large trees in manipulations that reduce wet season rainfall, occurring after 2–3 years of treatment. The goal of this study is to characterize the mechanisms that produce these contrasting ecosystem responses. A mechanistic model is developed based on the ecohydrological framework of TIN (Triangulated Irregular Network)-based Real Time Integrated Basin Simulator + Vegetation Generator for Interactive Evolution (tRIBS+VEGGIE). The model is used to test the roles of deep roots and soil capillary flux to provide water to the forest during the dry season. Also examined is the importance of “root niche separation,” in which roots of overstory trees extend to depth, where during the dry season they use water stored from wet season precipitation, while roots of understory trees are concentrated in shallow layers that access dry season precipitation directly. Observational data from the Tapajós National Forest, Brazil, were used as meteorological forcing and provided comprehensive observational constraints on the model. Results strongly suggest that deep roots with root niche separation adaptations explain both the observed resilience during seasonal drought and the vulnerability of canopy-dominant trees to extended deficits of wet season rainfall. These mechanisms appear to provide an adaptive strategy that enhances productivity of the largest trees in the face of their disproportionate heat loads and water demand in the dry season. A sensitivity analysis exploring how wet season rainfall affects the stability of the rainforest system is presented.

Description: We present direct numerical simulation results for both isothermal and density-stratified turbulent flow in an open channel into and around a 120° bend with a bulk Reynolds number of 7500 and Prandtl number of 1.5. The bend is sharp with a radius-to-channel breadth ratio of 1.5. The bulk Richardson number for the stratified flow is 2.4 based on overall channel depth. The gradient Richardson number (Rig) varies between 10 and 20 at the entrance to θ ≈ 60°, where θ is the angular location. Above θ ≈ 60 − 120, Rig ≈ 1. In isothermal flow, the well-known helical flow structure is observed. In stratified conditions, the vertical variation in relative strength of the outward-directed baroclinic pressure gradient and the centrifugal acceleration leads to a more complex circulation structure. In the near bed region and immediately above the interface, the centrifugal acceleration is greater, driving flow radially inward, while just below the density interface the baroclinic pressure gradient is greater, leading to outward-directed flow. This produces a four layer circulation structure with potentially significant implications for sediment erosion and transport. Additionally, this produces a complex dynamic at the density interface where the shear orientation varies through approximately 200° over the mixing layer depth.

Description: In the majority of large river systems, flow is regulated and/or otherwise affected by operational and management activities, such as ship locking. The effect of lock operation on sediment-water oxygen fluxes was studied within a 12.9 km long impoundment at the Saar River (Germany) using eddy-correlation flux measurements. The continuous observations cover a time period of nearly 5 days and 39 individual locking events. Ship locking is associated with the generation of surges propagating back and forth through the impoundment which causes strong variations of near-bed current velocity and turbulence. These wave-induced flow variations cause variations in sediment-water oxygen fluxes. While the mean flux during time periods without lock operation was 0.5 ± 0.1 g m−2 d−1, it increased by about a factor of 2 to 1.0 ± 0.5 g m−2 d−1 within time periods with ship locking. Following the daily schedule of lock operations, fluxes are predominantly enhanced during daytime and follow a pronounced diurnal rhythm. The driving force for the increased flux is the enhancement of diffusive transport across the sediment-water interface by bottom-boundary layer turbulence and perhaps resuspension. Additional means by which the oxygen budget of the impoundment is affected by lock-induced flow variations are discussed.

Description: Existing models of spatial relations do not consider that different concepts exist on different levels in a hierarchy and in turn that the spatial relations in a given scene are a function of the specific concepts considered. One approach to determining the existence of a particular spatial relation is to compute the corresponding high level concepts explicitly using map generalization before inferring the existence of the spatial relation in question. We explore this idea through the development of a model of the spatial relation “enters” that may exist between a road and a housing estate.

Description: Model errors are inevitable in any prediction exercise. One approach that is currently gaining attention in reducing model errors is by combining multiple models to develop improved predictions. The rationale behind this approach primarily lies on the premise that optimal weights could be derived for each model so that the developed multimodel predictions will result in improved predictions. A new dynamic approach (MM-1) to combine multiple hydrological models by evaluating their performance/skill contingent on the predictor state is proposed. We combine two hydrological models, “abcd” model and variable infiltration capacity (VIC) model, to develop multimodel streamflow predictions. To quantify precisely under what conditions the multimodel combination results in improved predictions, we compare multimodel scheme MM-1 with optimal model combination scheme (MM-O) by employing them in predicting the streamflow generated from a known hydrologic model (abcd model or VIC model) with heteroscedastic error variance as well as from a hydrologic model that exhibits different structure than that of the candidate models (i.e., “abcd” model or VIC model). Results from the study show that streamflow estimated from single models performed better than multimodels under almost no measurement error. However, under increased measurement errors and model structural misspecification, both multimodel schemes (MM-1 and MM-O) consistently performed better than the single model prediction. Overall, MM-1 performs better than MM-O in predicting the monthly flow values as well as in predicting extreme monthly flows. Comparison of the weights obtained from each candidate model reveals that as measurement errors increase, MM-1 assigns weights equally for all the models, whereas MM-O assigns higher weights for always the best-performing candidate model under the calibration period. Applying the multimodel algorithms for predicting streamflows over four different sites revealed that MM-1 performs better than all single models and optimal model combination scheme, MM-O, in predicting the monthly flows as well as the flows during wetter months.

Description: The potential for riverine drinking source water to become contaminated with pathogens is related to the production and transport of fecal waste from within the local catchment area. Identifying specific relationships between land-use types and fecal contamination in riverine water provides an indication of the risk associated with land-use change and helps to target mitigation measures toward land-use types of concern. Fecal coliform (FC) data from 42 riverine sites across British Columbia (BC), Canada, were examined in relation to land-use composition (including 16 land-use types) in the local catchment area. FC concentration significantly increased in relation to anthropogenic land-use impacts but was negatively associated with undisturbed and high-elevation land types. Regression tree analysis identified that highest FC concentrations occurred in catchments characterized by more than 12.5% agricultural land and more than 1.6% urban land. Furthermore, the risk of violation of the BC partial treatment raw drinking water quality guideline for FC concentration (100 CFU 100 mL−1) increased in relation to agricultural impacts. Additional factors, such as sewage treatment discharge, low dilution in smaller streams, and higher temperatures, were associated with higher FC concentration among sites with similar levels of agricultural development. These results identify land-use types that present the greatest threat to riverine contamination, namely agricultural and urban land, and indicate the proportion of such land use associated with high contamination. Land use should be managed and source water protection should be targeted in light of these results so as to minimize the risk of surface water exposure to fecal contaminants.

Description: We report on the calibration of the one-dimensional hydrodynamic lake model DYRESM to simulate the water temperature conditions of the pre-alpine Lake Ammersee (South-east Germany) which is representative of deep and large lakes in this region. Special focus is given to the calibration in order to reproduce the correct thermal distribution and stratification including the time of onset and duration of summer stratification. To ensure the application of the model to investigate the impact of climate change on lakes, an analysis of the model sensitivity under stepwise modification of meteorological input parameters (air temperature, wind speed, precipitation, global radiation, cloud cover, vapor pressure, and tributary water temperature) was conducted. The total mean error of the calibration results is –0.23 °C, the root mean square error amounts to 1.012 °C. All characteristics of the annual stratification cycle were reproduced accurately by the model. Additionally, the simulated deviations for all applied modifications of the input parameters for the sensitivity analysis can be differentiated in the high temporal resolution of monthly values for each specific depth. The smallest applied alteration to each modified input parameter caused a maximum deviation in the simulation results of at least 0.26 °C. The most sensitive reactions of the model can be observed through modifications of the input parameters air temperature and wind speed. Hence the results show that further investigations at Lake Ammersee, such as coupling the hydrodynamic model with chemo-dynamic models to assess the impact of changing climate on biochemical conditions within lakes, can be carried out using DYRESM. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In semi-arid and arid river basins, understanding the connectivity between rivers and alluvial aquifers is one of the key challenges for the management of groundwater resources. The type of connection present (gaining, losing-connected, transitional and losing-disconnected) was assessed at 12 sites along six Murray-Darling Basin river reaches. The assessments were made by measuring the hydraulic head in the riparian zone near the rivers to evaluate if the water tables intersected the riverbeds and by measuring fluid pressure (ψ) in the riverbeds. The rationale for the latter was that ψ will always be greater than or equal to zero under connected conditions (either losing or gaining) and always lesser than or equal to zero under losing-disconnected conditions. A mixture of losing-disconnected, losing-connected and gaining conditions was found among the 12 sites. The losing-disconnected sites all had a riverbed with a lower hydraulic conductivity than the underlying aquifer, usually in the form of a silty clay or clay unit 0.5 – 2 m in thickness. The riparian water tables were 6 to 25 m below riverbed level at the losing-disconnected sites but never lower than 1 m below riverbed level at the losing-connected ones. The contrast in water table depth between connected and disconnected sites was attributed to the conditions at the time of the study, when a severe regional drought had generated a widespread decline in regional water tables. This decline was apparently compensated near losing-connected rivers by increased infiltration rates, while the decline could not be compensated at the losing-disconnected rivers because the infiltration rates were already maximal there. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Abstact The distribution of groundwater inflows in a stream reach plays a major role in controlling the stream temperature, a vital component shaping the riverine ecosystem. In this study, the DTS system was installed in a small Danish lowland stream, Elverdamsåen, to assess the seasonal dynamics of groundwater inflow zones using high spatial (1 m) and temporal (3 minutes) resolution of water temperature measurements. Four simple criteria consisting of 30 min average temperature at 16:00, mean and standard deviation of diurnal temperatures, and the day-night temperature difference were applied to three DTS datasets representing stream temperature responses to the variable meteorological and hydrological conditions prevailing in summer, winter and spring. The standard deviation criterion was useful to identify groundwater discharge zones in summer and spring conditions, while the mean temperature criterion was better for the winter conditions. In total, 20 interactions were identified from the DTS datasets representing summer, 16 in winter and 19 in spring, albeit with only two interactions contributing in all three seasons. Higher baseflow to streamflow ratio, antecedent precipitation and presence of fractured clayey till in the stream reach were deemed as the vital factors causing apparent seasonal variation in the locations of upwelling zones, prompting use of DTS not only in preconceived scenarios of large diurnal temperature change but rather a long term deployment covering variable meteorological and hydrological scenarios. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Non-point source pollution is a key issue in integrated river basin management around the world, and has resulted in water contamination, aquatic ecology deterioration and eutrophication. Xin'anjiang Catchment is the key drinking water source area for Hangzhou City, China. A promising model (Soil and Water Assessment Tool: SWAT) was applied to assess the non-point source pollution and its effect on drinking water. Sensitivity analysis of model parameters was carried out using the SUFI-2 sensitivity technique. Water discharge, sediment, total nitrogen and total phosphorus loads processes from 2000 to 2010 were simulated and the spatial distributions of non-point source pollutants were evaluated at the catchment and administrative country levels. The results show that the hydrological parameters of SWAT were dominantly sensitive for non-point source pollution simulation, including CN2, RCHRG_DP, ALPHA_BF, SOL_AWC, ESCO and SOL_K, and the characteristic parameters of sub-basins (viz. HRU_SLP, SLSUBBSN). Also, water quality parameters (viz. CH_EROD, NPERCO, RSDCO and PPERCO, PHOSKD, etc.) have a significant effect on nutrients. The model performance was very satisfactory, especially for runoff, sediment, and TP simulation. The non-point source pollutant load increased from 2001 to 2010 in the whole catchment. Total nitrogen (TN) load increased from 3,428 tons (0.59 ton/km 2 ) to 7,315 tons (1.25 ton/km 2 ) and total phosphorus (TP) load increased from 299 tons (0.05 ton/km 2 ) to 867 tons (0.15 ton/km 2 ). The contribution of rice land was the largest, accounting for nearly 95%, followed by tea garden (3.56%), winter wheat (1.37%), forest (0.07%) and grassland (0.02%). Moreover, She County and Xiuning County contributed more than half of the non-point source pollutants. This study was expected to provide a method and reference for non-point source pollution quantification, and to support water quality management implementation in China. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The increased socio-economic relevance of flood risk assessment has led to the development of innovative methodologies for the hydraulic simulation of river and floodplain systems, and has promoted the development of new techniques for flood hazard and inundation mapping (e.g. Di Baldassarre et al., 2010; Vorogushyn et al., 2010). In particular, one-dimensional (1D) and two-dimensional (2D) hydraulic models have been used more and more as numerical tools (e.g. Aronica et al., 2002; Hesselink et al., 2003; Horritt et al., 2007; Pappenberger et al., 2005) as these models have proven to be able to effectively simulate river hydraulics and floodplain inundation at different levels of detail (e.g. Horritt & Bates, 2001, 2002). Flood inundation models appear also to be useful tools for the reconstruction and analysis of historical events (e.g. Di Baldassarre et al., 2009; Horritt et al., 2010), which can be very important to provide a comprehensive assessment of exposure to floods and to develop flood risk management plans as required by the recent Floods Directive 2007/60/EC (European Commission, 2007). Copyright © 2012 John Wiley & Sons, Ltd.

Description: In the absence of model deficiencies, simulation results at the correct parameter values lead to an unbiased description of observed data with remaining deviations due to observation errors only. However, this ideal cannot be reached in the practice of environmental modeling, because the required simplified representation of the complex reality by the model and errors in model input lead to errors that are reflected in biased model output. This leads to two related problems: First, ignoring bias of output in the statistical model description leads to bias in parameter estimates, model predictions and, in particular, in the quantification of their uncertainty. Second, as there is no objective choice of how much bias to accept in which output variable, it is not possible to design an “objective” model calibration procedure. The first of these problems has been addressed by introducing a statistical (Bayesian) description of bias, the second by suggesting the use of multiobjective calibration techniques that cannot easily be used for uncertainty analysis. We merge the ideas of these two approaches by using the prior of the statistical bias description to quantify the importance of multiple calibration objectives. This leads to probabilistic inference and prediction while still taking multiple calibration objectives into account. The ideas and technical details of the suggested approach are outlined and a didactical example as well as an application to environmental data are provided to demonstrate its practical feasibility and computational efficiency.

Description: The hydraulics of step-pool streams are characterized by rapidly varied flow at the step crest, a hydraulic jump, and gradually varied flow in the pool unit of the step-pool sequence. The flow characteristics at the step crests act as the hydraulic control for the water surface profile within the upstream pool unit. Using both field and flume investigations, we demonstrate the use of weir flow concepts for assessing and categorizing the hydraulic characteristics of natural step-crests in step-pool streams. We categorize the results of our investigations in terms of the crest-clast, planform, longitudinal, and instream wood geometries of the step crests. The broad-crested weir equation can be expressed as Q = C* g0.5Wh3/2, where Q is the flowrate, C* is a dimensionless discharge coefficient, W is the crest width, g is the acceleration of gravity, and h is the upstream flow depth above the step crest. Although the flow over a natural step is generally more complex than for an engineered weir, the results of our investigations indicate that the C*-value for simulated and natural steps increases linearly as a function of the upstream head (h), with C* values ranging from 0.15 to 0.97. As a result, the application of weir flow concepts to natural steps provides means for (1) indirectly estimating flow rates; (2) characterizing the hydraulics for individual steps; (3) defining external and/or internal boundary conditions at step crests for hydraulic model simulations of natural or restored step-pool streams; and (4) estimating the upstream pressure force acting on step-crest clasts.

Description: Surface ground-penetrating radar (GPR) techniques have been used by a number of previous researchers to characterize soil moisture content in the vadose zone. However, limited temporal sampling and low resolution near the surface in these studies greatly impedes the quantitative analysis of vertical soil moisture distribution and its associated dynamics within the shallow subsurface. To further examine the capacity of surface GPR, we have undertaken an extensive 26 month field study using concurrent high-frequency (i.e., 900 MHz) reflection profiling and common-midpoint (CMP) soundings to quantitatively monitor soil moisture distribution and dynamics within the shallow vadose zone. This unprecedented data set allowed us to assess the concurrent use of these techniques over two contrasting annual cycles of soil conditions. Reflection profiles provided high-resolution traveltime data between four stratigraphic reflection events while cumulative results of the CMP sounding data set produced precise depth estimates for those reflecting interfaces, which were used to convert interval-traveltime data into soil moisture. The downward propagation of major infiltration episodes associated with seasonal and transient events are well resolved by the GPR data. The use of CMP soundings permitted the determination of direct ground wave velocities, which provided high-resolution information along the air-soil interface. This improved resolution enabled better characterization of short-duration wetting/drying and freezing/thawing processes, and permitted better evaluation of the nature of the coupling between shallow and deep moisture conditions. The nature of transient infiltration pulses, evapotranspiration episodes, and deep drainage patterns observed in the GPR data series were further examined by comparing them with a vertical soil moisture flow simulation based on the variably saturated model, HYDRUS-1D. Using laboratory-derived soil hydraulic property information from soil samples and a number of simplifying assumptions about the upper and lower-boundary condition, we were able to achieve very good agreement between measured and simulated soil moisture profiles without model calibration; this is a strong indication of the overall quality of the GPR-derived soil moisture estimates. The only notable difference between simulated values and GPR water content estimates occurred during extended dry soil conditions near the surface.

Description: Typhoon rainfall characteristics over a mesoscale mountainous watershed (drainage area of 620 km2) located in eastern Taiwan were analyzed to fill the gaps in our knowledge concerning the linkage between typhoon track, rainfall patterns, and flood peak time. This study used spatially high-resolution radar-derived rainfall estimates from 38 storm events (∼2800 h) to investigate this linkage. The effect of spatial rainfall patterns on the timing of flood peak for the selected events was examined with the aid of a diffusive wave model. The results show that the typhoon rainfall was spatially aggregated and that the relative variations in the rainfall became smaller at higher rainfall rates. The maximum hourly rainfall was approximately twice the areal mean rainfall. Three major rainfall types were identified statistically, and different typhoon tracks appeared to have preferable rainfall types. This finding is presumably due to the interaction of the typhoon circulation and precipitation with the mountainous landscape. Flood lead times were derived for the different rainfall types, and it was found that differences in their lead times could be as large as ∼3 h over the studied mesoscale watershed. It is recommended that this empirical approach be incorporated into flood forecasting and warning systems.

Description: Quantification of precipitation extremes is important for flood planning purposes, and a common measure of extreme events is the T year return level. Extreme precipitation depths in Belgium are analyzed for accumulation durations ranging from 10 min to 30 days. Spatial generalized extreme value (GEV) models are presented by considering multisite data and relating GEV parameters to geographical/climatological covariates through a common regression relationship. Methods of combining data from several sites are in common use, and in such cases, there is likely to be nonnegligible intersite dependence. However, parameter estimation in GEV models is generally done with the maximum likelihood estimation method (MLE) that assumes independence. Estimates of uncertainty are adjusted for spatial dependence using methodologies proposed earlier. Consistency of GEV distributions for various durations is obtained by fitting a smooth function to the preliminary estimations of the shape parameter. Model quality has been assessed by various statistical tests and indicates the relevance of our approach. In addition, a methodology is applied to account for the fact that measurements have been made in fixed intervals (usually 09:00 UTC–09:00 UTC). The distribution of the annual sliding 24 h maxima was specified through extremal indices of a more than 110 year time series of 24 h aggregated 10 min rainfall and daily rainfall. Finally, the selected models are used for producing maps of precipitation return levels.

Description: During the Asian monsoon period, intense and precipitation commonly occurs for an extended period in accompaniment with a reduction in solar radiation. This suggests that wet surface evapotranspiration is important contributor to the total evapotranspiration. Therefore, investigating evapotranspiration over a wet canopy surface is critical to achieve a better understanding of water and energy cycles in Asia. In this study, we estimated surface resistances under wet conditions in a mixed forest influenced by the East Asian monsoon system. We showed that the surface resistance had a non-negligible magnitude of about 30  sm −1 even under wet conditions. We also found that the ratio between the actual and potential evapotranspiration depended on the friction velocity regardless of the time of day. Our analyses suggest that this dependency is tightly related to the underestimation of turbulent fluxes by the eddy-covariance system under wet surface conditions. Together, our findings suggest that the wet surface resistance, although small, should be considered in simulating evapotranspiration because the forest ecosystem is strongly coupled to the overlying atmosphere. This could significantly improve the shortcomings of evapotranspiration measurement in Asian forest canopies influenced by the monsoon system. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Forecast ensembles of hydrological and hydrometeorologial variables are prone to various uncertainties arising from climatology, model structure and parameters, and initial conditions at the forecast date. Postprocessing methods are usually applied to adjust the mean and variance of the ensemble without any knowledge about the uncertainty sources. This study initially addresses the drawbacks of a commonly used statistical technique, Quantile Mapping (QM), in bias correction of hydrologic forecasts. Then, an auxiliary variable, the failure index (γ) is proposed to estimate the ineffectiveness of the postprocessing method based on the agreement of adjusted forecasts with corresponding observations during an analysis period prior to the forecast date. An alternative postprocessor based on copula functions is then introduced such that marginal distributions of observations and model simulations are combined to create a multivariate joint distribution. A set of 2500 hypothetical forecast ensembles with parametric marginal distributions of simulated and observed variables are postprocessed with both QM and the proposed multivariate postprocessor. Deterministic forecast skills show that the proposed copula-based postprocessing is more effective than the QM method in improving the forecasts. It is found that the performance of QM is highly correlated with the failure index, unlike the multivariate postprocessor. In probabilistic metrics, the proposed multivariate postprocessor generally outperforms QM. Further evaluation of techniques is conducted for river flow forecast of Sprague River Basin in southern Oregon. Results show that the multivariate postprocessor performs better than the QM technique; it reduces the ensemble spread and is a more reliable approach for improving the forecast. Copyright © 2012 John Wiley & Sons, Ltd.

Description: In reservoirs used for geologic CO2 sequestration, brine films remaining on mineral surfaces can influence flow, diffusion, and reactions. We have investigated how the capillary (disjoining) potential influences the thickness of a KCsI2 brine film on both smooth and rough SiO2 surfaces [root mean square roughness (Rrms), 1.6 and 330 nm, respectively], under confinement with supercritical (sc) CO2. The thicknesses of brine films coating interior surfaces of SiO2 windows in a high-pressure cell were determined through synchrotron X-ray fluorescence of two tracer ions (I− and Cs+) at 7.8 MPa and 40°C (representative of conditions at about 0.78 km below the land surface), with scCO2 as the immiscible confining fluid. The measured area-averaged film thicknesses on the 330 nm Rrms silica surface ranged from 265 to 249 nm for capillary potentials measured within a narrow range from 0.18 to 3.7 kPa. Over this same range of potentials, film thicknesses measured on the smooth (1.6 nm Rrms) silica surface were about 2 nm, although equilibrium does not appear to have been achieved. The measured average brine film thicknesses were strongly controlled by surface roughness, with very weak variation in response to the fairly narrow range of tested capillary potentials.

Description: A well-established precept in forest hydrology is that any reduction of forest cover will always have a progressively smaller effect on floods with increasing return period. The underlying logic in snow environments is that during the largest snowmelt events the soils and vegetation canopy have little additional storage capacity and under these conditions much of the snowmelt will be converted to runoff regardless of the amount or type of vegetation cover. Here we show how this preconceived physical understanding, reinforced by the outcomes of numerous paired watershed studies, is indefensible because it is rationalized outside the flood frequency distribution framework. We conduct a meta-analysis of postharvest data at four catchments (3–37 km2) with moderate level of harvesting (33%–40%) to demonstrate how harvesting increases the magnitude and frequency of all floods on record (19–99 years) and how such effects can increase unchecked with increasing return period as a consequence of changes to both the mean (+11% to +35%) and standard deviation (−12% to +19%) of the flood frequency distribution. We illustrate how forest harvesting has substantially increased the frequency of the largest floods in all study sites regardless of record length and this also runs counter to the prevailing wisdom in hydrological science. The dominant process responsible for these newly emerging insights is the increase in net radiation associated with the conversion from longwave-dominated snowmelt beneath the canopy to shortwave-dominated snowmelt in harvested areas, further amplified or mitigated by basin characteristics such as aspect distribution, elevation range, slope gradient, amount of alpine area, canopy closure, and drainage density. Investigating first order environmental controls on flood frequency distributions, a standard research method in stochastic hydrology, represents a paradigm shift in the way harvesting effects are physically explained and quantified in forest hydrology literature.

Description: Surface transient storage (STS) has functional significance in stream ecosystems because it increases solute interaction with sediments. After volume, mean residence time is the most important metric of STS, but it is unclear how this can be measured accurately or related to other timescales and field-measureable parameters. We studied mean residence time of lateral STS in small streams over Reynolds numbers (Re) 5000–200,000 and STS width to length (W/L) aspect ratios between 0.2–0.75. Lateral STS have flow fields characterized by a shear layer spanning the length of the STS entrance, and one primary gyre and one or more secondary gyre(s) in the STS. The study's purpose was to define, measure, and compare residence timescales: volume to discharge ratio (Langmuir timescale); area under normalized concentration curve; and characteristic time of exponential decay, and to compare these timescales to field measureable parameters. The best estimate of STS mean residence time—primary gyre residence time—was determined to be the first characteristic time of exponential decay. An apparent mean residence time can arise, which is considerably larger than other timescales, if probes are placed within secondary gyre(s). The Langmuir timescale is the minimum mean residence time, and is linearly correlated to channel velocity and STS width. The lateral STS mean residence time can be predicted using a physically based hydromorphic timescale derived by Uijttewaal et al. (2001) with an entrainment coefficient of 0.031 ± 0.009 for the Re and W/L studied.

Description: Hydraulic stimulation of subsurface rocks is performed in developing geothermal and hydrocarbon reservoirs to create permeable zones and enhance flow and transport in low-permeability formations. Borehole fluid injection often induces measurable microearthquakes (MEQs). While the nature and source of the processes that lead to triggering of these events is yet to be fully understood, a major hypothesis has linked these events to an increase in pore pressure that decreases the effective compressional stress and causes sliding along preexisting cracks. Based on this hypothesis, the distribution of the resulting microseismicity clouds can be viewed as monitoring data that carry important information about the spatial distribution of hydraulic rock properties. However, integration of fluid-induced microseismicity events into prior rock permeability distributions is complicated by the discrete nature of the MEQ events, which is not amenable to well-established inversion methods. We use kernel density estimation to first interpret the MEQ data events as continuous seismicity density measurements and, subsequently, assimilate them to estimate rock permeability distribution. We apply the ensemble Kalman filter (EnKF) for microseimic data integration where we update a prior ensemble of permeability distributions to obtain a new set of calibrated models for prediction. The EnKF offers several advantages for this application, including the ensemble formulation for uncertainty assessment, convenient gradient-free implementation, and the flexibility to incorporate various failure mechanisms and additional data types. Using several numerical experiments, we illustrate the suitability of the proposed approach for characterization of reservoir hydraulic properties from discrete MEQ monitoring measurements.

Description: We propose a new runoff model including an outflow process that was applied to two adjacent basins (CL, TL) located in Lambir Hills National Park in north-central Sarawak, Malaysia. Rainfall, runoff, topography, and soil layer thickness were observed. About 19 % of annual runoff was observed in the CL basin (21.97 ha), whereas about 46 % was observed in the TL basin (23.25 ha). It was inferred that the CL basin has an outflow because of low base flow, small runoff peak, and excessive water loss. By incorporating the outflow process into the HYdrological CYcle MODEL (HYCYMODEL), good agreement between the data generated by the model and that observed was shown, with the exception of the data from the rainless period. Then, the fitting parameters for each basin were exchanged, except for the outflow parameter, and the characteristics of each basin were compared by calculating virtual runoff. As a result, the low base flow of the CL basin was estimated by the movement of the rainwater that escaped from the basin as deep percolation or lateral flow (11 % of rainfall). The potential of the CL basin for mitigating flood and drought appeared to be higher than that of the TL basin. This is consistent with the topographic characteristics of the CL basin, which has a gentler slope than the TL basin. Copyright © 2012 John Wiley & Sons, Ltd.

Description: For many basins, identifying changes to water quality over time and understanding current hydrologic processes are hindered by fragmented and discontinuous water-quality and hydrology data. In the coal mined region of the New River basin and Indian Fork sub-basin, muted and pronounced changes, respectively, to concentration-discharge relationships were identified using linear regression on log-transformed historical (1970s-1980s) and recent (2000s) water-quality and streamflow data. Changes to concentration-discharge relationships were related to coal mining histories and shifts in land use. Hysteresis plots of individual storms from 2007 (New River) and the fall of 2009 (Indian Fork) were used to understand current hydrologic processes in the basins. In the New River, storm magnitude was found to be closely related to the reversal of loop rotation in hysteresis plots; a peak-flow threshold of 25 cubic meters per second (m 3 /s) segregates hysteresis patterns into clockwise and counterclockwise rotational groups. Small storms with peak flow less than 25 m 3 /s often resulted in dilution of constituent concentrations in headwater tributaries like Indian Fork and concentration of constituents downstream in the mainstem of the New River. Conceptual two or three component mixing models for the basins were used to infer the influence of water derived from spoil material on water quality. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Thus far, measurements and estimations of actual evapotranspiration (ET) from high-altitude grassland ecosystems in remote areas like the Qinghai-Tibetan plateau are still insufficient. To address these issues, a comparison between the results of the eddy covariance (EC) measurements and the estimates, considering the Katerji and Perrier (KP), the Todorovic (TD) and the Priestley-Taylor (PT) models, was carried out over an alpine grassland (38 o 03’1.7'' N, 100 o 27’ 26'' E; 3032 m a.s.l.) during the growing seasons in 2008 and 2009. The results indicated that the KP model after a particularly simple calibration gave the most effective ET values in different time scales, the PT model slightly underestimate ET at night, and the TD model significantly overestimated ET at noon. In addition, the canopy resistance calculated by the TD model was completely different from that calculated using the inverted EC-measured data and the KP model, which may be due to some unrealistic assumptions made by the TD model. The KP parameters were a  = 0.17 and b  = 1.50 for the alpine grassland, and appeared to be interannually stable. However, the PT parameter showed some interannual variations ( α  = 0.83 and 0.74 for 2008 and 2009, respectively). Therefore, the KP model was preferred to estimate the actual ET at both hourly and daily time scales. The PT model, being the simplest approach and field condition dependent, was recommended when available weather data were rare. On the contrary, the TD model always overestimated the actual ET and should be avoided in case of the alpine grassland ecosystems. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Stochastic weather generators have evolved as tools for creating long time-series of synthetic meteorological data at a site for risk assessments in hydrologic and agricultural applications. Recently, their use has been extended as downscaling tools for climate change impact assessments. Non-parametric weather generators, which typically use a K-Nearest Neighbour (K-NN) resampling approach, require no statistical assumptions about probability distributions of variables and can be easily applied for multisite use. Two characteristics of traditional K-NN models result from resampling daily values: (a) temporal correlation structure of daily temperatures may be lost, and (b) no values lying below or exceeding historical observations can be simulated. Temporal correlation in simulated temperature data is important for hydrologic applications. Temperature is a major driver of many processes within the hydrologic cycle (for example, evaporation, snow melt, etc) that may affect flood levels. As such, a new methodology for simulation of climate data using the K-NN approach is presented (named KnnCAD Version 4). A block resampling scheme is introduced along with perturbation for the reshuffled daily temperature data to create 675 years of synthetic historical daily temperatures for the Upper Thames River basin in Ontario, Canada. The updated KnnCAD model is shown to adequately reproduce observed monthly temperature characteristics, temporal and spatial correlations while simulating reasonable values outside the range of observations. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Global change is predicted to increase temperature substantially in the North as well as altering run-off regimes with less synchronicity as the importance of snow melt declines. River biota and ecosystem processes will be influenced across all levels of organization, both in concert and individually. It is of vital importance that the impacts, and their likely magnitude, can be identified in order to deploy suitable adaptation strategies at the catchment scale. In this paper, we re-analyse 4 data sets from studies conducted in Greenland (66-69 o N), Iceland (64 o N), Sweden (60 o N) and Denmark (55-57 o N) to try and tease out the likely impacts of water temperature and hydrology in shaping the stream communities and ecosystem processes in high-latitude catchments. Water temperature was the environmental variable that best explained macroinvertebrate community composition across latitudes. In contrast, no significant relationship between macroinvertebrate community composition and measures of hydraulic stability (or nutrients) was found. We found a strong linear relationship between decay rate of leaf litter and water temperature (r 2  = 0.68; p 〈 0.0001) independent of latitudes. Our study suggests that temperature could be the primary driver of ecosystem change in future with northern catchments likely to be especially vulnerable. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The velocity field in a river flow cross-sectional area can be determined by applying entropy as done by Chiu (1987), who developed a two-dimensional model of flow velocity based on the knowledge of maximum velocity, u max , and the dimensionless entropic parameter, characteristic of the river site. This is appealing in the context of discharge monitoring, particularly for high floods, considering that u max occurs in the upper portion of flow area and can be easily sampled, unlike velocity in the lower portion of flow area. The simplified form of Chiu's entropy-based velocity model, proposed by Moramarco et al. (2004), has been found to be reasonably accurate for determining mean flow velocity along each vertical sampled in the flow area, but no uncertainty analysis has been reported for this simplified entropy-based velocity model. This study, therefore, performed uncertainty analysis of the simplified model following a procedure proposed by Misirli et al. (2003). The flow velocity measurements at the Rosciano River section along the Chiascio River, central Italy, carried out for a period spanning 20 years were used for this purpose. Results showed that the simplified entropy velocity model was able to provide satisfactory estimates of velocity profiles in the whole flow area and the 95% confidence bands for the computed estimated mean vertical velocity were quite representative of observed values. In addition, using these 95% confidence bands, it was possible to have an indication of the uncertainty in the determination of mean cross-sectional flow velocity as well. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Total evaporation ( ET ) is one of the major components of the water budget of a wetland. Very little research has been conducted on the loss of water to the atmosphere from different wetland vegetation types occurring in southern Africa. This study on the ET of taro (locally known as madumbe) and sedge within the Mbongolwane wetland was conducted to assess the potential impact of madumbe cultivation on the hydrology of the wetland. Sugarcane planted on the contributing catchment outside the wetland was the other crop examined. Two field campaigns were conducted in November 2009 and January 2010 during the growing season of the madumbe crop to quantify ET rates in the Mbongolwane wetland and from sugar cane in the surrounding catchment. Total evaporation was measured over two vegetation types in the wetland, namely: madumbe ( Colocasia esculenta ); sedge ( Cyperus latifolius ) with some reeds ( Phragmites australis ); and sugarcane in adjacent terrestrial areas. Total evaporation from the madumbes ranged from 1.0 to 6.0 mm day -1 . The daily average ET rates in November 2009 were 3.5 and 4.9 mm for the madumbe and sedge sites respectively and 4.0 mm for sugarcane grown in the catchment. The daily average ET rates in January 2010 were 3.3 and 3.7 mm for the madumbes and sedge sites, respectively and 2.4 mm for the sugarcane site. The daily ET was therefore lower at the madumbe site in November 2009 and in January 2010 compared to the sedge site. An average crop factor ( Kc ) of 0.6 was obtained from this study during the growth stage of the madumbes. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Much of what is known about groundwater circulation and geochemical evolution in carbonate platforms is based on platforms that are fully confined or unconfined. Much less is known about groundwater flow paths and geochemical evolution in partially confined platforms, particularly those supporting surface water. In north-central Florida, sea level rise and a transition to a wetter climate during the Holocene formed rivers in unconfined portions of the Florida carbonate platform. Focusing on data from the Santa Fe River basin, we show river formation has led to important differences in the hydrological and geochemical evolution of the Santa Fe River basin relative to fully confined or unconfined platforms. Runoff from the siliciclastic confining layer drove river incision and created topographic relief, reorienting the termination of local and regional groundwater flow paths from the coast to the rivers in unconfined portions of the platform. The most chemically evolved groundwater occurs at the end of the longest and deepest flow paths, which discharge near the center of the platform because of incision of the Santa Fe River at the edge of the confining unit. This pattern of discharge of mineralized water differs from fully confined or unconfined platforms where discharge of the most mineralized water occurs at the coast. Mineralized water flowing into the Santa Fe River is diluted by less evolved water derived from shorter, shallower flow paths that discharge to the river downstream. Formation of rivers shortens flow path lengths, thereby decreasing groundwater residence times and allowing freshwater to discharge more quickly to the oceans in the newly formed rivers than in platforms that lack rivers. Similar dynamic changes to groundwater systems should be expected to occur in the future as climate change and sea level rise develop surface water on other carbonate platforms and low lying coastal aquifer systems. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Sea ice dynamic and thermodynamic processes are important and highly variable elements of the marginal ice zone (MIZ). This study examines the detection and classification of statistically separable sea ice classes in the MIZ through a range of temporal and spatial scales. A helicopter-based laser system was used to obtain large-scale and a ship-based laser profiler to identify small-scale roughness types respectively. The analysis of variance (ANOVA) of surface height data from helicopter- and ship-based laser systems, active microwave (AMW) C-band backscattering data and passive microwave (PMW) (37 and 89 GHz) brightness temperature data reveal different classes that statistically differ from one another. We found significant statistical difference in variances in AMW data with six classes that differ in VV polarization, three classes in VH polarization and five classes in HH polarization in the MIZ (e.g. snow-covered first-year ice, ice rubble, pancake ice, frost flowers, melt pond, flooded ice, and ice edge) of southeastern Beaufort Sea. The PMW emission was not as effective at discrimination, yielding only one statistically separable class. The results can potentially be extended to satellite-based investigations of the MIZ at regional scales. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The USLE/RUSLE model was designed to predict long term (~20 years) average annual soil loss by accounting for the effects of climate, soil, topography and crops. The USLE/RUSLE model operates mathematically in two steps. The first step involves the prediction of soil loss from the “unit” plot, a bare fallow area 22.1 m long on a 9 % slope gradient with cultivation up and down the slope. Appropriate values of the factors accounting for slope length, gradient, crops and crop management and soil conservation practice are then used to adjust that soil loss to predict soil loss from areas that have conditions that are different from the unit plot. Replacing EI 30 , the USLE/RUSLE event erosivity index, by the product of the runoff ratio ( Q R ) and EI 30 can enhance the capacity of the model to predict short term soil loss from the unit plot if appropriate data on runoff is available. Replacing the EI 30 index by another index has consequences on other factors in the model. The USLE/RUSLE soil erodibility factor cannot be used when the erosivity factor is based on Q R EI 30 . Also, the USLE/RUSLE factors for slope length, slope gradient crops and crop management, and soil conservation practice cannot be used when runoff from other than the unit plot is used to calculate Q R . Here equations are provided to convert the USLE/RUSLE factors to values suitable for use when the erosivity factor is based on the Q R EI 30 index under these circumstances. At some geographic locations, non linear relationships exist between soil loss from bare fallow areas and the Q R EI 30 index. The effect of this on the slope length factor associated with the Q R EI 30 index is demonstrated using data from runoff and soil loss plots located at the Sparacia site, Sicily. Copyright © 2012 John Wiley & Sons, Ltd.

Description: A raster based glacier sub-model was successfully introduced in the distributed hydrological model FEST-WB to simulate the water balance and surface runoff of large Alpine catchments. The glacier model is based on temperature-index approach for melt, on linear reservoir for melt water propagation into the ice and on mass balance for accumulation; the initialization of the volume of ice on the basin was based on a formulation depending on surface topography. The model was first tested on a sub-basin of the Rhone basin (Switzerland), which is for 62% glaciated; the calibration and validation were based on comparison between simulated and observed discharge from 1999 to 2008. The model proved to be suitable to simulate the typical discharge seasonality of a heavily glaciated basin. The performance of the model was also tested by simulating discharge in the whole Swiss Rhone basin, in which glaciers contribution is not negligible, in fact in summer about the 40% of the discharge is due to glacier melt. The model allowed to take into account the volume of water coming from glaciers melt and its simple structure is suitable for analysis of the effects of climate change on hydrological regime of high mountain basins, with available meteorological forcing from current RCM. Copyright © 2012 John Wiley & Sons, Ltd.

Description: As part of this special issue on Hydrological Ensemble Prediction Systems (HEPS), this paper reports on the intercomparison experiment for post-processing techniques that has been initiated in 2011 by the international community on Hydrologic Ensemble Predictions (HEPEX). The design of this inter-comparison experiment and the data sets available are presented. The post-processing methods that have been applied to date are listed and example results are shown. It is expected that through the exchange and joint verification and analysis of the post-processing results, the inter-comparison experiment will contribute to a fast improvement and applicability of post-processing techniques. Readers are invited to join the inter-comparison experiment. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Biocrusts abound in southern Israel, covering the Hallamish dune field near Nizzana (NIZ) in the Negev (mean annual precipitation of 95 mm) and the coast of Nizzanim (NIM) near Ashdod (mean annual precipitation of 500 mm). While the hydrological response of the NIZ crust to natural rain events was thoroughly investigated, no data is available on the hydrological response of the NIM crust. Runoff was monitored in runoff plots during the years 2005-2008 and in addition, sprinkling experiments were carried out on NIM and NIZ crusts. For the evaluation of the possible factors that may control runoff initiation, fine content of the parent material, crust thickness, compressional strength, hydrophobicity, surface microrelief, organic matter, biomass (chlorophyll a and total carbohydrates) and the crust's species composition of NIM were studied and compared to that of NIZ. The data showed that in comparison to the NIZ crust that readily generated runoff, no runoff was produced by the NIM crust. This was so despite the fact that (a) Microculeus vaginatus predominated in both crusts, (b) the substantially higher rain intensities in NIM (c) the greater thickness and higher chlorophyll content and (d) the lower microrelief at NIM in comparison to NIZ. The lack of runoff in NIM was explained by its low amounts of exopolysaccharides that did not suffice to affectively clog the surface and in turn to facilitate runoff initiation. The absence of runoff and its consequences on the NIM ecosystem are discussed. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The Soil Conservation Service (SCS) curve number estimates of direct runoff from rainfall for semiarid catchments can be inaccurate. Investigation of the Walnut Gulch Experimental Watershed (Southeastern Arizona) and its 10 nested catchments determined that the inaccuracy is due to the original SCS ratio ( λ ) of 0.2 between initial abstraction and maximum potential retention. Sensitivity analyses indicate that runoff estimation can be very sensitive to the initial abstraction ratio, especially for relatively low rainfall amount and for watersheds covered by deep, coarse, and porous soil, conditions that dominate many semiarid watersheds worldwide. Changing the ratio of initial abstraction to the maximum potential retention to optimal values ranging from 0.01 to 0.53 for different Walnut Gulch catchments improved runoff estimates. The greater the channel area and the finer the soil, the smaller the initial abstraction ratio is. The variation of the initial abstraction ratio for the Walnut Gulch Experimental Watershed is due to the variation of maximum potential retention and initial abstraction, which are channel area and soil dependent parameters. The greater the channel area, the higher the maximum potential retention S is; and the coarser the soil, the larger the initial abstraction I a is. In addition, the effect of initial abstraction ratio on runoff estimation increases with decreasing curve number. Thus, impacts of initial abstraction ratio on runoff estimation should be considered, especially for semiarid watersheds where the curve number is usually low. Copyright © 2012 John Wiley & Sons, Ltd.

Description: The Future Midwestern Landscapes (FML) project is part of the U.S. Environmental Protection Agency's undertaken Ecosystem Services Research Program. The goal of the FML project is to quantify changes in ecosystem services across the region as a result of the growing demand for biofuels. Watershed models are an efficient way to quantify ecosystem services of water quality and quantity. By calibrating models we can better capture watershed characteristics before they are applied to make predictions. The Kaskaskia River watershed in Illinois was selected to investigate the effectiveness of different calibration strategies (single-site and multi-site calibrations) for streamflow, total suspended sediment (TSS) and total nitrogen (TN) loadings using the Soil and Water Assessment Tool (SWAT). Four USGS gauges were evaluated in this study. Single-site calibration was performed from downstream site to upstream site, and multi-site calibration was performed and fine-tuned based on the single-site calibration results. Generally, simulated streamflow and TSS were not much affected by different calibration strategies. However, when single-site calibration was performed at the most downstream site, the Nash-Sutcliffe Efficiency (NSE) values for TN ranged between -0.09 and 0.53 at the other sites; and when single-site calibration was performed at the most upstream site, the NSE values ranged between -8.38 and -0.07 for the other sites. The NSE values for TN were improved to 0.5 – 0.59 for all four sites when multi-site calibration was performed. The results of multi-site calibration and validation showed an improvement on model performance on TN and highlighted that multi-site calibrations are needed to assess the hydrological and water quality processes at various spatial scales. Copyright © 2012 John Wiley & Sons, Ltd.

Description: Snow over the Tibetan Plateau (TP) is an important water source of major Asian rivers and greatly influences water availability in the downstream areas. In this study, snow cover dynamics of the four characteristic lake basins, Cedo Caka, Selin Co, Nam Co, and Yamzhog Yumco during hydrological years 2001–2010 (September through August) are examined at the basin scale using the flexible multiday combined MODIS snow cover products. The time series of multiday, seasonal, and annual snow covered area (SCA), onset/disappearance dates of snow, snow covered days (SCD), peaks of maximum SCA, and snow cover index (SCI) for each hydrological year (HY) are examined. Results show there is no obvious trend of snow cover change in the examined period, although Nam Co basin has the greatest SCA in all four basins and in all years, and Cedo Caka and Selin Co basins show the smallest SCA in most of the years. Overall, the HY2007 shows a greater snow extent and HY2010 a smaller for the region, with exceptions for the Nam Co basin where the HY2003 is the greatest and for Cedo Caka basin where the HY2004 is the smallest. Statistical analysis between lake level changes and lake basin's SCA, precipitation and pan evaporation (ETpan) changes shows that (1) Cedo Caka's water level rise was highly correlated with the basin's SCA changes (r = 0.94, p = 0.063); (2) Selin Co's water level rise was significantly correlated with the basin's SCA, precipitation and ETpan changes (r = 0.99, p = 0.029); and (3) lake level changes of Nam Co and Yamzhog Yumco were correlated with their corresponding lake basin's SCA, precipitation and ETpan changes (r = 0.87 and r = 0.86, respectively), although insignificant at the 95% level. This could have been due to precipitation and ETpan data of a distant meteorological station for Nam Co lake basin and the complex hydrological processes in the Yamzhog Yumco basin. This study suggests that the examination of time series snow cover dynamics is important to evaluate the water budget of lake basins with snow as a major component of water balance.

Description: Recently evapotranspiration has been hypothesized to promote the secondary formation of calcium carbonate year-round on tree islands in the Everglades by influencing groundwater ions concentrations. However, the role of recharge and evapotranspiration as drivers of shallow groundwater ion accumulation has not been investigated. The goal of this study is to develop a hydrologic model that predicts the chloride concentrations of shallow tree island groundwater and to determine the influence of overlying biomass and underlying geologic material on these concentrations. Groundwater and surface water levels and chloride concentrations were monitored on eight constructed tree islands at the Loxahatchee Impoundment Landscape Assessment (LILA) from 2007–2010. The tree islands at LILA were constructed predominately of peat, or of peat and limestone, and were planted with saplings of native tree species in 2006 and 2007. The model predicted low shallow groundwater chloride concentrations when inputs of regional groundwater and evapotranspiration-to-recharge rates were elevated, while low evapotranspiration-to-recharge rates resulted in a substantial increase the chloride concentrations of the shallow groundwater. Modeling results indicated that evapotranspiration typically exceeded recharge on the older tree islands and those with a limestone lithology, which resulted in greater inputs of regional groundwater. A sensitivity analysis indicated the shallow groundwater chloride concentrations were most sensitive to alterations in specific yield during the wet season and hydraulic conductivity in the dry season. In conclusion the inputs of rainfall, underlying hydrologic properties of tree islands sediments and forest structure may explain the variation in ion concentration seen across Everglades tree islands. Copyright © 2012 John Wiley & Sons, Ltd.