Temperature sensitivity patterns of carbon and nitrogen processes in decomposition of boreal organic soils – Quantification in different compounds and molecule sizes based on a multifactorial experiment

Climate warming and organic matter decomposition are connected in a recursive manner; this recursion can be described by temperature sensitivity. We conducted a multifactorial laboratory experiment to quantify the temperature sensitivity of organic carbon (C) and nitrogen (N) decomposition processes of common boreal organic soils. We incubated 36 mor and 36 slightly decomposed Carex-Sphagnum peat samples in a constant moisture and ambient temperature for 6 months. The experiment included three temperature and two moisture levels and two food web manipulations (samples with and without fungivore enchytraeid worms). We determined the release of carbon dioxide (CO2) and dissolved organic carbon (DOC) in seven molecular size classes together with ammonium N and dissolved organic N in low molecular weight and high molecular weight fractions. The temperature sensitivity function Q10 was fit to the data. The C and N release rate was almost an order of magnitude higher in mor than in peat. Soil fauna increased the temperature sensitivity of C release. Soil fauna played a key role in N release; when fauna was absent in peat, the N release was ceased. The wide range of the studied C and N compounds and treatments (68 Q10 datasets) allowed us to recognize five different temperature sensitivity patterns. The most common pattern (37 out of 68) was a positive upwards temperature response, which was observed for CO2 and DOC release. A negative downward pattern was observed for extractable organic nitrogen and microbial C. Sixteen temperature sensitivity patterns represented a mixed type, where the Q10function was not applicable, as this does not allow changing the sign storage change rate with increasing or decreasing temperature. The mixed pattern was typically connected to intermediate decomposition products, where input and output fluxes with different temperature sensitivities may simultaneously change the storage. Mixed type was typical for N processes. Our results provide useful parameterization for ecosystem models that describe the feedback loop between climate warming, organic matter decomposition, and productivity of N-limited vegetation.

Introduction that the function of soil fauna responds to temperature change and is therefore indirectly linked with GHG balance. Although the importance of the interactions between soil fauna, microbes, and decomposition in GHG balance has been recognized [26,29], building these interactions into simulation models requires further studies to expose the hidden interactions between C and N processes. Model development will benefit especially from multifactorial decomposition experiments [18] that differentiate the effect of soil type, temperature, moisture, and soil food webs. Such studies are currently absent in the body of scientific literature.
The aim of this study was to quantify the temperature sensitivity of C and N release in boreal organic soils using a controlled laboratory experiment. We incubated the most common boreal organic soils, mor from upland forest sites and Carex-Sphagnum peat from drained forested peatland, and studied the release rate of CO 2 , high-molecular weight (HMW) DOC, lowmolecular weight (LMW) DOC, DOC in different relative molecule size classes, NH 4 -N, HMW-DON, and LMW-DON during the experiment. Incubation was conducted at +5˚C, +10˚C, and +20˚C and in two moisture levels. Furthermore, we incubated soil with the presence and absence of enchytraeid worms. Temperature sensitivity patterns were evaluated for each soil type and soil faunal treatment. This study continues the works presented by Laurén et al. [27] and Lappalainen et al. [5,28], who demonstrated the role of soil fauna and soil type on different C and N release rates but did not study the effect of abiotic factors, such as temperature and moisture on C and N release.

Study setup
We conducted a multifactorial laboratory experiment to study the temperature sensitivity of C and N release in decomposition of boreal organic soils. The study flowchart is presented in Fig  1. We studied storage change rates of 13 different C compounds and 11 N compounds; the abbreviations and descriptions of the studied compounds are summarized in Table 1.

