# Douglas fir stimulates nitrification in French forest soils

Bernhard Zeller, Arnaud Legout, Séverine Bienaimé, Bruno Gratia, Philippe Santenoise, Pascal Bonnaud, Jacques Ranger

<figure><img src="image_1.png" /></figure>

## To cite this version:

Bernhard Zeller, Arnaud Legout, Séverine Bienaimé, Bruno Gratia, Philippe Santenoise, et al. Douglas fir stimulates nitrification in French forest soils. Scientific Reports, 2019, 9 (1), pp.1-11. 10.1038/s41598-019-47042-6. [doi]

HAL Id: hal-02264369

https://hal.science/hal-02264369

Submitted on 6 Aug 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Douglas fir stimulates nitrification in French forest soils

Bernd Zeller<fnref n="1" />, Arnaud Legout<fnref n="1" />, Séverine Bienaimé<fnref n="1" />, Bruno Gratia<fnref n="1" />,², Philippe Santenoise<fnref n="1" />, Pascal Bonnaud<fnref n="1" /> & Jacques Ranger<fnref n="1" />

Douglas fir trees presumably stimulate nitrification in the soil. We studied in 21 French Douglas fir forests if and how nitrification is modulated by soil properties, past land use and current forest management. Soil (0–10 cm depth) was collected and initial concentrations of N-NH₄⁺ and N-NO₃⁻, potential net nitrogen mineralization (PNM) and net nitrification (PNN) rates and microbial biomass were measured. At 11 of the 21 sites, annual nitrate fluxes in the soil were measured using anion exchange resin bags. Soils contained between 2.3 to 29.4 mg N-NO₃⁻ kg soil⁻<fnref n="1" />. About 86% (±14%) of mineral N was nitrate. The proportion of nitrate increased to almost 100% during incubation. PNN varied from 0.10 mg N kg soil⁻<fnref n="1" /> day⁻<fnref n="1" /> to 1.05 mg N kg soil⁻<fnref n="1" /> day⁻<fnref n="1" /> (21 sites). Neither the initial nitrate concentration nor PNN was related to soil chemistry (pH, % C, %N, P, CEC), microbial biomass, texture, past land use or thinning. *In situ* net nitrate accumulation (NNA) estimated with resins beds varied from 4 to 100 kg N-NO₃⁻ ha⁻<fnref n="1" /> yr⁻<fnref n="1" /> (11 sites). It was positively correlated with base saturation, clay content, ELLENBERG N, temperature and negatively with soil organic N, C/N ratio and precipitation.

Received: 26 December 2018

Accepted: 5 July 2019

Published online: 23 July 2019

Trees are an integral part of nitrogen cycling in forests because trees take up nitrogen, grow, produce litter which decompose and releases nitrogen. Tree species influence nitrogen cycling through uptake of mineral and/or organic N, litter quality and quantity and interactions with microbes at the plant root interface<fnref n="1" />⁻³. Although, availability of mineral nitrogen is still a factor limiting the productivity of forest ecosystems in many biomes⁴⁻⁶. Especially, tree species with a high biomass increment and thus, a high demand for nitrogen and other nutrients, as Douglas fir (*Pseudotsuga menziesii* (Mirb.) Franco var. menziesii)) need to apply/develop efficient strategies to match these requirements.

Since the late 50's until the 80's significant areas in France were afforested with Douglas fir while during the same period broadleaf forests were converted into Douglas fir forests. Under optimal conditions annual increment of Douglas fir may reach 30 m³ ha⁻<fnref n="1" /> yr⁻<fnref n="1" /> which rank this species among the most productive coniferous species in Europe⁷,⁸. Since its discovery by the Scottish botanist David Douglas in 1824 this forest tree species had been introduced in Europe and contributes now to about 3% or 421000 ha of the French forest cover⁷,⁹,¹⁰. The most commonly cited reasons for this success are (i) rapid initial growth and high production levels¹¹ over a wide range of environments through effective coping strategies¹², (ii) efficient nutrient uptake strategies and biomass production¹³, (iii) absence of major pests and pathogens¹⁰ and at last (iv) wood products with rheological qualities that do not deteriorate with the radial growth rate¹⁴.

The impact of Douglas fir on biogeochemical cycles have been extensively studied since the 1960's in its native range¹⁵⁻¹⁸ or when introduced as alien species in Europe¹⁹. In Europe, some studies focused on its impact on the nitrogen cycle²⁰⁻²³. Among them, Jussy *et al.* at Vauxrenard (France) and Zeller *et al.* at Breuil (France), highlighted the dominance of nitrate as main mineral N form in the topsoil when this species was planted on previous farmland or broadleaf forest²⁴,²⁵. At both sites, soils were acid, nutrient poor with low amounts of plant available Ca and Mg. Such soils are representative for about half of the area of French Douglas fir forests²⁶. In addition, in a common garden experiment with 6 tree species, a fairly rapid and significant increase of the nitrate concentration in soil cores transferred 18 months earlier from 3 low nitrifying tree species into a Douglas fir stand happened²⁷. Based on this experiment Andrianarisoa *et al.* suggested that Douglas fir triggers nitrification through interactions with nitrifying bacteria at the root - soil interface²⁸. Although the precise mechanism(s) by which Douglas fir may stimulate nitrification and enhance soil nitrate concentrations remain still unknown, as well as a possible modulation by environmental factors (e.g. soil properties, past land use, current forest management).

