Different Land Use Intensities in Grassland Ecosystems Drive Ecology of Microbial Communities Involved in Nitrogen Turnover in Soil

Understanding factors driving the ecology of N cycling microbial communities is of central importance for sustainable land use. In this study we report changes of abundance of denitrifiers, nitrifiers and nitrogen-fixing microorganisms (based on qPCR data for selected functional genes) in response to different land use intensity levels and the consequences for potential turnover rates. We investigated selected grassland sites being comparable with respect to soil type and climatic conditions, which have been continuously treated for many years as intensely used meadows (IM), intensely used mown pastures (IP) and extensively used pastures (EP), respectively. The obtained data were linked to above ground biodiversity pattern as well as water extractable fractions of nitrogen and carbon in soil. Shifts in land use intensity changed plant community composition from systems dominated by s-strategists in extensive managed grasslands to c-strategist dominated communities in intensive managed grasslands. Along the different types of land use intensity, the availability of inorganic nitrogen regulated the abundance of bacterial and archaeal ammonia oxidizers. In contrast, the amount of dissolved organic nitrogen determined the abundance of denitrifiers (nirS and nirK). The high abundance of nifH carrying bacteria at intensive managed sites gave evidence that the amounts of substrates as energy source outcompete the high availability of inorganic nitrogen in these sites. Overall, we revealed that abundance and function of microorganisms involved in key processes of inorganic N cycling (nitrification, denitrification and N fixation) might be independently regulated by different abiotic and biotic factors in response to land use intensity.


Introduction
Soils provide a large number of ecosystem services including plant growth, carbon sequestration, degradation of xenobiotics and safeguarding of drinking water resources. Most of these functions are closely linked to the soil microbiome and its activity pattern [1,2,3]. Therefore many attempts were made to identify soil borne microbial communities as key drivers of ecological processes and describe factors that drive the abundance and diversity of selected functional communities [4]. Despite the high heterogeneity of soil microbes in time and space, it has become possible to figure out one general conclusions from these studies: Besides site-specific parameters, for example soil texture or climatic conditions, the type of land management and land use intensity has been identified as a major driver for microbial performance in soil [5,6,7,8].
Recently, the effects of land use changes have been studied mainly focusing on (i) conversion of grassland to forest or vice versa [9], (ii) alterations in tillage management [10,11,12], (iii) changes in crop rotation [13] or (iv) modifications in fertilizer quality [14]. However, studies addressing questions related to consequences of changes in land use intensity on the soil microbiome are rare, although in many parts of the world we are facing a tremendous increase in land use intensity, due to the demands of bioeconomy (production of food, feed, fuel and fiber).
This intensification is also frequently observed in grassland ecosystems. While in the past sites have been used extensively as pastures, nowadays up to four times per season, the same areas are managed as meadows for hay production and silage, entailing an intensive application of organic and inorganic fertilizers. Differences in intensity of agricultural practice like mowing, grazing and fertilization lead to changes in plant composition [15,16,17,13], microclimate, soil quality and hence to changes on macro-as well as micro-scale habitats. For some soil animals the impact of such changes is well known [18,19,20,21] but data on microbial communities in soil is rare. For example, [6] compared diversity pattern of microbial community involved in nitrogen fixation, denitrification and nitrification in grassland ecosystems under different management intensities. This study clearly demonstrated changes in diversity pattern of single functional groups involved in nitrogen transformation on low diverse grassland sites. However, this study did not address questions how land use intensity influences the abundance and activity of selected functional groups of microbes in soil and thus changes turnover processes and rates.
The aim of the present study was to characterize microbial communities responsible for key processes in the inorganic nitrogen cycle (nitrification, denitrification and N-fixation) in grasslands of different land use intensity and relate these results to the aboveground biodiversity of plants as well as important belowground properties (water extractable carbon and nitrogen fractions). Overall, we postulated that nitrogen cycling at extensively used sites is mainly driven by nitrogen fixation and internal nitrogen turnover is highly efficient resulting in low denitrification rates. In contrast, denitrification may play a more important role in nitrogen turnover at intensely used sites due to higher amounts of nitrogen available in soil and intensive plant growth (resulting in higher root exudation rates and increased microbial activity).