Study sites and soil sampling
The boreal landscape is characterized by a mosaic of upland forests and open or forested peatlands [30]. For this study, soil samples were collected from the middle boreal region in Sotkamo, eastern Finland, from sites that represent the most common upland and peatland site types in Finland [31,32,33], namely mesic upland forest and drained forested peatland with Carex-Sphagnum peat and medium fertility. The long-term (1981-2010) mean annual precipitation in the area was 591 mm, with about 40% falling as snow; the mean annual air temperature was +2.3˚C [34]. The upland site (Kangasvaara 63˚51'N /28˚58'E, altitude 220 m a.s.l.,) is dominated by old-growth Norway spruce (Picea abies (L.) Karsten) with Scots pine (Pinus sylvestris L.) mixture and has been described in detail by Finér et al. [35]. The mor layer thickness varied from 5 to 9 cm. The peatland site (Koivupuro 63˚52'N /28˚39'E, altitude 200 m a.s.l.) with Scots pine stand developed on peat layer with a depth of 1-5 m and with a slightly decomposed peat (H3-H4 on the von Post [36] scale of decomposition) in the top 40 cm of the profile. The Koivupuro site has been described earlier by Ahtiainen and Huttunen [37] and Haahti et al. [38].
The sampling was conducted along two parallel lines (15 m apart) along which the soil samples were collected. The microsite for the sampling was a topographically even spot located at least 1 m away from the nearest tree. A cylindrical core (diameter 20 cm, height 9-15 cm) was cut from the organic soil layer to fit into a plastic container (diameter 20 cm, height 20 cm). The mean dry mass for the incubated soil cores was 253 g for mor and 279 g for peat. All living The experiment included two soil types (S m = mor, S p = peat), three temperature levels (T 5 = +5˚C, T 10 = +10˚C, and T 20 = +20˚C), two moisture levels (θ 1 = lower moisture, θ 2 = higher moisture), and two food web manipulations (W 0 = without enchytraeid worms, W 1 = with worms) with three laboratory replicates in each treatment. Therefore, the total number of samples was 72. For moist soil samples, one plastic container represented one replicate. To ensure an adequate supply of soil solution, two containers of lower-moisture mor represented one replicate. During the incubation the same volume of soil solution from both containers was collected for chemical analyses. In addition, four parallel soil samples for each soil type were collected for analyses of extractable C and N contents and the other basic soil characteristics requiring a destructive sampling. Additional mor material was also collected for extraction of the enchytraeid worms that represent the soil faunal function in the experiment. Enchytraeids were inoculated into half of the soil containers as described below. The soil was stored at +4˚C prior to the worm extraction.

Preceding soil analyses and preparation for incubation
All samples were defaunated by exposing them to a freeze-thaw cycle that was repeated twice. The freeze-thaw cycle may have caused a momentary C and N flush, which was considered in the calculation by omitting the first measurement and thereafter considering the long-term rate of C and N release. Laurén et al. [27] and Lappalainen et al. [5] described the sample pretreatment and defaunation and discussed the consequences of the defaunation on soil microbes and fauna.
The moisture level in the samples was adjusted using the following procedure: gravitational water was first allowed to drain through a hole in the bottom of the containers for all samples. Soil samples representing lower moisture conditions (treatment θ 1 ) were then allowed to dry through evaporation at +4˚C. In peat, the drying was further enhanced using suction samplers (described below). For higher moisture treatment (θ 2 ), the mor samples were irrigated with deionized water whereas the peat samples were kept in the initial water content. We were able to generate a modest moisture difference between θ 1 and θ 2 treatments (Table 2), where the higher moisture level corresponds to approximately -3 kPa and the lower to -8 kPa matric potential [39,40]. Extending the moisture difference towards drier conditions would have limited the use of suction samplers in the incubation. In the present setup, the collection of the soil solution sample took 3 days in the θ 1 treatment for mor; the time tends to increase considerably towards drier conditions. A time domain reflectometer (TDR) was used to monitor the initial moisture. The water content during the experiment was kept constant by adding water into the container until the initial mass was achieved after soil solution sampling. For the incubation, the containers were placed in a dark growth chamber (GR77, Conviron Controlled Environments Ltd.) with constant ambient temperature (+5˚C, +10˚C, and +20˚C) and relative humidity of 80%. Half of the samples were inoculated with enchytraeids, which were extracted from the additional mor material using a wet funnel method [41]. At least 50 worms per container were inoculated at the beginning of the incubation and a further 50 worms were inoculated monthly to ensure the continuous presence of enchytraeid population in the containers. The number of enchytraeids inoculated during the experiment (per container) corresponds to a total of approximately 8 000 individuals m -2 , or 0.1 g m -2 in worm dry mass. The total number of enchytraeids in the soil containers was determined at the end of the incubation. At the beginning of the experiment, the soil CN ratio, pH, and contents of extractable C and N compounds were determined by means of parallel soil samples ( Table 2). The total C and N content in the soils was determined using a CHN analyzer (LECO CHN-2000). The soil pH was measured from a suspension of soil in H 2 O (1:2 v:v). C and N compounds were extracted by 0.5 M K 2 SO 4 solution and the extracts were filtered through a 0.45-μm filter. Total extractable N (TN ex ), extractable NH 4 -N (NH 4 -N ex ), and extractable NO 3 -N (NO 3 -N ex ) were analyzed by flow injection analysis (FIA-Star 5020 Analyzer, FOSS TECATOR). Extractable organic C (OC ex ) was analyzed by a total organic carbon analyzer (TOC 5000A, Shimadzu). The extractable organic N (ON ex ) content was calculated by subtracting NH 4 -N ex and NO 3 -N ex from TN ex . Extractable compounds are assumed to include both the compounds in the solution and the compounds adsorbed on solid surfaces. The microbial biomass C (C mic ) and N (N mic ) were determined using the fumigation-extraction method [42,43,44].

