Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

A simple parameterization of bulk canopy resistance from climatic variables for estimating hourly evapotranspiration

This paper examines a model for estimating canopy resistance r_c and reference evapotranspiration ET_o on an hourly basis. The experimental data refer to grass at two sites in Spain with semiarid and windy conditions in a typical Mediterranean climate. Measured hourly ET_o values were obtained over grass during a 4 year period between 1997 and 2000 using a weighing lysimeter (Zaragoza, northeastern Spain) and an eddy covariance system (Córdoba, southern Spain). The present model is based on the Penman–Monteith (PM) approach, but incorporates a variable canopy resistance r_c as an empirical function of the square root of a climatic resistance r* that depends on climatic variables. Values for the variable r_c were also computed according to two other approaches: with the r_c variable as a straight‐line function of r* (Katerji and Perrier, 1983, Agronomie 3 (6): 513–521) and as a mechanistic function of weather variables as proposed by Todorovic (1999, Journal of Irrigation and Drainage Engineering, ASCE 125 (5): 235–245). In the proposed model, the results show that r_c/r_a (where r_a is the aerodynamic resistance) presents a dependence on the square root of r*/r_a, as the best approach with empirically derived global parameters. When estimating hourly ET_o values, we compared the performance of the PM equation using those estimated variable r_c values with the PM equation as proposed by the Food and Agriculture Organization, with a constant r_c = 70 s m^?1. The results confirmed the relative robustness of the PM method with constant r_c, but also revealed a tendency to underestimate the measured values when ET_o is high. Under the semiarid conditions of the two experimental sites, slightly better estimates of ET_o were obtained when an estimated variable r_c was used. Although the improvement was limited, the best estimates were provided by the Todorovic and the proposed methods. The proposed approach for r_c as a function of the square root of r* may be considered as an alternative for modelling r_c, since the results suggest that the global coefficients of this locally calibrated relationship might be generalized to other climatic regions. It may also be useful to incorporate the effects of variable canopy resistances into other climatic and hydrological models. Copyright © 2005 John Wiley & Sons, Ltd.     

Runoff generation and routing on artificial slopes in a Mediterranean–continental environment: the Teruel coalfield,Spain

The aim of this study was to identify the mechanisms of runoff generation and routing and their controlling factors at the hillslope scale, on artificial slopes derived from surface coal mining reclamation in a Mediterranean–continental area. Rainfall and runoff at interrill and microcatchment scales were recorded for a year on two slopes with different substrata: topsoil cover and overburden cover. Runoff coefficient and runoff routing from interrill areas to microcatchment outlets were higher in the overburden substratum than in topsoil, and greater in the most developed rill network. Rainfall volume is the major parameter responsible for runoff response on overburden, suggesting that this substratum is very impermeable—at least during the main rainfall periods of the year (late spring and autumn) when the soil surface is sealed. In such conditions, most rainfall input is converted into runoff, regardless of its intensity. Results from artificial rainfall experiments, conducted 3 and 7 years after seeding, confirm the low infiltration capacity of overburden when sealed. The hydrological response shows great seasonal variability on the overburden slope in accordance with soil surface changes over the year. Rainfall volume and intensities (I₃₀, I₆₀) explain runoff at the interrill scale on the topsoil slope, where rainfall experiments demonstrated a typical Hortonian infiltration curve. However, no correlation was found at the microcatchment level, probably because of the loss of functionality of the only rill as ecological succession proceeded. The runoff generation mechanism on the topsoil slope is more homogeneous throughout the year. Runoff connectivity, defined as the ratio between runoff rates recorded at the rill network scale and those recorded at the interrill area scale in every rainfall event, was also greater on the rilled overburden slope, and in the most developed rill network. The dense rill networks of the overburden slope guarantee very effective runoff drainage, regardless of rainfall magnitude. Rills drain overland flow from interrill‐sealed areas, reducing the opportunity of reinfiltration in areas not affected by siltation. Runoff generation and routing on topsoil slopes are controlled by grass cover and soil moisture content, whereas on overburden slopes rill network density and soil moisture content are the main controlling factors. Copyright © 2002 John Wiley & Sons, Ltd.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

Performance of Jacobi preconditioning in Krylov subspace solution of finite element equations

