Soil Carbon Dioxide Respiration in Switchgrass Fields : Assessing Annual , Seasonal and Daily Flux Patterns

Quantifications of annual soil respiration in switchgrass systems are limited to the growing season or coarse-scale temporal sampling. This study evaluates daily and seasonal soil CO2 respiration in switchgrass croplands. Hourly measurements during a 12-month period were taken for soil CO2 flux, soil temperature, and soil moisture. Although both soil temperature and moisture were positively correlated with soil CO2 flux rates, soil temperature was the primary driver of soil respiration. During winter, lower soil temperatures corresponded with significant decreases in average daily CO2 flux rates, however, CO2 pulses associated with precipitation events increased flux rates up to three times the seasonal daily average. Soil temperature influenced both daily and seasonal flux patterns where the highest flux rates, up to 31.0 kg CO2 ha−1 h−1, were observed during the warmest hours of the day (13:00 to 15:00) and during the warmest season (Summer). Summer and Spring emissions combined accounted for 80.1% of annual flux, indicating that exclusion of non-growing season time periods may result in an underestimation of total annual CO2 efflux. Our results indicate that inclusion of the non-growing season and a fine-resolution temporal sampling approach provides more accurate quantifications of total annual CO2 emissions in switchgrass croplands.


Introduction
On a global scale, release of carbon dioxide (CO 2 ) through soil respiration is estimated to contribute more than ten times the amount of carbon (C) to the atmosphere than anthropogenic fossil fuel emissions [1,2].Consequently, changes in land use and management practices that affect soil respiration rates can have a significant impact on atmospheric CO 2 concentrations and soil C pools [3].In the United States, agriculture alone occupies 51.3% (1.16 billion acres) of the total land area [4] and is among the largest net sources of atmospheric C resulting from land use change [5].In recent years, the cultivation of bioenergy crops has been embraced by greenhouse gas (GHG) reduction efforts to not only mitigate atmospheric CO 2 emissions through the reduction of fossil fuel use, but to also promote C sequestration in agricultural soils [6,7].Evidence suggests that bioenergy crops can promote soil C capture and storage as much as 5.4 Mg C ha −1 yr −1 [7], however, a better understanding of soil CO 2 emissions, the primary mechanism of soil C loss in agricultural landscapes, is needed to determine the net C sequestration potential of these systems.Parkin and Kaspar [8] reported that careful temperature-based correction algorithms are necessary to accurately estimate annual CO 2 emissions from short term measurements.They suggested future research to assess the errors associated with C loss determined from short-term CO 2 flux measurements made at infrequent intervals.
Bioenergy crops are often regarded as C neutral or C sink systems due to their ability to offset fossil fuel GHG emissions, fix and deposit atmospheric C into soils, and recapture C released during utilization in crop production in subsequent growing seasons.Quantification of actual offsets generated by these crops has progressed to include many components of atmospheric C release and capture including soil and plant studies as well as complexities of the cultivation process [9,10].These assessments can be used to identify the most efficient bioenergy crops for meeting C offset and sequestration goals from various perspectives of production.In terms of GHG offsets, switchgrass has emerged as a highly effective bioenergy resource, reducing GHG emissions from fossil fuels by 115% compared to 85% and 40% by reed canary grass and corn, respectively [9].Studies assessing soil processes also indicate that although soil respiration in switchgrass fields can account for 44% of total plant biomass C loss, the increase in soil C content in these systems results in a greater net C uptake by soils than has been observed for some other biofuels [11,12].Soil CO 2 respiration studies assessing switchgrass croplands in the southeastern United States are, however, often limited either to data collected at a coarse temporal scale (weekly or biweekly) or to the growing season, including only warmer months.Scientists often try to estimate a total amount of C loss using the field measured data.The limited data for estimating a total annual soil C loss may mislead the results.
Extensive research has demonstrated that soil respiration rates are positively correlated with temperature [1,2,13,14] and soil moisture [15][16][17][18][19][20], both of which vary diurnally, seasonally, and episodically.Fine-resolution temporal data, i.e., multiple measurements per day, is required to assess episodic events, such as precipitation, to capture pulses in soil CO 2 respiration in response to soil wetting [16,20], which would be missed using a weekly to bi-weekly sampling approach [21].With respect to seasonal variation, although cooler soil temperatures are associated with a decrease in soil CO 2 rates, soil respiration has been shown to continue even in snowy, winter conditions [22].Omission of soil flux data during non-growing seasonal periods could thus significantly affect the net annual CO 2 flux estimated for the system.These limitations in soil CO 2 flux studies in switchgrass croplands create the potential to inaccurately estimate annual loss of C through soil respiration.
Although some studies do report switchgrass soil CO 2 flux rates based on fine-scale and long-term data, these studies have been conducted in regions outside of the southeastern United States [11,23,24], which differ significantly in annual climate regimes that affect soil temperature and soil moisture patterns.Applying these rates to switchgrass systems located in the southeastern U.S. may also misrepresent soil C losses within the region.Switchgrass production in East Tennessee has increased significantly through efforts to cultivate biofuel crops to offset fossil fuel production and use in the region.Consequently, geographically targeted studies are needed to accurately develop C sequestration values for switchgrass production in East Tennessee.
This study investigated diurnal and seasonal variation in soil CO 2 flux for a switchgrass cropland in East Tennessee, and assessed the role of soil moisture and soil temperature as controls on soil CO 2 flux rates.First, we hypothesized that soil CO 2 flux rates are higher than those reported for the Pacific Northwest and Central Plains regions of the United States due to the warmer climate and annual precipitation patterns in E TN.Further, we hypothesized that soil temperature was a more significant control on soil CO 2 flux rates than soil moisture due to the relatively consistent annual precipitation regime coupled with significant thermal changes throughout the seasons.A fine-scale, i.e., hourly, temporal sampling approach was used for one full year to provide a better understanding of soil respiration responses to diurnal, episodic, and seasonal conditions.

