Drought resistance and resilience of rhizosphere communities in forest soils from the cellular to ecosystem scale - insights from 13C pulse labeling.

The link between above- and belowground communities is a key uncertainty in drought and rewetting effects on forest carbon (C) cycle. In young beech model ecosystems and mature naturally dry pine forest exposed to 15-yr-long irrigation, we performed 13C pulse labeling experiments, one during drought and one 2 wk after rewetting, tracing tree assimilates into rhizosphere communities. The 13C pulses applied in tree crowns reached soil microbial communities of the young and mature forests one and 4 d later, respectively. Drought decreased the transfer of labeled assimilates relative to the irrigation treatment. The 13C label in phospholipid fatty acids (PLFAs) indicated greater drought reduction of assimilate incorporation by fungi (-85%) than by gram-positive (-43%) and gram-negative bacteria (-58%). 13C label incorporation was more strongly reduced for PLFAs (cell membrane) than for microbial cytoplasm extracted by chloroform. This suggests that fresh rhizodeposits are predominantly used for osmoregulation or storage under drought, at the expense of new cell formation. Two weeks after rewetting, 13C enrichment in PLFAs was greater in previously dry than in continuously moist soils. Drought and rewetting effects were greater in beech systems than in pine forest. Belowground C allocation and rhizosphere communities are highly resilient to drought.

Formation of mineral-associated organic matter in temperate soils is primarily controlled by mineral type and modified by land use and management intensity.

Formation of mineral-associated organic matter (MAOM) supports the accumulation and stabilization of carbon (C) in soil, and thus, is a key factor in the global C cycle. Little is known about the interplay of mineral type, land use and management intensity in MAOM formation, especially on subdecadal time scales. We exposed mineral containers with goethite or illite, the most abundant iron oxide and phyllosilicate clay in temperate soils, for 5 years in topsoils of 150 forest and 150 grassland sites in three regions across Germany. Results show that irrespective of land use and management intensity, more C accumulated on goethite than illite (on average 0.23 ± 0.10 and 0.06 ± 0.03 mg m-2 mineral surface respectively). Carbon accumulation across regions was consistently higher in coniferous forests than in deciduous forests and grasslands. Structural equation models further showed that thinning and harvesting reduced MAOM formation in forests. Formation of MAOM in grasslands was not affected by grazing. Fertilization had opposite effects on MAOM formation, with the positive effect being mediated by enhanced plant productivity and the negative effect by reduced plant species richness. This highlights the caveat of applying fertilizers as a strategy to increase soil C stocks in temperate grasslands. Overall, we demonstrate that the rate and amount of MAOM formation in soil is primarily driven by mineral type, and can be modulated by land use and management intensity even on subdecadal time scales. Our results suggest that temperate soils dominated by oxides have a higher capacity to accumulate and store C than those dominated by phyllosilicate clays, even under circumneutral pH conditions. Therefore, adopting land use and management practices that increase C inputs into oxide-rich soils that are under their capacity to store C may offer great potential to enhance near-term soil C sequestration.

Rhizosphere activity in an old-growth forest reacts rapidly to changes in soil moisture and shapes whole-tree carbon allocation.

Drought alters carbon (C) allocation within trees, thereby impairing tree growth. Recovery of root and leaf functioning and prioritized C supply to sink tissues after drought may compensate for drought-induced reduction of assimilation and growth. It remains unclear if C allocation to sink tissues during and following drought is controlled by altered sink metabolic activities or by the availability of new assimilates. Understanding such mechanisms is required to predict forests' resilience to a changing climate. We investigated the impact of drought and drought release on C allocation in a 100-y-old Scots pine forest. We applied 13CO2 pulse labeling to naturally dry control and long-term irrigated trees and tracked the fate of the label in above- and belowground C pools and fluxes. Allocation of new assimilates belowground was ca. 53% lower under nonirrigated conditions. A short rainfall event, which led to a temporary increase in the soil water content (SWC) in the topsoil, strongly increased the amounts of C transported belowground in the nonirrigated plots to values comparable to those in the irrigated plots. This switch in allocation patterns was congruent with a tipping point at around 15% SWC in the response of the respiratory activity of soil microbes. These results indicate that the metabolic sink activity in the rhizosphere and its modulation by soil moisture can drive C allocation within adult trees and ecosystems. Even a subtle increase in soil moisture can lead to a rapid recovery of belowground functions that in turn affects the direction of C transport in trees.

