Water Transport Model Research Articles

Multi-response modeling of reaction-diffusion to explain alpha-galactoside behavior during the soaking-cooking process in cowpea

A modelling approach was developed to better understand the behavior of the flatulence-causing oligosaccharides in cowpea seeds during isothermal water soaking-cooking process. Concentrations of verbascose, stachyose and raffinose were measured both in the seed and in the soaking water during the process (T=30, 60 and 95°C). A reaction-diffusion model was built for the three considered alpha-galactosides both in the seed and in the soaking water, together with a model of water transport in the seed. The model reproduced coupled reaction-diffusion of alpha-galactosides during the soaking-cooking process with a good fit. Produced, diffused and degraded alpha-galactoside fractions were identified by performing a mass balance. During soaking at 30°C, degradation predominated (maximum found for raffinose degradation rate constant of 3.22×10−4s−1) whereas diffusion predominated at higher temperatures (95°C).

Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone

AbstractHydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km2 catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.

A Comprehensive Biophysical Model of Ion and Water Transport in Plant Roots. I. Clarifying the Roles of Endodermal Barriers in the Salt Stress Response.

In this paper, we present a detailed and comprehensive mathematical model of active and passive ion and water transport in plant roots. Two key features are the explicit consideration of the separate, but interconnected, apoplastic, and symplastic transport pathways for ions and water, and the inclusion of both active and passive ion transport mechanisms. The model is used to investigate the respective roles of the endodermal Casparian strip and suberin lamellae in the salt stress response of plant roots. While it is thought that these barriers influence different transport pathways, it has proven difficult to distinguish their separate functions experimentally. In particular, the specific role of the suberin lamellae has been unclear. A key finding based on our simulations was that the Casparian strip is essential in preventing excessive uptake of Na+ into the plant via apoplastic bypass, with a barrier efficiency that is reflected by a sharp gradient in the steady-state radial distribution of apoplastic Na+ across the barrier. Even more significantly, this function cannot be replaced by the action of membrane transporters. The simulations also demonstrated that the positive effect of the Casparian strip of controlling Na+ uptake, was somewhat offset by its contribution to the osmotic stress component: a more effective barrier increased the detrimental osmotic stress effect. In contrast, the suberin lamellae were found to play a relatively minor, even non-essential, role in the overall response to salt stress, with the presence of the suberin lamellae resulting in only a slight reduction in Na+ uptake. However, perhaps more significantly, the simulations identified a possible role of suberin lamellae in reducing plant energy requirements by acting as a physical barrier to preventing the passive leakage of Na+ into endodermal cells. The model results suggest that more and particular experimental attention should be paid to the properties of the Casparian strip when assessing the salt tolerance of different plant varieties and species. Indeed, the Casparian strip appears to be a more promising target for plant breeding and plant genetic engineering efforts than the suberin lamellae for the goal of improving salt tolerance.

OpenSimRoot: widening the scope and application of root architectural models.

Summary OpenSimRoot is an open‐source, functional–structural plant model and mathematical description of root growth and function. We describe OpenSimRoot and its functionality to broaden the benefits of root modeling to the plant science community.OpenSimRoot is an extended version of simroot, established to simulate root system architecture, nutrient acquisition and plant growth. OpenSimRoot has a plugin, modular infrastructure, coupling single plant and crop stands to soil nutrient and water transport models. It estimates the value of root traits for water and nutrient acquisition in environments and plant species.The flexible OpenSimRoot design allows upscaling from root anatomy to plant community to estimate the following: resource costs of developmental and anatomical traits; trait synergisms; and (interspecies) root competition. OpenSimRoot can model three‐dimensional images from magnetic resonance imaging (MRI) and X‐ray computed tomography (CT) of roots in soil. New modules include: soil water‐dependent water uptake and xylem flow; tiller formation; evapotranspiration; simultaneous simulation of mobile solutes; mesh refinement; and root growth plasticity.OpenSimRoot integrates plant phenotypic data with environmental metadata to support experimental designs and to gain a mechanistic understanding at system scales.

Leaf Hydraulic Architecture and Stomatal Conductance: A Functional Perspective.

