Nitrogen deposition reduces plant diversity and alters ecosystem functioning: field-scale evidence from a nationwide survey of UK heathlands.

Findings from nitrogen (N) manipulation studies have provided strong evidence of the detrimental impacts of elevated N deposition on the structure and functioning of heathland ecosystems. Few studies, however, have sought to establish whether experimentally observed responses are also apparent under natural, field conditions. This paper presents the findings of a nationwide field-scale evaluation of British heathlands, across broad geographical, climatic and pollution gradients. Fifty two heathlands were selected across an N deposition gradient of 5.9 to 32.4 kg ha(-1) yr(-1). The diversity and abundance of higher and lower plants and a suite of biogeochemical measures were evaluated in relation to climate and N deposition indices. Plant species richness declined with increasing temperature and N deposition, and the abundance of nitrophilous species increased with increasing N. Relationships were broadly similar between upland and lowland sites, with the biggest reductions in species number associated with increasing N inputs at the low end of the deposition range. Both oxidised and reduced forms of N were associated with species declines, although reduced N appears to be a stronger driver of species loss at the functional group level. Plant and soil biochemical indices were related to temperature, rainfall and N deposition. Litter C:N ratios and enzyme (phenol-oxidase and phosphomonoesterase) activities had the strongest relationships with site N inputs and appear to represent reliable field indicators of N deposition. This study provides strong, field-scale evidence of links between N deposition--in both oxidised and reduced forms--and widespread changes in the composition, diversity and functioning of British heathlands. The similarity of relationships between upland and lowland environments, across broad spatial and climatic gradients, highlights the ubiquity of relationships with N, and suggests that N deposition is contributing to biodiversity loss and changes in ecosystem functioning across European heathlands.

Introduction

pan class="Species">Human activities associated with the production of energy, fertilisers and leguminous cropn>s have had a substantial effect on the global n>an class="Chemical">nitrogen (N) cycle [1]. Recent estimates suggest that anthropogenic contributions to the global reactive N pool have increased 10-fold since pre-industrial times and are predicted to double, relative to current day levels, by 2050 [2]. Deposition of N is spatially highly variable; in Europe, for example, rates currently range from ∼1 kg ha−1 yr−1 in the relatively pristine areas of northern Scandinavia to >50 kg ha−1 yr−1 in areas dominated by industry (e.g. Northern Italy) or intensive agripan class="Chemical">culture (e.g. Netherlands) [3].

The link between N deposition and changes in the structure and functioning of terrestrial ecosystems is now widely acknowledged, with N cited as one of the leading drivers of biodiversity loss at a global scale [4], [5]. pan class="Species">Manipn>ulation expn>eriments have demonstrated a wide range of plant respn>onses to N, including changes in plant phenology, physiology, biochemistry, productivity and assimilate allocation [6]–[8]. Resulting changes in compn>etition between neighbouring spn>ecies, as well as N-driven changes in soil pH and n>an class="Chemical">NH4 + concentrations and/or increased sensitivity to environmental stresses, are believed to be responsible for changes in plant community composition observed in many ecosystems [9]. Similarly, N-driven changes in the activity and composition of the soil microbial community [10], [11] have been linked with changes in key soil processes, such as rates of decomposition and nitrous oxide production [12]–[13], with implications for the functioning of affected ecosystems.

Although pan class="Species">many of the respn>onses repn>orted in N-pan class="Species">manipulation experiments are not specific to N and will be affected by factors such as climate and habitat management, several have been proposed as potential robust bioindicators of the effects and/or amounts of N deposition. These include foliar N accumulation, activities of enzymes involved in nutrient metabolism, soil C:N ratios and concentrations of N in soil solution and leachate [14], [15]. Relatively few studies have explicitly tested whether the potential bioindicators of (effects of) N deposition observed in manipulation experiments are apparent and/or widespread under natural, field conditions. Those studies which have done so have tended to focus on relationships between N deposition and either species richness or biochemical indicators, but rarely both. For example, Remke et al. [16] reported a decline in species richness of coastal dune systems with increasing N deposition, while in a recent analysis of UK Countryside Survey data, Maskell et al. [17] suggest that acidification and eutrophication associated with N deposition are responsible for biodiversity loss in UK grasslands and heathlands. Surveys focusing on biochemical indicators of N deposition identified in manipulation experiments (e.g. tissue N concentrations, soil nutrient availability, litter chemistry and/or enzyme activities) have provided evidence of field-scale changes in ecosystem functioning, particularly in relation to nutrient cycling, along gradients of ambient N deposition [18]–[20]. A small number of studies have demonstrated links between diversity loss, changes in ecosystem function and N deposition, including a regional survey of moorlands [18], and both a national [21] and pan-European [22] survey of acidic grasslands. To date, however, there has been no wide-scale evaluation of the relationships between community-level (compositional) bioindicators, biogeochemical (functional) indicators and N deposition in heathland ecosystems, despite the fact these have been shown to respond strongly to experimental N manipulation [23]–[26].

Atmospheric N ocpan class="Chemical">curs in both oxidised (NOy) and reduced forms (NHx); anthropn>ogenic emissions are associated principally with the combustion of fossil fuels (oxidised N) and intensive agripan class="Chemical">cultural practices (reduced N) processes [27]. Whilst ecosystem responses to increasing total N deposition have been extensively researched, far less is known about the sensitivity of soils and vegetation to the different forms of N [27]. Findings from the relatively few studies that have investigated the role of different forms of N have reported differences in the strengths of relationships between oxidised and reduced forms in terms of plant tissue chemistry [28] species richness and community composition [29], [30].