<footnote n="1">¹INRA Grand-EST Nancy, UR 1138 Biogéochimie des Ecosystèmes Forestiers, Route d'Amance, 54280, Champenoux, France. ²Office National de Forêts (ONF), Unité Territorial de Darney-Bains, 6 rue des Rochottes, 88260, Darney, France. Correspondence and requests for materials should be addressed to B.Z. (email: bernhard.zeller@inra.fr)</footnote>
<table><caption>**Table 1.** Soil type and soil chemical characteristics of the studied sites (soils 0–10 cm). Ca²⁺, Mg²⁺, K⁺ = exchangeable Cations (mg g⁻¹), CEC = cation exchange capacity (meq. 100 g⁻¹), S/T = base saturation (%), Al³⁺ (g kg⁻¹) and phosphorus as P₂O₅ (g kg⁻¹). In bold and *, sites (n = 11) where soil humidity, soil solution chemistry and nutrient fluxes are measured (solutions are collected once per month).</caption><thead><tr><th>Site</th><th>Soil type</th><th>C (%)</th><th>N (%)</th><th>pH</th><th>Sand (%)</th><th>Silt (%)</th><th>Clay (%)</th><th>Ca<sup>2+</sup> (mg g<sup>-1</sup>)</th><th>Mg<sup>2+</sup></th><th>K<sup>+</sup></th><th>CEC (meq. 100 g<sup>-1</sup>)</th><th>S/T</th><th>Al<sup>3+</sup></th><th>P<sub>2</sub>O<sub>5</sub></th></tr></thead><tbody><tr><td>ADI</td><td>Hyperdystric Cambisol</td><td>3.45</td><td>0.23</td><td>4.45</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ALL*</td><td>Andosols</td><td>10.10</td><td>0.72</td><td>4.67</td><td>34.9</td><td>46.0</td><td>19.1</td><td>3.30</td><td>0.58</td><td>0.43</td><td>9.47</td><td>53.1</td><td>17.3</td><td>1.13</td></tr><tr><td>BOU</td><td>Alfisols</td><td>3.43</td><td>0.17</td><td>3.86</td><td>23.5</td><td>63.2</td><td>13.3</td><td>1.16</td><td>0.31</td><td>0.17</td><td>5.97</td><td>35.9</td><td>1.2</td><td>0.09</td></tr><tr><td>CONT*</td><td>Eutric Cambisol</td><td>2.80</td><td>0.20</td><td>4.36</td><td>8.1</td><td>74.3</td><td>17.6</td><td>2.98</td><td>0.53</td><td>0.51</td><td>7.98</td><td>58</td><td>1.7</td><td>1.10</td></tr><tr><td>DOU34</td><td>Dystric Cambisol</td><td>4.01</td><td>0.33</td><td>4.56</td><td>22.9</td><td>34.9</td><td>31.2</td><td>2.46</td><td>0.43</td><td>0.16</td><td>5.36</td><td>56.9</td><td>24.0</td><td>0.03</td></tr><tr><td>GAI*</td><td>Dystric Cambisol</td><td>4.01</td><td>0.27</td><td>3.94</td><td>10.7</td><td>70.8</td><td>19.9</td><td>0.75</td><td>0.25</td><td>0.35</td><td>8.38</td><td>23.9</td><td>2.3</td><td>0.39</td></tr><tr><td>MEY*</td><td>Entic Podzols</td><td>9.40</td><td>0.66</td><td>3.96</td><td>52.2</td><td>23.7</td><td>24.1</td><td>0.09</td><td>0.13</td><td>0.27</td><td>9.90</td><td>10.7</td><td>3.9</td><td>0.21</td></tr><tr><td>MIG</td><td>Dystric Planosol</td><td>3.85</td><td>0.25</td><td>5.24</td><td>21.4</td><td>53.7</td><td>24.9</td><td>1.06</td><td>0.19</td><td>0.32</td><td>5.62</td><td>35.2</td><td>2.3</td><td>0.11</td></tr><tr><td>MONT*</td><td>Eutric Planosol</td><td>2.46</td><td>0.22</td><td>4.44</td><td>13.8</td><td>60.6</td><td>25.6</td><td>2.06</td><td>0.32</td><td>0.41</td><td>7.30</td><td>57.3</td><td>6.1</td><td>0.63</td></tr><tr><td>OLL*</td><td>Hyperdystric Cambisol</td><td>3.39</td><td>0.29</td><td>4.12</td><td>55.9</td><td>25.1</td><td>19.0</td><td>0.48</td><td>0.13</td><td>0.17</td><td>6.35</td><td>15.5</td><td>2.9</td><td>0.27</td></tr><tr><td>ORL-T*</td><td>Entic Podzols</td><td>10.98</td><td>0.74</td><td>4.13</td><td>62.0</td><td>21.7</td><td>16.3</td><td>0.09</td><td>0.13</td><td>0.28</td><td>8.02</td><td>9.9</td><td>6.5</td><td>0.20</td></tr><tr><td>ORL-F*</td><td>Entic Podzols</td><td>7.69</td><td>0.48</td><td>4.32</td><td>76.1</td><td>12.8</td><td>12.6</td><td>0.28</td><td>0.12</td><td>0.16</td><td>5.86</td><td>4.6</td><td>6.9</td><td>0.15</td></tr><tr><td>PAU*</td><td>Hyperdystric Cambisol</td><td>4.70</td><td>0.37</td><td>4.20</td><td>52.5</td><td>25.5</td><td>22.0</td><td>1.11</td><td>0.31</td><td>0.32</td><td>9.63</td><td>27.2</td><td>4.1</td><td>1.15</td></tr><tr><td>PEY-T</td><td>Entic Podzols</td><td>5.30</td><td>0.36</td><td>4.09</td><td>61.0</td><td>18.2</td><td>20.8</td><td>0.42</td><td>0.21</td><td>0.34</td><td>11.52</td><td>12.8</td><td>4.3</td><td>0.12</td></tr><tr><td>PEY-ESC</td><td>Entic Podzols</td><td>5.52</td><td>0.36</td><td>4.00</td><td>70.4</td><td>15.4</td><td>14.2</td><td>0.21</td><td>0.74</td><td>0.23</td><td>5.10</td><td>13.3</td><td>4.2</td><td>0.15</td></tr><tr><td>PIQ</td><td>Hyperdystric Cambisol</td><td>6.74</td><td>0.49</td><td>4.18</td><td>63.3</td><td>20.4</td><td>16.6</td><td>0.26</td><td>0.75</td><td>0.15</td><td>4.50</td><td>13.5</td><td>7.8</td><td>0.46</td></tr><tr><td>ROQ*</td><td>Eutric Cambisol</td><td>4.31</td><td>0.37</td><td>4.20</td><td>25.1</td><td>47.9</td><td>27.0</td><td>7.62</td><td>1.16</td><td>0.49</td><td>14.03</td><td>70.0</td><td>1.1</td><td>0.13</td></tr><tr><td>SAL-B</td><td>Calcaric Cambisol</td><td>10.08</td><td>0.70</td><td>6.71</td><td>24.7</td><td>16.8</td><td>58.5</td><td>32.2</td><td>4.14</td><td>0.86</td><td>40.21</td><td>96.4</td><td>2.9</td><td>0.67</td></tr><tr><td>SAL-H</td><td>Andosols</td><td>5.41</td><td>0.38</td><td>4.81</td><td>17.6</td><td>52.6</td><td>29.8</td><td>8.57</td><td>2.72</td><td>0.7</td><td>12.90</td><td>94.4</td><td>3.8</td><td>0.62</td></tr><tr><td>VIL*</td><td>Eutric Cambisol</td><td>2.58</td><td>0.19</td><td>4.04</td><td>14.1</td><td>54.6</td><td>31.3</td><td>1.72</td><td>0.34</td><td>0.39</td><td>8.92</td><td>35.7</td><td>4.4</td><td>0.14</td></tr><tr><td>VULC</td><td>Andosols</td><td>6.42</td><td>0.47</td><td>5.35</td><td>35.8</td><td>43.1</td><td>21.1</td><td>1.07</td><td>0.19</td><td>0.17</td><td>3.99</td><td>39.5</td><td>24.1</td><td>0.63</td></tr></tbody></table>

Consequences of high soil nitrate concentrations and resulting leaching of nitrate in Douglas fir forests is loss of cations (Al, Ca, Mg) and soil acidification¹⁹,²⁹. However, this impact depends on soil conditions e.g. under acidic conditions when nitrate is not immobilized probably in relation to deficient microbial activity or root uptake, a massive leaching of nitrate and monomeric Al³⁺ occurred in deep drainage solutions³⁰. Such a phenomenon potentially affects the quality of surface waters³¹.

Up to now, almost all studies about nitrate concentrations and potential nitrification were obtained in Douglas fir forests on acid, more or less nutrient poor soils. Thus, the question arises if this tree species affects soil nitrate concentrations and potential nitrification in a generalized way in less acid and more fertile soils. Our hypothesis is that in Douglas fir forest soils in France, nitrate is the dominant mineral N form and that potential nitrification rates were modulated by soil properties. In the light of this hypothesis our objectives were to (1) evaluate the amount and form of mineral N in the mineral soil (0–10 cm depth) in 21 Douglas fir forests in France, distributed all over the distribution area of this species, and assess net mineralization and nitrification potentials under controlled conditions, (2) compare these potentials to their expression under field conditions using ion-exchange resins and (3) tackle how environmental factors (e.g. climate, soil fertility, stand structure, tree growth and past land use) may modulate the availability of mineral nitrogen in the soil of Douglas fir forests.

## Results

**Soil chemical status of the experimental sites.** Amongst the studied sites, soil types and their soil physico-chemical characteristics varied broadly. Cambisols formed the biggest group, the other soil types like Podzols, Alfisols, Andosols or Planosols are less abundant (1–2 sites per soil type) (Table 1). In line with the soil type, soil pH, concentrations of soil organic C and soil organic N differed between the sites but remained in a range commonly observed in French forest soils. As for the soil type, the texture of the soils covered a large range, from sandy to loamy until clayey soils. Most other variables as, CEC, exchangeable cations like Al and Ca showed a similar pattern than the former parameters thus validating the selection of each site. Indeed, sites were thought to represent the whole panel of French soils now covered by Douglas fir forests.

