Impact of wildfire on source water contributions in Devil Creek,CA: evidence from end‐member mixing analysis

A geochemical and end‐member mixing analysis (EMMA) is undertaken in Devil Canyon catchment, located in southern California, to further understanding of watershed behaviour and source water contributions after an acute and extensive wildfire. Physical and chemical transformations in post‐fire watersheds are known to increase overland flow and decrease infiltration, mainly due to formation of a hydrophobic layer at, or near, the soil surface. However, less is known about subsurface flow response in burned watersheds. The current study incorporates EMMA to evaluate and quantify source water contributions before, and after, a catchment affected by wildfires in southern California during the fall of 2003. Pre‐ and post‐fire stream water data were available at several sampling sites within the catchment, allowing the identification of contributing water sources at varying spatial scales. Proposed end‐member observations (groundwater, overland flow, shallow subsurface flow) were also collected to constrain and develop the catchment mixing model. Post‐fire source water changes are more evident in the smaller and faster responding sub‐basin (interior sampling point). Early post‐fire storm events are dominated by overland flow with no significant soil water or groundwater flow contribution. Inter‐storm streamwater in this smaller basin shows an increase in groundwater and a decrease in soil water. In the larger, baseflow‐dominated system, source water components appear less affected by fire. A slight increase in lateral flow is observed with only a slight decrease in baseflow. Changes in the post‐fire flow regimes affect nutrient loading and chemical response of the basin. Relatively rapid recovery of the chaparral ecosystem is evidenced, with active re‐growth and evapotranspiration evidenced by the fourth post‐fire rainy season. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of integrated back‐propagation network and self organizing map for flood forecasting

Much of the nonlinearity and uncertainty regarding the flood process is because hydrologic data required for estimation are often tremendously difficult to obtain. This study employed a back‐propagation network (BPN) as the main structure in flood forecasting to learn and to demonstrate the sophisticated nonlinear mapping relationship. However, a deterministic BPN model implies high uncertainty and poor consistency for verification work even when the learning performance is satisfactory for flood forecasting. Therefore, a novel procedure was proposed in this investigation which integrates linear transfer function (LTF) and self‐organizing map (SOM) to efficiently determine the intervals of weights and biases of a flood forecasting neural network to avoid the above problems. A SOM network with classification ability was applied to the solutions and parameters of the BPN model in the learning stage, to classify the network parameter rules and to obtain the winning parameters. The outcomes from the previous stage were then used as the ranges of the parameters in the recall stage. Finally, a case study was carried out in Wu‐Shi basin to demonstrate the effectiveness of the proposal. Copyright © 2009 John Wiley & Sons, Ltd.     

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S_f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub‐reaches. Training set required for development of the fuzzy rule based model for each sub‐reach is obtained from discharge‐area relationship at its mean section. For highly heterogeneous sub‐reaches, optimized fuzzy rule based models are obtained by means of a neuro‐fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub‐reaches and also in a natural channel with a number of approximately uniform sub‐reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM outputs are comparable to those of the HEC‐RAS model. Copyright © 2009 John Wiley & Sons, Ltd.     

Effects of land‐use changes on hydrological processes in the middle basin of the Heihe River,northwest China

