Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Differences in plant function in phosphorus- and nitrogen-limited mangrove ecosystems

Mangrove ecosystems can be either nitrogen (N) or phosphorus (P) limited and are therefore vulnerable to nutrient pollution. Nutrient enrichment with either N or P may have differing effects on ecosystems because of underlying differences in plant physiological responses to these nutrients in either N- or P-limited settings. Using a common mangrove species, Avicennia germinans, in sites where growth was either N or P limited, we investigated differing physiological responses to N and P limitation and fertilization. We tested the hypothesis that water uptake and transport, and hydraulic architecture, were the main processes limiting productivity at the P-limited site, but that this was not the case at the N-limited site. We found that plants at the P-deficient site had lower leaf water potential, stomatal conductance and photosynthetic carbon-assimilation rates, and less conductive xylem, than those at the N-limited site. These differences were greatly reduced with P fertilization at the P-limited site. By contrast, fertilization with N at the N-limited site had little effect on either photosynthetic or hydraulic traits. We conclude that growth in N- and P-limited sites differentially affect the hydraulic pathways of mangroves. Plants experiencing P limitation appear to be water deficient and undergo more pronounced changes in structure and function with relief of nutrient deficiency than those in N-limited ecosystems.

Efffluent flux prediction in variably saturated soil zones within a septic tank-soil absorption trench

The treatment and hydraulic mechanisms in a septic tank-soil absorption system ( SAS) are highly influenced by the clogging layer or biomat zone which develops on bottom and lower sidewall surfaces within the trench. Flow rates through the biomat and sub-biomat zones are governed largely by the biomat hydraulic properties (resistance and hydraulic conductivity) and the unsaturated hydraulic conductivity of the underlying soil. One- and 2-dimensional models were used to investigate the relative importance of sidewall and vertical flow rates and pathways in SAS. Results of 1-dimensional modelling show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) reduce to a 1 order of magnitude variation in long-term flow rates. To increase the reliability of prediction of septic trench hydrology, HYDRUS-2D was used to model 2-dimensional flow. In the permeable soils, under high trench loading, effluent preferentially flowed in the upper region of the trench where no resistant biomat was present (the exfiltration zone). By comparison, flow was more evenly partitioned between the biomat zones and the exfiltration zones of the low permeability soil. An increase in effluent infiltration corresponded with a greater availability of exfiltration zone, rather than a lower resistance of biomat. Results of modelling simulations demonstrated the important role that a permeable A horizon may play in limiting surface surcharge of effluent under high trench hydraulic loading.

Long-term flow rates and biomat zone hydrology in soil columns receiving septic tank effluent

Soil absorption systems (SAS) are used commonly to treat and disperse septic tank effluent (STE). SAS can hydraulically fail as a result of the low permeable biomat zone that develops on the infiltrative surface. The objectives of this experiment were to compare the hydraulic properties of biomats grown in soils of different textures, to investigate the long-term acceptance rates (LTAR) from prolonged application of STE, and to assess if soils were of major importance in determining LTAR. The STE was applied to repacked sand, Oxisol and Vertisol soil columns over a period of 16 months, at equivalent hydraulic loading rates of 50, 35 and 8 L/m(2)/d, respectively Infiltration rates, soil matric potentials, and biomat hydraulic properties were measured either directly from the soil columns or calculated using established soil physics theory. Biomats 1 to 2 cm thick developed in all soils columns with hydraulic resistances of 27 to 39 d. These biomats reduced a 4 order of magnitude variation in saturated hydraulic conductivity (K.) between the soils to a one order of magnitude variation in LTAR. A relationship between biomat resistance and organic loading rate was observed in all soils. Saturated hydraulic conductivity influenced the rate and extent of biomat development. However, once the biomat was established, the LTAR was governed by the resistance of the biomat and the sub-biomat soil unsaturated flow regime induced by the biomat. Results show that whilst initial soil K. is likely to be important in the establishment of the biomat zone in a trench, LTAR is determined by the biomat resistance and the unsaturated soil hydraulic conductivity, not the K, of a soil. The results call into question the commonly used approach of basing the LTAR, and ultimately trench length in SAS, on the initial K, of soils. (c) 2006 Elsevier Ltd. All rights reserved.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Modelling trickle irrigation: Comparison of analytical and numerical models for estimation of wetting front position with time

