Simulating pan‐Arctic runoff with a macro‐scale terrestrial water balance model

A terrestrial hydrological model, developed to simulate the high‐latitude water cycle, is described, along with comparisons with observed data across the pan‐Arctic drainage basin. Gridded fields of plant rooting depth, soil characteristics (texture, organic content), vegetation, and daily time series of precipitation and air temperature provide the primary inputs used to derive simulated runoff at a grid resolution of 25 km across the pan‐Arctic. The pan‐Arctic water balance model (P/WBM) includes a simple scheme for simulating daily changes in soil frozen and liquid water amounts, with the thaw–freeze model (TFM) driven by air temperature, modelled soil moisture content, and physiographic data. Climate time series (precipitation and air temperature) are from the National Centers for Environmental Prediction (NCEP) reanalysis project for the period 1980–2001. P/WBM‐generated maximum summer active‐layer thickness estimates differ from a set of observed data by an average of 12 cm at 27 sites in Alaska, with many of the differences within the variability (1σ) seen in field samples. Simulated long‐term annual runoffs are in the range 100 to 400 mm year^?1. The highest runoffs are found across northeastern Canada, southern Alaska, and Norway, and lower estimates are noted along the highest latitudes of the terrestrial Arctic in North America and Asia. Good agreement exists between simulated and observed long‐term seasonal (winter, spring, summer–fall) runoff to the ten Arctic sea basins (r = 0·84). Model water budgets are most sensitive to changes in precipitation and air temperature, whereas less affect is noted when other model parameters are altered. Increasing daily precipitation by 25% amplifies annual runoff by 50 to 80% for the largest Arctic drainage basins. Ignoring soil ice by eliminating the TFM sub‐model leads to runoffs that are 7 to 27% lower than the control run. The results of these model sensitivity experiments, along with other uncertainties in both observed validation data and model inputs, emphasize the need to develop improved spatial data sets of key geophysical quantities (particularly climate time series) to estimate terrestrial Arctic hydrological budgets better. Copyright © 2003 John Wiley & Sons, Ltd.     

Wetlands with controlled drainage and sub‐irrigation systems—modelling of the water balance

Over the past centuries, the agricultural use of wetlands in Central Europe has required interference with the natural wetland water balance. Often this has consisted of drainage measures alone. In low‐precipitation areas, it has also involved the operation of combined drainage and sub‐irrigation systems. Model studies conducted as part of planning processes, or with a view to finding out the impact of changing climate conditions on the water balance of wetlands, must take these facts into account. For this reason, a water balance model has been devised for wetlands whose water balance is governed by water resources management systems. It is based on the WBalMo model system. Special modules were integrated into WBalMo to calculate the water balance of wetland areas (WABI module) and to regulate inflow partitioning within the wetland (REGINF module). When calculating the water balance, the WABI module takes into account precipitation and potential evapotranspiration, groundwater levels below surface, soil types, land‐use classes, inflows via the running water system, and data for target water levels. It provides actual evapotranspiration, discharge into the running water system, and groundwater levels in the area. The example of the Spreewald, a major wetland area in north‐eastern Germany, was used to design and test the WBalMo Spreewald model. The comparison of measured and calculated water balance parameters of the wetland area confirms the suitability of the model for water balance studies in wetlands with complex water resources management systems. The results reveal the strong influence of water management on the water balance of such areas. The model system has proved to be excellently suited for planning and carrying out water management measures aimed at the sustainable development of wetlands. Furthermore, scenario analyses can be used to assess the impact of global change on the water balance of wetlands. Copyright © 2006 John Wiley & Sons, Ltd.     

Lumped surface and sub‐surface runoff for erosion modeling within a small hilly watershed in northern Vietnam

