Search for: All records

In a complex ecohydrologic system, vegetation and soil variables combine to dictate heat fluxes, and these fluxes may vary depending on the extent to which drivers are linearly or nonlinearly interrelated. From a modeling and causality perspective, uncertainty, sensitivity, and performance measures all relate to how information from different sources “flows” through a model to produce a target, or output. We address how model structure, broadly defined as a mapping from inputs to an output, combines with source dependencies to produce a range of information flow pathways from sources to a target. We apply information decomposition, which partitions reductions in uncertainty into synergistic, redundant, and unique information types, to a range of model cases. Toy models show that model structure and source dependencies both restrict the types of interactions that can arise between sources and targets. Regressions based on weather data illustrate how different model structures vary in their sensitivity to source dependencies, thus affecting predictive and functional performance. Finally, we compare the Surface Flux Equilibrium theory, a land‐surface model, and neural networks in estimating the Bowen ratio and find that models trade off information types particularly when sources have the highest and lowest dependencies. Overall, this study extends an information theory‐based model evaluation framework to incorporate the influence of source dependency on information pathways. This could be applied to explore behavioral ranges for both machine learning and process‐based models, and guide model development by highlighting model deficiencies based on information flow pathways that would not be apparent based on existing measures.

 

At the edge of alpine and Arctic ecosystems all over the world, a transition zone exists beyond which it is either infeasible or unfavorable for trees to exist, colloquially identified as the treeline. We explore the possibility of a thermodynamic basis behind this demarcation in vegetation by considering ecosystems as open systems driven by thermodynamic advantage—defined by vegetation’s ability to dissipate heat from the earth’s surface to the air above the canopy. To deduce whether forests would be more thermodynamically advantageous than existing ecosystems beyond treelines, we construct and examine counterfactual scenarios in which trees exist beyond a treeline instead of the existing alpine meadow or Arctic tundra. Meteorological data from the Italian Alps, United States Rocky Mountains, and Western Canadian Taiga-Tundra are used as forcing for model computation of ecosystem work and temperature gradients at sites on both sides of each treeline with and without trees. Model results indicate that the alpine sites do not support trees beyond the treeline, as their presence would result in excessive CO$$_2$$${}_2$loss and extended periods of snowpack due to temperature inversions (i.e., positive temperature gradient from the earth surface to the atmosphere). Further, both Arctic and alpine sites exhibit negative work resulting in positive feedback between vegetation heat dissipation and temperature gradient, thereby extending the duration of temperature inversions. These conditions demonstrate thermodynamic infeasibility associated with the counterfactual scenario of trees existing beyond a treeline. Thus, we conclude that, in addition to resource constraints, a treeline is an outcome of an ecosystem’s ability to self-organize towards the most advantageous vegetation structure facilitated by thermodynamic feasibility.

 

Vegetation optimizes its geochemical environment for resource management via root exudation. We refer to the soil zone where biogeochemical behavior is significantly influenced, directly or indirectly, by root processes as the vegetation induced reactive zone (VIRZ). Root exudates react with VIRZ soil substrates creating temporally variable chemical environments through depth that extend below the rooting zone, impacting weathering, and releasing solutes and gases. We present a new framework, REWTCrunch, to capture VIRZ dynamics by integrating three modeling advances: the multicomponent reactive transport model CrunchFlow, the root exudation model REWT, and the multilayer canopy‐root ecohydrologic model MLCan. REWTCrunch's high‐resolution, process‐based simulation of root exudation, and the transport and transformation of carbon (C) and nutrients according to mass‐balanced and charge‐balanced reaction networks gives new insight into vertically resolved root‐soil‐microbe‐water interactions and their influence on solute fluxes at a daily timescale. We benchmark REWT and CrunchFlow, illustrate coupling mechanisms, and present REWTCrunch simulations for an agricultural site in the US Midwest. Results demonstrate root‐sourced reactive C can augment or reduce solute concentrations in the soil by several orders of magnitude. Silicate weathering products illustrate after‐harvest effects of plant C inputs in leaching patterns. Calcium simulations reveal the development of a stable weathering front. Aluminum concentrations are particularly responsive to root‐sourced reactivity, and analysis of leaching concentration versus leaching flux indicates hysteresis behavior. REWTCrunch significant improves our ability to simulate the link between root processes and soil biogeochemistry, thereby filling an important gap in the numerical simulation of root processes, weathering, and long‐term soil health.

 

Threshold changes in rainfall‐runoff generation commonly represent shifts in runoff mechanisms and hydrologic connectivity controlling water and solute transport and transformation. In watersheds with limited human influence, threshold runoff responses reflect interaction between precipitation event and antecedent soil moisture. Similar analyses are lacking in intensively managed landscapes where installation of subsurface drainage tiles has altered connectivity between the land surface, groundwater, and streams, and where application of fertilizer has created significant stores of subsurface nitrogen. In this study, we identify threshold patterns of tile‐runoff generation for a drained agricultural field in Illinois and evaluate how antecedent conditions—including shallow soil moisture, groundwater table depth, and the presence or absence of crops—control tile response. We relate tile‐runoff thresholds to patterns of event nitrate load observed across multiple storm events and evaluate how antecedent conditions control within‐event nitrate concentration‐discharge relationships. Our results demonstrate that an event tile‐runoff threshold emerges relative to the sum of gross precipitation and indices of antecedent shallow soil moisture and antecedent below‐tile groundwater moisture deficit, indicating that both shallow soil and below‐tile storages must be filled to generate significant runoff. In turn, event nitrate load shows a linear dependence on runoff for most time periods, suggesting that subsurface nitrate export and storage can be estimated using runoff threshold relationships and long‐term average nitrate concentrations. Finally, within‐event nitrate concentration‐discharge relationships are controlled by event size and the antecedent tile flow state because these factors dictate the sequence of flow path activation and tile connectivity over a storm event.

 