Upland and lowland heathlands are typically of low nutrient status and dominated by dwarf shrubs. However, their contrasting altitudes (lowlands typically <300m, uplands >300m [31]) result in characteristically different levels of soil organic matter, rates of nutrient cycling, plant species assemblages and pan class="Species">management regimes [32] which might be expn>ected to result in differing respn>onses to elevated N inpn>uts. Upn>land and lowland habitats also expn>erience different climates (typclass="Chemical">n>ically wetter and cooler in the former compared to the latter), which have the potential to modify their sensitivity to atmospn>heric N inpn>uts [33]. Indeed, although N pan class="Species">manipulation experiments have demonstrated similar qualitative responses to N in upland and lowland heathlands (e.g. increased Calluna productivity and reduced abundance of lower plants [24], [34], [35]), differences in terms of soil biogeochemical responses have been found [23], [26]. The ubiquity of heathland responses to N across a range of altitudinal and climatic gradients has not yet been tested. This study, therefore, presents the first comprehensive, nationwide survey of biological and biogeochemical indicators of N deposition across lowland and upland heathlands. It explicitly aims to evaluate relationships between N deposition and 1) plant community composition and species richness and 2) ecosystem functioning, as measured by soil nutrient stocks, availabilities and enzyme activities. Based on results from earlier manipulation experiments and surveys, we hypothesise that similar relationships with N deposition will be evident in upland and lowland heaths, and that reduced forms of N are having a greater impact than oxidised N on the diversity and functioning of UK heathlands.

Methods

Site selection and field surveying

Fifty two pan class="Species">Calluna vulgarisn> dominated heathland sites (representing building to mature phase stands, sensu Gimingham, [32]) were selected across a gradient of modelled N deposition (data provided by Ron Smith and Jane Hall, Centre for Ecology & Hydrology, Edinburgh, UK) spn>anning 5.9 kg N ha−1 yr−1 to 32.4 kg N ha−1 yr−1 (Figure 1). In addition, sites were chosen to be representative of the geographical, climatic and altitudinal (25 upland and 27 lowland heaths) variation encountered in heathland ecosystems across the United Kingdom.

Figure 1

Site survey locations in relation to total N deposition (2006).