Site Description
The study was conducted in an Alamo switchgrass plot at the University of Tennessee Plant Science Experiment Station (UT PSES) and at the University of Tennessee Holston Dairy Farm (UT HDF) both located in Knoxville, Tennessee (35 • 58 N, 83 • 56 W).Soil properties for both sites are listed in Table 1.Regional climatic conditions are humid subtropical with an average annual precipitation of 121.9 cm and weekly averages ranging from 3.0 cm to 5.5 cm.The mean annual temperature is 14 • C with the hottest month, July, ranging from 21 • C to 32 • C and the coldest month, January, ranging from −1 • C to 8 • C. The data were collected using two weather stations located in the UT PSES and HDF.The study was initiated at the UT PSES in well-established switchgrass plots planted in 1992.Crop management at this site, however, did not incorporate the use of nitrogen fertilizer which is outlined in University of Tennessee (UT) Extension recommendations and commonly practiced in the southeastern United States to increase crop yield.Flux measurements at the UT HDF site began in mid-November and continued at this location through to the completion of the experiment.UT HDF switchgrass plots were established in 2004 and managed according to UT Extension recommendations which included the application of 67 kg ha −1 y −1 of nitrogen annually.Thus, we ended up having a total of 14 month of data from the two study sites.
With the exception of the use of nitrogen fertilizer, site characteristics at both the UT PSES and UT HDF sites were comparable to reduce error resulting from differences in site environmental conditions.In addition to similar soil profiles and regional climatic conditions, both sites adhered to a similar cropping and harvesting regime in which grass was cut between 10 and 15 cm once per year in early November.Although the use of nitrogen fertilizer improves switchgrass crop production, this practice has been shown to have mixed results related to soil CO 2 flux [2,16,20].