**Initial mineral N content.** In the 84 soil samples (0–10 cm layer), collected at the 21 sites, the soil pH covered a range from 3.8 to 6.7, the C concentration varied from 2.5% C to 10.9% C, the N concentration varied from 0.2% N to 0.7% N and the C/N varied from 11.0 to 20.1 (Table 1). In the field fresh soils, immediately extracted after their sampling in the Douglas fir forests, the ammonium concentration varied from 0.7 to 12.8 mg N-NH₄⁺ kg soil⁻¹ and the nitrate concentration varied from 2.3 to 29.4 mg N-NO₃⁻ kg soil⁻¹ (Table 2). In all soils, the percentage of nitrate varied from 56% to 93% with 86% (±14%) as a mean. Nitrate was always much more abundant than ammonium at all sites.
<table><caption>**Table 2.** Initial nitrate concentration at T = 0 and percentage of nitrate in field fresh soils (n = 4), potential net nitrification rate (PNN) and percentage of nitrate in incubated soils (n = 4), net nitrate accumulation rate (NNA) measured with resin beds in the field (n = 30). Mean values and standard variation. In bold and *, sites (n = 11) where soil humidity, soil solution chemistry and nutrient fluxes are measured.</caption><thead><tr><th>Site</th><th>T = 0 (mg N-NO<sub>3</sub><sup>-</sup> kg soil)</th><th>% (N-NO<sub>3</sub><sup>-</sup>)</th><th>PNN (mg N-NO<sub>3</sub><sup>-</sup> kg soil<sup>-1</sup> day<sup>-1</sup>)</th><th>% (N-NO<sub>3</sub><sup>-</sup>)</th><th>NNA (kg N-NO<sub>3</sub><sup>-</sup> ha<sup>-1</sup> yr<sup>-1</sup>)</th></tr></thead><tbody><tr><td>ADI</td><td>22,18 (± 10,25)</td><td>88</td><td>0,37 (± 0,22)</td><td>95</td><td></td></tr><tr><td>ALL*</td><td>13,36 (± 6,27)</td><td>93</td><td>1,05 (± 0,45)</td><td>100</td><td>25,79 (± 5,12)</td></tr><tr><td>BOU</td><td>5,78 (± 1,90)</td><td>82</td><td>0,37 (± 0,14)</td><td>99</td><td></td></tr><tr><td>CONT*</td><td>19,99 (± 6,16)</td><td>99</td><td>0,69 (± 0,25)</td><td>100</td><td>36,39 (± 8,13)</td></tr><tr><td>DOU34</td><td>7,90 (± 1,81)</td><td>95</td><td>0,62 (± 0,29)</td><td>98</td><td></td></tr><tr><td>GAI*</td><td>10,30 (± 0,80)</td><td>75</td><td>0,63 (± 0,21)</td><td>98</td><td>42,09 (± 11,20)</td></tr><tr><td>MEY*</td><td>4,03 (± 1,37)</td><td>62</td><td>0,19 (± 0,11)</td><td>53</td><td>4,30 (± 2,27)</td></tr><tr><td>MIG</td><td>10,95 (± 1,87)</td><td>85</td><td>0,83 (± 0,10)</td><td>91</td><td></td></tr><tr><td>MONT*</td><td>29,44 (± 12,44)</td><td>99</td><td>0,44 (± 0,12)</td><td>99</td><td>78,18 (± 14,55)</td></tr><tr><td>OLL*</td><td>3,10 (± 1,59)</td><td>80</td><td>0,38 (± 0,11)</td><td>94</td><td>28,87 (± 10,38)</td></tr><tr><td>ORL-F*</td><td>2,31 (± 0,91)</td><td>80</td><td>0,64 (± 0,13)</td><td>96</td><td>10,53 (± 3,00)</td></tr><tr><td>ORL-T*</td><td>7,39 (± 1,76)</td><td>78</td><td>0,81 (± 0,07)</td><td>99</td><td>18,53 (± 4,95)</td></tr><tr><td>PAU*</td><td>6,75 (± 3,22)</td><td>94</td><td>0,60 (± 0,20)</td><td>99</td><td>37,86 (± 12,05)</td></tr><tr><td>PEY-F</td><td>10,60 (± 6,25)</td><td>96</td><td>0,58 (± 0,15)</td><td>98</td><td></td></tr><tr><td>PEY-T</td><td>11,93 (± 1,48)</td><td>91</td><td>0,53 (± 0,19)</td><td>99</td><td></td></tr><tr><td>PIQ</td><td>12,93 (± 2,02)</td><td>92</td><td>0,47 (± 0,04)</td><td>98</td><td></td></tr><tr><td>ROQ*</td><td>10,30 (± 8,13)</td><td>82</td><td>0,29 (± 0,05)</td><td>99</td><td>101,26 (± 23,24)</td></tr><tr><td>SAL-B</td><td>4,07 (± 0,73)</td><td>53</td><td>0,10 (± 0,15)</td><td>83</td><td></td></tr><tr><td>SAL-H</td><td>27,32 (± 12,81)</td><td>95</td><td>0,68 (± 0,12)</td><td>100</td><td></td></tr><tr><td>VII*</td><td>10,66 (± 5,38)</td><td>84</td><td>0,34 (± 0,12)</td><td>98</td><td>63,15 (± 12,24)</td></tr><tr><td>VULC</td><td>19,05 (± 7,46)</td><td>98</td><td>0,66 (± 0,27)</td><td>99</td><td></td></tr></tbody></table>

**Potential net mineralization (PNM) and potential net nitrification (PNN) rate.** The daily PNN rate varied by a factor of 10 among the sites, from 0.10 mg N-NO₃⁻ kg⁻¹ day⁻¹ (SAL-B) till 1.05 mg N-NO₃⁻ kg⁻¹ day⁻¹ (ALL) (Table 2). After three weeks and at the end of the incubation period of six weeks, about 90–98% of mineral N was present as nitrate, except at one site (53% at MEY). We observed a highly significant linear relationship between the rates of PNM and PNN (r² = 0.97, P < 0.05) (Fig. 1A). The few values, observed below the 1:1 line, all correspond to the same site (MEY) with the lowest percentage of nitrate. Both, PNM and PNN varied among the 21 studied sites, but PNN was neither related to the soil organic N concentration, the pH, the base saturation (S/T ratio) or the clay content (Fig. 1B, C, E, F). Neither a trend nor any relationship appeared between the PNN and other site-specific variables like former land use and current forest management (Fig. 1G, H). Obviously, the species effect on PNN rates is much stronger than site-specific parameters as soil fertility or past land use. Exclusion of the data from the two sites with the lowest PNN and the lowest or highest pH (MEY and SAL-B) revealed a positive correlation between PNN and soil microbial biomass C and N (r² = 0.43, P < 0.05) (Fig. 1D).

**In situ net nitrate accumulation (resin beds).** As for PNN, a large variation of nitrate accumulation in resins was observed among the 11 sites (Table 2). Average annual net nitrate accumulation (NNA) varied by a factor of 25 among the sites, from 4.3 kg N-NO₃⁻ ha⁻¹ yr⁻¹ at MEY to 101.3 kg N-NO₃⁻ ha⁻¹ yr⁻¹ at ROQ. Nevertheless, beside these two extremes, NNA varied in a narrower range among the remaining 9 sites (10.5 kg N-NO₃⁻ ha⁻¹ yr⁻¹ to 78.2 kg N-NO₃⁻ ha⁻¹ yr⁻¹). NNA was positively correlated to the base saturation (S/T ratio, r² = 0.69; P < 0.05) (Fig. 2A), the clay content (r² = 0.44; P < 0.05) (Fig. 2B), the Ellenberg N (r² = 0.53; P < 0.05, Fig. 2E) and the number of days with a temperature >5 °C (r² = 0.32; P < 0.05, Fig. 2H). NNA was negatively correlated to soil organic N (r² = 0.51; P < 0.05) (Fig. 2C), the C/N ratio of the mineral soil 0–10 cm (r² = 0.45; P < 0.05) (Fig. 2D) and mean annual precipitation (r² = 0.31; P < 0.05) (Fig. 2G). NNA was higher in former coniferous forests than in former broadleaf forests and former farmland or pasture (Fig. 2F). Excluding the data from the site MEY, with the lowest PNN and NNA, revealed a negative relationship between NNA and PNN (r² = 0.58; P < 0.05) (Fig. 2I).