Developing models to predict on‐site soil erosion and off‐site sediment transport at the agricultural watershed scale represent an on‐going challenge in research today. This study attempts to simulate the daily discharge and sediment loss using a distributed model that combines surface and sub‐surface runoffs in a small hilly watershed (< 1 km²). The semi‐quantitative model, Predict and Localize Erosion and Runoff (PLER), integrates the Manning–Strickler equation to simulate runoff and the Griffith University Erosion System Template equation to simulate soil detachment, sediment storage and soil loss based on a map resolution of 30 m × 30 m and over a daily time interval. By using a basic input data set and only two calibration coefficients based, respectively, on water velocity and soil detachment, the PLER model is easily applicable to different agricultural scenarios. The results indicate appropriate model performance and a high correlation between measured and predicted data with both Nash–Sutcliffe efficiency (Ef) and correlation coefficient (r²) having values > 0.9. With the simple input data needs, PLER model is a useful tool for daily runoff and soil erosion modeling in small hilly watersheds in humid tropical areas. Copyright © 2013 John Wiley & Sons, Ltd.     

Gypsies in the palace: experimentalist's view on the use of 3‐D physics‐based simulation of hillslope hydrological response

As a fundamental unit of the landscape, hillslopes are studied for their retention and release of water and nutrients across a wide range of ecosystems. The understanding of these near‐surface processes is relevant to issues of runoff generation, groundwater–surface water interactions, catchment export of nutrients, dissolved organic carbon, contaminants (e.g. mercury) and ultimately surface water health. We develop a 3‐D physics‐based representation of the Panola Mountain Research Watershed experimental hillslope using the TOUGH2 sub‐surface flow and transport simulator. A recent investigation of sub‐surface flow within this experimental hillslope has generated important knowledge of threshold rainfall‐runoff response and its relation to patterns of transient water table development. This work has identified components of the 3‐D sub‐surface, such as bedrock topography, that contribute to changing connectivity in saturated zones and the generation of sub‐surface stormflow. Here, we test the ability of a 3‐D hillslope model (both calibrated and uncalibrated) to simulate forested hillslope rainfall‐runoff response and internal transient sub‐surface stormflow dynamics. We also provide a transparent illustration of physics‐based model development, issues of parameterization, examples of model rejection and usefulness of data types (e.g. runoff, mean soil moisture and transient water table depth) to the model enterprise. Our simulations show the inability of an uncalibrated model based on laboratory and field characterization of soil properties and topography to successfully simulate the integrated hydrological response or the distributed water table within the soil profile. Although not an uncommon result, the failure of the field‐based characterized model to represent system behaviour is an important challenge that continues to vex scientists at many scales. We focus our attention particularly on examining the influence of bedrock permeability, soil anisotropy and drainable porosity on the development of patterns of transient groundwater and sub‐surface flow. Internal dynamics of transient water table development prove to be essential in determining appropriate model parameterization. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaporation and energy balance of a wet grassland at Tadham Moor on the Somerset Levels

Wet grasslands are important both for their conservation value and for their important hydrological function. Evaporation is an important component of the water balance of a wetland. Where water is limited rainfall and/or surface (or sub‐surface) inflows are required to balance the summer evaporation and thus become a significant factor in the maintenance and environmental health of a wetland. This study presents an almost complete year of measurements of the water and energy balance of a wet grassland in the Somerset Levels in southwest England. The majority of the radiant energy at this site goes into evaporation. There is a strong seasonality of the controls on evaporation; the roughness length varies by a factor of 10 between winter and summer. The surface resistance to evaporation is low, close to zero, during the winter when the water table is at, or just below the surface. In the summer the water table drops to 80 cm below the surface; there is no sign of soil water stress on the evaporation but a clear effect of the senescence of grass during seed head production and of the subsequent harvest. There is clear evidence that water for evaporation is provided by the drainage ditches through sub‐surface flow. Standard evaporation formulae—such as the Penman‐Monteith equation with constant and standard parameters—provide a reasonable simulation of the total evaporation at this site, although they miss some of the seasonal detail. Copyright © 2008 John Wiley & Sons, Ltd.     

A method to detect and characterize sub‐daily flow fluctuations

Hydro peaking causes an important environmental impact on running water ecosystems. Many affected rivers have a poor ecological status. In rivers affected by hydro peaking, the flow conditions are highly complex and difficult to grasp. To develop a general framework for detecting and characterizing sub‐daily flow fluctuations, we analysed more than 500 Austrian hydrographs, covering the whole range from unimpacted to heavily impacted rivers. Different fluctuation types could be identified according to the potential source: e.g. sub‐daily flow fluctuations caused by hydro peaking, rainfall or snow and glacier melt. Additionally, the term ‘hydro fibrillation’ was established, to indicate frequently occurring artificial fluctuations with comparably low intensities. An automatic procedure was used to detect frequency and intensity of each flow fluctuation. Using variables based on duration curves of flow fluctuation rates (ramping rates), amplitudes, flow ratios, durations and daily numbers of fluctuations, a predictive model (linear discriminant analysis) was fitted to classify hydrographs into predominant fluctuation regimes. This is the basis for a detailed investigation of present sub‐daily flow regimes and to analyse the differences between the regimes. Based on the results, we finally propose a general framework that enables a standardized assessment of flow fluctuations regarding event intensities and/or event timing. The proposed framework offers a standardized selection of particular flow fluctuations referring to increase and decrease events separately. The selection of specific flow fluctuations can be defined with respect to several research questions (e.g. ecologically relevant fluctuations), which offers a wide range of applications. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and temporal variability of groundwater dynamics in a sub‐Mediterranean mountain catchment

The temporal and spatial dynamics of groundwater was investigated in a small catchment in the Spanish Pyrenees, which was extensively used for agriculture in the past. Analysis of the water table fluctuations at five locations over a 6‐year period demonstrated that the groundwater dynamics had a marked seasonal cycle involving a wetting‐up period that commenced with the first autumn rainfall events, a saturation period during winter and spring and a drying‐down period from the end of spring until the end of the summer. The length of the saturation period showed great interannual variability, which was mainly influenced by the rainfall and evapotranspiration characteristics. There was marked spatial variability in the water table, especially during the wetting‐up period, which could be related to differences in slope and drainage area, geomorphology, soil properties and local topography. Areas contributing to runoff generation were identified within the catchment by field mapping of moisture conditions. Areas contributing to infiltration excess runoff were correlated with former cultivated fields affected by severe sheetwash erosion. Areas contributing to saturation excess runoff were characterized by a marked spatial dynamics associated with catchment wetness conditions. The saturation spatial pattern, which was partially related to the topographic index, was very patchy throughout the catchment, suggesting the influence of other factors associated with past agricultural activities, including changes in local topography and soil properties. The relationship between water table levels and stream flow was weak, especially during the wetting‐up period, suggesting little connection between ground water and the hydrological response, at least at some locations. The results suggest that in drier and human‐disturbed environments, such as sub‐Mediterranean mountains, saturation patterns cannot be represented only by the general topography of the catchment. They also suggest that groundwater storage and runoff is not a succession of steady‐state flow conditions, as assumed in most hydrological models. Copyright © 2013 John Wiley & Sons, Ltd.     

Water contents of clinopyroxenes from sub‐arc mantle peridotites

One poorly constrained reservoir of the Earth's water budget is that of clinopyroxene in metasomatised, mantle peridotites. This study presents reconnaissance Sensitive High‐Resolution, Ion Microprobe–Stable Isotope (SHRIMP–SI) determinations of the H₂O contents of (dominantly) clinopyroxenes in rare mantle xenoliths from four different subduction zones, i.e. Mexico, Kamchatka, Philippines, and New Britain (Tabar‐Feni island chain) as well as one intra‐plate setting (western Victoria). All of the sub‐arc xenoliths have been metasomatised and carry strong arc trace element signatures. Average measured H₂O contents of the pyroxenes range from 70 ppm to 510 ppm whereas calculated bulk H₂O contents range from 88 ppm to 3 737 ppm if the variable presence of amphibole is taken into account. In contrast, the intra‐plate, continental mantle xenolith from western Victoria has higher water contents (3 447 ppm) but was metasomatised by alkali and/or carbonatitic melts and does not carry a subduction‐related signature. Material similar to the sub‐arc peridotites can either be accreted to the base of the lithosphere or potentially be transported by convection deeper into the mantle where it will lose water due to amphibole breakdown.     

Urban fluvial flood modelling using a two‐dimensional diffusion‐wave treatment,part 2: development of a sub‐grid‐scale treatment