Soil Properties
Chemical soil properties for both UT PSES and UT HDF are listed in Table 1.The laser diffraction method (LDM) was employed to confirm soil classification based on particle size and ratio.Soil type classification was consistent within the first 30 cm of the soil profile with UT PSES characterized by a loam and UT HDF by sandy clay.Chemical soil analyses included soil pH, cation-exchange capacity (CEC) and soil organic carbon (SOC) content.CEC was determined using the soil testing protocol developed by the Clemson soil testing with mineral extraction performed using a Perkin-Elmer 5300 Dual View Inductively Coupled Plasma spectrometer.SOC was determined using the dry combustion method described in Sparks et al. [25].Six replicate soil samples from the four depths at each location were collected.All soil chemistry analyses were performed at the University of Tennessee Soil, Plant and Pest Center in Nashville, TN.

Soil Respiration Measurements
Daily CO 2 flux measurements were taken once per hour from 00:00 to 23:00 throughout the duration of the experiment with the exception of the cool season (December to March) in which measurements were taken once every two hours due to solar power limitations during this season.Flux rates determined from hourly data were used to assess annual, seasonal, and daily soil flux patterns.Soil CO 2 flux was measured using the LICOR LI-8100 Automated Soil CO 2 Flux system comprised of an infrared gas analyzer control unit and multiplexer capable of supporting continuous flux measurements for up to eight LICOR 8100 long-term closed chambers.A total of six chambers were established at the UT PSES study site.Three chambers were positioned atop switchgrass crowns to represent covered areas (Cover) and the other three chambers were positioned in the interspaces between crowns to represent bare soil (Bare).
The LICOR LI-8100 long-term CO 2 flux chambers were installed to implement the closed chamber method for estimating soil respiration of CO 2 .In this method, the LI-8100 Analyzer control unit estimates soil respiration rates by measuring the increase or decrease (flux) of gaseous CO 2 concentrations in the chamber headspace over time.The observation time for all measurements was 90 s to minimize chamber CO 2 concentration changes during analysis and a 30 s dead band was programmed to allow for equilibration of the chamber pressure upon closure.All measurement protocols were programmed using a LICOR 8100 palm wireless controller linked with the LI-8100 Analyzer control unit.
Each flux chamber has an internal volume of 4076.1 cm 3 with an exposed soil area of 317.8 cm 2 .The whole system dimensions are 48.3 cm long by 38.1 cm wide by 33 cm high.To reduce error from lateral diffusion of CO 2 in the soil column, PVC soil collars measuring 20.3 cm in diameter were inserted to a depth of 3-5 cm and extended approximately 6-10 cm above the soil surface.A double gasket system sealed the chamber outside the soil collar and between the chamber and mounting plate to minimize CO 2 leaks and wind effects [26].Between measurements, the chamber head moved away from atop the measurement collar to minimize perturbations to the natural microclimate.Vegetation height was maintained at 5 cm in each chamber to allow for automated chamber movement, though we do recognize that clipping vegetation for maintenance may result in an underestimation of soil CO 2 flux rates due to decreased soil C accumulation [27].Soil temperature and moisture was also recorded at the time of each measurement.Soil temperature was monitored using Omega T-handled Type E thermocouples with 6.4 mm diameter and 250 mm immersion length inserted to a depth of 250 mm.Volumetric soil moisture was measured using 5 cm long EH2O Model EC-5 dielectric sensors.
A cubic polynomial regression analysis was applied to hourly measurements to evaluate diurnal patterns in soil CO 2 flux with changes in soil moisture and temperature [28] for both Cover and Bare plots separately.Seasonal and annual flux patterns were determined using daily average flux rates calculated from hourly measurements.Seasonal assessments were represented as Spring (20 March to 20 June), Summer (21 June to 22 September), Fall (23 September to 21 December), and Winter (22 December to 19 March).To determine annual flux patterns, averaged daily CO 2 flux values were summed with Cover and Bare sites representing 75% and 25%, respectively, of the total land area.

