Microbial community dynamics in soil aggregates shape biogeochemical gas fluxes from soil profiles – upscaling an aggregate biophysical model

Microbial communities inhabiting soil aggregates dynamically adjust their activity and composition in response to variations in hydration and other external conditions. These rapid dynamics shape signatures of biogeochemical activity and gas fluxes emitted from soil profiles. Recent mechanistic models of microbial processes in unsaturated aggregate‐like pore networks revealed a highly dynamic interplay between oxic and anoxic microsites jointly shaped by hydration conditions and by aerobic and anaerobic microbial community abundance and self‐organization. The spatial extent of anoxic niches (hotspots) flicker in time (hot moments) and support substantial anaerobic microbial activity even in aerated soil profiles. We employed an individual‐based model for microbial community life in soil aggregate assemblies represented by 3D angular pore networks. Model aggregates of different sizes were subjected to variable water, carbon and oxygen contents that varied with soil depth as boundary conditions. The study integrates microbial activity within aggregates of different sizes and soil depth to obtain estimates of biogeochemical fluxes from the soil profile. The results quantify impacts of dynamic shifts in microbial community composition on CO₂ and N₂O production rates in soil profiles in good agreement with experimental data. Aggregate size distribution and the shape of resource profiles in a soil determine how hydration dynamics shape denitrification and carbon utilization rates. Results from the mechanistic model for microbial activity in aggregates of different sizes were used to derive parameters for analytical representation of soil biogeochemical processes across large scales of practical interest for hydrological and climate models.     

Characterising the biophysical,economic and social impacts of soil carbon sequestration as a greenhouse gas removal technology

To limit warming to well below 2°C, most scenario projections rely on greenhouse gas removal technologies (GGRTs); one such GGRT uses soil carbon sequestration (SCS) in agricultural land. In addition to their role in mitigating climate change, SCS practices play a role in delivering agroecosystem resilience, climate change adaptability and food security. Environmental heterogeneity and differences in agricultural practices challenge the practical implementation of SCS, and our analysis addresses the associated knowledge gap. Previous assessments have focused on global potentials, but there is a need among policymakers to operationalise SCS. Here, we assess a range of practices already proposed to deliver SCS, and distil these into a subset of specific measures. We provide a multidisciplinary summary of the barriers and potential incentives towards practical implementation of these measures. First, we identify specific practices with potential for both a positive impact on SCS at farm level and an uptake rate compatible with global impact. These focus on: (a) optimising crop primary productivity (e.g. nutrient optimisation, pH management, irrigation); (b) reducing soil disturbance and managing soil physical properties (e.g. improved rotations, minimum till); (c) minimising deliberate removal of C or lateral transport via erosion processes (e.g. support measures, bare fallow reduction); (d) addition of C produced outside the system (e.g. organic manure amendments, biochar addition); (e) provision of additional C inputs within the cropping system (e.g. agroforestry, cover cropping). We then consider economic and non‐cost barriers and incentives for land managers implementing these measures, along with the potential externalised impacts of implementation. This offers a framework and reference point for holistic assessment of the impacts of SCS. Finally, we summarise and discuss the ability of extant scientific approaches to quantify the technical potential and externalities of SCS measures, and the barriers and incentives to their implementation in global agricultural systems.     

Soil pH as the chief modifier for regional nitrous oxide emissions: New evidence and implications for global estimates and mitigation

Nitrous oxide (N₂O) is a greenhouse gas that also plays the primary role in stratospheric ozone depletion. The use of nitrogen fertilizers is known as the major reason for atmospheric N₂O increase. Empirical bottom‐up models therefore estimate agricultural N₂O inventories using N loading as the sole predictor, disregarding the regional heterogeneities in soil inherent response to external N loading. Several environmental factors have been found to influence the response in soil N₂O emission to N fertilization, but their interdependence and relative importance have not been addressed properly. Here, we show that soil pH is the chief factor explaining regional disparities in N₂O emission, using a global meta‐analysis of 1,104 field measurements. The emission factor (EF) of N₂O increases significantly (p < .001) with soil pH decrease. The default EF value of 1.0%, according to IPCC (Intergovernmental Panel on Climate Change) for agricultural soils, occurs at soil pH 6.76. Moreover, changes in EF with N fertilization (i.e. ΔEF) is also negatively correlated (p < .001) with soil pH. This indicates that N₂O emission in acidic soils is more sensitive to changing N fertilization than that in alkaline soils. Incorporating our findings into bottom‐up models has significant consequences for regional and global N₂O emission inventories and reconciling them with those from top‐down models. Moreover, our results allow region‐specific development of tailor‐made N₂O mitigation measures in agriculture.     