Mechanistic modeling of water transport from petiole to stomata provides new perspectives on optimality in vascular and mesophyll transport properties.

Automatic plant-based water status monitoring in grapevine using an improved water transport and storage model

As all plants, grapevines (Vitis vinifera L.) need water to function properly. A certain level of drought stress might however be beneficial as it improves the composition of the grape. Earlier research has shown that differences in water status result in wines with different appearance, aroma, flavour and colour. Nevertheless, the level and timing is of utmost importance. Therefore, an adequate monitoring of the water status is crucial for improving grape (and wine) quality. It is internationally recognised that this should be based on plant measurements, because only then information is gained about the actual plant water status. Mechanistic models are promising for this purpose and allow a deeper understanding of the underlying mechanisms. In this study we use a dynamic water transport and storage model that links sap flow, or whole plant water consumption, and stem diameter variations in order to simulate stem water potential. This variable is considered as one of the best indicators for water status. We aimed at improving the model to perform well under both wet and pronounced dry conditions and evaluated it for real-time water status monitoring. To this end, the former constant flow resistance in the xylem has been replaced by a dynamic resistance depending on measured soil water potential and combines the resistance experienced in the soil, roots and stem. Furthermore, also the radial flow resistance (between xylem and storage tissues), originally implemented as a constant value, has been replaced by an equation. The improved model is able to accurately simulate the plant water status during both wet and pronounced dry conditions. The model seems very promising to apply as an automatic plant-based water status monitoring system and may be a tool to improve grape and wine quality.

Geometric and fractal analysis of dynamic cracking patterns subjected to wetting-drying cycles

Desiccation cracks detrimentally affect the hydraulic and mechanical properties of soils. Soil cracking dynamic subjected to wetting-drying (W-D) cycles aids in the understanding of soil shrinking-swelling mechanism and the regularity of water movement in unsaturated soils with variable solid phase. In the present paper, we experimentally investigate the geometric and fractal features of dynamic cracking patterns with the soils undergoing four W-D cycles. Crack patterns were recorded by a digital camera mounted above the soil specimen. Digital image processing and morphological algorithms were employed to yield the geometric and fractal indexes that depict crack patterns (incl. crack ratio, length, mean width, similarity degree, box-counting fractal and multifractal dimensions). The analysis of dynamic crack patterns suggest an inconsistency exhibited in the crack morphology towards the drying (DR) process and backwards the wetting (WT) process, which could be possibly attributed to the different suction stress between soil particles during moisture desorption and sorption, closely analogous to the hysteresis phenomenon of soil-water characteristic curve. The gradient of crack propagating and dissolving related to water content reaches different peak values during an individual W-D cycle. High repeatability is found across the stabilized crack patterns in the four drying cycles and the cracking degree reaches equilibrium with the increasing number of W-D cycles. The crack patterns in the subsequent drying processes are significantly dominated by the first cracking. The similarity degree between the first and respectively the second, third, fourth crack patterns decreases with the increasing W-D cycles, suggesting the soil cracked in the previously cracking vicinity, and saturated wetting/submerge hardly healed the failure zones that caused the repeatable cracking. Attempt of fractal and multifractal methods applied to the crack networks corroborates the applicability of the scaling method to the crack patterns. The analysis indicates that soil cracks progress in hierarchical networks in the first drying cycle, while the gradational crack segments nearly simultaneously emerged or vanished in the subsequent drying or wetting cycles. The study reveals the discrepancy of crack morphology between the crack formation and dissolution during W-D cycles, and demonstrates the repeatability degree and hierarchy of cracks across W-D cycles, which would aid in quantifying physical processes in cracked soil shrinking/swelling and modelling of water transport in unsaturated soils of dynamic structure.

Evaluation of Drip Irrigation System Uniformity on Cotton Yield in an Arid Region using a Two‐Dimensional Soil Water Transport and Crop Growth Coupling Model