Repeated measurements during the incubation
Repeated measurements included monitoring of CO 2 efflux and C and N storage in the soil solution (Fig 1, Table 1). Six measurement events were conducted with 3-to 7-week intervals during the incubation period of 179 to 185 days. The slight difference in the incubation time results from the time-consuming sampling and measuring sequence and the large number of containers in the experiment. CO 2 efflux from the soil containers was measured using a static chamber method with an infrared gas analyzer (ADC LCA-2, the ADC Bioscientific Ltd.) just before each soil solution sampling. The equipment and the procedure have been described in detail by Lappalainen et al. [5].
Soil solution was trapped from the containers using suction samplers (see [5] for procedural details). Three sampling tips (MacroRhizon with syringe, Eijkelkamp, length 9 cm, diameter 4.5 mm) were inserted into the soil and approximately 35 to 100 ml of soil solution was extracted in each sampling event.
The extracted soil solution was divided into two subsamples, one for determining C and N content and the other for measuring pH. The first subsample was further divided into two parts; one was ultrafiltered through a membrane with a nominal molecular weight limit of 1 kDa (Amicon Stirred Cell model 8400, pressure 1.5-2 bar) and one remained unfiltered. Prior to the ultrafiltration, the solution was diluted with 2 M KCl (sample:KCl of 3:1) to prevent flocculation and retention of the molecules on the membrane. In the filtration, two thirds of the load volume were allowed to pass through the membrane. The ultrafiltered fraction represents the low molecular weight fraction of the DOM.
The total dissolved N (TN sol ), NH 4 -N sol , NO 3 -N sol , and DOC were determined from both the ultrafiltered and the unfiltered samples (TOC-5000A, FIA-Star 5000 Analyzer). DON content was calculated as the difference between TN sol and mineral N (NH 4 -N sol + NO 3 -N sol ). Although NO 3 -N sol was included in the calculations, it is not presented in the figures and tables because the amount was marginal throughout the experiment. The relative molecularsize distributions of DOC were determined from soil solution samples by using size-exclusion chromatography analysis (HPLC, Agilent Technologies, USA) as described in Lappalainen et al. (2018). Seven different relative molecular-size classes were distinguished on the basis of the peaks in the chromatography results. Size classes 1 and 2 represented high-molecularweight DOC (HMW-DOC), i.e. > 1 kDa, and classes 3-7, LMW-DOC. Ultrafiltration and HPLC were used in parallel as together they provide a more precise picture of the molecularsize distribution in soil solution; HPLC provides the relative molecular-size distribution where the actual molecular size remains unknown, whereas the compounds passed through the ultrafiltration membrane are known to be < 1 kDa.