Irrigation practices that are profligate in their use of water have come under closer scrutiny by water managers and the public. Trickle irrigation has the propensity to increase water use efficiency but only if the system is designed to meet the soil and plant conditions. Recently we have provided a software tool, WetUp (http://www.clw.csiro.au/products/wetup/), to calculate the wetting patterns from trickle irrigation emitters. WetUp uses an analytical solution to calculate the wetted perimeter for both buried and surface emitters. This analytical solution has a number of assumptions, two of which are that the wetting front is defined by water content at which the hydraulic conductivity (K) is I mm day(-1) and that the flow occurs from a point source. Here we compare the wetting patterns calculated with a 2-dimensional numerical model, HYDRUS2D, for solving the water flow into typical soils with the analytical solution. The results show that the wetting patterns are similar, except when the soil properties result in the assumption of a point source no longer being a good description of the flow regime. Difficulties were also experienced with getting stable solutions with HYDRUS2D for soils with low hydraulic conductivities. (c) 2005 Elsevier Ltd. All rights reserved.

Changes in the properties of a Vertisol and responses of wheat after compaction with harvester traffic

Soil compaction has been recognised as the greatest problem in terms of damage to Australia's soil resource. Compaction by tractor and harvester tyres, related to trafficking of wet soil, is one source of the problem. In this paper an array of soil properties was measured before and immediately after the application of a known compaction force to a wet Vertisol, A local grain harvester was used on soil that was just trafficable; a common scenario at harvest. The primary aim was to determine the changes in various soil properties in order to provide a benchmark against which the effectiveness of future remedial treatments could be evaluated. A secondary aim was a comparison of the measurements' efficiency to assess a soil's structural degradation status. Also assessed was the subsequent effect of the applied compaction on wheat growth and yield in the following cropping season. Nine of the soil properties measured gave statistically significant differences as a result of the soil compaction. Differences were mostly restricted to the top 0.2 m of the soil. The greatest measured depth of effect was decreased soil porosity to 0.4 m measured from intact soil clods. There was 72% emergence of the wheat crop planted into the compact soil and 93% in the uncompact soil. Wheat yield, however, was not affected by the compaction. This may demonstrate that wheat, growing on a full profile of stored soil water as did the current crop, may be little affected by compaction, Also, wheat may have potential to facilitate rapid repair of the damage in a Vertisol such as the current soil by drying the topsoil between rainfall events so increasing shrinking and swelling cycles. If this is true, then sowing a suitable crop species in a Vertisol may be a better option than tillage for repairing compaction damage by agricultural traffic. (C) 2000 Elsevier Science B.V. All rights reserved.

Prediction of steady-state flux through variably saturated zones within a septic absorption trench

This paper describes effluent flow dynamics within a septic absorption system and the prediction of flow through the biomat and sub-biomat zone. Using soil hydraulic properties in a one dimensional model we demonstrate how soil hydraulic properties interact with biomat resistances to determine long-term acceptance rate (LTAR). The LTAR is a key parameter used in the Australian and New Zealand Standard AS1547:2000 to calculate the area of trench required to ensure trenches are not overloaded. Results show that several orders of magnitude variation in saturated hydraulic conductivity (Ks) collapse to a one order of magnitude variation in LTAR. These results are calculated from a model using basic flow theory, allowing LTAR to be estimated for any combination of biomat resistance and soil hydraulic properties. To increase the reliability of prediction of septic trench hydrology, HYDRUS 2D was used to model two dimensional flow. For more permeable soils, the exfiltration zone above sidewall biomat growth is shown to be a key pathway for excess effluent flow.

Turbulence in Hydraulic Jumps: Experimental Measurements

A hydraulic jump is the transition from a supercritical open channel flow to a subcritical regime. It is characterised by a highly turbulent flow with macro-scale vortices, some kinetic energy dissipation and a bubbly two-phase flow structure. New air-water flow measurements were performed in hydraulic jump flows for a range of inflow Froude numbers. The experiments were conducted in a large-size facility using two types of phase-detection intrusive probes: i.e., single-tip and double-tip conductivity probes. These were complemented by some measurements of free-surface fluctuations using ultrasonic displacement meters. The present study was focused on the turbulence characteristics of hydraulic jumps with partially-developed inflow conditions. The void fraction measurements showed the presence of an advective diffusion shear layer in which the void fractions profiles matched closely an analytical solution of the advective diffusion equation for air bubbles. The present results highlighted some influence of the inflow Froude number onto the air bubble entrainment process. At the largest Froude numbers, the advected air bubbles were more thoroughly dispersed vertically, and larger amount of air bubbles were detected in the turbulent shear layer. In the air-water mixing layer, the maximum void fraction and bubble count rate data showed some longitudinal decay function in the flow direction. Such trends were previously reported in the literature. The measurements of interfacial velocity and turbulence level distributions provided new information on the turbulent velocity field in the highly-aerated shear region. The present data suggested some longitudinal decay of the turbulence intensity. The velocity profiles tended to follow a wall jet flow pattern. The air–water turbulent time and length scales were deduced from some auto- and cross-correlation analyses based upon the method of CHANSON (2006,2007). The results provided the integral turbulent time and length scales of the eddy structures advecting the air bubbles in the developing shear layer. The experimental data showed that the auto-correlation time scale Txx was larger than the transverse cross-correlation time scale Txz. The integral turbulence length scale Lxz was a function of the inflow conditions, of the streamwise position (x-x1)/d1 and vertical elevation y/d1. Herein the dimensionless integral turbulent length scale Lxz/d1 was closely related to the inflow depth: i.e., Lxz/d1 = 0.2 to 0.8, with Lxz increasing towards the free-surface. The free-surface fluctuations measurements showed large turbulent fluctuations that reflected the dynamic, unsteady structure of the hydraulic jumps. A linear relationship was found between the normalized maximum free-surface fluctuation and the inflow Froude number.

Nanocomposite Polymer Electrolyte Membranes: Methanol Crossover and Conductivity

Commercial Nafion® 117 membranes were successfully modified by in-situ reactions (sol-gel of TEOS and/or polymerization of aniline) within Nafion structures. Water-methanol permeability and proton conductivity were investigated in order to determine the potential performance of these membranes for DMFC systems. Silica-polyaniline modification resulted in 84% methanol crossover reduction, from 2.45x10^-5 cm2.s^-1 for conventional Nafion membranes to 3.71x10^-6 cm2.s^-1 for the modified silica-polyaniline composite membrane at 75 degrees C. In addition, conductivity was not hindered, as the polyaniline-Nafion membrane increased from 12.2 to 15 mS.cm^-1 as compared to Nafion, while a reduction of 11% was observed for silica-polyaniline-Nafion composite membrane. The results in this work strongly suggest the potential of polyaniline nanocomposites to enhance the performance of DMFCs.

Research Quality, Publications and Impact in Hydraulic Engineering into the 21st Century. Publish or Perish, Commercial versus Open Access, Internet versus Libraries

One of the main objectives of the first International Junior Researcher and Engineer Workshop on Hydraulic Structures is to provide an opportunity for young researchers and engineers to present their research. But a research project is only completed when it has been published and shared with the community. Referees and peer experts play an important role to control the research quality. While some new electronic tools provide further means to disseminate some research information, the quality and impact of the works remain linked with some thorough expert-review process and the publications in international scientific journals and books. Importantly unethical publishing standards are not acceptable and cheating is despicable.

Bubbly Flow Structure in Hydraulic Jump

In an open channel, a hydraulic jump is the rapid transition from super- to sub-critical flow associated with strong turbulence and air bubble entrainment in the mixing layer. New experiments were performed at relatively large Reynolds numbers using phase-detection probes. Some new signal analysis provided characteristic air-water time and length scales of the vortical structures advecting the air bubbles in the developing shear flow. An analysis of the longitudinal air-water flow structure suggested little bubble clustering in the mixing layer, although an interparticle arrival time analysis showed some preferential bubble clustering for small bubbles with chord times below 3 ms. Correlation analyses yielded longitudinal air-water time scales Txx*V1/d1 of about 0.8 in average. The transverse integral length scale Z/d1 of the eddies advecting entrained bubbles was typically between 0.25 and 0.4, irrespective of the inflow conditions within the range of the investigations. Overall the findings highlighted the complicated nature of the air-water flow

Hydraulic Design of Stepped Spillways and Downstream Energy Dissipators for Embankment Dams

In recent years, the design flows of many dams were re-evaluated, often resulting in discharges larger than the original design. In many cases, the occurrence of the revised flows could result in dam overtopping because of insufficient storage and spillway capacity. An experimental study was conducted herein to gain a better understanding of the flow properties in stepped chutes with slopes typical of embankment dams. The work was based upon a Froude similitude in large-size experimental facilities. A total of 10 configurations were tested including smooth steps, steps equipped with devices to enhance energy dissipation and rough steps. The present results yield a new design procedure. The design method includes some key issues not foreseen in prior studies : e.g., gradually varied flow, type of flow regime, flow resistance. It is believed that the outcomes are valid for a wide range of chute geometry and flow conditions typical of embankment chutes.

Diseño hidráulico de vertedores escalonados con pendientes moderadas: Metodología basada en un estudio experimental

Durante las últimas tres décadas el interés y diversidad en el uso de canales escalonados han aumentado debido al desarrollo de nuevas técnicas y materiales que permiten su construcción de manera rápida y económica (Concreto compactado con rodillo CCR, Gaviones, etc.). Actualmente, los canales escalonados se usan como vertedores y/o canales para peces en presas y diques, como disipadores de energía en canales y ríos, o como aireadores en plantas de tratamiento y torrentes contaminados. Diversos investigadores han estudiado el flujo en vertedores escalonados, enfocándose en estructuras de gran pendiente (  45o) por lo que a la fecha, el comportamiento del flujo sobre vertedores con pendientes moderadas (  15 a 30o) no ha sido totalmente comprendido. El presente artículo comprende un estudio experimental de las propiedades físicas del flujo aire-agua sobre canales escalonados con pendientes moderadas, típicas en presas de materiales sueltos. Un extenso rango de gastos en condiciones de flujo rasante se investigó en dos modelos experimentales a gran escala (Le = 3 a 6): Un canal con pendiente 3.5H:1V (  16o) y dos alturas de escalón distintas (h = 0.1 y 0.05 m) y un canal con pendiente 2.5H:1V (  22o) y una altura de escalón de h = 0.1 m. Los resultados incluyen un análisis detallado de las propiedades del flujo en vertedores escalonados con pendientes moderadas y un nuevo criterio de diseño hidráulico, el cual está basado en los resultados experimentales obtenidos. English abstract: Stepped chutes have been used as hydraulic structures since antiquity, they can be found acting as spillways and fish ladders in dams and weirs, as energy dissipators in artificial channels, gutters and rivers, and as aeration enhancers in water treatment plants and polluted streams. In recent years, new construction techniques and materials (Roller Compacted Concrete RCC, rip-rap gabions, etc.) together with the development of the abovementioned new applications have allowed cheaper construction methods, increasing the interest in stepped chute design. During the last three decades, research in stepped spillways has been very active. However, studies prior to 1993 neglected the effect of free-surface aeration. A number of studies have focused since on steep stepped chutes (  45o) but the hydraulic performance of moderate-slope stepped channels is not yet totally understood. This study details an experimental investigation of physical air-water flow properties down moderate slope stepped spillways conducted in two laboratory models: the first model was a 3.15 m long stepped chute with a 15.9o slope comprising two interchangeable step heights (h = 0.1 m and h = 0.05 m); the second model was a 3.3 m long, stepped channel with a 21.8o slope (h = 0.1 m). A broad range of discharges within transition and skimming flow regimes was investigated. Measurements were conducted using a double tip conductivity probe. The study provides new, original insights into air-water stepped chute flows not foreseen in prior studies and presents a new design criterion for chutes with moderate slopes based on the experimental results.

Experimental Analysis of Froude Number Effect on Air Entrainment in the Hydraulic Jump

A hydraulic jump is characterized by strong energy dissipation and mixing, large-scale turbulence, air entrainment, waves and spray. Despite recent pertinent studies, the interaction between air bubbles diffusion and momentum transfer is not completely understood. The objective of this paper is to present experimental results from new measurements performed in rectangular horizontal flume with partially-developed inflow conditions. The vertical distributions of void fraction and air bubbles count rate were recorded for inflow Froude number Fr1 in the range from 5.2 to 14.3. Rapid detrainment process was observed near the jump toe, whereas the structure of the air diffusion layer was clearly observed over longer distances. These new data were compared with previous data generally collected at lower Froude numbers. The comparison demonstrated that, at a fixed distance from the jump toe, the maximum void fraction Cmax increases with the increasing Fr1. The vertical locations of the maximum void fraction and bubble count rate were consistent with previous studies. Finally, an empirical correlation between the upper boundary of the air diffusion layer and the distance from the impingement point was provided.