AbstractAppropriate drip irrigation system uniformity can maintain a higher crop yield and reduce the initial investment of the system. In this study, investigating the optimal drip system uniformity for cotton, a two‐dimensional soil water transport and crop growth coupling model was calibrated and validated by soil water dynamics, crop growth indices and cotton yield obtained from mulched drip irrigation field experiments in an arid region of Xinjiang Uygur Autonomous Region, China. The model performed well with the normalized root mean square error of 2.6–22.6% for soil moisture, 6.3–14.1% for leaf area index, 4.9–7.2% for above‐ground biomass, and with the difference ranging from 2.5 to 8.2% for seed cotton yield (SCY). Next, SCYs were simulated under five Christiansen uniformity coefficients (CUs) of 55, 65, 75, 85 and 95% with three irrigation levels of 50, 75 and 100% full irrigation treatment. The predicted SCY was found to decrease as CU decreased under different irrigation treatments, and the more obvious effect of CU on SCY at 75 and 100% full irrigation treatment was obtained. In addition, SCY increased by 2.8–8.0% when CU increased from 55 to 75%, and increased by 0.1–4.6% when CU increased from 75 to 95%. Therefore, a target drip irrigation uniformity of CU = 75% was recommended to obtain a higher cotton yield, while maintaining a lower installation cost in the arid region. Copyright © 2017 John Wiley & Sons, Ltd.

Outside-Xylem Vulnerability, Not Xylem Embolism, Controls Leaf Hydraulic Decline during Dehydration.

Leaf hydraulic supply is crucial to maintaining open stomata for CO2 capture and plant growth. During drought-induced dehydration, the leaf hydraulic conductance (Kleaf) declines, which contributes to stomatal closure and, eventually, to leaf death. Previous studies have tended to attribute the decline of Kleaf to embolism in the leaf vein xylem. We visualized at high resolution and quantified experimentally the hydraulic vulnerability of xylem and outside-xylem pathways and modeled their respective influences on plant water transport. Evidence from all approaches indicated that the decline of Kleaf during dehydration arose first and foremost due to the vulnerability of outside-xylem tissues. In vivo x-ray microcomputed tomography of dehydrating leaves of four diverse angiosperm species showed that, at the turgor loss point, only small fractions of leaf vein xylem conduits were embolized, and substantial xylem embolism arose only under severe dehydration. Experiments on an expanded set of eight angiosperm species showed that outside-xylem hydraulic vulnerability explained 75% to 100% of Kleaf decline across the range of dehydration from mild water stress to beyond turgor loss point. Spatially explicit modeling of leaf water transport pointed to a role for reduced membrane conductivity consistent with published data for cells and tissues. Plant-scale modeling suggested that outside-xylem hydraulic vulnerability can protect the xylem from tensions that would induce embolism and disruption of water transport under mild to moderate soil and atmospheric droughts. These findings pinpoint outside-xylem tissues as a central locus for the control of leaf and plant water transport during progressive drought.

Field trial and modelling of different strategies for remediation of soil salinity and sodicity in the Lower Murray irrigation areas

Rising saline groundwatertables and drought in the Lower Murray Reclaimed Irrigation Area (LMRIA) has created soil salinity and sodicity conditions, which has resulted in a decline in agricultural production. A 2 month field experiment was conducted at Mobilong irrigation area using a randomised block design with trial plots (4m2) in each of the three blocks containing the following six treatments: (1) control (not irrigated), (2) irrigation (River Murray water) only, (3) gypsum application (1.5 kg/m2) and irrigation, (4) limestone application (1.5 kg/m2) and irrigation, (5) seawater (source of dissolved Ca2+) application (100mm depth) and irrigation and (6) acid (pH 3) drainage (to dissolve CaCO3 in soil to release Ca2+) application (100mm irrigation water depth) and irrigation. Soil electrical conductivity (EC), pH, exchangeable cations (Ca2+, Mg2+, Na+, K+), exchangeable sodium percentage (ESP) and saturated hydraulic conductivity were measured. The decrease in EC was greatest in the irrigation only treatment followed by the gypsum and limestone treatments. At the end of the trial, the EC in the irrigation only treatment was about one-quarter of the control, in which the EC increased. The seawater and acid drainage treatments reduced the soil EC relative to the control but the EC was 2–3dS/m higher than the irrigation, gypsum and limestone treatments at the end of the trial. The gypsum, irrigation only and limestone treatments approximately doubled the exchangeable Ca relative to the control and exchangeable Na was reduced. Unsaturated water and solute transport model (HYDRUS-UNSATCHEM) simulations were able to represent the general trends in the field results. The results suggest that only River Murray water irrigation and drainage could be effective for the management of soil salinity and sodicity in the LMRIA but further research is required to establish the threshold electrolyte concentration to prevent soil dispersion.