Hydraulic redistribution is the transport of water from wet to dry soil layers, upward or downward, through plant roots. Often in savanna and woodland ecosystems, deep‐rooted trees, and shallow‐rooted grasses coexist. The degree to which these different species compete for or share soil‐water derived from precipitation or groundwater, as well as how these interactions are altered by hydraulic redistribution, is unknown. We use a multilayer canopy model and field observations to examine how the presence of deep, but tree‐root accessible, groundwater impacts seasonal patterns of hydraulic redistribution, and interaction between coexisting vegetation species in a semiarid riparian woodland (US‐CMW). Based on the simulation, trees absorb moisture at the water table (∼10 m depth) and release it in the shallow soil depth (0–3 m) during the dry pre‐monsoon season. We observed the occurrence of a new convergent hydraulic redistribution pattern during the monsoon season, where moisture is transported from both the near‐surface (0–0.5 m) and the water table to intermediate soil layers (1–5 m) through tree roots. We found that hydraulic redistribution demonstrates a growth facilitation effect at this site, supporting 49% of growing season tree transpiration and 14% of the grass transpiration. Compared to a similarly structured upland savanna without accessible groundwater, the riparian site shows an increased amount of hydraulically redistributed water and more facilitative water use between coexisting grasses and trees. These results shed light on the linkage between accessible groundwater and the role of hydraulic redistribution on the interaction between deep‐rooted and shallow‐rooted vegetation.

 

At the biosphere–atmosphere interface, nonlinear interdependencies among components of an ecohydrological complex system can be inferred using multivariate high frequency time series observations. Information flow among these interacting variables allows us to represent the causal dependencies in the form of a directed acyclic graph (DAG). We use high frequency multivariate data at 10 Hz from an eddy covariance instrument located at 25 m above agricultural land in the Midwestern US to quantify the evolutionary dynamics of this complex system using a sequence of DAGs by examining the structural dependency of information flow and the associated functional response. We investigate whether functional differences correspond to structural differences or if there are no functional variations despite the structural differences. We base our analysis on the hypothesis that causal dependencies are instigated through information flow, and the resulting interactions sustain the dynamics and its functionality. To test our hypothesis, we build upon causal structure analysis in the companion paper to characterize the information flow in similarly clustered DAGs from 3-min non-overlapping contiguous windows in the observational data. We characterize functionality as the nature of interactions as discerned through redundant, unique, and synergistic components of information flow. Through this analysis, we find that in turbulence at the biosphere–atmosphere interface, the variables that control the dynamic character of the atmosphere as well as the thermodynamics are driven by non-local conditions, while the scalar transport associated with CO2 and H2O is mainly driven by short-term local conditions. 

Eddy covariance measurements quantify the magnitude and temporal variability of land-atmosphere exchanges of water, heat, and carbon dioxide (CO 2 ) among others. However, they also carry information regarding the influence of spatial heterogeneity within the flux footprint, the temporally dynamic source/sink area that contributes to the measured fluxes. A 25 m tall eddy covariance flux tower in Central Illinois, USA, a region where drastic seasonal land cover changes from intensive agriculture of maize and soybean occur, provides a unique setting to explore how the organized heterogeneity of row crop agriculture contributes to observations of land-atmosphere exchange. We characterize the effects of this heterogeneity on latent heat ( LE ), sensible heat ( H ), and CO 2 fluxes ( F c ) using a combined flux footprint and eco-hydrological modeling approach. We estimate the relative contribution of each crop type resulting from the structured spatial organization of the land cover to the observed fluxes from April 2016 to April 2019. We present the concept of a fetch rose, which represents the frequency of the location and length of the prevalent upwind distance contributing to the observations. The combined action of hydroclimatological drivers and land cover heterogeneity within the dynamic flux footprint explain interannual flux variations. We find that smaller flux footprints associated with unstable conditions are more likely to be dominated by a single crop type, but both crops typically influence any given flux measurement. Meanwhile, our ecohydrological modeling suggests that land cover heterogeneity leads to a greater than 10% difference in flux magnitudes for most time windows relative to an assumption of equally distributed crop types. This study shows how the observed flux magnitudes and variability depend on the organized land cover heterogeneity and is extensible to other intensively managed or otherwise heterogeneous landscapes. 

The migration timing of Pacific salmon in the Columbia River basin is subject to multiple influences related to climate, human water resource management, and lagged effects such as oceanic conditions. We apply an information theory-based approach to analyze drivers of adult Chinook salmon migration within the spring and fall spawning seasons and between years based on salmon counts at dams along the Columbia and Snake Rivers. Time-lagged mutual information and information decomposition measures, which characterize lagged and nonlinear dependencies as reductions in uncertainty, are used to detect interactions between salmon counts and lagged streamflows, air and water temperatures, precipitation, snowpack, climate indices and downstream salmon counts. At a daily timescale, these interdependencies reflect migration timing and show differences between fall and spring run salmon, while dependencies based on variables at an annual resolution reflect long-term predictability. We also highlight several types of joint dependencies where predictability of salmon counts depends on the knowledge of multiple lagged sources. This study illustrates how co-varying human and natural drivers could propagate to influence salmon migration timing or overall returns, and how nonlinear types of dependencies between variables enhance predictability of a target. This information-based framework is broadly applicable to assess driving factors in other types of complex water resources systems or species life cycles.