Research Article |
Corresponding author: Michele Freppaz ( michele.freppaz@unito.it ) Academic editor: Alessandro Campanaro
© 2019 Michele Freppaz, Davide Viglietti, Raffaella Balestrini, Michele Lonati, Nicola Colombo.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Freppaz M, Viglietti D, Balestrini R, Lonati M, Colombo N (2019) Climatic and pedoclimatic factors driving C and N dynamics in soil and surface water in the alpine tundra (NW-Italian Alps). In: Mazzocchi MG, Capotondi L, Freppaz M, Lugliè A, Campanaro A (Eds) Italian Long-Term Ecological Research for understanding ecosystem diversity and functioning. Case studies from aquatic, terrestrial and transitional domains. Nature Conservation 34: 67-90. https://doi.org/10.3897/natureconservation.34.30737
|
In alpine tundra the interannual and seasonal variability of C and N forms in soil and lake water during the short snow-free season could be significant and related to climatic and pedoclimatic variables. The hypothesis that not only the climatic and pedoclimatic parameters recorded during the summer season but also the ones measured during the previous snow-covered season could contribute to explaining the C and N dynamics in soil and surface water was tested along 10 snow-free seasons in 3 sites in the alpine tundra in the north-western Italian Alps (LTER site Istituto Mosso). Among the considered parameters, the snow cover duration (SCD) exerted a primary control on soil N-NH4+, DOC, Cmicr, Nmicr and DOC:DON ratio, with an inverse relationship. A long SCD might cause the consumption of all the subnival substrata by the soil microorganisms, determining a C starvation during the subsequent snow-free season. An opposite trend was observed for the lake water, where a longer SCD corresponded to a higher content of inorganic N forms. Among the pedoclimatic indices, the number of soil freeze/thaw cycles (FTC) recorded during the snow-covered season had a positive relation with most of soil C and N forms and N-NO3− in lake water. Only the soil DON showed an inverse pattern, and this result is consistent with the hypothesis that FTC released soil DON, subsequently decomposed and mineralized. Only N-NO3− had a significant intraseasonal variability, reaching the highest values in September both in soil and water, revealing a significant slowdown of the contribution of soil N immobilization processes.
LTER, snow cover duration, soil temperature, freeze/thaw cycles, leaching, N-NO3-
The alpine tundra is a high-mountain environment located above the tree-line, which occurs across a wide range of latitudes and landscapes with common properties such as: a short growing season, extended periods with air temperature below freezing, and long periods with snow-covered soils (
Several studies performed in other ecotones, such as the boreal forests, demonstrated that the climatic conditions of the preceding winter are important driving factors of the C cycling during the growing season (e.g.
Although organic matter decomposition may be slow, the decomposition below the snowpack may still constitute a significant proportion of the total decomposition because the snow-covered season is long (
Snowmelt- and rainfall-driven leachate of nitrate is a key hydrochemical feature in montane catchments and has been studied extensively in order to understand its underlying mechanisms and ecological consequences (
In subarctic streams and rivers, higher inorganic N concentrations have been measured during the late growing season and fall (
Although much is known in the alpine tundra about litter decomposition, soil C and N cycling and microbial communities under the snowpack and during the period of snow melting (e.g.
The research area (Long Term Ecological Research [LTER] site Angelo Mosso Scientific Institute) is located in NW Italy (Piemonte Region), close to the Monte Rosa Massif (4634 m a.s.l.), along the border with the Valle d’Aosta Region (Fig.
Localization of the study area in Italy, Cimalegna Lake basin and lake, soil sampling sites (1, 3 and 5), and the automatic weather station (AWS). Rock and soil land cover categories refer to the selected basin area of Cimalegna Lake.
All the three study sites are ascribable to the ‘Siliceous alpine and boreal grasslands’ (habitat 6150, according to the EU Habitat Directive), but a large between-site difference in plant species composition was observed according to contrasting extremes of exposure and snow cover duration. Site 1 was a typical snow-bed community belonging to the Salicion herbaceae phytosociological alliance. Site 3 was an alpine microthermal Carex curvula-dominated grassland, ascribable to the Caricion curvulae alliance. Site 5, located at the lowest altitude, was dominated by Agrostis schraderiana Bech. (Festucion variae alliance).
Among the research sites, site 1 was located within the basin of the alpine lake Cimalegna, as shown in Figure
Although the Alpine Permafrost Index Map (APIM) (
Air temperature, liquid precipitation (during the snow-free season), and snowfall have been continuously recorded since 2005 by an Automatic Weather Station (AWS) located at 2901 m a.s.l. and belonging to the Italian Army (Comando Truppe Alpine – Ufficio Meteomont) (Fig.
To assess the climatic conditions in the area, several indices were extracted from the AWS data (listed and described in Table
Indices used to assess the influence of climatic and pedoclimatic conditions on C and N forms in soil and water (adapted from
Index | Term | Definition | Unit |
---|---|---|---|
Climatic index | |||
Cumulative snowfall | CS | Cumulative daily fresh snow calculated for each hydrological year (1 October to 30 September) | cm |
Heavy precipitation days | HPD | Number of days, between samplings, when daily liquid precipitation >10 mm. For the first sampling the considered period is between melt-out day and sampling day | days |
Very heavy precipitation days | VHPD | Number of days, between samplings, when daily liquid precipitation >20 mm. For the first sampling the considered period is between melt-out day and sampling day | days |
Consecutive wet days | CWD | Maximum number of consecutive days, between samplings, when precipitation >1mm. For the first sampling the considered period is between melt-out day and sampling day | days |
Consecutive dry days | CDD | Maximum number of consecutive days, between samplings, when precipitation <1mm. For the first sampling the considered period is between melt-out day and sampling day | days |
Pedoclimatic index | |||
Snow cover duration* | SCD | Sum of “snow-covered days” in each hydrological year | days |
Melt-out day of snow* | MOD | Date of complete snowmelt (indicated as day of the year - DOY) | DOY |
Duration of soil freezing | DSF | Cumulative number of days, from October 1 to the melt-out day, when mean daily soil temperature <0 °C | days |
Soil freeze/thaw cycles | FTC | Number of soil freeze/thaw cycles in each hydrological year | number |
Mean soil temperature during soil freezing | MTF | Mean daily soil temperature when the soil is frozen, from October 1 to the melt-out day | °C |
Mean soil temperature during the snow-covered season | MTSC | Mean daily soil temperature when the soil is snow-covered | °C |
Mean soil temperature during the snow-free season | MTSF | Mean daily soil temperature between samplings. For the first sampling the considered period is between melt-out day and sampling day | °C |
Intensity of soil freezing | ISF | Minimum soil temperature when soil is frozen | °C |
Each soil study site consisted of three plots, each 9 m2. From 2008 until 2010, once a year at the end of the snow-free season (September), and from to 2011 to 2017, monthly during the snow-free season, three topsoil samples (A horizon, 0–10 cm depth) were collected, which in turn consisted of three subsamples in each plot. Soil samples were homogenized by sieving at 2 mm within 24 h of collection. At each sampling time, subsamples were dried at 100 °C overnight in order to obtain the gravimetric water content. An aliquot of 20 g of fresh soil was extracted with 100 mL K2SO4 0.5 M, whereas a 10 g aliquot was subjected to chloroform fumigation for 18 h before extraction with 50 mL K2SO4 0.5 M. Dissolved organic carbon (DOC) was determined with 0.45 μm membrane, which filtered K2SO4 extracts (extractable DOC) with a TOC analyzer (Elementar, Vario TOC, Hanau, Germany). Microbial carbon (Cmicr) was calculated from the difference in DOC between fumigated and non-fumigated samples corrected by a recovery factor of 0.45 (
Lake water was sampled approximately at the same time that soil was sampled, with a total of 79 samples analyzed in the time-span 2008–2017. The lake was sampled from 3 points on the shore with no vegetation, at ca. 10-cm depth. Shore samples were assumed not to be significantly different from mid-lake samples because the investigated lake was small and shallow (
In order to assess the interannual variability in C and N forms in soils and lake water in the time-span 2008–2017, only the samples collected at the end of the snow-free season (September) were considered (during the first 3 years the sampling was carried out only during September), while for the evaluation of the intraseasonal variability, the monthly samples (July, August, September, October) were considered. We tested significant differences among years and months by one-way ANOVA and Bonferroni post hoc test (p < 0.05). Data were previously tested for homoscedasticity (Levene’s test) and for normality (Kolmogorov–Smirnov test), and transformed when necessary.
In order to describe the variation of climatic and pedoclimatic indices over the 10 studied years, a Principal Component Analysis (PCA) was performed using all the variables described in Table
The influence of climatic and pedoclimatic variables on C and N forms in soil was evaluated by fitting generalized linear models (GLMs). Soil C and N forms were used as dependent variables and a gamma distribution was used because normality (tested with Kolmogorov-Smirnoff test) was not met (
All the statistical analyses were performed using SPSS v.19 (
The mean annual air temperature measured at the AWS in the time-span 2008–2017 was -2.3 °C, with mean daily values ranging from a minimum of -21.1 °C (16 January 2017) to a maximum of +12.4 °C (24 August 2016) (Fig.
a Air temperature, snow depth and rain recorded at the Automatic Weather Station (AWS) from 1 October 2007 to 30 September 2017 (daily mean values) b soil temperature recorded in the topsoil (A horizon – 10 cm depth) at study sites 1, 3, and 5 (daily mean values).
Climatic indices derived from the AWS data. CS (cumulative snowfall), HPD (heavy precipitation days), VHPD (very heavy precipitation days), CWD (consecutive wet days), and CDD (consecutive dry days). All indices are listed and described in Table
Index | Unit | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | Min | Max | Mean | St.dev |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CS | cm | 605 | 1099 | 879 | 752 | 890 | 947 | 820 | 820 | 756 | 608 | 605 | 1099 | 817 | 150 |
HPD | days | 12 | 5 | 5 | 3 | 8 | 5 | 2 | 5 | 6 | 6 | 2 | 12.0 | 5.7 | 2.8 |
VHPD | days | 4 | 2 | 2 | 4 | 4 | 3 | 0 | 7 | 7 | 3 | 0 | 7 | 4 | 2 |
CWD | days | 7 | 3 | 6 | 4 | 4 | 7 | 5 | 7 | 9 | 6 | 3 | 9 | 6 | 2 |
CDD | days | 18 | 9 | 5 | 9 | 6 | 5 | 23 | 12 | 9 | 16 | 5 | 23 | 11 | 6 |
Soil temperature measured at all study sites during the snow-covered season was generally close to 0 °C (Fig.
Site-specific pedoclimate indices measured at sites 1, 3, and 5 between 2008 and 2017. SCD (snow cover duration*), MOD (melt-out day of snow*), DSF (duration of soil freezing), FTC (soil freeze/thaw cycles), MTF (mean soil temperature during soil freezing), MTSC (mean soil temperature during the snow-covered season), MTSF (mean soil temperature during the snow-free season), and ISF (intensity of soil freezing). All indices are listed and described in Table
Index | Unit | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | Min | Max | Mean | St.dev |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Site 1 | |||||||||||||||
SCD | days | 249 | 272 | 271 | 255 | 241 | 284 | 286 | 250 | 283 | 262 | 241 | 286 | 265 | 16 |
MOD | DOY | 181 | 221 | 192 | 187 | 163 | 211 | 205 | 186 | 193 | 175 | 163 | 221 | 191 | 17 |
DSF | days | 205 | 246 | 269 | 208 | 235 | 217 | 223 | 15 | 102 | 124 | 15 | 269 | 184 | 79 |
FTC | number | 3 | 1 | 3 | 1 | 2 | 1 | 1 | 1 | 1 | 1 | 1 | 3 | 2 | 1 |
MTF | °C | -0.3 | -0.1 | -0.1 | -0.4 | -0.5 | -0.1 | 0.0 | 0.0 | -0.3 | -0.1 | -0.5 | 0.0 | -0.2 | 0.2 |
MTSC | °C | -0.3 | 0.0 | -0.3 | -0.2 | -0.5 | -0.1 | 0.0 | 0.1 | 0.0 | 0.1 | -0.5 | 0.1 | -0.1 | 0.2 |
MTSF | °C | 5.8 | 8.5 | 7.6 | 7.2 | 6.8 | 5.1 | 4.5 | 9.0 | 7.7 | 6.4 | 4.5 | 9.0 | 6.9 | 1.4 |
ISF | °C | -2.6 | -0.1 | -0.8 | -0.3 | -1.1 | -0.3 | -0.1 | 0.0 | -0.5 | -0.1 | -2.6 | 0.0 | -0.6 | 0.8 |
Site 3 | |||||||||||||||
SCD | days | 222 | 246 | 222 | 219 | 223 | 247 | 231 | 210 | 274 | 228 | 210 | 274 | 232 | 19 |
MOD | DOY | 182 | 174 | 179 | 139 | 162 | 184 | 162 | 148 | 188 | 146 | 139 | 188 | 166 | 18 |
DSF | days | 167 | 125 | 211 | 154 | 58 | 184 | 125 | 59 | 122 | 57 | 57 | 211 | 126 | 55 |
FTC | number | 4 | 1 | 2 | 2 | 2 | 1 | 1 | 1 | 1 | 1 | 1 | 4 | 2 | 1 |
MTF | °C | -1.0 | -0.1 | -1.1 | -0.3 | -0.9 | -1.6 | -0.5 | -0.1 | -0.9 | -1.6 | -1.6 | -0.1 | -0.8 | 0.6 |
MTSC | °C | -0.1 | 0.1 | -0.8 | -0.1 | 0.1 | -1.1 | -0.3 | 0.2 | -0.6 | 0.2 | -1.1 | 0.2 | -0.3 | 0.4 |
MTSF | °C | 8.1 | 10.0 | 9.1 | 8.5 | 8.3 | 6.6 | 7.2 | 10.0 | 6.5 | 9.1 | 6.5 | 10.0 | 8.3 | 1.3 |
ISF | °C | -2.1 | -0.1 | -3.4 | -1.0 | -2.5 | -8.3 | -1.9 | -0.2 | -7.4 | -4.1 | -8.3 | -0.1 | -3.1 | 2.8 |
Site 5 | |||||||||||||||
SCD | days | 219 | 244 | 235 | 219 | 216 | 253 | 235 | 222 | 248 | 216 | 253 | 232 | 14 | |
MOD | DOY | 183 | 178 | 164 | 139 | 160 | 182 | 161 | 157 | 165 | 139 | 183 | 165 | 14 | |
DSF | days | 190 | 156 | 176 | 147 | 163 | 160 | 181 | 0 | 68 | 0 | 190 | 138 | 63 | |
FTC | number | 4 | 1 | 1 | 1 | 1 | 1 | 1 | 0 | 1 | 0 | 4.0 | 1.2 | 1.1 | |
MTF | °C | -0.7 | -0.1 | -0.1 | -0.1 | -0.5 | -0.1 | -0.1 | -0.3 | -0.7 | -0.1 | -0.2 | 0.2 | ||
MTSC | °C | -0.1 | 0.0 | 0.0 | 0.0 | 0.0 | 0.2 | 0.0 | 0.3 | 0.2 | -0.1 | 0.3 | 0.1 | 0.1 | |
MTSF | °C | 9.2 | 11.3 | 10.8 | 9.3 | 9.7 | 7.8 | 8.4 | 10.8 | 6.5 | 6.5 | 11.3 | 9.3 | 1.6 | |
ISF | °C | -1.9 | -0.3 | -0.3 | -0.3 | -0.1 | -0.1 | -0.8 | -1.4 | -1.9 | -0.1 | -0.7 | 0.7 |
The PCA revealed the distribution of the years in three different groups (Fig.
Principal Component Analysis showing the variation of climatic (red arrows) and pedoclimatic (black arrows) indices over the 10 studied years. Blue points represent the average coordinates of different sampling date and replicates within the same year. Descriptions for each climate and pedoclimate index can be found in Table
On the interannual basis, the N-NO3– concentration was significantly higher in years 2008, 2009 and 2010, while the lowest values were recorded in years 2014, 2016 and 2017 (Fig.
a Mean extractable soil nitrate (N-NO3–) concentrations (mg kg–1) recorded in September at the three study sites in the time-span 2008–2017 (n = 87, in 2017 the site 5 was not sampled) b mean seasonal concentration N-NO3– (mg kg–1) in the 3 study sites in years 2008–2017 c mean water nitrate (N-NO3–) (mg L–1) concentration recorded in September at the Cimalegna Lake during the years 2008–2017 (n = 30) d Mean seasonal concentration (mg L-1) of N-NO3– at the Cimalegna Lake in years 2008–2017 (n = 30). Upper-case letters represent significant differences between years and months (p < 0.05).
Mean interannual and seasonal concentrations of N-NH4+, N-NO3–, DOC, TDN, Cmicr, Nmicr, DON (mg kg–1), and values of C:Nmicr and DOC:DON at the 3 study sites in years 2008–2017. Letters represent significant differences between years and months (p < 0.05), letters are not reported when differences are not significant (p > 0.05).
Parameter | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | Min | Max | Mean | St.dev |
Interannual | ||||||||||||||
N-NH4+ | 9.9 (B) | 9.7 (B) | 21.2 (A) | 6.9 (B) | 6.9 (B) | 5.4 (B) | 4.5 (B) | 5.1 (B) | 2.2 (B) | 4.4 (B) | 2.2 | 21.2 | 7.6 | 5.3 |
N-NO3- | 2.6 (B) | 7.3 (A) | 5.9 (A) | 0.4 (BC) | 1.0 (BC) | 0.9 (BC) | 0.2 (C) | 1.2 (B) | 0.2 (C) | 0.2 (C) | 0.2 | 7.3 | 2.0 | 2.5 |
DOC | 324.2 (B) | 308.2 (B) | 576.5 (A) | 251.7 (B) | 335.7 (B) | 302.6 (B) | 212.0 (B) | 365.6 (AB) | 164.9 (B) | 180.7 (B) | 164.9 | 576.5 | 302.2 | 117.9 |
TDN | 30.8 (AB) | 32.7 (AB) | 24.0 (B) | 8.2 (C) | 32.5 (AB) | 61.9 (A) | 39.7 (AB) | 44.5 (AB) | 34.8 (AB) | 36.4 (AB) | 8.2 | 61.9 | 34.5 | 13.8 |
Cmicr | 2178.4 (B) | 1150.4 (BC) | 1769.1 (BC) | 2106.5 (AB) | 3029.8 (A) | 872.3 (C) | 1001.1 (C) | 944.7 (C) | 729.2 (C) | 760.1 (C) | 729.2 | 3029.8 | 1454.2 | 777.1 |
Nmicr | 228.4 (A) | 110.6 (AB) | 102.7 (B) | 67.3 (B) | 98.5 (B) | 64.4 (B) | 135.7 (AB) | 78.2 (B) | 44.6 (B) | 35.4 (B) | 35.4 | 228.4 | 96.6 | 55.6 |
DON | 18.3 (B) | 16.1 (B) | 6.7 (B) | 4.4 (B) | 32.4 (AB) | 55.5 (A) | 35.0 (AB) | 38.8 (AB) | 32.3 (AB) | 31.7 (AB) | 4.4 | 55.5 | 27.1 | 15.7 |
C:Nmicr | 12.2 (BC) | 11.0 (C) | 19.4 (BC) | 55.5 (B) | 42.9 (A) | 25.2 (BC) | 8.1 (C) | 12.2 (BC) | 16.7 (BC) | 21.1 (BC) | 8.1 | 55.5 | 22.4 | 15.3 |
DOC:DON | 17.6 (B) | 24.5 (A) | 41.5 (A) | 33.7 (A) | 11.0 (BC) | 5.2 (C) | 6.5 (BC) | 11.7 (BC) | 5.5 (C) | 7.7 (BC) | 5.2 | 41.5 | 16.5 | 12.8 |
Jul | Aug | Sep | Oct | Min | Max | Mean | St.dev | |||||||
Seasonal | ||||||||||||||
N-NH4+ | 6.0 | 6.8 | 7.8 | 6.2 | 6.0 | 7.8 | 6.7 | 0.8 | ||||||
N-NO3- | 0.6 (A) | 0.6 (A) | 1.9 (B) | 1.0 (A) | 0.6 | 1.9 | 1.0 | 0.7 | ||||||
DOC | 299.9 | 271.6 | 303.7 | 258.2 | 258.2 | 303.7 | 283.4 | 22.0 | ||||||
TDN | 30.8 | 32.7 | 24.0 | 8.2 | 8.2 | 32.7 | 23.9 | 11.1 | ||||||
Cmicr | 1476.9 | 1066.6 | 1482.7 | 1276.7 | 1066.6 | 1482.7 | 1325.7 | 197.5 | ||||||
Nmicr | 104.3 | 107.7 | 100.8 | 90.5 | 90.5 | 107.7 | 100.8 | 7.4 | ||||||
DON | 30.6 | 39.5 | 29.2 | 31.3 | 29.2 | 39.5 | 32.7 | 4.6 | ||||||
C:Nmicr | 17.4 | 21.5 | 21.7 | 29.9 | 17.4 | 29.9 | 22.6 | 5.3 | ||||||
DOC:DON | 17.7 | 11.0 | 12.5 | 12.7 | 11.0 | 17.7 | 13.5 | 2.9 |
All soil C and N forms were positively correlated with CS, with the exception of the microbial C:N ratio that was inversely correlated (Tab.
Results of generalized linear models (GLMs) showing the effects of climatic and pedoclimatic variables on soil C and N forms. Explanatory variables were standardized (Z-scores) to allow for analysis of effect size by scrutinizing model parameters (β coefficients). p-values are also shown.
Predictor | N-NH4+ | N-NO3– | DOC | TDN | Cmicr | Nmicr | DON | C:Nmicr | DOC:DON | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
β | p | β | p | β | p | β | p | β | p | β | p | β | p | β | p | β | p | |
Climatic index | ||||||||||||||||||
CS | 0.19 | 0.000 | 0.43 | 0.000 | 0.25 | 0.000 | 0.33 | 0.000 | 0.11 | 0.004 | 0.26 | 0.000 | 0.31 | 0.000 | -0.13 | 0.014 | -0.25 | 0.001 |
HPD | -0.02 | 0.772 | 0.37 | 0.000 | -0.02 | 0.575 | 0.13 | 0.029 | 0.08 | 0.134 | 0.11 | 0.155 | 0.14 | 0.083 | 0.04 | 0.566 | -0.35 | 0.002 |
VHPD | 0.07 | 0.225 | 0.10 | 0.340 | 0.06 | 0.141 | -0.14 | 0.019 | 0.10 | 0.077 | -0.06 | 0.477 | -0.18 | 0.018 | 0.08 | 0.306 | -0.02 | 0.835 |
CWD | -0.14 | 0.018 | -0.20 | 0.053 | -0.05 | 0.226 | -0.02 | 0.783 | -0.13 | 0.037 | -0.21 | 0.016 | -0.03 | 0.722 | 0.13 | 0.159 | -0.05 | 0.711 |
CDD | 0.03 | 0.473 | -0.04 | 0.516 | 0.01 | 0.621 | 0.01 | 0.106 | -0.07 | 0.100 | 0.08 | 0.111 | 0.04 | 0.486 | -0.22 | 0.000 | -0.32 | 0.000 |
Pedoclimatic index | ||||||||||||||||||
SCD | -0.38 | 0.000 | -0.32 | 0.000 | -0.39 | 0.000 | -0.09 | 0.042 | -0.45 | 0.000 | -0.33 | 0.000 | -0.06 | 0.307 | -0.21 | 0.001 | -0.62 | 0.000 |
DSF | 0.14 | 0.010 | 0.00 | 0.646 | 0.02 | 0.656 | -0.28 | 0.000 | 0.15 | 0.009 | -0.13 | 0.117 | -0.31 | 0.000 | 0.52 | 0.000 | 0.45 | 0.000 |
FTC | 0.13 | 0.015 | 0.20 | 0.014 | 0.12 | 0.001 | -0.07 | 0.164 | 0.12 | 0.007 | 0.23 | 0.002 | -0.17 | 0.011 | -0.09 | 0.238 | 0.05 | 0.614 |
MTF | -0.02 | 0.685 | 0.06 | 0.481 | -0.03 | 0.481 | -0.21 | 0.001 | -0.04 | 0.523 | 0.18 | 0.032 | -0.22 | 0.004 | -0.31 | 0.000 | 0.37 | 0.002 |
MTSC | 0.00 | 0.991 | -0.16 | 0.082 | -0.03 | 0.509 | 0.04 | 0.484 | 0.01 | 0.210 | -0.13 | 0.116 | 0.00 | 0.693 | 0.33 | 0.000 | 0.09 | 0.445 |
MTSF | 0.05 | 0.277 | -0.07 | 0.351 | -0.01 | 0.753 | -0.05 | 0.285 | 0.02 | 0.620 | -0.03 | 0.606 | -0.01 | 0.902 | 0.00 | 0.968 | -0.15 | 0.177 |
On an interannual basis, the maximum N-NO3– concentration in water was recorded in 2008, while the lowest was measured in 2011 and 2012 (Fig.
Mean interannual and seasonal concentrations of N-NH4+, N-NO3–, DOC, TDN, DON (mg L–1), and values of DOC:DON at the Cimalegna Lake. Letters represent significant differences between years and months (p < 0.05), letters are not reported when differences are not significant (p > 0.05).
Parameter | 2008 | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | 2015 | 2016 | 2017 | Min | Max | Mean | St.dev |
Interannual | ||||||||||||||
N-NH4+ | 0.2 (B) | 0.0 (C) | 0.3 (A) | 0.0 (C) | 0.0 (C) | 0.0 (C) | 0.0 (C) | 0.0 (C) | 0.0 (C) | 0.0 (C) | 0.0 | 0.3 | 0.1 | 0.1 |
N-NO3- | 0.6 (A) | 0.2 (C) | 0.1 (DE) | 0.1 (E) | 0.1 (E) | 0.1 (D) | 0.0 (F) | 0.1 (DE) | 0.4 (B) | 0.1 (DE) | 0.0 | 0.6 | 0.2 | 0.2 |
DOC | 2.5 (A) | 1.6 (B) | 0.7 (C) | 1.2 (BC) | 1.5 (B) | 1.6 (B) | 1.1 (BC) | 0.7 | 2.5 | 1.5 | 0.6 | |||
TDN | 0.7 (A) | 0.4 (B) | 0.2 (CD) | 0.2 (CD) | 0.3 (C) | 0.5 (B) | 0.2 (D) | 0.2 | 0.7 | 0.4 | 0.2 | |||
DON | 0.6 (A) | 0.3 (B) | 0.1 (CD) | 0.2 (C) | 0.3 (B) | 0.1 (D) | 0.1 (D) | 0.1 | 0.6 | 0.2 | 0.2 | |||
DOC:DON | 4.3 (B) | 5.6 (B) | 6.8 (B) | 6.0 (B) | 6.1 (B) | 19.7 (A) | 18.3 (A) | 4.3 | 19.7 | 9.5 | 6.5 | |||
Jul | Aug | Sep | Oct | Min | Max | Mean | St.dev | |||||||
Seasonal | ||||||||||||||
N-NH4+ | 0.0 | 0.0 | 0.1 | 0.0 | 0.0 | 0.1 | 0.0 | 0.0 | ||||||
N-NO3- | 0.1 (AB) | 0.1 (A) | 0.2 (B) | 0.1 (AB) | 0.1 | 0.2 | 0.1 | 0.1 | ||||||
DOC | 1.7 | 1.8 | 1.4 | 1.5 | 1.4 | 1.8 | 1.6 | 0.2 | ||||||
TDN | 0.5 | 0.4 | 0.4 | 0.4 | 0.4 | 0.5 | 0.4 | 0.1 | ||||||
DON | 0.3 | 0.3 | 0.2 | 0.2 | 0.2 | 0.3 | 0.3 | 0.1 | ||||||
DOC:DON | 6.3 | 6.6 | 10.5 | 6.7 | 6.3 | 10.5 | 7.5 | 2.0 |
Among the selected climatic indices, CS was inversely correlated with N-NH4+ (Tab.
Results of generalized linear models (GLMs) showing the effects of climatic, pedoclimatic and soil variables on water C and N forms. Explanatory variables were standardized (Z-scores) to allow for analysis of effect size by scrutinizing model parameters (β coefficients). p-values are also shown.
Predictor | N-NH4+ | N-NO3– | DOC | TDN | DON | DOC:DON | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
β | p | β | p | β | p | β | p | β | p | β | p | |
Climatic index | ||||||||||||
CS | -0.32 | 0.023 | 0.15 | 0.342 | -0.02 | 0.760 | -0.02 | 0.756 | -0.09 | 0.460 | 0.07 | 0.536 |
HPD | 0.07 | 0.684 | 0.24 | 0.201 | 0.03 | 0.744 | -0.07 | 0.352 | 0.12 | 0.419 | -0.11 | 0.359 |
VHPD | -0.72 | 0.000 | 0.01 | 0.950 | -0.01 | 0.831 | 0.00 | 0.966 | -0.19 | 0.008 | 0.11 | 0.115 |
CWD | 0.46 | 0.003 | 0.28 | 0.060 | 0.13 | 0.038 | 0.12 | 0.061 | 0.09 | 0.434 | -0.01 | 0.883 |
CDD | -0.19 | 0.161 | -0.38 | 0.000 | 0.08 | 0.168 | -0.01 | 0.924 | 0.07 | 0.466 | -0.11 | 0.173 |
Pedoclimatic index | ||||||||||||
SCD | 1.07 | 0.000 | 0.59 | 0.029 | -0.06 | 0.689 | -0.02 | 0.922 | -0.25 | 0.379 | 0.35 | 0.160 |
DSF | -0.24 | 0.120 | -0.37 | 0.106 | -0.02 | 0.878 | 0.11 | 0.307 | 0.26 | 0.154 | -0.38 | 0.029 |
FTC | 0.27 | 0.372 | 0.59 | 0.046 | -0.25 | 0.313 | -0.54 | 0.028 | -0.69 | 0.134 | 0.31 | 0.447 |
MTF | 0.21 | 0.156 | -0.38 | 0.023 | -0.21 | 0.022 | -0.51 | 0.000 | -0.44 | 0.004 | 0.13 | 0.359 |
MTSF | 0.12 | 0.424 | -0.24 | 0.071 | 0.03 | 0.675 | -0.06 | 0.392 | -0.15 | 0.193 | 0.30 | 0.003 |
Soil parameter | ||||||||||||
S_N-NH4+ | 0.16 | 0.404 | 0.47 | 0.023 | -0.12 | 0.228 | -0.06 | 0.528 | -0.19 | 0.211 | 0.23 | 0.124 |
S_N-NO3– | -0.25 | 0.047 | 0.43 | 0.000 | -0.13 | 0.210 | 0.04 | 0.712 | 0.21 | 0.223 | -0.48 | 0.001 |
S_DOC | 0.49 | 0.006 | 0.25 | 0.078 | 0.17 | 0.007 | 0.22 | 0.001 | 0.27 | 0.010 | -0.12 | 0.229 |
S_Cmicr | 0.43 | 0.009 | 0.28 | 0.253 | 0.13 | 0.288 | -0.02 | 0.886 | -0.35 | 0.111 | 0.64 | 0.001 |
S_Nmicr | -0.38 | 0.082 | -0.49 | 0.046 | 0.04 | 0.727 | 0.00 | 0.984 | 0.43 | 0.041 | -0.60 | 0.001 |
S_DON | -0.24 | 0.075 | -0.15 | 0.356 | -0.17 | 0.026 | -0.09 | 0.232 | -0.08 | 0.560 | -0.17 | 0.129 |
S_C:Nmicr | 0.00 | 0.982 | -0.30 | 0.021 | 0.00 | 0.966 | -0.09 | 0.158 | 0.18 | 0.116 | -0.23 | 0.049 |
S_DOC:DON | -0.26 | 0.053 | -0.42 | 0.009 | 0.08 | 0.299 | -0.05 | 0.478 | 0.03 | 0.831 | -0.08 | 0.477 |
Along the study period, the climate and pedoclimate conditions showed a great interannual variability, with some extreme meteorological events. For example, in 2008, the little and delayed snowpack accumulation and consequently the low insulation exerted by the snowpack caused a high number of soil FTC in all sites. In high-elevation ecosystems, the most frequent periods for FTC are spring and fall, when the soil cannot be covered by a consistent snowpack. Sometimes they can also occur throughout the winter, due to little snowpack accumulation and/or the wind action that causes a snow removal exposing the soil to cold air temperature (
In our research sites, the resulting FTC number recorded during the snow-covered season significantly and positively correlated with the concentrations of most of the soil C and N forms in the subsequent growing season (cf.
The duration and intensity of soil freezing during the snow-covered season was positively related to an increase in the soil and microbial C:N ratio, suggesting the prevalence of fungi, characterized by a higher C:N ratio in comparison to bacteria.
In the time-span considered in this study, the interannual variability of the SCD and melt-out day was marked and it was possible to discriminate years with short and long snow-cover duration, which corresponded to early and late melt-out days, respectively. According to the conceptual model of
Among the climatic indices, the cumulative snowfall was positively related to all the soil C and N forms. We assume that a higher cumulative snowfall in the study area corresponded to a higher soil water content and nutrient inputs into the soil during the spring melting of snow, enhancing the microbial activity and the C and N transformations. In the same area,
In contrast to what was found by
A number of studies have demonstrated that some physical features of the catchment strongly influence the chemical composition of surface water, and may control ecosystem responses to global perturbations, such as changes in climate (
In our study the resultant N-NO3– content in lake water positively correlated with the soil inorganic N forms. This concurs with the findings of
As reported for the soil matrix, the number of FTC had an important control on the N-NO3– content in lake water and, as reported by
The SCD had a first order control on N-NH4+ and N-NO3– concentration in water, but an opposite pattern was observed in comparison to soil C and N forms. A longer SCD caused an increase of both N-NH4+ and N-NO3– concentrations in the lake water, which could be related to the reduction of the microbial nutrient immobilization processes in soil. In our study DOC concentration in water was positively related to the DOC concentration in soil, revealing how allochthonous DOC could represent a large fraction of the total DOC in lakes (
In the LTER site Istituto Mosso the C and N forms analyzed in soil and water for a decade showed a significant interannual variability, while a seasonal change was observed only for N-NO3– both in soil and water lake, with the greatest values recorded in early fall, probably due to the slowdown of biological-mediated processes of N immobilization.
Both the climatic and pedoclimatic conditions recorded during the snow-free and snow-covered season significantly influenced the C and N forms in soil and water. A little and delayed snowpack accumulation caused a high number of soil freeze/thaw cycles, which resulted in a high nitrate content both in soil and water. The longer the snow cover duration, the lower are all the soil C and N forms measured during the subsequent snow-free season, with the exception of DON. An opposite trend was observed for the lake water, where a longer snow-cover duration caused a higher content of inorganic N forms, probably due to a reduction in soil N immobilization potential.
This study was supported by NextData Data-LTER-Mountain Project. This research has been also partially developed in the framework of the European Regional Development Fund in Interreg Alpine Space project Links4Soils (ASP399): Caring for Soil-Where Our Roots Grow (http://www.alpine-space.eu/projects/links4soils). Thanks to the Comando Truppe Alpine - Ufficio Meteomont for the data from the AWS Col d’Olen and to Monterosa 2000 and Monterosa S.p.A. (MonterosaSki) for the logistic support.