**Multivariate analysis.** A Principal Component Analysis (PCA) has been used to summarize the potential relationship between site characteristics (soil and stand properties) and net nitrification (PNN) or net nitrate accumulation rate (NNA) (Fig. 3, Table S1). Base saturation (S/T), sum of base cations (S), silt and C<sub>mic</sub>/N<sub>mic</sub> are opposed to C/N<sub>horA</sub>, N<sub>mic</sub> and soil organic nitrogen (SON) on the principal component 1 (PC1). C/N<sub>horO</sub>, duration of the growing season (GW) are opposed to aluminum extracted by the Tamm method (AlTamm), pH, soil organic nitrogen (SON) and phosphorus (P) on the principal component 2 (PC2). The principal plan PC1-PC2 did not allow to discriminate the PNN groups, whereas at least the high NNA group differentiates from the others NNA groups. These findings corroborate previous results presented in Figs 1, 2, that emphasized strong relationships
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Relationships between potential net nitrification (PNN) and (A) potential net mineralization (PNM), (B) pH, (C) soil organic N, (D) microbial biomass N (except MEY and SAL_B), (E) clay content, (F) base saturation (S/T), (G) past land use (1 = conifers, 2 = broadleaf, 3 = pasture or farmland) and (H) forest management (years since last thinning).</figcaption></figure>

between net nitrate accumulation rate and some site characteristics, whereas no clear relationship was observed between net nitrification and site characteristics.

**Review of PNN and PNM in coniferous, broadleaf and Douglas fir forest soils.** In the selected studies (n = 9), PNM and PNN rates were obtained by incubating mineral soils under fairly similar incubation conditions than ours (temperature, soil water content). Over the panel of published data, large variations in PNM and PNN among coniferous, broadleaf and Douglas fir forests were observed (Fig. 4). In Douglas fir forests, there is a strong positive correlation between PNM and PNN (r² = 0.95, p < 0.05) which follows the 1:1 line, indicating that PNN = PNM. A similar positive correlation is observed for coniferous forests (r² = 0.69, p < 0.05). Compared to Douglas fir, PNN is lower than PNM in coniferous forests, but nitrate is still the dominant mineral N form. Although a large variability, in broadleaf forests, a positive relation between PNN and PNM is observed (r² = 0.32, p < 0.05). Compared to Douglas fir and coniferous forests, PNN is much lower than PNM in broadleaf forests. Mean rates of PNM are 0.54 mg N kg⁻¹ day⁻¹, 0.48 mg N kg⁻¹ day⁻¹ and 0.84 mg N kg⁻¹ day⁻¹ in Douglas fir forests, coniferous forests and broadleaf forest, respectively. Mean rates of PNN are 0.55 mg N-NO₃⁻ kg⁻¹ day⁻¹, 0.40 mg N-NO₃⁻ kg⁻¹ day⁻¹ and 0.48 mg N-NO₃⁻ kg⁻¹ day⁻¹ in Douglas fir forests, coniferous forests and broadleaf forest, respectively.

Discussion

**Dominance of nitrate in the topsoil.** Our hypothesis that nitrification is predominant in French Douglas fir forests soils is supported by the initial ammonium and nitrate concentrations measured in the topsoil at the 21 sites.

Firstly, although large inter-site variations of the initial soil nitrate concentration (3–30 mg N kg⁻¹) were observed, we show that a dominant share of the mineral N in the 0–10 cm soil (on average 86% ± 14%) was present as nitrate. As this percentage was neither related to soil properties nor forest management, we suppose that the main reason for such an imbalance between ammonium and nitrate is the tree species itself. Climate events (drought or enhanced soil temperature) may reinforce or weaken the amount of residual soil nitrate but this factor
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Relationships between annual net nitrate accumulation (NNA) and (A) base saturation (S/T), (B) clay content, (C) soil organic N, (D) C/N of the 0–10 cm soil), (E) Ellenberg N, (F) past land use, (G) mean annual precipitation, (H) number of days with T > 5°C and (I) PNN (except Mey). Data obtained at 11 sites.</figcaption></figure>

will not basically modify this trait. In addition, the average autumnal nitrate concentration in Douglas fir (12 mg N·NO₃⁻·kg soil⁻¹) is definitively at the upper boundary compared to nitrate concentrations measured in other temperate forests soils²⁵,³²,³³ and the Ellenberg N values of the forest floor vegetation confirms the long-lasting dominance of nitrate in the studied Douglas fir forests. Usually, the leftover of nitrate in late autumn is a sign that the production of nitrate exceeds the uptake by the trees and the understory vegetation³⁴. Our results are consistent with earlier findings at the Breuil common garden experiment²⁷,³⁰,³⁵, where apparently higher soil and soil solution nitrate concentrations were measured in Douglas fir compared to other deciduous or coniferous species. In addition, Douglas fir increased rapidly the soil nitrate concentration in soil cores transferred into the Douglas fir plots from other less nitrifying plots whereas the opposite effect was observed in soil cores from Douglas fir transferred into less nitrifying plots²⁷.

Secondly, we show that the potential net nitrification (PNN) is high at all sites, except one and is probably only limited by the bottleneck's depolymerization of soil organic N and ammonification. These findings illustrate the high potential of nitrifying bacteria to convert ammonium into nitrate and corroborate the dominance of nitrate in French Douglas fir soils. We may suppose that under non-limiting conditions, the potential nitrification rate would be even higher than observed. Positive or negative effects of different tree species on net and gross N mineralization and nitrification had been observed in a certain number of common gardens experiments³⁶⁻⁴⁰. In the above-mentioned studies, net mineralization was close to net nitrification in Douglas fir soils, a factor that approves the dominance of nitrate under this species. Nevertheless, old growth Douglas fir forests in the US did not really fit into this schema, because net nitrification is here much lower⁴¹. In those older forests the buildup of thicker humus layers indicates modifications of the soil chemistry and reduction of biological activity with a switch towards a higher abundance of fungi. For instance, increase in forest floor thickness along a chronosequence of beech stands resulted in a sharp decrease in nitrification⁴².

**Modulation of soil nitrate contents in Douglas fir forests.** Litter (needle and roots) is the most important source of recent soil organic matter in the studied Douglas fir forests and the thin organic layers in all studied stands indicated a highly efficient team constituted by fauna, fungi and bacteria that decomposes litter thus creating a continuous flux of soil organic matter into the topsoil. Some authors also reported a positive effect of Douglas fir on the humus C:N ratio, soil pH and base saturation¹⁰,²⁹. Nevertheless, the absence of any relation between PNN and soil organic C and N suggest that soil organic matter is probably not a principal driver of nitrification at our sites. Moreover, the potential net mineralization and net nitrification measured at the 21 sites
<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) distribution of the 17 parameters (soil chemistry, stand characteristics) (b) distribution of the sites according to the PNN group, (c) distribution of the sites according to the NNA group.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Relationships between PNM and PNN in French Douglas fir forests soils (red dots), in coniferous forests soils (black dots) (White *et al.*, 1988; Gower and Son 1992; Menyailo *et al.*, 2002; Nugroho *et al.*, 2006; Colman and Schimel 2013; Urukawa *et al.*, 2016) and in broadleaf forest soils (grey dots) (Morris and Boerner 1998; Verchot *et al.*, 2001; Gower and Son; Menyailo *et al.*, 2002; Andrianarisoa *et al.*, 2009; Colman and Schimel 2013; Urukawa *et al.*, 2016)</figcaption></figure>

cope fairly little with soil chemical parameters. These findings concur with the results of two meta-analysis with numerous tree species<sup>43,44</sup> where a minor part of the observed variability of net mineralization and nitrification was explained by climatic and soil chemical parameters (33% by MAP + %C + %N + clay content) but both studies supposed that biological parameters like soil organic matter quality and microbial community composition strongly increase this percentage.