Destructive analyses at the end of incubation
At the end of the experiment, the volume and the mass of the soil in the containers were measured and the soil was cut vertically into four similar sectors. The soil volume was calculated from the height of the soil core and the diameter of the cylinder-shaped container. The mass of the first sector was determined before and after drying at 105˚C to estimate the volumetric water content. The dried sector was then further used to determine the loss on ignition (LOI) at 550˚C. The second sector was cut horizontally into 5-cm slices, from where enchytraeids were extracted using the wet funnel method. The third sector was used for analyses of OC ex , TN ex , N mic , and C mic . For lower moisture treatment, where two containers formed together one replicate, the mean over the containers represented the soil characteristics for the replicate.

Data processing and statistical analyses
We selected incubation time such that the change in C and N storage would be linear in time.
The studies by Conant et al. [45] and Hamdi et al. [46] suggest that C release rate at the beginning of incubation is highest but stabilizes after the labile C has been consumed (approximately after 50 days). Our previous experiments [5,27,28] with the same soil types show that the compound storage in the repeated measurements change linearly between 50 to 180 days of incubation. Linearity allowed us to use a single value, the rate of storage change (γ), as a response variable in the further analyses, thus facilitating a balanced comparison of different treatments and compounds. Due to linearity, the slight difference in the incubation time is not likely to affect the temperature sensitivity results.
The repeated measurements were processed to derive the rate of storage change in the following manner: for soil solution samples we kept record for the quantity of the substances that were removed in the soil solution extraction and considered that in the further data processing. We first added the removed quantities to the substance storage dataset and then standardized the datasets by dividing them with the dry mass of the organic fraction of the sample. We next omitted the first sampling occasion (time < 50 days) to avoid the effect of the C and N flush and to linearize the data. The original dataset is available as S1 Dataset. We then fitted a linear mixed effect model (Eq 1) to the remaining storage-time dataset. All statistical analyses were performed using the lme package in R version 3.5.3. The rate of change (μg g -1 day -1 ) of C and N fractions was expressed as the slope (γ ijkm ) of the line where the x-coordinate represents elapsed time and the y-coordinate represents storage [5,27]. To make the CO 2 release data directly comparable to the soil solution data, we constructed a cumulative release dataset by assuming that the observed CO 2 release continues with constant rate until the next measurement. The CO 2 dataset was thereafter treated similarly to the soil solution data. The application of the linear mixed effect model allows statistical comparison of the treatment effects and their interactions.
where Q ijkmn is the measured compound content per mass unit of soil (μg g -1 dry mass) in the soil type i (S m , S p ), moisture level j (θ 1 , θ 2 ), temperature k (T 5 , T 10 , T 20 ), and worm treatment m (W 0 , W 1 ), replicate n, and sampling event p. In the fixed part of the model, α is the intercept term for soil type i, moisture level j, temperature k, and worm treatment m. The intercept term includes worm-temperature, worm-moisture, worm-soil type, temperature-moisture, temperature-soil type, and moisture-soil type interaction terms. γ ijkm represents the release rate of the compound (μg g -1 d -1 ) and t ijkmnp is the time (d, days) from the beginning of the experiment. The release rate term includes the interactions as the intercept term does. The random component includes the intercept a ijkm , the slope c ijkmn , and the residual term e ijkmnp . The rate of change for the storages requiring destructive sampling (ON ex , NH 4 -N ex , OC ex , N mic , C mic ) was calculated as the difference between the end stage and a mean storage measured from parallel samples at the beginning of the experiment. We used the variance components application of mixed linear models (Eq 2): where Q ijkmn is the measured compound content at the end of the experiment (μg g -1 dry mass) in the soil type i, moisture level j, temperature k, worm treatment m, and replicate n. Q Mi is the mean compound content at the beginning of the experiment (μg g -1 ) for the soil type i measured from parallel samples and t is the duration of the experiment (d, days). The calculated change rate (the right-hand side of the equation) for each treatment is directly comparable to γ ijkm in Eq 1. A positive change rate refers to an increasing compound pool size. The application of this method allows evaluation of the significance of the treatment effects and their interactions.
We studied the temperature sensitivity of the C and N processes by assessing the release rates (Eqs 1 and 2) obtained in different temperatures. In cases where the response was positive or negative throughout the studied temperature range, we fitted a Q 10 function into the dataset (Eq 3) using a non-linear least-square algorithm in Python [47].
Where k T is the measured release rate at incubation temperature T, k ref is the measured release rate at reference temperature T ref (here 10˚C), and Q 10 is a fitted temperature sensitivity parameter. Q 10 function could not be used in cases where the change rate turned from positive to negative or from negative to positive within the temperature range as it does not allow crossing the x-axis. The implications of this will be discussed later.