(http://pollutantdeposition.defra.gov.uk/pollutant-maps).

Sites were visited during the summer of 2009 and grid reference and altitudinal details were recorded at each site. pan class="Species">Vascular and non-pan class="Species">vascular plant species were recorded and relative abundance was estimated by percentage cover in five randomly placed quadrats (2 m×2 m) at each site. Canopy height was measured to the nearest mm at the centre of each of the four 1 m2 quadrants comprising the 2 m×2 m quadrat, and averaged. 5–10 g of pan class="Species">Calluna shoots were collected into paper bags and air dried for later chemical analyses. Three soil cores (15 cm depth ×4 cm diameter) were extracted from the centre of the fifth quadrat for chemical analyses; these were stored in an ice-filled cool-box before transfer to controlled temperature rooms upon return to the laboratory, typically within 1–3 days. Five samples of fresh Calluna litter (5–10g) were also collected from each quadrat for enzyme assays. The litter was air dried and stored in paper bags.

Chemical analyses

To determine extractable soil pan class="Chemical">ammonium and pan class="Chemical">nitrate concentrations, 40 g of fresh soil was mixed with 150 ml of 1M pan class="Chemical">KCl and shaken for one hour on a rotary shaker before filtering. The filtrate was measured for inorganic N (NH4 + and NO3 −) using a SKALAR san flow++ segmented flow analyser (Skalar UK (Ltd), York, UK). Litter phosphomonoesterase (PME) activity was assayed using the artificial substrate p-nitrophenyl phosphate (p-NPP), in accordance with the method outlined by Johnson et al. [36]. Samples were centrifuged for 5 minutes at 3000 rpm and the resulting liquid was analysed colourimetrically (using a Perkin Elmer Lambda 3 UV spectrophotometer, Perkin Elmer, USA) at a wavelength of 460 nm. Paired non- reactive controls were prepared for each of the samples by adding the PNP and toluene solution immediately after the stopping chemicals to prevent chemical reaction. The readings were compared with pigment densities from a calibration curve obtained using six standards of p-nitrophenol prepared within the absorbance range of the samples.

Total N and pan class="Chemical">carbon (C) concentrations were determined using a FlashEA 1112 CN analyser (Thermo Scientific, UK). pan class="Species">Calluna shoot and litter samples were dried at 80°C for 24 hours before being ground into a fine powder with a ball mill for two minutes. 15–20 µg of milled material was weighed into tin capsules and combusted. Standard samples were prepared using certified reference materials. Soil pH (pan class="Chemical">H2O) was measured by mixing 10g field moist soil with 25ml deionised water, then allowing the mixture to stand for 30 minutes; pH was then recorded using a CORNING 220 pH meter (Corning Inc., USA).

Litter phenol oxidase activity was measured by vortex-mixing 1 g of ground pan class="Species">Calluna litter with 9 ml deionised class="Chemical">n>an class="Chemical">water. 3 ml of this solution was then added to 4.5 ml water and 7.5 ml of 1 M dihydroxyphenylalanine (DOPA) and then shaken for 9 minutes. After shaking, the mixture was centrifuged at 6000 rpm for 5 minutes. A sub-sample was then taken, filtered (0.2 µm) and absorbance measured at 460 nm using a Perkin Elmer Lambda 3 UV spectrophotometer, Perkin Elmer, USA. Soil moisture content was determined gravimetrically by weighing a 10g sample of fresh soil, drying this at 105°C for 24 hours, then re-weighing.

Data analyses

Relationships between response and explanatory variables (N deposition, rainfall and temperature) were analysed using the R statistical package, version 2.12.2 (R development core team, 2011). The climate data used were based upon UK 5 km2 gridded data sets from the UK Meteorological Office. Variables representing total annual precipitation and temperature were used, the latter represented by data on growing degree days (pan class="Chemical">GDD, sum of degree days above 5°C in a given year). Both precipn>itation and n>an class="Chemical">GDD were averaged over the period 1997–2006. Since total deposition is the product of oxidised and reduced N deposition these three explanatory variables could not be included in the same statistical model. In order to separate the combined effects of total N deposition from the individual influences of oxidised and reduced N forms, two separate models were, therefore, run for each response variable. Model 1  =  total N deposition, rainfall and temperature; Model 2  =  reduced N, oxidised N, rainfall and temperature. Percentage cover data for vegetation were arc-sine transformed prior to analysis. Log or square root transformations were used to transform other non-normally distributed data prior to modelling. Relationships were tested using multiple regression analysis or, for count data (e.g. species richness), generalised linear models (GLM) with poisson errors. Relationships with species richness were investigated firstly for all species recorded at a site and subsequently for the different plant functional groups present (lichens, mosses, grasses and forbs). Where curvature in the data was evident, models were fitted with a quadratic function. The most complicated models were fitted (to include interactions between explanatory variables, plus quadratic terms to test for non-linearity of responses) and then simplified using likelihood ratio deletion tests until the minimum adequate model was attained. Evidence of overdispersal was corrected by re-fitting the model with quasi-empirical structures, followed by F tests as outlined by Crawley [37]. Three analyses were performed for each parameter to encompass (a) all sites, (b) upland sites only and (c) lowland sites only.

Results

Higher plant responses

A significant decline in overall species richness was observed across all sites in relation to increasing levels of total N deposition (F1,50 = 27.5, P<0.001, Figure 2). The biggest losses were seen at the lower end of the N deposition gradient, with an average of 13 species lost per site between 5–10 kg N ha−1 yr−1, compared to an average of only 3 species lost as N increased from 10 to 20 kg N ha−1 yr−1. Beyond this, pan class="Disease">declines in spn>ecies richness plateaued, indicating a reduction in sensitivity as N loading increased. Tempn>erature was also a significant factor across all sites, both on its own and in interaction with N depn>osition (F1,50 = 4.4, P<0.05, Table 1: model 1). Spn>ecies richness was seen to decline by 50 % across all sites around a tempn>erature threshold of 1000 growing degree days (n>an class="Chemical">GDD) and total N inputs of 15 kg N ha−1 yr−1, and declined by a further 50 % when GDD exceeded 2000; this relationship was significant in the lowland-only dataset and as a not quite significant trend for the uplands.

Figure 2

Relationships between N deposition and plant species richness and abundance.

A decline of species richness across all plant groups correlates with increasing N deposition (A) and where N deposition and temperatures are higher (B). Whilst graminoid species richness declines in relation to increased N deposition inputs (C), graminoid abundance increases (D), indicating the promotion of fewer species at the higher end of the N gradient.

Table 1

Summary of optimal models for plant species richness in relation to N deposition and climate variables.

	Model 1 (Total N deposition + environmental variables)									Model 2 (oxidized/reduced N + environmental variables)
Response variable	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P
Total species	(−) N	51	***	(−) N	23	***	(−) N	26	**	(−) NOy	51	**	(−) NOy	24	***	(−) NOy	26	**
	(−) gdd		***				(−) gdd		**	(−) gdd		***	(−) NHx		*	(−) gdd		**
	N x gdd		*							(−) NHx		***
Total lichen species	(−) N	51	**	(−) N	24	**	-	-	-	(−) NHx	51	***	(−) NOy	24	**	(−) NHx	24	*
	N x gdd		*	N x gdd		**				(−) NOy		**	(−) NHx		*
										(−) gdd		***	NHx x gdd		**
													NHx x rain		*
Total bryophyte species	(−) gdd	51	***	-	-	-	(-) gdd	26	**	(-) gdd	51	***	-	-	-	(−) NOy	26	**
																(−) gdd		*
																(−) rain		*
Total graminoid species	(−) N	51	**	(−) N	24	*	(−) N	26	*	(−) NOy	51	*	(−) NHx	24	*	(−) NOy	26	*
	(−) gdd		**				(−) gdd		**	(−) NHx		*				(−) gdd		***
										(−) gdd		**				NHx x gdd		*
Total forb species	(−) N	51	***	(−) N	24	**	(−) N	26	**	(−) NOy	51	***	(−) NOy	24	***	(−) NOy	26	***
				N x rain		*				(−) NHx		*	(−) NHx		*	NOy x rain		*
										(−) gdd		*				NHx x rain		*

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; gdd  =  growing degree days.

A decline of species richness across all n class="Chemical">plant groups correlates with increasing N deposition (A) and where N deposition and temperatures are higher (B). Whilst pan class="Chemical">graminoid spn>ecies n>an class="Disease">richness declines in relation to increased N deposition inputs (C), graminoid abundance increases (D), indicating the promotion of fewer species at the higher end of the N gradient.

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; pan class="Chemical">gdd  =  growing degree days.

The number of forb species also decreased with increasing total n class="Chemical">N deposition (F1,50 = 21.2, P<0.001, Table 1: model 1; Figure 2b), with the steepest rate of decline seen with modest increases of N at the lower end of the deposition gradient (4 species lost on average per site between 5–10 kg N ha−1 yr−1 with negligible losses after this). Detailed analysis of the role of N forms indicates that, for pan class="Chemical">forbs across all sites, stronger relationshipn>s were found with oxidised N (F1,50 = 16.0, P<0.001, Table 1: model 2). Significant interactions between N deposition and rainfall indicate that forb diversity decline is exacerbated in wetter areas receiving higher levels of N deposition. In the upn>lands, the relationshipn> was significant with total N deposition (F1,19 = 6.9, P<0.05, Table 1: model 1), whilst in the lowlands, NOy and NHx both had significant interactions with rain (F1,20 = 4.6, P<0.05, F1,19 = 6.2, P<0.05 respn>ectively) Table 1: model 2).