The effects of land‐use changes on the runoff process in the midstream plain of this arid inland river basin are a key factor in the rational allocation of water resources to the middle and lower reaches. The question is whether and by how much increasingly heavy land use impacts the hydrological processes in such an arid inland river basin. The catchment of the Heihe River, one of the largest inland rivers in the arid region of northwest China, was chosen to investigate the hydrological responses to land‐use change. Flow duration curves were used to detect trends and variations in runoff between the upper and lower reaches. Relationships among precipitation, upstream runoff, and hydrological variables were identified to distinguish the effects of climatic changes and upstream runoff changes on middle and downstream runoff processes. The quantitative relation between midstream cultivated land use and various parameters of downstream runoff processes were analysed using the four periods of land‐use data since 1956. The Volterra numerical function relation of the hydrological non‐linear system response was utilized to develop a multifactor hydrological response simulation model based on the three factors of precipitation, upstream runoff, and cultivated land area. The results showed that, since 1967, the medium‐ and high‐coverage natural grassland area in the midstream region has decreased by 80·1%, and the downstream runoff has declined by 27·32% due to the continuous expansion of the cultivated land area. The contribution of cultivated land expansion to the impact on the annual total runoff is 14–31%, on the annual, spring and winter base flow it is 44–75%, and on spring and winter discharge it is 23–64%. Once the water conservation plan dominated by land‐use structural adjustments is implemented over the next 5 years, the mean annual discharge in the lower reach could increase by 8·98% and the spring discharge by 26·28%. This will significantly alleviate the imbalance between water supply and demand in both its quantity and temporal distribution in the middle and lower reaches. Copyright © 2006 John Wiley & Sons, Ltd.     

The effects of rock fragments on hydrologic and hydraulic responses along a slope

A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd.     

Estimation of potential evapotranspiration over the Yellow River basin: reference crop evaporation or Shuttleworth–Wallace?

Potential evapotranspiration (PET) is a key input to hydrological models. Its estimation has often been via the Penman–Monteith (P–M) equation, most recently in the form of an estimate of reference evapotranspiration (RET) as recommended by FAO‐56. In this paper the Shuttleworth–Wallace (S–W) model is implemented to estimate PET directly in a form that recognizes vegetation diversity and temporal change without reference to experimental measurements and without calibration. The threshold values of vegetation parameters are drawn from the literature based on the International Geosphere–Biosphere Programme land cover classification. The spatial and temporal variation of the LAI of vegetation is derived from the composite NOAA‐AVHRR normalized difference vegetation index (NDVI) using a method based on the SiB2 model, and the Climate Research Unit database is used to provide the required meteorological data. All these data inputs are publicly and globally available. Consequently, the implementation of the S–W model developed in this study is applicable at the global scale, an essential requirement if it is to be applied in data‐poor or ungauged large basins. A comparison is made between the FAO‐56 method and the S–W model when applied to the Yellow River basin for the whole of the last century. The resulting estimates of RET and PET and their association with vegetation types and leaf area index (LAI) are examined over the whole basin both annual and monthly and at six specific points. The effect of NDVI on the PET estimate is further evaluated by replacing the monthly NDVI product with the 10‐day product. Multiple regression relationships between monthly PET, RET, LAI, and climatic variables are explored for categories of vegetation types. The estimated RET is a good climatic index that adequately reflects the temporal change and spatial distribution of climate over the basin, but the PET estimated using the S–W model not only reflects the changes in climate, but also the vegetation distribution and the development of vegetation in response to climate. Although good statistical relationships can be established between PET, RET and/or climatic variables, applying these relationships likely will result in large errors because of the strong non‐linearity and scatter between the PET and the LAI of vegetation. It is concluded that use of the implementation of the S–W model described in this study results in a physically sound estimate of PET that accounts for changing land surface conditions. Copyright © 2006 John Wiley & Sons, Ltd.     

Statistical characterization of impact of system variables on temporal dynamics of groundwater in highly weathered regoliths