Annual CO 2 Flux
The total annual soil CO 2 emission for Cover and Bare sites were 81.17 Mg CO 2 ha −1 y −1 and 51.69 Mg CO 2 ha −1 y −1 , respectively.These values were used to determine an average weighted annual CO 2 emission of 73.80 Mg CO 2 ha −1 y −1 for the switchgrass plot, where Cover is representative of 75% and Bare 25% of the total land surface area.This equates to a loss of 20.14 Mg C ha −1 y −1 .The total annual soil CO 2 flux measured here for switchgrass is significantly greater than estimations of annual soil CO 2 efflux reported for other perennial grasses [29,30].Compared to other crop grasses, annual loss of CO 2 -C flux for switchgrass is five times greater than has been observed for barley, 4.0 Mg C ha −1 y −1 and more than two times greater than reports for crops such as wheatgrass, 8.65 Mg C ha −1 y −1 [30].In some instances, however, annual soil flux under annual and perennial crop grasses have been shown to range between 14.8 Mg CO 2 -C ha −1 y −1 and 20.3 Mg CO 2 -C ha −1 y −1 during the growing season alone [20].This may be due to differences in climate, environmental conditions and species types.The warm, humid climate in the East Tennessee region may provide conditions conducive to higher average soil respiration rates compared to other switchgrass croplands in other regions.This is supported by lower estimated annual flux rates reported from other regions such as North Dakota (10.8 Mg CO 2 -C ha −1 y −1 ) [11], Iowa (10.2 Mg CO 2 -C ha −1 y −1 ) [24], and Louisiana (3.72 Mg CO 2 -C ha −1 y −1 ) [31].Consequently, although it has been suggested that switchgrass croplands in the East Tennessee region function as C sinks or C neutral, our results suggest that loss of soil C through CO 2 respiration may exceed an annual SOC accumulation rate of 7.55 Mg C ha −1 y −1 [32].These findings demonstrate the need for a more detailed analysis of total net ecosystem exchange of C in East Tennessee switchgrass croplands to better understand the role of these systems in local and regional C budgets.

Seasonal CO 2 Flux
Seasonal flux trends for total cropland area show a positive correlation between average seasonal temperatures and soil CO 2 respiration rates.The highest seasonal CO 2 fluxes were observed during the warmest season (summer) and the lowest fluxes during the coolest (winter) (Figure 1).Average daily CO 2 respiration rates increased and decreased throughout the Spring and Fall seasons, respectively.Increasing ambient and soil temperatures in the Spring stimulate both plant growth and microbial activity previously dormant during the Winter.The creation of plant litter and movement of new roots through the soil supply fresh carbon and oxygen to soil microbes, resulting in a steady increase in soil CO 2 respiration during the Spring transition from Winter to Summer.Similar soil temperature ranges were observed during the Spring and Fall transition, however, higher soil flux rates were observed during the Spring.Lower flux rates in the Fall coincide with generation of organic material with matured plants and decreasing ambient and soil temperatures, which can slow respiration from microbial and autotrophic processes [1,14].Further, average daily Spring soil moisture was greater than in the Fall at 25.8% and 19.2%, respectively.Greater soil moisture content has been shown to correlate with higher soil C content [16,20] which may also contribute to the greater flux rates during the Spring.Summer alone accounted for approximately 47.4% of annual CO 2 emissions while Spring and Summer combined was representative of 80.1% of the annual flux.The Fall and Winter seasons contributed least to the overall annual soil CO 2 flux, 13.3% and 6.6% of respectively.However, representing close to 20% of total annual emissions, exclusion of these seasons may result in a significant underestimation of total annual soil CO 2 respiration.
Average respiration rates at Cover sites were greater than at Bare sites for every season (Figure 1), with the greatest flux rate offsets between site types observed during Spring and Summer.Previous studies have shown that vegetated soils can exhibit higher rates of soil respiration which can be correlated with larger microbial communities [33,34].Results of this study are consistent with these findings, with higher soil CO 2 respiration under Cover due to greater plant biomass, including switchgrass crowns and root structure, as well as increased microbial activity due to in situ litter supplied by the switchgrass crowns [33,34].Reduced flux variation between Cover and Bare during Winter was mostly likely the result of reduced microbial activity due to cooler soil temperatures, in spite of higher soil moisture content during this season.
has been shown to correlate with higher soil C content [16,20] which may also contribute to the greater flux rates during the Spring.Summer alone accounted for approximately 47.4% of annual CO2 emissions while Spring and Summer combined was representative of 80.1% of the annual flux.The Fall and Winter seasons contributed least to the overall annual soil CO2 flux, 13.3% and 6.6% of respectively.However, representing close to 20% of total annual emissions, exclusion of these seasons may result in a significant underestimation of total annual soil CO2 respiration.Average respiration rates at Cover sites were greater than at Bare sites for every season (Figure 1), with the greatest flux rate offsets between site types observed during Spring and Summer.Previous studies have shown that vegetated soils can exhibit higher rates of soil respiration which can be correlated with larger microbial communities [33,34].Results of this study are consistent with

Daily CO 2 Flux
Daily CO 2 flux patterns followed the patterns of diurnal soil temperature shifts.During all seasons except winter, the highest and lowest daily flux rates were observed during the warmest and coolest hours of the day, respectively (Figure 2).these findings, with higher soil CO2 respiration under Cover due to greater plant biomass, including switchgrass crowns and root structure, as well as increased microbial activity due to in situ litter supplied by the switchgrass crowns [33,34].Reduced flux variation between Cover and Bare during Winter was mostly likely the result of reduced microbial activity due to cooler soil temperatures, in spite of higher soil moisture content during this season.

Daily CO2 Flux
Daily CO2 flux patterns followed the patterns of diurnal soil temperature shifts.During all seasons except winter, the highest and lowest daily flux rates were observed during the warmest and coolest hours of the day, respectively (Figure 2).During Spring, Summer and Fall the highest average daily fluxes occurred between 13:00 and 16:00 while the lowest flux rates were measured between 00:00 and 09:00.The greatest diurnal variation in hourly flux rates was observed during the summer with an average daily range of 6.79 (±0.77) kg CO2 ha −1 h −1 .The highest average hourly flux rates for the total switchgrass cropland area were measured during June and ranged from 20.6 kg CO2 ha −1 h −1 at 05:00 and 31.0 kg CO2 ha −1 h −1 at 15:00.The smallest daily variation in hourly flux rates was observed during the winter with an average daily range of 1.39 (±0.21) kg CO2 ha −1 h −1 .Diurnal flux patterns were consistent with the patterns between soil temperature and soil CO2 respiration discussed with respect to seasonal trends.Higher temperatures during the daylight hours corresponded with daily peak soil respiration rates.Conversely, the decrease in CO2 flux throughout the evening correlated well with decreasing daily soil temperatures which slows biological processes.In addition to lower average soil temperatures, reduced diurnal variations in soil temperature during this season resulted in more consistent conditions from day to night and contributed to reduced variation in diurnal flux rate trends.Reduced daily variations in soil temperature under wetter soil conditions combined with cool season temperatures have been shown to be correlated with decreased soil respiration variability [35].The During Spring, Summer and Fall the highest average daily fluxes occurred between 13:00 and 16:00 while the lowest flux rates were measured between 00:00 and 09:00.The greatest diurnal variation in hourly flux rates was observed during the summer with an average daily range of 6.79 (±0.77) kg CO 2 ha −1 h −1 .The highest average hourly flux rates for the total switchgrass cropland area were measured during June and ranged from 20.6 kg CO 2 ha −1 h −1 at 05:00 and 31.0 kg CO 2 ha −1 h −1 at 15:00.The smallest daily variation in hourly flux rates was observed during the winter with an average daily range of 1.39 (±0.21) kg CO 2 ha −1 h −1 .Diurnal flux patterns were consistent with the patterns between soil temperature and soil CO 2 respiration discussed with respect to seasonal trends.Higher temperatures during the daylight hours corresponded with daily peak soil respiration rates.Conversely, the decrease in CO 2 flux throughout the evening correlated well with decreasing daily soil temperatures which slows biological processes.In addition to lower average soil temperatures, reduced diurnal variations in soil temperature during this season resulted in more consistent conditions from day to night and contributed to reduced variation in diurnal flux rate trends.Reduced daily variations in soil temperature under wetter soil conditions combined with cool season temperatures have been shown to be correlated with decreased soil respiration variability [35].The cool, wet winters of the East Tennessee climate are consistent with these conditions, as are the results of this study.A cubic polynomial regression model showed that soil temperature accounted for 75% and 69% (p < 0.0001) of soil CO2 flux variability for Cover and Bare sites, respectively (Figure 5).The positive relationship between higher flux rates and increases in soil temperature corresponded with previous studies in which increased rates of microbial and plant activity with warmer soil temperatures resulted in increased soil CO2 respiration [1,13,36].When applied to determine the impact of soil moisture on soil CO2 respiration, the regression model accounted for only 6% of Cover flux variability and 9% for Bare (p < 0.0001).Figure 3 shows the annual CO2 flux as a function of soil moisture.These results were inconsistent with some previous studies that demonstrated a more significant positive  A cubic polynomial regression model showed that soil temperature accounted for 75% and 69% (p < 0.0001) of soil CO2 flux variability for Cover and Bare sites, respectively (Figure 5).The positive relationship between higher flux rates and increases in soil temperature corresponded with previous studies in which increased rates of microbial and plant activity with warmer soil temperatures resulted in increased soil CO2 respiration [1,13,36].When applied to determine the impact of soil moisture on soil CO2 respiration, the regression model accounted for only 6% of Cover flux variability and 9% for Bare (p < 0.0001).Figure 3 shows the annual CO2 flux as a function of soil moisture.These results were inconsistent with some previous studies that demonstrated a more significant positive A cubic polynomial regression model showed that soil temperature accounted for 75% and 69% (p < 0.0001) of soil CO 2 flux variability for Cover and Bare sites, respectively (Figure 5).The positive Soil Syst.2018, 2, 13 8 of 11 relationship between higher flux rates and increases in soil temperature corresponded with previous studies in which increased rates of microbial and plant activity with warmer soil temperatures resulted in increased soil CO 2 respiration [1,13,36].When applied to determine the impact of soil moisture on soil CO 2 respiration, the regression model accounted for only 6% of Cover flux variability and 9% for Bare (p < 0.0001).Figure 3 shows the annual CO 2 flux as a function of soil moisture.These results were inconsistent with some previous studies that demonstrated a more significant positive correlation between soil moisture and soil CO 2 flux [16,17,19,20].Sainju et al. [20], for instance, showed that irrigation in some croplands increased CO 2 flux rates by up to a 50% due to a correlated increase in soil moisture.A notable difference in our study, however, is that the highest average soil moisture (~34%) was observed during Spring but CO 2 flux did not follow the soil moisture pattern.(Figure 4).Conversely, the warmer months were characterized by less frequent and less intense precipitation events resulting in a lower average soil moisture content, 21.5% during Summer, and more consistent soil moisture conditions.Figure 6 shows the annual precipitation and soil moisture.In all seasons, soil moisture increased immediately following precipitation events which simultaneously corresponded with episodic CO 2 pulse events (Figure 7).Although the average daily soil CO 2 flux during Winter was 60 kg CO 2 ha −1 , precipitation-induced pulses stimulated increases in daily flux rates of up to 178 kg CO 2 ha −1 .
Further, decreased CO 2 flux rates were observed at soil moisture contents of 5% or below.These results suggest that while soil moisture is less significantly correlated to soil CO 2 respiration than soil temperature in this system, episodic wetting and drying events do influence smaller scale flux rate variations.
Our results indicate that soil temperature is a primary driver of soil CO 2 flux rates in switchgrass croplands in East Tennessee.As the interactive effects of soil temperature and moisture influence soil biological activity [35], the lower variability in soil moisture during warmer months may have enhanced heterotophic and autotrophic respiration sensitivity to soil temperature variations during the peak soil CO 2 flux seasons (Spring and Summer).Still, the combination of both soil temperature and moisture increased the regression model's ability to explain flux variance for Cover and Bare sites to 82.7% and 81.5% (p < 0.0001), respectively.This suggests that while soil temperature plays a more significant role in soil respiration in this system, the interaction between soil temperature and moisture should be considered in flux rate models.

Conclusions
In this study, the majority (~93%) of the annual soil CO 2 efflux was observed during the Spring, Summer and Fall seasons.During the Winter, episodic precipitation events induced pulses of CO 2 respiration indicating that significant respiration continues throughout the entire year.The growing season for switchgrass excludes Winter and portions of Spring and Fall suggesting that studies limited to only this time period may underestimate annual CO 2 by more than 10%.Soil temperature was the primary driver of soil respiration rates during all seasons except Winter, during which soil moisture was more influential on average daily flux rates.CO 2 flux rates at Cover sites were consistently greater than flux rates observed at Bare sites; however, significant soil respiration at Bare sites demonstrates the importance of including bare ground interspaces in soil CO 2 flux assessments for switchgrass croplands.The total annual CO 2 efflux reported here is greater than has been reported for other perennial grasses and greater than CO 2 flux data reported for switchgrass in other regions in the United States.Greater annual CO 2 flux results from the wet, warm subtropical-temperate climate of the region.Our results suggest that geographically focused studies, which include the growing and non-growing seasons, are critical to estimations of annual soil CO 2 -C soil losses needed for the development of more accurate carbon offset values for switchgrass production as a bioenergy crop in Tennessee.

Figure 1 .
Figure 1.Average seasonal CO2 flux rates (kg CO2 ha −1 ) and soil temperature.Spring was dated 20 March to 20 June, Summer 21 June to 22 September, Fall 23 September to 21 December, and Winter 22 December to 19 March.

Figure 1 .
Figure 1.Average seasonal CO 2 flux rates (kg CO 2 ha −1 ) and soil temperature.Spring was dated 20 March to 20 June, Summer 21 June to 22 September, Fall 23 September to 21 December, and Winter 22 December to 19 March.

Figure 2 .
Figure 2. Daily variations in CO2 flux for each season.

Figure 2 .
Figure 2. Daily variations in CO 2 flux for each season.

3. 4 . 11 Figure 3 .
Figures 3 and 4 show the annual CO 2 flux along with soil temperature and moisture.Overall, CO 2 flux corresponds well with soil temperature change, while soil moisture was not a primary factor controlling CO 2 flux in our study.Soil Syst.2018, 2, x FOR PEER REVIEW 7 of 11

Figure 4 .
Figure 4. Variations of daily CO2 flux and soil moisture.

Figure 4 .
Figure 4. Variations of daily CO2 flux and soil moisture.

Figure 4 .
Figure 4. Variations of daily CO 2 flux and soil moisture.

Figure 5 .
Figure 5.A cubic polynomial regression analysis of daily CO2 flux relative to daily soil temperature.Data are representative of individual measurements collected throughout the entire study period.

Figure 6 .
Figure 6.Annual variations of precipitation and soil moisture.

Figure 5 .
Figure 5.A cubic polynomial regression analysis of daily CO 2 flux relative to daily soil temperature.Data are representative of individual measurements collected throughout the entire study period.

Figure 5 .
Figure 5.A cubic polynomial regression analysis of daily CO2 flux relative to daily soil temperature.Data are representative of individual measurements collected throughout the entire study period.

Figure 6 .
Figure 6.Annual variations of precipitation and soil moisture.

Figure 6 .
Figure 6.Annual variations of precipitation and soil moisture.

Figure 5 .
Figure 5.A cubic polynomial regression analysis of daily CO2 flux relative to daily soil temperature.Data are representative of individual measurements collected throughout the entire study period.

Figure 6 .
Figure 6.Annual variations of precipitation and soil moisture.

Figure 7 .
Figure 7. Sporadic CO2 flux response to soil moisture.Figure 7. Sporadic CO 2 flux response to soil moisture.

Figure 7 .
Figure 7. Sporadic CO2 flux response to soil moisture.Figure 7. Sporadic CO 2 flux response to soil moisture.

Table 1 .
Soil properties of the study site.