Nitrification had also been related to the tree species<sup>39,45,46</sup>, past land use<sup>32</sup> and forest management practices<sup>18,29,47-49</sup>. In the present study, recognized drivers of nitrification in forest ecosystems like past land use and
forest management practices did not mark net nitrification and nitrate accumulation in resins in the Douglas fir stands. Thus, other intrinsic parameters seem responsible for the predominance of nitrate in French Douglas fir soils. Different kinds of interactions at different levels (molecular, diversity) between the tree species and microbes and/or soil parameters are appropriate to affect soil N transformations and nitrification⁵⁰,⁵¹.

Some tree species like Norway spruce or birch produce monoterpenes or diterpenes as agents to down regulate nitrifying bacteria⁵²⁻⁵⁴ while the absence of monoterpenes was observed in Douglas fir stands⁵². Moreover, Nugroho *et al.* suggested that Douglas fir could even up regulate the number of nitrifying bacteria: although the mechanism is not well understood, a positive correlation between AOB's like *Nitrosospira* and nitrate production was observed³⁹. Andrianarisoa *et al.* revealed that interactions at the root-soil-interface leading to a rapid up regulation of the number of ammonium oxidizing bacteria (AOB) in soils cores transferred from two low nitrifying species (Norway spruce, Nordmann fir) into a Douglas fir plot²⁸. According to these studies there is some evidence that Douglas fir stimulate AOB's, but how is still an open question. Nugroho *et al.* supposed that tree species regulate AOB's through their effects on the soil C/N ratio, but Zeller *et al.* observed a clear difference between different species in net and gross nitrification in a common garden experiment where soils had a similar C/N ratio²⁵,³⁹. As previously detailed, neither the humus nor the soil C/N ratio was a primary determinant of PNN or NNA in our study. As at our sites climate and soil chemistry varied largely, the high PNN and soil nitrate concentrations suggests a straight stimulating effect of Douglas fir on nitrifying bacteria³⁹. In addition, Douglas fir is a neophyte in Europe, young trees were often inoculated with north American mycorrhizal strains as Laccaria laccata S238N⁵⁵. Probably, Douglas fir fine roots are not familiar to most of the ectomycorrhizal fungi (ECM) present in European forest soils which may affect the interactions between microorganisms and fine roots. In the light of such a scenario it could be possible that with ongoing time the ectomycorrhizal diversity increases thus inducing a switch towards a more conservative N cycling.

**Ex situ versus in situ approach and implication for forest management.** Ion exchange resins (IER) have been recognized as a relevant tool for the assessment of the availability of mineral nitrogen in soils⁵⁶. IER have the big advantage to investigate under real field conditions (temperature, moisture, soil structure, biotic factors) soil mineral nitrogen fluxes with a very limited initial soil disturbance⁵⁷. The availability of mineral nitrogen in the nitrate form, as measured by the resins bags is the expression of the difference between the production by nitrifying bacteria, the immobilization by microorganisms and the vegetation and finally loss by leaching of excess nitrate. The negative correlation between NNA and PNN reflects the complexity of the N cycle in Douglas fir forests and underlines that the potential is not expressed *in situ*. Contrary to PNN, NNA was related to the S/T ratio, the clay content and the soil organic N concentration. Interestingly, less nitrate was accumulated in soils with higher soil organic C and N contents, which is somewhat surprising and opposite to observations made in other forests³⁴. Sites characterized by higher soil fertility (S/T ratio, Ca concentration) sustained nitrification but also loss of nitrate. Such a behavior suggest/indicate that even under conditions of optimal growth (nutrients, water) for Douglas fir a substantial proportion of nitrate is probably transferred into deeper soil layers and/or lost to the ground water⁵⁸. The greater amount of clay at these sites could result in more fine porosity, which would be a suitable habitat for nitrification for much of the year. Leaching of excess nitrate will result in loss of cations, acidification and *in fine* deterioration of drinking water quality²⁹⁻³¹,⁵⁹. At least, two forest management measures could reduce nitrate leaching: mixed stands composed of beech and Douglas fir²⁹,⁶⁰ and fertilization with Ca-Mg-P for acid nutrient poor soils, because such a treatment reduced the nitrate concentration significantly in soil solutions collected at 0.6 m soil depth in a so treated Douglas fir stand³⁵. Nevertheless, further experiments are necessary to evaluate the impact of fertilization on soil nitrate fluxes in other acid contexts and along a soil fertility gradient.

## Conclusion

All results converge to the same conclusion, that namely Douglas fir, an alien tree species in Europe, has a strong encouraging/stimulating effect on nitrification and top soil N-NO₃⁻ concentrations either *in situ* and *ex situ*. This effect comes first into sight in the field fresh soil samples (collected in autumn) from all 21 sites. Potential net nitrification rates (PNN) in incubated soils but also the net nitrate accumulation (NNA) in the field corroborate the domination of nitrate in French Douglas fir soils. Net nitrate accumulation is higher in more fertile soils and at warmer sites. This suggests that even under optimal conditions for tree growth (high biomass increment) an excess of nitrate remains in the soil with a peak in autumn. As nitrate is highly mobile in the soil profile, leaching loss of nitrate and cations may affect surface and groundwater quality, as well as the sustainability of soils by an acidification process. Future studies are needed to assess the environmental impact of Douglas fir by evaluating properly the leaching losses of nitrate at the plot scale and at larger scales (watersheds etc.) in monospecific Douglas fir forests. Another perspective, which seems the most difficult, is to understand the mechanism(s) by which Douglas fir stimulates the activity of nitrifying bacteria. Lastly, further research should focus on how forest management measures may mitigate nitrate fluxes, probably by the set up of mixed stands and through liming/fertilization.

## Material and Methods

**Sites.** All over France, 21 Douglas fir forests (private and public forests) were selected from a larger number of sites (Table 3, Fig. S1). Those sites are thought to represent all areas in France where Douglas fir had been introduced since the last century and where Douglas fir covers a significant surface of the forested area. Soil properties (physical and chemical parameters) were one important criterion for the final selection of a site in order to represent the whole panel of soil types and soil fertility. Additional criteria were vegetation cover (Ellenberg N) past land use and current forest management practices⁶¹. Climatic parameters like annual rainfall, air temperature and length of the growing season were recorded but not included into the site selection process. Among these 21 sites,
a sub-network of 11 sites was set up with the objective to monitor continuously soil moisture, fluxes of nutrients and litter fall.

**Potential net N-mineralization and nitrification rates.** At each of the 21 sites, mineral soil samples (0–10 cm depth) were collected from 4 random subplots situated at 0°, 90°, 180° and 270° about 3 m around an anchor point. The more or less thin humus layer was carefully removed, then all soil (0–10 cm depth) was collected from a surface of 0.25 m², put into a big box, passed through a 4 mm sieve, mixed and finally a representative sample of about 3 kg soil per subplot was shipped in a cooled container to the laboratory. Mineral N concentrations were measured in all soil samples immediately after their arrival in the laboratory. Then, the soils were stored at 4 °C, until the soil sampling campaign was completed and all soils were available for the laboratory incubation. This short storage period (3 weeks) was necessary in order to setup and start the incubation experiment with all soil samples at the same time. Potential net mineralization and nitrification rates were obtained at 20 °C and 60% water holding capacity (WHC). To meet the later criteria for each site the water holding capacity was measured in all soils before the set-up of the incubations and, if necessary, deionized water was added. Soil from each of the four subplots per site (300 g) were put into 2L Mason jars, pre - incubated during 10 days, followed by an incubation period of 42 days. During the whole period, jars were opened twice a week for about 15 minutes, the soil mixed with a spoon to assist the evacuation of accumulated CO2. Subsamples were collected from each jar at the end of the pre - incubation period (T0) and after 21 days (T1) and 42 days (T2) of incubation. Each sample of 40 g of soil was extracted with 200 ml of 0.5 M K₂SO₄, shaken during 1 hour, filtered and finally the concentrations of N-NH₄⁺ and N-NO₃⁻ in the extracts were measured using continuous flow analysis (SKALAR San + +). Potential net N mineralization (PNM = Nmin final - Nmin initial) was the amount of total inorganic N accumulated during the incubation period and PNN (PNN = N-NO₃⁻ final - N-NO₃⁻ initial) was the amount of N-NO₃⁻ formed from the nitrification of the NH₄⁺ already present in the soil core before the incubation as well as that mineralized during the incubation. Both were calculated as mg N kg⁻¹ soil d⁻¹. All concentrations and rates are presented on a dry weight basis. Microbial biomass C and N was measured after 3 weeks of incubation using the fumigation extraction method. Concentration of soil organic C and N were measured using a CHN analyzer (Carlo Erba). In addition, basic soil parameters like soil texture, pH, concentration of P and the CEC were measured at the soil analysis laboratory (LAS) of INRA in Arras.