Research permits
The study areas were located on state-owned land allocated to research and managed by Metsähallitus. Research activities in these areas are regulated in Partnership Agreement (2200/12 01 01 05 01/2018) between Metsähallitus and Natural Resources Institute Finland (Luke). The agreement allows the field activities done in this study including measurement of trees and soil sampling. No endangered or protected species were threatened in the field work.

The effect of enchytraeid worms
Enchytraeid worms increased the release rate of mineral N in both soil types and the effect was pronounced in the highest incubation temperature (Table 3, Figs 2 and 3). Worms enhanced CO 2 efflux especially at high temperatures and moisture levels. The effect of enchytraeid worms on CO 2 efflux was greater in mor than in peat ( Table 3).

The effect of temperature
Increasing temperature enhanced the release of CO 2 and DOC both in mor and in peat. Furthermore, the release of TN sol , NH 4 -N sol , and DON increased with increasing temperature in mor. Interestingly, a contrasting pattern for N dynamics was observed for peat. In peat without worms, the amount of N in soil solution decreased; the decrease rate was most pronounced at the highest temperature. Examining the change rate in various molecular-size DOC pools revealed that the release of DOC in all size classes increased with increasing temperature (Fig 4,  Table 4). However, this response was more pronounced in mor than in peat. The largest molecular-size fraction (DOC size1 ) and the smallest fraction (DOC size5-7 ) were practically absent in mor. Overall, the release of DOC in size classes 2-4 was smaller in peat than in mor. The destructive analyses showed that the microbial N and C decreased during the experiment and the decrease was most pronounced in the highest temperature and in the presence of worms ( Table 5).

The effect of soil moisture
The moisture range in the study was rather small and therefore moisture alone did not have a significant effect on the C and N release. However, at the highest temperature, the larger water content increased CO 2 release and decreased DOC release. As moisture played a minor role in the results, we combined the moisture levels (θ 1 and θ 2 ) for Figs 2-6.