Experimental Setup and Sampling
Experiments were carried out in the frame of the German Biodiversity Exploratories [22], which form an ideal platform for such type of studies, as here for the first time gradients in land use intensity were defined on a large scale in three regions in Germany. For the present study soil samples was taken in 2008 in the southernmost Exploratory ''Schwä bische Alb'' which covers more than 45,000 ha of the state of Baden-Württemberg in SW-Germany. The mean annual precipitation in this area ranged from 938-963 mm, whereas the annual mean temperature was around 7uC. All sampled soils have been described as a Rendzic Leptosols with a clayey or loamy texture and a pH-value between 5.7 and 6.9. A more detailed soil description is given in table 1. Field work permits were given by the responsible state environmental offices of the state of Baden Württemberg (according to 1 72 BbgNatSchG).
From each site five sampling replicates were taken from the upper 10 cm; each replicate consist of five pooled bulk soil cores (d = 5.5 cm) taken with a soil auger. All samples were cooled directly after sampling for DNA based analyses at 220uC and at 4uC for measurements of enzyme activities and soil parameter.

Plant Diversity
At all investigated sites the vegetation has been recorded on an area of 4 6 4 m. Plants were identified on taxa level and their percentage cover was estimated separately for the shrub layer (0-5 m woody species) and the herbaceous layer (including phanerophyte seedlings), respectively (data here not shown). The ecological strategy type of each vascular plant according to [23] were determined using the ''Biolflor'' data base (www.biolflor.de). This concept describing the general limits to ecology and evolution based on the trade-off that organisms face when the resources they gain from the environment are allocated between either growth, maintenance or regeneration -known as the universal three-way trade-off.: (C) the survival of the individual using traits that maximise resource acquisition and resource control in consistently productive niches, (S) individual survival via maintenance of metabolic performance in variable and unproductive niches, or (R) rapid gene propagation via rapid completion of the lifecycle and regeneration in niches where events are frequently lethal to the individual. Based on this assumption the csr triangle has been defined, in which each plant could be integrated according to its individual life type strategies. Furthermore, site conditions were characterized by calculating the mean indicator values based on plant species composition according to [24]. Ellenberg indicator values represent a well-established method for bio-indication of a range of environmental parameters.

Chemical Analyses of Soils
CaCl 2 extracts of the soil samples were prepared to analyze carbon and nitrogen contents in water extractable fraction. Therefore soil was mixed with 0.01 M CaCl 2 at a ratio of 1:2 [25] and homogenized by overhead shaken for 2 h followed by a filtration step. Water extractable organic carbon (WEOC) and total nitrogen bound (TNb) were measured using a DIMA-TOC 100 (DIMATEC Analysentechnik, Essen, Germany), nitrate and ammonium by using the commercial kits Nanocolor Nu50 (NO 3 -N) and Nu3 (NH 4 +-N) (Macherey-Nagel, Düren, Germany).

Potential Nitrification and Denitrification Rate in Soils
Potential nitrification rates (PNR) were determined via micro titer plate assay as described by [26]. This method is based on the quantification of nitrite formed by ammonia oxidation. In brief, 2.5 g of fresh soil (three technical replicates for each soil sample), 50 ml sodium chlorate (1.5 M) and 10 ml ammonium sulfate (10 mM) were shaken for 5 h at 20uC whereas the added sodium chlorate should repress the oxidation of nitrite to nitrate. To stop the oxidation of ammonium 2.5 ml potassium chloride was added and incubated for 20 min. After short centrifugation 150 ml of the hydrous supernatant, 90 ml ammonium chloride (0.2 M) and 60 ml Griess-reagent (0.002 M naphthylenediamine dihydrochloride, 0.06 M sulphanilamide and 2.5 M phosphoric acid) were pipetted into the micro titer plates and the color change were detected via plate spectrometer at a absorbance wavelength of 540 nm (SpectraMax 340, MWG BIOTECH, Ebersberg, Germany). Nitrite concentrations of non-incubated soil samples served as control.
Potential denitrification activity (DEA) was determined using the anaerobic slurry technique as described by [27], slightly modified by [28]. Soil slurries, containing 10 g field-fresh soil and 10 ml of a mixture of 1 mM glucose/1 mM potassium nitrate (KNO 3 ), were flushed with helium in airtight serum bottles. Subsequent acetylene was added to the slurries and the bottles were shaken via overhead shaker at 25uC for two hours. Every hour gas samples of headspace atmosphere were taken with a gastight syringe, and N 2 O concentrations measured using a gas chromatograph (GC-14B, Shimadzu, Japan).