Methane and carbon dioxide emissions from inland waters in India – implications for large scale greenhouse gas balances

Inland waters were recently recognized to be important sources of methane (CH₄) and carbon dioxide (CO₂) to the atmosphere, and including inland water emissions in large scale greenhouse gas (GHG) budgets may potentially offset the estimated carbon sink in many areas. However, the lack of GHG flux measurements and well‐defined inland water areas for extrapolation, make the magnitude of the potential offset unclear. This study presents coordinated flux measurements of CH₄ and CO₂ in multiple lakes, ponds, rivers, open wells, reservoirs, springs, and canals in India. All these inland water types, representative of common aquatic ecosystems in India, emitted substantial amounts of CH₄ and a major fraction also emitted CO₂. The total CH₄ flux (including ebullition and diffusion) from all the 45 systems ranged from 0.01 to 52.1 mmol m^?2 d^?1, with a mean of 7.8 ± 12.7 (mean ± 1 SD) mmol m^?2 d^?1. The mean surface water CH₄ concentration was 3.8 ± 14.5 μm (range 0.03–92.1 μm ). The CO₂ fluxes ranged from ?28.2 to 262.4 mmol m^?2 d^?1 and the mean flux was 51.9 ± 71.1 mmol m^?2 d^?1. The mean partial pressure of CO₂ was 2927 ± 3269 μatm (range: 400–11 467 μatm). Conservative extrapolation to whole India, considering the specific area of the different water types studied, yielded average emissions of 2.1 Tg CH₄ yr^?1 and 22.0 Tg CO₂ yr^?1 from India's inland waters. When expressed as CO₂ equivalents, this amounts to 75 Tg CO₂ equivalents yr^?1 (53–98 Tg CO₂ equivalents yr^?1; ± 1 SD)_, with CH₄ contributing 71%. Hence, average inland water GHG emissions, which were not previously considered, correspond to 42% (30–55%) of the estimated land carbon sink of India. Thereby this study illustrates the importance of considering inland water GHG exchange in large scale assessments.     

Soil–atmosphere exchange of CH4, CO,N2O and NOx and the effects of land‐use change in the semiarid Mallee system in Southeastern Australia

The semiarid and arid zones cover a quarter of the global land area and support one‐fifth of the world's human population. A significant fraction of the global soil–atmosphere exchange for climatically active gases occurs in semiarid and arid zones yet little is known about these exchanges. A study was made of the soil–atmosphere exchange of CH₄, CO, N₂O and NO_x in the semiarid Mallee system, in north‐western Victoria, Australia, at two sites: one pristine mallee and the other cleared for approximately 65 years for farming (currently wheat). The mean (± standard error) rates of CH₄ exchange were uptakes of ?3.0 ± 0.5 ng(C) m^?2 s^?1 for the Mallee and ?6.0 ± 0.3 ng(C) m^?2 s^?1 for the Wheat. Converting mallee forest to wheat crop increases CH₄ uptake significantly. CH₄ emissions were observed in the Mallee in summer and were hypothesized to arise from termite activity. We find no evidence that in situ growing wheat plants emit CH₄, contrary to a recent report. The average CO emissions of 10.1 ± 1.8 ng(C) m^?2 s^?1 in the Mallee and 12.6 ± 2.0 ng(C) m^?2 s^?1 in the Wheat. The average N₂O emissions were 0.5 ± 0.1 ng(N) m^?2 s^?1 from the pristine Mallee and 1.4 ± 0.3 ng(N) m^?2 s^?1 from the Wheat. The experimental results show that the processes controlling these exchanges are different to those in temperate systems and are poorly understood.     

The production‐ecological sustainability of cassava,sugarcane and sweet sorghum cultivation for bioethanol in Mozambique