This paper develops and tests a sub‐grid‐scale wetting and drying correction for use with two‐dimensional diffusion‐wave models of urban flood inundation. The method recognizes explicitly that representations of sub‐grid‐scale topography using roughness parameters will provide an inadequate representation of the effects of structural elements on the floodplain (e.g. buildings, walls), as such elements not only act as momentum sinks, but also have mass blockage effects. The latter may dominate, especially in structurally complex urban areas. The approach developed uses high‐resolution topographic data to develop explicit parameterization of sub‐grid‐scale topographic variability to represent both the volume of a grid cell that can be occupied by the flow and the effect of that variability upon the timing and direction of the lateral fluxes. This approach is found to give significantly better prediction of fluvial flood inundation in urban areas than traditional calibration of sub‐grid‐scale effects using Manning's n. In particular, it simultaneously reduces the need to use exceptionally high values of n to represent the effects of using a coarser mesh process representation and increases the sensitivity of model predictions to variation in n. Copyright © 2005 John Wiley & Sons, Ltd.     

Pixel‐oriented land use classification in energy balance modelling

Mass and energy transfer between soil, vegetation and atmosphere is the process that allows to maintain an adequate energy and water balance in the earth–atmosphere system. However, the evaluation of the energy balance components, such as the net radiation and the sensible and latent heat fluxes, is characterized by significant uncertainties related to both the dynamic nature of heat transfer processes and surfaces heterogeneity. Therefore, a detailed land use classification and an accurate evaluation of vegetation spatial distribution are required for an accurate estimation of these variables. For this purpose, in the present article, a pixel‐oriented supervised classification was applied to obtain land use maps of the Basilicata region in Southern Italy by processing three Landsat TM and ETM+ satellite images. An accuracy analysis based on the overall accuracy index and the agreement Khat of Cohen coefficient showed a good performance of the applied classification methodology and a good quality of the obtained maps. Subsequently, these maps were used in the application of a simplified two‐source energy balance model for estimating the actual evapotranspiration at a regional scale. The comparison between the simulations made by applying the simplified two‐source energy balance model and the measurements of evapotranspiration at a lysimetric station located in the study area showed the applicability and the validity of the proposed methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Flow behaviour,suspended sediment transport and transmission losses in a small (sub‐bank‐full) flow event in an Australian desert stream

The behaviour of a discrete sub‐bank‐full flow event in a small desert stream in western NSW, Australia, is analysed from direct observation and sediment sampling during the flow event and from later channel surveys. The flow event, the result of an isolated afternoon thunderstorm, had a peak discharge of 9 m³/s at an upstream station. Transmission loss totally consumed the flow over the following 7·6 km. Suspended sediment concentration was highest at the flow front (not the discharge peak) and declined linearly with the log of time since passage of the flow front, regardless of discharge variation. The transmission loss responsible for the waning and eventual cessation of flow occurred at a mean rate of 13.2% per km. This is quite rapid, and is more than twice the corresponding figure for bank‐full flows estimated by Dunkerley (1992) on the same stream system. It is proposed that transmission losses in ephemeral streams of the kind studied may be minimized in flows near bank‐full stage, and be higher in both sub‐bank‐full and overbank flows. Factors contributing to enhanced flow loss in the sub‐bank‐full flow studied included abstractions of flow to pools, scour holes and other low points along the channel, and overflow abstractions into channel filaments that did not rejoin the main flow. On the other hand, losses were curtailed by the shallow depth of banks wetted and by extensive mud drapes that were set down over sand bars and other porous channel materials during the flow. Thus, in contrast with the relatively regular pattern of transmission loss inferred from large floods, losses from low flows exhibit marked spatial variability and depend to a considerable extent on streamwise variations in channel geometry, in addition to the depth and porosity of channel perimeter sediments. Copyright © 1999 John Wiley & Sons, Ltd.     

Evaluation of performance of non‐invasive upgrade strategy for beam–column sub‐assemblages of poorly designed structures under seismic type loading