Nucleic Acid Extraction and Quantification
Genomic DNA was extracted from 0.5 g bulk soil (wet weight) using FastDNA Spin Kit for soil (MP Biomedicals, Germany) according to manufacturer's protocol. Quality and Quantity of DNA extracts were determined with Nanodrop 1000 Spectrophotometer (Peqlab, Germany). A260/A280 ratios were approx 1.9. As the used kit is known to highly influence A230 due to the used binding matrix which has its absorption maximum at 230 nm, we used qPCR dilution series to exclude inhibition effects and did not calculate A230/260 ratios.

Real-time PCR Assay
Quantitative real-time PCR (qPCR) was conducted on a 7300 Real-Time PCR System (Applied Biosystems, Germany) using SybrGreen as fluorescent dye to estimate the abundance of ammonia oxidizing bacteria and archaea, nitrogen fixing microbes as well as nitrite reducers. The marker genes used as well as the corresponding conditions for qPCR are listed in table 2.
In a pre-experiment a dilution of the DNA of 1:64 turned out to be best suited avoid inhibition of PCR, e.g. by co-extracted humic substances (data not shown). Quantitative real-time PCR was performed in 96-well plates (Applied Biosystems). Reaction mixtures with total volume of 25 ml were composed of 12.5 ml of Power SybrGreen PCR Master Mix (Applera, Germany), 15 mg bovine serum albumine (Sigma-Aldrich, Germany), gene specific concentrations of the forward and reverse primer (Metabion Germany) (see table 2) and 2 ml of diluted DNA extract. In the case of the genes nirS and nirK dimethyl sulfoxide (DMSO) was additionally used (final concentration of 0.3 M). The PCR reaction was started with a hot start of 94uC for 10 min according to the manufacturer's instruction of the SybrGreen master mix, followed by 40 cycles with specific temperature profile according to the gene targeted (see table 2). Data collection was performed at each elongation step. The purity of the amplicons was checked by melting curve analysis and the presence of a unique band of the expected size in a 1.5% agarose gel stained with ethidium bromide. Standard curves was obtained using serial dilutions of plasmids DNA (10 6 -10 1 gene copies/ml) containing the respective cloned gene (see table 2). The amplification efficiencies were calculated from the formula Eff = [10 (-1/slope) -1] and resulted in values from 97% to 83%.