We present an approach for providing quantitative insight into the production‐ecological sustainability of biofuel feedstock production systems. The approach is based on a simple crop‐soil model and was used for assessing feedstock from current and improved production systems of cassava for bioethanol. Assessments were performed for a study area in Mozambique, a country considered promising for biomass production. Our focus is on the potential role of smallholders in the production of feedstock for biofuels. We take cassava as the crop for this purpose and compare it with feedstock production on plantations using sugarcane, sweet sorghum and cassava as benchmarks. Production‐ecological sustainability was defined by seven indicators related to resource‐use efficiency, soil quality, net energy production and greenhouse gas (GHG) emissions. Results indicate that of the assessed systems, sugarcane performed better than cassava, although it requires substantial water for irrigation. Targeted use of nutrient inputs improved sustainability of smallholder cassava. Cassava production systems on more fertile soils were more sustainable than those on less fertile soils; the latter required more external inputs for achieving the same output, affecting most indicators negatively and reducing the feasibility for smallholders. Cassava and sweet sorghum performed similarly. Cassava production requires much more labour per hectare than production of sugarcane or sweet sorghum. Production of bioethanol feedstock on cultivated lands was more sustainable and had potential for carbon sequestration, avoiding GHG emissions from clearing natural vegetation if new land is opened.     

Comparing the consequences of natural selection,adaptive phenotypic plasticity,and matching habitat choice for phenotype–environment matching,population genetic structure,and reproductive isolation in meta‐populations

Organisms commonly experience significant spatiotemporal variation in their environments. In response to such heterogeneity, different mechanisms may act that enhance ecological performance locally. However, depending on the nature of the mechanism involved, the consequences for populations may differ greatly. Building on a previous model that investigated the conditions under which different adaptive mechanisms (co)evolve, this study compares the ecological and evolutionary population consequences of three very different responses to environmental heterogeneity: matching habitat choice (directed gene flow), adaptive plasticity (associated with random gene flow), and divergent natural selection. Using individual‐based simulations, we show that matching habitat choice can have a greater adaptive potential than plasticity or natural selection: it allows for local adaptation while protecting genetic polymorphism despite global mating or strong environmental changes. Our simulations further reveal that increasing environmental fluctuations and unpredictability generally favor the emergence of specialist genotypes but that matching habitat choice is better at preventing local maladaptation by individuals. This confirms that matching habitat choice can speed up the genetic divergence among populations, cause indirect assortative mating via spatial clustering, and hence even facilitate sympatric speciation. This study highlights the potential importance of directed dispersal in local adaptation and speciation, stresses the difficulty of deriving its operation from nonexperimental observational data alone, and helps define a set of ecological conditions which should favor its emergence and subsequent detection in nature.     

Sensitivity of Soil Restoration Indicators during Páramo Succession in the High Tropical Andes: Chronosequence and Permanent Plot Approaches

The aim of this study was to analyze indicators of soil restoration during secondary succession in a heterogeneous valley in the high tropical Andes. A combination of chronosequence and permanent plot methods was used to detect changes in this heterogeneous matrix. Thirty-six plots with different fallow times (1–9 years) and four noncultivated plots with natural vegetation (páramo) were sampled twice in a 3-year interval (1996 and 1999). The following soil properties were determined: total C and N, pH, exchangeable bases, cation exchange capacity, and microbial biomass N (MB-N). Using the chronosequence approach, successional increases in soil pH and Mg were detected, pointing to these variables as indicators of soil restoration during the fallow period. Comparing the noncultivated páramo with the fallow plots, a significant decrease in MB-N was found, suggesting that this is a sensitive agricultural disturbance indicator. The permanent plot analysis failed to detect successional trends in any of the study variables, probably as a result of a lack of sensitivity of the indicators used within the 3-year interval. Nevertheless, a strong acidification was detected by the permanent plot method when fallow plots were cultivated. We conclude that the size of important soil components such as total soil organic matter or microbial biomass is not a sensitive soil restoration indicator in these heterogeneous mountain systems but that other integrative variables such as pH could be more sensitive to successional changes in key soil processes (e.g., nitrification or nutrient losses).     

Conversion of coastal wetlands,riparian wetlands,and peatlands increases greenhouse gas emissions: A global meta‐analysis

Land‐use/land‐cover change (LULCC) often results in degradation of natural wetlands and affects the dynamics of greenhouse gases (GHGs). However, the magnitude of changes in GHG emissions from wetlands undergoing various LULCC types remains unclear. We conducted a global meta‐analysis with a database of 209 sites to examine the effects of LULCC types of constructed wetlands (CWs), croplands (CLs), aquaculture ponds (APs), drained wetlands (DWs), and pastures (PASs) on the variability in CO₂, CH₄, and N₂O emissions from the natural coastal wetlands, riparian wetlands, and peatlands. Our results showed that the natural wetlands were net sinks of atmospheric CO₂ and net sources of CH₄ and N₂O, exhibiting the capacity to mitigate greenhouse effects due to negative comprehensive global warming potentials (GWPs; ?0.9 to ?8.7 t CO₂‐eq ha^?1 year^?1). Relative to the natural wetlands, all LULCC types (except CWs from coastal wetlands) decreased the net CO₂ uptake by 69.7%?456.6%, due to a higher increase in ecosystem respiration relative to slight changes in gross primary production. The CWs and APs significantly increased the CH₄ emissions compared to those of the coastal wetlands. All LULCC types associated with the riparian wetlands significantly decreased the CH₄ emissions. When the peatlands were converted to the PASs, the CH₄ emissions significantly increased. The CLs, as well as DWs from peatlands, significantly increased the N₂O emissions in the natural wetlands. As a result, all LULCC types (except PASs from riparian wetlands) led to remarkably higher GWPs by 65.4%?2,948.8%, compared to those of the natural wetlands. The variability in GHG fluxes with LULCC was mainly sensitive to changes in soil water content, water table, salinity, soil nitrogen content, soil pH, and bulk density. This study highlights the significant role of LULCC in increasing comprehensive GHG emissions from global natural wetlands, and our results are useful for improving future models and manipulative experiments.     

Testing for Genetic Association in the Presence of Linkage and Gene–Covariate Interactions

In order to study family‐based association in the presence of linkage, we extend a generalized linear mixed model proposed for genetic linkage analysis (Lebrec and van Houwelingen (2007), Human Heredity 64 , 5–15) by adding a genotypic effect to the mean. The corresponding score test is a weighted family‐based association tests statistic, where the weight depends on the linkage effect and on other genetic and shared environmental effects. For testing of genetic association in the presence of gene–covariate interaction, we propose a linear regression method where the family‐specific score statistic is regressed on family‐specific covariates. Both statistics are straightforward to compute. Simulation results show that adjusting the weight for the within‐family variance structure may be a powerful approach in the presence of environmental effects. The test statistic for genetic association in the presence of gene–covariate interaction improved the power for detecting association. For illustration, we analyze the rheumatoid arthritis data from GAW15. Adjusting for smoking and anti‐cyclic citrullinated peptide increased the significance of the association with the DR locus.     

The effects of species pool,dispersal and competition on the diversity–productivity relationship

Aim The diversity–productivity relationship is a controversial issue in ecology. Diversity is sometimes seen to increase with productivity but a unimodal relationship has often been reported. Competitive exclusion was cited initially to account for the decrease of diversity at high productivity. Subsequently, the roles of evolutionary history (species pool size) and dispersal rate have been acknowledged. We explore how the effects of species pool, dispersal and competition combine to produce different diversity–productivity relationships. Methods We use a series of simulations with a spatially explicit, individual‐based model. Following empirical expectations, we used four scenarios to characterize species pool size along the productivity gradient (uniformly low and high, linear increase and unimodal). Similarly, the dispersal rate varied along the productivity gradient (uniformly low and high, and unimodal). We considered both neutral communities and communities with competitive exclusion. Results and main conclusions Our model predicts that competitive interactions will result in unimodal diversity–productivity relationships. The model often predicts unimodal patterns in neutral communities as well, although the decline in richness at high productivity is less than in competing communities. A positive diversity–productivity relationship is simulated for neutral communities when the species pool size increases with productivity and the dispersal rate is high. This scenario is probably more widespread in nature than the others since positive diversity–productivity relationships have been observed more frequently than previously expected, especially in the tropics and for woody species. Our simulated effects of species pool, dispersal and competition on diversity patterns can be linked to empirical observations to uncover mechanisms behind the diversity–productivity relationship.     