Beam–column sub‐assemblages are the one of the most vulnerable structural elements to the seismic loading and may lead to devastating consequences. In order to improve the performance of the poorly/under‐designed building structures to the critical loading scenarios, introduction of steel bracing at the RC beam–column joint is found to be one of the modern and implementable techniques. In the present work, a diagonal metallic single haunch/bracing system is introduced at the beam–column joints to provide an alternate load path and to protect the joint zone from extensive damage because of brittle shear failure. In this paper, an investigation is reported on the evaluation of tae influence of different parameters, such as angle of inclination, location of bracing and axial stiffness of the single steel bracing on improving the performance through altering the force transfer mechanism. Numerical investigations on the performance of the beam–column sub‐assemblages have been carried out under cyclic loading using non‐linear finite element analysis. Experimentally validated numerical models (both GLD and upgraded specimen) have been further used for evaluating the performance of various upgrade schemes. Cyclic behaviour of reinforcement, concrete modelling based on fracture energy, bond‐slip relations between concrete and steel reinforcement have been incorporated. The study also includes the numerical investigation of crack and failure patterns, ultimate load carrying capacity, load displacement hysteresis, energy dissipation and ductility. The findings of the present study would be helpful to the engineers to develop suitable, feasible and efficient upgrade schemes for poorly designed structures under seismic loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Estimates of spatial variation in evaporation using satellite‐derived surface temperature and a water balance model

Evaporation dominates the water balance in arid and semi‐arid areas. The estimation of evaporation by land‐cover type is important for proper management of scarce water resources. Here, we present a method to assess spatial and temporal patterns of actual evaporation by relating water balance evaporation estimates to satellite‐derived radiometric surface temperature. The method is applied to a heterogeneous landscape in the Krishna River basin in south India using 10‐day composites of NOAA advanced very high‐resolution radiometer satellite imagery. The surface temperature predicts the difference between reference evaporation and modelled actual evaporation well in the four catchments (r² = 0·85 to r² = 0·88). Spatial and temporal variations in evaporation are linked to vegetation type and irrigation. During the monsoon season (June–September), evaporation occurs quite uniformly over the case‐study area (1·7–2·1 mm day^?1), since precipitation is in excess of soil moisture holding capacity, but it is higher in irrigated areas (2·2–2·7 mm day^?1). In the post‐monsoon season (December–March) evaporation is highest in irrigated areas (2·4 mm day^?1). A seemingly reasonable estimate of temporal and spatial patterns of evaporation can be made without the use of more complex and data‐intensive methods; the method also constrains satellite estimates of evaporation by the annual water balance, thereby assuring accuracy at the seasonal and annual time‐scales. Copyright © 2007 John Wiley & Sons, Ltd.     

Variable power‐law scaling of hillslope Hortonian rainfall–runoff processes

Hydrological studies focused on Hortonian rainfall–run‐off scaling have found that the run‐off depth generally declines with the plot length in power‐law scaling. Both the power‐law proportional coefficient and the scaling exponent show great variability for specific conditions, but why and how they vary remain unclear. In the present study, the scaling of hillslope Hortonian rainfall–run‐off processes is investigated for different rainfall, soil infiltration, and hillslope surface characteristics using the physically based cell‐based rainfall‐infiltration‐run‐off model. The results show that both temporally intermittent and steady rainfalls can result in prominent power‐law scaling at the initial stage of run‐off generation. Then, the magnitude of the power‐law scaling decreases gradually due to the decreasing run‐on effect. The power‐law scaling is most sensitive to the rainfall and soil infiltration parameters. When the ratio of rainfall to infiltration exceeds a critical value, the magnitude of the power‐law scaling tends to decrease notably. For different intermittent rainfall patterns, the power‐law exponent varies in the range of ?1.0 to ?0.113, which shows an approximately logarithmic increasing trend for the proportional coefficient as a function of the run‐off coefficient. The scaling is also sensitive to the surface roughness, soil sealing, slope angle, and hillslope geometry because these factors control the run‐off routing and run‐on infiltration processes. These results provide insights into the variable scaling of the Hortonian rainfall–run‐off process, which are expected to benefit modelling of large‐scale hydrological and ecological processes.     

Modification and testing of a one‐dimensional energy and mass balance model for supraglacial snowpacks

