Biotic and Abiotic Factors Interact to Regulate Northern Peatland Carbon Cycling

Understanding the spatio-temporal variability of controls on peatland carbon (C) cycling is essential to project the effects of future environmental change. While there is understanding of individual drivers of C cycling, the effect of multiple drivers, including interactions, remains poorly understood. Using a spatially and temporally explicit sampling framework, we examined the effects of biotic and abiotic controls on key indicators of peatland functioning: ecosystem respiration (R eco), photosynthesis (P cal), net ecosystem exchange (NEE), methane (CH4) fluxes, and pore water dissolved organic carbon concentration ([DOC]). Measurements were made over 12 months in a blanket peatland hosting a wind farm in Scotland, UK. Overall, we found that (i) season and plant functional type (PFT) explained most variation in R eco and P cal, (ii) PFT and spatial location within the wind farm, which integrates several peat properties, were dominant predictors of CH4 fluxes, and (iii) season and location within the wind farm correlated with pore water [DOC]. Examination of predictors indicated that interactions, between and within biotic and abiotic factors, explained a significant amount of variation in greenhouse gas fluxes and [DOC]. These findings indicate that combinations of biotic and abiotic factors could mediate or exacerbate the effects of future environmental change on peatland C cycling. Given this, studies of C cycling need to capture the spatial and temporal variance of biotic and abiotic factors and their interactions to project the likely impacts of environmental change.


