Dynamics of soil biogeochemical gas emissions shaped by remolded aggregate sizes and carbon configurations under hydration cycles

Multiscale Porosity Microfluidics to Study Bacterial Transport in Heterogeneous Chemical Landscapes

Microfluidic models are proving to be powerful systems to study fundamental processes in porous media, due to their ability to replicate topologically complex environments while allowing detailed, quantitative observations at the pore scale. Yet, while porous media such as living tissues, geological substrates, or industrial systems typically display a porosity that spans multiple scales, most microfluidic models to date are limited to a single porosity or a small range of pore sizes. Here, a novel microfluidic system with multiscale porosity is presented. By embedding polyacrylamide (PAAm) hydrogel structures through in-situ photopolymerization in a landscape of microfabricated polydimethylsiloxane (PDMS) pillars with varying spacing, micromodels with porosity spanning several orders of magnitude, from nanometers to millimeters are created. Experiments conducted at different porosity patterns demonstrate the potential of this approach to characterize fundamental and ubiquitous biological and geochemical transport processes in porous media. Accounting for multiscale porosity allows studies of the resulting heterogeneous fluid flow and concentration fields of transported chemicals, as well as the biological behaviors associated with this heterogeneity, such as bacterial chemotaxis. This approach brings laboratory studies of transport in porous media a step closer to their natural counterparts in the environment, industry, and medicine.

Consider the Anoxic Microsite: Acknowledging and Appreciating Spatiotemporal Redox Heterogeneity in Soils and Sediments

Reduction-oxidation (redox) reactions underlie essentially all biogeochemical cycles. Like most soil properties and processes, redox is spatiotemporally heterogeneous. However, unlike other soil features, redox heterogeneity has yet to be incorporated into mainstream conceptualizations of soil biogeochemistry. Anoxic microsites, the defining feature of redox heterogeneity in bulk oxic soils and sediments, are zones of oxygen depletion in otherwise oxic environments. In this review, we suggest that anoxic microsites represent a critical component of soil function and that appreciating anoxic microsites promises to advance our understanding of soil and sediment biogeochemistry. In sections 1 and 2, we define anoxic microsites and highlight their dynamic properties, specifically anoxic microsite distribution, redox gradient magnitude, and temporality. In section 3, we describe the influence of anoxic microsites on several key elemental cycles, organic carbon, nitrogen, iron, manganese, and sulfur. In section 4, we evaluate methods for identifying and characterizing anoxic microsites, and in section 5, we highlight past and current approaches to modeling anoxic microsites. Finally, in section 6, we suggest steps for incorporating anoxic microsites and redox heterogeneities more broadly into our understanding of soils and sediments.

Regulated bacterial interaction networks: A mathematical framework to describe competitive growth under inclusion of metabolite cross-feeding

When bacterial species with the same resource preferences share the same growth environment, it is commonly believed that direct competition will arise. A large variety of competition and more general 'interaction' models have been formulated, but what is currently lacking are models that link monoculture growth kinetics and community growth under inclusion of emerging biological interactions, such as metabolite cross-feeding. In order to understand and mathematically describe the nature of potential cross-feeding interactions, we design experiments where two bacterial species Pseudomonas putida and Pseudomonas veronii grow in liquid medium either in mono- or as co-culture in a resource-limited environment. We measure population growth under single substrate competition or with double species-specific substrates (substrate 'indifference'), and starting from varying cell ratios of either species. Using experimental data as input, we first consider a mean-field model of resource-based competition, which captures well the empirically observed growth rates for monocultures, but fails to correctly predict growth rates in co-culture mixtures, in particular for skewed starting species ratios. Based on this, we extend the model by cross-feeding interactions where the consumption of substrate by one consumer produces metabolites that in turn are resources for the other consumer, thus leading to positive feedback in the species system. Two different cross-feeding options were considered, which either lead to constant metabolite cross-feeding, or to a regulated form, where metabolite utilization is activated with rates according to either a threshold or a Hill function, dependent on metabolite concentration. Both mathematical proof and experimental data indicate regulated cross-feeding to be the preferred model to constant metabolite utilization, with best co-culture growth predictions in case of high Hill coefficients, close to binary (on/off) activation states. This suggests that species use the appearing metabolite concentrations only when they are becoming high enough; possibly as a consequence of their lower energetic content than the primary substrate. Metabolite sharing was particularly relevant at unbalanced starting cell ratios, causing the minority partner to proliferate more than expected from the competitive substrate because of metabolite release from the majority partner. This effect thus likely quells immediate substrate competition and may be important in natural communities with typical very skewed relative taxa abundances and slower-growing taxa. In conclusion, the regulated bacterial interaction network correctly describes species substrate growth reactions in mixtures with few kinetic parameters that can be obtained from monoculture growth experiments.

Reactive Nitrogen Hotspots Related to Microscale Heterogeneity in Biological Soil Crusts

Biocrusts covering drylands account for major fractions of terrestrial biological nitrogen fixation and release large amounts of gaseous reactive nitrogen (N_r) as nitrous acid (HONO) and nitric oxide (NO). Recent investigations suggested that aerobic and anaerobic microbial nitrogen transformations occur simultaneously upon desiccation of biocrusts, but the spatio-temporal distribution of seemingly contradictory processes remained unclear. Here, we explore small-scale gradients in chemical concentrations related to structural characteristics and organism distribution. X-ray microtomography and fluorescence microscopy revealed mixed pore size structures, where photoautotrophs and cyanobacterial polysaccharides clustered irregularly in the uppermost millimeter. Microsensor measurements showed strong gradients of pH, oxygen, and nitrite, nitrate, and ammonium ion concentrations at micrometer scales in both vertical and lateral directions. Initial oxygen saturation was mostly low (∼30%) at full water holding capacity, suggesting widely anoxic conditions, and increased rapidly upon desiccation. Nitrite concentrations (∼6 to 800 μM) and pH values (∼6.5 to 9.5) were highest around 70% WHC. During further desiccation they decreased, while emissions of HONO and NO increased, reaching maximum values around 20% WHC. Our results illustrate simultaneous, spatially separated aerobic and anaerobic nitrogen transformations, which are critical for N_r emissions, but might be impacted by future global change and land management.