A one‐dimensional energy and mass balance snow model (SNTHERM) has been modified for use with supraglacial snowpacks and applied to a point on Haut Glacier d'Arolla, Switzerland. It has been adapted to incorporate the underlying glacier ice and a site‐specific, empirically derived albedo routine. Model performance was tested against continuous measurements of snow depth and meltwater outflow from the base of the snowpack, and intermittent measurements of surface albedo and snowpack density profiles collected during the 1993 and 2000 melt seasons. Snow and ice ablation was simulated accurately. The timing of the daily pattern of meltwater outflow was well reproduced, although magnitudes were generally underestimated, possibly indicating preferential flow into the snowpack lysimeter. The model was used to assess the quantity of meltwater stored temporally within the unsaturated snowpack and meltwater percolation rates, which were found to be in agreement with dye tracer experiments undertaken on this glacier. As with other energy balance studies on alpine valley glaciers, the energy available for melt was dominated by net radiation (64%), with a sizable contribution from sensible heat flux (36%) and with a negligible latent heat flux overall, although there was more complex temporal variation on diurnal timescales. A basic sensitivity analysis indicated that melt rates were most sensitive to radiation, air temperature and snowpack density, indicating the need to accurately extrapolate/interpolate these variables when developing a spatially distributed framework for this model. Copyright © 2008 John Wiley & Sons, Ltd.     

Energy balance above a boreal coniferous forest: a difference in turbulent fluxes between snow‐covered and snow‐free canopies

To evaluate the interactive effects of snow and forest on turbulent fluxes between the forest surface and the atmosphere, the surface energy balance above a forest was measured by the eddy correlation method during the winter of 1995–1996. The forest was a young coniferous plantation comprised of spruce and fir. The study site, in Sapporo, northern Japan, had heavy and frequent snowfalls and the canopy was frequently covered with snow during the study period. A comparison of the observed energy balance above the forest for periods with and without a snow‐covered canopy and an analysis using a single‐source model gave the following results: during daytime when the canopy was covered with snow, the upward latent heat flux was large, about 80% of the net radiation, and the sensible heat flux was positive but small. On the other hand, during daytime when the canopy was dry and free from snow, the sensible heat flux was dominant and the latent heat flux was minor, about 10% of the net radiation. To explain this difference of energy partition between snow‐covered and snow‐free conditions, not only differences in temperature but also differences in the bulk transfer coefficients for latent heat flux were necessary in the model. Therefore, the high evaporation rate from the snow‐covered canopy can be attributed largely to the high moisture availability of the canopy surface. Evaporation from the forest during a 60‐day period in midwinter was estimated on a daily basis as net radiation minus sensible heat flux. The overall average evaporation during the 60‐day period was 0·6 mm day⁻¹, which is larger than that from open snow fields. Copyright © 1999 John Wiley & Sons, Ltd.     

Incipient weathering rind development on introduced machine‐polished granite discs in an Arctic alpine environment,northern Scandinavia

Weathering rinds, zones of alteration on the exterior surfaces of rock outcrops and coarse unconsolidated surficial debris are widely used by geomorphologists and Quaternary geologists as indicators of the relative age of landforms and landscapes. Additionally they provide unique insights into the earliest stages of rock and mineral weathering, yet the origin of these alteration zones is relatively poorly understood. This lack of understanding applies especially to the initial stages of rind formation. The study reported in this paper has two principal objectives. The first is to use lightly polished granite discs inserted in soil profiles under several different plant communities in an Arctic alpine environment for a period of four or five years to investigate the nature of incipient weathering rind development. The second is to investigate the factors responsible for spatial variability in the nature and rates of rind formation. Incipient weathering rind development on the outer edges of the granite discs is observable and measurable over a period of time as short as four years in the mild Arctic alpine environment of Swedish Lapland. The earliest stages of rind development involve the development of a porous structure consisting of a combination of pits and fractures which have been solutionally enlarged and modified. Solution appears to be preferentially concentrated on the surfaces of feldspars and, to a lesser extent, quartz. In addition, iron oxides are present along grain boundaries and in grain interiors and are interpreted to have been derived from the oxidation of ferromagnesian minerals. Spatial variability in weathering rind development appears to be particularly driven by differences in moisture but is not related to soil pH. Copyright © 2005 John Wiley & Sons, Ltd.