INTRODUCTION
Peat soils are globally important carbon (C) stores, with approximately one-quarter of the world soil C stock, 547 Gt C, stored in northern peatlands (Jobbá gy and Jackson 2000; Yu and others 2010). The C balance of a peatland is dependent on fluxes of C to and from the system, primarily the influx and efflux of greenhouse gases (GHG), carbon dioxide (CO 2 ), and methane (CH 4 ), and the production and export of dissolved organic and inorganic carbon (DOC, DIC) (Kutsch and others 2009). These fluxes reflect biological processes by which C is assimilated and processed. CO 2 is released through decomposition of organic matter, plant respiration, and the oxidation of CH 4 (positive fluxes, together hereafter termed ecosystem respiration) and taken up through photosynthesis (negative fluxes). The balance of ecosystem respiration and photosynthesis is net ecosystem exchange and indicates if C is being sequestered (negative flux) or released (positive flux). CH 4 is produced during anaerobic decomposition but can be oxidized into CO 2 by methanotrophs. Methanotrophs can also oxidize atmospheric CH 4 , leading to negative fluxes. DOC can be produced under aerobic and anaerobic conditions and exported from peatlands by drainage channels, but also acts as a substrate for microbes, with CO 2 efflux from streams comprising an important component of the overall C balance (Billett and Harvey 2013). Measurements of GHG fluxes represent the contributions of the sampling plots, whereas [DOC] values are a result of production (including mobilization into the dissolved phase), microbial degradation, input of water from precipitation at the sampling point, and inflow and outflow of water and DOC from upslope and downslope. Consequently the pore water [DOC] values give an indication of both potential DOC runoff [although discharge is known to dominate fluxes (Armstrong and others 2010)], and substrate available for microbial decomposition which may feed back to GHG fluxes.
Characteristics of the biotic communities, both plant and soil, influence the assimilation and processing of C in peatlands. The composition and structure of the vegetation regulates C cycling as different plant functional types (PFTs) have specific traits that dictate how they cycle C (Díaz and others 2007;De Deyn and others 2008;Ward and others 2009;Armstrong and others 2012). Traits can be ecophysiological (for example, their nitrogen (N) use efficiency), biochemical (for example, how they influence litter nutrient content), demographical (for example, diameter of growth), morphological (for example, root architecture), or phenological (for example, flowering time) (Lavorel and others 2007). Many plant functional traits are associated with growth rate, with fast-growing plants typically characterized by a greater photosynthetic capacity, shorter lifespan, lower dry matter content, and greater and higher quality litter inputs than slower growing species (Aerts and Chapin III 1999). Peatlands are dominated by shrubs, sedges, and mosses, with sedges relatively fast growing and shrubs relatively slow growing. Soil faunal and microbial communities are also central to peatland C cycling, stimulating decomposition and nutrient supply, with consequent effects on productivity and plant community composition (Bardgett and others 2008;Weedon and others 2012;Bragazza and others 2014).
Biological C cycling processes respond to abiotic conditions, with climate established as a dominant control (Freeman and others 2001;Dorrepaal and others 2009). As a consequence, much research has considered the response of both GHG and pore water [DOC] to temperature and water table (including their interactions) (Bubier and others 2003;Holden 2005;Ise and others 2008;Clark and others 2009). In northern latitudes, there are strong seasonal controls over GHG exchange and [DOC]: warmer temperatures generally enhance C cycling in summer (Ward and others 2007;Clark and others 2008) while seasonally lower water tables decrease photosynthesis rates and net CH 4 emissions, and increase respiration and DOC production (Blodau and others 2004). However, increases in DOC production may not be apparent in pore and stream water DOC samples until autumn when the water table commonly rises causing an 'autumn flush' of DOC produced over the summer (Grayson and Holden 2012).
Other abiotic properties also influence the biological processing of C in peatlands. Physical properties, such as peat depth and bulk density, affect C cycling through controlling substrate availability and the diffusion rates of water, compounds, and gas through the peat profile (Dorrepaal and others 2009;Levy and others 2012). The chemical composition of the peat, including nutrient availability and the level of secondary metabolites, also exerts a significant influence on both productivity and decomposition (Bragazza and others 2006). For example, 7 years of N addition decreased C:N, increased bacterial biomass, and stimulated soil enzyme activities in an Italian alpine peatland (Bragazza and others 2012).
There are known interactions between and within biotic and abiotic controls, for example, the thermal conditions of peat are affected by the water table depth (Kettridge and Baird 2008). Further, litter quality, a plant functional trait and thus a biotic property, influences the chemical composition of the peat (Bragazza and others 2013), an abiotic property. Moreover, physical and chemical peat properties may influence plant and microbial community composition and thus indirectly affect peatland C cycling (De Deyn and others 2008). However, these interactions are poorly resolved as initial studies examining the effect of environmental change on peatland C cycling focused on single drivers (Freeman and others 1992). More recently the result of multiple drivers changing simultaneously has been elucidated by examining more than one driver and the interactions between them (Ellis and others 2009), but knowledge of interactions between abiotic and biotic factors remains poorly understood (Ward and others 2013). Moreover, many studies imposed treatments and thus results may be indicative of the short-term transition state rather than longer term effects (Bubier and others 2007).
Interactive effects, both between and within biotic and abiotic properties, maybe synergistic or antagonistic. For example, plant community composition can modulate the effects of temperature warming, necessitating projections of future peatland C sink strength to take into account any concurrent changes in vegetation (Ward and others 2013). Further, nutrient addition can cause a reduction in photosynthesis despite an increase in shrub biomass and leaf area, as a result of reduced moss cover and increased shading by shrubs and litter accumulation (Bubier and others 2007). Therefore, there is a critical need to examine the role of biotic and abiotic factors and their interactions play in the functioning of natural peatland ecosystems. The understanding generated will provide valuable insight into the potential implications of future environmental change. Further, it will provide a firmer basis for inference of the effect of environmental change on peatlands with different biotic and abiotic characteristics, and the implications of natural or managed vegetation change.
The overarching aim of this study was to determine the role of abiotic and biotic factors as regulators of peatland GHG fluxes and pore water [DOC]. To achieve this, we established a factorial field experiment and subsequently measured GHG fluxes (that is, R eco , P cal , NEE, and CH 4 ) and pore water [DOC] for 12 months to assess how ecosystem C cycling processes varied temporally and spatially. First, we examined the effects of season, PFT, and site, which integrates abiotic and biotic factors that vary spatially across the peatland. Second, biotic and abiotic metrics were used to explain variances in carbon cycling, including peat temperature, water table depth, PFT biomass (that is, shrubs, sedges and mosses), peat depth, bulk density, and peat C and N contents. We hypothesized that (1) PFT explains more of the variation in ecosystem CO 2 and CH 4 fluxes and [DOC] than spatial variation 'site'; and (2) interactions, both between and within biotic and abiotic factors, explain a significant proportion of the variance in GHG fluxes and [DOC] in northern peatlands.