Scale and uncertainty in modeled soil organic carbon stock changes for US croplands using a process‐based model

Process‐based model analyses are often used to estimate changes in soil organic carbon (SOC), particularly at regional to continental scales. However, uncertainties are rarely evaluated, and so it is difficult to determine how much confidence can be placed in the results. Our objective was to quantify uncertainties across multiple scales in a process‐based model analysis, and provide 95% confidence intervals for the estimates. Specifically, we used the Century ecosystem model to estimate changes in SOC stocks for US croplands during the 1990s, addressing uncertainties in model inputs, structure and scaling of results from point locations to regions and the entire country. Overall, SOC stocks increased in US croplands by 14.6 Tg C yr^?1 from 1990 to 1995 and 17.5 Tg C yr^?1 during 1995 to 2000, and uncertainties were ±22% and ±16% for the two time periods, respectively. Uncertainties were inversely related to spatial scale, with median uncertainties at the regional scale estimated at ±118% and ±114% during the early and latter part of 1990s, and even higher at the site scale with estimates at ±739% and ±674% for the time periods, respectively. This relationship appeared to be driven by the amount of the SOC stock change; changes in stocks that exceeded 200 Gg C yr^?1 represented a threshold where uncertainties were always lower than ±100%. Consequently, the amount of uncertainty in estimates derived from process‐based models will partly depend on the level of SOC accumulation or loss. In general, the majority of uncertainty was associated with model structure in this application, and so attaining higher levels of precision in the estimates will largely depend on improving the model algorithms and parameterization, as well as increasing the number of measurement sites used to evaluate the structural uncertainty.     

The coevolution of juvenile play–fighting and adult competition

Although play–fighting is widespread among juvenile mammals, its adaptive significance remains unclear. It has been proposed that play is beneficial for developing skills to improve success in adult contests (motor‐training hypothesis), but the links between juvenile play–fighting and adult aggression are complex and not well understood. In this theoretical study, we investigate the coevolution between juvenile play–fighting and adult aggression using evolutionary computer simulations. We consider a simple life history with two sequential stages: a juvenile phase in which individuals play–fight with other juveniles to develop their fighting skills; and an adult phase in which individuals engage in potentially aggressive contests over access to resources and ultimately mating opportunities, leading to reproductive success. The simulations track genetic evolution in key traits affecting adult contests, such as the level of aggression, as well as juvenile investment in play–fighting, capturing the coevolutionary feedbacks between juvenile and adult decisions. We find that coevolution leads to one of two outcomes: a high‐play, high‐aggression situation with highly aggressive adult contests preceded by a prolonged period of juvenile play–fighting to improve fighting ability, or a low‐play, low‐aggression situation in which adult contests are resolved without fighting and there is minimal investment in play–fighting before individuals mature. Which of these outcomes is favoured depends on the mortality costs and on the type of societal structure: societies with strong reproductive skew, favouring monopolization of resources, show high levels of adult aggression and high investment in juvenile play–fighting, whereas societies with low reproductive skew have both low adult aggression and low levels of play–fighting. A review of empirical evidence, particularly in the primate genus Macaca, highlights some limitations of our model and suggests that other, complementary functional explanations are needed to account for the full range of competitive and cooperative forms of play–fighting. Our study illustrates the power of evolutionary simulations to shed light on the long‐standing puzzle of animal play.     

A trans‐national monarch butterfly population model and implications for regional conservation priorities

1. The monarch has undergone considerable population declines over the past decade, and the governments of Mexico, Canada, and the United States have agreed to work together to conserve the species. 2. Given limited resources, understanding where to focus conservation action is key for widespread species like monarchs. To support planning for continental‐scale monarch habitat restoration, we address the question of where restoration efforts are likely to have the largest impacts on monarch butterfly (Danaus plexippus Linn.) population growth rates. 3. We present a spatially explicit demographic model simulating the multi‐generational annual cycle of the eastern monarch population, and use the model to examine management scenarios, some of which focus on particular regions of North America. 4. Improving the monarch habitat in the north central or southern parts of the monarch range yields a slightly greater increase in the population growth rate than restoration in other regions. However, combining restoration efforts across multiple regions yields population growth rates above 1 with smaller simulated improvements in habitat per region than single‐region strategies. 5. S ynthesis and applications: These findings suggest that conservation investment in projects across the full monarch range will be more effective than focusing on one or a few regions, and will require international cooperation across many land use categories.     