Phase Saturation Control on Mixing-Driven Reactions in 3D Porous Media

Transported chemical reactions in unsaturated porous media are relevant to environmental and industrial applications. Continuum scale models are based on equivalent parameters derived from analogy with saturated conditions and cannot appropriately account for incomplete mixing. It is also unclear how the third dimension controls mixing and reactions. We obtain three-dimensional (3D) images by magnetic resonance imaging using an immiscible nonwetting liquid as a second phase and a fast irreversible bimolecular reaction. We study the impact of phase saturation on the dynamics of mixing and the reaction front. We quantify the temporally resolved effective reaction rate and describe it using the lamellar theory of mixing, which explains faster than Fickian (t^0.5) rate of product formation by accounting for the deformation of the mixing interface between the two reacting fluids. For a given Péclet, although stretching and folding of the reactive front enhance as saturation decreases, enhancing the product formation, the product formation is larger as saturation increases. After breakthrough, the extinction of the reaction takes longer as saturation decreases because of the larger nonmixed volume behind the front. These results are the basis for a general model to better predict reactive transport in unsaturated porous media not achievable by the current continuum paradigm.

Microscale pH variations during drying of soils and desert biocrusts affect HONO and NH3 emissions

Microscale interactions in soil may give rise to highly localised conditions that disproportionally affect soil nitrogen transformations. We report mechanistic modelling of coupled biotic and abiotic processes during drying of soil surfaces and biocrusts. The model links localised microbial activity with pH variations within thin aqueous films that jointly enhance emissions of nitrous acid (HONO) and ammonia (NH₃) during soil drying well above what would be predicted from mean hydration conditions and bulk soil pH. We compared model predictions with case studies in which reactive nitrogen gaseous fluxes from drying biocrusts were measured. Soil and biocrust drying rates affect HONO and NH₃ emission dynamics. Additionally, we predict strong effects of atmospheric NH₃ levels on reactive nitrogen gas losses. Laboratory measurements confirm the onset of microscale pH localisation and highlight the critical role of micro-environments in the resulting biogeochemical fluxes from terrestrial ecosystems.

The mechanisms that determine the composition of nitrogen gas emissions from soil remain unclear. A biocrust mechanistic model was developed to resolve puzzling dynamics of nitrous acid and ammonia emissions from drying soil pointing to previously unknown microscale pH zonation in thinning water films that affect soil biogeochemical fluxes.

Soil aggregates as biogeochemical reactors and implications for soil-atmosphere exchange of greenhouse gases-A concept

Soil-atmosphere exchange significantly influences the global atmospheric abundances of carbon dioxide (CO₂ ), methane (CH₄ ), and nitrous oxide (N₂ O). These greenhouse gases (GHGs) have been extensively studied at the soil profile level and extrapolated to coarser scales (regional and global). However, finer scale studies of soil aggregation have not received much attention, even though elucidating the GHG activities at the full spectrum of scales rather than just coarse levels is essential for reducing the large uncertainties in the current atmospheric budgets of these gases. Through synthesizing relevant studies, we propose that aggregates, as relatively separate micro-environments embedded in a complex soil matrix, can be viewed as biogeochemical reactors of GHGs. Aggregate reactivity is determined by both aggregate size (which determines the reactor size) and the bulk soil environment including both biotic and abiotic factors (which further influence the reaction conditions). With a systematic, dynamic view of the soil system, implications of aggregate reactors for soil-atmosphere GHG exchange are determined by both an individual reactor's reactivity and dynamics in aggregate size distributions. Emerging evidence supports the contention that aggregate reactors significantly influence soil-atmosphere GHG exchange and may have global implications for carbon and nitrogen cycling. In the context of increasingly frequent and severe disturbances, we advocate more analyses of GHG activities at the aggregate scale. To complement data on aggregate reactors, we suggest developing bottom-up aggregate-based models (ABMs) that apply a trait-based approach and incorporate soil system heterogeneity.

A moisture function of soil heterotrophic respiration that incorporates microscale processes

Soil heterotrophic respiration (HR) is an important source of soil-to-atmosphere CO₂ flux, but its response to changes in soil water content (θ) is poorly understood. Earth system models commonly use empirical moisture functions to describe the HR–θ relationship, introducing significant uncertainty in predicting CO₂ flux from soils. Generalized, mechanistic models that address this uncertainty are thus urgently needed. Here we derive, test, and calibrate a novel moisture function, f_m, that encapsulates primary physicochemical and biological processes controlling soil HR. We validated f_m using simulation results and published experimental data, and established the quantitative relationships between parameters of f_m and measurable soil properties, which enables f_m to predict the HR–θ relationships for different soils across spatial scales. The f_m function predicted comparable HR–θ relationships with laboratory and field measurements, and may reduce the uncertainty in predicting the response of soil organic carbon stocks to climate change compared with the empirical moisture functions currently used in Earth system models.

Empirical moisture functions that describe the relationship between soil heterotrophic respiration and moisture introduce considerable uncertainty in soil CO₂ flux predictions. Here, the authors derive a process-based moisture function by incorporating mechanisms that control soil heterotrophic respiration.