Measuring and modelling water transport on Skaftafellsheiði, Iceland

Areas with thick basaltic aquifers are used for drinking water supply and irrigation purposes, such as the Columbia River Basalt group in northwest USA and the Deccan Traps in India. However, rainfall-runoff processes in these basaltic areas are poorly understood. Cooling joints can transport large amounts of water, but – due to their limited porosity – they are vulnerable for over-abstraction. On Iceland, in the small Skaftafellsheiði basaltic catchment (4 x 6 km), field data were collected in 2014 and 2015. Two small streams discharge the rain surplus. Precipitation was measured at various elevations on the ridge. Also, the discharge of the streams was measured. A groundwater flow model was constructed in order to get more insight in the physical properties of the basalt aquifer and its rainfall-runoff properties. The field experiments showed that precipitation increases linearly with surface elevation. On average, the precipitation at 800 m+msl was almost double, relative to the precipitation at 200 m+msl. Calculated ETpot was rather high, due to the 19 potential sun hours per day during the Icelandic summer. Field experiments revealed quick discharge response on rainfall events, but also rather constant base flow during dryer periods. This indicates a limited infiltration capacity, but also a considerable storage capacity in the subsequent layers. The peat layer is believed to be the dominant storage/reservoir. Peat, regolith and an organic layer formed the top layer in the GMS-Modflow groundwater model. The thick basaltic aquifer was split in a series of model layers. Best results were obtained by using a decreasing hydraulic conductivity to depth. The transient model overestimated the groundwater levels at the outlet, but managed to reproduce the wet/dry conditions in the catchment rather well. This indicates that it is possible to model complex basaltic aquifers, by taking a large Representative Elementary Volume (REV) as starting point

Applying Process-Based Models for Subsurface Flow Treatment Wetlands: Recent Developments and Challenges

To date, only few process-based models for subsurface flow treatment wetlands have been developed. For modelling a treatment wetland, these models have to comprise a number of sub-models to describe water flow, pollutant transport, pollutant transformation and degradation, effects of wetland plants, and transport and deposition of suspended particulate matter. The two most advanced models are the HYDRUS Wetland Module and BIO-PORE. These two models are briefly described. This paper shows typical simulation results for vertical flow wetlands and discusses experiences and challenges using process-based wetland models in relation to the sub-models describing the most important wetland processes. It can be demonstrated that existing simulation tools can be applied for simulating processes in treatment wetlands. Most important for achieving a good match between measured and simulated pollutant concentrations is a good calibration of the water flow and transport models. Only after these calibrations have been made and the effect of the influent fractionation on simulation results has been considered, should changing the parameters of the biokinetic models be taken into account. Modelling the effects of wetland plants is possible and has to be considered when important. Up to now, models describing clogging are the least established models among the sub-models required for a complete wetland model and thus further development and research is required.

Leaf water stable isotopes and water transport outside the xylem.

How water moves through leaves, and where the phase change from liquid to vapour occurs within leaves, remain largely mysterious. Some time ago, we suggested that the stable isotope composition of leaf water may contain information on transport pathways beyond the xylem, through differences in the development of gradients in enrichment within the various pathways. Subsequent testing of this suggestion provided ambiguous results and even questioned the existence of gradients in enrichment within the mesophyll. In this review, we bring together recent theoretical developments in understanding leaf water transport pathways and stable isotope theory to map a path for future work into understanding pathways of water transport and leaf water stable isotope composition. We emphasize the need for a spatially, anatomically and isotopically explicit model of leaf water transport.

A Coupled Groundwater-Surface Water Modeling Framework for Simulating Transition Zone Processes.