Long-term manipulation of mean climatic conditions alters drought effects on C- and N-cycling in an arable soil.

Climate is changing and predicted future scenarios include both changes in long-term mean climatic conditions and intensification of extreme events such as drought. Drought can have a major impact on soil functional processes; soil microorganisms, key to these processes, depend on water and temperature dynamics. Consequently, feedback mechanisms regarding microbially mediated carbon and nitrogen cycling in soils may be affected. There are indications that microbial exposure to increasingly unfavorable environmental conditions influences their stress responses. Here, the long-term field experiment Hohenheim Climate Change (HoCC) provided a research platform to explore how microbial exposure to long-term reduced water availability and soil warming modifies microbially driven soil processes, especially gas fluxes from soil, both during drought and after rewetting. The HoCC experiment is an agroecosystem in which the soil microbiome has been exposed to reduced annual mean precipitation and elevated temperature since 2008. Treatment levels were chosen based on a realistic future climate scenario. In June 2019, we exposed this system to a drought period of four weeks. We found that even after 11 years, warming remained a driver of CO2 and N2 O fluxes across the different soil moisture conditions in our drought experiment. Importantly, however, microbial exposure to long-term reduced water availability limited the stimulatory effect of warming on gas fluxes during drought and after rewetting. Our results were neither related to a legacy effect within overall microbial biomass carbon levels nor a shift towards enhanced fungal abundance. We found no indications that extracellular enzyme activities or microbial substrate availability explained the gas flux dynamics observed in our drought experiment. Our study indicates that soil warming promotes gaseous C and N loss even under extreme drought conditions. We suspect, however, that a shift in microbial function following long-term water limitation can hamper the enhancing effect of warming on soil gas fluxes.

Temperature response of soil respiration largely unaltered with experimental warming.

The respiratory release of carbon dioxide (CO2) from soil is a major yet poorly understood flux in the global carbon cycle. Climatic warming is hypothesized to increase rates of soil respiration, potentially fueling further increases in global temperatures. However, despite considerable scientific attention in recent decades, the overall response of soil respiration to anticipated climatic warming remains unclear. We synthesize the largest global dataset to date of soil respiration, moisture, and temperature measurements, totaling >3,800 observations representing 27 temperature manipulation studies, spanning nine biomes and over 2 decades of warming. Our analysis reveals no significant differences in the temperature sensitivity of soil respiration between control and warmed plots in all biomes, with the exception of deserts and boreal forests. Thus, our data provide limited evidence of acclimation of soil respiration to experimental warming in several major biome types, contrary to the results from multiple single-site studies. Moreover, across all nondesert biomes, respiration rates with and without experimental warming follow a Gaussian response, increasing with soil temperature up to a threshold of ∼25 °C, above which respiration rates decrease with further increases in temperature. This consistent decrease in temperature sensitivity at higher temperatures demonstrates that rising global temperatures may result in regionally variable responses in soil respiration, with colder climates being considerably more responsive to increased ambient temperatures compared with warmer regions. Our analysis adds a unique cross-biome perspective on the temperature response of soil respiration, information critical to improving our mechanistic understanding of how soil carbon dynamics change with climatic warming.

Offsetting global warming-induced elevated greenhouse gas emissions from an arable soil by biochar application.

Global warming will likely enhance greenhouse gas (GHG) emissions from soils. Due to its slow decomposability, biochar is widely recognized as effective in long-term soil carbon (C) sequestration and in mitigation of soil GHG emissions. In a long-term soil warming experiment (+2.5 °C, since July 2008) we studied the effect of applying high-temperature Miscanthus biochar (0, 30 t/ha, since August 2013) on GHG emissions and their global warming potential (GWP) during 2 years in a temperate agroecosystem. Crop growth, physical and chemical soil properties, temperature sensitivity of soil respiration (Rs ), and metabolic quotient (qCO2 ) were investigated to yield further information about single effects of soil warming and biochar as well as on their interactions. Soil warming increased total CO2 emissions by 28% over 2 years. The effect of warming on soil respiration did not level off as has often been observed in less intensively managed ecosystems. However, the temperature sensitivity of soil respiration was not affected by warming. Overall, biochar had no effect on most of the measured parameters, suggesting its high degradation stability and its low influence on microbial C cycling even under elevated soil temperatures. In contrast, biochar × warming interactions led to higher total N2 O emissions, possibly due to accelerated N-cycling at elevated soil temperature and to biochar-induced changes in soil properties and environmental conditions. Methane uptake was not affected by soil warming or biochar. The incorporation of biochar-C into soil was estimated to offset warming-induced elevated GHG emissions for 25 years. Our results highlight the suitability of biochar for C sequestration in cultivated temperate agricultural soil under a future elevated temperature. However, the increased N2 O emissions under warming limit the GHG mitigation potential of biochar.