pan class="Chemical">Graminoid spn>ecies richness declined significantly across all sites with increasing total N inpn>uts (F1,50 = 9.3, P<0.01, Figure 2), as well as with increasing tempn>erature (F1,50 = 10.9, P<0.01). Losses associated with increasing N were highest at the lower end of the depn>osition range (an average of 4 spn>ecies lost per site between 5–10 kg ha−1 yr−1) with a slower rate of loss as inpn>uts increased further. Additionally, each incremental tempn>erature increase of 500 growing degree days correspn>onded to an average loss of 1 spn>ecies per site. Across all sites, both forms of N (oxidised and reduced) were associated with n>an class="Chemical">graminoid species diversity decline. However, when upland and lowland datasets were considered independently, stronger relationships with NOy were found in the lowland sites, and with NHx in the upland sites. Graminoid abundance, by contrast, increased with increasing N (F1,50 = 15.8, P<0.001), indicating a shift towards more grass dominated systems with greater N deposition (Figure 2); reduced N was the most significant driver of this response (F1,50 = 19.5, P<0.001, Table 2: model 2). Increases in grass abundance were, however, found to be chiefly characterised by a few species, specifically Molinia caerulea and Deschampsia flexuosa, the latter being significantly more abundant at sites receiving greater reduced N inputs (F2,49 = 9.2, P<0.01, Table 2: model 1).

Table 2

Summary of optimal models for vegetation responses in relation to N deposition and climate variables.

	Model 1 (total N deposition + environmental variables)									Model 2 (oxidized/reduced N + environmental variables)
Response variable	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P
Graminoid cover (%)	(+) N	50	***	(+) N	23	**	(+) N	22	**	(+) NHx	50	***	(+) NHx	23	**	(+) NHx	25	**
							gdd x rain		**
Bryophyte cover (%)	(−) gdd	50	***	(−) gdd	23	*	(−) gdd	25	**	(−) gdd	50	***	(−) gdd	25	*	(−) gdd	25	**
Lichen cover (%)	(−) N	47	**	(−) N	20	*	(−) N	25	*	(−) NHx	47	**	(−) NOy	19	*	(−) NHx	25	**
	(+) rain		*	(+) rain		**				rain x ht		***	(+) rain		***
	rain x ht		**	rain x ht		***							NOy x ht		*
													rain x ht		***
Canopy height (cm)	(+) N	49	***	(+) N	22	***	(+) gdd	25	*	(+) NHx	49	***	(+) NHx	22	***	(+) gdd	25	*
	(+) gdd		***	(+) gdd		***				(+) gdd		**	(+) gdd		***
Pleurozium (%)	(−) N	48	***	(+) rain	22	**	(−) N	24	***	(−) NOy	47	**	(+) rain	22	**	(−) NHx	24	**
	(+) rain		***	N x rain		***	(−) gdd		**	(−) NHx		***	NHx x rain		***	(−) gdd		**
	N x rain		***							(+) rain		***
										NHx x rain		***
Brachythecium (%)	(+) N	50	*	-	-	-	-	-	-	(+) NHx	50	*	-	-	-	-	-	-
Hylocomnium (%)	(−) N	49	**	(−) gdd	23	*	(−) gdd	25	*	(−) NHx	48	**	(−) gdd	23	*	(−) NOy	23	*
	(−) gdd		**							(−) gdd		**				(−) gdd		*
Deschampsia (%)	(+) N	50	***	(+) N	23	**	-	-	-	(+) NHx	49	***	-	-	-	(+) NHx	25	*
										(−) gdd		*

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; gdd  =  growing degree days; ht  =  Calluna canopy height.

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; pan class="Chemical">gdd  =  growing degree days; ht  =  n>an class="Species">Calluna canopy height.

Across both models and all sites, significant positive relationships were observed between pan class="Species">Calluna canopn>y height and tempclass="Chemical">n>erature (F2,49 = 12.6, P<0.001, Table 1: model 1). Whilst height increased in a linear fashion in relation to total N deposition across all sites (F2,49 = 15.3, P<0.001, Table 2: model 1), this relationshipn> was largely driven by NHy in the upn>lands (F2,22 = 24.8, P<0.001, Table 2: model 2).

Lower plant responses

Significant negative relationships between lichen species richness and total N deposition were found across the combined (all sites) dataset (F1,50 = 7.8, P<0.01, Table 1: model 1, Figure 3), driven largely by upland sites. An average of 5 lichen species per site was recorded in areas receiving less than 10 kg N ha−1 yr−1. This declined to an average of 1.4 species at input rates of 20 kg ha−1 yr−1, indicating that even modest rates of N loading are associated with a considerable loss of lichen species. Significant interactions indicate that the negative effects of N deposition on lichen richness were greater at sites with high temperature (F1,48 = 4.8, P<0.05, Table 1: model 1), partipan class="Chemical">cularly in the upn>lands (F1,20 = 8.9, P<0.01, Table 1: model 1).