Study Site
This study was undertaken at Black Law Wind Farm, Scotland (55°46¢01¢¢N 03°44¢20¢¢W, altitude 250-320 m). The wind farm covers 18.6 km 2 of blanket bog, grassland, and plantation forestry. The blanket bog is typical of UK peatlands in terms of microtopography and vegetation [National Vegetation Community M19 Calluna vulgaris -Eriophorum vaginatum blanket mire (Rodwell 2000)]. Four sites, coded from A to D, along a 3.5 km SW to NE transect, were established in areas of blanket bog. The sites were selected to be broadly visually similar, but given the inherent spatial variability of the blanket bog, sites A and B had thinner peat, and were less hummocky and shrubby compared with sites C and D. Within all four sites, the dominant PFT of mosses, sedges and shrubs were found, with the most prevalent species being Sphagnum spp. (hollow), Calluna vulgaris (hummock), and Eriophorum angustifolium (hummock). Within each site, four blocks were established, each dominated by one of the three PFTs (n = 12 at each site; n = 48 in total). A wooden sampling platform was installed at the plots to minimize disturbance during site visits. At each sampling plot, NEE, R eco , CH 4 flux, pore water [DOC], and water table depth were measured monthly for 1 year commencing from May 2011; peat temperature was recorded every 30 min at three of the four blocks and vegetation biomass and peat properties were characterized once. The vegetation cover at Black Law was very patchy, thus there was a minimum buffer area (that is, the same PFT) of 5 cm around the monitoring equipment (plots were at least 40 cm by 40 cm).
Vegetation Biomass Above-ground biomass within the GHG sampling plots was estimated in August (time of peak biomass) 2011 by selecting a similar plot nearby (to avoid destructive harvesting of the GHG sampling plots), harvesting all of the vegetation (including other PFTs present), oven drying it at 50°C, and weighing the dry mass of each PFT.
Peat Properties Peat temperature was recorded at 5 cm depth using Onset Hobo Pendant temperature loggers (Onset, USA). C and N contents of peat surrounding each block were evaluated from a 0.1 g homogenized and oven-dried (65°C) subsample of a 15 cm long, 5 cm diameter peat core using a LECO Truspec CN Analyzer (LECO, USA) and C:N was calculated (Carter and Gregorich 2007). Dry bulk density was calculated from a core of the same size, with the peat dried at 105°C (Carter and Gregorich 2007). There may be varia-Biotic and Abiotic Factors to Regulate Carbon Cycling tion in temperature, C, N, and bulk density with depth. Average peat depth for each plot was calculated from triplicate probed peat depths within 100 cm of each plot.
Greenhouse Gas Fluxes CO 2 and CH 4 (hereafter abbreviated to GHG) fluxes were measured using static chambers, as detailed in Ward and others (2007), at each of the 48 plots. Six weeks prior to the first measurement, 30 cm diameter plastic collars were carefully inserted to the litter-peat interface, minimizing root damage. If there was not an airtight seal, peat from the surrounding area was used to create a seal on the outside of the collar. R eco and NEE were measured using an EGM-4 portable infrared gas analyzer (PP Systems, USA) and CH 4 by taking samples for analysis in the laboratory using a Perkin Elmer Autosystem XL Gas Chromatograph, as detailed in Ward and others (2013). R eco and CH 4 were measured using dark chambers and NEE using a transparent chamber. The measurements were quality controlled and those affected by poor chamber seals, leaking vials or instrument malfunction removed. Fluxes, as mg CO 2 -C or CH 4 -C m -2 h -1 , were calculated following Levy and others (2012).
Photosynthesis (P cal ) was calculated by subtracting R eco from NEE: positive NEE indicates land-toatmosphere exchange dominates and negative indicates atmosphere-to-land, that is, the peat is a source and sink, respectively. Extreme CH 4 fluxes, ascribed to measurement error (Christiansen and others 2011) or ebullition events (Baird and others 2004), were excluded using thresholds of -1.15 to 11.55 mg m 2 h -1 CH 4 as prescribed by Levy and others (2012). Consequently, our results are only applicable to non-ebullition CH 4 fluxes.
Hydrology Water table depth was measured within 30 min of GHG measurements using a dip well (2.2 cm diameter plastic tubing with perforations along its length below ground). Water table values above 0 cm represent a saturated litter layer or water ponding on the surface. Pore water samples (peat surface to 15 cm depth) were collected from a second dip well (2.2 cm diameter perforated plastic tubes with sealed bases). The pore water dip wells were emptied and the refill sampled after 24 h. The samples were stored at 4°C, filtered to 0.7 lm within 36 h, and analyzed for [DOC] within 2 months. The samples were analyzed for [DOC] by high temperature oxidation using a Thermalox Total C analyzer. Potassium phthalate was used to produce standards and the precision of sample analysis was less than 0.5 mg l -1 C.