13C assimilation as well as functional gene abundance and expression elucidate the biodegradation of glyphosate in a field experiment.

Glyphosate (N-phosphonomethylglycine; GLP) and its main metabolite AMPA (aminomethylphosphonic acid), are frequently detected in relatively high concentrations in European agricultural topsoils. Glyphosate has a high sorption affinity, yet it can be detected occasionally in groundwater. We hypothesized that shrinkage cracks occurring after dry periods could facilitate GLP transport to greater depths where subsoil conditions slow further microbial degradation. To test this hypothesis, we simulated a heavy rainfall event (HRE) on a clay-rich arable soil. We applied 2.1 kg ha-1 of 100% 13C3, 15N-labeled GLP one day before the simulated rainfall event. Microbial degradation of translocated GLP over a 21-day period was assessed by quantifying 13C incorporation into phospholipid fatty acids. Microbial degradation potential and activity were determined by quantifying the abundance and expression of functional genes involved in the two known degradation pathways of GLP; to AMPA (goxA) or sarcosine (sarc). We confirmed that goxA transcripts were elevated in the range of 4.23 x 105 copy numbers g-1 soil only one day after application. The increase in AMPA associated with a rise in goxA transcripts and goxA-harboring microorganisms indicated that the degradation pathway to AMPA dominated. Based on 13C-enrichment 3 h after the HRE, fungi appeared to initiate glyphosate degradation. At later time points, Gram+-bacteria proved to be the main degraders due to their higher 13C-incorporation. Once GLP reached the subsoil, degradation continued but more slowly. By comparing GLP distribution and its microbial degradation in macropores and in the bulk soil, we demonstrated different time- and depth-dependent GLP degradation dynamics in macropores. This indicates the need for field studies in which soil properties relevant to GLP degradation are related to limiting environmental conditions, providing a realistic assessment of GLP fate in soils.

Heavy rainfall following a summer drought stimulates soil redox dynamics and facilitates rapid and deep translocation of glyphosate in floodplain soils.

We present field data on the effects of heavy rainfall after drought on the mobility of glyphosate and redox conditions in a clayey floodplain soil. By applying glyphosate together with deuterated water as conservative tracer in combination with time resolved in situ redox potential measurements, the spatial and temporal patterns of water infiltration and pesticide transport as well as the concomitant changes of the redox conditions were revealed. Our findings demonstrate that shrinkage cracks in dry soils can serve as effective transport paths for atmospheric oxygen, water and glyphosate. The rain intensity of a typical summer storm event (approx. 25 mm within one hour) was sufficient to translocate deuterated water and glyphosate to the subsoil (50 cm) within 2 hours. Soil wetting induced partial closure of the shrinkage cracks and stimulated microbial activity resulting in pronounced dynamics of in situ soil redox conditions. Redox potentials in 40 to 50 cm depth dropped permanently to strongly reducing conditions within hours to days but fluctuated between reducing and oxidizing conditions in 10 to 30 cm depth. Our findings highlight the close link between the presence of macropores (shrinkage cracks), heavy rainfall after drought, redox dynamics and pesticide translocation to the subsoil and thus call for further studies addressing the effects of dynamic redox conditions as a limiting factor for glyphosate degradation.

Small-scale diversity and succession of fungi in the detritusphere of rye residues.

Transport of litter carbon in the detritusphere might determine fungal abundance and diversity at the small scale. Rye residues were applied to the surface of soil cores with two different water contents and incubated at 10 degrees C for 2 and 12 weeks. Fungal community structure was analysed by constructing clone libraries of 18S rDNA and subsequent sequencing. Litter addition induced fungal succession in the adjacent soil and decreased detectable fungal diversity mainly due to the huge supply of substrates. Ergosterol content and N-acetyl-glucosaminidase activity indicated fungal growth after 2 weeks. Simultaneously, the structure of the fungal community changed, with Mortierellaceae proliferating during the initial phase of litter decomposition. Ergosterol measurements were unable to detect this early fungal growth because Mortierellaceae do not produce ergosterol. In the late phase during decomposition of polymeric substrates, like cellulose and chitin, the fungal community was dominated by Trichocladium asperum. Water content influenced community composition only during the first 2 weeks due to its influence on transport processes in the detritusphere and on competition between fungal species. Our results underline the importance of species identification in understanding decomposition processes in soil.