**In situ accumulation of ammonium and nitrate (resin bed method).** In order to evaluate the *in-situ* N mineralization and nitrification we set up a batch of ion exchange resin traps (anionic and cationic) at 10 cm soil depth and measured over one year at 4, 8, and 12 months the amounts of ammonium and nitrate captured by the resins.

DOWEX 21 K (MERCK- CEC = 1.95 meq g⁻¹ dry resin, 20–50 mesh) were used to catch anions (nitrate and others) Anionic resins were first rinsed under a slow stream of distilled water and then saturated with Cl⁻ from 1 M HCl (hydrochloric acid) at a ratio of 1000 ml HCl per 100 g of dry resin. Then, 20 grams of wet resin representing about 36–40 meq for anionic were introduced into traps made of PVC rings (Ø 50 mm and 10 mm in height), closed, from both sides by a nylon mesh (20 µm mean mesh size), allowing free percolation of water. At each of the 11 sub-sites, 6 sets of 10 resin traps were set up (total = 660 traps). After opening a pit of about 3 m long, 0,4 m large and 0,15 m deep, the traps were carefully introduced horizontally into the undisturbed soil (at a depth of 10 cm and with a minimum distance of 25 cm between two traps). Each trap was marked with dyed rope and the trench carefully refilled with the excavated soil, slight compacted and finally the organic layer was restored. In total 6 independent pits were set up per site and two pits were sampled every four months between June 2012 and June 2013.

Each resin trap was first rinsed three times with distilled water (10 ml applied with a syringe) for eliminating soil particles. Then each resin trap was weighed and stored at 4 °C until extraction. N-NO₃⁻ was extracted with 1M NaCl at a ratio of 4 g resin per 40 ml of NaCl, after a one-hour batch contact on a rotary shaker. Concentration of N-NO₃⁻ in the desorbed extracts was measured by colorimetry using a continuous flow analyzer (SKALAR San + +). Net nitrate accumulation rate (NNA) measured by resins is the amount of NO₃⁻-N fixed by the resins during the time of incubation in the field, divided by the number of days between installation and collection. We computed only the mean annual rates of NNA. Net nitrate accumulation was expressed as kg N-NO₃⁻ ha⁻¹ y⁻¹.

**Statistics.** Single and multiple regression analysis were used to detect significant relationship among potential net nitrification and soil and forest management parameters. The same approach was used for the net nitrate accumulation rate. A standardized Principal Component Analysis (sPCA) has been used to analyze the data collected from the 21 sites. This data comprised 17 quantitative variables based on soil physical and chemical properties (*pH*, *silt*, *Clay*, *SOC*, *SON*, *Nmic*, *Cmic*, *Nmic*, *S*, *S:T*, *AlTamm*, *P*, *C:NHorO*, *C:NHorA*) and site characteristics (*Ellen*, *GW*, *G*, *TSLT*). The global idea of sPCA method is to project the standardized data onto a lower dimensional linear space generated by uncorrelated variables, named principal components (PC), such as the variance of the projected data is maximized. The first step is to select a minimal number of PCs while keeping a large part of the explained variance in the new PC subspace. For this, we decided (1) to retain a number of PCs explaining at least 50% of the total variance and (2) to stop introducing new PCs when the fraction of the explained variance decrease sharply. Once the PC subspace considered, the second step is to detect outliers into this subspace. For the outlier detection, we are calculated the Hotelling's T² statistic (Jackson *et al.* 1991, Villegas *et al.* 2010). The T² statistic is distributed according to a F distribution and the upper 95% confidence limit of this statistic is determined as the maximal accept limit. In case where outliers are detected, a second application of the sPCA method is required after deleting these ones. The retained PCs are interpreted from the correlation circle and the contribution of individual quantitative variables. Moreover, the quality of representation (squared cosine) is calculated
<table><caption>**Table 3.** Description of the study sites (mean annual precipitation, altitude, age and number of the trees, past land use, G = basal area, SV = days with a daily mean temperature > 5,8 °C, Ellenberg N = indicator of N availability). Sites in bold and with * belongs to the network of sub-sites with continuous monitoring of soil solution chemistry, soil humidity and litterfall.</caption><thead><tr><th>Site</th><th>Code</th><th>Precipitation (mm per year)</th><th>Altitude (m)</th><th>SV (days)</th><th>Age (years)</th><th>Trees per ha</th><th>G (m² ha<sup>-1</sup>)</th><th>Ellenberg (N)</th><th>Past land use</th></tr></thead><tbody><tr><td>Adinfert</td><td>ADI</td><td>681</td><td>117</td><td>208</td><td>37</td><td>374</td><td>34</td><td>4.2</td><td>Coniferous Forest (CF)</td></tr><tr><td>Allagnat*</td><td>ALL*</td><td>917</td><td>1005</td><td>133</td><td>68</td><td>211</td><td>53</td><td>6.0</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Bouillancourt</td><td>BOU</td><td>780</td><td>146</td><td>211</td><td>43</td><td>368</td><td>28</td><td>4.3</td><td>Coniferous Forest (CF)</td></tr><tr><td>Wailly*</td><td>CONT*</td><td>706</td><td>111</td><td>208</td><td>42</td><td>204</td><td>36</td><td>6.1</td><td>Coniferous Forest (CF)</td></tr><tr><td>Verreries de Moussans</td><td>DOU34</td><td>1150</td><td>685</td><td>221</td><td>41</td><td>416</td><td>46</td><td>5.6</td><td>Coniferous Forest (CF)</td></tr><tr><td>Gaillefontaine*</td><td>GAI*</td><td>888</td><td>232</td><td>208</td><td>32</td><td>230</td><td>39</td><td>5.4</td><td>Farmland (pasture or field) (FL)</td></tr><tr><td>Meymac*</td><td>MEY*</td><td>1125</td><td>916</td><td>152</td><td>51</td><td>200</td><td>18</td><td>2.7</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Azy-le-Vif</td><td>MIG</td><td>730</td><td>302</td><td>152</td><td>38</td><td>278</td><td>14</td><td>5.9</td><td>Coniferous Forest (CF)</td></tr><tr><td>Montapas*</td><td>MONT*</td><td>996</td><td>288</td><td>153</td><td>46</td><td>472</td><td>36</td><td>6.9</td><td>Coniferous Forest (CF)</td></tr><tr><td>Olliergues*</td><td>OLL*</td><td>1030</td><td>643</td><td>185</td><td>31</td><td>450</td><td>37</td><td>6.2</td><td>Farmland (pasture or field) (FL)</td></tr><tr><td>Pérols sur Vézère*</td><td>ORL-T*</td><td>1326</td><td>840</td><td>162</td><td>48</td><td>429</td><td>46</td><td>2.5</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Pérols sur Vézère*</td><td>ORL-F*</td><td>1326</td><td>841</td><td>162</td><td>48</td><td>425</td><td>52</td><td>4.5</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Anglès*</td><td>PAU*</td><td>1424</td><td>802</td><td>176</td><td>41</td><td>300</td><td>39</td><td>7.0</td><td>Farmland (pasture or field) (FL)</td></tr><tr><td>Saint Amans le Petit</td><td>PEY-T</td><td>1250</td><td>586</td><td>172</td><td>43</td><td>1333</td><td>53</td><td>4.5</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Saint Amans le Petit</td><td>PEY-ESC</td><td>1250</td><td>586</td><td>172</td><td>43</td><td>285</td><td>40</td><td>4.5</td><td>Deciduous forest or heather (DF)</td></tr><tr><td>Le Bez</td><td>PIQ</td><td>1000</td><td>649</td><td>184</td><td>73</td><td>217</td><td>51</td><td>6.2</td><td>Coniferous Forest (CF)</td></tr><tr><td>Montreuil en Auge*</td><td>ROQ*</td><td>788</td><td>79</td><td>222</td><td>40</td><td>304</td><td>42</td><td>5.6</td><td>Coniferous Forest (CF)</td></tr><tr><td>Sallèdes</td><td>SAL-B</td><td>680</td><td>743</td><td>201</td><td>42</td><td>388</td><td>36</td><td>4.8</td><td>Coniferous Forest (CF)</td></tr><tr><td>Sallèdes</td><td>SAL-H</td><td>720</td><td>817</td><td>196</td><td>42</td><td>353</td><td>33</td><td>4.2</td><td>Coniferous Forest (CF)</td></tr><tr><td>Felleries*</td><td>VIL*</td><td>822</td><td>195</td><td>213</td><td>40</td><td>338</td><td>41</td><td>4.6</td><td>Coniferous Forest (CF)</td></tr><tr><td>Ceyssat</td><td>VULC</td><td>900</td><td>959</td><td>134</td><td>71</td><td>201</td><td>53</td><td>5.7</td><td>Deciduous forest or heather (DF)</td></tr></tbody></table>

