Responses of dissolved organic carbon to freeze–thaw cycles associated with the changes in microbial activity and soil structure
Arctic warming accelerates snowmelt, exposing soil surfaces with shallow or no snow cover to freeze–thaw cycles (FTCs) more frequently in early spring and late autumn. FTCs influence Arctic soil C dynamics by increasing or decreasing the amount of dissolved organic carbon (DOC); however, mechanism-based explanations of DOC changes that consider other soil biogeochemical properties are limited. To understand the effects of FTCs on Arctic soil responses, we designed microcosms with surface organic soils from Alaska and investigated several soil biogeochemical changes for seven successive temperature fluctuations of freezing at −9.0 ± 0.3 ∘C and thawing at 6.2 ± 0.3 ∘C for 12 h each. FTCs significantly changed the following soil variables: soil CO2 production (CO2), DOC and total dissolved nitrogen (TDN) contents, two DOC quality indices (SUVA254 and A365 A254), microaggregate (53–250 µm) distribution, and small-sized mesopore (0.2–10 µm) proportion. Multivariate statistical analyses indicated that the FTCs improved soil structure at the scale of microaggregates and small-sized mesopores, facilitating DOC decomposition by soil microbes and changes in DOC quantity and quality by FTCs. This study showed that FTCs increased soil CO2 production, indicating that FTCs affected DOC characteristics without negatively impacting microbial activity. Soil microaggregation enhanced by FTCs and the subsequent increase in microbial activity and small-sized pore proportion could promote DOC decomposition, decreasing the DOC quantity. This study provides a mechanism-based interpretation of how FTCs alter DOC characteristics of the organic soil in the active layer by incorporating structural changes and microbial responses, improving our understanding of Arctic soil C dynamics.
Arctic tundra soils store approximately 1300 ± 200 Pg of soil organic carbon (SOC) in permafrost (Čapek et al., 2015), which accounts for approximately 30 % of the global SOC pool (Xu et al., 2009). Recently, Arctic warming, which occurs 4 times faster than global warming (Rantanen et al., 2022), has enhanced permafrost thaw, causing the previously stored SOC to be released into greenhouse gases (CO2 and CH4) and/or leaching dissolved organic carbon (DOC) (Estop-Aragonés et al., 2020). In particular, DOC released from the active layer could be further decomposed by soil microorganisms into CO2 and CH4, leading to a positive feedback on permafrost thawing (Foster et al., 2016; Yi et al., 2015). Permafrost thaw also influences the Arctic watershed through the export of terrestrial-derived DOC into the surrounding lakes and seas (Al-Houri et al., 2009). The exported DOC from the active layer can horizontally migrate along the unfrozen vicinity between frozen layers; in addition, it can infiltrate the deeper active layer and upper permafrost during the thawing phase (Ban et al., 2016; Han et al., 2018). Thus, the measurement of quantitative and qualitative DOC changes are necessary for understanding the response of permafrost C dynamics to Arctic warming (Xu et al., 2009; Foster et al., 2016; Perez-Mon et al., 2020).
Moreover, increased temperatures in Arctic regions accelerate snowmelt (Henry, 2008; Førland et al., 2011; Kreyling et al., 2008) and cause rainfall instead of snowfall (Henry, 2013; IPCC, 2014), leading to the absence of snow cover on the soil surface (Callaghan et al., 1998; Heal et al., 1998). In Arctic regions, snow plays a key role in protecting tundra soils against dramatic temperature changes caused by harsh climates (Royer et al., 2021). Exposed soil surfaces lacking snow cover are likely to undergo more frequent freeze–thaw cycles (FTCs) in the early spring and late autumn because they are directly influenced by the diurnal fluctuations of atmospheric temperature (Kreyling et al., 2008; Henry, 2013; Freppaz et al., 2007). Climate models project that the air temperature of the Arctic may continue to rise owing to climate change, thereby enhancing the occurrence of FTCs in permafrost soils within the near future (Henry, 2008; Groffman et al., 2011; Grogan et al., 2004).
Numerous studies have reported influences of FTCs on the labile soil C content, which is strongly related to microbial activity in Arctic tundra soils (Sawicka et al., 2010; Schimel and Clein, 1996; Larsen et al., 2002; Männistö et al., 2009; Lipson and Monson, 1998; Perez-Mon et al., 2020; Grogan et al., 2004; Foster et al., 2016; Schimel and Mikan, 2005; Sjursen et al., 2005; Yi et al., 2015). FTCs have been reported to increase the amount of DOC in a few tundra and non-tundra soils, attributed to a decrease in microbial utilization of DOC due to cell lysis generally occurring below −7 to −11 ∘C of freezing temperature (D. Gao et al., 2018, 2021; Song et al., 2017; Schimel and Clein, 1996; Larsen et al., 2002). In contrast, several studies have shown negligible changes or decreases in the amount of DOC by FTCs, without negative responses of microbial biomass, community, and enzymatic activities (Männistö et al., 2009; Lipson and Monson, 1998; Perez-Mon et al., 2020). This was interpreted as soil microorganisms in the Arctic tundra having already adapted to extreme temperature fluctuations for a long period of time; thus, FTCs could not inhibit microbial DOC utilization for growth and activity compared to non-tundra regions, such as forest, grassland, and cropland (D. Gao et al., 2018, 2021; Song et al., 2017). These controversial results suggest that further research and evidence of DOC changes by FTCs are required for an improved mechanism-based understanding of tundra soil C dynamics in the early spring and late autumn.
FTCs can indirectly affect the Arctic tundra DOC dynamics through soil structural changes such as the fragmentation, rearrangement, and aggregation of soil particles (Matzner and Borken, 2008; Zhang et al., 2016). Owing to the phase transitions in soil water during FTCs, soil matrix cracks and the physical degradation of soil aggregates have been reported in previous studies (Oztas and Fayetorbay, 2003; Wang et al., 2012; Hall and André, 2003). In contrast, several researchers have found that FTCs enhance soil aggregate stability (Lehrsch, 1998) and small-sized aggregate formation (50–250 and 500–1000 µm) (Li and Fan, 2014). Changes in soil aggregate distribution by FTCs likely affect the soil pore volume and spatial distribution (Lu et al., 2021; Al-Houri et al., 2009; Oztas and Fayetorbay, 2003; Viklander, 1998), leading to alterations in soil water retention and DOC release (Matzner and Borken, 2008; Song et al., 2017; D. Gao et al., 2018; Feng et al., 2007). Additionally, these soil structural changes may affect microbe-mediated soil C mineralization and utilization by improving the soil water and nutrient distribution (Athmann et al., 2013; Liang et al., 2019; Sander and Gerke, 2007). However, the linkage of how structural changes caused by FTCs, such as the formation of aggregates and pores with specific sizes, affect DOC changes has not been well understood.
This study aimed to identify the effects of FTCs on Arctic tundra DOC dynamics using surface organic soils from the Alaskan tundra undergoing temperature fluctuation during the early spring. We designed two parallel microcosms simulating the FTCs of the study site for two different purposes. One set of microcosms was established for destructive sampling to investigate the temporal changes in soil CO2 production, soil enzyme activities, DOC characteristics, and soil aggregate size distribution. The other set of soil core incubation was prepared with repacked soil to measure pore size distribution (PSD) using soil water retention curves. We tested the following hypotheses: (1) FTCs alter DOC quantity and quality without decreasing the activities of soil microbes previously adapted to temperature fluctuations in the Arctic, and (2) soil aggregate distribution influenced by FTCs changes DOC characteristics by enhancing microbial activities and altering specific-sized soil pore proportion.
2.1 Site description and soil preparation
Soil samples for microcosm incubation were collected from the moist acidic tundra in Council (64.51∘ N, 163.39∘ W) on the Seward Peninsula in Northwest Alaska. The average temperature and precipitation over the past 30 years (1981–2020) are −3.1 ∘C and 258 mm (Alaska Climate Research Center). In the early spring (April to May), the minimum and maximum temperatures are −8.5 and 7.1 ∘C, respectively (Alaska Climate Research Center). This site is in the discontinuous permafrost zone and tussock tundra dominated by cotton grasses (Eriophrum vaginatum), blueberries (Vaccinium uliginosum), lichen, and moss (Sphagnum spp.) (Kim et al., 2016).
Soil sampling was performed at three random points under similar vegetation compositions. Each point was within a distance of approximately 100 m from each other. At the time of sampling (early July 2010), the active layer depth was approximately 25–30 cm, measured by a steel rod (1 m). Soil samples were acquired by hammering a stainless steel pipe (7.6 cm diameter × 50 cm long) into the partially or well-degraded organic layer (Oe), mixed with soil minerals, after removing the litter layer (Oi) on the surface. The soil samples were stored at −20 ∘C before initiating microcosm incubation. The frozen soil was thawed at < 4 ∘C, and the surface organic soils were passed through a 2 mm sieve and homogenized by hand. Figure 1 shows the soil sampling and preparation procedure. Soil textural analysis was conducted by a wet sieving and pipette method (Kim et al., 2022). Soil bulk density (BD) was determined by calculating the soil dry weight contained in the soil core volume. Volumetric water content (VWC) in the soil was measured using a portable sensor with an accuracy of ±0.01 cm3 cm−3 (ProCheck Decagon Devices, Washington, USA). Total carbon (C) and nitrogen (N) contents were determined through combustion (950 ∘C) with an elemental analyzer (vario MAX cube; Elementar vario MAX cube; Elementar, Langenselbold, Germany). The basic soil properties are summarized in Table S1 in the Supplement.
2.2 Soil incubation with freeze–thaw cycles
Soil incubation was conducted with two parallel sets of microcosms: one for destructive sampling and the other for monitoring soil PSD changes (Fig. 1). The destructive sampling set was established using a 380 mL polypropylene bottle to investigate the soil biogeochemical properties influenced by FTCs. The other microcosm set was created by reconstructing the small-sized soil core (5 cm diameter × 5 cm long) to compare PSD alterations under incubation conditions with/without the impact of FTCs. We established FTC and CON as experimental groups. The FTC is a treatment with seven successive temperature fluctuations of freezing at −9.0 ± 0.3 ∘C and thawing at 6.2 ± 0.3 ∘C for 12 h each (Fig. S1 in the Supplement). The meta-analysis and several other studies showed that soil carbon dynamics (Gao et al., 2021) and total porosity (Liu et al., 2021; Ma et al., 2021) responded to 3–10 FTCs; thus, seven successive FTCs were adopted in this study. CON is a control group that maintained an average temperature of −2 ∘C without any fluctuations (Fig. S1). The freezing and thawing temperatures in the FTC treatment corresponded to the early spring conditions observed at the study site. We ensured that our FTC treatment was adequate for complete freezing and thawing of the soil based on previous studies conducted under similar conditions (Freppaz et al., 2007; Larsen et al., 2002; Song et al., 2017; Han et al., 2018). Three replicates were used for the CON and FTC treatments of each microcosm set. In all incubation sets, the initial soil BD was adjusted to 0.72 g cm−3, similar to field-soil conditions (Fig. 1; Table S1). For the destructive sampling set, we used 260 g of homogenized fresh soil (154.7 g dry weight) to a soil volume of 215 cm3. The microcosm set using the small-sized cores was established with 120 g of homogenized fresh soil (71.5 g dry weight) in a 99 cm3 volume. The VWC for all incubation soils was also standardized by spraying water using a pipette to 0.50 cm3 cm−3 (70 % water-filled pore space), a similar level to field soils (Fig. 1; Table S1).
2.3 Soil analyses
All soil analyses, except for the PSD measurement, were conducted using the first incubation set for destructive sampling. Soil CO2 production (CO2), widely accepted as a proxy for overall microbial activity (Kim and Yoo, 2021; Maikhuri and Rao, 2012; Davidson et al., 1998; Kuzyakov and Domanski, 2000; Lipson and Schmidt, 2004; Raich and Schlesinger, 1992), was estimated by daily measurement of the CO2 flux from soil incubation during the entire incubation period. We collected gas samples from the headspace through a septum using 10 mL syringes (BD Luer-Lok tip, BD Company, Franklin Lakes, NJ, USA) before and after sealing the incubation bottle for 60 min. The gas samples were analyzed using gas chromatography (Agilent 7890A, Santa Clara, CA, USA) with a hydrogen flame ionization detector to determine the CO2 concentration. The CO2 flux was calculated based on changes in headspace concentration over 60 min using the following Eq. (1) (Troy et al., 2013):
where dGasdt is the change in the CO2 concentration before and after sealing the incubation bottle for 60 min, V and A are the volume and area of the incubation bottle, P is the atmospheric pressure (1 atm), MW is the molecular weight of CO2 (44.01 g mol−1), R is a gas constant (0.082 atm L mol−1 K−1), and T is the absolute temperature during gas collection (293 K). In addition, we calculated the mean CO2 production ( by averaging the daily measured CO2 flux during the entire incubation period.
Several soil biogeochemical properties were measured at the end of seven successive FTCs. Soil extracellular enzyme activity was determined: two oxidases (peroxidase and phenol-oxidase) and four hydrolases (β-glucosidase, cellobiase, N-acetyl-β-glucosaminidase, and aminopeptidase) involved in soil C and N cycling (Liao et al., 2022) were identified. These enzyme activities were measured by fluorometric assays using L-3,4-dihydroxyphenylalanine (L-DOPA) solution for oxidases and methylumbelliferyl (MUF)-linked substrates for hydrolases (Kwon et al., 2013).
To quantify the available C and N in soils, we measured the DOC and total dissolved nitrogen (TDN) content via water extraction. After adding 40 mL of distilled water, 20 g of fresh soil was shaken for 1 h, centrifuged, and filtered with a 0.45 µm filter to obtain supernatant. The supernatants were measured using a Multi N/C 3100 analyzer (Analytik Jena, Jena, Thüringen, Germany). The filtered samples were also used to estimate DOC qualitative indices. The specific ultraviolet absorbance at 254 nm (SUVA254), which allowed for the estimation of DOC aromaticity, was calculated using UV absorbance at 254 nm (A254) divided by DOC concentration (mg C L−1) and the path length (m) of the UV cuvette of the spectrometer (Eppendorf, Hamburg, Germany) (Lim et al., 2021). The ratio of A254 to A365 (A254 A365) was used as a proxy that is negatively related to the molecular weight of the DOC compounds (Berggren et al., 2018). For NH-N and NO-N content analysis, 5 g of fresh soil was shaken with a 2 M KCl solution for 1 h, centrifuged, and filtered through Whatman #42 paper. The filtrates were analyzed using an auto-analyzer (QuAAtro, SEAL Analytical GmbH, Norderstedt, Schleswig-Holstein, Germany).
Soil aggregate fractionation was performed using density separation and a subsequent wet sieving method at the end of incubation (Kim et al., 2021; Yoo et al., 2017). Then, 20 g of air-dried and 2 mm sieved soil was mixed with 35 mL of distilled water for 30 min. The soil–water mixture was left overnight and then centrifuged at 3200 rpm for 10 min. The supernatant, which was a free-light fraction (< 1.0 g cm−3), was collected using precombusted glass microfiber filters (GF/A). The heavy fraction was wet-sieved using 1000, 250, and 53 µm sieves to separate water-stable aggregates into four size classes: mega-aggregates (1000–2000 µm), macroaggregates (250–1000 µm), microaggregates (53–250 µm), and mineral-associated fractions (< 53 µm). The wet sieving procedure was performed by manually shaking each sieve 100 times for more than 2 min. All aggregate fractions remaining on the GF/A filters and sieves were transferred to an aluminum dish, dried in an oven at 60 ∘C for a week, and then weighed.
To estimate the PSD on the post-incubation core soils (5 cm diameter × 5 cm long), soil water release curves were generated by the Hydrus-1D model equipped with van Genuchten soil-hydraulic equations, which can be applied to organic and mineral soil (Šimůnek et al., 2013; Dettmann et al., 2014). The modeling procedure requires the van Genuchten parameters, calculated using volumetric water content at field capacity and wilting point (Kameyama et al., 2012; Likos et al., 2014). The volumetric water content at field capacity was measured by the soil water content at a matric potential of −33 kPa using a sand box (Eijkelkamp Agrisearch Equipment, Santa Barbara, CA, USA) (Yoo et al., 2020), after saturating soil core samples. As the matric potential used in a sandbox did not sufficiently cover the entire soil water release curves, the volumetric water content at the wilting point (−1500 kPa) was calculated using the pedotransfer function from soil carbon content and bulk density (da Silva and Kay, 1997). Lastly, the PSD was estimated from the matric potential corresponding to each pore size using the Young–Laplace equation (Kim et al., 2021).
Note: The asterisks and * indicate significant differences between treatments at the p< 0.05 and p< 0.10 levels, respectively. The numbers in parentheses are standard errors (n=3).
2.4 Statistical analyses
Analysis of variance (ANOVA) was performed using the generalized linear model (GLM) procedure (SAS 9.4, SAS Institute Inc., Cary, NC, USA) to compare the measurement data between the CON and FTC treatments. Least square means were used to assess significant differences among treatments at p < 0.05. Following ANOVA, we performed a principal component analysis (PCA) using the “FactoMineR” package in RStudio 4.2.1 (Rstudio Inc., Boston, MA, USA) to verify whether soil variables with significant responses could discriminate the FTC soil from the CON soil. The Pearson's correlation analysis was conducted using the CORR procedure (SAS 9.4) to examine the relationship between soil variables. Finally, multiple linear regression (MLR) analyses were employed to illustrate the mechanisms by which the FTC-influenced soil biogeochemical variables contributed to changes in DOC characteristics. Based on the results of ANOVA, PCA, and the Pearson's correlation analysis, we extracted all plausible interactions among soil physical and biogeochemical variables significantly affected by FTCs and refined these interactions by finding the best-fitting regression sets from MLR. The MLR analyses were generated using SigmaPlot 13.0.
3.1 Soil biogeochemical and structural changes by freeze–thaw cycles
The quantity and quality of DOC in the soil solution were altered by FTCs, as presented in Table 1. The FTC soil exhibited lower DOC and TDN contents by 29 % and 35 %, respectively, compared to the CON soil (p< 0.001). As proxies for DOC quality, SUVA254 was higher (p=0.002), but A254 A365 was lower (p=0.016) in the FTC soil than in the CON soil. The increase in SUVA254 indicates an increase in the aromaticity of DOC (Lim et al., 2021), while the decrease in A254 A365 reflects an increase in the molecular weight of DOC (Berggren et al., 2018). In contrast, no significant changes in inorganic N (-N and -N) content were determined with FTCs (p>0.100).
The mean CO2 production ( in the FTC soil was 42.54 mg m−2 h−1, 12 times higher than that in the CON soil (3.65 mg m−2 h−1, p=0.004), as shown in Table 2. This is because the in FTC soil was significantly higher than in CON soil from the early stages of FTCs (p < 0.05) and remained consistently higher until the end of the incubation (Fig. S2). Conversely, no significant differences (p>0.10) were observed between treatments in all types of microbial extracellular enzyme activities (Table 2).
Note: The asterisks and * indicate significant differences between treatments at the p < 0.05 and p < 0.10 levels, respectively. The numbers in parentheses are standard errors (n=3).
The FTC resulted in minor differences in the mass proportion of microaggregates (53–250 µm) and mineral-associated fractions (< 53 µm), which account for an average of 35 % and 34 % of the total in each soil (Table 3). The mass proportion of microaggregates marginally increased by 17 % in the FTC soil compared to that in the CON soil (p=0.066). Although the mineral-associated fractions were insignificantly reduced by FTCs (p=0.257), the reduction level (18 %) corresponded to the increased distribution of microaggregates by FTCs. Moreover, FTCs caused a significant difference in the PSD, particularly in the small-sized mesopores (0.2–10 µm), which accounted for 44 %–45 % of the total soil pores (Table 4), as estimated using water retention curves (Fig. S3). Despite the small magnitude of difference, the proportion of small-sized mesopores in the FTC soil exhibited a statistically significant increase compared to that in the CON soil (p=0.024).
3.2 Influencing variables deriving dissolved organic carbon changes by freeze–thaw cycles
The PCA was used to further identify relationships among seven soil variables that showed significant responses to FTCs: DOC and TDN contents, SUVA254, A254 A365, , microaggregates, and small-sized mesopores (Tables 1, 2, 3, and 4). The first two principal components (PCs) accounted for 93.9 % of the total variance, with PC1 clearly clustering the FTC and CON treatments (Fig. 2). PC1 exhibited positive correlations with SUVA254, , microaggregates, and small-sized mesopores, while it showed negative correlations with DOC, TDN, and A254 A365. In Fig. 2, the microaggregate was nearly perpendicular to the DOC and TDN contents, SUVA254, and A254 A365, indicating a weak or no correlation between them. This is consistent with the results of the Pearson's correlation analysis (Fig. 3). The DOC and TDN contents showed strong correlations with and small-sized mesopores (p < 0.05) but had a weaker correlation with microaggregates at a significance level of p < 0.10. Furthermore, there were no significant correlations between microaggregates and proxies for DOC quality, including SUVA254 and A254 A365 (p>0.10). Lastly, a conceptual diagram (Fig. 4) was created using MLR analyses (Table 5) to depict the relationships between soil structural properties, microbial activity, and DOC quantity and quality as influenced by FTCs. In Table 5, and small-sized mesopores were the best-fitting variables for explaining the contents of DOC (adjusted R2=0.911, p=0.012) and TDN (adjusted R2=0.869, p=0.022). The variance inflation factors (VIFs) resulting from the MLR analyses were < 10, indicating no collinearity between and small-sized mesopores as independent variables. The variables for best representing SUVA254 and A254 A365 were (adjusted R2=0.703, p=0.023) and small-sized mesopores (adjusted R2=0.878, p=0.004), respectively. The addition of microaggregates reduced the adjusted R2 of the best-fitting regression for explaining DOC, TDN, SUVA, and A254 A365. Microaggregates correlated with (adjusted R2=0.557, p=0.054) and small-sized mesopores (adjusted R2=0.618, p=0.039). As a result, we speculated that microaggregates indirectly, rather than directly, affected the quantitative and qualitative DOC variables through its correlation with and small-sized mesopores, as illustrated in Fig. 4.
Note: R2 denotes coefficient of determination; denotes adjusted R2; VIF denotes variance inflation factor; is the mean CO2 production. The asterisks and * indicate significant differences between treatments at the p < 0.05 and p < 0.10 levels, respectively.
4.1 Effects of freeze–thaw cycles on dissolved organic carbon associated with microbial activities
The seven successive FTCs reduced soil DOC and TDN contents compared to the non-treated condition, aligning with the expectation that the quantitative characteristics of DOC were significantly affected by FTCs (Table 1). These results indicate that FTCs can accelerate the microbial decomposition of labile organic matter (Grogan et al., 2004; Han et al., 2018; Foster et al., 2016; Gao et al., 2021). A proxy for overall microbial activity, CO2 production, remained high throughout the incubation under the influence of FTCs (Fig. S2; Table 2). The main reason for our findings likely results from the soil microbial characteristics in the Arctic tundra. In other words, soil microorganisms have already adapted to the frequent temperature fluctuations in early spring and late autumn in the Arctic tundra (Perez-Mon et al., 2020; Koponen and Bååth, 2016; Walker et al., 2006; Song et al., 2017). Soil microbes in the Arctic tundra could survive at temperatures below −7 to −11 ∘C (Lipson et al., 2000; Männistö et al., 2009; Lipson and Monson, 1998), the general threshold for microbial cell lysis in non-tundra environments (D. Gao et al., 2018, 2021; Song et al., 2017). The microbes that can survive under these freezing conditions actively play a role in decomposing available DOC in the surface organic layer during thaw phases in FTCs. In addition, the top organic layer was composed of a higher quality plant-derived organic matter compared to the underlying mineral layer in the moist acidic tundra in Council, Alaska, which is the same as our study site (White et al., 2004). Thus, the biologically labile DOC could be available in the surface organic layer (L. Gao et al., 2018). Hence, decreases in DOC associated with activated microbial activities following FTCs suggest that responses of the DOC in the organic layer to FTCs would be crucial in affecting the tundra C cycle under Arctic warming. More frequent FTCs and a longer thawing length in tundra soils with warming could enhance soil C availability in the active layer of the Arctic terrestrial ecosystems, leading to a high risk of CO2 being released into the atmosphere (Estop-Aragonés et al., 2020).
Meanwhile, FTCs did not significantly change the activities of extracellular enzymes (Table 2), which are released by soil microbes to obtain C and N from recalcitrant soil organic matter such as cellulose, chitin, polypeptides, and lignin (Sinsabaugh, 2010; Liao et al., 2022). However, the enzyme activities in this study were measured under laboratory conditions with sufficient substrate supplies and suitable environment; therefore, these potential activities may not properly reflect actual microbial enzyme activities under the field conditions. Despite this inherent limitation, we argue that there is a non-significance in measured enzyme activity caused by FTCs, as soil microbes preferentially utilize simple compounds that do not require enzymes for degradation in DOC decomposition enhanced by FTCs (Foster et al., 2016; Gao et al., 2021; Perez-Mon et al., 2020). These results were evidenced by the different DOC quality between the FTC and CON soils (Table 1). The DOC quality indices, SUVA254 and A254 A365, significantly differed between the FTC and CON soils, indicating that complex substrates with high aromaticity and molecular weight remained in the dissolved organic matter after successive FTCs (Berggren et al., 2018; Yang et al., 2019).
A series of multivariate analyses indicated the relationships between DOC characteristics and soil microbial activities influenced by FTCs. The results in the PCA and Pearson's correlation analysis showed that soil CO2 production, influenced by FTCs, was closely related to the quantitative and qualitative changes in DOC (Figs. 2 and 3). Furthermore, as shown in Table 5, MLR analyses identified that those relationships are direct. These results eventually confirmed the first hypothesis that FTCs can change DOC quantity and quality without inhibiting soil microbial activities previously adapted to temperature fluctuations in the Arctic.
4.2 Effects of freeze–thaw cycles on dissolved organic carbon associated with soil structural properties
FTCs caused an increase in microaggregates (53–250 µm) and a corresponding decrease in mineral-associated fractions (< 53 µm), despite low significance levels (Table 3). In other words, the formation of microaggregates by FTCs is likely enhanced by the binding of smaller-sized aggregates rather than the breakdown of larger-sized ones. This could be related to the ice formation in the soil as the ambient temperature drops to −9.0 ± 0.3 ∘C during the FTCs (Fig. S1). Under the freezing phase of FTCs, soil water gradually freezes, but at a microscale level, it can still form thin films of unfrozen water on the surfaces of soil particles. These unfrozen water films can intensively contain dissolved solutes that are charged or excluded during icing, which contribute to stabilizing soil structure (Sletten, 1988; Zhang et al., 2016). These characteristics allow the unfrozen water films to function as binding agents between soil particles, enhancing soil microaggregation after seven successive FTCs. This is further supported by the fact that the extent of the decrease in mineral-associated fractions, despite no significance, is comparable to that of the increase in microaggregates (Table 3).
Soil microaggregates enhanced by FTCs affected DOC quantity and quality mainly through changes in the microbe-mediated mechanism rather than the direct pathway (Fig. 4). Because soil structural dynamics derived from FTCs can be a critical process for soil quality and function (Rabot et al., 2018), soil aggregate formation can improve soil structural stability and govern nutrient cycling and water retention, resulting in enhanced microbial activity (Bird et al., 2000; Yoo et al., 2017; Kim et al., 2021). Consequently, our findings suggest that soil structural improvement by FTCs, at the microaggregate scale, contributes to DOC decomposition by soil microbes, thereby resulting in reduced DOC content and increased DOC aromaticity in the FTC soil (Fig. 4).
Furthermore, the quantitative and qualitative changes in DOC can be attributed to the formation of specific-sized pores by soil microaggregates enhanced by FTCs. We found that the significant difference in soil PSD affected by FTCs was in the small-sized mesopores (Table 5). These pores were strongly related to soil microaggregate formation (Figs. 3 and 5). Our findings indicate that the increase in microaggregate formation by FTCs likely created the corresponding small-sized mesopores through the rearrangement and formation of soil pores (Peng et al., 2015; Dal Ferro et al., 2012; Zaffar and Lu, 2015). Furthermore, such pores are able to hold water surrounding the soil particles (Jim and Ng, 2018; Kim et al., 2021), potentially contributing to water-film development on the soil particle surfaces. As previously mentioned, these water films can have electrical charges and condensed solutes during the freezing periods of the FTCs, serving as binding materials for soil microaggregation (Zhang et al., 2016). Thus, increase in mesopores by enhanced soil microaggregation may permit the dissolved solutes in the soil pore water to be adsorbed and occluded to the soils, thereby decreasing the DOC content in the soil solution of the FTC treatment (Fig. 4).
Multivariate analyses confirmed the second hypothesis about the quantitative and qualitative characteristics of DOC associated with soil structural changes by FTCs. The FTCs led to an increase in soil microaggregate formation and consequent changes in soil microbial activity and pore distribution, accelerating DOC decomposition and decreasing its content in the soil solution (Figs. 2 and 3; Table 5). Our findings contribute to a mechanism-based understanding of the effect of FTCs on DOC properties through systematic measurements on soil biogeochemical properties.
This study demonstrated the responses of organic soils in the Arctic tundra to FTCs, focusing on the changes in DOC characteristics associated with microbial activity and soil physical structure. We found that the following seven variables differed significantly after FTCs: soil CO2 production (CO2), DOC and TDN contents, two DOC quality indices (SUVA254 and A365 A254), microaggregate (53–250 µm) distribution, and small-sized mesopore (0.2–10 µm) proportion. Multivariate statistical analyses, including PCA, Pearson's correlation, and MLR, contributed to the mechanism-based interpretation of how FTCs altered DOC quantity and quality mediated by the changes in microbial activity and soil physical structure. As a result, FTCs altered the DOC quantity and quality with higher CO2, indicating that FTCs affected DOC characteristics without negatively impacting microbial activity. In addition, soil microaggregation enhanced by FTCs and the subsequent increase in soil CO2 production and small-sized pore distribution could promote DOC decomposition, eventually decreasing the DOC content in the soil solution. In conclusion, we elucidated the effects of FTCs on DOC characteristics in the Arctic organic soils of active layer by incorporating soil structural changes and microbial responses. Further study is required to determine how the deeper active layer or ice-rich permafrost thaw under warming would affect the permafrost C dynamics with FTCs.
All data can be provided by the corresponding author upon request.
The supplement related to this article is available online at: https://doi.org/10.5194/tc-17-3101-2023-supplement.
YJK and JYJ planned the campaign; YJK and JK performed the measurement and analyzed the data; YJK wrote the manuscript draft; JYJ and JK reviewed and edited the paper; JYJ acquired the financial support for the project leading to this publication.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This study was supported by the National Research Foundation of Korea funded by the Korean Government (grant no. NRF-2021M1A5A1075508/KOPRI-PN23012). We are grateful to Sungjin Nam, a researcher at the Korea Polar Research Institute, for performing fieldwork and soil sampling.
This research has been supported by the National Research Foundation of Korea (grant no. NRF-2021M1A5A1075508/KOPRI-PN23012).
This paper was edited by Hanna Lee and reviewed by Liam Heffernan and one anonymous referee.
Al-Houri, Z. M., Barber, M. E., Yonge, D. R., Ullman, J. L., and Beutel, M. W.: Impacts of frozen soils on the performance of infiltration treatment facilities, Cold Reg. Sci. Technol., 59, 51–57, https://doi.org/10.1016/j.coldregions.2009.06.002, 2009.
Athmann, M., Kautz, T., Pude, R., and Köpke, U.: Root growth in biopores-evaluation with in situ endoscopy, Plant Soil, 371, 179–190, https://doi.org/10.1007/s11104-013-1673-5, 2013.
Ban, Y., Lei, T., Liu, Z., and Chen, C.: Comparison of rill flow velocity over frozen and thawed slopes with electrolyte tracer method, J. Hydrol., 534, 630–637, https://doi.org/10.1016/j.jhydrol.2016.01.028, 2016.
Berggren, M., Klaus, M., Selvam, B. P., Ström, L., Laudon, H., Jansson, M., and Karlsson, J.: Quality transformation of dissolved organic carbon during water transit through lakes: contrasting controls by photochemical and biological processes, Biogeosciences, 15, 457–470, https://doi.org/10.5194/bg-15-457-2018, 2018.
Bird, N. R. A., Perrier, E., and Rieu, M.: The water retention function for a model of soil structure with pore and solid fractal distributions, Eur. J. Soil Sci., 51, 55–63, https://doi.org/10.1046/j.1365-2389.2000.00278.x, 2000.
Callaghan, T., Körner, C., Heal, O., Lee, S., and Cornelissen, J.: Scenarios for ecosystem responses to global change, Global Change in Europe's Cold Regions, Ecosystem Research Report, 65–134, http://urn.kb.se/resolve?urn=urn:nbn:se:polar:diva-4698 (last access: 17 July 2023), 1998.
Čapek, P., Diáková, K., Dickopp, J. E., Bárta, J., Wild, B., Schnecker, J., Alves, R. J. E., Aiglsdorfer, S., Guggenberger, G., Gentsch, N., Hugelius, G., Lashchinsky, N., Gittel, A., Schleper, C., Mikutta, R., Palmtag, J., Shibistova, O., Urich, T., Richter, A., and Šantrůčková, H.: The effect of warming on the vulnerability of subducted organic carbon in arctic soils, Soil Biol. Biochem., 90, 19–29, https://doi.org/10.1016/j.soilbio.2015.07.013, 2015.
Dal Ferro, N., Delmas, P., Duwig, C., Simonetti, G., and Morari, F.: Coupling X-ray microtomography and mercury intrusion porosimetry to quantify aggregate structures of a cambisol under different fertilisation treatments, Soil Tillage Res., 119, 13–21, https://doi.org/10.1016/j.still.2011.12.001, 2012.
da Silva, A. and Kay, B.: Estimating the least limiting water range of soils from properties and management, Soil Sci. Soc. Am. J., 61, 877–883, 1997.
Davidson, E. A., Belk, E., and Boone, R. D.: Soil water content and temperature as independent or confounded factors controlling soil respiration in a temperate mixed hardwood forest, Glob. Change Biol., 4, 217–227, https://doi.org/10.1046/j.1365-2486.1998.00128.x, 1998.
Dettmann, U., Bechtold, M., Frahm, E., and Tiemeyer, B.: On the applicability of unimodal and bimodal van Genuchten-Mualem based models to peat and other organic soils under evaporation conditions, J. Hydrol., 515, 103–115, https://doi.org/10.1016/j.jhydrol.2014.04.047, 2014.
Estop-Aragonés, C., Olefeldt, D., Abbott, B. W., Chanton, J. P., Czimczik, C. I., Dean, J. F., Egan, J. E., Gandois, L., Garnett, M. H., Hartley, I. P., Hoyt, A., Lupascu, M., Natali, S. M., O'Donnell, J. A., Raymond, P. A., Tanentzap, A. J., Tank, S. E., Schuur, E. A. G., Turetsky, M., and Anthony, K. W.: Assessing the Potential for Mobilization of Old Soil Carbon After Permafrost Thaw: A Synthesis of 14C Measurements From the Northern Permafrost Region, Global Biogeochem. Cy., 34, 1–26, https://doi.org/10.1029/2020GB006672, 2020.
Feng, X., Nielsen, L. L., and Simpson, M. J.: Responses of soil organic matter and microorganisms to freeze-thaw cycles, Soil Biol. Biochem., 39, 2027–2037, https://doi.org/10.1016/j.soilbio.2007.03.003, 2007.
Førland, E. J., Benestad, R., Hanssen-Bauer, I., Haugen, J. E., and Skaugen, T. E.: Temperature and Precipitation Development at Svalbard 1900–2100, Adv. Meteorol., 2011, 1–14, https://doi.org/10.1155/2011/893790, 2011.
Foster, A., Jones, D. L., Cooper, E. J., and Roberts, P.: Freeze–thaw cycles have minimal effect on the mineralisation of low molecular weight, dissolved organic carbon in Arctic soils, Polar Biol., 39, 2387–2401, https://doi.org/10.1007/s00300-016-1914-1, 2016.
Freppaz, M., Williams, B. L., Edwards, A. C., Scalenghe, R., and Zanini, E.: Simulating soil freeze/thaw cycles typical of winter alpine conditions: Implications for N and P availability, Appl. Soil Ecol., 35, 247–255, https://doi.org/10.1016/j.apsoil.2006.03.012, 2007.
Gao, D., Zhang, L., Liu, J., Peng, B., Fan, Z., Dai, W., Jiang, P., and Bai, E.: Responses of terrestrial nitrogen pools and dynamics to different patterns of freeze-thaw cycle: A meta-analysis, Glob. Change Biol., 24, 2377–2389, https://doi.org/10.1111/gcb.14010, 2018.
Gao, D., Bai, E., Yang, Y., Zong, S., and Hagedorn, F.: A global meta-analysis on freeze-thaw effects on soil carbon and phosphorus cycling, Soil Biol. Biochem., 159, 108283, https://doi.org/10.1016/j.soilbio.2021.108283, 2021.
Gao, L., Zhou, Z., Reyes, A. V., and Guo, L.: Yields and Characterization of Dissolved Organic Matter From Different Aged Soils in Northern Alaska, J. Geophys. Res.-Biogeo., 123, 2035–2052, https://doi.org/10.1029/2018JG004408, 2018.
Groffman, P. M., Hardy, J. P., Fashu-Kanu, S., Driscoll, C. T., Cleavitt, N. L., Fahey, T. J., and Fisk, M. C.: Snow depth, soil freezing and nitrogen cycling in a northern hardwood forest landscape, Biogeochemistry, 102, 223–238, https://doi.org/10.1007/s10533-010-9436-3, 2011.
Grogan, P., Michelsen, A., Ambus, P., and Jonasson, S.: Freeze-thaw regime effects on carbon and nitrogen dynamics in sub-arctic heath tundra mesocosms, Soil Biol. Biochem., 36, 641–654, https://doi.org/10.1016/j.soilbio.2003.12.007, 2004.
Hall, K. and André, M. F.: Rock thermal data at the grain scale: Applicability to granular disintegration in cold environments, Earth Surf. Process. Landf., 28, 823–836, https://doi.org/10.1002/esp.494, 2003.
Han, Z., Deng, M., Yuan, A., Wang, J., Li, H., and Ma, J.: Vertical variation of a black soil's properties in response to freeze-thaw cycles and its links to shift of microbial community structure, Sci. Total Environ., 625, 106–113, https://doi.org/10.1016/j.scitotenv.2017.12.209, 2018.
Heal, O., Broll, G., Hooper, D., McConnel, J., Webb, N., and Wookey, P.: Impacts of global change on tundra soil biology, Global Change in Europe's Cold Regions, Ecosystem Research Report, 65–134, http://urn.kb.se/resolve?urn=urn:nbn:se:polar:diva-4680 (last access: 17 July 2023), 1998.
Henry, H. A. L.: Climate change and soil freezing dynamics: Historical trends and projected changes, Climatic Change, 87, 421–434, https://doi.org/10.1007/s10584-007-9322-8, 2008.
Henry, H. A. L.: Plant and Microbe Adaptations to Cold in a Changing World, Proceedings from Plant and Microbe Adaptations to Cold, Springer New York, NY, 17–28, https://doi.org/10.1007/978-1-4614-8253-6, 2013.
IPCC: Climate Change 2014: Synthesis Report, Contribution of Working group to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Team, C. W., Pachauri, R. K., and Meyer, L. A., Geneva, Switzerland, https://doi.org/10.1177/0002716295541001010, 2014.
Jim, C. Y. and Ng, Y. Y.: Porosity of roadside soil as indicator of edaphic quality for tree planting, Ecol. Eng., 120, 364–374, https://doi.org/10.1016/j.ecoleng.2018.06.016, 2018.
Kameyama, K., Miyamoto, T., Shiono, T., and Shinogi, Y.: Influence of Sugarcane Bagasse-derived Biochar Application on Nitrate Leaching in Calcaric Dark Red Soil, J. Environ. Qual., 41, 1131–1137, https://doi.org/10.2134/jeq2010.0453, 2012.
Kim, H. M., Lee, M. J., Jung, J. Y., Hwang, C. Y., Kim, M., Ro, H. M., Chun, J., and Lee, Y. K.: Vertical distribution of bacterial community is associated with the degree of soil organic matter decomposition in the active layer of moist acidic tundra, J. Microbiol., 54, 713–723, https://doi.org/10.1007/s12275-016-6294-2, 2016.
Kim, Y. J. and Yoo, G.: Suggested key variables for assessment of soil quality in urban roadside tree systems, J. Soils Sediments, 21, 2130–2140, https://doi.org/10.1007/s11368-020-02827-5, 2021.
Kim, Y. J., Hyun, J., Yoo, S. Y., and Yoo, G.: The role of biochar in alleviating soil drought stress in urban roadside greenery, Geoderma, 404, 115223, https://doi.org/10.1016/j.geoderma.2021.115223, 2021.
Kim, Y. J., Laffly, D., Kim, S. eun, Nilsen, L., Chi, J., Nam, S., Lee, Y. B., Jeong, S., Mishra, U., Lee, Y. K., and Jung, J. Y.: Chronological changes in soil biogeochemical properties of the glacier foreland of Midtre Lovénbreen, Svalbard, attributed to soil-forming factors, Geoderma, 415, 115777, https://doi.org/10.1016/j.geoderma.2022.115777, 2022.
Koponen, H. T. and Bååth, E.: Soil bacterial growth after a freezing/thawing event, Soil Biol. Biochem., 100, 229–232, https://doi.org/10.1016/j.soilbio.2016.06.029, 2016.
Kreyling, J., Beierkuhnlein, C., Pritsch, K., Schloter, M., and Jentsch, A.: Recurrent soil freeze-thaw cycles enhance grassland productivity, New Phytol., 177, 938–945, https://doi.org/10.1111/j.1469-8137.2007.02309.x, 2008.
Kuzyakov, Y. and Domanski, G.: Carbon input by plants into the soil. Review, J. Plant Nutr. Soil Sci., 163, 421–431, https://doi.org/10.1002/1522-2624(200008)163:4<421::AID-JPLN421>3.0.CO;2-R, 2000.
Kwon, M. J., Haraguchi, A., and Kang, H.: Long-term water regime differentiates changes in decomposition and microbial properties in tropical peat soils exposed to the short-term drought, Soil Biol. Biochem., 60, 33–44, https://doi.org/10.1016/j.soilbio.2013.01.023, 2013.
Larsen, K. S., Jonasson, S., and Michelsen, A.: Repeated freeze-thaw cycles and their effects on biological processes in two arctic ecosystem types, Appl. Soil Ecol., 21, 187–195, https://doi.org/10.1016/S0929-1393(02)00093-8, 2002.
Lehrsch, G. A.: Freeze-thaw cycles increase near-surface aggregate stability, Soil Sci., 163, 63–70, https://doi.org/10.1097/00010694-199801000-00009, 1998.
Li, G. Y. and Fan, H. M.: Effect of Freeze-Thaw on Water Stability of Aggregates in a Black Soil of Northeast China, Pedosphere, 24, 285–290, https://doi.org/10.1016/S1002-0160(14)60015-1, 2014.
Liang, A., Zhang, Y., Zhang, X., Yang, X., McLaughlin, N., Chen, X., Guo, Y., Jia, S., Zhang, S., Wang, L., and Tang, J.: Investigations of relationships among aggregate pore structure, microbial biomass, and soil organic carbon in a Mollisol using combined non-destructive measurements and phospholipid fatty acid analysis, Soil Tillage Res., 185, 94–101, https://doi.org/10.1016/j.still.2018.09.003, 2019.
Liao, X., Kang, H., Haidar, G., Wang, W., and Malghani, S.: The impact of biochar on the activities of soil nutrients acquisition enzymes is potentially controlled by the pyrolysis temperature: A meta-analysis, Geoderma, 411, 115692, https://doi.org/10.1016/j.geoderma.2021.115692, 2022.
Likos, W. J., Lu, N., and Godt, J. W.: Hysteresis and Uncertainty in Soil Water-Retention Curve Parameters, J. Geotech. Geoenviron., 140, 1–11, https://doi.org/10.1061/(asce)gt.1943-5606.0001071, 2014.
Lim, A. G., Loiko, S. V., Kuzmina, D. M., Krickov, I. V., Shirokova, L. S., Kulizhsky, S. P., Vorobyev, S. N., and Pokrovsky, O. S.: Dispersed ground ice of permafrost peatlands: Potential unaccounted carbon, nutrient and metal sources, Chemosphere, 266, 128953, https://doi.org/10.1016/j.chemosphere.2020.128953, 2021.
Lipson, D. A. and Monson, R. K.: Plant-microbe competition for soil amino acids in the alpine tundra: Effects of freeze-thaw and dry-rewet events, Oecologia, 113, 406–414, https://doi.org/10.1007/s004420050393, 1998.
Lipson, D. A. and Schmidt, S. K.: Seasonal changes in an alpine soil bacterial community in the Colorado Rocky Mountains, Appl. Environ. Microbiol., 70, 2867–2879, https://doi.org/10.1128/AEM.70.5.2867-2879.2004, 2004.
Lipson, D. A., Schmidt, S. K., and Monson, R. K.: Carbon availability and temperature control the post-snowmelt decline in alpine soil microbial biomass, Soil Biol. Biochem., 32, 441–448, https://doi.org/10.1016/S0038-0717(99)00068-1, 2000.
Liu, B., Ma, R., and Fan, H.: Evaluation of the impact of freeze-thaw cycles on pore structure characteristics of black soil using X-ray computed tomography, Soil Tillage Res., 206, 104810, https://doi.org/10.1016/j.still.2020.104810, 2021.
Lu, Y., Liu, S., Zhang, Y., Wang, L., and Li, Z.: Hydraulic conductivity of gravelly soils with various coarse particle contents subjected to freeze–thaw cycles, J. Hydrol., 598, 126302, https://doi.org/10.1016/j.jhydrol.2021.126302, 2021.
Ma, R., Jiang, Y., Liu, B., and Fan, H.: Effects of pore structure characterized by synchrotron-based micro-computed tomography on aggregate stability of black soil under freeze-thaw cycles, Soil Tillage Res., 207, 104855, https://doi.org/10.1016/j.still.2020.104855, 2021.
Maikhuri, R. K. and Rao, K. S.: Soil quality and soil health: A review, Int. J. Ecol. Environ. Sci., 38, 19–37, 2012.
Männistö, M. K., Tiirola, M., and Häggblom, M. M.: Effect of freeze-thaw cycles on bacterial communities of Arctic tundra soil, Microb. Ecol., 58, 621–631, https://doi.org/10.1007/s00248-009-9516-x, 2009.
Matzner, E. and Borken, W.: Do freeze-thaw events enhance C and N losses from soils of different ecosystems? A review, Eur. J. Soil Sci., 59, 274–284, https://doi.org/10.1111/j.1365-2389.2007.00992.x, 2008.
Oztas, T. and Fayetorbay, F.: Effect of freezing and thawing processes on soil aggregate stability, Catena, 52, 1–8, https://doi.org/10.1016/S0341-8162(02)00177-7, 2003.
Peng, X., Horn, R., and Hallett, P.: Soil structure and its functions in ecosystems: Phase matter & scale matter, Soil Tillage Res., 146, 1–3, https://doi.org/10.1016/j.still.2014.10.017, 2015.
Perez-Mon, C., Frey, B., and Frossard, A.: Functional and Structural Responses of Arctic and Alpine Soil Prokaryotic and Fungal Communities Under Freeze-Thaw Cycles of Different Frequencies, Front. Microbiol., 11, 1–14, https://doi.org/10.3389/fmicb.2020.00982, 2020.
Rabot, E., Wiesmeier, M., Schlüter, S., and Vogel, H. J.: Soil structure as an indicator of soil functions: A review, Geoderma, 314, 122–137, https://doi.org/10.1016/j.geoderma.2017.11.009, 2018.
Raich, J. and Schlesinger, W.: The global carbon dioxdie flux in soil respiration and its relationship to vegetation and climate, Tellus B, 44, 81–99, 1992.
Rantanen, M., Karpechko, A. Y., Lipponen, A., Nordling, K., Hyvärinen, O., Ruosteenoja, K., Vihma, T., and Laaksonen, A.: The Arctic has warmed nearly four times faster than the globe since 1979, Commun. Earth Environ., 3, 1–10, https://doi.org/10.1038/s43247-022-00498-3, 2022.
Royer, A., Picard, G., Vargel, C., Langlois, A., Gouttevin, I., and Dumont, M.: Improved Simulation of Arctic Circumpolar Land Area Snow Properties and Soil Temperatures, Front. Earth Sci., 9, 1–19, https://doi.org/10.3389/feart.2021.685140, 2021.
Sander, T. and Gerke, H. H.: Preferential Flow Patterns in Paddy Fields Using a Dye Tracer, Vadose Zone J., 6, 105–115, https://doi.org/10.2136/vzj2006.0035, 2007.
Sawicka, J. E., Robador, A., Hubert, C., Jørgensen, B. B., and Brüchert, V.: Effects of freeze-thaw cycles on anaerobic microbial processes in an Arctic intertidal mud flat, ISME J., 4, 585–594, https://doi.org/10.1038/ismej.2009.140, 2010.
Schimel, J. P. and Clein, J. S.: Microbial response to freeze-thaw cycles in tundra and taiga soils, Soil Biol. Biochem., 28, 1061–1066, https://doi.org/10.1016/0038-0717(96)00083-1, 1996.
Schimel, J. P. and Mikan, C.: Changing microbial substrate use in Arctic tundra soils through a freeze-thaw cycle, Soil Biol. Biochem., 37, 1411–1418, https://doi.org/10.1016/j.soilbio.2004.12.011, 2005.
Šimůnek, J., van Genuchten, T., and Kodešová, R.: HYDRUS Software Applications to Subsurface Flow and Contaminant Transport Problems, in: Proceedings of the 4th International Conference, 404, ISBN 978-80-213-2380-3, 2013.
Sinsabaugh, R. L.: Phenol oxidase, peroxidase and organic matter dynamics of soil, Soil Biol. Biochem., 42, 391–404, https://doi.org/10.1016/j.soilbio.2009.10.014, 2010.
Sjursen, H., Michelsen, A., and Holmstrup, M.: Effects of freeze-thaw cycles on microarthropods and nutrient availability in a sub-Arctic soil, Appl. Soil Ecol., 28, 79–93, https://doi.org/10.1016/j.apsoil.2004.06.003, 2005.
Sletten, R. S.: The formation of pedogenic carbonates on Svalbard: The influence of cold temperatures and freezing, Permafrost, 5, 467–472, 1988.
Song, Y., Zou, Y., Wang, G., and Yu, X.: Altered soil carbon and nitrogen cycles due to the freeze-thaw effect: A meta-analysis, Soil Biol. Biochem., 109, 35–49, https://doi.org/10.1016/j.soilbio.2017.01.020, 2017.
Troy, S. M., Lawlor, P. G., O' Flynn, C. J., and Healy, M. G.: Impact of biochar addition to soil on greenhouse gas emissions following pig manure application, Soil Biol. Biochem., 60, 173–181, https://doi.org/10.1016/j.soilbio.2013.01.019, 2013.
Viklander, P.: Permeability and volume changes in till due to cyclic freeze/thaw, Can. Geotech. J., 35, 471–477, https://doi.org/10.1139/t98-015, 1998.
Walker, V. K., Palmer, G. R., and Voordouw, G.: Freeze-thaw tolerance and clues to the winter survival of a soil community, Appl. Environ. Microbiol., 72, 1784–1792, https://doi.org/10.1128/AEM.72.3.1784-1792.2006, 2006.
Wang, E., Cruse, R. M., Chen, X., and Daigh, A.: Effects of moisture condition and freeze/thaw cycles on surface soil aggregate size distribution and stability, Can. J. Soil Sci., 92, 529–536, https://doi.org/10.4141/CJSS2010-044, 2012.
White, D. M., Garland, D. S., Ping, C. L., and Michaelson, G.: Characterizing soil organic matter quality in arctic soil by cover type and depth, Cold Reg. Sci. Technol., 38, 63–73, https://doi.org/10.1016/j.coldregions.2003.08.001, 2004.
Xu, C., Guo, L., Dou, F., and Ping, C. L.: Potential DOC production from size-fractionated Arctic tundra soils, Cold Reg. Sci. Technol., 55, 141–150, https://doi.org/10.1016/j.coldregions.2008.08.001, 2009.
Yang, X. Y., Chang, K. H., Kim, Y. J., Zhang, J., and Yoo, G.: Effects of different biochar amendments on carbon loss and leachate characterization from an agricultural soil, Chemosphere, 226, 625–635, https://doi.org/10.1016/j.chemosphere.2019.03.085, 2019.
Yi, Y., Kimball, J. S., Rawlins, M. A., Moghaddam, M., and Euskirchen, E. S.: The role of snow cover affecting boreal-arctic soil freeze-thaw and carbon dynamics, Biogeosciences, 12, 5811–5829, https://doi.org/10.5194/bg-12-5811-2015, 2015.
Yoo, G., Kim, H., and Choi, J. Y.: Soil Aggregate Dynamics Influenced by Biochar Addition using the 13C Natural Abundance Method, Soil Sci. Soc. Am. J., 81, 612–621, https://doi.org/10.2136/sssaj2016.09.0313, 2017.
Yoo, S. Y., Kim, Y. J., and Yoo, G.: Understanding the role of biochar in mitigating soil water stress in simulated urban roadside soil, Sci. Total Environ., 738, 139798, https://doi.org/10.1016/j.scitotenv.2020.139798, 2020.
Zaffar, M. and Lu, S. G.: Pore size distribution of clayey soils and its correlation with soil organic matter, Pedosphere, 25, 240–249, https://doi.org/10.1016/S1002-0160(15)60009-1, 2015.
Zhang, Z., Ma, W., Feng, W., Xiao, D., and Hou, X.: Reconstruction of Soil Particle Composition During Freeze-Thaw Cycling: A Review, Pedosphere, 26, 167–179, https://doi.org/10.1016/S1002-0160(15)60033-9, 2016.