There is an identified need for fully representing groundwater-surface water transition zone (i.e., the sediment zone that connects groundwater and surface water) processes in modeling fate and transport of contaminants to assist with management of contaminated sediments. Most existing groundwater and surface water fate and transport models are not dynamically linked and do not consider transition zone processes such as bioturbation and deposition and erosion of sediments. An interface module is developed herein to holistically simulate the fate and transport by coupling two commonly used models, Environmental Fluid Dynamics Code (EFDC) and SEAWAT, to simulate surface water and groundwater hydrodynamics, while providing an enhanced representation of the processes in the transition zone. Transition zone and surface water contaminant processes were represented through an enhanced version of the EFDC model, AQFATE. AQFATE also includes SEDZLJ, a state-of-the-science surface water sediment transport model. The modeling framework was tested on a published test problem and applied to evaluate field-scale two- and three-dimensional contaminant transport. The model accurately simulated concentrations of salinity from a published test case. For the field-scale applications, the model showed excellent mass balance closure for the transition zone and provided accurate simulations of all transition zone processes represented in the modeling framework. The model predictions for the two-dimensional field case were consistent with site-specific observations of contaminant migration. This modeling framework represents advancement in the simulation of transition zone processes and can help inform risk assessment at sites where contaminant sources from upland areas have the potential to impact sediments and surface water.

Prediction of mass transport in cracked-unsaturated concrete by mesoscale lattice model

The ingress of mass in a harsh chloride environment leads to the corrosion of reinforcing steel bars and deterioration of the structural performance. This paper deals with a two-dimensional (2D) numerical model for water and chlorides transport in the cracked-unsaturated concrete on the mesoscale level, in which concrete is modeled as a composite material with impermeable coarse aggregates embedded in the porous matrix which is separated by vulnerable interfacial transition zone (ITZ). Coupled with the unsaturated flow theory and the cubic law of water through a single crack, the transport coefficients for water and chloride ions in a single crack are developed by treating the crack as a parallel-plate configuration. The lattice network developed on the basis of Voronoi tessellation is presented to investigate mass transport process in the cracked-unsaturated concrete subjected to drying-wetting (D/W) cycles. The above proposed transport model coupled with the mesoscale lattice approach is validated by comparison with the available experimental findings from the literature. The numerical results indicated that the computational models are able to well represent the mass (water and chloride) movement within the cracked-unsaturated concrete. Furthermore, the cyclic D/W action and crack width/length within concrete are crucial for the transport properties of unsaturated concrete.

Analysis of the dynamic impact of Three-Gorge Project on regime of soil water and salt in Yangtze River Delta

The topography is low and flat in the estuarine region of the Yangtze River, where waterlog and salt-affected soil are all through a latent threat to the agricultural development for long time. After the Three-Gorge Project is built up, the Three-Gorge Reservoir adjusts seasonably runoff will impact on the dynamic of soil water and salt in the region. These influence relations are very complicated. Based on the theory of solute transport and the principle of soil water dynamics, the model of soil water and salt transport was created, which involved root water extraction, adsorption---desorption, mobile---immobile region. To overcome numerical dispersion, the paper solved the convection---dispersion equation by Bresler algorithm, the term of time derivative used the approximate difference formula of quadratic differential coefficient based on central difference. Then, the impact of Three-Gorge Project on the regime of soil water and salt in the estuarine region by numerical simulation was studied in. The numerical simulation indicates that electric conductivity of soil one meter thick of each observed section will go down by 0.005---0.174 mS/cm as discharge increases in spring after building the Three-Gorge Reservoir, and up by 0.004---0.416 mS/cm as discharge decreases in autumn---winter. The elementary predicted result of numerical simulation shows that salification will be bigger in autumn than desalinization in spring. The impact degree of the project on the regime of soil salt will be gradually weakened along the river up from the estuary, the nearer to the sea, the bigger the impact, and from cross profiles, the nearer to the river, the bigger the impact. Especially within the range of 4 km of the river, it is more possible that the soil will produce salification or secondary salinization.

A numerical model for water and heat transport in freezing soils with nonequilibrium ice‐water interfaces

