Next Article in Journal
Pollination Is Sufficient, Even with Low Bee Diversity, in Pumpkin and Winter Squash Fields
Next Article in Special Issue
Indicative Value of the Dominant Plant Species for a Rapid Evaluation of the Nutritional Value of Soils
Previous Article in Journal
Assessing Soil Organic Carbon in Soils to Enhance and Track Future Carbon Stocks
Previous Article in Special Issue
Soil Amendment with Biosolids and Inorganic Fertilizers: Effects on Biochemical Properties and Oxidative Stress in Basil (Ocimum basilicum L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating Heathland Restoration Belowground Using Different Quality Indices of Soil Chemical and Biological Properties

1
Department of Sustainable Land Management & Soil Research Centre, School of Agriculture, Policy and Development, University of Reading, Reading, Berkshire RG6 6AR, UK
2
Department of Protection of the Soil, Plant, Water System, Institute of Natural Resources and Agrobiology of Seville, 41012 Seville, Spain
3
Department of Life & Environmental Sciences, Faculty of Science & Technology, Bournemouth University, Dorset BH12 5BB, UK
4
Department of Geography and Environmental Science, University of Reading, Reading, Berkshire RG6 6DW, UK
5
School of Water, Energy and Environment, Cranfield University, Cranfield MK43 0AL, Bedfordshire, UK
6
Department of Geosciences and Natural Resource Management, University of Copenhagen, Rolighedsvej 23, DK-1958 Frederiksberg C, Denmark
*
Author to whom correspondence should be addressed.
Agronomy 2020, 10(8), 1140; https://doi.org/10.3390/agronomy10081140
Submission received: 25 June 2020 / Revised: 30 July 2020 / Accepted: 1 August 2020 / Published: 5 August 2020
(This article belongs to the Special Issue Soil Quality Evaluation Using Biological Properties)

Abstract

:
Reversion of agricultural land to heathland and acid grassland is a priority for the conservation of these rare habitats. Restoration processes require interventions to reverse the effects of fertilization and acidity amelioration undertaken during decades of agricultural production. Belowground assessments of restoration success are few, and we have examined the utility of soil indices as a rationalized tool for land managers and restoration practitioners to assess the efficacy of restoration practice. To achieve this, we assessed a large number of variables, many of which might be near redundant, that could be optimized for such indices. We used a 14-year field experiment contrasting acidified pasture (treated with elemental sulphur), control (untreated) pasture, and adjacent native heathland and acid grassland sites. Based on biotic and abiotic parameters, several ‘heathland restoration indices’ (resembling soil quality indices) were generated using a minimum dataset identified through principal component analysis and a linear scoring system. For comparison we also conducted alternative analyses of all parameters, using nonmetric multidimensional scaling plots and analyses of similarity (ANOSIM). Use of heathland restoration indices showed that elemental sulphur application had changed the soil chemical conditions, along with the vegetation assemblage, to be comparable to that of native acid grassland, but not the belowground biology. ANOSIM on full datasets confirmed this finding. An index based on key variables, rather than an analysis of all biotic and abiotic parameters, can be valuable to land managers and stakeholders in acid grassland and heathland restoration.

1. Introduction

Lowland heathland and acid grasslands are both listed as priority habitats in the European Commission Habitats Directive [1] and the UK Biodiversity Action Plan [2,3]. These landscapes are of particular importance as they now provide a last refuge for many rare and endangered fauna, but have seen rapid decline with agricultural intensification since the beginning of the 19th century. In England, for example, 72% of heathland was lost between 1830 and 1980 [4], largely as a result of improvement of these systems for agricultural production. The remaining ~58,000 ha of lowland heathland in the UK represents approximately 20% of the total global area of this habitat [5]. Therefore, particular interest has been paid for decades to methods that revert improved agricultural land back to heathland or acid grassland systems [2].
Soil conditions required for heathland or acid grassland ecosystems to establish include low pH and low nutrient availability [6,7,8]. Improved agricultural land has typically undergone a long-term regime of heavy fertilization and acidity amelioration in order to maximize yield for modern agricultural production [9]. These activities can cause long-lasting changes to soil chemical characteristics that leave soils unsuitable for the growth of heathland flora [10].
The re-establishment of heathland and acid grassland plant communities, therefore, needs to begin belowground through the modification of soil chemistry toward conditions found in existing heaths and grasslands [11]. One such restoration technique involves chemical amendment of improved agricultural land through the application of elemental sulphur as an acidifying agent [9,12]. The reduction of soil pH to <5 is usually a prerequisite to re-establishing acid grassland and heathland communities [10]. Therefore, successful establishment of heathland and acid grassland requires a reversal of the changes imposed on the soil by decades of fertilization and liming [13,14,15]. Consequently, heathlands and acid grasslands would be considered poor quality in the context of agricultural production. Therefore, restoration activities on agricultural land would require a perceived reduction in conventional soil quality (at least in terms of chemical fertility) presuming a good quality soil is intended for optimized agricultural productivity.
Statistically based soil quality indices, which use principal component analysis (PCA) to create a minimum dataset for linear scoring functions, are commonly used in the quantification of soil quality. Soil quality indices (SQIs), calculated in this manner, have been used to assess the effects of management on soil quality, such as land use change [16] and application of soil amendments [17,18], and agricultural systems, such as rice paddy [19] and vegetable production [20]. Variables that have been included in these indices include soil physico-chemical variables (e.g., pH, soil organic matter content, extractable nutrients, particle size distribution) and biological (e.g., microbial biomass, respiration, enzyme activity). Consolidation of multivariate datasets into a single index makes results of complicated datasets more accessible and easily communicated to practitioners and stakeholders. The use of statistical tools to generate a minimum dataset not only avoids data redundancy (i.e., lots of variables may be correlated or confound with one another, rendering the use of all of them unnecessary), but also reduces the opportunity for bias, compared to methods of soil quality assessment that use expert opinion [21].
Soil conditions required for a heathland ecosystem to establish include a low pH and low nutrient availability. Therefore, the variables that are considered to be ‘less is better’ or ‘more is better’ in an SQI scoring system to assess agricultural land would be reversed when applied to an SQI developed for heathland restoration. In response, we have replaced the term Soil Quality Index (SQI) with Heathland Restoration Index (HRI) to reflect reversal of soil conditions from those suitable for agricultural purposes toward those that are beneficial for heathland restoration (i.e., a soil with a high HRI may score poorly in an SQI for agricultural land).
Based on a field experiment on the Isle of Purbeck (Dorset, UK), a study developed an HRI to evaluate the efficacy of long-term chemical acidification treatments for the reversion of agricultural soil to heathland and acid grassland systems using a field experiment [22]. The work concluded that the application of elemental sulphur to acidify pasture was successful in creating an edaphic environment that was comparable with nearby native heathland [22]. A long-term impact on the plant community composition was also observed, with an increase in forbs and heather species in the acidified pasture, compared to the control [22,23]. As the restoration technique was the application of chemical amendments, only soil chemical parameters, such as extractable nutrients and pH, were used in the development of this HRI. An important limitation of an HRI based entirely on chemical characteristics is that it excludes consideration of the influence of, and interaction with, soil biota. Our work will build on this original assessment, on the same field site, to encompass biological variables, as well as chemical variables.
Soil biota, such as ‘decomposers’ and ‘nutrient transformers’ (e.g., microorganisms), ‘biocontrollers’ (e.g., nematodes), and ‘ecosystem engineers’ (e.g., earthworms), are an important indicator of heathland habitat quality [24,25]. Moreover, management activities and perturbations may promote a response in soil biological properties more rapidly than chemical and physical attributes [26]. Therefore, if the biological characteristics of the chemically treated soils were also to be considered when developing an HRI, it is likely that a different HRI would be observed. If so, an index based solely on chemical variables could be inadequate as a tool to provide an early indication of restoration success belowground. Alternatively, if the inclusion of biological variables in the development of an HRI has no impact on the overall structure of the resulting HRI, then this will also provide valuable insights concerning how the soil system should be assessed during reversion to heathland and acid grassland. Biological indicators often require more resources to measure than chemical indicators. So, if an HRI based on chemical variables alone is sufficient to inform restoration managers of the success of reversion to heathland, this outcome would have clear financial and logistical benefits. Such a finding would also result in an HRI that satisfies the criteria, laid out by Menta [27], that indicators of soil health or quality should: (1) Be sensitive to variations of soil management, (2) have good correlation with the beneficial soil functions, (3) be helpful in revealing ecosystem processes, (4) have comprehensibility and utility for land managers, and (5) be cheap and easy to measure.
The current paper extends the development of an HRI to incorporate biological parameters, including micro- and macrofauna abundance and microbial respiration, in order to assess soil quality in a heathland-restoration context. The objectives of this study were: (1) To establish whether the acidification treatments that alter soil chemistry also affect the aboveground (autotrophic) community, (2) to investigate whether soil chemistry (or an HRI based on chemical variables) also leads to changes in belowground communities (or an HRI based on biological variables), (3) to ascertain whether an HRI based on biological variables can distinguish between treatments more accurately than an HRI based only on chemical variables, and (4) to determine whether any HRI based on a minimum dataset and a linear scoring system identifies acidification treatment effects in the same manner as multivariate techniques such as PCA or analysis of similarity. Our results encompassed numerous chemical and biological variables, which we did not consider as univariate data as this was beyond the scope of our objectives. Here we will focus on the separation of acidifying treatments in terms of multivariate analyses and heathland restoration indices. Detailed discussion of individual chemical (e.g., pH, extractable elements) and biological variables (microbial community, decomposition, earthworm biomass, nematode abundance, etc.) can be found in Tibbett et al. [22]. Comparisons of different soil quality index evaluation methods, such as expert opinion versus statistically modeled minimum datasets and linear versus nonlinear scoring have taken place elsewhere [17,20,21,28]. However, our research was interested in the effects of the initial variables included in the PCA on the resultant index, with particular reference to biological versus chemical variables in a heathland restoration context.

2. Materials and Methods

2.1. Site Description

We compared four environments: (1) Chemically acidified pasture, (2) control pasture, (3) native heathland, and (4) native acid grassland. A detailed description of the field experiment utilized can be found in Diaz et al. [23], Green el al. [6], and Tibbett et al. [9,22]. Briefly, experimental plots treated with different acidifying agents were established in 1999 on the Isle of Purbeck, Dorset, UK (50.657425° N, −2.067085° W). The plots were on agricultural pasture of two contiguous farms, farmed as one consistent enterprise since the 1950s. During the 1950s and 1960s an aggressive policy of agricultural improvement of heathland through the addition of rock phosphate, manure, and chalk marl altered the soil conditions considerably. Compared to the surrounding (unimproved) acidic grassland and heath, soil pH had been increased from 3.6 to 6.6. In addition, rock phosphate had been applied over this period to elevate total phosphorus (Olsen P) levels by one-third, from 1000 mg/kg in the surrounding landscape to 1550 mg/kg in the experimental field (see Tibbett et al. [9] for more details). At the start of the experiment in the vegetation was an improved pasture in which Lolium perenne and Agrostis capillaris were the most constant and abundantly occurring grasses (see Diaz et al. [23] for a more detailed vegetation assessment).
The chemically acidified pasture treatment consisted of 10 replicate 50- × 50-m plots of improved agricultural pasture that were treated with 2000 kg ha−1 of elemental sulphur in May 2000 as an acidifying agent, and an additional 1600 kg ha−1 in 2001. Clippings of Calluna vulgaris were sown in the control and acidified pasture one year after the final application of elemental sulphur. The control pasture treatment consisted of 10 replicates of 50- × 50-m plots that received no chemical amendments. In addition, four replicate sites, representative of native heathland and native acid grassland, were selected in the nearby Middlebere Heath and Scotland Heath, respectively, in order to put the acidified pasture into context.

2.2. Soil Sampling

Fourteen years after initial application of elemental sulphur amendments, soil samples were collected using a gauge auger at 0–5-cm, 5–10-cm, and 10–15-cm depths (after removal of the litter layer). Twenty-five composite soil samples were taken from each plot following a standard ‘W’ sampling pattern, as outlined in ISO 18400-104 [29] and homogenized, resulting in three samples for each plot, one for each depth. All soils were sieved to 2 mm. The 0–5-cm depth layer was then split into two subsamples, the first for microbial was stored at 4 °C until analysis (described below). The second subsample was air-dried for chemical analysis. Chemical analysis was also conducted on the 2-mm sieved, air-dried, 5-10-cm, and 10-15-cm layer.

2.3. Vegetation Survey

Each plot was surveyed 15 years after initial application of elemental sulphur amendments on the control pasture and acidified pasture treatments. The percentage cover of all plant species was recorded in 10 randomly located 2-m × 2-m quadrats in each plot. Plant species were then classified as either: (1) Characteristic of acid/calcifugous grassland and heathlands (AH), (2) characteristic of mesotrophic grassland (M), or (3) intermediate species which can occur in both mesotrophic and acid grasslands (In). Classifications were according to the National Vegetation Classification [30,31].

2.4. Chemical Variables

Soil pH of 2.5:1 water-soil slurry was measured after shaking for 15 min at 120 rpm [32]. Exchangeable Al3+, Ca2+, Na+, Mg2+, K+ and Mn2+ and extractable Fe were determined in a 10:1 soil to 1 M ammonium nitrate (NH4NO3) extraction [33]. Soil available P was extracted with a 0.1 M H2SO4 solution [34]. Results and methods of chemical analysis are presented in Tibbett et al. [22].

2.5. Biological Variables

Biological variables encompassed both microbial analysis and soil fauna. Soil microbes react to changes in their environment by adjusting their (1) activity, (2) biomass, and (3) community structure [35]. Consequently, microbial analysis of the soil included microbial biomass determination of C, N, and P; microbial respiration; and microbial community composition. Microbial techniques were applied only to the top 0–5-cm soil depth.
Faunal groups selected as indicators needed to: (1) form a dominant group and occur in all soil types, (2) have a high abundance and high biodiversity, and (3) play an important role in soil functioning, e.g., in food webs [35]. These factors, combined with ease of collection and identification, led to our selection of the following soil fauna variables: Abundance and biomass of juvenile and adult epigeic, endogeic, and anecic earthworms; abundance and feeding group of nematodes; abundance of rotifers; and abundance of tardigrades. Methods for each of these variables are described below.

2.5.1. Soil Microbial Biomass

Microbial biomass was determined by the fumigation–extraction method. Two sets of field-fresh soils were bottled with a mass equivalent to 12.5 g for C and N extraction, and 5 g for P extraction of oven-dried sample. One set of the samples was fumigated with ethanol-free chloroform (CHCl3) for 24 h ± 1 h, for comparison.
Microbial biomass C and N were extracted in fumigated and nonfumigated samples with 50 mL of 0.5 M potassium sulphate (K2SO4) solution [36]. Microbial biomass P was extracted in fumigated and nonfumigated samples with 100 mL of sodium hydrogen carbonate (NaHCO3) reagent [37]. Briefly, samples were shaken for 30 min and then filtered and stored at −15 °C until analysis. The determination of microbial biomass C and N consisted of the removal of inorganic C in the form of CO2 through potassium persulphate (K2S2O8) acidification and ultraviolet (UV) light irradiation. The solution was measured at 550 nm in a spectrophotometer.
For the determination of microbial biomass P, 5 mL soil extract, 1 mL 1.5 M H2SO4, 20 mL 0.15% m/v ammonium molybdate reagent, and 5 mL ascorbic acid (C6H8O6) were mixed and left for 30 min for color development. The absorbance was measured at 880 nm in the spectrophotometer and microbial P, in mg L−1, determined from a standard graph.
Final microbial biomass C, N, and P were calculated with the conversion coefficients for carbon (KEC = 0.45), nitrogen (KEN = 0.45), and phosphorus (KEP = 0.40) [38,39,40,41].

2.5.2. Community-Level Physiological Profiles

Microbial respiration for different carbon sources was determined using the multiple substrate-induced respiration (MSIR) MicroResp™ method [42,43] as a convenient, sensitive, and rapid analysis for detection of community-level physiological profiles (CLPP). The MicroResp™ colorimetric method [44] is based on the measurement of the rate of CO2 generated by soil microorganisms in response to addition of different carbon substrates. The resultant profiles reflect the functional ability of the soil community to degrade the added substrates. After incubation of fresh soils (3–5 days at 25 °C), 0.25–0.48 g of soil was added to a 96-deep-well microtiter plate. Each soil was prepared in 32 wells (4 replicates).
The carbon sources selected were: Two amino acids, ɣ-Aminobutyric acid (GABA) and L-Alanine (ALA); one carbohydrate, α-D-Glucose (Glu); three carboxylic acids, L-Malic acid (MalA), Citric acid (CitA), and α-Ketoglutaric acid (aKG); and water. Carbon sources were added to soils in deep-well plates and then sealed with the detection plates containing cresol red.
The detection plates were read at 570-nm absorbance at the beginning of the incubation (time 0) and immediately after 6-h incubation (time 6). Carbon dioxide development, expressed as µg C g−1 h−1, was assessed using SoftMax Pro 6 Software.

2.5.3. Phospholipid Fatty Acid Analysis

Fatty acid profiles, indicative of microbial community phenotypic profile, were attained by phospholipid fatty acid (PLFA) analysis. The lipid extraction was accomplished by the modified Bligh and Dyer method [45] with a water replacement by citrate buffer [46] in a final 0.8:1:2 (v/v/v) citrate buffer:chloroform:methanol reagent solution. The reagent was added to 5 g of ground, freeze-dried soil. The resultant upper organic layer was mixed with 4 mL of chloroform and 4 mL of 0.15 M citrate buffer. Under a stream of N2 at <37 °C, the lower organic layer was dried. Solid phase extraction, or silicic acid column chromatography, was used as the lipid fractionation method [47]. Phospholipid fractions were eluted with 8 mL of methanol and dried under a stream of N2 at <37 °C. The phospholipid fraction was reconstituted and 0.2 M methanolic potassium hydroxide (KOH/MeOH) in order to hydrolyze the lipids. The solution was incubated for 30 min and 1 M acetic acid was added to stop the reaction. Finally, 5 mL of 4:1 (v/v) hexane:chloroform and 3 mL of deionized water were added. The top layer was filtered through sodium sulphate (Na2SO4) and rinsed with hexane. The solution was dried under a stream of N2 at room temperature. The dried samples were reconstituted with hexane, peak separated, and detected gas chromatography flame ionization detector (GC-FID) with an Agilent Technologies 6890N. The software used for analysis was Agilent G2070 ChemStation for GC systems and results were expressed as mol%.
Signature PLFAs were classified as biomarkers in five microbial groups: Gram-positive bacteria, gram-negative bacteria, arbuscular mycorrhizal fungi, ectomycorrhizal fungi, and actinobacteria. Twelve of the 33 fatty acids analyzed have been reported as biomarkers of specific microorganisms to differentiate the bacterial or fungal dominance in the soils (Table 1). The PLFA 18:1ω9c was not included in any group for the total biomarkers’ calculation as it has been found in both fungi [48,49,50] and bacteria [51,52,53].

2.5.4. Soil Fauna

Earthworms, nematodes, rotifers, and tardigrades were collected in each of the plots 16 years after initial application of elemental sulphur amendments. For each plot, a 20-cm × 20-cm × 20-cm soil block was removed using a flat shovel and placed in trays in the field for hand sorting. Earthworms were carefully removed, counted, and placed in a subsample of soil to be transported back to the lab for classification. Specimens were rinsed, blotted dry, and weighed individually (fresh weight). Earthworms were also recorded as juvenile or adult, and allocated to ecological groups: Epigeic, endogeic, or anecic, following Sherlock [66].
Soil from the block was homogenized prior to collecting ~100 g subsamples which were sieved to 3.35 mm and stored at 4 °C. Nematode, tardigrade, and rotifer extractions were conducted on these subsamples using a modified Baermann funnel technique, substituting extraction trays for funnels (samples collected after 24 and 72 h). Nematodes, rotifers, and tardigrades in the extracts were counted while alive on a Leitz Wilovert inverted microscope at 4x magnification. The two sampling times (after 24 and 72 h) were counted separately and then combined. Abundance was expressed as the number of individuals per 100 g soil dry weight equivalent. Gravimetric soil moisture content was determined by drying at 105 °C [32]. In addition, 100 nematodes per sample were also identified to feeding-group level (bacterial feeder, fungal feeder, plant parasite, omnivore, or predator) as described by Yeates et al. [67,68].

2.6. Heathland Restoration Index (HRI)

Five different HRIs were created in this study (Table 2): (1) An HRI based on chemical variables only (HRIchem), (2) an HRI based on biological variables only (HRIbio), and (3) a combined HRI encompassing both chemical and biological indicators (HRIcomb). Biological indicators were then partitioned to create a biological HRI for: (4) Microbes (HRImic) and (5) soil micro- and macrofauna (HRIfau) (Table 2).
Calculations of the HRIs were based on a soil quality index (SQI) described by Andrews et al. [20] and Tibbett et al. [22]. Each variable was initially subjected to Kruskal-Wallis one-way ANOVA on ranks. Only variables that showed statistically significant differences between treatments at p < 0.05 were then further analyzed by PCA. Variables were analyzed using PCA in order to generate a minimum dataset based on the so-called ‘highly weighted factors’. Highly weighted factors are defined as the variable with the highest eigenvector in each principal component (PC) that has an eigenvalue >1, plus any variables with eigenvectors within 10% of the highest eigenvector, or >0.4 in PCs with an eigenvalue >1. A Pearson’s product moment correlation coefficient (r2) was conducted pair-wise on variables highlighted to be highly weighted factors. When r2 > 0.60, only the variable with the highest eigenvector was retained for the minimum dataset; when r2 < 0.60, both variables were retained for the minimum dataset. Variables within the minimum dataset were then subjected to a linear scoring system based on whether ‘less is better’ (native heathland sites are lower than the control plots) or ‘more is better’ (native heathland sites are higher than the control plots). The higher the HRI, the closer the sample resembles a heathland environment. For ‘less is better’ variables, the lowest observed value was divided by each observation (e.g., the lowest observed value received a score of 1) and for ‘more is better’ indicators, each observation was divided by the highest observed value (e.g., the highest observed value received a score of 1).
The HRI was then calculated as follows:
H R I = i = 1 n W i S i
where S is the score of the indicator variable and W the weighted factor derived from the PCA:
W = V % C V %
where V is the proportion variance (%) explained by the principal component in which the variable has the highest eigenvalue (or is within 10% of the highest eigenvector, or above >0.4). CV is the cumulative variance (%) explained by all principal components with eigenvalue >1. The variables input into the PCA depended on the HRI in question (Table 2). The weighed factors, based on the PCA output, means that each HRI will have a different value for the highest possible HRI that can be assigned to a sample. We elected not to standardize the different HRIs to one another because, at this stage, we were interested in treatment effects within the different HRIs (chemical only, biological only, etc.).

2.7. Statistical Analysis

The effect of treatments (chemically acidified pasture, control pasture, native heathland, and native acid grassland) on a single characteristic, such as HRI or pH, was assessed with the use of a one-way ANOVA and Tukey’s post hoc testing. This analysis used Minitab Version 19, after being tested for homogeneity of variances with a Levene’s test. Multivariate analysis, including PCA, nonmetric multidimensional scaling (nMDS), cluster analysis, and analysis of similarity (ANOSIM), was conducted in Primer Version 6. An nMDS of plant community was conducted using a Bray Curtis resemblance matrix based on the square root of plant species’ abundance, with 1000 restarts and one-way ANOSIM, with treatment as a factor (1000 permutations). Soil chemical and biological variables were analyzed using PCA, cluster analysis, and ANOSIM (based on Euclidean distance resemblance matrices) to examine multivariate differences between treatments. Data were normalized prior to Euclidean distance resemblance matrices and were generated using Primer Version 6.

3. Results

3.1. Soil pH

Fourteen years after the first application of the acidifying agent, the acidified pasture maintained a significantly lower soil pH than the control pasture in all depth increments (Figure 1). This reduced soil pH in the acidified pasture was comparable with the neighboring acid grassland environment at all depths, and heathland from 5 cm downwards (Figure 1).

3.2. Plant Community

Nonmetric multidimensional scaling of the vegetation assemblage of the acidified pasture and the control pasture showed that acidified pasture resulted in a plant community that is significantly different to the control pasture (Figure 2). This difference is significant according to ANOSIM (p < 0.05). Vectors representing heathland and acid grassland species, such as Calluna vulgaris (labeled as AH 3 in Figure 2), Erica cinereal (AH 4), Erica tetralix (AH 5), and Agrostis capillaris (AH 1), resulted in acidified pasture sites starting to cluster above the control pasture, which is dominated by more mesotrophic grassland species.

3.3. Chemical Heathland Restoration Index (HRIchem)

The soil chemical conditions in the acidified pasture moved in a trajectory that resulted in an HRIchem that is significantly higher than the control pasture and statistically comparable to acid grassland (Figure 3a). The chemical environment in the acidified pasture has not, however, matched the native heathland in terms of the HRIchem.
Principal component analysis put forward seven variables for the minimum dataset of the HRIchem. These vectors are labelled in Figure 3b. In PC1, the acidified pasture showed a chemical composition closer to acid grassland and between the control pasture and heathland soils. In PC2, the acidified pasture deviated from the other treatments (Figure 3b).
The trends observed in the HRIchem were in concordance with ANOSIMchem (Table 3). The acidified pasture was significantly different from the control pasture and heathland, but comparable with the acid grassland (Table 3).

3.4. Microbial Heathland Restoration Index (HRImic)

Contrary to the HRIchem, the HRImic was not significantly different for the acidified pasture when compared to the control (Figure 4a). That said, native acid grassland was also comparable to the control pasture. In terms of HRImic, only the heathland sites had a significantly different HRImic than the other treatments (Figure 4a). The PCA revealed that, although the heathland is more dissimilar to the other treatments due to clear separation along PC1 (Figure 4b), the acidified pasture is separated from the control pasture on PC2, but it seems to have become more dissimilar to the acid grassland and heathland environment (Figure 4b). These observations are reiterated in the ANOSIMmic (Table 4), where there is a significant difference between the acidified pasture and all other treatments. The native heathland and the acid grassland environments are the only two treatments that are not significantly different to one another, according to ANOSIMmic.

3.5. Soil Fauna Heathland Restoration Index (HRIfau)

The acidified pasture had a greater HRIfau, compared to the control pasture, but not to a degree that resulted in a significant difference (Figure 5a). However, the acidified pasture HRIfau was not significantly different to the heathland, suggesting an upward trajectory of the acidified pasture treatment toward the native heathland environment, in terms of HRIfau (Figure 5a). Using multivariate analysis (Figure 5b), the acidified pasture is more similar to the heathland resulting in there being no significant difference between them, according to ANOSIMfau (Table 5). However, unlike HRIfau, there was no significant difference between the heathland and the control pasture according to ANOSIMfau.

3.6. Biological Heathland Restoration Index (HRIbio)

When microbial and soil faunal parameters were combined, the resulting HRIbio mirrored the trends seen in HRImic, i.e., there were no significant differences between the control pasture, acidified pasture, and acid grassland. However, the heathland HRIbio is significantly higher than all other treatments (Figure 6a). This is also apparent in the separation in the clustering of treatments in the PCA (Figure 6b) since the acidified pasture is dissimilar from the control pasture on PC2, resulting in a significant difference according to ANOSIMbio (Table 6). However, the acidified pasture treatments are still separated from the acid grassland or heathland treatments on PC1, reflected in Table 6 with acidified pasture being significantly different to all other treatments according to ANOSIMbio.

3.7. Combined Chemical and Biological Heathland Restoration Index (HRIcomb)

When all chemical, microbial, and soil fauna variables were included to generate a combined HRI (HRIcomb), the acidified pasture was not significantly different from the control pasture treatment (Figure 7a). However, it was a little higher than the control treatment, so that it was also not significantly different from the acid grassland environment, suggesting an upwards trajectory. Multivariate analysis (ANOSIMcomb, Table 7) indicates that the acidified pasture is significantly different from all other treatments. However, closer examination of Figure 7b indicates that, although the acidified pasture treatment is separated from the control pasture on PC2, it is still separated from the heathland and acid grassland environments on PC1.

3.8. Data Summary

Table 8 summarizes the significant differences observed between the control pasture and the acidified pasture, and the acidified pasture and heathland. According to ANOSIM, regardless of the variables used (chemical, microbial faunal, or a combination), the acidified pasture treatment was significantly different to the control pasture. Using an HRI, however, only resulted in a significant difference between the acidified pasture and the control pasture when the index was based on chemical variables alone (HRIchem). With the exception of soil fauna, there was a significant difference between the acidified pasture and the heathland in both HRI and ANOSIM (Table 8).

4. Discussion

4.1. Aboveground Response to Changes in Soil Chemical Properties

Application of elemental sulphur has resulted in the chemical conditions in the acidified pasture being significantly different to the control pasture, both in terms of multivariate analysis of similarity (ANOSIMchem) and Heathland Restoration Index (HRIchem). Soil pH and phosphorus have been found to drive aboveground floral community composition in heathland and grassland systems [69,70,71] and low pH and P concentrations are characteristic of the native heathland in this study [22].
Consequently, the plant community in the acidified pasture is starting to shift toward one containing more ericaceous species and other species typical of UK heathlands, such as Agrostis capillaris (common bent) [72]. As a result, the vegetational assemblage of the acidified pasture is significantly different to the control pasture according to ANOSIM. The influence of acidity on aboveground plant communities in heathlands has been well documented in other locations [71,73] and this study is in concordance with those findings.
For surficial soil pH (0–5 cm), the acidified pasture is equivalent to that observed in the acid grassland; however, the pH has not reduced enough to be comparable with the heathland. For soil chemistry, using ANOSIMchem or HRIchem, there is no significant difference between the acidified pasture and the acid grassland, but we do observe a significant difference between the acidified pasture and the heathland. It is clear that the elemental sulphur amendments have reduced the soil pH so that it now reflects an acid grassland environment. However, the soil pH of the acidified pasture, 14 years after application of elemental sulphur, is more than one pH unit higher than it was just three years after application, when it would have been comparable to nearby heathland [23]. This suggests that we have passed the peak of influence of elemental sulphur application and the acidification process is beginning to reverse in the acidified pasture.

4.2. Belowground Response to Changes in Soil Chemical Properties

Management activities and perturbations are thought to induce a response in soil biological properties more rapidly than a response in soil chemical or physical properties [26]. However, no HRI that included biological variables (HRImic, HRIfau, HRIbio, or HRIcomb) resulted in a significant difference between the acidified pasture and the control pasture. This finding suggests that either (1) belowground biology, in this system, is slower to respond to adjusted soil chemistry than the aboveground plant community, (2) belowground biology is less affected by changes in soil chemistry compared to aboveground vegetation, or (3) belowground biology rebounded more rapidly if the acidified pasture soils were beyond the peak effects of the elemental sulphur application (discussed above), i.e., biological parameters respond more rapidly to perturbation (i.e., acidifying elemental sulphur application), but also revert more rapidly after the perturbation.
Whether changes in plant community composition lead to changes in belowground soil microbial communities or whether the belowground microbial communities facilitate a change in vegetation is a quandary for restoration ecologists [74]. The data we present on the shifts in plant community composition and the HRImic suggest that the former is the case in heathland reversion from improved pasture. Vegetation structure was found by Vogels et al. [75] to play an important role in determining density and species richness of the aboveground faunal groups in heathlands; so, perhaps the same is true for belowground soil fauna and microorganisms. However, clippings of Calluna vulgaris were sown in the acidified pasture one year after the final application of elemental sulphur [23], facilitating aboveground community change. Early-stage presence, or inoculation, of heathland belowground communities can reinforce the development of a heathland system through the interactions they form [76]. No such introduction of microbial communities was undertaken in our study.
Return of species in restoration does not only depend on the habitat quality, e.g., the appropriate soil conditions, but also on the capacity of the species to recolonize the habitat [77]. It has been suggested that plant communities condition the soil through rhizodeposition and exudates, which modify the physical, chemical, and biological environment in soil, thereby impacting on the microbial community [78]. However, some studies had contrasting results and found that plants have no impact on shaping the microbial community in restoration projects [79,80]. Alteration in the aboveground community does not automatically initiate a response in the belowground community, and the presence of aboveground or belowground heathland species alone does not necessarily lead to restored plant–soil interactions [81,82]
Soil amended chemically with elemental sulphur, although similar in pH, will differ from a soil that has acidified naturally as a result of build-up of ericaceous plant litter. Tibbett et al. [22] observed that soil carbon in acidified pasture was not significantly different to the control pasture, and was much lower than carbon held in heathland soils (data taken the same year as our soil chemical data). This similarity between acidified pasture and control pasture would have implications for the soil biology and may provide a reason why differences in biologically based HRIs have not been observed between the acidified pasture and the control pasture.
Some studies have shown that only the presence of both above- (vegetation) and belowground (microbiota and fauna) heathland species leads to a fast assembly toward the target ecosystem [81,83]. Therefore, in terms of assessing the efficacy of application of elemental sulphur to acidify pasture, sole reliance on chemical indicators to quantify the success of a restoration may lead to misleading conclusions.

4.3. HRI Based on a Minimum Dataset and Linear Scoring System Compared to Analysis of Similarity

The use of an HRI based on a minimum dataset, as opposed to multivariate analysis of a large dataset, will have obvious advantages to practitioners. Specifically, resources required for analysis of soil samples will be reduced with a minimum dataset. The findings using HRIs and analysis of similarities based on belowground biology concurred, showing that acidified pasture was not yet comparable with heathland and acid grassland systems. However, the HRIs failed to reflect differences in belowground biology between the control and acidified pasture, whereas multivariate analysis of similarity identified there were significant differences. On closer inspection of the multivariate analysis, we noted that rather than being similar to an acid grassland or a heathland (which would have resulted in a similar HRI), the acidified pasture shifted toward a different state entirely in terms of belowground biology. From this we might ask: Is the multivariate dissimilarity between the control pasture and the acidified pasture of consequence to land managers? Aboveground ericaceous species are being supported in the acidified pasture; so, depending on the motivations of the restoration project, belowground biology may not concern land managers if heathland plant species are dominant. Therefore, in the interest of directness, simplicity, and practicality for land managers and stakeholders, the use of an index may be more appropriate than comprehensive multivariate outputs, and will at least bring belowground processes into consideration. Further research will be needed to explore the motivations of conservation projects, and whether plants growing on acidified pasture can support the rare and endangered aboveground fauna that are the likely drivers for conservation efforts. To scientists, however, the new state of acidified pasture will be of interest, particularly with regard to the sustainability of the heathland vegetation in the absence of a comparable heathland soil biological community. This will require a more detailed analysis of variables that are contributing to the differences between acidified pasture, heathland, and acid grassland systems, to allow a more functional understanding of the processes and mechanisms in play.

5. Conclusions

Lowland heathland and acid grasslands are both listed as priority habitats in conservation efforts in the UK. Reversion of agricultural land back to heathland and acid grassland environments requires intervention to reverse the effects of the heavy fertilization and acidity amelioration undertaken during agricultural conversion. Application of elemental sulphur to acidify pasture resulted in chemical conditions in the acidified pasture being comparable to acid grassland. This, in turn, impacted aboveground plant communities toward more ericaceous species and other species typical of UK acid grassland and heathlands. However, heathland restoration indices that included belowground biological variables were not significantly different between acidified pasture and control pasture. When assessing the efficacy of heathland restoration, biological variables provide different insights compared to abiotic variables into soil responses to exogenous chemical modification. Multivariate analysis of similarity and heathland restoration indices concurred, showing that belowground biology in acidified pasture was not yet comparable to that of heathland or acid grassland. We, therefore, propose that the use of a simplified heathland restoration index, rather than an analysis that includes all variables, can be valuable to stakeholders in heathland restoration.

Author Contributions

Conceptualization, M.T., A.D., and S.D.; methodology, all authors; formal analysis, all authors; resources, M.T. and K.R.-R.; data curation, M.T., A.D., S.D., and T.F.; writing—original draft preparation, S.D.; writing—review and editing, all authors; visualization, S.D.; supervision, M.T.; project administration, M.T.; funding acquisition, M.T. All authors have read and agreed to the published version of the manuscript

Funding

This project was part of the RECARE project (http://www.recare-project.eu/), which received funding from the European Union’s Seventh Framework Programme for research, technological development, and demonstration under grant agreement no. 603498.

Acknowledgments

The authors would like to thank Undergraduate Research Opportunities Programme (UROP) students Dewi Lloyd and Jhojan-Fernando Rios-Orjuela for assistance with field work activities and the National Trust for the use of the field site and support for this project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Commission. Council Directive 92/43/EEC of 21 May 1992 on the conservation of natural habitats and of wild fauna and flora. Off. J. Eur. Union 1992. [CrossRef]
  2. HMSO. Biodiveristy: The UK Action Plan; HMSO: London, UK, 1994. [Google Scholar]
  3. UKBAP. UK Biodiversity Action Plan Priority Habitat Descriptions 2008 (Updated 2011); Department of Environment: London, UK, 2011. [Google Scholar]
  4. Clarke, C.T. Role of soils in determining sites for lowland heathland reconstruction in England. Restor. Ecol. 1997, 5, 256–264. [Google Scholar] [CrossRef]
  5. Newton, A.C.; Stewart, G.B.; Myers, G.; Diaz, A.; Lake, S.; Bullock, J.M.; Pullin, A.S. Impacts of grazing on lowland heathland in north-west Europe. Biol. Conserv. 2009, 142, 935–947. [Google Scholar] [CrossRef]
  6. Green, I.; Stockdale, J.; Tibbett, M.; Diaz, A. Heathland restoration on former agricultural land: Effects of artificial acidification on the availability and uptake of toxic metal cations. Water Air. Soil Pollut. 2007, 178, 287–295. [Google Scholar] [CrossRef]
  7. Mitchell, R.J.; Auld, M.H.D.; Hughes, J.M.; Marrs, R.H. Estimates of nutrient removal during heathland restoration on successional sites in Dorset, southern England. Biol. Conserv. 2000, 95, 233–246. [Google Scholar] [CrossRef]
  8. Farrell, M.; Healey, J.R.; Godbold, D.L.; Nason, M.A.; Tandy, S.; Jones, D.L. Modification of Fertility of Soil Materials for Restoration of Acid Grassland Habitat. Restor. Ecol. 2011, 19, 509–519. [Google Scholar] [CrossRef]
  9. Tibbett, M.; Diaz, A. Are sulfurous soil amendments (S0, Fe(II)SO4, Fe(III)SO4) an effective tool in the restoration of heathland and acidic grassland after four decades of rock phosphate fertilization? Restor. Ecol. 2005, 13, 83–91. [Google Scholar] [CrossRef]
  10. Walker, K.J.; Pywell, R.F.; Warman, E.A.; Fowbert, J.A.; Bhogal, A.; Chambers, B.J. The importance of former land use in determining successful re-creation of lowland heath in southern England. Biol. Conserv. 2004, 116, 289–303. [Google Scholar] [CrossRef]
  11. Kleijn, D.; Bekker, R.M.; Bobbink, R.; De Graaf, M.C.C.; Roelofs, J.G.M. In search for key biogeochemical factors affecting plant species persistence in heathland and acidic grasslands: A comparison of common and rare species. J. Appl. Ecol. 2008, 45, 680–687. [Google Scholar] [CrossRef]
  12. Owen, K.M.; Marrs, R.H.; Snow, C.S.R.; Evans, C.E. Soil acidification—The use of sulphur and acidic plant materials to acidify arable soils for the recreation of heathland and acidic grassland at Minsmere, UK. Biol. Conserv. 1999, 87, 105–121. [Google Scholar] [CrossRef]
  13. Davis, J.; Lewis, S.; Putwain, P. The re-creation of dry heathland and habitat for a nationally threatened butterfly at Prees Heath Common Reserve, Shropshire. Asp. Appl. Biol. 2011, 108, 247–254. [Google Scholar] [CrossRef]
  14. Lawson, C.S.; Ford, M.A.; Mitchley, J.; Warren, J.M. The establishment of heathland vegetation on ex-arable land: The response of Calluna vulgaris to soil acidification. Biol. Conserv. 2004, 116, 409–416. [Google Scholar] [CrossRef]
  15. Pywell, R.F.; Webb, N.R.; Putwain, P.D. Soil fertility and its implications for the restoration of heathland on farmland in Southern Britain. Biol. Conserv. 1994, 70, 169–181. [Google Scholar] [CrossRef]
  16. Raiesi, F. A minimum data set and soil quality index to quantify the effect of land use conversion on soil quality and degradation in native rangelands of upland arid and semiarid regions. Ecol. Indic. 2017, 75, 307–320. [Google Scholar] [CrossRef]
  17. Masto, R.E.; Chhonkar, P.K.; Singh, D.; Patra, A.K. Alternative soil quality indices for evaluating the effect of intensive cropping, fertilisation and manuring for 31 years in the semi-arid soils of India. Environ. Monit. Assess. 2008, 136, 419–435. [Google Scholar] [CrossRef] [PubMed]
  18. Romaniuk, R.; Giuffre, L.; Romero, R. A soil quality index to evaluate the vermicompost amendments effects on soil properites. J. Environ. Prot. (Irvine, Calif) 2011, 2, 502–510. [Google Scholar] [CrossRef] [Green Version]
  19. Biswas, S.; Hazra, G.C.; Purakayastha, T.J.; Saha, N.; Mitran, T.; Roy, S.S.; Basak, N.; Mandal, B. Establishment of critical limits of indicators and indices of soil quality in rice-rice cropping systems under different soil orders. Geoderma 2017, 292, 34–48. [Google Scholar] [CrossRef]
  20. Andrews, S.S.; Karlen, D.L.; Mitchell, J.P. A comparison of soil quality indexing methods for vegetable production systems in Northern California. Agric. Ecosyst. Environ. 2002, 90, 25–45. [Google Scholar] [CrossRef]
  21. Mukherjee, A.; Lal, R. Comparison of soil quality index using three methods. PLoS ONE 2014, 9, e105981. [Google Scholar] [CrossRef] [Green Version]
  22. Tibbett, M.; Gil-Martínez, M.; Fraser, T.; Green, I.D.; Duddigan, S.; De Oliveira, V.H.; Raulund-Rasmussen, K.; Sizmur, T.; Diaz, A. Long-term acidification of pH neutral grasslands affects soil biodiversity, fertility and function in a heathland restoration. Catena 2019, 180, 401–415. [Google Scholar] [CrossRef]
  23. Diaz, A.; Green, I.; Tibbett, M. Re-creation of heathland on improved pasture using top soil removal and sulphur amendments: Edaphic drivers and impacts on ericoid mycorrhizas. Biol. Conserv. 2008, 141, 1628–1635. [Google Scholar] [CrossRef]
  24. Fagúndez, J. Heathlands confronting global change: Drivers of biodiversity loss from past to future scenarios. Ann. Bot. 2013, 111, 151–172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Pulleman, M.; Creamer, R.; Hamer, U.; Helder, J.; Pelosi, C.; Peres, G.; Rutgers, M. Soil biodiversity, biological indicators and soil ecosystem services-an overview of European approaches. Curr. Opin. Environ. Sustain. 2012, 4, 529–538. [Google Scholar] [CrossRef]
  26. Paz-Ferreiro, J.; Fu, S. Biological Indices for Soil Quality Evaluation: Perspectives and Limitations. Land Degrad. Dev. 2016, 27. [Google Scholar] [CrossRef]
  27. Menta, C. Soil Fauna Diversity—Function, Soil Degradation, Biological Indices, Soil Restoration. In Biodiversity Conservation and Utilization in A Diverse World; IntechOpen: London, UK, 2012; pp. 59–94. [Google Scholar]
  28. Guo, L.; Sun, Z.; Ouyang, Z.; Han, D.; Li, F. A comparison of soil quality evaluation methods for Fluvisol along the lower Yellow River. Catena 2017, 152, 135–143. [Google Scholar] [CrossRef]
  29. International Organization for Standardization (ISO). ISO18400-104:Soil Quality—Sampling—Part 104: Strategies; International Organization for Standardization: Geneva, Switzerland, 2018. [Google Scholar]
  30. Rodwell, J.S. (Ed.) Mires and Heaths. In British Plant Communities; Cambridge University Press: Cambridge, UK, 1998; Volume 2. [Google Scholar]
  31. Rodwell, J.S. (Ed.) Grasslands and Montane Communities. In British Plant Communities; Cambridge University Press: Cambridge, UK, 1998; Volume 3. [Google Scholar]
  32. Rowell, D.L. Soil Science: Methods & Applications; Taylor & Francis: London, UK, 1994. [Google Scholar]
  33. Stuanes, A.O.; Ogner, G.; Opem, M. Ammonium nitrate as extractant for soil exchangeable cations, exchangeable acidity and aluminum. Commun. Soil Sci. Plant Anal. 1984, 15, 773–778. [Google Scholar] [CrossRef]
  34. Sørensen, N.K.; Bülow-Olsen, A. Metode 14. Fosofortallet Pt. Plantedirektoratet. In Flles Arb. Jordbundsanalyser; Landbrugsministeriet: Lyngby, Denmark, 1994; pp. 1–4. [Google Scholar]
  35. Schloter, M.; Dilly, O.; Munch, J.C. Indicators for evaluating soil quality. Agric. Ecosyst. Environ. 2003, 98, 255–262. [Google Scholar] [CrossRef]
  36. Vance, E.D.; Brookes, P.C.; Jenkinson, D.S. An extraction method for measuring soil microbial biomass, C. Soil Biol. Biochem. 1987, 19, 703–707. [Google Scholar] [CrossRef]
  37. Brookes, P.C.; Powlson, D.S.; Jenkinson, D.S. Measurement of microbial biomass phosphorus in soil. Soil Biol. Biochem. 1982, 14, 319–329. [Google Scholar] [CrossRef]
  38. Hedley, M.J.; Stewart, J.W.B. Method to measure microbial phosphate in soils. Soil Biol. Biochem. 1982, 14, 377–385. [Google Scholar] [CrossRef]
  39. Wu, J.; Joergensen, R.G.; Pommerening, B.; Chaussod, R.; Brookes, P.C. Measurement of soil microbial biomass C by fumigation-extraction-an automated procedure. Soil Biol. Biochem. 1990, 22, 1167–1169. [Google Scholar] [CrossRef]
  40. Jenkinson, D.S.; Brookes, P.C.; Powlson, D.S. Measuring soil microbial biomass. Soil Biol. Biochem. 2004, 36, 5–7. [Google Scholar] [CrossRef]
  41. Griffiths, B.S.; Spilles, A.; Bonkowski, M. C:N:P stoichiometry and nutrient limitation of the soil microbial biomass in a grazed grassland site under experimental P limitation or excess. Ecol. Process 2012, 1, 1–11. [Google Scholar] [CrossRef] [Green Version]
  42. Chapman, S.J.; Campbell, C.D.; Artz, R.R.E. Assessing CLPPs using MicroResp™—A comparison with biolog and multi-SIR. J. Soils Sediments 2007, 7, 406–410. [Google Scholar] [CrossRef]
  43. Creamer, R.E.; Stone, D.; Berry, P.; Kuiper, I. Measuring respiration profiles of soil microbial communities across Europe using MicroResp™ method. Appl. Soil Ecol. 2016, 97, 36–43. [Google Scholar] [CrossRef]
  44. Campbell, C.D.; Chapman, S.J.; Cameron, C.M.; Davidson, M.S.; Potts, J.M. A rapid microtiter plate method to measure carbon dioxide evolved from carbon substrate amendments so as to determine the physiological profiles of soil microbial communities by using whole soil. Appl. Environ. Microbiol. 2003, 69, 3593–3599. [Google Scholar] [CrossRef] [Green Version]
  45. Bligh, E.G.; Dyer, W.J. A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 1959, 37, 911–917. [Google Scholar] [CrossRef] [Green Version]
  46. Frostegård, A.; Bååth, E. The use of phospholipid fatty acid analysis to estimate bacterial and fungal biomass in soil. Biol. Fertil. Soils 1996, 22, 59–65. [Google Scholar] [CrossRef]
  47. Tunlid, A.; White, D.C. Biochemical analysis of biomass, community structure, nutritional status, and metabolic activity of microbial communities in soil. In Soil Biochemistry; Stotzky, G., Bollag, J.M., Eds.; Marcel Dekker Inc.: New York, NY, USA, 1992; Volume 7, pp. 229–262. [Google Scholar]
  48. Bååth, E. The use of neutral lipid fatty acids to indicate the physiological conditions of soil fungi. Microb. Ecol. 2003, 45, 373–383. [Google Scholar] [CrossRef]
  49. Bååth, E.; Anderson, T.H. Comparison of soil fungal/bacterial ratios in a pH gradient using physiological and PLFA-based techniques. Soil Biol. Biochem. 2003, 35, 955–963. [Google Scholar] [CrossRef]
  50. Kaur, A.; Chaudhary, A.; Kaur, A.; Choudhary, R.; Kaushik, R. Phospholipid fatty acid—A bioindicator of environment monitoring and assessment in soil ecosystem. Curr. Sci. 2005, 89, 1103–1112. [Google Scholar]
  51. Zelles, L. Identification of single cultured micro-organisms based on their whole-community fatty acid profiles, using an extended extraction procedure. Chemosphere 1999, 39, 665–682. [Google Scholar] [CrossRef]
  52. De Deyn, G.B.; Quirk, H.; Oakley, S.; Ostle, N.; Bardgett, R.D. Rapid transfer of photosynthetic carbon through the plant-soil system in differently managed species-rich grasslands. Biogeosciences 2011, 8, 1131–1139. [Google Scholar] [CrossRef] [Green Version]
  53. Pankhurst, L.J.; Whitby, C.; Pawlett, M.; Larcombe, L.D.; McKew, B.; Deacon, L.J.; Morgan, S.L.; Villa, R.; Drew, G.H.; Tyrrel, S.; et al. Temporal and spatial changes in the microbial bioaerosol communities in green-waste composting. FEMS Microbiol. Ecol. 2012, 79, 229–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Harwood, J.L.; Russell, N.J. Lipids in Plants and Microbes; George Allen & Unwin Ltd.: London, UK, 1984. [Google Scholar]
  55. Zogg, G.P.; Zak, D.B.; Ringelberg, N.W.; MacDonald, N.W.; Pregitzer, K.S.; White, D.C. Compositional and functional shifts in microbial communities due to soil warming. Soil Sci. Soc. Am. J. 1997, 61, 475–481. [Google Scholar] [CrossRef] [Green Version]
  56. Parkes, R.J.; Taylor, J. The relationship between fatty acid distributions and bacterial respiratory types in contemporary marine sediments. Estuar. Coast. Shelf. Sci. 1983, 16. [Google Scholar] [CrossRef]
  57. O’Leary, W.M.; Wilkinson, S.G. Gram-positive bacteria. In Microbial Lipids; Ratledge, C., Wilkinson, S.G., Eds.; Academic Press Ltd.: London, UK, 1988; Volume 1, pp. 117–185. [Google Scholar]
  58. Wilkinson, S.G. Gram-negative bacteria. In Microbial Lipids; Ratledge, C., Wilkinson, S.G., Eds.; Academic Press Ltd.: London, UK, 1988; Volume 1, pp. 299–457. [Google Scholar]
  59. Olsson, P.A. Signature fatty acids provide tools for determination of the distribution and interaction of mycorrhizal fungi. FEMS Microbiol. Ecol. 1999, 29, 303–310. [Google Scholar] [CrossRef]
  60. Federle, T.W. Microbial distribution in soil—New techniques. In Perspectives in Microbial Ecology; Megusar, F., Gantar, M., Eds.; Slovene Society for Microbiology: Ljubljana, Slovenia, 1986; pp. 493–498. [Google Scholar]
  61. Vestal, J.R.; White, D.C. Lipid analysis in microbial ecology: Quantitative approaches to the study of microbial communities. Biogeoscience 1989, 39, 535–541. [Google Scholar] [CrossRef]
  62. Frostegard, A.; Tunlid, A.; Baath, E. Phospholipid fatty acid composition, biomass, and activity of microbial communities from two soil types experimentally exposed to different heavy metals. Appl. Environ. Microbiol. 1993, 59, 3605–3617. [Google Scholar] [CrossRef] [Green Version]
  63. Pankhurst, C.E.; Yu, S.; Hawke, B.G.; Harch, B.D. Capacity of fatty acid profiles and substrate utilization patterns to describe differences in soil microbial communities associated with increased salinity or alkalinity at three locations in South Australia. Biol. Fertil. Soils 2001, 33, 204–217. [Google Scholar] [CrossRef]
  64. Lechevalier, M.P. Lipids in bacterial taxonomy: A taxonomists view. Crit. Rev. Microbiol. 1977, 5, 109–210. [Google Scholar] [CrossRef] [PubMed]
  65. Kroppenstedt, R.M. Fatty acid and menaquinone analysis of actinomycetes and related organisms. In Chemical Methods Bacteria System; Goodfellow, M., Minnikin, D.E., Eds.; Elsevier Science & Technology Books: London, UK, 1985; pp. 173–199. [Google Scholar]
  66. Sherlock, E. Key to the Earthworms of Britain and Ireland; Field Studies Council: Shrewsbury, UK, 2012. [Google Scholar]
  67. Yeates, G.W. Nematodes as soil indicators: Functional and biodiversity aspects. Biol. Fertil. Soils 2003, 37, 199–210. [Google Scholar] [CrossRef]
  68. Yeates, G.W.; Bongers, T.; De Goede, R.G.M.; Freckman, D.W.; Georgieva, S.S. Feeding Habits in Soil Nematode Families and Genera-An Outline for Soil Ecologists. J. Nematol. 1993, 25, 315–331. [Google Scholar] [PubMed]
  69. Riesch, F.; Stroh, H.G.; Tonn, B.; Isselstein, J. Soil pH and phosphorus drive species composition and richness in semi-natural heathlands and grasslands unaffected by twentieth-century agricultural intensification. Plant Ecol. Divers. 2018, 11, 239–253. [Google Scholar] [CrossRef]
  70. Manning, P.; Putwain, P.D.; Webb, N.R. The role of soil phosphorus sorption characteristics in the functioning and stability of lowland heath ecosystems. Biogeochemistry 2006, 81, 205–217. [Google Scholar] [CrossRef]
  71. Green, I.; Evans, D.; Diaz, A. Modifying soil chemistry to enhance heathland recreation: A use for sulphur captured during oil refining. Int. J. Plant Soil Sci. 2015, 6, 272–282. [Google Scholar] [CrossRef] [Green Version]
  72. Owen, K.M.; Marrs, R.H. The use of mixtures of sulfur and bracken litter to reduce pH of former arable soils and control ruderal species. Restor. Ecol. 2001, 9, 397–409. [Google Scholar] [CrossRef]
  73. Roem, W.J.; Klees, H.; Berendse, F. Effects of nutrient addition and acidification on plant species diversity and seed germination in heathland. J. Appl. Ecol. 2002, 39, 937–948. [Google Scholar] [CrossRef]
  74. Harris, J. Soil microbial communities and restoration ecology: Facilitators or followers? Science 2009, 325, 573–574. [Google Scholar] [CrossRef]
  75. Vogels, J.J.; Verberk, W.C.E.P.; Lamers, L.P.M.; Siepel, H. Can changes in soil biochemistry and plant stoichiometry explain loss of animal diversity of heathlands? Biol. Conserv. 2017, 212, 432–447. [Google Scholar] [CrossRef]
  76. Radujković, D.; van Diggelen, R.; Bobbink, R.; Weijters, M.; Harris, J.; Pawlett, M.; Vicca, S.; Verbruggen, E. Initial soil community drives heathland fungal community trajectory over multiple years through altered plant–soil interactions. New Phytol. 2020, 225, 2140–2151. [Google Scholar] [CrossRef] [PubMed]
  77. De Graaf, M.C.C.; Verbeek, P.J.M.; Bobbink, R.; Roelofs, J.G.M. Restoration of species-rich dry heaths: The importance of appropriate soil conditions. Acta. Bot. Neerl. 1998, 47, 89–111. [Google Scholar]
  78. Hart, S.C.; DeLuca, T.H.; Newman, G.S.; MacKenzie, M.D.; Boyle, S.I. Post-fire vegetative dynamics as drivers of microbial community structure and function in forest soils. For. Ecol. Manag. 2005, 220, 166–184. [Google Scholar] [CrossRef]
  79. Orozco-Aceves, M.; Standish, R.J.; Tibbett, M. Soil conditioning and plant-soil feedbacks in a modified forest ecosystem are soil-context dependent. Plant Soil 2015, 390, 183–194. [Google Scholar] [CrossRef]
  80. Orozco-Aceves, M.; Standish, R.J.; Tibbett, M. Long-term conditioning of soil by plantation eucalypts and pines does not affect growth of the native jarrah tree. For. Ecol. Manag. 2015, 338, 92–99. [Google Scholar] [CrossRef]
  81. Van der Bij, A.U.; Weijters, M.J.; Bobbink, R.; Harris, J.A.; Pawlett, M.; Ritz, K. Facilitating ecosystem assembly: Plant-soil interactions as a restoration tool. Biol. Conserv. 2018, 220, 272–279. [Google Scholar] [CrossRef]
  82. Kardol, P.; Wardle, D.A. How understanding aboveground-belowground linkages can assist restoration ecology. Trends Ecol. Evol. 2010, 25, 670–679. [Google Scholar] [CrossRef]
  83. Lane, M.; Hanley, M.E.; Lunt, P.; Knight, M.E.; Braungardt, C.B.; Ellis, J.S. Chronosequence of former kaolinite open cast mines suggests active intervention is required for the restoration of Atlantic heathland. Restor. Ecol. 2020, 28, 661–667. [Google Scholar] [CrossRef]
Figure 1. Soil pH at 0–5-cm, 5–10-cm, and 10–15-cm depths. Measured in 2.5:1 water–soil slurry. Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different, within each depth, according to one-way ANOVA and Tukey’s post hoc testing (p > 0.05); n = 10 for acidified pasture and control pasture; n = 4 for heathland and acid grassland.
Figure 1. Soil pH at 0–5-cm, 5–10-cm, and 10–15-cm depths. Measured in 2.5:1 water–soil slurry. Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different, within each depth, according to one-way ANOVA and Tukey’s post hoc testing (p > 0.05); n = 10 for acidified pasture and control pasture; n = 4 for heathland and acid grassland.
Agronomy 10 01140 g001
Figure 2. Nonmetric multidimensional scaling of plant community on acidified and control pasture. Based on Bray-Curtis resemblance matrix, dashed-line clusters encircle plots ≥ 50% similarity. Species vectors’ key: AH 1—Agrostis capillaris; AH 2—Anthoxanthum odoratum; AH 3—Calluna vulgaris; AH 4—Erica cinereal; AH 5—Erica tetralix; AH 6—Molinia caerulea; AH 7—Ulex europaeus; In 1—Cynosurus cristatus; In 2—Holcus lanatus; In 3—Juncus articulates; In 4—Juncus effuses; In 5—Lotus corniculatus; In 6—Luzula campestris/multiflora; In 7—Rumex acetosella; M 1—Achillea millefolium; M 2—Agrostis stolonifera; M 3—Bellis perennis; M 4—Centaurium erythraea; M 5—Cerastium fontanum; M 6—Cirsium arvense; M 7—Cirsium palustre; M 8—Crepis capillaris; M 9—Elymus repens; M 10—Festuca rubra; M 11—Hypochoeris radicata; M 12—Leontodon autumnalis; M 13—Lolium perenne; M 14—Plantago lanceolate; M 15—Ranunculus repens; M 16—Rubus fructicosus; M 17—Senecio jacobaea; M 18—Taraxacum officinale agg.; M 19—Trifolium dubium; M 20—Trifolium pratense; M 21—Trifolium repens; M22—Vicia sativa. Groupings of species (AH, In, or M) according to National Vegetation Classification [30,31].
Figure 2. Nonmetric multidimensional scaling of plant community on acidified and control pasture. Based on Bray-Curtis resemblance matrix, dashed-line clusters encircle plots ≥ 50% similarity. Species vectors’ key: AH 1—Agrostis capillaris; AH 2—Anthoxanthum odoratum; AH 3—Calluna vulgaris; AH 4—Erica cinereal; AH 5—Erica tetralix; AH 6—Molinia caerulea; AH 7—Ulex europaeus; In 1—Cynosurus cristatus; In 2—Holcus lanatus; In 3—Juncus articulates; In 4—Juncus effuses; In 5—Lotus corniculatus; In 6—Luzula campestris/multiflora; In 7—Rumex acetosella; M 1—Achillea millefolium; M 2—Agrostis stolonifera; M 3—Bellis perennis; M 4—Centaurium erythraea; M 5—Cerastium fontanum; M 6—Cirsium arvense; M 7—Cirsium palustre; M 8—Crepis capillaris; M 9—Elymus repens; M 10—Festuca rubra; M 11—Hypochoeris radicata; M 12—Leontodon autumnalis; M 13—Lolium perenne; M 14—Plantago lanceolate; M 15—Ranunculus repens; M 16—Rubus fructicosus; M 17—Senecio jacobaea; M 18—Taraxacum officinale agg.; M 19—Trifolium dubium; M 20—Trifolium pratense; M 21—Trifolium repens; M22—Vicia sativa. Groupings of species (AH, In, or M) according to National Vegetation Classification [30,31].
Agronomy 10 01140 g002
Figure 3. Chemical variables (a) Heathland Restoration Index (HRIchem). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIchem was 1.92 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIchem. Vectors of minimum dataset variables used in the HRIchem labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤4.5.
Figure 3. Chemical variables (a) Heathland Restoration Index (HRIchem). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIchem was 1.92 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIchem. Vectors of minimum dataset variables used in the HRIchem labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤4.5.
Agronomy 10 01140 g003
Figure 4. Microbial variables (a) Heathland Restoration Index (HRImic). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRImic is 1.24 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRImic. Vectors of minimum dataset variables used in the HRImic labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤4.
Figure 4. Microbial variables (a) Heathland Restoration Index (HRImic). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRImic is 1.24 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRImic. Vectors of minimum dataset variables used in the HRImic labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤4.
Agronomy 10 01140 g004
Figure 5. Soil fauna variables (a) Heathland Restoration Index (HRIfau). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIfau is 1.00 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIfau. Vectors of minimum dataset variables used in the HRIfau labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤2.
Figure 5. Soil fauna variables (a) Heathland Restoration Index (HRIfau). Error bars represent standard error and the same lowercase letters in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIfau is 1.00 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIfau. Vectors of minimum dataset variables used in the HRIfau labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤2.
Agronomy 10 01140 g005
Figure 6. Biological variables, microbial, and soil fauna variables combined (a) Heathland Restoration Index (HRIbio). Error bars represent standard error and the same lowercase letter in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIbio is 1.63 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIbio. Vectors of minimum dataset variables used in the HRIbio labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤ 6.
Figure 6. Biological variables, microbial, and soil fauna variables combined (a) Heathland Restoration Index (HRIbio). Error bars represent standard error and the same lowercase letter in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIbio is 1.63 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIbio. Vectors of minimum dataset variables used in the HRIbio labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤ 6.
Agronomy 10 01140 g006
Figure 7. All available variables combined (a) Heathland Restoration Index (HRIcomb). Error bars represent standard error and the same lowercase letter in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIcomb is 3.02 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIcomb. Vectors of minimum dataset variables used in the HRIcomb labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤ 6.
Figure 7. All available variables combined (a) Heathland Restoration Index (HRIcomb). Error bars represent standard error and the same lowercase letter in bars denote treatments that are not significantly different (p > 0.05) according to one-way ANOVA and Tukey’s post hoc testing. The highest possible value for HRIcomb is 3.02 (i.e., if a sample scored 1 for every variable). (b) Principal component analysis for HRIcomb. Vectors of minimum dataset variables used in the HRIcomb labeled. See Table 2 for list of variables used. Dashed-line clusters surround samples with a Euclidean distance ≤ 6.
Agronomy 10 01140 g007
Table 1. Phospholipid fatty acids (PLFAs) or biomarkers in specific microbial groups.
Table 1. Phospholipid fatty acids (PLFAs) or biomarkers in specific microbial groups.
Microbial GroupSignature PLFAs or BiomarkersReferences
Total bacteria15:0, 15:0i, 15:0ai, 16:0i, 17:0i, 17:0ai, 16:1 ω7c, 17:0c, 19:0c[47,54,55]
Gram-positive bacteria15:0i, 15:0ai, 16:0i, 17:0i, 17:0ai[56,57]
Gram-negative bacteria16:1ω7c, 17:0c, 19:0c[51,58]
Arbuscular mycorrhizal fungi16:1ω5[59]
Ectomycorrhizal fungi18:2ω6,9[46,59,60,61,62,63]
Actinobacteria18:0 (10Me)[64,65]
Table 2. Variables included in principal component analysis in order to generate a minimum dataset for five different variations of Heathland Restoration Index: chemical variables only (HRIchem); soil micro- and macrofauna variables only (HRIfau); soil microbial variables only (HRImic); soil biological variables (microbes an fauna) combined (HRIbio); and all variables (chemical and biological) combined (HRIcomb).
Table 2. Variables included in principal component analysis in order to generate a minimum dataset for five different variations of Heathland Restoration Index: chemical variables only (HRIchem); soil micro- and macrofauna variables only (HRIfau); soil microbial variables only (HRImic); soil biological variables (microbes an fauna) combined (HRIbio); and all variables (chemical and biological) combined (HRIcomb).
Chemical HRI (HRIchem)Soil Micro- and Macrofauna HRI (HRIfau)Soil Microbes HRI (HRImic)Biological HRI (HRIbio)Combined HRI (HRIcomb)
0-5 cm: Al; Ca; Fe; K;Earthworm Abundance and Mass:Total Bacteria AbundanceAll variables listed in HRIfau and HRImicAll variables listed in HRIchem and HRIbio
Mg; Mn; Na;Total (juvenile and adult; all ecological groups)G+:G-togethertogether
P; pHAdult (all ecological groups)Total Actinobacteria Abundance
Juvenile (all ecological groups)Total Fungi Abundance34 Variables58 variables
5-10 cm: Al; Ca; Fe; K;Total epigeic (juvenile and adult)Total Arbuscular Mycorrhizal Fungi Abundance
Mg; Mn; Na;Adult epigeicBacteria:Fungi7 Variables in minimum dataset:15 variables in the minimum dataset
P; pHJuvenile epigeicMicrobial Biomass C, N and P
Total endogeic (juvenile and adult) Gram positive bacteria: Gram negative bacteria ratio (G+:G-)G+:G
10-15 cm: Al; Ca; Fe;Adult endogeicMicrobial Respiration (MicroResp™)Total Actinobacteria abundanceTotal Actinobacteria abundance
K; Mg; Mn;Juvenile endogeicα-Ketoglutaric acidTotal Fungi AbundanceBacteria:Fungi
Na; P; pHTotal anecic (juvenile and adult)L-AlanineTotal Arbuscular Mycorrhizal FungiMicroresp α-D-Glucose
Adult anecicCitric acidMicroresp Citric AcidMicrobial Biomass N
Juvenile anecicɣ-Aminobutyric acidJuvenile Endogeic AbundanceJuvenile Endogeic Mass
24 VariablesNematode Abundance:α-D-GlucoseTotal Tardigrade AbundanceTotal Tardigrade Abundance
Total (all functional groups)L-Malic acid 0–5 cm pH
7 * Variables inPlant ParasiteWater 0–5 cm Ca
minimum datasetBacterial Feeder 0–5 cm Fe
Fungal Feeder14 Variables 5–10 cm Mg
Predator 5–10 cm P
Omnivore3 * Variables in minimum dataset 10–15 cm Al
Total Rotifer Abundance 10–15 cm Mn
Total Tardigrade Abundance 10–15 cm Na
24 Variables
2 * Variables in minimum dataset
* Variables that are part of the minimum dataset are in bold underline.
Table 3. Analysis of similarity matrix based on chemical variables (ANOSIMchem). Pair-wise tests of Euclidean distance of treatment resemblance based on chemical variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have a p-value <0.05 and are shown in bold.
Table 3. Analysis of similarity matrix based on chemical variables (ANOSIMchem). Pair-wise tests of Euclidean distance of treatment resemblance based on chemical variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have a p-value <0.05 and are shown in bold.
ANOSIM Pair-Wise Test p-Value
HeathlandAcid GrasslandAcidified Pasture
Heathland---
Acid Grassland0.029--
Acidified Pasture0.0040.378-
Control Pasture0.0020.0210.002
Table 4. Analysis of similarity matrix based on microbial variables (ANOSIMmic). Pair-wise tests of Euclidean distance of treatment resemblance based on microbial variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have a p-value < 0.05 and are shown in bold.
Table 4. Analysis of similarity matrix based on microbial variables (ANOSIMmic). Pair-wise tests of Euclidean distance of treatment resemblance based on microbial variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have a p-value < 0.05 and are shown in bold.
ANOSIM Pair-Wise Test p-Value
HeathlandAcid GrasslandAcidified Pasture
Heathland---
Acid Grassland0.057--
Acidified Pasture0.0030.013-
Control Pasture0.0010.0040.002
Table 5. Analysis of similarity matrix based on soil fauna variables (ANOSIMfau). Pair-wise tests of Euclidean distance of treatment resemblance based on soil fauna variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have p-value < 0.05 and are shown in bold.
Table 5. Analysis of similarity matrix based on soil fauna variables (ANOSIMfau). Pair-wise tests of Euclidean distance of treatment resemblance based on soil fauna variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have p-value < 0.05 and are shown in bold.
ANOSIM Pair-Wise Test p-Value
HeathlandAcid GrasslandAcidified Pasture
Heathland---
Acid Grassland0.086--
Acidified Pasture0.0680.016-
Control Pasture0.2430.2290.017
Table 6. Analysis of similarity matrix based on biological variables, microbial, and soil fauna variables combined (ANOSIMbio). Pair-wise tests of Euclidean distance of treatment resemblance based on biological variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have p-value < 0.05 and are shown in bold.
Table 6. Analysis of similarity matrix based on biological variables, microbial, and soil fauna variables combined (ANOSIMbio). Pair-wise tests of Euclidean distance of treatment resemblance based on biological variables (see Table 2 for list of variables used). Treatments that are significantly different from one another have p-value < 0.05 and are shown in bold.
ANOSIM Pair-Wise Test p-Value
HeathlandAcid GrasslandAcidified Pasture
Heathland---
Acid Grassland0.057--
Acidified Pasture0.0010.006-
Control Pasture0.0030.0840.001
Table 7. Analysis of similarity matrix based on all available variables combined (ANOSIMcomb). Pair-wise tests of Euclidean distance of treatment resemblance based on all variables (see Table 2 for list of variables used.). Treatments that are significantly different from one another have p-value <0.05 and are shown in bold.
Table 7. Analysis of similarity matrix based on all available variables combined (ANOSIMcomb). Pair-wise tests of Euclidean distance of treatment resemblance based on all variables (see Table 2 for list of variables used.). Treatments that are significantly different from one another have p-value <0.05 and are shown in bold.
ANOSIM Pair-Wise Test p-Value
HeathlandAcid GrasslandAcidified Pasture
Heathland---
Acid Grassland0.029--
Acidified Pasture0.0010.048-
Control Pasture0.0010.0600.001
Table 8. Summary of pair-wise comparisons of acidified pasture compared to control pasture and heathland, according to Heathland Restoration Index and analysis of similarity. * Heathland restoration indices marked with a Agronomy 10 01140 i001 resulted in a significant difference in treatments according to one-way ANOVA and Tukey’s post hoc testing (p < 0.05). See part a of Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 for full results. ** Analysis of similarity pair-wise comparisons marked with a Agronomy 10 01140 i002 resulted in a significant difference in treatments based on Euclidean matrix of similarity (p < 0.05).
Table 8. Summary of pair-wise comparisons of acidified pasture compared to control pasture and heathland, according to Heathland Restoration Index and analysis of similarity. * Heathland restoration indices marked with a Agronomy 10 01140 i001 resulted in a significant difference in treatments according to one-way ANOVA and Tukey’s post hoc testing (p < 0.05). See part a of Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 for full results. ** Analysis of similarity pair-wise comparisons marked with a Agronomy 10 01140 i002 resulted in a significant difference in treatments based on Euclidean matrix of similarity (p < 0.05).
Significant Difference in Pair-wise Comparison? ( Agronomy 10 01140 i003 if p < 0.05)
Heathland Restoration Index *Analysis of Similarity **
Acidified Pasture vs. Control PastureHRIchem Agronomy 10 01140 i004ANOSIMchem Agronomy 10 01140 i005
HRImic ANOSIMmic Agronomy 10 01140 i006
HRIfau ANOSIMfau Agronomy 10 01140 i007
HRIbio ANOSIMbio Agronomy 10 01140 i008
HRIcomb ANOSIMcomb Agronomy 10 01140 i009
Acidified Pasture vs. HeathlandHRIchem Agronomy 10 01140 i010ANOSIMchem Agronomy 10 01140 i011
HRImic Agronomy 10 01140 i012ANOSIMmic Agronomy 10 01140 i013
HRIfau ANOSIMfau
HRIbio Agronomy 10 01140 i014ANOSIMbio Agronomy 10 01140 i015
HRIcomb Agronomy 10 01140 i016ANOSIMcomb Agronomy 10 01140 i017

Share and Cite

MDPI and ACS Style

Duddigan, S.; Gil-Martínez, M.; Fraser, T.; Green, I.; Diaz, A.; Sizmur, T.; Pawlett, M.; Raulund-Rasmussen, K.; Tibbett, M. Evaluating Heathland Restoration Belowground Using Different Quality Indices of Soil Chemical and Biological Properties. Agronomy 2020, 10, 1140. https://doi.org/10.3390/agronomy10081140

AMA Style

Duddigan S, Gil-Martínez M, Fraser T, Green I, Diaz A, Sizmur T, Pawlett M, Raulund-Rasmussen K, Tibbett M. Evaluating Heathland Restoration Belowground Using Different Quality Indices of Soil Chemical and Biological Properties. Agronomy. 2020; 10(8):1140. https://doi.org/10.3390/agronomy10081140

Chicago/Turabian Style

Duddigan, Sarah, Marta Gil-Martínez, Tandra Fraser, Iain Green, Anita Diaz, Tom Sizmur, Mark Pawlett, Karsten Raulund-Rasmussen, and Mark Tibbett. 2020. "Evaluating Heathland Restoration Belowground Using Different Quality Indices of Soil Chemical and Biological Properties" Agronomy 10, no. 8: 1140. https://doi.org/10.3390/agronomy10081140

APA Style

Duddigan, S., Gil-Martínez, M., Fraser, T., Green, I., Diaz, A., Sizmur, T., Pawlett, M., Raulund-Rasmussen, K., & Tibbett, M. (2020). Evaluating Heathland Restoration Belowground Using Different Quality Indices of Soil Chemical and Biological Properties. Agronomy, 10(8), 1140. https://doi.org/10.3390/agronomy10081140

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop