Productivity and Seasonality Drive Total Soil Respiration in Semi-Arid Juniper Woodlands ( Juniperus thurifera L., Southern Spain)

: We analyzed the relationship between forest productivity (joint effect of forest maturity and soil quality) and total soil respiration ( µ mol CO 2 m − 2 s − 1 ) in semi-arid juniper woodlands (young woodlands growing in abandoned farmlands with deeper soils, and mature woodlands in lithic soils), and investigated the seasonal variation in soil CO 2 efﬂux as a function of soil temperature and the soil water content. We measured the soil CO 2 efﬂux from twelve cylinders in the soil over a three-year period using LI-6400 equipment. The results show that, in the more productive site (young woodland), soil CO 2 efﬂux was higher due to greater respiration, mainly in the driest periods. Soil respiration followed a seasonal trend, being higher in spring and decreasing in cold periods. In both juniper woodlands and especially in the older forest, the CO 2 efﬂux rates were low (<2.5 for Q 10 ), typical of slow-growing species. Soil respiration was controlled by soil temperature without drought and in the temperate-warm season, whereas respiration showed sensitivity to soil water content in periods when edaphic humidity was low (but only in the more productive, young forest, which seemed to show better adaptation to drought), and under high soil moisture (soil water > 25%) for both woodlands, coinciding with warm temperatures in the spring. This period also corresponded to the highest CO 2 efﬂux recorded in both woodlands. The accumulation of organic C seems to also be important to maintain elevated soil respiration in summer, especially in young woodlands. Thus, apart from microclimatic conditions, factors related to productivity regulate respiratory activity.


Introduction
The soil-to-atmosphere CO 2 efflux constitutes a large part of the terrestrial C cycle and is the second-largest cause of terrestrial carbon flux [1,2]. The soil CO 2 efflux is defined as total soil respiration because it is mainly produced by underground respiration [3]. Total soil respiration is the sum of autotrophic respiration by plant roots and heterotrophic respiration by soil microbes and organisms [4,5], with the contribution of root respiration to total soil respiration being 0.54, 0.41, 0.56, 0.22, and 0.50 for the boreal evergreen needleleaf, temperate deciduous broadleaf, temperate evergreen needleleaf, temperate evergreen broadleaf, and tropical evergreen broadleaf forests, respectively [6].
Mill. is only found in Spain and the French Pyrenees [30], and due to its rarity and slow growth, it has been included in catalogs of endangered species in Spain. The phytoclimatic area of this evergreen conifer is correlated with severe climates [31], such as those predicted in climate change scenarios [32]. For this reason, this species is utilized for forest restoration in sites with low productivity, such as semi-arid areas [26]. Thermophilic juniper woodlands occupy the driest areas of the distribution of Juniperus thurifera in Spain [30]. The species here constitute relict forests [33] that can be classified into two levels of maturity and productivity [34]: (i) old or mature woodlands that generally contain trees that are over 150 years old that grow in stony soils with a low tree density (<100 trees ha −1 ), and (ii) young woodlands, with trees that are less than 100 years old growing in abandoned farmlands (i.e., better soils with a higher water storage capacity). These forests have higher growth rates, productivity, and tree density [34,35].
After analyzing our study area, we hypothesized that total soil respiration would differ among the woodlands, and in addition to forest productivity, these differences would also be controlled by microclimatic factors. Consequently, the specific aims were the following: (i) to analyze the effects of the juniper woodland type (young and mature woodlands, effects of forest "productivity") on soil CO 2 efflux, (ii) to study the seasonal variation in soil CO 2 efflux in the two juniper woodlands under a semi-arid Mediterranean climate ("seasonality"), and (iii) to evaluate the effects of microclimatic conditions (soil temperature and soil water content) on soil CO 2 efflux for both juniper woodlands throughout the seasons.

Study Area
This research was carried out in the Spanish juniper woodlands of the El Campo de Montiel Mountains (Southern Spain; Figure 1). summers typical of semi-arid climates [29]. The Spanish juniper (Juniperus thurifera L., C pressaceae) is an endemic dioecious species of the Western Mediterranean. J. thurifera v hispanica Mill. is only found in Spain and the French Pyrenees [30], and due to its rari and slow growth, it has been included in catalogs of endangered species in Spain. T phytoclimatic area of this evergreen conifer is correlated with severe climates [31], su as those predicted in climate change scenarios [32]. For this reason, this species is utiliz for forest restoration in sites with low productivity, such as semi-arid areas [26]. Therm philic juniper woodlands occupy the driest areas of the distribution of Juniperus thurife in Spain [30]. The species here constitute relict forests [33] that can be classified into tw levels of maturity and productivity [34]: (i) old or mature woodlands that generally co tain trees that are over 150 years old that grow in stony soils with a low tree density (<1 trees ha −1 ), and (ii) young woodlands, with trees that are less than 100 years old growi in abandoned farmlands (i.e., better soils with a higher water storage capacity). These fo ests have higher growth rates, productivity, and tree density [34,35].
After analyzing our study area, we hypothesized that total soil respiration wou differ among the woodlands, and in addition to forest productivity, these differenc would also be controlled by microclimatic factors. Consequently, the specific aims we the following: (i) to analyze the effects of the juniper woodland type (young and matu woodlands, effects of forest "productivity") on soil CO2 efflux, (ii) to study the season variation in soil CO2 efflux in the two juniper woodlands under a semi-arid Mediterrane climate ("seasonality"), and (iii) to evaluate the effects of microclimatic conditions (s temperature and soil water content) on soil CO2 efflux for both juniper woodlan throughout the seasons.

Study Area
This research was carried out in the Spanish juniper woodlands of the El Campo Montiel Mountains (Southern Spain; Figure 1).   [30]. UTM coordinates, ETRS89. These juniper woodlands are thermophilic according to [30]. The climate is semi-arid cold (BSk type, KÖPPEN [36]) with extreme temperatures of between 43 and −21 • C (data from "El Bonillo" station: 38 • 57 N, 01 • 09 W, 1068 m.a.s.l., period of 30 years, located within 7 km of the research area, and owned by the State Meteorological Agency of Spain). The average annual temperature is 12.8 • C, and annual rainfall corresponds to a dry ombroclimate (450 mm). Soils in which the total soil CO 2 efflux was measured were selected from the two types of thermophilic juniper woodland (Table 1, and Figure 2): (i) mature woodland with an average age of 170 years, growing in Leptosols (soil depth < 15 cm), and (ii) young woodland with a mean age of 66 years, growing in abandoned farmlands with deeper soils (Cambisols, with a mean depth of 42 cm). In each type of juniper woodland, a representative and homogenous experimental area ("plots", approximately 1 hectare) was selected. The areas were less than 1 km apart, ensuring that both woodlands are growing under the same climatology, as well as similar altitude (1050 m.a.s.l.), slope (gentle slopes, 5 • ), and exposure (paramo formation, flat surfaces). Referring to vegetation characteristics (Table 1, vegetation variables have been collected from the study of Garcia [33]), biomass stocks and productivity were scarce in both woodlands. However, the young site has higher productivity. Juniper trees in the young woodland area have slow growth, and trees in the mature woodland area represent very-slow-growing species. Litterfall and the LAI are also higher in young woodlands due to the superior productivity (joint effect of age and quality of soil). The two juniper woodlands contained lime-rich substrates (Jurassic dolomites) with a basic pH (Table 1, soil characteristics). The soils in the young woodland area had a greater clay content and water storage capacity (137.1 vs. 20.5 mm). The bulk density was slightly greater in the mature woodland. The average annual temperature is 12.8 °C, and annual rainfall corresponds to a dry ombroclimate (450 mm). Soils in which the total soil CO2 efflux was measured were selected from the two types of thermophilic juniper woodland (Table 1, and Figure 2): (i) mature woodland with an average age of 170 years, growing in Leptosols (soil depth < 15 cm), and (ii) young woodland with a mean age of 66 years, growing in abandoned farmlands with deeper soils (Cambisols, with a mean depth of 42 cm). In each type of juniper woodland, a representative and homogenous experimental area ("plots", approximately 1 hectare) was selected. The areas were less than 1 km apart, ensuring that both woodlands are growing under the same climatology, as well as similar altitude (1050 m.a.s.l.), slope (gentle slopes, 5°), and exposure (paramo formation, flat surfaces). Referring to vegetation characteristics (Table 1, vegetation variables have been collected from the study of Garcia [33]), biomass stocks and productivity were scarce in both woodlands. However, the young site has higher productivity. Juniper trees in the young woodland area have slow growth, and trees in the mature woodland area represent veryslow-growing species. Litterfall and the LAI are also higher in young woodlands due to the superior productivity (joint effect of age and quality of soil). The two juniper woodlands contained lime-rich substrates (Jurassic dolomites) with a basic pH ( Table 1, soil characteristics). The soils in the young woodland area had a greater clay content and water

Measurement of Soil CO 2 Efflux
Total soil CO 2 efflux (µmol CO 2 m −2 s −1 ) was measured in situ on 17 sampling days, spread over the 4 seasons (2 sampling days in winter, 3 in summer and autumn, and 9 in spring, the most important growing season). Within each woodland area, six metallic collars (cylinder: 10.5 cm in diameter and 5 cm in height) were installed on the soil between the trunk and the projection limit of the crown in six juniper trees ( Figure 3). Thus, a total of 12 metal cylinders were installed on the soil for the respiration measurement (6 juniper trees in each of the mature and young forests).  The total soil CO2 efflux was measured using LI-6400 CO2 measurement equipment (LI-COR, Inc., Lincoln, NE, USA), to which an LI-6400-09 chamber was attached (closed system). For each CO2 efflux measurement, three repeated measurements were performed, and the mean was the value used for data analysis. Thus, the total dataset included 204 measures. The measurements were conducted from 9:00 to 15:00 to minimize the effects of daytime soil temperature fluctuations and rainless days [10]. Measurements were performed by reducing the CO2 concentration inside the chamber and then letting it increase to an upper concentration limit. Before starting each measurement, the ambient CO2 concentration was used as the target. The air CO2 concentration ranged between 370 and 421 p.p.m. during the research. The air flow (µmol air s −1 ) pumping down was fixed depending on the efflux level (200-700 for low and high efflux levels, respectively). The rate of change of CO2 inside the chamber was calculated in three consecutive measurement cycles.

Microclimatic Conditions and Organic Matter in Soils: Measurement of the Soil Temperature (Ts, °C), Soil Water Content (Sw, %), and Soil Organic C (%)
The soil temperature (°C) was measured at a depth of 7 cm below the ground surface next to each collar using a temperature probe. Temperature measurements were recorded with HOBO sensors and data loggers (Onset Computers, Bourne, MA, USA). The soil water content (volumetric moisture, %) was measured with WatchDog sensors and data loggers (Spectrum Technologies, Aurora, IL., USA) close to the collars (at the depth of 7 cm). Finally, to determine the soil organic C, four samples of topsoil (0-10 cm) were collected in each woodland type in each season for the two years (in the middle of the season) using a soil core (5 cm inside diameter). To minimize the soil variability, each soil sample was made up of five random subsamples collected between the trunk and canopy edge. Soil samples were taken from the upper 10 cm layer after the litter had been removed (this layer presents the greatest microbial activity because of the shallowness of soils). In the laboratory, soil samples were sieved (<2 mm) and stored at 4 °C until analysis. The total soil organic carbon (%) content was determined by oxidation with K2CrO7 in an acid medium and titration of the excess dichromate with (NH4)2Fe(SO4)2 [39]. Samples were analyzed by assessing three replicates in the laboratory, and the mean value was used as the sampling data. The total soil CO 2 efflux was measured using LI-6400 CO 2 measurement equipment (LI-COR, Inc., Lincoln, NE, USA), to which an LI-6400-09 chamber was attached (closed system). For each CO 2 efflux measurement, three repeated measurements were performed, and the mean was the value used for data analysis. Thus, the total dataset included 204 measures. The measurements were conducted from 9:00 to 15:00 to minimize the effects of daytime soil temperature fluctuations and rainless days [10]. Measurements were performed by reducing the CO 2 concentration inside the chamber and then letting it increase to an upper concentration limit. Before starting each measurement, the ambient CO 2 concentration was used as the target. The air CO 2 concentration ranged between 370 and 421 p.p.m. during the research. The air flow (µmol air s −1 ) pumping down was fixed depending on the efflux level (200-700 for low and high efflux levels, respectively). The rate of change of CO 2 inside the chamber was calculated in three consecutive measurement cycles. The soil temperature ( • C) was measured at a depth of 7 cm below the ground surface next to each collar using a temperature probe. Temperature measurements were recorded with HOBO sensors and data loggers (Onset Computers, Bourne, MA, USA). The soil water content (volumetric moisture, %) was measured with WatchDog sensors and data loggers (Spectrum Technologies, Aurora, IL, USA) close to the collars (at the depth of 7 cm). Finally, to determine the soil organic C, four samples of topsoil (0-10 cm) were collected in each woodland type in each season for the two years (in the middle of the season) using a soil core (5 cm inside diameter). To minimize the soil variability, each soil sample was made up of five random subsamples collected between the trunk and canopy edge. Soil samples were taken from the upper 10 cm layer after the litter had been removed (this layer presents the greatest microbial activity because of the shallowness of soils). In the laboratory, soil samples were sieved (<2 mm) and stored at 4 • C until analysis. The total soil organic carbon (%) content was determined by oxidation with K 2 CrO 7 in an acid medium and titration of the excess dichromate with (NH 4 ) 2 Fe(SO 4 ) 2 [39]. Samples were analyzed by assessing three replicates in the laboratory, and the mean value was used as the sampling data.

Data Analysis
The influence of main factors on soil respiration (Rs, µmol CO 2 m −2 s −1 ) was analyzed using a nested linear mixed model (LMM) with woodland type as the fixed factor (two levels: young and mature woodland), measurement day (date) as the repeated fixed factor (17 levels, to study the seasonality in Rs), and collar nested within woodland as the random factor. Random variance components were calculated using restricted maximum likelihood estimation (RMLE). Two variance structures were considered (compound symmetry and continuous AR1). Bayesian Information Criterion (BIC) and Akaike Information Criterion (AIC) were used as criteria for the selection of the best model. Mean values were compared using Fisher's LSD post-hoc test (α = 0.05). Linear mixed models were made using the nlme package in R [40]. It the model, "woodland type" (W) was a confounded effect of both maturity (age) and soil quality (depth and water retention capacity). For this, "woodland" represents the productivity of the juniper sites.
To analyze the relationships among the total soil CO 2 efflux (Rs, µmol CO 2 m −2 s −1 ), soil temperature (Ts, • C), and/or soil water content (Sw, %) throughout the seasons included in the study period, the following bivariate model was fitted [41]: In Model 1, the effect of the juniper woodland type on the soil CO 2 efflux was included through the dummy variable "W" (W = 0 for mature woodland, and W = 1 for young woodland). Model 1 shows the exponential relationship between cell respiration and temperature [42] and adds the influence of soil moisture on this relationship. In the model, Ts was standardized by 10 • C (Ts/10) to obtain the Q 10 value [43]. When the soil water content (Sw) was not significant and was therefore removed from Model 1, we obtained a univariate model to describe the sensitivity of Rs to soil temperature, described by the coefficient Q 10 (temperature coefficient): where α and α are the coefficients of temperature sensitivity of soil respiration for mature and young woodlands, respectively. The temperature coefficient, Q 10 (rate of variation of Rs when the soil temperature increased by 10 • C), was calculated from the univariate function as [41,43]: Q 10 = e (α+α W) , α for mature woodland and α + α for young woodland, if W was significant. Thus, significant equations in function of the woodland type could be formulated in their exponential forms [5,44]: Rs = Rr × Q 10 Ts−Tr 10 , where Rr is the respiration at the reference temperature (Tr), which was fixed at 0 • C in our study. Models and Q 10 values were estimated for the short term, i.e., seasons ("seasonality" of Q 10 ), and the long term (yearly, i.e., the annual whole data for all seasons without division by years [16]) for each woodland area.
The selection of significant parameters in Model 1 was carried out by means of logarithmic transformation (linearization of the exponential function [44,45]) and stepwise regression [46]: log(Rs) This corresponds with the comparison of two regression lines as well as the study of significant differences in Q 10 between the two woodlands during the seasons [44]. In all regressions fitted, coefficients with a value of p < 0.05 were considered significant [46]. The linear regressions were assessed according to the F ratio (p < 0.05) and the adjusted R 2 (%). Values with absolute DIFT > (2× n p ) were considered influential points (p is the number of coefficients, and n is the number of data [47]) and were removed. Regressions were fitted using Statgraphics Centurion XVIII ® software (Statgraphics Technologies, Inc., The Plains, VA, USA). Regression models were assessed for seasonal and annual periods.  Table 2 shows the main statistics for the selected linear mixed model with the lowest AIC (482.30) and BIC (654.77). The continuous AR1 variance was included in this model. The "woodland type" (W, productivity) was shown to significantly affect soil CO 2 efflux (F = 7.8; p = 0.006; Table 2) because the mean soil respiration was higher in the young woodland (2.03 ± 0.10 vs. 1.73 ± 0.10 µmol CO 2 m −2 s −1 ). The interaction between the fixed factors also significantly influenced Rs (p = 0.028). Thus, throughout the study period, the soil respiration response differed in function of seasonality for each woodland type.

Seasonal Variation of Soil CO 2 Efflux in the Juniper Woodlands
The "date" factor (D) was shown to have a significant effect (F = 21.9; p < 0.001) on the response variable, Rs (Table 2); thus, the total soil respiration exhibited a clear seasonal variation throughout the study period ( Figure 4). Overall, the maximum soil respiration rate for the two woodlands occurred in the spring months (more pronounced in May 2018, under a higher soil temperature and water content; Figure 5), the most important period for vegetation activity, whereas the minimum values were registered in the winter months (January and February), the colder season ( Figure 4).  Table 2 shows the main statistics for the selected linear mixed model with the lowest AIC (482.30) and BIC (654.77). The continuous AR1 variance was included in this model. The "woodland type" (W, productivity) was shown to significantly affect soil CO2 efflux (F = 7.8; p = 0.006; Table 2) because the mean soil respiration was higher in the young woodland (2.03 ± 0.10 vs. 1.73 ± 0.10 µmol CO2 m −2 s −1 ). The interaction between the fixed factors also significantly influenced Rs (p = 0.028). Thus, throughout the study period, the soil respiration response differed in function of seasonality for each woodland type.

Seasonal Variation of Soil CO2 Efflux in the Juniper Woodlands
The "date" factor (D) was shown to have a significant effect (F = 21.9; p < 0.001) on the response variable, Rs (Table 2); thus, the total soil respiration exhibited a clear seasonal variation throughout the study period ( Figure 4). Overall, the maximum soil respiration rate for the two woodlands occurred in the spring months (more pronounced in May 2018, under a higher soil temperature and water content; Figure 5), the most important period for vegetation activity, whereas the minimum values were registered in the winter months (January and February), the colder season ( Figure 4). The soil CO2 efflux as a function of the woodland type was found to be clearly significant in summer months, the most important drought period. In this case, the CO2 efflux in soils was higher for the young woodland area (Figure 4). Additionally, in spring of 2017, there were significant differences between measures. In general, the soil respiration  significantly decreased from spring to autumn in both sites. In the cold seasons, winter and autumn, there were no significant differences between juniper woodland types.

Microclimatic Conditions and Organic Matter: Soil Temperature (Ts, °C), Soil Water Content (Sw, %), and Soil Organic C (%)
Seasonal changes in soil temperature (Ts, °C) and soil water content (Sw, %) are shown in Figure 5. We observed that the temporal variation of the soil temperature was similar for mature and young woodland areas. However, Figure 5 shows that the soil moisture was often lower in the mature woodland because soils have less depth and water retention capacity ( Table 1). The microclimatic variables were negatively correlated, in such a way that when the soil temperature increased, the soil moisture was progressively depleted.
The soil temperature ( Figure 5) ranged from −1.3 ± 0.8 °C and −0.4 ± 0.3 °C (mean ± standard error) for young and mature woodlands in winter, respectively, to 27.8 ± 1.7 °C and 24.3 ± 2.9 °C (summer months, young and mature woodlands, respectively). The seasonal behavior of the soil temperature indicated a positive correlation between this variable and the soil CO2 efflux: greater soil CO2 efflux was obtained with higher soil temperature values, i.e., in spring and summer.
The soil water content ( Figure 5) varied from 0.9% ± 0.1% and 1.7% ± 0.2% (27 July 2017, in summer, for mature and young woodland, respectively) to 27.9% ± 1.9% and 27.1% ± 1.9% (27 April and 10 May 2018, in spring, for mature and young woodland, respectively), with a peak of 27.6% ± 1.9% on 17 January for the young site. Due to the inherent variability that defines this semi-arid climate, mainly the soil water content in spring showed a great variation between years: in spring of 2018, the soil humidity was much higher than that recorded in the previous springs. Contrary to what happened with the soil temperature, the seasonal variation in soil moisture was not always positively correlated with the soil CO2 efflux: when the soil moisture was maximal, the soil efflux was minimal (e.g., winter), and the low value of moisture in summer did not correspond to the lower soil CO2 efflux rates. In summer, the peak of soil C could also contribute to increase the respiration ( Figure 6). The soil CO 2 efflux as a function of the woodland type was found to be clearly significant in summer months, the most important drought period. In this case, the CO 2 efflux in soils was higher for the young woodland area (Figure 4). Additionally, in spring of 2017, there were significant differences between measures. In general, the soil respiration significantly decreased from spring to autumn in both sites. In the cold seasons, winter and autumn, there were no significant differences between juniper woodland types.

Microclimatic Conditions and Organic Matter: Soil Temperature (Ts, • C), Soil Water Content (Sw, %), and Soil Organic C (%)
Seasonal changes in soil temperature (Ts, • C) and soil water content (Sw, %) are shown in Figure 5. We observed that the temporal variation of the soil temperature was similar for mature and young woodland areas. However, Figure 5 shows that the soil moisture was often lower in the mature woodland because soils have less depth and water retention capacity ( Table 1). The microclimatic variables were negatively correlated, in such a way that when the soil temperature increased, the soil moisture was progressively depleted.
The soil temperature ( Figure 5) ranged from −1.3 ± 0.8 • C and −0.4 ± 0.3 • C (mean ± standard error) for young and mature woodlands in winter, respectively, to 27.8 ± 1.7 • C and 24.3 ± 2.9 • C (summer months, young and mature woodlands, respectively). The seasonal behavior of the soil temperature indicated a positive correlation between this variable and the soil CO 2 efflux: greater soil CO 2 efflux was obtained with higher soil temperature values, i.e., in spring and summer.
The soil water content ( Figure 5) varied from 0.9% ± 0.1% and 1.7% ± 0.2% (27 July 2017, in summer, for mature and young woodland, respectively) to 27.9% ± 1.9% and 27.1% ± 1.9% (27 April and 10 May 2018, in spring, for mature and young woodland, respectively), with a peak of 27.6% ± 1.9% on 17 January for the young site. Due to the inherent variability that defines this semi-arid climate, mainly the soil water content in spring showed a great variation between years: in spring of 2018, the soil humidity was much higher than that recorded in the previous springs. Contrary to what happened with the soil temperature, the seasonal variation in soil moisture was not always positively correlated with the soil CO 2 efflux: when the soil moisture was maximal, the soil efflux was minimal (e.g., winter), and the low value of moisture in summer did not correspond to the lower soil CO 2 efflux rates. In summer, the peak of soil C could also contribute to increase the respiration (Figure 6). Finally, the measures of total soil organic C (SOC; Figure 6) were in the range cited by Lal [48] for forest soils. Overall, the value of SOC was slightly greater in mature than in young woodlands (average annual value of 3.27% vs. 2.58%, respectively). Despite the higher litterfall in the young site, the existence of a greater amount of organic C in the mature woodland area is noteworthy and could be a consequence of the greater cumulative deposition of organic matter over the years [35]. SOC also experienced seasonal variations, reflecting a general trend with the peak in the summer (more pronounced in the young juniper woodland) and decreases in autumn and winter. At the season level, significant differences in SOC between sites were obtained in the months of spring. In spring 2016 and 2017 (the soil respiration was only measured until 10 May 2018), soil organic C was higher in the older woodland, probably due to the greater accumulation of C that did not return to the atmosphere by respiration, and a fraction of the annual C input was stabilized in soils [49]. As a consequence, the soil organic C content was more stable in the mature forest across the seasons. Table 3 shows the results of adjusting the regression model of soil CO2 efflux (Rs) as a function of soil temperature (Ts, °C), soil water content (Sw, %), and the woodland type, "W" (bivariate, full model: Model 1). Parameters showing that the dummy variable (W) was significant (this represents two independent equations for each woodland) resulted only in fitted regressions for summer, spring 2017, and for all seasons or yearly (Table 3).   Finally, the measures of total soil organic C (SOC; Figure 6) were in the range cited by Lal [48] for forest soils. Overall, the value of SOC was slightly greater in mature than in young woodlands (average annual value of 3.27% vs. 2.58%, respectively). Despite the higher litterfall in the young site, the existence of a greater amount of organic C in the mature woodland area is noteworthy and could be a consequence of the greater cumulative deposition of organic matter over the years [35]. SOC also experienced seasonal variations, reflecting a general trend with the peak in the summer (more pronounced in the young juniper woodland) and decreases in autumn and winter. At the season level, significant differences in SOC between sites were obtained in the months of spring. In spring 2016 and 2017 (the soil respiration was only measured until 10 May 2018), soil organic C was higher in the older woodland, probably due to the greater accumulation of C that did not return to the atmosphere by respiration, and a fraction of the annual C input was stabilized in soils [49]. As a consequence, the soil organic C content was more stable in the mature forest across the seasons. Table 3 shows the results of adjusting the regression model of soil CO 2 efflux (Rs) as a function of soil temperature (Ts, • C), soil water content (Sw, %), and the woodland type, "W" (bivariate, full model: Model 1). Parameters showing that the dummy variable (W) was significant (this represents two independent equations for each woodland) resulted only in fitted regressions for summer, spring 2017, and for all seasons or yearly ( Table 3).