Plant litter enhances degradation of the herbicide MCPA and increases formation of biogenic non-extractable residues in soil.

Amendment of soils with plant residues is common practice for improving soil quality. In addition to stimulated microbial activity, the supply of fresh soluble organic (C) from litter may accelerate the microbial degradation of chemicals in soils. Therefore, the aim of this study was to test whether the maize litter enhances degradation of 4-chloro-2-methylphenoxyacetic acid (MCPA) and increases formation of non-toxic biogenic non-extractable residues (bioNERs). Soil was amended with 13C6-MCPA and incubated with or without litter addition on the top. Three soil layers were sampled with increasing distance from the top: 0-2 mm, 2-5 mm and 5-20 mm; and the mass balance of 13C6-MCPA transformation determined. Maize litter promoted microbial activity, mineralization of 13C6-MCPA and bioNER formation in the upper two layers (0-2 and 2-5 mm). The mineralization of 13C6-MCPA in soil with litter increased to 27% compared to only 6% in the control. Accordingly, maize addition reduced the amount of extractable residual MCPA in soil from 77% (control) to 35% of initially applied 13C6-MCPA. While non-extractable residues (NERs) were <6% in control soil, litter addition raised NERs to 21%. Thereby, bioNERs comprised 14% of 13C6-MCPA equivalents. We found characteristic differences of bioNER formation with distance to litter. While total NERs in soil at a distance of 2-5 mm were mostly identified as 13C-bioNERs (97%), only 45-46% of total NERs were assigned to bioNERs in the 0-2 and 5-20 mm layers. Phospholipid fatty acid analysis indicated that fungi and Gram-negative bacteria were mainly involved in MCPA degradation. Maize-C particularly stimulated fungal activity in the adjacent soil, which presumably facilitated non-biogenic NER formation. The plant litter accelerated formation of both non-toxic bioNERs and non-biogenic NERs. More studies on the structural composition of non-biogenic NERs with toxicity potential are needed for future recommendations on litter addition in agriculture.

Biodegradation of Pesticides at the Limit: Kinetics and Microbial Substrate Use at Low Concentrations.

The objective of our study was to test whether limited microbial degradation at low pesticide concentrations could explain the discrepancy between overall degradability demonstrated in laboratory tests and their actual persistence in the environment. Studies on pesticide degradation are often performed using unrealistically high application rates seldom found in natural environments. Nevertheless, biodegradation rates determined for higher pesticide doses cannot necessarily be extrapolated to lower concentrations. In this context, we wanted to (i) compare the kinetics of pesticide degradation at different concentrations in arable land and (ii) clarify whether there is a concentration threshold below which the expression of the functional genes involved in the degradation pathway is inhibited without further pesticide degradation taking place. We set up an incubation experiment for four weeks using 14C-ring labeled 2-methyl-4-chlorophenoxyacetic acid (MCPA) as a model compound in concentrations from 30 to 20,000 µg kg-1 soil. To quantify the abundance of putative microorganisms involved in MCPA degradation and their degradation activity, tfdA gene copy numbers (DNA) and transcripts (mRNA) were determined by quantitative real-time PCR. Mineralization dynamics of MCPA derived-C were analyzed by monitoring 14CO2 production and 14C assimilation by soil microorganisms. We identified two different concentration thresholds for growth and activity with respect to MCPA degradation using tfdA gene and mRNA transcript abundance as growth and activity indices, respectively. The tfdA gene expression started to increase between 1,000 and 5,000 µg MCPA kg-1 dry soil, but an actual increase in tfdA sequences could only be determined at a concentration of 20,000 µg. Accordingly, we observed a clear shift from catabolic to anabolic utilization of MCPA-derived C in the concentration range of 1,000 to 5,000 µg kg-1. Concentrations ≥1,000 µg kg-1 were mainly associated with delayed mineralization, while concentrations ≤1,000 µg kg-1 showed rapid absolute dissipation. The persistence of pesticides at low concentrations cannot, therefore, be explained by the absence of functional gene expression. Nevertheless, significant differences in the degradation kinetics of MCPA between low and high pesticide concentrations illustrate the need for studies investigating pesticide degradation at environmentally relevant concentrations.

