Linking physiological processes with mangrove forest structure: phosphorus deficiency limits canopy development, hydraulic conductivity and photosynthetic carbon gain in dwarf Rhizophora mangle

Spatial gradients in mangrove tree height in barrier islands of Belize are associated with nutrient deficiency and sustained flooding in the absence of a salinity gradient. While nutrient deficiency is likely to affect many parameters, here we show that addition of phosphorus (P) to dwarf mangroves stimulated increases in diameters of xylem vessels, area of conductive xylem tissue and leaf area index (LAI) of the canopy. These changes in structure were consistent with related changes in function, as addition of P also increased hydraulic conductivity (K-s), stomatal conductance and photosynthetic assimilation rates to the same levels measured in taller trees fringing the seaward margin of the mangrove. Increased xylem vessel size and corresponding enhancements in stern hydraulic conductivity in P fertilized dwarf trees came at the cost of enhanced midday loss of hydraulic conductivity and was associated with decreased assimilation rates in the afternoon. Analysis of trait plasticity identifies hydraulic properties of trees as more plastic than those of leaf structural and physiological characteristics, implying that hydraulic properties are key in controlling growth in mangroves. Alleviation of P deficiency, which released trees from hydraulic limitations, reduced the structural and functional distinctions between dwarf and taller fringing tree forms of Rhizophora mangle.

Variation in hydraulic conductivity of mangroves: influence of species, salinity, and nitrogen and phosphorus availability

We investigated how species identity and variation in salinity and nutrient availability influence the hydraulic conductivity of mangroves. Using a fertilization study of two species in Florida, we found that stem hydraulic conductivity expressed on a leaf area basis (K-leaf) was significantly different among species of differing salinity tolerance, but was not significantly altered by enrichment with limiting nutrients. Reviewing data from two additional sites (Panama and Belize), we found an overall pattern of declining leaf-specific hydraulic conductivity (K-leaf) with increasing salinity. Over three sites, a general pattern emerges, indicating that native stem hydraulic conductivity (K-h) and K-leaf are less sensitive to nitrogen (N) fertilization when N limits growth, but more sensitive to phosphorus (P) fertilization when P limits growth. Processes leading to growth enhancement with N fertilization are probably associated with changes in allocation to leaf area and photosynthetic processes, whereas water uptake and transport processes could be more limiting when P limits growth. These findings suggest that whereas salinity and species identity place broad bounds on hydraulic conductivity, the effects of nutrient availability modulate hydraulic conductivity and growth in complex ways.

Effect of cultivation on soil C contents and saturated hydraulic conductivity

The irrigation of pasture with saline, Na-contaminated industrial wastewater typically results in an increase in soil ESP. From current knowledge (derived largely from cultivated agricultural soils), although these sodic soils are likely to remain stable whilst irrigated with effluent (due to the effluent’s large electrolyte concentration), during rainfall periods of low electrolyte concentration these soils would be expected to disperse. However, effluent irrigated pasture soils have been observed to maintain their structure even during intense rainfall events. Three soil types were collected (Sodosol, Vertosol and Dermosol), each with a cultivated/non-cultivated pair. The soils were equilibrated with various SAR solutions and then leached with deionised water to allow the measurement of saturated hydraulic conductivity (Ksat). At low SARs, Ksat tended to be greater in non-cultivated than cultivated soils and is attributable to a loss of structure associated with cultivation. In addition, as SAR increased, the reduction in relative Ksat tended to be significantly greater in cultivated than non-cultivated soils. The relatively rapid saturated hydraulic conductivity in the non-cultivated soils at large SARs is due to a greater aggregate stability due to greater soil C content. For the sustainable disposal of saline effluent, it is therefore necessary to ensure that soils remain undisturbed and preferably under pasture, thus maximising soil structural stability and hydraulic conductivity.

Hydraulic properties of layered soils influence survival of Rhodes grass (Chloris gayana Kunth.) during water stress

Survival of vegetation on soil-capped mining wastes is often impaired during dry seasons due to the limited amount of water stored in the shallow soil capping. Growth and survival of Rhodes grass (Chloris gayana) during soil drying on various layered capping sequences constructed of combinations of topsoil, subsoil, seawater-neutralised residue sand and low grade bauxite was determined in a glasshouse. The aim was to describe the survival of Rhodes grass in terms of plant and soil water relationships. The soil water characteristic curve and soil texture analysis was a good predictor of plant survival. The combination of soil with a high water holding capacity and low soil water diffusivity (e.g. subsoil with high clay contents) with soil having a high water holding capacity and high diffusivity (e.g. residue sand) gave best survival during drying down (up to 88 days without water), whereas topsoil and low grade bauxite were unsuitable (plants died within 18-39 days). Clayey soil improved plant survival by triggering a water stress response during peak evaporative water demand once residue sand dried down and its diffusivity fell below a critical range. Thus, for revegetation in seasonally dry climates, soil capping should combine one soil with low diffusivity and one or more soils with high total water holding capacity and high diffusivity.

