Holocene land-cover reconstructions for studies on land cover-climate feedbacks

The major objectives of this paper are: (1) to review the pros and cons of the scenarios of past anthropogenic land cover change (ALCC) developed during the last ten years, (2) to discuss issues related to pollen-based reconstruction of the past land-cover and introduce a new method, REVEALS (Regional Estimates of VEgetation Abundance from Large Sites), to infer long-term records of past land-cover from pollen data, (3) to present a new project (LANDCLIM: LAND cover – CLIMate interactions in NW Europe during the Holocene) currently underway, and show preliminary results of REVEALS reconstructions of the regional land-cover in the Czech Republic for five selected time windows of the Holocene, and (4) to discuss the implications and future directions in climate and vegetation/land-cover modeling, and in the assessment of the effects of human-induced changes in land-cover on the regional climate through altered feedbacks. The existing ALCC scenarios show large discrepancies between them, and few cover time periods older than AD 800. When these scenarios are used to assess the impact of human land-use on climate, contrasting results are obtained. It emphasizes the need for methods such as the REVEALS model-based land-cover reconstructions. They might help to fine-tune descriptions of past land-cover and lead to a better understanding of how long-term changes in ALCC might have influenced climate. The REVEALS model is demonstrated to provide better estimates of the regional vegetation/land-cover changes than the traditional use of pollen percentages. This will achieve a robust assessment of land cover at regional- to continental-spatial scale throughout the Holocene. We present maps of REVEALS estimates for the percentage cover of 10 plant functional types (PFTs) at 200 BP and 6000 BP, and of the two open-land PFTs "grassland" and "agricultural land" at five time-windows from 6000 BP to recent time. The LANDCLIM results are expected to provide crucial data to reassess ALCC estimates for a better understanding of the land suface-atmosphere interactions.

Transient simulations of the carbon and nitrogen dynamics in northern peatlands: from the Last Glacial Maximum to the 21st century

The development of northern high-latitude peatlands played an important role in the carbon (C) balance of the land biosphere since the Last Glacial Maximum (LGM). At present, carbon storage in northern peatlands is substantial and estimated to be 500 ± 100 Pg C (1 Pg C = 1015 g C). Here, we develop and apply a peatland module embedded in a dynamic global vegetation and land surface process model (LPX-Bern 1.0). The peatland module features a dynamic nitrogen cycle, a dynamic C transfer between peatland acrotelm (upper oxic layer) and catotelm (deep anoxic layer), hydrology- and temperature-dependent respiration rates, and peatland specific plant functional types. Nitrogen limitation down-regulates average modern net primary productivity over peatlands by about half. Decadal acrotelm-to-catotelm C fluxes vary between −20 and +50 g C m−2 yr−1 over the Holocene. Key model parameters are calibrated with reconstructed peat accumulation rates from peat-core data. The model reproduces the major features of the peat core data and of the observation-based modern circumpolar soil carbon distribution. Results from a set of simulations for possible evolutions of northern peat development and areal extent show that soil C stocks in modern peatlands increased by 365–550 Pg C since the LGM, of which 175–272 Pg C accumulated between 11 and 5 kyr BP. Furthermore, our simulations suggest a persistent C sequestration rate of 35–50 Pg C per 1000 yr in present-day peatlands under current climate conditions, and that this C sink could either sustain or turn towards a source by 2100 AD depending on climate trajectories as projected for different representative greenhouse gas concentration pathways.

Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

Grazing response patterns indicate isolation of semi-natural European grasslands