AbstractA one‐dimensional numerical model of heat and water transport in freezing soils is developed by assuming that ice‐water interfaces are not necessarily in equilibrium. The Clapeyron equation, which is derived from a static ice‐water interface using the thermal equilibrium theory, cannot be readily applied to a dynamic system, such as freezing soils. Therefore, we handled the redistribution of liquid water with the Richard's equation. In this application, the sink term is replaced by the freezing rate of pore water, which is proportional to the extent of supercooling and available water content for freezing by a coefficient, β. Three short‐term laboratory column simulations show reasonable agreement with observations, with standard error of simulation on water content ranging between 0.007 and 0.011 cm3 cm−3, showing improved accuracy over other models that assume equilibrium ice‐water interfaces. Simulation results suggest that when the freezing front is fixed at a specific depth, deviation of the ice‐water interface from equilibrium, at this location, will increase with time. However, this deviation tends to weaken when the freezing front slowly penetrates to a greater depth, accompanied with thinner soils of significant deviation. The coefficient, β, plays an important role in the simulation of heat and water transport. A smaller β results in a larger deviation in the ice‐water interface from equilibrium, and backward estimation of the freezing front. It also leads to an underestimation of water content in soils that were previously frozen by a rapid freezing rate, and an overestimation of water content in the rest of the soils.

One-dimensional coupled model for landfill gas and water transport in layered unsaturated soil cover systems

Cover systems are used to prevent water infiltration into a waste body. They also play an important role in controlling landfill gas transport from the waste body to the atmosphere. It is important to assess the flux of landfill gas at the surface of a cover system by considering the coupled effects of rainwater infiltration and gas transport in the cover soils. We have developed a 1D mathematical model for coupled transient gas and water transport in unsaturated cover soils. The coupled model was solved by the finite element method. Results obtained by the proposed model agreed well with experimental data. Based on the proposed solution, the influences of gas pressure, gas permeability, and the thickness of the cover soils on soil gas concentration profiles were investigated. The difference in soil gas concentration reached up to 31% as the thickness of cover increased from 1 to 2 m. Gas concentration at a depth of 0.2 m decreased by 6% as the amplitude of atmospheric gas pressure fluctuation increased from 20 to 100 Pa. The gas concentration increased by only 3% when gas permeability increased by a factor of 2 for a relatively long period of gas migration (e.g., 60 h) under the given conditions. Results suggest that both diffusion and advection should be considered when estimating gas transport in unsaturated cover soils. The numerical model can be used in the design of cover systems in relation to gas breakthrough time, breakthrough concentration, and flux.

Computer-aided identification of the water diffusion coefficient for maize kernels dried in a thin layer

Abstract Uncertainties in mathematical modelling of water transport in cereal grain kernels during drying and storage are mainly due to implementing unreliable values of the water diffusion coefficient and simplifying the geometry of kernels. In the present study an attempt was made to reduce the uncertainties by developing a method for computer-aided identification of the water diffusion coefficient and more accurate 3D geometry modelling for individual kernels using original inverse finite element algorithms. The approach was exemplified by identifying the water diffusion coefficient for maize kernels subjected to drying. On the basis of the developed method, values of the water diffusion coefficient were estimated, 3D geometry of a maize kernel was represented by isoparametric finite elements, and the moisture content inside maize kernels dried in a thin layer was predicted. Validation of the results against experimental data showed significantly lower error values than in the case of results obtained for the water diffusion coefficient values available in the literature.

Role of Aquaporins in a Composite Model of Water Transport in the Leaf.

Water-transport pathways through the leaf are complex and include several checkpoints. Some of these checkpoints exhibit dynamic behavior that may be regulated by aquaporins (AQPs). To date, neither the relative weight of the different water pathways nor their molecular mechanisms are well understood. Here, we have collected evidence to support a putative composite model of water pathways in the leaf and the distribution of water across those pathways. We describe how water moves along a single transcellular path through the parenchyma and continues toward the mesophyll and stomata along transcellular, symplastic and apoplastic paths. We present evidence that points to a role for AQPs in regulating the relative weight of each path in the overall leaf water-transport system and the movement of water between these paths as a result of the integration of multiple signals, including transpiration demand, water potential and turgor. We also present a new theory, the hydraulic fuse theory, to explain effects of the leaf turgor-loss-point on water paths alternation and the subsequent reduction in leaf hydraulic conductivity. An improved understating of leaf water-balance management may lead to the development of crops that use water more efficiently, and responds better to environmental changes.