for each variable in order to evaluate which variables are well represented onto different principal plans. Thus, a variable with a contribution superior to the uniform contribution and a high quality of representation has a large impact for the considered PC. The dispersion of sites into principal plan is then illustrated for the variables PNN and NNA transformed in qualitative variables: the variable PNN was divided into 4 groups (low <0.4; medium [0.4; 0.6]; high [0.6; 0.7]; very high ≥0.7) while the variable NNA was divided in 3 groups (low <30; medium [30; 50]; high ≥50).

## References

1. Knops, J. M. H., Bradley, K. L. & Wedin, D. A. Mechanisms of plant species impacts on ecosystem nitrogen cycling. *Ecology Letters* **5**, 454–466 (2002).

2. Laughlin, D. C., Hart, S. C., Kaye, J. P. & Moore, M. M. Evidence for indirect effects of plant diversity and composition on net nitrification. *Plant and Soil* **330**, 435–445 (2010).

3. Meier, I. C., Finzi, A. C. & Phillips, R. P. Root exudates increase N availability by stimulating microbial turnover of fast-cycling N pools. *Soil Biology and Biochemistry* **106**, 119–128 (2017).

4. Le Bauer, D. S. & Treseder, K. K. Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. *Ecology* **89**, 371–379 (2008).

5. Perakis, S. S. & Sinkhorn, E. R. Biogeochemistry of a temperate forest nitrogen gradient. *Ecology* **92**, 1481–1491 (2011).

6. Rennenberg, H. & Dannenmann, M. Nitrogen nutrition of trees in temperate forests - the significance of nitrogen availability in the pedosphere and atmosphere. *Forests* **6**, 2820–2835 (2015).

7. Hermann, R. K. & Lavender, D. P. Douglas-fir planted forests. *New Forests* **17**, 53–70 (1999).

8. Waring, R. H. Why is the productivity of Douglas-fir higher in New Zealand than in its native range in the Pacific Northwest, USA? *Forest Ecology and Management* **255**, 4040–4046 (2008).

9. Debazac, E. F. Manuel des conifères, 2nd Edition (Editions ENGREF 1991).

10. Schmid, M., Pautasso, M. & Holdenrieder, O. Ecological consequences of Douglas fir (Pseudotsuga menziesii) cultivation in Europe. *European Journal of Forest Research* **133**, 13–29 (2014).

11. Oswald, H. & Parde, J. 1984 Le Douglas en France: Sylviculture et Production. *Revue Forestière Française* **26**, 56–68 (1984).

12. Du, B. *et al*. Foliar nitrogen metabolism of adult Douglas-fir trees is affected by soil water availability and varies little among provenances. *PloS one* **13**, e0194684 (2018).

13. Ranger, J. *et al*. Nutrient budgets for a rotation of Douglas-fir plantation in the Beaujolais (France) based on a chronosequence study. *Forest Ecology and Management* **171**, 3–16 (2002).
14. Nepveu, G. & Blachon, J. L. Largeur de cernes et aptitudes à l'usage en structure de quelques confières: Douglas - fir, Pin sylvestre, Pin maritime, Epicéa de Sitka, Epicéa commun et Sapin pectiné. *Revue Forestière Française* **6**, 497-506 (1989).

15. Wright, E. & Bollen, W. B. Microflora of Douglas-fir forest soil. *Ecology* **42**, 825-828 (1961).

16. Van Migroet, H. & Cole, D. W. Acidification source in red alder and Douglas - fir soils - Importance of nitrification. *Soil Science Society of America Journal* **49**, 1274-1279 (1985).

17. Horner, J. D., Cates, R. G. & Gosz, J. R. Tannin, nitrogen, and cell wall composition of green vs. senescent Douglas-fir foliage. *Oecologia* **72**, 515-519 (1987).

18. White, C. S., Gosz, J. R., Horner, J. D. & Moore, D. I. Seasonal, annual, and treatment-induced variation in available nitrogen pools and nitrogen-cycling processes in soils of two Douglas-fir stands. *Biology and Fertility of Soils* **6**, 93-99 (1988).

19. Marques, R. & Ranger, J. Nutrient dynamics in a chronosequence of Douglas-fir (*Pseudotsuga menziesii* (Mirb.) Franco) stands on the Beaujolais Mounts (France). 1: Qualitative approach. *Forest Ecology and Management* **91**(2-3), 255-277 (1997).

20. Gijsman, A. J. Nitrogen nutrition of Douglas - fir (*Pseudotsuga menziesii*) on strongly acid sandy soil. 1. growth, nutrient uptake and ionic balance. *Plant and Soil* **126**, 53-61 (1990a).

21. Gijsman, A. J. Nitrogen nutrition of Douglas-fir (*Pseudotsuga menziesii*) on strongly acid sandy soil. 2. Proton excretion and rhizosphere pH. *Plant and Soil* **126**, 63-70 (1990b).

22. Tietema, A., De Boer, W., Riemer, L. & Verstraten, J. M. Nitrate production in nitrogen-satured acid forest soils: vertical distribution and characteristics. *Soil Biology and Biochemistry* **24**, 235-240 (1992).

23. Tietema, A. et al. Nitrogen cycling in acid forest soils subject to increased atmospheric nitrogen input. *Forest Ecology and Management* **57**, 29-44 (1993).

24. Jussy, J. H., Colin-Belgrand, M. & Ranger, J. Production and root uptake of mineral nitrogen in a chronosequence of Douglas-fir (*Pseudotsuga menziesii*) in the Beaujolais Mounts. *Forest Ecology and Management* **128**, 197-209 (2000).

25. Zeller, B. et al. Influence of tree species on gross and net N transformation in forest soils. *Annals of Forest Science* **67**, 51-158 (2007).

26. https://inventaire-forestier.ign.fr.

27. Andrianarisoa, S. et al. Dynamic control of nitrification by tree species. *Ecosystems* **13**, 1171-1187 (2010).

28. Andrianarisoa, S., Zeller, B., Siegenfuhr, H., Bienaimé, S. & Dambrine, E. Root colonization and nitrate availability in forest soils. *Soil Use & Management* **33**, 45-53 (2017).

29. Prietzel, J. & Bachmann, S. Changes in soil organic C and N stocks after forest transformation from Norway spruce and Scots pine into Douglas fir, Douglas fir/spruce, or European beech stands at different sites in Southern Germany. *Forest Ecology and Management* **269**, 134-148 (2012).

30. Legout, A., Van Der Heijden, G., Jaffrain, J., Boudot, J.-P. & Ranger, J. Tree species effects on solution chemistry and major element fluxes: A case study in the Morvan (Breuil, France). *Forest Ecology and Management* **378**, 244-258 (2016).