Rhizosphere bacteria affected by transgenic potatoes with antibacterial activities compared with the effects of soil, wild-type potatoes, vegetation stage and pathogen exposure.

A greenhouse experiment was performed to analyze a potential effect of genetically modified potatoes expressing antibacterial compounds (attacin/cecropin, T4 lysozyme) and their nearly isogenic, nontransformed parental wild types on rhizosphere bacterial communities. To compare plant transformation-related variations with commonly accepted impacts caused by altered environmental conditions, potatoes were cultivated under different environmental conditions, for example using contrasting soil types. Further, plants were challenged with the blackleg pathogen Erwinia carotovora ssp. atroseptica. Rhizosphere soil samples were obtained at the stem elongation and early flowering stages. The activities of various extracellular rhizosphere enzymes involved in the C-, P- and N-nutrient cycles were determined as the rates of fluorescence of enzymatically hydrolyzed substrates containing the highly fluorescent compounds 4-methylumbelliferone or 7-amino-4-methyl coumarin. The structural diversity of the bacterial communities was assessed by 16S rRNA-based terminal restriction fragment length polymorphism analysis, and 16S rRNA gene clone libraries were established for the flowering conventional and T4 lysozyme-expressing Desirée lines grown on the chernozem soil, each line treated with and without E. carotovora ssp. atroseptica. Both genetic transformation events induced a differentiation in the activity rates and structures of associated bacterial communities. In general, T4 lysozyme had a stronger effect than attacin/cecropin. In comparison with the other factors, the impact of the genetic modification was only transient and minor, or comparable to the dominant variations caused by soil type, plant genotype, vegetation stage and pathogen exposure.

Evidence for the importance of litter as a co-substrate for MCPA dissipation in an agricultural soil.

Environmental controls of 2-methyl-4-chlorophenoxyacetic acid (MCPA) degradation are poorly understood. We investigated whether microbial MCPA degraders are stimulated by (maize) litter and whether this process depends on concentrations of MCPA and litter. In a microcosm experiment, different amounts of litter (0, 10 and 20 g kg(-1)) were added to soils exposed to three levels of the herbicide (0, 5 and 30 mg kg(-1)). The treated soils were incubated at 20 °C for 6 weeks, and samples were taken after 1, 3 and 6 weeks of incubation. In soils with 5 mg kg(-1) MCPA, about 50 % of the MCPA was dissipated within 1 week of the incubation. Almost complete dissipation of the herbicide had occurred by the end of the incubation with no differences between the three litter amendments. At the higher concentration (30 mg kg(-1)), MCPA endured longer in the soil, with only 31 % of the initial amount being removed at the end of the experiment in the absence of litter. Litter addition greatly increased the dissipation rate with 70 and 80 % of the herbicide being dissipated in the 10 and 20 g kg(-1) litter treatments, respectively. Signs of toxic effects of MCPA on soil bacteria were observed from related phospholipid fatty acid (PLFA) analyses, while fungi showed higher tolerance to the increased MCPA levels. The abundance of bacterial tfdA genes in soil increased with the co-occurrence of litter and high MCPA concentration, indicating the importance of substrate availability in fostering MCPA-degrading bacteria and thereby improving the potential for removal of MCPA in the environment.

The impact of chemical pollution on the resilience of soils under multiple stresses: A conceptual framework for future research.