In situations where the water table fluctuates during the rainy season the characterization of the impact of system variables on the temporal dynamics of the groundwater (GW) is essential to improve the understanding at catchment or regional scale behaviour of GW. In this study the appropriateness of the statistical parameters; mean, median, the 80^th percentile (PC80), coefficient of variation (CV), correlation coefficient (r), and multiple regression models were assessed to characterize the impact of system variables on the temporal dynamics of hydraulic head relative to ground surface (HH) during rainy seasons. The study was conducted from 1999 to 2003 in the wet tropical Johnstone River catchment (JRC) in north‐east Queensland, Australia. Piezometer wells were installed at 32 sites under cropping to 5–90 m depth on different soil types, landscape positions, and varying proximity to surface water bodies (i.e. four system variables). The HH was measured, at least at 10–15 day intervals during 1–5 consecutive rainy seasons. The HH in the 32 wells fluctuated throughout each of the five rainy seasons. The mean HH averaged over the seasons ranged from 1·1 to 17·2 m across the wells, the median from 0·9 to 17·3 m, and the PC80 from 0·3 to 16·1 m. The temporal behaviour of HH characterization by mean of means of HH, the mean of medians of HH, and the mean of PC80 of HH, indicated the HH can be classified to belong to three different groups for each one of these parameters. The impact of the system variables on temporal dynamics, explored using multiple regression procedure, indicated that the model for median was marginally better than mean. The CV was found to be most appropriate parameter to characterize the impact of GW system variable (aquifer type), a component of the system variables, on temporal dynamics. The interactions of GW (i) belonging to different GW system and (ii) at shoulder with footslope in a landscape were best characterized by simple linear correlations. Copyright © 2007 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.     

Spatial variability of snowmelt timing from AMSR‐E and SSM/I passive microwave sensors,Pelly River,Yukon Territory,Canada

Spring snow melt run‐off in high latitude and snow‐dominated drainage basins is generally the most significant annual hydrological event. Melt timing, duration, and flow magnitude are highly variable and influence regional climate, geomorphology, and hydrology. Arctic and sub‐arctic regions have sparse long‐term ground observations and these snow‐dominated hydrologic regimes are sensitive to the rapidly warming climate trends that characterize much of the northern latitudes. Passive microwave brightness temperatures are sensitive to changes in the liquid water content of the snow pack and make it possible to detect incipient melt, diurnal melt‐refreeze cycles, and the approximate end of snow cover on the ground over large regions. Special Sensor Microwave Imager (SSM/I) and Advanced Microwave Scanning Radiometer for EOS (AMSR‐E) passive microwave brightness temperatures (T_b) and diurnal amplitude variations (DAV) are used to investigate the spatial variability of snowmelt onset timing (in two stages, ‘DAV onset’ and ‘melt onset’) and duration for a complex sub‐arctic landscape during 2005. The satellites are sensitive to small percentages of liquid water, and therefore represent ‘incipient melt’, a condition somewhat earlier than a traditional definition of a melting snowpack. Incipient melt dates and duration are compared to topography, land cover, and hydrology to investigate the strength and significance of melt timing in heterogeneous landscapes in the Pelly River, a major tributary to the Yukon River. Microwave‐derived melt onset in this region in 2005 occurred from late February to late April. Upland areas melt 1–2 weeks later than lowland areas and have shorter transition periods. Melt timing and duration appear to be influenced by pixel elevation, aspect, and uniformity as well as other factors such as weather and snow mass distribution. The end of the transition season is uniform across sensors and across the basin in spite of a wide variety of pixel characteristics. Copyright © 2007 John Wiley & Sons, Ltd.     

Coupled HM analysis using zero‐thickness interface elements with double nodes—Part II: Verification and application

In a companion Part I of this paper (Int. J. Numer. Anal. Meth. Geomech. 2008; DOI: 10.1002/nag.735 ), a coupled hydro‐mechanical (HM) formulation for geomaterials with discontinuities based on the finite element method (FEM) with double‐node, zero‐thickness interface elements was developed and presented. This Part II paper includes the numerical solution of basic practical problems using both the staggered and the fully coupled approaches. A first group of simulations, based on the classical consolidation problem with an added vertical discontinuity, is used to compare both the approaches in terms of accuracy and convergence. The monolithic or fully coupled scheme is also used in an application example studying the influence of a horizontal joint in the performance of a reservoir subject to fluid extraction. Results include a comparison with other numerical solutions from the literature and a sensitivity analysis of the mechanical parameters of the discontinuity. Some simulations are also run using both a full non‐symmetric and a simplified symmetric Jacobian matrix. On top of verifying the model developed and its capability to reflect the conductivity changes of the interface with aperture changes, the results presented also lead to interesting observations of the numerical performance of the methods implemented. Copyright © 2008 John Wiley & Sons, Ltd.