This paper examines the performance of the Jacobi preconditioner when used with two Krylov subspace iterative methods. The number of iterations needed for convergence was shown to be different for drained, undrained and consolidation problems, even for similar condition number. The differences were due to differences in the eigenvalue distribution, which cannot be completely described by the condition number alone. For drained problems involving large stiffness ratios between different material zones, ill‐conditioning is caused by these large stiffness ratios. Since Jacobi preconditioning operates on degrees‐of‐freedom, it effectively homogenizes the different spatial sub‐domains. The undrained problem, modelled as a nearly incompressible problem, is much more resistant to Jacobi preconditioning, because its ill‐conditioning arises from the large stiffness ratios between volumetric and distortional deformational modes, many of which involve the similar spatial domains or sub‐domains. The consolidation problem has two sets of degrees‐of‐freedom, namely displacement and pore pressure. Some of the eigenvalues are displacement dominated whereas others are excess pore pressure dominated. Jacobi preconditioning compresses the displacement‐dominated eigenvalues in a similar manner as the drained problem, but pore‐pressure‐dominated eigenvalues are often over‐scaled. Convergence can be accelerated if this over‐scaling is recognized and corrected for. Copyright © 2002 John Wiley & Sons, Ltd.     

Water temperature and conductivity variability as indicators of groundwater behaviour in complex aquifer systems in the south‐east of Spain

Groundwater temperature and electrical conductivity logs were applied in this study to gain an understanding of the hydrogeological behaviour of a complex system. This system, compartmentalized by a large number of fractures, which may even cause geothermal anomalies, is that of the Campo de Dalías in the south‐east of the Iberian Peninsula. In addition to this complexity one must consider the different levels with aquifer properties and their varying lithologies, which in some cases are in direct contact with the sea. Vertical sections based upon the various logs allow lateral correlation to be made between the anomalies detected, and provide information on the existence of both ascending and descending vertical flows, and the inflow of sea water into the system. Analysis of environmental isotopes (¹⁸O and D) indicate the areas of recharge with the greatest flow velocity. No seasonal variation in the isotopes was observed, which suggests that the residence time of the groundwater in the aquifer is several years. Radiocarbon dating gave estimates of generally long residence times, although in some sectors it is possible to detect the current recharge based on its tritium content, from rainwater infiltration or seawater intrusion. Copyright © 2002 John Wiley & Sons, Ltd.     

Biodegradation and transport of benzene,toluene, and xylenes in a simulated aquifer: comparison of modelled and experimental results

Both laboratory experiments and numerical modelling were conducted to study the biodegradation and transport of benzene–toluene–xylenes (BTX) in a simulated semi‐confined aquifer. The factors incorporated into the numerical model include advection, hydrodynamic dispersion, adsorption, and biodegradation. The various physico‐chemical parameters required by the numerical model were measured experimentally. In the experimental portion of the study, BTX compounds were introduced into the aquifer sand. After the contaminants had been transported through the system, BTX concentrations were measured at 12 equally spaced wells. Subsequently, microorganisms obtained from the activated sludge of a sewage treatment plant and cultured in BTX mixtures were introduced into the aquifer through the 12 sampling wells. The distribution data for BTX adsorption by the aquifer sand form a nonlinear isotherm. The degree of adsorption by the sand varies, depending on the composition of the solute. The degradation time, measured from the time since the bacteria were added to the aquifer until a specific contaminant was no longer detectable, was 35–42 h for BTX. The dissolved oxygen, after degradation by BTX compounds and bacteria, was consumed by about 40–60% in the entire simulated aquifer; thus the aerobic conditions were maintained. This study provides insights for the biodegradation and transport of BTX in aquifers by numerical modelling and laboratory experiments. Experimental and numerical comparisons indicate that the results by Monod degradation kinetics are more accurate than those by the first‐order degradation kinetics. Copyright © 2002 John Wiley & Sons, Ltd.     

Image registration for the underwater inspection using the maximum a posteriori technique

Describes an image registration method for underwater inspection tasks. A remotely operated vehicle equipped with a video camera and a scanning sonar is used as the testbed vehicle. Each image of the underwater scene is saved along with the video camera's position and orientation. The images are then combined to create a large composite picture of the underwater structure being inspected. This method is based upon a maximum a posteriori estimation technique and provides smooth and robust estimates of image shifts. Our results demonstrate the feasibility of this highly promising underwater inspection procedure.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.