Statistical Analysis
All results were related on the basis of one gram of soil dry weight (g dw -1 ). Using the R environment for statistical computing (http://www.R-project.org) univariate and multivar- iate methods were applied to test the effect of land use and season. Assumptions were tested by checking the residuals (equal variance in the groups, no outliers, normal distribution). To account for the existence of replicate groups (n = 5) in each case linear mixed-effect models (lme in package nlme) were fitted. Two factor lmes with pairwise comparisons were used to construct a table of p-values. The integrative multivariate data analysis was based on a between group analysis (between in package ade4). To investigate differences in plant diversity data a one factor ANOVA (lm in R) and subsequent pairwise t-tests with adjusted pvalues (pairwise.t.test with method holm) were applied.

Plant Diversity
The number of plant species on the investigated sites ranged between 17 and 58. Significantly lower numbers were observed on the intensively managed plots (24 respectively 25 for IM and IP) compared to the EP plots where 49 different species were found (p = 0.019). The coverage of vascular plant was between 102 and 170%, with highest plant coverage on the IM plots. The plant species composition showed a clear dependency on the land use intensity level (table S1). While grasses were dominant on the intensely used meadows (58 to 105% of coverage), the pastures (EP and IP) were more colonized by herbs (84 to 107% of coverage) indicating that grasses are more adapted to the disturbance caused by mowing. Considering the different ecological strategy types of the collected plants we found a clear trend for plants more linked to r-strategists on the extensively used sites (up to 96% of coverage) to c-strategists on the intensely used sites (up to 99% of the coverage).
The analysis of plants with indicator values for abiotic sites properties according to Ellenberg revealed a strong response of the plant community structure towards nutrient availability and soil water content (table 3). Most plant species from EP sites could be classified as indicators for nitrogen deficiency, whereas plants of intensely used sites were classified as indicators of moderate to strong nitrogen concentrations in soil. Additionally plants grown on the EP plots were better adapted to water shortage and matched well with plants being an indicator for moderate water shortage compared to plants obtained from the IP and IM plots which are indicators for moderate to well water conditions in soils.

Labile Soil Nitrogen and Carbon Pools
Soil ammonium content ranged from 0.19 to 4.1 mg N g 21 dw and was up to ten times higher in spring than in summer (p = 0.0444). In July values were consistently low (0.19 and 0.53 mg N g 21 dw) in dependent from the site, whereas in April (0.62 to 4.1 mg N g 21 dw) differences between sites were more pronounced, with a tendency for higher values in soil samples from EP and IM compared to IP (p = 0.2976). Nitrate concentrations were in the range of 1.61 and 54.5 mg N g 21 dw in April and 1.

Potential Nitrification and Denitrification Rates
The highest PNR values were measured in samples derived from IM sites. However transformation rates between the three investigated sites (AEG1 -3) were highly variable with values from 0.28 up to 2.55 mgNO 2 2 N g 21 dw h 21 ( figure 1 and table 4). PNR in samples from plots characterized as IP ranged between 0.26 and 1.91 mgNO 2 2 N g 21 dw h 21. and rates were neither significant different to IM nor to EP (p = 0.2944) plots. Again a high variation of the values, comparing the three sites under investigation (AEG4 -6), was visible. The lowest PNR was measured in samples derived from EP (0.0620.31 mgNO 2 2 N g 21 dw h 21 ), whereas a significant difference between EP and IM was observed (p = 0.0327). Overall a significant decrease of PNR from intensely to extensively managed sites was found (p = 0.0350). The sampling time point did not influence PNR in the investigated samples (p = 0.4882).
In contrast DEA showed an interdependency between the sampling time point and the land use intensity (p = 0.0346; figure 1  and table 4)). In spring the values for DEA in samples derived from IM and EP were up to the factor of three times higher (ranging from 0.

Abundances of Functional Genes
Gene copy numbers of bacterial amoA genes (AOB) increased significantly from the spring to the summer sampling at all intensively managed sites (IM and IP), ranging in April from 4.4610 6 to 5.5610 7 copies g -1 and in July from 2.2610 7 to 2.1 10 8 copies g -1 (figure 2). Gene copy numbers in samples derived from the extensively used plots were in the range of 7.7610 5 and 8.7610 6 g 21 and hence significantly lower compared to the intensively managed sites (table 4). In addition, no clear seasonal effect could be described for samples derived from EP sites, as on two plots (AEG7 and AEG8) a decrease from April to July was observed whereas on the third one an increase was found. Overall, AOB community size was influenced by season (p = 0.0293) and by land use intensity (p = 0.0004). Gene copy numbers for the archaeal amoA gene (AOA) ranged from 3.0610 6 and 2.5610 8 g 21 in April and from 8.5610 6 and 4.5610 8 copies g 21 in July; no statistically significant seasonal effect were proven (p = 0.5872), but a significant increase from intensely to extensively used plots was observed (see table 4). As for AOB also for AOA an increase from April to July of amoA gene abundances on the intensely used plots and a decrease on the extensively used sites was visible determined. Considering the ratios of AOA and AOB except the site AEG6 (AOA:AOB ratio ,1) AOA : AOB ratios between 2.1 and 16 were found indicating a dominance from AOA over AOB. However in soil samples from the intensely used sites the ratio decreased significantly from April to July, whereas in samples from EP in two plots no changes were observed between the two sampling time points and a increase from 2.2 to 16 was found in the 3 rd plot (AEG7).
Gene copy numbers for nirK ranged in April from 8.9610 7 to 3.7610 8 copies g 21 and increased in July ranging from 5.0610 8 to 3.4610 9 copies g 21 (p,0.0001; figure 2 and table 4). In April all plots of one land use category showed comparable abundance levels for nirK, whereas in July AEG3 (belonging to IM) revealed up to ten times lower copy numbers than AEG1 Table 3. Soil parameters and plant analysis.   April and between 4.0610 6 and 2.9610 7 copies g dw -1 in July. Also land use had a significant influence on the number of nirS gene copies. However nirS gene copy numbers were up to the factor of 500 lower compared to nirK. Gene copy numbers for nifH ranged from 1.8610 6 to 6.9610 7 copies g dw 21 (figure 2). The data did not reveal significant differences between the sampling time points (p = 0.639) but between the land use intensities (p = 0.0006). IP sites show significant higher nifH copy numbers (from 1.9610 7 to 6.9610 7 copies g dw 21 ) than the other two land use categories. The lowest values were detected on the EP plots, where quantities between 1.4610 6 and 9.0610 6 copies g dw 21 were measured.

Integrative Data Analysis
Based on all data obtained a between group analysis (BGA) was performed using the mean values of all five replicates from each site (Figure 3). The BGA revealed fertilization and other types of soil disturbance (e.g. soil compaction by machines or animals) are interconnected and explain the differences between the intensely (IM and IP) and the extensively (EP) managed grassland sites to a large extent (PCA1:43% PCA2:24%). Based on the BGA data, it can be assumed that the effects of fertilization were more dominant compared to the influence of grazing and mowing. Moreover, clear differences between the sampling time points appeared (samples from July located in the lowerpart, from April in the upper part). Overall abundance of nifH, amoA AOA and AOB genes as well as potential nitrification activity were more related to land use intensity, whereas ammonium content and dissolved organic carbon as well as DNA content and abundance of nirS and nirK genes were affected more by time.

Discussion
Plant Diversity [29] was among the first, who described a close link between land use intensity of grasslands and plant species composition. These observations have been confirmed in the present study. These results can be discussed on the basis of life strategies of plants and soil microorganisms: Whereas more plant species could be assigned to the exploitation group of so called c-strategists on the intensively managed sites, with high assimilation and growth rates if high nutrient levels were present in soil, the extensively managed plots were more characterized by ''csr-plants'' characterized by higher tolerance level to disturbance and nutrient stress but slower growth rates [23]. In addition to differences in species composition we found a significant reduction of plant species richness in response to land use intensification which has been also found in several other studies [30,31,32,33]. These losses were comparable to an over regional inventory by [34] who revealed that increased fertilization by 90-130 kg N ha 21 is linked to a 50% higher productivity and a decrease of species richness by 30%.

Land Use Intensity and Nitrification
Several studies which investigated the impacts of different management regimes on the nitrifying potential in soil [35,6,36] described higher PNR in response to higher intensities of land use of grassland ecosystems. These results were confirmed in our study where all intensely managed meadows showed higher PNR values compared to the extensively managed plots. Interestingly despite an increase in amoA gene copy numbers in samples from the IP and IM from spring to summer no comparable increase in PNR activities was observed. This could be the result of the fact that during summer the soil pH mainly in the rhizosphere decreased due to higher exudation rates of the plants, which might reduces the ratio of ammonia to ammonium in soil. As ammonia oxidizing microbes are only able to use ammonia for oxidation [37], this fact may explain that there is no direct link between PNR and amoA gene copies when data from the spring and summer sampling were compared. AOA harbor urease genes allowing AOA to use urea as alternative substrate, when ammonia is limited in soil [38]. The observed increase of AOA between spring and summer might be related to that urea utilization and a mixotrophic lifestyle of AOA. The relatively low PNR in summer (in relation to the abundance of AOA and AOB) might be also a fact of nitrification inhibitors which are excreted by many plants into soil to increase their competitiveness towards ammonia against microbes [39]. This ability of plants to reduce activity of AOA and AOB becomes even more important in those sites which have been extensively managed and nitrogen contents are very low. Obviously plants on these plots were indeed very competitive for ammonia uptake as not only PNR rates did not change between spring and summer but also the abundance of AOA and AOB did not change.
In nearly all samples a higher number of amoA AOA than amoA AOB copies (AOA/AOB ratios were between 2 and 16) was found. These results are in accordance with some other recent studies indicating dominance from AOA over AOB in soils [40,41,14]. Up to now it is still unclear if these higher abundance values of AOA might be related to a significant contribution to ammonification in soils [42,43,44,45,46]. However considering the correlation of different amoA genes and the nitrification potential (R (AOA) = 0.836, p (AOA) ,0.001 and R (AOB) = 0.680, p (AOB ) ,0.001) our dataset delivers one more hint that AOA plays an important role at least for ammonia oxidation in soil, although based on molecular data still the oxidation of hydroxylamine to nitrite has not been proven for AOA. Overall it seems that AOA and AOB are similarly affected by land use intensification and it has been crystallized that nitrogen availability is the major driver for performance and abundance of the nitrifying community.

Land Use Intensity and Denitrification
Microbial communities in soil are generally stimulated by plant growth especially in the rhizosphere by root exudates, resulting in an increase of anoxic habitats in the rhizosphere and the need for microbes to use alternative electron acceptors like nitrate inducing denitrification [47,48]. This was confirmed in our study by reduced DEA activities in soil samples derived from EP, where plant biomass was significantly lower. However in our study DEA was not stimulated in summer compared to spring, which is surprising on the first glance but might be a result of a low nitrate concentrations found mainly in soils from intensively managed plots IP and IM. Furthermore the differences in water content and the corresponding differences in the redox potential between spring and summer sampling might have influenced DEA activity independent from land use intensities.
To describe denitrifying communities, the two nitrite reductase genes (nirK and nirS) were quantified. Overall, many prokaryotes from soil are able to denitrify and the proportion of denitrifiers within the soil microbial community was considered between 10% and 60% of the total bacterial and archaeal microbiota reaching values between 10 5 and 10 9 copies g 21 [27,49,50,51,52,53]. The values found in our study were in a similar range. In all samples gene copy numbers for nirK genes dominated over nirS genes (nirS/ nirK ratio 0.002-0.02). This is in accordance with data from other studies, which indicate nirK phylotypes being more dominant in the rhizosphere or in soils characterized by an intensive root development [54], or high nutrient contents due to climatic conditions (e.g. freezing and thawing, [55]). The relatively high WEOC (water extractable organic carbon) values confirmed this hypothesis (Table 3). A strong time effect (p,0.0001) was found for both types of nitrite reductase genes in all soil samples, with values up to 20 times higher in summer, which correlated to the increased overall microbial biomass values obtained in soil samples from summer compared to spring (data not shown). The increase in abundance of the nitrite reductase over time is in contrast to the potential measurements of denitrification, but could be easily explained as all denitrifiers are facultative anaerobes and can shift their metabolism from respiration to denitrification based on the conditions present in soil. Therefore not surprisingly this missing link between denitrifying activity and genetic potential was described in some studies before [56,57].
However several studies have indicated that the plant species composition can act as a major driver for denitrifiers in soil [58] probably by providing different amounts and compositions of root exudates [59,60] mainly in summer which may explain the significant copy number increase of both genes in our study. Overall it can be concluded that abundance of denitrifiers is mainly driven by the amount and quality of carbon present in soil and consequently increase or decrease of denitrifiers is strongly linked to development of overall microbial biomass, whereas the induction of genes involved in denitrification is linked to the presence of nitrate and the redox conditions present in soil.

Land Use Intensity and Nitrogen Fixation
Surprisingly the highest amount of nifH copies was found at IP sites and not as we assumed at sites, which have been used extensively (EP) were nitrogen is limited. Despite no differences were found in the WEOC, this might be related to the fact that in soil the amount of energy equivalents needed for nitrogen fixation is too low in soil samples derived from EP [61], confirming earlier data obtained by [62] from glacier fore fields, where also low abundance of nifH as well as low potential activity for nitrogen fixation was found at the youngest sites, which had the lowest supply with nitrogen and lacked of energy equivalents. [63] as well as [64] showed in their study that intensive grazing has positive influences on nitrogen fixing communities in grassland soil, due to the high amounts of excrements which provide good sources for energy generation.
However it must be stated that the occurrence of legumes was not correlated to the nifH gene abundance, so we assume that the obtained results just reflect the occurrence of free living nitrogen fixing microbes, as we did not include nodule analysis in our study.

Conclusion
Obviously there is a strong link between land use intensity and microbial community structure in soil. Some of these effects might be direct effects and a response to e.g. increased fertilization regimes, other might be indirect effects and mediated by changes in the plant diversity and biomass (e.g. through changes in exudation). However the described data just reflect potentials for certain processes. If the same response pattern can be confirmed for the activation of certain functional traits remains an open question for future research focusing more on mRNA based studies with a more dense net of sampling time points.

Supporting Information
Table S1 List of plant species with species categorization according to ecological strategy type (Grime 1997) and Ellenberg values for moisture (M) and nutrients (N) (Ellenberg et. al. 2001). Different ecological strategy types as given in the table are competitors (c), competitive ruderals (cr), stress-tolerant competitors (cs), stress-tolerant ruderals (sr), csr-plants (competition is restricted by the combined effects of stress and disturbance), stresstolerants (s) and ruderals (r). Ellenberg indicator values are normalized values on an ordinal scale from low to high values (1 to 9) whereas 1 means high tand 9 low tolerance to nutrient-and moisture stress. NA shows plants which are not categoriesed by Ellenberg and x signs these with a large ecolocial amplitude. Occurence of each plant species on the different plots is given in percent of coverage. (DOC)