Figure 3

Relationships between N deposition and lower plant species richness and abundance.

Declines in lichen species richness are related to increasing N deposition (A) whereas bryophyte richness per se is principally influenced by temperature (B). At a species level, however, significant species-specific relationships between N deposition and bryophytes were evident (C, D).

pan class="Disease">Declines in lichen spn>ecies richness are related to increasing N depn>osition (A) whereas bryopn>hyte richness per se is principn>ally influenced by tempn>erature (B). At a spn>ecies level, however, significant spn>ecies-spn>ecific relationshipn>s between N depn>osition and n>an class="Species">bryophytes were evident (C, D).

Considered independently, patterns in lichen abundance at upland and lowland sites were best explained by the different forms of N. The relationship with oxidised N was statistically stronger in upland sites (F5,19 = 4.5, P<0.05, Table 2: model 2), whilst the relationship was stronger with reduced N in the lowlands (F1,25 = 8.9, P<0.01, Table 2: model 2). Climatic conditions were also important (partipan class="Chemical">cularly in the wetter, upn>land areas), with lichen abundance generally greater in areas receiving more rainfall (F4,47 = 4.4, P<0.05, Table 2: model 1). Lichen cover was also generally lower in wetter areas with taller n>an class="Species">Calluna canopies (F4,47 = 9.1, P<0.001, Table 2: model 1).

Temperature was the main driver influencing bryophyte species richness and abundance across all sites (F1,50 = 21.3, P<0.001, Table 1: model 1, Figure 3), with less species recorded at the warmer end of the temperature gradient (<1000 pan class="Chemical">GDD  =  10 spn>ecies (on average) per site; > 2000 n>an class="Chemical">GDD  =  5 species (on average) per site). This overall relationship was driven by the strong decline at lowland sites experiencing the warmest conditions, with a non-significant relationship in the uplands dataset. Whilst a general downward trend in bryophyte diversity was apparent in relation to N deposition across all sites, the only significant N signal was found in the lowland sites in the form of oxidised N (F1,26 = 11.2, P<0.01, Table 1: model 2). Stronger relationships with N were, however, found at the species level: the characteristic heathland mosses Hylocomium splendens (Figure 3) and Pleurozium schreberii were seen to severely decline with increasing N (Table 2) whilst, conversely, the nitrophilous species [38] Brachythecium rutabulum (Figure 3) was more abundant at sites receiving higher N inputs (F1,50 = 6.1, <0.05).

Calluna tissue chemistry

The main factor affecting pan class="Species">Calluna tissue N concentrations across all sites was temperature (Figure 4). A significant quadratic relationshipn> was found between (increasing) tempn>erature and (higher) tissue N concentrations upn> to a threshold equivalent to 1500 n>an class="Chemical">GDD (F1,48 = 5.4, P<0.05, Table 3: Model 1). When upland and lowland datasets were considered separately, this relationship was only significant for the lowlands (F2,23 = 5.5, P<0.05).

Figure 4

The influence of N deposition and climate on heathland plant and soil chemistry.

Calluna foliar N concentrations were principally related to changes in temperature (A). However, positive relationships were seen between Calluna litter N concentrations and total N deposition (B). Soil carbon concentrations were higher in areas receiving greater oxidised N inputs and rainfall (C), and soil C:N ratios were also positively related to levels of oxidised N deposition (D).

Table 3

Summary of optimal models for plant and soil biogeochemical responses in relation to N deposition and climatic variables.

	Model 1 (Total N deposition + environmental variables)									Model 2 (oxidized/reduced N + environmental variables)
Response variable	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P	All sites	d.f.	P	Uplands	d.f.	P	Lowlands	d.f.	P
Soil pH	(−) N	49	**	(−) N	20	***	-	-	-	(−) NOy	47	***	(−) NOy	23	***	-	-	-
	(+) gdd		**	gdd x rain		*				(+) gdd		*
				N x rain		*				NHx x rain		*
Extractable NO3	gdd x moist	38	**	(+) moist	12	***	(+) moist	24	**	(+) moist	36	***	(+) NHx	10	*	(+) moist		*
										gdd x moist		*	(+) moist		*
Extractable NHx	(+) gdd	44	**	(+) N	11	**	-	-	-	(+) gdd	44	**	(+) NOy	10	**	-	-	-
				(+) moist		*							(+) gdd		*
Total soil N	(+) rain	45	***	N x rain	19	**	-	-	-	NOy x rain	44	***	(+) rain	18	*	(−) NOy	23	*
Total soil C	(+) rain	45	***	(−) rain	17	*	-	-	-	NOy x rain	44	***	NOy x rain	19	**	-	-	-
				N x rain		*
				gdd x rain		*
Soil C:N ratio	(−) rain	43	*	-	-	-	-	-	-	(−) rain	43	*	-	-	-	(+) NOy	22	**
Litter Ph-Oxidase	(−) gdd	44	*	-	-	-	gdd x rain	19	**	(−) gdd	44	*	(+) NHx	19	**	NHx x gdd	19	**
	gdd x rain		*							gdd x rain		*	NOy x rain		**
Litter PME	(−) rain	47	***	(+) N	21	*	-	-	-	(+) NHx	46	*	(−) rain	21	*	(+) NOy	23	*
	rain x N		**	(−) rain		*				(−) rain		***
										NHx x rain		**
Calluna litter N	(+) N	23	*	(+) N	2	*	-	-	-	-	-	-	(+) NHx	2	*	-	-	-
				(−) gdd		*							(−) gdd		*