Identifying drivers of species diversity is a major challenge in understanding and predicting the dynamics of species-rich semi-natural grasslands. In particular in temperate grasslands changes in land use and its consequences, i.e. increasing fragmentation, the on-going loss of habitat and the declining importance of regional processes such as seed dispersal by livestock, are considered key drivers of the diversity loss witnessed within the last decades. It is a largely unresolved question to what degree current temperate grassland communities already reflect a decline of regional processes such as longer distance seed dispersal. Answering this question is challenging since it requires both a mechanistic approach to community dynamics and a sufficient data basis that allows identifying general patterns. Here, we present results of a local individual- and trait-based community model that was initialized with plant functional types (PFTs) derived from an extensive empirical data set of species-rich grasslands within the `Biodiversity Exploratories' in Germany. Driving model processes included above- and belowground competition, dynamic resource allocation to shoots and roots, clonal growth, grazing, and local seed dispersal. To test for the impact of regional processes we also simulated seed input from a regional species pool. Model output, with and without regional seed input, was compared with empirical community response patterns along a grazing gradient. Simulated response patterns of changes in PFT richness, Shannon diversity, and biomass production matched observed grazing response patterns surprisingly well if only local processes were considered. Already low levels of additional regional seed input led to stronger deviations from empirical community pattern. While these findings cannot rule out that regional processes other than those considered in the modeling study potentially play a role in shaping the local grassland communities, our comparison indicates that European grasslands are largely isolated, i.e. local mechanisms explain observed community patterns to a large extent.

(Table 1) Species richness and diversity of tundra plants in herbivore exclosure and control plots near Barrow, Alaska

To determine the role lemmings play in structuring plant communities and their contribution to the 'greening of the Arctic', we measured plant cover and biomass in 50 + year old lemming exclosures and control plots in the coastal tundra near Barrow, Alaska. The response of plant functional types to herbivore exclusion varied among land cover types. In general, the abundance of lichens and bryophytes increased with the exclusion of lemmings, whereas graminoids decreased, although the magnitude of these responses varied among land cover types. These results suggest that sustained lemming activity promotes a higher biomass of vascular plant functional types than would be expected without their presence and highlights the importance of considering herbivory when interpreting patterns of greening in the Arctic. In light of the rapid environmental change ongoing in the Arctic and the potential regional to global implications of this change, further exploration regarding the long-term influence of arvicoline rodents on ecosystem function (e.g. carbon and energy balance) should be considered a research priority.

Initial vegetation, and spore density, diversity and evenness of different arbuscular mycorrhizal fungi morphotaxa 29,36 and 39 months after experimental vegetation removal, central Argentina

1. Dominant plant functional types (PFTs) are expected to be primary determinants of communities of other above- and below-ground organisms. Here, we report the effects of the experimental removal of different PFTs on arbuscular mycorrhizal fungi (AMF) communities in a shrubland ecosystem in central Argentina. 2. On the basis of the biomass-ratio hypothesis and plant resource use strategy theory, we expected the effect of removal of PFTs on AMF colonization and spores to be proportional to the biomass removed and to be stronger when more conservative PFTs were removed. The treatments applied were: undisturbed control (no plant removed), disturbed control (mechanical disturbance), no shrub (removal of deciduous shrubs), no perennial forb (removal of perennial forbs), no graminoid (removal of graminoids) and no annual forb (removal of annual forbs). AMF colonization was assessed after 5,17 and 29 months. Total density of AMF spores, richness and evenness of morpho-taxa, and AMF functional groups were quantified after 5,17,29,36 and 39 months. 3. Five months after the initial removal we found a significant reduction in total AMF colonization in all plots subjected to PFT removals and in the disturbed control plots, as compared with the undisturbed controls. This effect disappeared afterwards and no subsequent effect on total colonization and colonization by arbuscules was observed. In contrast, a significant increase in colonization by vesicles was observed in months 17 and 29, mainly in no graminoid plots. In general, treatments did not significantly affect AMF spores in the soil. On the other hand, no annual forb promoted transient (12-18 months) higher ammonia availability, and no shrub promoted lower nitrate availability in the longer term (24-28 months). 4. Synthesis. Our experiment, the first to investigate the effects of the removal of different PFTs on AMF communities in natural ecosystems, indicates that AMF communities are resilient to changes in the soil and in the functional composition of vegetation. Furthermore, it does not provide consistent evidence in support of the biomass-ratio hypothesis or differential trait-based direct or indirect effects of different PFTs on AMF in this particular system.

Tropical climate and vegetation simulations during the Heinrich event 1 using an Earth System Model of Intermediate Complexity (EMIC) - the University of Victoria Earth System-Climate Model (UVic ESCM)

Abrupt climate changes from 18 to 15 thousand years before present (kyr BP) associated with Heinrich Event 1 (HE1) had a strong impact on vegetation patterns not only at high latitudes of the Northern Hemisphere, but also in the tropical regions around the Atlantic Ocean. To gain a better understanding of the linkage between high and low latitudes, we used the University of Victoria (UVic) Earth System-Climate Model (ESCM) with dynamical vegetation and land surface components to simulate four scenarios of climate-vegetation interaction: the pre-industrial era, the Last Glacial Maximum (LGM), and a Heinrich-like event with two different climate backgrounds (interglacial and glacial). We calculated mega-biomes from the plant-functional types (PFTs) generated by the model to allow for a direct comparison between model results and palynological vegetation reconstructions. Our calculated mega-biomes for the pre-industrial period and the LGM corresponded well with biome reconstructions of the modern and LGM time slices, respectively, except that our pre-industrial simulation predicted the dominance of grassland in southern Europe and our LGM simulation resulted in more forest cover in tropical and sub-tropical South America. The HE1-like simulation with a glacial climate background produced sea-surface temperature patterns and enhanced inter-hemispheric thermal gradients in accordance with the "bipolar seesaw" hypothesis. We found that the cooling of the Northern Hemisphere caused a southward shift of those PFTs that are indicative of an increased desertification and a retreat of broadleaf forests in West Africa and northern South America. The mega-biomes from our HE1 simulation agreed well with paleovegetation data from tropical Africa and northern South America. Thus, according to our model-data comparison, the reconstructed vegetation changes for the tropical regions around the Atlantic Ocean were physically consistent with the remote effects of a Heinrich event under a glacial climate background.

Stable hydrogen and carbon isotope records and concentration data for four long chain n-alkanes from core GIK16160-3 derived from the Zambezi delta in the Mozambique Channel during the last 37,000 years

Hydrogen isotope values (dD) of sedimentary terrestrial leaf wax such as n-alkanes or n-acids have been used to map and understand past changes in rainfall amount in the tropics because dD of precipitation is commonly assumed as the first order controlling factor of leaf wax dD. Plant functional types and their photosynthetic pathways can also affect leaf wax dD but these biological effects are rarely taken into account in paleo studies relying on this rainfall proxy. To investigate how biological effects may influence dD values we here present a 37,000-year old record of dD and stable carbon isotopes (d13C) measured on four n-alkanes (n-C27, n-C29, n-C31, n-C33) from a marine sediment core collected off the Zambezi River mouth. Our paleo d13C records suggest that each individual n-alkanes had different C3/C4 proportional contributions. n-C29 was mostly derived from a C3 dicots (trees, shrubs and forbs) dominant vegetation throughout the entire record. In contrast, the longer chain n-C33 and n-C31 were mostly contributed by C4 grasses during the Glacial period but shifted to a mixture of C4 grasses and C3 dicots during the Holocene. Strong correlations between dD and d13C values of n-C33 (correlation coefficient R2 = 0.75, n = 58) and n-C31 (R2 = 0.48, n = 58) suggest that their dD values were strongly influenced by changes in the relative contributions of C3/C4 plant types in contrast to n-C29 (R2 = 0.07, n = 58). Within regions with variable C3/C4 input, we conclude that dD values of n-C29 are the most reliable and unbiased indicator for past changes in rainfall, and that dD and d13C values of n-C31 and n-C33 are sensitive to C3/C4 vegetation changes. Our results demonstrate that a robust interpretation of palaeohydrological data using n-alkane dD requires additional knowledge of regional vegetation changes from which nalkanes are synthesized, and that the combination of dD and d13C values of multiple n-alkanes can help to differentiate biological effects from those related to the hydrological cycle.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Potential Arctic tundra vegetation shifts in response to changing temperature, precipitation and permafrost thaw

Over the past decades, vegetation and climate have changed significantly in the Arctic. Deciduous shrub cover is often assumed to expand in tundra landscapes, but more frequent abrupt permafrost thaw resulting in formation of thaw ponds could lead to vegetation shifts towards graminoid-dominated wetland. Which factors drive vegetation changes in the tundra ecosystem are still not sufficiently clear. In this study, the dynamic tundra vegetation model, NUCOM-tundra (NUtrient and COMpetition), was used to evaluate the consequences of climate change scenarios of warming and increasing precipitation for future tundra vegetation change. The model includes three plant functional types (moss, graminoids and shrubs), carbon and nitrogen cycling, water and permafrost dynamics and a simple thaw pond module. Climate scenario simulations were performed for 16 combinations of temperature and precipitation increases in five vegetation types representing a gradient from dry shrub-dominated to moist mixed and wet graminoid-dominated sites. Vegetation composition dynamics in currently mixed vegetation sites were dependent on both temperature and precipitation changes, with warming favouring shrub dominance and increased precipitation favouring graminoid abundance. Climate change simulations based on greenhouse gas emission scenarios in which temperature and precipitation increases were combined showed increases in biomass of both graminoids and shrubs, with graminoids increasing in abundance. The simulations suggest that shrub growth can be limited by very wet soil conditions and low nutrient supply, whereas graminoids have the advantage of being able to grow in a wide range of soil moisture conditions and have access to nutrients in deeper soil layers. Abrupt permafrost thaw initiating thaw pond formation led to complete domination of graminoids. However, due to increased drainage, shrubs could profit from such changes in adjacent areas. Both climate and thaw pond formation simulations suggest that a wetter tundra can be responsible for local shrub decline instead of shrub expansion.

Is the maximum velocity of carboxylation (Vcmax) well adjusted for deciduous shrubs in DGVMs: a case study for the Caatinga Biome in Brazil.

Vcmax is the rate of maximum velocity of carboxylation of plants and is considered one of the most critical parameters for changes in vegetation in face of global changes and it has a direct impact on gross primary productivity. Physiological processes are considered the main sources of uncertainties in dynamic global vegetation models (DGVMs). The Caatinga biome, in the semiarid region of northeastern Brazil, is extremely important due to its biodiversity and endemism. In a field work realized in an area of preserved Caatinga forest, measurements of carbon assimilation (in response to light and CO2) were performed on 11 individuals of a native species. These results of Vcmax measurements in Caatinga were compared with parameterization of models, revealing that Vcmax is not well adjusted in several DGVMs. Also, the values obtained in the Caatinga field experiments were very close to empirical values obtained in the Northern hemisphere (Austria). These ecophysiological measurements can contribute in understanding of this biome

Canopy-scale biophysical controls of transpiration and evaporation in the Amazon Basin.

Canopy and aerodynamic conductances (gC and gA) are two of the key land surface biophysical variables that control the land surface response of land surface schemes in climate models. Their representation is crucial for predicting transpiration (λET) and evaporation (λEE) flux components of the terrestrial latent heat flux (λE), which has important implications for global climate change and water resource management. By physical integration of radiometric surface temperature (TR) into an integrated framework of the Penman?Monteith and Shuttleworth?Wallace models, we present a novel approach to directly quantify the canopy-scale biophysical controls on λET and λEE over multiple plant functional types (PFTs) in the Amazon Basin. Combining data from six LBA (Large-scale Biosphere-Atmosphere Experiment in Amazonia) eddy covariance tower sites and a TR-driven physically based modeling approach, we identified the canopy-scale feedback-response mechanism between gC, λET, and atmospheric vapor pressure deficit (DA), without using any leaf-scale empirical parameterizations for the modeling. The TR-based model shows minor biophysical control on λET during the wet (rainy) seasons where λET becomes predominantly radiation driven and net radiation (RN) determines 75 to 80 % of the variances of λET. However, biophysical control on λET is dramatically increased during the dry seasons, and particularly the 2005 drought year, explaining 50 to 65 % of the variances of λET, and indicates λET to be substantially soil moisture driven during the rainfall deficit phase. Despite substantial differences in gA between forests and pastures, very similar canopy?atmosphere "coupling" was found in these two biomes due to soil moisture-induced decrease in gC in the pasture. This revealed the pragmatic aspect of the TR-driven model behavior that exhibits a high sensitivity of gC to per unit change in wetness as opposed to gA that is marginally sensitive to surface wetness variability. Our results reveal the occurrence of a significant hysteresis between λET and gC during the dry season for the pasture sites, which is attributed to relatively low soil water availability as compared to the rainforests, likely due to differences in rooting depth between the two systems. Evaporation was significantly influenced by gA for all the PFTs and across all wetness conditions. Our analytical framework logically captures the responses of gC and gA to changes in atmospheric radiation, DA, and surface radiometric temperature, and thus appears to be promising for the improvement of existing land?surface?atmosphere exchange parameterizations across a range of spatial scales.

Canopy-scale biophysical controls of transpiration and evaporation in the Amazon Basin.

Canopy and aerodynamic conductances (gC and gA) are two of the key land surface biophysical variables that control the land surface response of land surface schemes in climate models. Their representation is crucial for predicting transpiration (?ET) and evaporation (?EE) flux components of the terrestrial latent heat flux (?E), which has important implications for global climate change and water resource management. By physical integration of radiometric surface temperature (TR) into an integrated framework of the Penman?Monteith and Shuttleworth?Wallace models, we present a novel approach to directly quantify the canopy-scale biophysical controls on ?ET and ?EE over multiple plant functional types (PFTs) in the Amazon Basin. Combining data from six LBA (Large-scale Biosphere-Atmosphere Experiment in Amazonia) eddy covariance tower sites and a TR-driven physically based modeling approach, we identified the canopy-scale feedback-response mechanism between gC, ?ET, and atmospheric vapor pressure deficit (DA), without using any leaf-scale empirical parameterizations for the modeling. The TR-based model shows minor biophysical control on ?ET during the wet (rainy) seasons where ?ET becomes predominantly radiation driven and net radiation (RN) determines 75 to 80?% of the variances of ?ET. However, biophysical control on ?ET is dramatically increased during the dry seasons, and particularly the 2005 drought year, explaining 50 to 65?% of the variances of ?ET, and indicates ?ET to be substantially soil moisture driven during the rainfall deficit phase. Despite substantial differences in gA between forests and pastures, very similar canopy?atmosphere "coupling" was found in these two biomes due to soil moisture-induced decrease in gC in the pasture. This revealed the pragmatic aspect of the TR-driven model behavior that exhibits a high sensitivity of gC to per unit change in wetness as opposed to gA that is marginally sensitive to surface wetness variability. Our results reveal the occurrence of a significant hysteresis between ?ET and gC during the dry season for the pasture sites, which is attributed to relatively low soil water availability as compared to the rainforests, likely due to differences in rooting depth between the two systems. Evaporation was significantly influenced by gA for all the PFTs and across all wetness conditions. Our analytical framework logically captures the responses of gC and gA to changes in atmospheric radiation, DA, and surface radiometric temperature, and thus appears to be promising for the improvement of existing land?surface?atmosphere exchange parameterizations across a range of spatial scales.

Plant regeneration functional groups modulate the response to fire of soil enzyme activities in a Mediterranean shrubland

Soil enzymes are critical to soil nutrient cycling function but knowledge on the factors that control their response to major disturbances such as wildfires remains very limited. We evaluated the effect of fire-related plant functional traits (resprouting and seeding) on the resistance and resilience to fire of two soil enzyme activities involved in phosphorus and carbon cycling (acid phosphatase and β-glucosidase) in a Mediterranean shrublands in SE Spain. Using experimental fires, we compared four types of shrubland microsites: SS (vegetation patches dominated by seeder species), RR (patches dominated by resprouter species), SR (patches co-dominated by seeder and resprouter species), and IP (shrub interpatches). We assessed pre- and post-fire activities of the target soil enzymes, available P, soil organic C, and plant cover dynamics over three years after the fire. Post-fire regeneration functional groups (resprouter, seeder) modulated both pre- and post-fire activity of acid phosphatase and β-glucosidase, with higher activity in RR and SR patches than in SS patches and IP. However, we found no major differences in enzyme resistance and resilience between microsite types, except for a trend towards less resilience in SS patches. Fire similarly reduced the activity of both enzymes. However, acid phosphatase and β-glucosidase showed contrasting post-fire dynamics. While β-glucosidase proved to be rather resilient to fire, fully recovering three years after fire, acid phosphatase showed no signs of recovery in that period. Overall, the results indicate a positive influence of resprouter species on soil enzyme activity that is very resistant to fire. Long-lasting decrease in acid phosphatase activity probably resulted from the combined effect of P availability and post-fire drought. Our results provide insights on how plant functional traits modulate soil biochemical and microbiological response to fire in Mediterranean fire-prone shrublands.

Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities

The controls on aboveground community composition and diversity have been extensively studied, but our understanding of the drivers of belowground microbial communities is relatively lacking, despite their importance for ecosystem functioning. In this study, we fitted statistical models to explain landscape-scale variation in soil microbial community composition using data from 180 sites covering a broad range of grassland types, soil and climatic conditions in England. We found that variation in soil microbial communities was explained by abiotic factors like climate, pH and soil properties. Biotic factors, namely community- weighted means (CWM) of plant functional traits, also explained variation in soil microbial communities. In particular, more bacterial-dominated microbial communities were associated with exploitative plant traits versus fungal-dominated communities with resource-conservative traits, showing that plant functional traits and soil microbial communities are closely related at the landscape scale.