skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Minimum Hydraulic Resistance Uncertainty and the Development of a Connectivity‐Based Iterative Sampling Strategy
Abstract The presence of well‐connected paths is commonly observed in spatially heterogeneous porous formations. Channels consisting of high hydraulic conductivity (K) values strongly affect fate and transport of dissolved species in the subsurface environment. Several studies have established a correlation between connectivity properties of the spatially variableK‐field and solute first arrival times. However, due to limited knowledge of the spatial structure of theK‐field, connectivity metrics are subject to uncertainty. In this work, we utilize the concept of the minimum hydraulic resistance and least resistance path to evaluate the connectivity of aK‐field in a stochastic framework. We employ a fast graph theory‐based algorithm to alleviate the computational burden associated with stochastic computations in order to investigate both the impact of the hydrogeological structural conceptualization and domain dimensionality (2‐D vs. 3‐D) on the uncertainty of the minimum hydraulic resistance. Finally, we propose an iterative data acquisition strategy that can be utilized to identify the least resistance path (which is linked to preferential flow channels) in real sites. A synthetic benchmark test is presented, showing the advantages of the proposed sampling strategy when compared to a regular sampling strategy. By using the iterative data sampling strategy, we were able to reduce first arrival time uncertainty by 47% (when compared to the regular sampling strategy), while maintaining site characterization efforts constant.  more » « less
Award ID(s):
1654009
PAR ID:
10456898
Author(s) / Creator(s):
 ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Water Resources Research
Volume:
55
Issue:
7
ISSN:
0043-1397
Page Range / eLocation ID:
p. 5593-5611
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Summary Wood anatomical traits shape a xylem segment’s hydraulic efficiency and resistance to embolism spread due to declining water potential. It has been known for decades that variations in conduit connectivity play a role in altering xylem hydraulics. However, evaluating the precise effect of conduit connectivity has been elusive. The objective here is to establish an analytical linkage between conduit connectivity and grouping and tissue‐scale hydraulics.It is hypothesized that an increase in conduit connectivity brings improved resistance to embolism spread due to increased hydraulic pathway redundancy. However, an increase in conduit connectivity could also reduce resistance due to increased speed of embolism spread with respect to pressure. We elaborate on this trade‐off using graph theory, percolation theory and computational modeling of xylem. The results are validated using anatomical measurements ofAcerbranch xylem.Considering only species with vessels, increases in connectivity improve resistance to embolism spread without negatively affecting hydraulic conductivity. The often measured grouping index fails to capture the totality of the effect of conduit connectivity on xylem hydraulics.Variations in xylem network characteristics, such as conduit connectivity, might explain why hypothesized trends among woody species, such as the ‘safety‐efficiency’ trade‐off hypothesis, are weaker than expected. 
    more » « less
  2. ABSTRACT Identifying the physiological mechanisms by which plants are adapted to drought is critical to predict species responses to climate change. We measured the responses of leaf hydraulic and stomatal conductances (Kleafandgs, respectively) to dehydration, and their association with anatomy, in seven species of CaliforniaCeanothusgrown in a common garden, including some of the most drought‐tolerant species in the semi‐arid flora. We tested for matching of maximum hydraulic supply and demand and quantified the role of decline ofKleafin driving stomatal closure. AcrossCeanothusspecies, maximumKleafandgswere negatively correlated, and bothKleafandgsshowed steep declines with decreasing leaf water potential (i.e., a high sensitivity to dehydration). The leaf water potential at 50% decline ingswas linked with a low ratio of maximum hydraulic supply to demand (i.e., maximumKleaf:gs). This sensitivity ofgs, combined with low minimum epidermal conductance and water storage, could contribute to prolonged leaf survival under drought. The specialized anatomy of subg.Cerastesincludes trichomous stomatal crypts and pronounced hypodermis, and was associated with higher water use efficiency and water storage. Combining our data with comparative literature of other California species, species of subg. Cerastesshow traits associated with greater drought tolerance and reliance on leaf water storage relative to other California species. In addition to drought resistance mechanisms such as mechanical protection and resistance to embolism, drought avoidance mechanisms such as sensitive stomatal closure could contribute importantly to drought tolerance in dry‐climate adapted species. 
    more » « less
  3. Abstract A coupling of random walkers on the same finite graph, who take turns sequentially, is said to be anavoidance couplingif the walkers never collide. Previous studies of these processes have focused almost exclusively on complete graphs, in particular how many walkers an avoidance coupling can include. For other graphs, apart from special cases, it has been unsettled whether even two noncolliding simple random walkers can be coupled. In this article, we construct such a coupling on (i) anyd‐regular graph avoiding a fixed subgraph depending ond; and (ii) any square‐free graph with minimum degree at least three. A corollary of the first result is that a uniformly random regular graph onnvertices admits an avoidance coupling with high probability. 
    more » « less
  4. Abstract To increase diversity and realism, surface bidirectional scattering distribution functions (BSDFs) are often modelled as consisting of multiple layers, but accurately evaluating layered BSDFs while accounting for all light transport paths is a challenging problem. Recently, Guoet al. [GHZ18] proposed an accurate and general position‐free Monte Carlo method, but this method introduces variance that leads to longer render time compared to non‐stochastic layered models. We improve the previous work by presenting two new sampling strategies,pair‐product samplingandmultiple‐product sampling. Our new methods better take advantage of the layered structure and reduce variance compared to the conventional approach of sequentially sampling one BSDF at a time. Ourpair‐product samplingstrategy importance samples the product of two BSDFs from a pair of adjacent layers. We further generalize this tomultiple‐product sampling, which importance samples the product of a chain of three or more BSDFs. In order to compute these products, we developed a new approximate Gaussian representation of individual layer BSDFs. This representation incorporates spatially varying material properties as parameters so that our techniques can support an arbitrary number of textured layers. Compared to previous Monte Carlo layering approaches, our results demonstrate substantial variance reduction in rendering isotropic layered surfaces. 
    more » « less
  5. A numerical inverse method called FlowPaths is presented to solve for the hydraulic conductivity field of an isotropic heterogeneous porous medium from a known specific discharge field (and constant-head boundary conditions). This method makes possible a new approach to reactive transport experiments, aimed at understanding the dynamic spatial and temporal evolution of hydraulic conductivity, which simultaneously record the evolving reaction and the evolving flow geometry. This inverse method assumes steady, two-dimensional flow through a square matrix of grid blocks. A graph-theoretical approach is used to find a set of flow paths through the porous medium using the known components of the specific discharge, where every vertex is traversed by at least one path from the upstream high-head boundary to the downstream low-head boundary. Darcy’s law is used to create an equation for the unknown head drop across each edge. Summation of these edge equations along each path through the network generates a set of linearly independent head-drop equations that is solved directly for the hydraulic conductivity field. FlowPaths is verified by generating 12,740 hydraulic conductivity fields of varying size and heterogeneity, calculating the corresponding specific discharge field for each, and then using that specific discharge field to estimate the underlying hydraulic conductivity field. When estimates from FlowPaths are compared to the simulated hydraulic conductivity fields, the inverse method is demonstrated to be accurate and numerically stable. Accordingly, within certain limitations, FlowPaths can be used in field or laboratory applications to find hydraulic conductivity from a known velocity field. 
    more » « less