A new analytical solution for water table fluctuations in coastal aquifers with sloping beaches

The Boussinesq equation appears as the zeroth-order term in the shallow water flow expansion of the non-linear equation describing the flow of fluid in an unconfined aquifer. One-dimensional models based on the Boussinesq equation have been used to analyse tide-induced water table fluctuations in coastal aquifers. Previous analytical solutions for a sloping beach are based on the perturbation parameter, epsilon(N) = alphaepsilon cot beta (in which beta is the beach slope, alpha is the amplitude parameter and epsilon is the shallow water parameter) and are limited to tan(-1) (alphaepsilon) much less than beta less than or equal to pi/2. In this paper, a new higher-order solution to the non-linear boundary value problem is derived. The results demonstrate the significant influence of the higher-order components and beach slope on the water table fluctuations. The relative difference between the linear solution and the present solution increases as 6 and a increase, and reaches 7% of the linear solution. (C) 2003 Elsevier Ltd. All rights reserved.

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated. (c) 2004 Elsevier Ltd. All rights reserved.

Boron nutrition and chilling tolerance of warm climate crop species

Background Field observations and glasshouse studies have suggested links between boron (B)-deficiency and leaf damage induced by low temperature in crop plants, but causal relationships between these two stresses at physiological, biochemical and molecular levels have yet to be explored. Limited evidence at the whole-plant level suggests that chilling temperature in the root zone restricts B uptake capacity and/or B distribution/utilization efficiency in the shoot, but the nature of this interaction depends on chilling tolerance of species concerned, the mode of low temperature treatment (abrupt versus gradual temperature decline) and growth conditions (e.g. photon flux density and relative humidity) that may exacerbate chilling stress. Scope This review explores roles of B nutrition in chilling tolerance of continual root or transient shoot chills in crop species adapted to warm season conditions. It reviews current research on combined effects of chilling temperature (ranging from > 0 to 20 degrees C) and B deficiency on growth and B nutrition responses in crop species differing in chilling tolerance. Conclusion For subtropical/tropical species (e.g. cucumber, cassava, sunflower), root chilling at 10-17 degrees C decreases B uptake efficiency and B utilization in the shoot and increases the shoot : root ratio, but chilling-tolerant temperate species (e.g. oilseed rape, wheat) require much lower root chill temperatures (2-5 degrees C) to achieve the same responses. Boron deficiency exacerbates chilling injuries in leaf tissues, particularly under high photon flux density. Suggested mechanisms for B x chilling interactions in plants are: (a) chilling-induced reduction in plasmalemma hydraulic conductivity, membrane fluidity, water channel activity and root pressure, which contribute to the decrease in root hydraulic conductance, water uptake and associated B uptake; (b) chilling-induced stomatal dysfunction affecting B transport from root to shoot and B partitioning in the shoot; and (c) B deficiency induced sensitivity to photo-oxidative damage in leaf cells. However, specific evidence for each of the mechanisms is still lacking. Impacts of B status on chilling tolerance in crop species have important implications for the management of B supply during sensitive stages of growth, such as early growth after planting and early reproductive development, both of which can coincide with the occurrence of chilling temperatures in the field.

Estimation of water retention curve from mercury intrusion porosimetry and van Genuchten model

The water retention curve (WRC) is a hydraulic characteristic of concrete required for advanced modeling of water (and thus solute) transport in variably saturated, heterogeneous concrete. Unfortunately, determination by a direct experimental method (for example, measuring equilibrium moisture levels of large samples stored in constant humidity cells) is a lengthy process, taking over 2 years for large samples. A surrogate approach is presented in which the WRC is conveniently estimated from mercury intrusion porosimetry (MIP) and validated by water sorption isotherms: The well-known Barrett, Joyner and Halenda (BJH) method of estimating the pore size distribution (PSD) from the water sorption isotherm is shown to complement the PSD derived from conventional MIP. This provides a basis for predicting the complete WRC from MIP data alone. The van Genuchten equation is used to model the combined water sorption and MIP results. It is a convenient tool for describing water retention characteristics over the full moisture content range. The van Genuchten parameter estimation based solely on MIP is shown to give a satisfactory approximation to the WRC, with a simple restriction on one. of the parameters.

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.