Data Analysis
For statistical analysis, the temperature closest to the mid-point of the GHG measurement time was selected but for [DOC] the average temperature over the 24 h preceding collection was calculated (as the sample was an integration of water movement in the previous 24 h). The sampling dates were categorized into spring (March-May), summer (June-August), autumn (September-November), and winter (December-February) (Jenkins and others 2008).
To test our hypothesis that PFT explains more of the variance in ecosystem CO 2 and CH 4 fluxes and [DOC] than the integrated effect of location on the peatland, repeated measures mixed effects models were used to explain the importance of season, PFT, and site, and their interactions for total biomass, peat temperature, water table depth, bulk density, peat depth, C, N, C:N, R eco , P cal , NEE, CH 4 , and [DOC]. This analysis was chosen as it provides a broad scale overview of the relative importance of groups of factors that vary with season (climate), PFT, and site (peat biological, and chemical and physical properties, some of which were not measured in this study). The distributions of response variables were assessed and the data transformed to promote normality as required: R eco , P cal , CH 4 , and [DOC] were logarithmically transformed, and the remainder were not transformed. To enable logarithmic transformation (that is, to eradicate negative values), P cal was multiplied by -1, 35 was added to each CH 4 measurement, and 1000 added to each NEE measurement. The distribution of residuals and residual versus fitted data plots were constructed for each model. Individual drivers and interactions with a P value <0.05 were deemed as significant. The effect size of the individual drivers and interactions were inferred from v 2 values.
To test our hypothesis that interactions, both between and within biotic and abiotic factors, explain a significant proportion of the variance in GHG fluxes and [DOC] in northern peatlands, generalized linear models (GLM) with a log-link and gamma distribution were used to assess the role of the biotic (biomass of each PFT) and abiotic (temperature, water table depth, bulk density, peat depth, C, N, and C:N) factors and their interactions in driving GHG fluxes and [DOC]. This approach was chosen to explore the importance of interactions over peatland C cycling. To allow interpretation of both the individual drivers and interactions, including their relative effect sizes, the explanatory variables were transformed (centered and standardized following Schielzeth (2010)): the coeffi-cients are indicative of effect size. As before, the P cal data were multiplied by -1, 35 were added to each CH 4 measurement, and 1000 added to each NEE measurement to allow the GLM models to run with a log-link. To avoid model over-specification given the large number of explanatory variables, the modeling was undertaken in three stages: (1) check for a relationship between the response variables and explanatory variables in turn and exclude nonsignificant factors from model development; (2) build a GLM model with the remaining explanatory variables and two-way interactions, ensuring that no explanatory variables were correlated (threshold of r = 0.80 used); and (3) step-wise removal of insignificant individual drivers (p > 0.05) and interactions. Model fit plots were assessed and pseudo-R 2 was calculated from the observed and predicted values. To assess the importance of interactions, the coefficients (absolute values) were summed and compared. To explore the interactions, plots of the predicted GHG flux or [DOC] using the mean of the explanatory variables and one standard deviation above and below were drawn. Stata13 was used for all data analysis (Sta-taCorp. 2013).