Soils are faced with man-made chemical stress factors, such as the input of organic or metal-containing pesticides, in combination with non-chemical stressors like soil compaction and natural disturbance like drought. Although multiple stress factors are typically co-occurring in soil ecosystems, research in soil sciences on this aspect is limited and focuses mostly on single structural or functional endpoints. A mechanistic understanding of the reaction of soils to multiple stressors is currently lacking. Based on a review of resilience theory, we introduce a new concept for research on the ability of polluted soil (xenobiotics or other chemical pollutants as one stressor) to resist further natural or anthropogenic stress and to retain its functions and structure. There is strong indication that pollution as a primary stressor will change the system reaction of soil, i.e., its resilience, stability and resistance. It can be expected that pollution affects the physiological adaption of organisms and the functional redundancy of the soil to further stress. We hypothesize that the recovery of organisms and chemical-physical properties after impact of a follow-up stressor is faster in polluted soil than in non-polluted soil, i.e., polluted soil has a higher dynamical stability (dynamical stability=1/recovery time), whereas resilience of the contaminated soil is lower compared to that of not or less contaminated soil. Thus, a polluted soil might be more prone to change into another system regime after occurrence of further stress. We highlight this issue by compiling the literature exemplarily for the effects of Cu contamination and compaction on soil functions and structure. We propose to intensify research on effects of combined stresses involving a multidisciplinary team of experts and provide suggestions for corresponding experiments. Our concept offers thus a framework for system level analysis of soils paving the way to enhance ecological theory.

Impacts of temperature increase and change in precipitation pattern on crop yield and yield quality of barley.

Spring barley was grown in a field experiment under moderately elevated soil temperature and changed summer precipitation (amount and frequency). Elevated temperature affected the performance and grain quality characteristics more significant than changes in rainfall. Except for the decrease in thousand grain weight, warming had no impacts on aboveground biomass and grain yield traits. In grains, several proteinogenic amino acids concentrations were increased, whereas their composition was only slightly altered. Concentration and yield of total protein remained unaffected under warming. The concentrations of total non-structural carbohydrates, starch, fructose and raffinose were lower in plants grown at high temperatures, whereas maltose was higher. Crude fibre remained unaffected by warming, whereas concentrations of lipids and aluminium were reduced. Manipulation of precipitation only marginally affected barley grains: amount reduction increased the concentrations of several minerals (sodium, copper) and amino acids (leucine). The projected climate changes may most likely affect grain quality traits of interest for different markets and utilisation requirements.

Succession of bacterial and fungal 4-chloro-2-methylphenoxyacetic acid degraders at the soil-litter interface.

Phenoxyacetic acids can be degraded by diverse soil microorganisms. Nevertheless, we miss information about the succession of 4-chloro-2-methylphenoxyacetic acid (MCPA) degraders in micro-environments of soils as well as specific functions of different microbial groups during MCPA degradation. We studied MCPA degradation at the soil-litter interface in a microcosm experiment and followed the succession of different degrader populations by quantifying the abundance of 16S rRNA genes as well as, the fungal ITS fragment and the functional genes tfdA (in total and divided into three classes) and cadA. Adjacent to the litter layer, a dynamic depletion zone of MCPA indicated that the litter effect on MCPA degradation depends on substrate availability and the affected soil volume. The increase of the tfdA class III and cadA genes was linked to MCPA mineralisation. Total abundance of tfdA genes was dominated by class I MCPA degraders and did not reflect MCPA degradation potential of the soil. Litter addition induced the development of pioneer and late-stage fungal communities, which were probably both involved in MCPA degradation. The results underline the importance of the ecological behaviour of different degrader populations for the understanding of herbicide degradation in soils.

Influence of land-use intensity on the spatial distribution of N-cycling microorganisms in grassland soils.

A geostatistical approach using replicated grassland sites (10 m × 10 m) was applied to investigate the influence of grassland management, i.e. unfertilized pastures and fertilized mown meadows representing low and high land-use intensity (LUI), on soil biogeochemical properties and spatial distributions of ammonia-oxidizing and denitrifying microorganisms in soil. Spatial autocorrelations of the different N-cycling communities ranged between 1.4 and 7.6 m for ammonia oxidizers and from 0.3 m for nosZ-type denitrifiers to scales >14 m for nirK-type denitrifiers. The spatial heterogeneity of ammonia oxidizers and nirS-type denitrifiers increased in high LUI, but decreased for biogeochemical properties, suggesting that biotic and/or abiotic factors other than those measured are driving the spatial distribution of these microorganisms at the plot scale. Furthermore, ammonia oxidizers (amoA ammonia-oxidizing archaea and amoA ammonia-oxidizing bacteria) and nitrate reducers (napA and narG) showed spatial coexistence, whereas niche partitioning was found between nirK- and nirS-type denitrifiers. Together, our results indicate that spatial analysis is a useful tool to characterize the distribution of different functional microbial guilds with respect to soil biogeochemical properties and land-use management. In addition, spatial analyses allowed us to identify distinct distribution ranges indicating the coexistence or niche partitioning of N-cycling communities in grassland soil.