Soil carbon and nitrogen cycling and storage throughout the soil profile in a sweetgum plantation after 11 years of CO2‐enrichment

Increased partitioning of carbon (C) to fine roots under elevated [CO₂], especially deep in the soil profile, could alter soil C and nitrogen (N) cycling in forests. After more than 11 years of free‐air CO₂ enrichment in a Liquidambar styraciflua L. (sweetgum) plantation in Oak Ridge, TN, USA, greater inputs of fine roots resulted in the incorporation of new C (i.e., C with a depleted δ¹³C) into root‐derived particulate organic matter (POM) pools to 90‐cm depth. Even though production in the sweetgum stand was limited by soil N availability, soil C and N contents were greater throughout the soil profile under elevated [CO₂] at the conclusion of the experiment. Greater C inputs from fine‐root detritus under elevated [CO₂] did not result in increased net N immobilization or C mineralization rates in long‐term laboratory incubations, possibly because microbial biomass was lower in the CO₂‐enriched plots. Furthermore, the δ¹³CO₂ of the C mineralized from the incubated soil closely tracked the δ¹³C of the labile POM pool in the elevated [CO₂] treatment, especially in shallower soil, and did not indicate significant priming of the decomposition of pre‐experiment soil organic matter (SOM). Although potential C mineralization rates were positively and linearly related to total SOM C content in the top 30 cm of soil, this relationship did not hold in deeper soil. Taken together with an increased mean residence time of C in deeper soil pools, these findings indicate that C inputs from relatively deep roots under elevated [CO₂] may increase the potential for long‐term soil C storage. However, C in deeper soil is likely to take many years to accrue to a significant fraction of total soil C given relatively smaller root inputs at depth. Expanded representation of biogeochemical cycling throughout the soil profile may improve model projections of future forest responses to rising atmospheric [CO₂].     

Long‐term no‐till and stover retention each decrease the global warming potential of irrigated continuous corn

Over the last 50 years, the most increase in cultivated land area globally has been due to a doubling of irrigated land. Long‐term agronomic management impacts on soil organic carbon (SOC) stocks, soil greenhouse gas (GHG) emissions, and global warming potential (GWP) in irrigated systems, however, remain relatively unknown. Here, residue and tillage management effects were quantified by measuring soil nitrous oxide (N₂O) and methane (CH₄) fluxes and SOC changes (ΔSOC) at a long‐term, irrigated continuous corn (Zea mays L.) system in eastern Nebraska, United States. Management treatments began in 2002, and measured treatments included no or high stover removal (0 or 6.8 Mg DM ha^?1 yr^?1, respectively) under no‐till (NT) or conventional disk tillage (CT) with full irrigation (n = 4). Soil N₂O and CH₄ fluxes were measured for five crop‐years (2011–2015), and ΔSOC was determined on an equivalent mass basis to ~30 cm soil depth. Both area‐ and yield‐scaled soil N₂O emissions were greater with stover retention compared to removal and for CT compared to NT, with no interaction between stover and tillage practices. Methane comprised <1% of total emissions, with NT being CH₄ neutral and CT a CH₄ source. Surface SOC decreased with stover removal and with CT after 14 years of management. When ΔSOC, soil GHG emissions, and agronomic energy usage were used to calculate system GWP, all management systems were net GHG sources. Conservation practices (NT, stover retention) each decreased system GWP compared to conventional practices (CT, stover removal), but pairing conservation practices conferred no additional mitigation benefit. Although cropping system, management equipment/timing/history, soil type, location, weather, and the depth to which ΔSOC is measured affect the GWP outcomes of irrigated systems at large, this long‐term irrigated study provides valuable empirical evidence of how management decisions can impact soil GHG emissions and surface SOC stocks.