RESULTS
The mixed effects models, based on season, PFT, and site, explained between 37 and 55 % of the variation in [DOC] and GHG fluxes, with the exception of NEE, for which only 21 % of the variation was explained ( Table 1). The GLM models, based on biotic and abiotic factors and their interactions, explained between 33 and 59 % of the variance in GHG fluxes and [DOC], again except for NEE for which only 10 % was explained (Table 2).
Below we describe our results in relation to our hypotheses that (1) PFT explains more of the variation in ecosystem CO 2 and CH 4 fluxes and [DOC] than the integrated effect of properties that vary with location on the peatland and (2) interactions, both between and within biotic and abiotic factors, explain a significant amount of the variance in GHG fluxes and [DOC] in northern peatlands.
The Importance of Season, PFT, and the Integrated Effect of Spatially Varying Properties The biotic and abiotic metrics varied with season, PFT, and site, with water table varying with all three and biomass with PFT and site ( Table 3). As expected, given the northern latitude of Black Law, season explained the most variation in temperature and water table, with the warmest temperatures and lowest water tables observed in summer, coldest temperatures in winter, and highest water tables in autumn (Tables 3, 4). Site explained significant variability in bulk density, water table, peat depth, biomass, peat C and N contents, and C:N and, with the exception of water table, site explained more of the variance in these factors than season or PFT (Table 3). Site A had the shallowest peat, lowest peat C content (which was similar at the other sites), highest bulk density and C:N, and intermediate water  The average values of each GHG and [DOC] are given in Table 5. n = number of data points; R 2 = pseudo R 2 calculated from the square of the correlation between the observed and actual values; Df = degrees of freedom; v 2 = chi 2 value; P = P value.
the lowest bulk density, and intermediate peat depth (Tables 3, 4). Biomass and water table were the only factors that varied significantly with PFT, although for water table the effects of site and season were greater (Table 3). We found that shrubs had the greatest biomass and mosses the lowest and the water table was deepest under shrubs and shallowest under mosses (Table 4).
On average the dominant PFT constituted 65, 84, and 84 % of the total biomass in the shrub-, sedge-, and moss-dominated plots, respectively. Within plots they dominated, sedges and mosses comprised a minimum of 65-70 % of the total biomass, respectively, and with the exception of four plots (within which the shrubs died off and thus the biomass was not included), the shrub dominated plots had a minimum shrub biomass of 60 % of the total.
Unsurprisingly season explained most variance in R eco , P cal , and NEE: R eco and P cal were lowest in winter and highest in summer and, on balance, most C was sequestered in the autumn (that is, NEE was negative) (Tables 1, 5; Figure 1). For    The average values of biotic and abiotic properties are given in Table 4. n = number of data points; R 2 = pseudo R 2 calculated from the square of the correlation between the observed and actual values; Df = degrees of freedom; v 2 = chi 2 value; P = P value.  [DOC], season also dominated, with concentrations generally increasing from winter through to autumn (Tables 1, 5; Figure 1). In contrast season explained the least variance in CH 4 fluxes, with CH 4 peaking in autumn (Table 1, 5; Figure 1). Season also interacted with PFT for all models except CH 4 and [DOC] and interacted with the variation in properties integrated by site for all models except CH 4 (Table 1). For example, R eco increased in magnitude and variability in shrub and sedge plots in the summer, while rates from moss plots remained relatively low throughout the year (Figure 1). PFT explained more variation in R eco and P cal than the variation in properties integrated by site; R eco and P cal were significantly greater and more variable from shrubs and sedges than mosses (Table 1, 5; Figure 1). Mosses were generally a sink of CO 2 , sedges were the largest sink (although they were predominantly weak sources during winter) and shrubs were the largest source (especially during the summer months) (Table 5, Figure 1). PFT explained the least variance in CH 4 fluxes, with greater and more variable CH 4 fluxes from sedges compared with shrubs and mosses (  (Table 1). CO 2 fluxes were broadly similar across sites with some significant differences driven by greater variability during the spring and summer (Figure 1; Table 5). Specifically, R eco was greater at site B than site C with differences most apparent in the summer, P cal was higher (that is, greater CO 2 uptake but lower values as P cal is negative) at site A than C, most notably during the spring and less C was sequestered at site B compared with sites A and C, especially during the summer (Table 5; Figure 1). Sites A and B had significantly lower CH 4 fluxes than C and D (Tables 4, 5; Figure 1). Finally, [DOC] generally increased from site A to D although site B had similar [DOC] to sites C and D ( Table 5).

The Importance of Interactions
Interactions explained more of the variance in P cal , NEE, and CH 4 than individual drivers and contributed notably to R eco and [DOC], as indicated by the sum of the main and interactive effect coefficients ( Table 2).The majority of the interactive drivers for R eco included temperature, with the response rate (that is, change in R eco for a unit increase in temperature) increasing with temperature at a greater rate at deeper water table depths, shallower peat depths, lower moss biomass, and higher shrub biomass (Table 2; Figure 2A). Out of the interactive drivers, the shrub-moss interaction had the greatest effect for R eco ; ecosystem respiration was positively related to shrub biomass but increased to a greater extent where moss biomass was also high ( Figure 2B). Temperature did not interact with as many drivers for P cal ; instead the majority of interactions featured a biotic variable ( Table 2). The shrub-sedge interaction was dominant with the effect size (as indicated by the coefficient) double that of any other interactions (Table 2). P cal decreased with increasing shrub Table 5. Average ± Standard Deviation R eco (mg C m -2 h -1 ), P cal (mg C m -2 h -1 ), NEE (mg C m -2 h -1 ), CH 4 (mg C m -2 h -1 ), and [DOC] (mg l -1 ) for the Different PFT, Sites (A-D), and Seasons  The upper whiskers on these plot represent the data point that is £ 75th percentile + (75th percentile -25th percentile) 9 1.5 and the lower whisker the data point that is £ 25th percentile -(75th percentile -25th percentile) 9 1.5 biomass at high sedge biomass, increased at low sedge biomass, and remained constant with mean sedge biomass ( Figure 2C). As for R eco , the rate of increase in P cal with temperature was greater at shallower peat depths and also at lower C:N ( Figure 2D). For CH 4 , as for P cal , the shrub-sedge interaction explained the most variance (and dominated over the individual drivers) ( Table 2). CH 4 flux decreased with increasing shrub biomass and at a greater rate when sedge biomass was higher. However, at sites with the lowest sedge biomass, CH 4 flux increased with increasing shrub biomass ( Figure 2E). The other dominant interactions for CH 4 included peat C and bulk density (Table 2). CH 4 decreased with bulk density at a faster rate if Figure 2. Interaction plots for the R eco , P cal , CH 4 , and [DOC] GLM models. The plots are for the transformed data and thus 0 = the mean of the biotic or abiotic property, -1 is one standard deviation below the mean and 1 is one standard deviation above the mean, and P cal has been multiplied by -1. SD = standard deviation; PD = peat depth; Mo = moss biomass; Se = sedge biomass; CN = C:N; BD = bulk density peat C was lower or if C:N or sedge biomass was higher. CH 4 fluxes did not vary with peat C if sedge or shrub biomass was intermediate. However, at higher sedge and lower shrub biomass, CH 4 fluxes decreased with increasing peat C and vice versa.
Compared to GHG fluxes, interactions were not as important for [DOC], both in terms of the number and effect size ( Table 2). The dominant interaction was temperature-bulk density, with [DOC] increasing with rising temperatures at a greater rate at higher bulk densities ( Figure 2F).
[DOC] also increased with increasing peat C but at lower rates at the higher bulk densities.

DISCUSSION
Our field experiment was designed to examine the individual and interactive effects of critical biotic and abiotic drivers of plant peat C cycling across a blanket peatland located in Scotland, UK. Results demonstrate the high spatio-temporal variability in peatland C cycling (that is, R eco , P cal , NEE, CH 4 fluxes, and pore water [DOC]) with season, dominant PFT, and locations across the peatland (Table 1). Our first hypothesis, that PFT would explain more of the variance in GHG fluxes and [DOC] than the integrated effect of location within the peatland (site), was supported for R eco and P cal (Table 1). However, site explained more variation in CH 4 fluxes and DOC than PFT (Table 1). Our second hypothesis that interactions explained a significant proportion of the variance in GHG fluxes and [DOC] was supported (Tables 1, 2).

The Importance of Season, PFT, and the Integrated Effect of Spatially Varying Properties
The northern latitude of Black Law wind farm promoted a strong seasonal cycle, with both temperature and water table dominated by seasonal effects (Tables 1, 3). Consequently, CO 2 fluxes and [DOC] showed strong seasonality: the lowest CO 2 fluxes were observed during the winter months, when low temperatures suppressed photosynthesis and decomposition processes, and the highest fluxes occurred during summer, when biological processes were promoted by warmer temperatures (Clark and others 2005;Bonnett and others 2006). Pore water [DOC] did not peak until autumn, attributable to the relatively high temperatures in summer and autumn that promoted decomposition and mobilization into the dissolved phase in autumn when the water table rose (Tables 4, 5) (Grayson and Holden 2012). CH 4 also peaked in autumn, which coincided with relatively high water tables and temperature (including at depth) (Tables 4,5). The relatively high autumnal water tables and temperatures may have promoted greater methanogen activity (CH 4 production) and limited methanotroph activity (CH 4 oxidation), leading to the higher CH 4 fluxes (Moore and others 2011;Dijkstra and others 2012;van Winden and others 2012). Further, long CH 4 residence times may have promoted the accumulation of a large CH 4 pool by autumn (Kellner and others 2005;Strack and others 2008).
PFT which, as hypothesized, dominated over the integrated effect of site for CO 2 fluxes (but not CH 4 or [DOC]), has been shown to influence peatland C cycling due to the direct and indirect effects of different traits (for example, physiological, phenological, and ecological) others 2009, 2013). In this study, there were significant effects of PFT on CO 2 fluxes with R eco and P cal rates lower in moss-dominated plots than from shrub and sedge plots (Table 5). This is unsurprising due to their lower biomass, less efficient photosynthesis, lower litter quality, and the production of secondary metabolites which are known to inhibit decomposition (Turetsky 2003;Marschall and Proctor 2004;Toet and others 2006). However, water table depth may also be playing a role: depth to the water table was significantly greater at shrub plots compared with moss plots (Table 4). The R eco and P cal rates of moss plots were less responsive to seasonal change compared with sedges and shrubs, as mosses are active throughout the year, whereas vascular plants have a defined growing season ( Figure 1). The dominance of PFT indicates the pivotal role of plant communities in influencing peatland C balances and the potential implications of both natural and managed vegetation change.
PFT explained the second largest amount of variation in the CH 4 fluxes (site dominated), with higher fluxes associated with sedges (Tables 1, 5). Sedges possess aerenchyma (a morphological trait that shrubs and mosses do not possess), which act as a conduit between the soil and atmosphere. Aerenchyma prevent oxidation of CH 4 , which typically occurs in aerobic layers of peat, although they can also transport oxygen into the anoxic zone, and thus promote CH 4 oxidation in the rhizosphere (Levy and others 2012;Bhullar and others 2013). Further, differences in root exudate quality (for example, organic acid concentrations) between the different PFTs may contribute to the higher CH 4 fluxes from sedges (Strö m and others 2012). As the water table depth was not significantly different under sedges compared with mos-ses and shrubs (Table 4), we infer that oxidation potential within the peat profile did not play a role in the differences in CH 4 fluxes.
There was no relationship between PFT and [DOC], despite the influence of plant functional traits on physical and biological conditions that regulate [DOC] (Armstrong and others 2012). At Black Law the vegetation was a very patchy mosaic of shrubs, sedges, and mosses, and thus the physical and biological effects of the PFTs at a sampling point would be less apparent compared with other peatlands at which the vegetation occurred in larger patches and [DOC] was found to be significantly correlated with PFT (Armstrong and others 2012).
The integrated effect of site was significant for all GHG fluxes and [DOC], although less dominant than season and PFT for R eco and P cal , and less dominant than season for [DOC] (Table 1). Within site differences can be more important than temporal (inter-annual) controls on C cycling in peatlands and thus understanding of the underlying mechanisms is crucial (Waddington and Roulet 2000). Spatial differences in C cycling will be influenced by variation in the driving biotic and abiotic properties. At Black Law, both relatively static (bulk density, peat depth, C, N, and C:N) and dynamic properties (water table and biomass) varied significantly with site (Tables 3, 4). Further, site will have incorporated other unmeasured static and dynamic spatially variable characteristics (for example, microbial community metrics (Schmidt and others 2007)), known to influence C cycling. Resolving the role and interaction between these static and dynamic, biotic, and abiotic characteristics is necessary to infer the effects of environmental change on peatland C cycling, especially for [DOC] and CH 4 fluxes.