31. Gauthier, A. et al. Dynamics of the water extractable organic carbon pool during mineralization in soils from a Douglas fir plantation and an oak-beech forest - an incubation experiment. *Plant and Soil* **330**, 465-479 (2010).

32. Jussy, J. H., Koerner, W., Dambrine, E., Dupouey, J. L. & Benoit, M. Influence of former agricultural land use on net nitrate production in forests soils. *European Journal of Soil Sciences* **53**, 367-374 (2004).

33. Persson, T. et al. Soil nitrogen turnover - mineralisation, nitrification and denitrification in European forest soils in *Carbon and nitrogen cycling in European forest ecosystems* (ed. Schulze, E. D.) 297-311 (Springer, 2000).

34. Hackl, E., Zechmeister-Boltenstern, S., Bodrossy, L. & Sessitsch, A. Comparison of diversities and compositions of bacterial populations inhabiting natural forest soils. *Applied and Environmental Microbiology* **70**, 5057-5065 (2004).

35. Ranger, J., Bonnaud, P., Bouriaud, O., Gelhaye, D. & Picard, J. F. Effects of the clear-cutting of a Douglas-fir (*Pseudotsuga menziesii* (Mirb.) Franco) plantation on chemical soil fertility. *Annals of Forest Science* **65**, 1-13 (2008).

36. Gower, S. T. & Son, Y. Differences in soil and leaf litterfall nitrogen dynamics for five forest plantations. *Soil Science Society of America Journal* **56**, 1959-1966 (1992).

37. Menyailo, O. V., Hungate, B. A. & Zech, W. The effect of single tree species on soil microbial activities related to C and N cycling in the Siberian artificial afforestation experiment. *Plant and Soil* **242**, 183-196 (2002).

38. Malchair, S. & Carnol, M. Microbial biomass and C and N transformations in forest floors under European beech, sessile oak, Norway spruce and Douglas-fir at four temperate forest sites. *Soil Biology and Biochemistry* **41**, 831-839 (2009).

39. Nugroho, R. H., Röling, W. E. M., Laverman, A. M. & Verhoef, H. A. Net nitrification rate and presence of Nitrospira cluster 2 in acid coniferous forest soils appear to be tree species specific. *Soil Biology and Biochemistry* **35**, 1166-1171 (2006).

40. Trum, F., Titeux, H., Ranger, J. & Delvaux, B. Influence of tree species on carbon and nitrogen transformation patterns in forest floor profiles. *Annals of Forest Science* **68**, 837-847 (2011).

41. Hart, S. C., Nason, G., Myrold, D. D. & Perry, D. A. Dynamics of gross nitrogen transformations in an old-growth forest: The carbon connection. *Ecology* **75**, 880-891 (194).

42. Trap, J., Bureau, F., Vinceslas-Akpa, M., Chevalier, R. & Aubert, M. Changes in soil N mineralization and nitrification pathways along a mixed forest chronosequence. *Forest Ecology and Management* **258**, 1284-1292 (2009).

43. Colman, B. P. & Schimel, J. P. Drivers of microbial respiration and net N mineralization at the continental scale. *Soil Biology and Biochemistry* **60**, 65-76 (2013).

44. Urakawa, R. et al. Factors contributing to soil nitrogen mineralization and nitrification rates of forest soils in the Japanese archipelago. *Forest Ecology and Management* **361**, 382-396 (2016).

45. De Boer, W. & Kowalchuk, G. A. Nitrification in acid soils: micro-organisms and mechanisms. *Soil Biology and Biochemistry* **33**, 853-866 (2001).

46. Laughlin, D. C. Nitrification is linked to dominant leaf traits rather than functional diversity. *Journal of Ecology* **99**, 1091-1099 (2011).

47. Prietzel, J., Wagoner, G. L. & Harrison, R. B. Long-term effects of repeated urea fertilization in Douglas-fir stands on forest floor nitrogen pools and nitrogen mineralization. *Forest Ecology and Management* **193**, 413-426 (2004).

48. Strahm, B. D. et al. Soil solution nitrogen concentrations and leaching rates as influenced by organic matter retention on a highly productive Douglas-fir site. *Forest Ecology and Management* **218**, 74-88 (2005).

49. Bai, S. H. et al. Effects of forest thinning on soil-plant carbon and nitrogen dynamics. *Plant and Soil* **411**, 437-449 (2017).

50. Steinlein, T. Invasive alien plants and their effects on native microbial soil communities in *Progress in Botany* (eds Lüttge, U., Beyschlag, W., Francis, D., Cushman, J.) 123-125 (Springer 2013).

51. Prescott, C. E. & Grayston, S. J. Tree species influence on microbial communities in litter and soil: current knowledge and research needs. *Forest Ecology and Management* **309**, 19-27 (2013).

52. White, C. S. Monoterpenes: their effects on ecosystem nutrient cycling. *Journal of Chemical Ecology* **20**, 1381-1406 (1994).

53. Priha, O. & Smolander, A. Nitrogen transformations in soil under *Pinus sylvestris*, *Picea abies* and *Betula pendula* at two forest sites. *Soil Biology and Biochemistry* **31**, 965-977 (1999).

54. Adamczyk, S., Kiikkilä, O., Kitunen, V. & Smolander, A. Potential response of soil processes to diterpenes, triterpenes and tannins: Nitrification, growth of microorganisms and precipitation of proteins. *Applied Soil Ecology* **67**, 47-52 (2013).

55. Selosse, M. A., Jacquot, D., Bouchard, D., Martin, F. & Le Tacon, F. Temporal persistence and spatial distribution of an American inoculant strain of the ectomycorrhizal basidiomycete *Laccaria bicolor* in a French forest plantation. *Molecular Ecology* **7**, 561-573 (1998).
56. Binkley, D. & Matson, P. Ion Exchange resin bag method for assessing forest soil nitrogen availability. *Soil Science Society of America Journal* **47**, 1050–1052 (1983).

57. Carlyle, J. C. & Malcolm, D. C. The use of ion exchange resin bags to assess N availability beneath pure spruce and larch + spruce stands growing on a deep peat soil. *Plant and Soil* **93**, 123–127 (1986).

58. Perakis, S. S. et al. Coupled nitrogen and calcium cycles in forests of the Oregon Coast Range. *Ecosystems* **9**, 63–74 (2006).

59. Mareschal, L., Turpault, M. P., Bonnaud, P. & Ranger, J. Relationship between the weathering of clay minerals and the nitrification rate: a rapid tree species effect. *Biogeochemistry* **112**, 293–309 (2013).

60. Thurm, E. A. & Pretzsch, H. Improved productivity and modified tree morphology of mixed versus pure stands of European beech (Fagus sylvatica) and Douglas-fir (Pseudotsuga menziesii) with increasing precipitation and age. *Annals of Forest Science* **73**, 1047–1061 (2016).

61. Ellenberg, H. 1986 Vegetation Mitteleuropas mit den Alpen. 4th Edition (ULMER 1986).

## Acknowledgements

We gratefully acknowledge Laurent Bisonnier, Bernard Catry, Jean François Sinet, Noémie Havet, David Daupin, Louis Gilet, Christel Joseph, Marc Lafaye, Magali Maviel, Yann Mozziconacci, Bruno Vanstaël for sample collection and assistance in the forest. Marie Clothilde Maintenant assisted at the soil sampling Tourde-France, sample analysis and prepared the resin traps. Gilles Nourrisson contributed to the set-up of the resin traps in the forest. This project received financial support from ONF – RDI and the MEAG company. The UR BEF is supported by the French National Research Agency through the Cluster of Excellence ARBRE (ANR-11-LABX-0002-01). The experiment was supported by the mobile lab (M-POETE) of ANAEE-France. ANAEE-France is an infrastructure from the French Investment for the Future (Investissements d'Avenir) program, overseen by the French National Research Agency (ANR-11-INBS-0001).

## Author Contributions

B.Z., J.R., B.G., A.L. designed the study. B.Z., S.B., P.B. contributed toward conducting the research. B.Z., J.R., A.L., P.S. analyzed the data. All authors contributed toward composing the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-019-47042-6.

**Competing Interests:** The authors declare no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019