Calluna litter C:N	(−) N	23	**	-	-	-	-	-	-	(−) NHx	23	**	-	-	-	(+) NHx	18	*
Calluna foliar N	(−) gdd	48	*	-	-	-	(+) gdd	20	*	(+) gdd	48	*	-	-	-	NHx x gdd	22	*
							N x rain		*

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; gdd  =  growing degree days; moist  =  soil moisture; PME  =  phosphomonoesterase.

pan class="Disease">Calluna foliar N concentrations were principally related to changes in tempn>erature (A). However, positive relationshipn>s were seen between n>an class="Disease">Calluna litter N concentrations and total N deposition (B). Soil carbon concentrations were higher in areas receiving greater oxidised N inputs and rainfall (C), and soil C:N ratios were also positively related to levels of oxidised N deposition (D).

(+/−)  =  direction of response; x  =  interaction; P values  =  * <0.05, **<0.001, ***<0.001; pan class="Chemical">gdd  =  growing degree days; moist  =  soil moisture; PME  =  phosphomonoesterase.

Soil pH

Soil pH ranged from 3.5 to 5.4, decreasing significantly in relation to total n class="Chemical">N deposition (F2,49 = 12.0, P<0.01, Table 3: model 1) and increasing at sites with higher temperatures (F2,49 = 7.5, P<0.01, Table 3: model 1) across all sites. The relationship with N appears to be driven by oxidised N (F4,47 = 13.7, P<0.001, Table 3: model 2), partipan class="Chemical">cularly in the upn>land sites.

Extractable soil nutrients

A positive relationship between reduced N deposition and extractable soil pan class="Chemical">nitrate (n>an class="Chemical">NO3 −) concentrations was found for upland sites (F3,10 = 9.6, P<0.05, Table 3: model 2). Across the combined, all sites dataset, however, interactions between temperature and soil moisture (F2,38 = 6.9, P<0.01) suggest that extractable NO3 − concentrations are greater in areas with clement climatic conditions (i.e. warm and wet).

Extractable soil pan class="Chemical">ammonium (pan class="Chemical">NH4 +) concentrations also increased with increasing total N in the upland sites (F2,11 = 11.9, P<0.01, Table 3: model 1), although this pattern was not significant in the lowlands. As with NO3 −, significant relationships with climate indicate that extractable NH4 + concentrations are higher in warmer areas (F4,9 = 8.6, P<0.05, Table 3: model 2); this relationship appears to be driven by the uplands dataset, where greater soil NH4 + availability is also associated with higher rainfall (F4,9 = 18.1, P<0.01, Table 3: model 2).

Total soil nitrogen and carbon concentrations

Significant interactions between total N deposition and rainfall were found for total soil N (F2,19 = 10.7, P<0.001, Table 3: model 1) and C (F4,17 = 6.5, P<0.05, Table 3: model 1) in upland sites, with the highest concentrations ocpan class="Chemical">curring at wetter sites receiving high N inpn>uts. Model 2 analyses show that oxidised N is driving this relationshipn> (N: F2,44 = 23.9, P<0.001; C: F2,44 = 9.3, P<0.001, Table 3, Figure 4). At lowland sites, a negative relationshipn> was seen between oxidised N and total soil N concentrations, although this was only significant at the 5% level (F1,23 = 6.3, P<0.05, Table 3: model 2). Climate is clearly also an impn>ortant driver of soil C and N, with significantly higher concentrations at sites receiving more rainfall (N, F1,45 = 15.1, P<0.001; C: F1,45 = 13.3, P<0.001, Table 3: model 1). Indeed, at upn>land sites, an interaction between tempclass="Chemical">n>erature and precipn>itation indicates that soil N and C concentrations are highest where rainfall and tempn>erature are both high (N: F3,18 = 9.7, P<0.01, model 2; C: F4,17 = 5.8, P<0.05, model 1). At lowland sites, soil C:N ratios were highest at sites receiving greater oxidised N inpn>uts (F2,22 = 9.0, P<0.01, Table 3: model 2, Figure 4).

Litter chemistry, enzyme activity and nutrient turnover

Across all sites, litter N concentrations were significantly higher in areas of higher N deposition (F1,23 = 5.3, P<0.05, Table 3: model 1, Figure 4). This relationship was driven by data from upland sites, with the strongest patterns seen in relation to reduced N (F2,2 = ,34.0, P<0.05, Table 3: model 2). Litter C:N ratios were lower at sites with the highest N deposition (F1,23 = 8.6, P<0.01, Table 3: model 1); this relationship was strongest in the lowlands for reduced N (F1,18 = 5.5, P<0.05, Table 3: model 2), with a similar (non-significant) trend apparent for the uplands.

Overall, litter phenol oxidase activity was more strongly related to climate than to levels of N deposition. Across the full (all sites) data set, enzyme activities were lower at warmer sites (F2,44 = 6.1, P<0.05, Table 3: model 1); a significant interaction was also apparent between temperature and rainfall, with the lowest activities recorded from areas where higher temperatures combined with greater quantities of rainfall (F2,44 = 6.7, P<0.05, Table 3: model 1). In upland areas, phenol-oxidase activity was, however, higher at sites receiving greater inputs of reduced N (F3,19 = 8.3, P<0.01, Table 3: model 2, Figure 5), as well as those with high levels of both rainfall and oxidised N deposition (F3,19 = 12.2, P<0.01, Table 3: model 2).

Figure 5

Relationships between N deposition and litter enzyme activity.

