Finger-like protrusions that form along fluid–fluid displacement fronts in porous media are often excited by hydrodynamic instability when low-viscosity fluids displace high-viscosity resident ...fluids. Such interfacial instabilities are undesirable in many natural and engineered displacement processes. We report a phenomenon whereby gradual and monotonic variation of pore sizes along the front path suppresses viscous fingering during immiscible displacement, that seemingly contradicts conventional expectation of enhanced instability with pore size variability. Experiments and porescale numerical simulations were combined with an analytical model for the characteristics of displacement front morphology as a function of the pore size gradient. Our results suggest that the gradual reduction of pore sizes act to restrain viscous fingering for a predictable range of flow conditions (as anticipated by gradient percolation theory). The study provides insights into ways for suppressing unwanted interfacial instabilities in porous media, and provides design principles for new engineered porous media such as exchange columns, fabric, paper, and membranes with respect to their desired immiscible displacement behavior.
Soil bacterial communities are central to ecosystem functioning and services, yet spatial variations in their composition and diversity across biomes and climatic regions remain largely unknown. We ...employ multivariate general additive modeling of recent global soil bacterial datasets to elucidate dependencies of bacterial richness on key soil and climatic attributes. Although results support the well-known association between bacterial richness and soil pH, a hierarchy of novel covariates offers surprising new insights. Defining climatic soil water content explains both, the extent and connectivity of aqueous micro-habitats for bacterial diversity and soil pH, thus providing a better causal attribution. Results show that globally rare and abundant soil bacterial phylotypes exhibit different levels of dependency on environmental attributes. Surprisingly, the strong sensitivity of rare bacteria to certain environmental conditions improves their predictability relative to more abundant phylotypes that are often indifferent to variations in environmental drivers.
Spatial self-organization is a hallmark of surface-associated microbial communities that is governed by local environmental conditions and further modified by interspecific interactions. Here, we ...hypothesize that spatial patterns of microbial cell-types can stabilize the composition of cross-feeding microbial communities under fluctuating environmental conditions. We tested this hypothesis by studying the growth and spatial self-organization of microbial co-cultures consisting of two metabolically interacting strains of the bacterium Pseudomonas stutzeri. We inoculated the co-cultures onto agar surfaces and allowed them to expand (i.e. range expansion) while fluctuating environmental conditions that alter the dependency between the two strains. We alternated between anoxic conditions that induce a mutualistic interaction and oxic conditions that induce a competitive interaction. We observed co-occurrence of both strains in rare and highly localized clusters (referred to as "spatial jackpot events") that persist during environmental fluctuations. To resolve the underlying mechanisms for the emergence of spatial jackpot events, we used a mechanistic agent-based mathematical model that resolves growth and dispersal at the scale relevant to individual cells. While co-culture composition varied with the strength of the mutualistic interaction and across environmental fluctuations, the model provides insights into the formation of spatially resolved substrate landscapes with localized niches that support the co-occurrence of the two strains and secure co-culture function. This study highlights that in addition to spatial patterns that emerge in response to environmental fluctuations, localized spatial jackpot events ensure persistence of strains across dynamic conditions.
Abstract
Evidence suggests that the response of rainfed crops to dry or wet years is modulated by soil texture. This is a central tenet for certain agronomic operations in water-limited regions that ...rely on spatial distribution of soil texture for guiding precision agriculture. In contrast, natural vegetation in climatic equilibrium evolves to form a dynamic assemblage of traits and species adapted to local climatic conditions, primarily precipitation in water-limited regions. For undisturbed landscapes, we hypothesize that natural vegetation responds to rainfall anomalies irrespectively of local soil texture whereas rainfed crops are expected to respond to texture-mediated plant available water. Earth system models (ESMs) often quantify vegetation response to drought and water stress based on traditional agronomic concepts despite fundamental differences in composition and traits of natural vegetation and crops. We seek to test the hypothesis above at local and regional scales to differentiate natural vegetation and rainfed crops response to rainfall anomalies across soil types and better link them to water and carbon cycles. We employed field observations and remote sensing data to systematically examine the response of natural and rainfed cropped vegetation across biomes and scales. At local scales (field to ∼0.1 km), we used crop yields from literature data and natural vegetation productivity as gross primary productivity (GPP) from adjacent FLUXNET sites. At regional scales (∼10
2
km), we rely exclusively on remote-sensing-based GPP. Results confirm a lack of response of natural vegetation productivity to soil texture across biomes and rainfall anomalies at all scales. In contrast, crop yields at field scale exhibit correlation with soil texture in dry years (in agreement with conventional agronomic practices). These results support the hypothesis that natural vegetation is decoupled from soil texture, whereas rainfed crops retain dependency on soil texture in dry years. However, the observed correlation of crops with soil texture becomes obscured at larger scales by spatial variation of topography, rainfall, and uncertainty in soil texture and GPP values. The study provides new insights into what natural vegetation’s climatic equilibrium might mean and reveals the role of scale in expressing such sensitivities in ESMs.
Leaves within a canopy may experience rapid and extreme fluctuations in ambient conditions. A shaded leaf, for example, may become exposed to an order of magnitude increase in solar radiation within ...a few seconds, due to sunflecks or canopy motions. Considering typical time scales for stomatal adjustments, (2 to 60 minutes), the gap between these two time scales raised the question whether leaves rely on their hydraulic and thermal capacitances for passive protection from hydraulic failure or over-heating until stomata have adjusted. We employed a physically based model to systematically study effects of short-term fluctuations in irradiance on leaf temperatures and transpiration rates. Considering typical amplitudes and time scales of such fluctuations, the importance of leaf heat and water capacities for avoiding damaging leaf temperatures and hydraulic failure were investigated. The results suggest that common leaf heat capacities are not sufficient to protect a non-transpiring leaf from over-heating during sunflecks of several minutes duration whereas transpirative cooling provides effective protection. A comparison of the simulated time scales for heat damage in the absence of evaporative cooling with observed stomatal response times suggested that stomata must be already open before arrival of a sunfleck to avoid over-heating to critical leaf temperatures. This is consistent with measured stomatal conductances in shaded leaves and has implications for water use efficiency of deep canopy leaves and vulnerability to heat damage during drought. Our results also suggest that typical leaf water contents could sustain several minutes of evaporative cooling during a sunfleck without increasing the xylem water supply and thus risking embolism. We thus submit that shaded leaves rely on hydraulic capacitance and evaporative cooling to avoid over-heating and hydraulic failure during exposure to typical sunflecks, whereas thermal capacitance provides limited protection for very short sunflecks (tens of seconds).
Saturated hydraulic conductivity (Ksat) is a key soil hydraulic parameter for representing infiltration and drainage in land surface models. For large scale applications, Ksat is often estimated from ...pedotransfer functions (PTFs) based on easy‐to‐measure soil properties like soil texture and bulk density. The reliance of PTFs on data from uniform arable lands and the omission of soil structure limits the applicability of texture‐based predictions of Ksat in vegetated lands. To include effects of terrain, climate, and vegetation in the derivation of a new global Ksat map at 1 km resolution, we harness technological advances in machine learning and availability of remotely sensed surrogate information. For model training and testing, a global compilation of 6,814 geo‐referenced Ksat measurements from the literature was used. The accuracy assessment based on spatial cross‐validation shows a concordance correlation coefficient (CCC) of 0.16 and a root mean square error (RMSE) of 1.18 for log10 Ksat values in cm/day (CCC = 0.79 and RMSE = 0.72 for non‐spatial cross‐validation). The generated maps of Ksat represent spatial patterns of soil formation processes more distinctly than previous global maps of Ksat based on easy‐to‐measure soil properties. The validation of the model indicates that Ksat could be modeled without bias using Covariate‐based GeoTransfer Functions (CoGTFs) that harness spatially distributed surface and climate attributes, compared to soil information based PTFs. The relatively poor performance of all models in the validation (low CCC and high RMSE) highlights the need for the collection of additional Ksat values to train the model for regions with sparse data.
Plain Language Summary
The soil saturated hydraulic conductivity (Ksat) defines how fast water infiltrates into and percolates through the soil. To model water flow at large scales, accurate maps of Ksat are needed. Usually, Ksat is not measured directly but deduced from well‐known basic soil properties (e.g., soil texture, bulk density). However, these estimates neglect the influence of vegetation and climate on formation of soil structures that control Ksat. To improve global predictions of Ksat, we use a new spatially referenced Ksat data collection and apply machine learning to exploit correlations between Ksat and other properties (e.g., soil information, terrain, climate, and vegetation). These correlations are then implemented at global scale using maps of all relevant properties (so called “environmental covariates”) that were measured by remote sensing. We call this new approach to predictive Ksat mapping “Covariate‐based GeoTransfer Function” (CoGTF) to highlight differences with other maps that neglect spatial correlation with soil formation processes and that are based only on soil data (so called “pedotransfer functions”, PTFs). We show that the new maps based on CoGTF perform better than approaches based on PTFs.
Key Points
Climate, vegetation, and terrain affect spatial patterns of saturated hydraulic conductivity (Ksat)
The effect of these environmental covariates on Ksat is quantified using remote sensing data and machine learning
We introduce Covariate‐based GeoTransfer Functions to improve Ksat predictions based on pedotransfer functions
Resource patchiness and aqueous phase fragmentation in soil may induce large differences local growth conditions at submillimeter scales. These are translated to vast differences in bacterial age ...from cells dividing every thirty minutes in close proximity to plant roots to very old cells experiencing negligible growth in adjacent nutrient poor patches. In this study, we link bacterial population demographics with localized soil and hydration conditions to predict emerging generation time distributions and estimate mean bacterial cell ages using mechanistic and heuristic models of bacterial life in soil. Results show heavy-tailed distributions of generation times that resemble a power law for certain conditions, suggesting that we may find bacterial cells of vastly different ages living side by side within small soil volumes. Our results imply that individual bacteria may exist concurrently with all of their ancestors, resulting in an archive of bacterial cells with traits that have been gained (and lost) throughout time-a feature unique to microbial life. This reservoir of bacterial strains and the potential for the reemergence of rare strains with specific functions may be critical for ecosystem stability and function.
Individual contributions of capillarity and adsorptive surface forces to the matric potential are seldom differentiated in determination of soil water characteristic (SWC) curves. Typically, ...capillary forces dominate at the wet end, whereas adsorptive surface forces dominate at the dry end of a SWC where water is held as thin liquid films. The amount of adsorbed soil water is intimately linked to soil specific surface area (SA) and plays an important role in various biological and transport processes in arid environments. Dominated by van der Waals adsorptive forces, surface-water interactions give rise to a nearly universal scaling relationship for SWC curves at low water contents. We demonstrate that scaling measured water content at the dry end by soil specific surface area yields remarkable similarity across a range of soil textures and is in good agreement with theoretical predictions based on van der Waals interactions. These scaling relationships are important for accurate description of SWC curves in dry soils and may provide rapid and reliable estimates of soil specific surface area from SWC measurements for matric potentials below -10 MPa conveniently measured with the chilled-mirror dew point technique. Surface area estimates acquired by fitting the scaling relationship to measured SWC data were in good agreement with SA data measured by standard methods. Preliminary results suggest that the proposed method could provide reliable SA estimates for natural soils with hydratable surface areas smaller than 200 m2/g.
Abstract
Earthworm activity modifies soil structure and promotes important hydrological ecosystem functions for agricultural systems. Earthworms use their flexible hydroskeleton to burrow and expand ...biopores. Hence, their activity is constrained by soil hydromechanical conditions that permit deformation at earthworm’s maximal hydroskeletal pressure (≈200kPa). A mechanistic biophysical model is developed here to link the biomechanical limits of earthworm burrowing with soil moisture and texture to predict soil conditions that permit bioturbation across biomes. We include additional constraints that exclude earthworm activity such as freezing temperatures, low soil pH, and high sand content to develop the first predictive global map of earthworm habitats in good agreement with observed earthworm occurrence patterns. Earthworm activity is strongly constrained by seasonal dynamics that vary across latitudes largely due to soil hydromechanical status. The mechanistic model delineates the potential for earthworm migration via connectivity of hospitable sites and highlights regions sensitive to climate.
We quantify mechanical processes common to soil penetration by earthworms and growing plant roots, including the energetic requirements for soil plastic displacement. The basic mechanical model ...considers cavity expansion into a plastic wet soil involving wedging by root tips or earthworms via cone-like penetration followed by cavity expansion due to pressurized earthworm hydroskeleton or root radial growth. The mechanical stresses and resulting soil strains determine the mechanical energy required for bioturbation under different soil hydro-mechanical conditions for a realistic range of root/earthworm geometries. Modeling results suggest that higher soil water content and reduced clay content reduce the strain energy required for soil penetration. The critical earthworm or root pressure increases with increased diameter of root or earthworm, however, results are insensitive to the cone apex (shape of the tip). The invested mechanical energy per unit length increase with increasing earthworm and plant root diameters, whereas mechanical energy per unit of displaced soil volume decreases with larger diameters. The study provides a quantitative framework for estimating energy requirements for soil penetration work done by earthworms and plant roots, and delineates intrinsic and external mechanical limits for bioturbation processes. Estimated energy requirements for earthworm biopore networks are linked to consumption of soil organic matter and suggest that earthworm populations are likely to consume a significant fraction of ecosystem net primary production to sustain their subterranean activities.