The Importance of Interactions
The statistical modeling demonstrated that biotic and abiotic factors had individual and interactive effects on peatland C cycling and that their relative importance varied among R eco , P cal , CH 4 fluxes, and [DOC] ( Table 2). Although these models indicated that the effect sizes of the individual drivers were generally greater than those of the interactions, together the interactions explained more variance in P cal , NEE, and CH 4 than individual drivers, and the variation in R eco and [DOC] explained by the interactions was not inconsequential (Table 2). Consequently, our second hypothesis, that interactions would explain a significant proportion of the variance in GHG fluxes and [DOC] was sup-ported. Our analysis also highlighted thresholds in interactions. For example, P cal decreased with increasing shrub biomass when sedge biomass was high, was constant when sedge biomass was average, and increased when sedge biomass was low. Such thresholds further indicate the complexity of C cycling in peatlands and the need to develop understanding of interactive effects.
Interactions with temperature were key for CO 2 fluxes ( Table 2), suggesting that there is potential for the response to global warming to be exacerbated, or mediated, by the concomitant variation in biotic and abiotic conditions. The interactions with temperature may also explain why temperature sensitivities are highly variable between studies/ treatments (Davidson and Janssens 2006). There were synergistic effects between temperature and shrubs for R eco and temperature and sedge for P cal ( Table 2). The synergistic effect between temperature and shrubs is concerning with regard to peat C storage given the observed and projected increase in both in some systems (Sturm and others 2001). Further, the interactions between temperature and water table, peat depth, shrub and moss biomass for R eco , and temperature and peat depth and C:N for P cal , suggest that the sensitivity of R eco and P cal to these factors may increase under a warmer climate and that response between peatlands and different sites within a peatland may diverge.
For all GHG fluxes, there was either an interaction between shrub and sedge biomass or shrub and moss biomass, with the effect sizes, as indicated by the coefficient size, among the largest for the interactions (Table 2). This indicates the key role that PFT plays in peatland C cycling and also indicates the importance of plant community composition. The centrality of PFT may reflect direct effects, such as photosynthesis rates given leaf area, or indirect effects, such as different decomposition rates due to litter and root exudation characteristics and microbial communities associated with each PFT (Crow and Wieder 2005; Moore and others 2007; Bardgett and van der Putten 2014). The biomass of the PFTs also interacted with a range of abiotic factors, including climatic and peat property variables (Table 2). Previously, plant removal experiments have shown that warming effects on R eco did not differ between PFTs (Ward and others 2013). In contrast, we found that the R eco response to temperature increased at higher shrub and lower moss biomass, although our analysis is based on seasonal temperature changes, whereas Ward and others had a warming treatment ( Table 2). The interactions between temperature and shrub and moss biomass likely indicate the indirect effects of PFT on C cycling. For example, the weaker response of R eco to increasing temperature at high moss biomass may reflect the inhibiting role of secondary metabolites and low N, both of which are associated with Sphagnum (Turetsky and others 2008). Further, microbial taxa are associated with specific PFTs and their response to temperature, and thus decomposition rates, may vary (Bardgett and van der Putten 2014).
The importance and nature of the interactions indicates that under future environmental change, components of the system, which together determine the C balance of the peatland, may respond in different ways. PFT biomass interacted with both biotic and abiotic factors, to affect GHG fluxes and [DOC] demonstrating the centrality of vegetation to peatland C cycling. CO 2 fluxes were dominated by interactions including temperature, thus signaling an immediate response to climate warming. In contrast, interactions including peat physical and chemical properties were more important for CH 4 and [DOC] suggesting that immediate effects could occur in response to land management that affects peat physical and chemical properties (including fertilization and compaction through vehicle access). In addition there may also be longer term effects on CH 4 fluxes and [DOC] in response to indirect feedbacks to environmental change (including change in litter quality). These varying responses to drivers that change over different time scales and interact make projection of the composite effect of environmental change on peatland C cycling complex.

CONCLUSIONS
Resolving the relative role of biotic and abiotic factors, including interactions, on globally important northern peatlands is central to predicting their contribution to mediating or exacerbating climate change. The number and strength of the interactions in the models used here and the presence of thresholds indicates the need for multivariate experiments within which the roles of interactions are explicitly tested. Without understanding the interactions between and within biotic and abiotic factors, it is difficult to predict the effects of future change on peatland C cycling, or translate results between sites with different characteristics. The interactions highlighted in this field experiment, within which biotic and abiotic factors varied naturally, suggest future research should be focussed on temperature interactions for CO 2 fluxes, bulk density, and peat C concentration for CH 4 and [DOC], and PFT for all components. Such studies will include spatially and temporally extensive assessment of C cycling in field environments, but also controlled laboratory studies examining the interactions highlighted in this study.