Increasing reduced N deposition was positively correlated with increased litter phenol-oxidase activity in upland sites (A) and litter phosphomonoesterase (PME) activity across all sites (B).

As for phenol-oxidase, litter PME activity across all sites was significantly related to climate, with activity decreasing as precipitation levels increased (F2,47 = 15.0, P<0.001, Table 3: model 1). A positive relationship between total N deposition and litter PME activity was, however, apparent for upland sites (F2,21 = 4.3, P<0.05, Table 3: model 1). Analysis of N forms revealed that this was driven by reduced N and that this relationship was also significant for the all sites dataset (F3,46 = 4.3, P<0.05, model 2, Figure 5).

Discussion

Above-ground responses

A decline in plant species diversity in relation to increasing atmospheric N deposition has been observed in a variety of ecosystems ranging from acid and calcareous grasslands [39], [40] to boreal and temperate forests [41], [42] and even pan class="Disease">arctic tundra [43]. Measurement of plant spn>ecies richness can provide a reliable indication of habitat quality [44], [45] and is considered to be a fundamental tool in assessing community and wider-scale diversity [46]. Where spn>ecies losses are countered by increases in a few, dominant and/or invasive spn>ecies, the risk to habitats are homogenisation with class="Chemical">n>an class="Disease">declines in genetic-, ecosystem- and landscape-scale diversity. The findings from this survey provide clear evidence that elevated N deposition significantly reduces the floristic diversity of heathland ecosystems, either solely or in combination with climatic drivers – particularly increasing temperature – a major driver of global biodiversity loss in both terrestrial and non-terrestrial ecosystems [4], [47]. The nature of the relationship between total N deposition and species richness suggests similar patterns in both upland and lowland systems, namely that the greatest proportional reduction of species richness in response to increasing N inputs is seen at sites currently experiencing the lowest deposition loads. It appears that both oxidised and reduced forms of N are contributing to heathland biodiversity loss, although, at the functional group level there were more significant relationships between species richness and reduced N, compared to oxidised N, as hypothesised. This is in line with findings from similar surveys elsewhere for grassland [48) and heathland [19], providing further support for an additional impact associated with the acidifying properties of reduced N and/or direct toxicity of ammonium ions, over and above that of N-enrichment per se. Given the broad cover of spatial, climatic and pollution indices, the magnitude of diversity loss that has been observed amongst forb, graminoid, bryophyte and lichen functional groups in this study not only validates experimental findings from many N manipulation studies [39], [49] but, more profoundly, provides compelling evidence that N deposition is driving biodiversity loss amongst UK heathlands at a national scale.

In N-limited environments, increased N availability frequently leads to increased plant growth and/or expansion of nitrophilous species, potentially leading to competitive exclusion of slower growing habitat-specialist species adapted to low levels of N availability [40]. Effects of N on soil pH, partipan class="Chemical">cularly in its reduced form [50], can also result in the decline of n>an class="Species">sensitive plant species, particularly those with narrow pH niches or low tolerance of high NH4 + ion concentrations [51], [52]. In the current study, Calluna plants were taller and soil pH lower in areas receiving high N inputs, suggesting that reduced understorey light availability and/or unfavourable soil pH may be mechanisms contributing to biodiversity loss in British heathlands. Whilst bryophyte cover was not related per se to N deposition, shifts within bryophyte communities, notable at the species level, were significantly related to N inputs. For example, the characteristic heathland moss species Hylocomnium splendens and Pleurozium schreberii were seen to decline as N deposition increased, whilst Brachythecium rutabulum, a nitrophilous species more commonly associated with base rich, eutrophicated habitats [53] was more abundant at sites with higher N deposition. Similar species shifts have been seen in boreal forest systems [41], [54] with H. splendens declining, and being replaced by B. rutabulum, at relatively modest input rates of 10 kg N ha−1 yr−1.

Although pan class="Chemical">graminoid abundance increased significantly with increasing N deposition, a concurrent decrease in graminoid diversity, coupled with sharp pan class="Disease">declines in the abundance of typical heathland grass species (e.g. Festuca ovina, Danthonia decumbens) indicates that grass invasion at higher N deposition loads is dominated by only a small number of species (specifically Deschampsia flexuousa or Molinia caerulea). These findings support earlier observations of D. flexuousa and M. caerulea increases in high N deposition areas in the Netherlands [55] and an associated conversion of 35 % of Dutch heathlands to grasslands in the 1980s [56]. These species have also been found to increase in Danish bog systems exposed to N levels of 10–15 kg ha−1 yr−1, to the detriment of characteristic ombrotrophic vegetation [57].

N deposition is associated with negative impacts on lichen communities in pan class="Species">many ecosystems [58]–[60], a phenomenon often linked to increased compclass="Chemical">n>etition for light as a result of N-driven stimulation in higher plant growth [9], [61]. However, direct n>an class="Disease">toxicity of NO2 −, pan class="Chemical">NO3 − and NH4 + is also likely to be a contributing factor in the decline of non-vascular plants [16], [62], [63]. Whilst lichen abundance was negatively related to both reduced and oxidised forms of N in this study, declines in both the cover and species richness of lichens were statistically stronger with reduced N. As with most higher plants, lichens can effectively assimilate NH4 + at modest levels [64] rapidly converting it into amino acids to mitigate its cytotoxicity. However, for acidophytic species (such as those found in heathlands) this assimilation capacity is limited and toxic levels of unassimilated NH4 + can accumulate in the lichen thallus [65]. It seems likely, therefore, that both competitive exclusion and potentially also direct toxicity effects are contributing to the decline in abundance and diversity of lichens observed with increasing N deposition in the current study.