Dissolved and extractable N
The change rate of NH 4 -N ex followed the pattern of NH 4 -N sol , but the absolute values of the NH 4 -N ex were generally 2 to 10 times higher than for the dissolved fraction (Figs 5 and 6, Tables 3 and 5). An interesting pattern was exposed when comparing signs of the change rate components between extractable organic nitrogen (ON ex ) and dissolved organic nitrogen (DON). In most cases, when ON ex components were positive, the respective DON component was negative and vice versa (Tables 3 and 5), indicating exchange processes between the adsorbed and dissolved pools. There was also a similar inverse connection between the components of NH 4 -N and ON ex . Table 3. Change rate (μg g -1 day -1 ) in dissolved N and C pools and CO 2 release during the incubation calculated using mixed linear effect models (Eq 1). Parameter γ represents the change rate for mor (S m ), no worms (W 0 ), lower moisture status (θ 1 ) at 5˚C temperature (T 5

Temperature sensitivity
The experiment revealed the following five different temperature sensitivity patterns: positive upwards (Q 10 > 1), positive downwards (0 < Q 10 < 1), negative upwards (0 < Q 10 < 1), negative downwards (Q 10 > 1), and mixed responses (Fig 6). A negative response was obtained for negative k ref values in Eq 3. The data allowed us to describe 68 temperature response curves (Figs 2-6). A positive upwards pattern was the most common and occurred in 37 cases (Table 6). This pattern was found for CO 2 and DOC in all treatments, NH 4 -N sol and DON in mor, and NH 4 -N sol in peat where worms were present. The second most common was the mixed response (16 cases), where the sign of the response changed from positive to negative or from negative to positive with increasing temperature. It is noteworthy that the Q 10 function cannot be applied for the mixed type. A negative response, where the release rate was always negative, was typically found for destructive analyses (NH 4 -N ex in peat without worms, C mic and ON ex in mor). C release in mor was considerably higher when worms were present (Fig 2). DOC release drastically increased when temperature rose from 10˚C to 20˚C; at the highest temperature the worms had the largest effect on DOC release. Enchytraeids increased the temperature sensitivity of C release in peat. The temperature sensitivity of total N release was higher than that of total carbon in mor. The difference in C release with and without worms emerged with increasing temperature. The new finding in our study was that mineral N release in peat was totally dependent on worms; there was no net release of N in samples without worms. On the contrary, N keeps accumulating in samples in the absence of worms. Thus, a negative temperature sensitivity emerged; immobilization increased with increasing temperature.

Temperature sensitivity patterns
This study presents a unique, multifactorial experiment consisting of a broad set of abiotic and biotic factors that control the temperature sensitivity of organic C and N decomposition processes of boreal organic soils. Traditionally, incubation studies use small volumes of sieved Table 4. Change rate (μg g -1 day -1 ) in different DOC molecular-size pools during the incubation (Eq 1). Parameter γ represents the change rate for mor (S m ), no worms (W 0 ), lower moisture status (θ 1 ) at 5˚C temperature (T 5 ). The other treatment effects and the interaction effects (from W 1 to T 20 & S p ) are expressed as a difference with respect to γ. Molecular size decreases from size class 1 to size class 7. DOC size1 and DOC size2 belong to HMW-DOC. Temperature sensitivity patterns of C and N processes in organic soil decomposition soil. In contrast, we used large, intact soil cores that are more likely to reflect processes in field conditions. There are only few studies where temperature sensitivity of N [9,17,18] or DOC processes [48] have been studied; to our knowledge this is the first study where temperature sensitivity of different molecular size fractions of DOC and DON have been quantified. The wide range of the studied C and N compounds (altogether 17 compounds in two soils and two worm treatments equaling 68 different combinations), including intermediate and end products of decomposition, allowed us to recognize different temperature sensitivity patterns ( Fig  7A and 7B). The most common pattern (37 out of 68) was a positive temperature response where the storage change rate was positive throughout the studied temperature range and where the rate increased with increasing temperature (Table 6, Fig 7A, solid red line). This pattern was observed for CO 2 and DOC release (Figs 2 and 3); Q 10 values ranged from 1.7 to 3.6. This falls within the range reported for temperate, boreal, and arctic organic soils [12,49,50]. ON ex and C mic (e.g. Fig 5B and 5E) represented a negative temperature response pattern where the storage change rate was negative throughout the studied temperature range and decreased with increasing temperature (Fig 7A, solid blue line). This emerges because higher temperature is likely to promote the biodegradation of ON ex into NH 4 -N and thereby decrease the storage of ON ex . The Q 10 value also has a meaningful interpretation in these cases; the negative sign is only a matter of perspective, since a decrease in one storage is an increase in another.

DOC
Sixteen out of 68 studied temperature sensitivity patterns represented the mixed type ( Fig  7B, dotted light blue line). In these cases, the sign of the storage change rate flipped from negative to positive or vice versa within the studied temperature range. Q 10 function is not applicable in these cases as it does not allow changing the sign. The mixed pattern was typically connected to intermediate decomposition products, where input and output fluxes with different temperature sensitivities may change the storage (Fig 7B and 7C). The mixed type was especially typical for N processes. This highlights the complex nature of N temperature response and may clarify results from several other studies that have reported variable or unclear temperature responses for N mineralization [13,14,17,18,51].

Temperature sensitivity patterns in context of kinetic theory
Conant et al. [52] identified three sets of processes that potentially affect the temperature sensitivity of decomposition: depolymerization of biochemically complex compounds; changes in microbial enzyme production while acclimating to warmer temperature; and processes that limit the availability of substrates for the decomposition process (e.g. adsorption and desorption).
Temperature sensitivity of decomposition and biodegradation processes can be described through kinetic theory as reviewed by Davidson and Janssens [53]. The basis of the theory can be expressed through temperature-driven Arrhenius kinetics combined with Michaelis-Menten enzyme kinetics. Arrhenius kinetics describes the relative reaction rate as affected by temperature, whereas Michaelis-Menten kinetics describes the reaction rate as affected by enzyme substrate availability. When the maximum reaction rate (V max ) and Michaelis-Menten constant K m (representing the substrate concentration at which the reaction rate is half of V max ) have their own temperature sensitivities [53], the Michaelis-Menten kinetics can be rewritten Release rate of CO 2 and rates of change of C and N storage in soil solution collected from mor samples. Samples (n = 36) were incubated under constant ambient temperature of 5, 10, and 20˚C with (worms 1) and without enchytraeid worms (worms 0). Q 10 values (Eq 3) and the temperature response patterns (see Fig 6) are presented for each compound. The observed compounds were total carbon (Total C), total nitrogen (TN sol ), carbon dioxide (CO 2 ), ammonium nitrogen (NH 4 -N sol ), dissolved organic carbon (DOC), dissolved organic nitrogen (DON), low-molecular-weight DOC (LMW-DOC), and lowmolecular-weight DON (LMW-DON). https://doi.org/10.1371/journal.pone.0223446.g002 Temperature sensitivity patterns of C and N processes in organic soil decomposition Temperature sensitivity patterns of C and N processes in organic soil decomposition in the form where Q 10,Vmax and Q 10,Km are the temperature sensitivities of V max and K m , respectively, and (S) is substrate concentration. Now, when the concentration exceeds K m , the overall temperature sensitivity is governed by the temperature sensitivity of V max . If, on the other hand, the substrate concentration is limiting the reaction, the overall temperature sensitivity reflects the combined effect of the temperature sensitivities of both V max and K m . It is notable that in the situation where the substrate is limiting the reaction, the affinity of enzyme to substrate decreases with increasing temperature and reaction slows down, even though V max increases. This provides an explanation for Q 10 < 1 (Fig 7A, dotted red and blue line) and results in positive downwards and negative upwards temperature sensitivity patterns. The complex temperature sensitivity patterns in N processes could be partly explained by enzyme chemistry. The activities of some N-acquiring enzymes have been found to increase strongly with temperature (protease, urease), while others (amidase) have been found to be insensitive to temperature and controlled instead by C availability [54]. Especially in N-poor peat, amidase activity can be detected without a labile C source probably because microbes are N limited [26]. Furthermore, the temperature sensitivity of some C-acquiring enzymes (betaglucosidase) have been found to respond differently to pH than N-acquiring enzymes  Table 1 for abbreviations.
https://doi.org/10.1371/journal.pone.0223446.g004 Table 5. Change rate (μg g -1 day -1 ) in microbial and extractable N and C pools during the incubation calculated using mixed random effect models (Eq 2). Parameter γ represents the change rate for mor (S m ), no worms (W 0 ), lower moisture status (θ 1 ) at 5˚C temperature (T 5  Temperature sensitivity patterns of C and N processes in organic soil decomposition (N-Acetyl-Glucosaminidase; NAGase decomposing chitin). Consequently, N release can be either less or more temperature sensitive than C release depending on pH [55], which may partly explain the different responses of mor and peat in our study and the variable responses also in previous studies.

Role of enchytraeid worms
Enchytraeids overall enhanced the release of N and CO 2 . The effect of enchytraeid worms on CO 2 efflux was greater in mor than in peat (Table 3). DOC release was also enhanced by the enchytraeids (Tables 3 and 4) and was clearest at the highest temperature ( Fig 2E). The enhanced mineralization and release of dissolved organic compounds is consistent with earlier findings [5,27,56]. Enchytraeids played a crucial role in N cycling especially in peat. In peat without worms, the N amount in soil solution decreased; the decrease rate was most pronounced at the highest temperature. This emphasizes the role of soil fauna in the nutrient release of N-poor peat soils. Enchytraeids can stimulate decomposition by processing organic matter into a more labile form [57] or by excreting N-rich compounds [58]. In the absence of worms, during the incubation N appears to be immobilized in peat, either by adsorption on soil particles or by assimilation into microbial biomass. The presence of enchytraeids increased the temperature sensitivity of C release. To our knowledge, the effect of enchytraeids on the temperature sensitivity of C and N release has not been previously reported.
Microbial N and C decreased during the experiment in mor and the decrease was most pronounced at the highest temperature (reflecting the higher loss of labile C at a higher incubation temperature [49]) and in the presence of worms (reflecting the higher grazing pressure on microbial biomass). The concomitant change in C availability and quality (loss of labile C) and microbial community (decrease in biomass and decrease especially in fungal biomass) has been previously shown in incubation studies [12].
Enchytraeids are common in forest soils and their population may fluctuate according to season [22] and due to different disturbances, such as forest management practices [59]. In this experiment we compared the presence and absence of enchytraeids. Even though the absence of enchytraeids in field conditions is unlikely, our results provide a useful estimate for the maximum range for how the fluctuating enchytraeid population may reflect C and N processes at the ecosystem level as has been demonstrated by Lappalainen et al. [5].

Mor and peat
The soluble storage change rates for C and N were almost an order of magnitude higher in mor than in peat (Figs 2 and 3), although the microbial biomass at the beginning of the experiment was almost equal in both soil types (Table 2). However, the CN ratio in the microbial biomass was drastically different (8.8 for mor and 22.8 for peat). This may reflect a higher N limitation in peat than in mor, which may slow down the biological activity in soil. The CN ratio of the decomposing material and the Q 10 values for peat were higher than for mor (Figs  2A and 3A), which is consistent with the CQT hypothesis [7]. However, the N release in peat responded less to increasing temperature in peat than in mor. DOC release increases with temperature, and smallest and largest molecule size was missing from mor, which is consistent with results shown by Lappalainen et al. [5].

Fig 5. Rate of change in K 2 SO 4 -extractable C and N pools in mor samples (n = 36).
Incubation was conducted at 5, 10, and 20˚C with (worms 1) and without enchytraeid worms (worms 0). The measured pools were ammonium (NH 4 -N ex ), microbial C, microbial N, total organic C (OC ex ), and total organic N (ON ex ). See Table 1  Temperature sensitivity patterns of C and N processes in organic soil decomposition Incubation was conducted at 5, 10, and 20˚C with (worms 1) and without enchytraeid worms (worms 0). The measured pools were ammonium (NH 4 -N ex ), microbial C, microbial N, total organic C (OC ex ), and total organic N (ON ex ). See Table 1 for abbreviations.  Temperature sensitivity patterns of C and N processes in organic soil decomposition Implications to global warming Q 10 is a robust and computationally efficient approach that enables upscaling experimental decomposition studies into the ecosystem level and beyond to a global scale. Q 10 applies easily available input, such as vegetation type and temperature scenarios, making it a suitable approach e.g. for regional, spatially explicit modeling studies. In this study, we observed higher temperature sensitivity for peat compared to mor, suggesting that peatlands are more vulnerable to climate warming than mor-covered upland soils. From a C balance point of view, the primary production of vegetation is, however, equally important and largely depends on the N mineralization from decomposition of organic material [6,10]. Our results indicate that in peatlands much less N is released per unit mass of released C. In mor, however, the temperature sensitivity for N release was higher than C release (Fig 2A and 2B). Consequently, it is likely that primary production in peatlands does not increase to the same extent as in upland sites. From a modeling point of view, the feedback loop from decomposition through N release to primary production of vegetation is essential. As pointed out by Novem Auyeueng et al. [18], ecosystem models that assume that warming will consistently increase N mineralization rates and inputs of plant-available N may overestimate the increase in terrestrial productivity and the magnitude of an important negative feedback to climate change. However, Q 10 can be made more dynamic when we account for the dependency of temperature sensitivity of N on environmental factors such as rainfall and edaphic conditions [10]. There are limitations to the Q 10 approach; thus more mechanistic models, which can account for kinetic theory, soil pH, substrate characteristics, and soil physical conditions such as hydrology and oxygen supply are needed to increase process-based understanding. However, such models tend to be computationally demanding and require input that is not available from large areas.