Seasonal Effects of Soil Temperature and Soil Water Content on Soil CO 2 Efflux
As a result of the above, Table 4 reflects the models fitted for each woodland type and as a function of seasons. This table also summarizes the microclimatic conditions (Ts and Sw) for the study period. No differences in soil temperature between juniper sites were found within a season (p > 0.05). Therefore, the differences in the respiratory process between woodlands should be attributed to other factors, such as those related to site or soil water content. Table 3. Significant parameters (±standard error), model significance (p), and goodness of fit (adjusted R 2 ) in the bivariate model (Model 1), fitted for seasonal and annual periods. The model explains the soil respiration (Rs, µmol CO 2 m −2 s −1 ) as a function of soil temperature (Ts, • C), soil water content (Sw, %), and woodland type (W): Rs = exp (β 0 +β 0 W)+(α+α W)×Ts Sw (β 1 +β 1 W) . n.s. = non-significant (p > 0.05, 95% probability). In the model, Ts is standardized by 10 • C (Ts/10).

Years
Seasons  Table 4. Fitted models for each woodland type and season period. The values of the temperature coefficient, Q 10 , were calculated when the exponential relationship between soil respiration (Rs) and soil temperature (Ts) was significant (Model 2). In this case, Q 10 was calculated by applying: Q 10 = e (α+α W) (α for mature woodland, and α + α for young woodland, if W was significant). *: Q 10 coefficient depends on soil water. Average values of Ts (soil temperature, • C) and Sw (soil water content, %) within each adjusting period followed by different letters reflect significant differences (Fisher's LSD test, 95% probability, α = 0.05). The full bivariate model (Model 1, soil temperature and soil water content as independent variables) was only significant in spring 2018 for both woodlands, and spring 2017 in the young site (Tables 3 and 4). In addition, Model 1 reached the highest goodness of fit in spring 2018 (R 2 = 94.1%; Table 3). In this model fitted for spring 2018, the high weight of the "soil water" variable (1.09) reflected that the exceptional CO 2 efflux rate was mainly due to high edaphic moisture levels.

Years
The soil water content seems to control the fitted models of soil respiration under two microclimatic conditions (Table 4): (i) under drought or low edaphic humidity (summer 2016 and spring 2017), and (ii) under high levels of soil moisture content (Sw > 25%) coinciding with warm-temperate temperatures (15-20 • C) in the growing season (conditions in spring of 2018). These microclimatic conditions also coincided with the season of the highest CO 2 efflux recorded in both woodlands (especially in the mature site; Figure 4).
The univariate model (Model 2, soil respiration as a function of soil temperature) was only significant in spring 2016 (with moderate soil water content levels) and yearly or throughout all seasons (Tables 3 and 4). Univariate models explained no more than 37% of the variance in CO 2 efflux (minor compared to when Sw was included in the models). Model 2 was the same for young and mature woodland areas in spring 2016, whereas when the model was fitted for all seasons (annual data in 2016, 2017, and all data pooled), two models resulted as a function of the woodland type. Thus, the model explained different annual rates of soil respiration for both types of juniper woodland as a function of only the soil temperature: for a certain increase in soil temperature (Ts), the variation in soil CO 2 efflux was significantly greater in the young juniper woodland, the most productive site (for example 1.40 vs. 1.20 in 2017). In these cases, the intra-annual effect of soil water content was confounded with the evolution of the soil temperature. Consequently, the global temperature coefficient, Q 10 (rate of respiration increases when the temperature increases by 10 • C), was significantly higher for the young woodland (1.40) than for the mature one (1.20). Regardless of the reference temperature (0 • C), Q 10 was invariant when the annual data were fitted by utilizing the exponential form.
Finally, for winter and autumn in the two woodland types, and in summer 2017 also for both woodlands (when soil moisture was minimal), the resulting models showed that soil respiration was independent of the soil temperature and soil water (Tables 3 and 4), thus representing a constant function or value (e β0 ). For this, the regression was not significant (p > 0.05; Table 3). The constant value indicated a higher level of respiration from summer to winter, and in the colder season, the estimated soil respiration was minimal (0.64 µmol CO 2 m −2 s −1 ). Models of soil respiration with low, constant values could indicate minimal activity in soil, due to low or high temperatures or a severe lack of humidity, conditions that would induce junipers and soil microbes to slow down their activity. However, the response of each woodland type to the mentioned edaphic conditions was different: under the microclimatic conditions of autumn and winter (no lack of moisture), no differences in soil respiration between sites were found, whereas in the driest season of the study period, summer 2017, the soil respiration was higher in the young site, which seems to indicate that this woodland could have a better response to most severe droughts.

Effects of Juniper Woodland Type and Seasonality on Soil CO 2 Efflux
Due to the strong relationship between woodland type (confounding effect of productivity and maturity) and seasonality, we discuss the effects of these factors together. Globally, the total soil CO 2 efflux measured at both juniper woodland sites fits the previously reported range (0.5-5 µmol CO 2 m −2 s −1 ) for forest ecosystems [50,51]. However, the annual mean soil CO 2 efflux appeared low (2.03 and 1.73 µmol CO 2 m −2 s −1 , for young and mature junipers, respectively) when comparing it with values from other biomes. For example, it was lower than the mean soil respiration for evergreen needleleaf forests (2.4 [50]). The results are an indicator that Spanish juniper woodlands have low productivity, which is in accordance with a previous study that focused on stem respiration in juniper trees [44].
We highlight that our results indicated a great level of seasonality in soil CO 2 efflux in both juniper woodlands. It is expected that the highest level of soil respiration will be obtained in spring, the season of the greatest growth activity [52], whereas the low soil respiration activity during cold seasons could result from minimal microbial respiration and maintenance respiration from the roots of juniper trees [53]. This semi-arid climate, with a warm and dry season alternating with a cold, wet season, seemed to be a decisive factor in driving the seasonal behavior of soil activity [13].
Further, the overall mean annual soil CO 2 efflux was higher in the young woodland. This is because the soil at this site maintained more respiratory activity in periods with a lack of soil humidity (summer, or dry spring); in the other seasons, there were no significant differences between sites, except for the driest spring (2017). Several factors can explain this difference in soil respiration in drought conditions, but the most important could be: (i) the higher productivity (growth rate) of the young woodland area, (ii) the higher quality of the soil (soil depth and water storage capacity) at this site, and (iii) the better adaptation or response to drought in younger juniper sites. In fact, the latter two factors are related to tree physiology and growth rates, which partly define forest production.
Referring to productivity, parameters that define the growth rate of a forest site are correlated with soil respiration [3]. The LAI is the main factor correlated with soil respiration in forest ecosystems, due to the higher accumulation of organic matter (litterfall) and Net Primary Production [8,50]. Soil C in Mediterranean ecosystems is one of the main components of the soil efflux because the C pool provides the substrate for heterotrophic organisms [54]. Spanish Juniper is a conifer species that deposits a large proportion of litter in the late spring, or early summer, to reduce transpiration when soil water content reaches a minimum [33]. Although we have only measured the soil C (the values of litterfall of Table 1 only describe the site productivity), the peak of soil organic C in summer (which was more pronounced in the young woodland) reflected an increase in deposition of litter. This could be due to the labile C released in the first stage of litter decomposition in soils [55]. Previous studies (for example, García Morote et al. [35]) also revealed the increase of the organic C pool in the season of abscission (late spring or early summer). Therefore, soil organic C could also contribute to greater respiration in young woodland areas in summer compared to mature ones. In relation to this physiological adaptation, it was found that β-glucosidase activity and basal respiration increased in summer in juniper woodlands due to the peak of organic matter accumulation [35]. The study presented by Hibbard et al. [50] confirmed that soil respiration is positively related to previous litterfall, and Chiang et al. [56] showed that soil respiration increased immediately after litterfall input to the soil.
In accordance with Jia et al. [57] and Allison et al. [58], maturity is another factor that influences soil respiration and microbial activity. Gustavo et al. [59] identified a negative correlation between stand maturity and total soil respiration, but whereas root respiration declined with age, heterotrophic respiration was correlated with soil C, and this factor slightly increased in mature forests [35]. Thus, in our study, differences in the soil CO 2 efflux under drought conditions (mainly in summer, or dry spring) may have resulted from increased root respiration in young trees, as cited previously [41]. Other studies also found that the contribution of root-rhizosphere respiration increased during the drought, due to an increasing supply of recent photosynthates to the belowground system [60].
Regardless of the soil microclimatic conditions, our results also highlight an important link between tree physiology and soil processes. Given the scarcity of water resources in the mature site (with the worst soil; lower water storage capacity, in mm), it is likely that older junipers adapt to episodes of water stress by reducing growth and autotrophic respiration. This response was also observed for the stem respiration of mature juniper trees [44]. On the contrary, in the site with better soil (young woodland), vegetal activity (and autotrophic respiration) could be continued in dry periods due to the greater reserves of water and nutrients. In addition, soils in the study area have a bedrock of fissured and karstified dolomite, which permits juniper trees to obtain water from the deeper soil profiles in the summer [35], because Juniperus thurifera is a species that inverts a high percentage of biomass in its roots (belowground biomass above 30% [34]). Research in areas with similar Mediterranean climates has confirmed that deep-rooted trees have greater water availability in dry seasons [61].

Effects of Soil Temperature and Soil Water on Soil CO 2 Efflux
The response of soil respiration to microclimatic conditions has not been consistent across all ranges because there are multiple confounding factors in different seasons [62,63]. In our research, the temperature sensitivity of respiration decreased under colder, wet conditions (winter and autumn), suggesting that seasonality is explained by the soil temperature without water limitations, in the temperate-growing season (spring), or when the effect of soil water was confounded (annual data). On the contrary, soil moisture was found to be significant when it is a limiting factor for soil. It is known that, for ecosystems without drought stress, soil temperature is the most important factor in determining the seasonal variation of soil respiration [41]. However, soil respiration, like other physiological processes of plants, usually responds to the most limiting factor [64], and soil moisture can become a limiting factor, especially as temperature and evaporative demand increase [14]. Especially in arid and semi-arid ecosystems, soil water can be the major factor that limits respiration sensitivity to temperature, particularly in summer [53]. In this sense, previous authors have shown that soil temperature is the main driver of soil respiration under high soil moisture conditions (above 15%), whereas under low levels, the soil water content becomes a better predictor [41,54]. In accordance with Reichstein et al. [65], soil respiration was shown to have a decreasing sensitivity to temperature in response to drought. In this semi-arid climate, high soil temperatures coincide with low moisture availability (except for spring 2018), thus confounding the effects of temperature and soil water on soil CO 2 efflux.
However, our study also reflects the joint effect of high soil temperatures with very high soil moisture content (>25%) in a growing season (spring 2018). Under these microclimatic conditions, the highest respiration rates in both sites were recorded. The inherent interannual variability of the Mediterranean climate also determined differences in soil respiration and their relationships with edaphics conditions, for the same season between years, and mainly in the spring. Previous studies have also found that the relationship between CO 2 efflux and microclimatic conditions depends on the magnitude of soil moisture, for example above 15% was the better predictor of the total soil respiration [54]. Some researches fixed the optimal soil moisture content for total soil respiration within a range of 15-20% [66]. Other studies showed that maximum soil respiration occurred when the soil moisture content reached the field capacity, for example 31% in beech [67]. In this sense, we hypothesize that when junipers have high levels of available soil water under warm temperatures, soil respiration could be acceleretad, as cited by Inclán et al. [54] in other mediterranean forests. These conditions would also be favorable to increase microbial activity [68,69].
Finally, the global Q 10 of total soil respiration calculated for the two juniper sites fell within the range of 1.0-5.0 [70] reported for European and North-American forest ecosystems. However, in connection with the soil efflux rates measured, the Q 10 value was low compared with those reported for coniferous forest [70]. The Q 10 quotient may vary depending on the season, growth rate, or acclimatization situation [10,71]. Based on our results (when data were analyzed throughout all seasons), the Q 10 coefficient in juniper woodlands also depends on the site. Thus, we conclude that soil respiration rates also vary within the same species as a function of ecological requirements (for example soil quality), i.e., productivity, as has been cited for respiration at the stem level [44,72].

Conclusions
Our results confirmed the hypothesis that, in semi-arid juniper woodlands, there is a high correlation between the total soil respiration and productivity (quality site). Seasonality was also found to be a decisive factor in the variations in soil respiration as functions of microclimatic conditions in each woodland area. Total soil respiration was positively correlated to soil water content under low moisture conditions, i.e., in dry spring, and summer (the season in which soil organic C peaks). Under these conditions, young woodland appears to be better adapted to drought conditions, due to the higher respiration rates in most dry seasons. Our study has also shown that the effect of soil temperature on soil respiration can depend on changes in soil moisture content: soil respiration was maximum when soil moisture content reached the highest values (>25%) in both sites in the growing season (spring).
In accordance with the soil efflux rates measured, the Q 10 coefficient was found to be dependent on the site and seasons, although an exponential response of soil respiration as a function of soil temperature is only possible if soil moisture is not limiting and in temperate seasons. Pooling the data yearly makes it possible to calculate a "mean" Q 10 coefficient to represent the mean rate of soil activity (it was significantly higher for the young woodland than for the mature one), but seasonal information is lost, and the effects of soil moisture and soil C can be minimized at the season level (and confounded yearly).
Therefore, research in semi-arid sites that attempts to model respiration or upscale data to the ecosystem level should analyze data by phenophase.
In conclusion, our study represents an advance in knowledge on the variation in soil respiration in classic juniper forests (mature and young woodlands) growing in a semi-arid climate (at their distribution limits). The uniqueness of this juniper species is reflected by low soil respiratory rates (especially in the older stands), typical of species linked to scarcely productive woodlands, relict species that must be protected by law.