Below-ground responses

At all sites in this study, litter N concentrations increased, and C:N ratios decreased, with increasing N inputs suggesting greater availability of pan class="Chemical">organic N for soil microorganisms and thus faster rates of nutrient cycling. Microbial enzyme activity provides a useful proxy for decompn>osition rates and has been frequently used to monitor the functional respn>onses of microbial populations to global change perturbations [10], [66]. Our results show that overall the activity of phenol-oxidase – an enzyme respn>onsible for breakdown of the structural compn>onents of plant material [67] – in litter was greater at sites receiving higher inpn>uts of reduced N. This relationshipn> was seen in both upn>land and lowland systems, the latter, however, stronger in areas of greater rainfall. Significant interactions between oxidised N and rainfall in upn>land systems were also apn>parent, possibly reflecting the wider spn>an of oxidised N inpn>uts (3–15 kg ha−1 yr−1) compn>ared to lowland sites (3–9 kg ha−1 yr−1). The pattern of increased litter/soil enzyme activity with increasing N depn>osition supn>ports results from a recent survey of biochemical indicators in lowland heathlands [19]. However, contrasting results have been repn>orted previously for moorlands [18] and tempn>erate forests [66], [67].

The positive relationship between N deposition – partipan class="Chemical">cularly reduced N - and litter PME activity provides the first field scale evidence of an increase in P depan class="Species">mand across both upland and lowland heathland sites receiving high rates of N deposition. These findings are in line with those from manipulation experiments across a range of ecosystems [68], [69], and with survey evidence of increased plant uptake of P across transects of increasing N deposition in lowland heathland [15], [70].

The main factors regulating decomposition are moisture and temperature [6], [71]; this is also apparent in our survey results with enzyme activities highest in warmer, wetter areas. Climate change predictions for warmer, drier summers and cooler, wetter winters will, therefore, have implications for the activities of enzymes involved in the decomposition process, with likely consequences for nutrient availabilities and uptake. Concentrations of soil extractable pan class="Chemical">NO3 and n>an class="Chemical">NH4 were significantly higher with increasing N deposition in the upland sites, although no such relationship was evident in the lowlands. The strong influence of climate on rates of nutrient turnover and availability [72], [73], and generally lower annual rainfall in lowland (600–1300 mm) compared to upland (800–2500 mm) sites, may explain the lack of relationship with N in lowlands and the generally higher abundance of litter in upland sites. This may be particularly the case in this survey since sampling was carried out in the summer when soil moisture levels are at their lowest and plant and microbial immobilisation of N are at their highest [74]. Caution must be exercised when making assumptions based on ‘snapshot’ measurements of mineral pools (particularly across broad temporal and spatial gradients) that are inherently subject to large inter-seasonal and inter-annual variations [75]. Furthermore, recent N-focused surveys carried out by Jones & Power [19] and Stevens et al. [76] found that soil type, underlying geology, climatic and altitudinal variation were important drivers of soil and plant responses, influencing the nature of relationships with N deposition.

Higher soil C:N ratios at sites with higher N inputs, specifically in the lowlands, indicate a shift towards immobilisation with increasing N deposition. Lower litter C:N ratios at high N sites contrast with this pattern, although may be explained by the previously mentioned N-driven stimulation at the early stage of decomposition where labile, energy rich substrates are abundant, followed by a slowdown in later stages as the relative proportions of recalcitrant compounds (such as pan class="Chemical">lignin) progressively increase [77]. Whilst many studies of N deposition impacts have found that total soil N content increased with elevated N inputs [69], [78], [79], no evidence of N or C accumulation in relation to elevated N as a sole driver was found in this survey. This may reflect the depth (0–15cm) of the soil core sampled, which not only results in high levels of between site variability in depth of organic layer included in the sample (from <1cm to 15 cm), but may also be too deep to capture nutrient signals from recent decades that are present in the upper horizons.

Across all sites, the most significant driver influencing soil total N and C content was precipitation. Since climatic factors are known to strongly influence the production of plant material (and, subsequently, plant residues) [80], the dominant influence of rainfall is not surprising. Indeed, our results suggest that climatic variables are important drivers of pan class="Species">many ecosystem processes and, therefore, are likely to influence respn>onses to N depn>osition. Given that climate models predict altered precipn>itation and tempn>erature patterns for the northern hemispn>here [81], and that widespn>read exceedance of critical N loads is predicted across Europn>e for the foreseeable future [2], [82], the impn>lications of global change for the biodiversity, structure and function of these internationally impn>ortant heathland ecosystems are very considerable.

Conclusion

The findings from this survey provide novel support to a growing body of evidence demonstrating a strong link between N deposition and an increase in nitrophilous plant species across a wide range of habitats [17], [22]. An associated decline in the diversity of both higher and lower plants across broad spatial and climatic gradients – in both upland and lowland environments - highlights the ubiquity of relationships with N, providing clear field evidence of widespread N-driven changes in the composition, structure and diversity of British heathlands. There is also strong evidence to suggest that pan class="Chemical">current rates of N depn>osition are altering many biogeochemical processes, particularly rates of nutrient cycling and availability, and that these relationships are influenced by climatic conditions. Litter N concentrations, litter enzyme activity (phenol-oxidase and PME) and, for the uplands, soil extractable N concentrations, all appear to have consistent relationships with total N deposition and may, therefore, be useful biochemical indicators of, and responses to, N deposition. Given future predictions for climate change and continued wide scale exceedance of critical N loads across Europe, further effects of these global change phenomena on heathland diversity and associated ecosystem services can be expected.