An evaluation of empirical regression models for predicting temporal variations in soil respiration in a cool-temperate deciduous broad-leaved forest
An evaluation of empirical regression models for predicting temporal variations in soil respiration in a cool-temperate deciduous broad-leaved forest
Journal of Ecology and Environment. 2010. Jun, 33(2): 165-173
Copyright ©2010, The Ecological Society of Korea
This is an Open Access article distributed under the terms of theCreative Commons Attribution Non-Commercial License( permits unrestrictednon-commercial use, distribution, and reproduction in any medium,provided the original work is properly cited.
  • Received : February 17, 2011
  • Published : June 01, 2010
Export by style
Cited by
About the Authors
Na-Yeon, Lee
Present adress: National Park Research Institute, Korea National Park Service, Namwon 590-811, Korea

Soil respiration (R S ) is a critical component of the annual carbon balance of forests, but few studies thus far have attempted to evaluate empirical regression models in R S . The principal objectives of this study were to evaluate the relationship between R S rates and soil temperature (ST) and soil water content (SWC) in soil from a cool-temperate deciduous broad-leaved forest, and to evaluate empirical regression models for the prediction of R S using ST and SWC. We have been measuring R S , using an open-flow gas-exchange system with an infrared gas analyzer during the snowfree season from 1999 to 2001 at the Takayama Forest, Japan. To evaluate the empirical regression models used for the prediction of R S , we compared a simple exponential regression (flux = ae bt : Eq. [1]) and two polynomial multiple-regression models (flux = ae bt × (θν ? c) × (d ? θν) f : Eq. [2] and flux = ae bt × (1 ? (1 ? (θν/c)) 2 ): Eq. [3]) that included two variables (ST: t and SWC: θν) and that utilized hourly data for R S . In general, daily mean R S rates were positively well-correlated with ST, but no significant correlations were observed with any significant frequency between the ST and R S rates on periods of a day based on the hourly R S data. Eq. (2) has many more site-specific parameters than Eq. (3) and resulted in some significant underestimation. The empirical regression, Eq. (3) was best explained by temporal variations, as it provided a more unbiased fit to the data compared to Eq. (2). The Eq. (3) (ST × SWC function) also increased the predictive ability as compared to Eq. (1) (only ST exponential function), increasing the R 2 from 0.71 to 0.78.
Soil respiration (RS) is an important component of the annual carbon (C) balance of forests, as RS is the second most important C flux after gross primary production and is the largest respiratory C flux at the ecosystem scale. Interest in the rate-controlling factors of RS is growing because of the potential for climate change to increase the flux of CO 2 from ecosystems to the atmosphere (Raich and Potter 1995). Information regarding the relationship between the RS rate and the variables that drive it is also necessary for the development of the models of value for the assessment of climate-change effects and for the interpretation of the processes involved.
The R S rate is characterized by sizeable seasonal fluctuations owing to seasonal changes in root and microbial respiration. Although several biotic and abiotic factors influence root and microbial activity, the control exerted by soil temperature (ST), and in some cases soil water content (SWC), is generally dominant. In the absence of
water stress, variations in RS account for the majority of the seasonal and diurnal variations in RS rates. In areas in which water stress occurs frequently, RS rates may not be correlated with ST, but rather with its SWC (Rout and Gupta 1989).
Microbial communities and plant root systems are particularly sensitive to changes in ST (Killham 1994). Both specific R S rates and microbial and root biomass are affected positively by the evaluated ST, and in the majority of ecosystems a positive correlation between the R S rate and ST is observed (Killham 1994), resulting in sizeable seasonal fluctuations in the flux rates. There is, however, currently no consensus regarding the exact form of the relationships (Lloyd and Taylor, 1994). The R S rate has been modeled using linear (Witkamp 1966, Anderson 1973), power (Kucera and Kirkham 1971), and sigmoid (Schlentner and van Cleve 1985, Janssens et al. 1999, Matteucci et al. 2000) relationships with ST. Exponential relationships, most especially the Q 10 relationship, are more frequently utilized to predict the R S rates from ST (Raich and Potter 1995, Boone et al. 1998, Davidson et al. 1998).
However, ST is frequently correlated poorly with the RS rate, and SWC may be the best predictor of RS rates (Rout and Gupta 1989, Holt et al. 1990, Keith et al. 1997). SWC may negatively affect RS rates when it becomes either very high (poor aeration and reduced CO2 diffusivity) or very low (desiccation stress) (Janssens et al. 2003). Even in ecosystems characterized by abundant rainfall, the RS rate generally peaks in spring and after rainfall events (Matteucci et al. 2000, Lee et al. 2002). Field observations have shown a more variable relationship between the SWC and RS rates. The response of the RS rate to SWC has been successfully described using linear (Rout and Gupta 1989, Holt et al. 1990), exponential (Keith et al. 1997, Davidson et al. 1998), and polynomial functions (Bunnell et al. 1977, Schlentner and van Cleve 1985, Carlyle and Than 1988). However, Pinol et al. (1995) found SWC to be significant at one of their sites. Davidson et al. (1998) also suggested that SWC and ST controlled RS as independent or confounded factors in a temperate mixed hardwood forest. Thus, these empirical models are site-specific, and only a few attempts have been made thus far to apply them to other locations.
Moreover, Hanson et al. (1993) and Janssens et al. (2003) determined, when including a moisture response function, that the empirical model fit the data better and more parsimoniously explained the temporal variability. However, only a relatively few studies have thus far evaluated empirical regression models, including SWC, in RS. The principal objectives of this study were to evaluate the relationship between RS rates versus ST and SWC using hourly and daily data collected for 1999 and 2001 in a cool-temperate deciduous broad-leaved forest, and to determine the empirical regression models for the prediction of RS.
- Site description
The study site is located on the northeast slope of Mt. Norikura, Gifu Prefecture, in central Japan (36°08’ N, 137°26’ E, 1,430 m asl). The area is classified as cool-temperate and is located in a transitional zone between the Pacific side and the Japan Sea side of Japan. From 1980 to 2002, the annual mean air temperature was 6.1°C and the annual mean precipitation was 2,275 mm (Takayama Experimental Station, Gifu University; 1,342 m asl). The total rainfall varied from year to year, with a rainy season in June and a typhoon season in September. The site is covered with snow from December to April. The annual mean snowfall is approximately 600 cm, and the maximum snow depth is 100 to 180 cm (data from Takayama Research Station).
Approximately 50-year-old secondary, deciduous broad-leaved forest composed mainly of Quercus crispula Blume and Betula ermanii Cham. grow at the site. In areas with brown forest soils (Dystric Cambisols) (Food and Agriculture Organization 1999), the forest floor vegetation consisted of Sasa senanensis Rehd., a bamboo. The trees produce leaf buds at the beginning of June, and the leaves begin to senesce in early October. The leaf area index of the forest is approximately 3.5 between the end of June and the end of September (Yamamoto et al. 1999).
The study site has also been one of the long-term flux monitoring sites in the AsiaFlux network established in 1993, with the cooperation of the National Institute of Advanced Industrial Science and Technology and Gifu University. The site description in this study has been previously detailed by Lee et al. (2006).
- Methodology
The soil CO 2 flux was measured at four points in the study site from September 1999 to September 2001. The soil respiration rate was determined via an open dynamic chamber technique similar to that described by Lee et al. (2002). During the spring growing season and the summer, measurements were taken every 2 weeks. During
the autumn, measurements were made once per month, with just a few accidental interruptions (e.g. the electricity failed during a rainstorm). The open-flow IRGA method was utilized for the flux measurements. Ambient air (50 cm from the soil surface) was passed through the chambers at a rate of 0.5 L/min. The carbon dioxide concentrations of the air being pumped into and withdrawn from each chamber in turn (N = 4) were measured via IRGA. One measurement cycle of 25 minutes was repeated for 24?48 hours. Methods for the soil respiration measurements have been detailed by Lee et al. (2002). The measuring system setup was comprised of one reference line and four identical sample lines, an IRGA in absolute mode (Li-6252; Li-Cor Biosciences, Inc., Lincoln, NE, USA) and a data logger (Model NR250; Keyence, Tokyo, Japan). The reference and sample lines were pneumatically independent of one another. In this system, the ambient air entered the chambers at a flow rate of 1.2-1.6 L/min during the summer, and 0.8 L/min during the winter. One measurement cycle of 25 minutes was repeated for 24-48 hours. During the CO 2 analysis, air flow to the IRGA was maintained at a rate 0.5 L/min lower than the chamber flow-through rate, and water vapor was removed prior to analysis with a perma-pure drier (MP-110?72P, ZBJ02502?72P; Fuji Electric Co, Ltd., Tokyo, Japan). The data logger recorded the electric output from the IRGA. Four chambers were placed 4-5 m apart. Each chamber consisted of two parts. The lower body, which featured two ports for the inlet and outlet of air, was a 15-cm-high PVC cylinder with an internal diameter of 21 cm. When the bottom edge of the cylinder was pushed 4 cm into the soil, the bottom circumscribed an area of 346 cm 2 . The upper part, a 2.5-cm-high PVC lid, was placed on the top of the body immediately prior to the commencement of the measurements. Surface litter was contained within the chambers. However, great care was taken to avoid live or standing dead vegetation.
ST was measured in each chamber during the flux measurement at 5-10 cm below the top of the forest litter layer, which was approximately 4-8 cm below the surface of the mineral soil, depending on local variations in the thickness of the litter layer. Continuous ST measurements were conducted at two points near the flux measuring points at 30 minutes intervals (1, 10, 20, and 50 cm below the top of the litter layer) using copper-constantan thermocouples over the entirety of the study period. SWC was measured every 30 minutes by time-domain reflectometry (Model CS612; Campbell Scientific Inc., Logan, UT, USA). Two 15-cm long rod-balanced probes were placed horizontally at soil depths of 15 and 40 cm below the top of the litter layer, as described by Davidson et al. (1998).
It is likely that ST and/or SWC in the surface soil layers fluctuates more widely than in deeper soil layers. In a previous study conducted at in our site, Mo et al. (2005) reported that the temperature of the upper soil may better predict daily soil carbon efflux than the temperature of the deeper soil, as the upper 20 cm of the soil could contribute greatly to the RS rate, owing to the fact that it contained more than 90% of the roots (Lee et al. 2003) and a greater proportion of the microbial biomass (Mishima 2002). In the present study, therefore, ST at 10 cm depth and SWC at 15 cm depth were used to derive the empirical regression models of the RS rates.
Table 1 shows the list of ST and SWC response functions of the R S rate that have been applied in the empirical models. In order to evaluate the ST response function, we fitted an exponential relationship (Eq. [1]) using the data of hourly and daily mean R S rates. To evaluate the SWC dependence, we selected the two SWC response functions, linear (Eq. [4]) and convex-shape relationships (Eq. [5]). To evaluate the empirical regression models used to predict R S , we conducted comparisons with a simple exponential regression (Eq. [1]) and two polynomial multiple-regression models (Eq. [2] and [3]) that included two variables (ST and SWC) that utilized hourly
List of temperature and moisture response functions of soil respiration rate that have been applied in empirical models
Lager Image
List of temperature and moisture response functions of soil respiration rate that have been applied in empirical models
data for R S . Regressions using ST are denoted by t, and SWC is denoted by θν.
The significance for all statistical analysis was accepted at P = 0.05. All statistical analyses were conducted using StatView ver. 5.0 (SAS Institute Inc., Cary, NC, USA).
- Observed RS rates during 1999-2001
From September 1999 through September 2001, the R S rates increased from May to August and then decreased during autumn (Fig. 1). In 1999, the R S rate was 600-900 mg CO 2 m ?2 h ?1 in early September, peaked at 700-1,300 mg CO 2 m ?2 h ?1 in late September, and decreased to 30-300 mg CO 2 m ?2 h ?1 in November ( Fig. 1 a). In 2000, the R S rate was 20-200 mg CO 2 m ?2 h ?1 in late April, increased to 200-400 mg CO 2 m ?2 h ?1 in May, peaked at 400-900 mg CO 2 m ?2 h ?1 in June and August, and decreased to 200 mg CO 2 m ?2 h ?1 in November ( Fig. 1 b). In 2001, the R S rate was 180-230 mg CO 2 m ?2 h ?1 in May, peaked at 700-1,300 mg CO 2 m ?2 h ?1 in June, and decreased to 700-1,100 mg CO 2 m ?2 h ?1 in September ( Fig. 1 c). The R S rates tended to be lower in 2000 than other years, possibly because of the slightly lower soil temperatures in 2000 (Lee et al. 2003).
In the Japanese cedar (Cryptomeria japonica D. Don) forest in our study watershed, Lee et al. (2008) reported that the temporal variation in RS from soil and snow surfaces was measured continuously with an open dynamic chamber system. They investigated that the mean values of RS rate obtained by all six chambers ranged between 255 and 1,529 mg CO2 m?2 h?1, with lower values in mid-May and at the end of November and higher values between the beginning of July and the middle of October. In a previous study conducted at our site, Mo et al. (2005) also reported that the daily average RS rate was 245 mg CO2 m?2 h?1 in late April, increased to 455 mg CO2 m?2 h?1 in mid May, peaked at 642 mg CO2 m?2 h?1in late August, and decreased to 232 mg CO2 m?2 h?1 in November. Therefore, our values were comparable to the values reported by Mo et al. (2005) and Lee et al. (2008).
- ST versus RS
A temperature function used for the prediction of R S was derived on the basis of 3 years of the Takayama Forest study. According to the daily mean data, Fig. 2a shows the relationship between the R S rates (g C m ?2 day ?1 ) and ST was expressed as follows:
Lager Image
This regression accounted for 69% of the variation in the respiration measurements, primarily representing the seasonal changes in R S at this site (N = 66). According to the hourly data, Fig. 2 b shows the relationship between
Lager Image
Seasonal changes in hourly soil CO2 efflux (RS, mg CO2 m?2 h?1) in 1999 (a), in 2000 (b), and in 2001 (c).
Lager Image
Relationship between soil CO2 efflux (RS) and soil temperature daily-based data (a) and hourly-based data (b). The solid and dashed lines represent regression curves for daily and hourly data, respectively.
the R S rates (mg CO 2 m ?2 h ?1 ) and ST was expressed as follows:
Lager Image
Fig. 2 b shows that the relationship between RS rates and ST on hourly data was also highly significant (R2 = 0.71, N = 2,050). These results demonstrated that ST was the primary factor driving temporal variations in RS.
In general, the daily mean RS rates were positively well correlated with ST, but no significant correlations were noted to exist between the ST and RS rates on the periods of a day based on the hourly RS data. In our study site, no significant correlations were noted between the ST and hourly RS rates recorded on a rainy day (Lee et al. 2002). They indicated that post-rainfall increases in the RS rate represent approximately 16-21% of the annual soil carbon flux.
Fig. 3 shows that at ST > 15°C, there was considerable scatter in residuals both above and below the predicted RS line for the function of ST (using the Flux hourly (ST) equation). It appears likely that values above and below the predicted RS line affected the SWC and/or rainfall conditions. Generally, it is fairly clear that ST was the main factor when ST < 15ºC (Figs. 1 and 3). It has been reported that when the ST was not limited, the SWC affects RS. Murayama et al. (2003) also reported that the variation in ecosystem respiration was relatively great when ST > 15ºC; this was attributed to SWC and rainfall events occurring at our study site. They suggested that a single exponential or Q10 functional relationship of biological respiration to ST does not appear to be an adequate representation of biological respiration processes over the entire growing season, at least those occurring at this study site. Lloyd and Taylor (1994) also reported that the assumption of an exponential Q10 relation between RS and ST is invalid and systematically results in underestimated fluxes at low temperatures, and overestimated fluxes at high temperatures.
Long-term continuous monitoring of R S rates in situ is also urgently necessary to clarify the major sources of uncertainty in our ecosystem carbon balance. Opens and closes were automatic, and the measurements are based on an open-flow dynamic method (AOCC), providing continuous and semicontinuous measurements for investigating the temporal variations in R S (e.g., Suh et al. 2006, Lee et al. 2008). The major advantages of continuous chamber systems are the availability of a series of data gathered over a long period of time, and measurements conducted under relatively undisturbed condi
Lager Image
Differences (observed – predicted RS) in hourly soil CO2 efflux (RS). At soil temperature (ST) > 15°C, considerable scatter in residuals both above and below the predicted RS line for the function of ST (using the Flux hourly (ST) equation).
Differences (observed ? predicted R S ) in hourly soil CO 2 efflux (RS). At soil temperature (ST) > 15°C, considerable scatter in residuals both above and below the predicted R S line for the function of ST (using the Flux hourly (ST) equation)
Lager Image
Relationships between soil CO2 efflux (RS) and soil water content in daily-based data (a) and hourly-based data (b). The solid line represents the regression curve constructed from the whole dataset (with drought stress). The dashed line represents the linear function (without drought stress). The list of regressions (Eq. [4] and [5]) that were fitted to the data from daily means and hourly soil respiration rates and the obtained parameter values. Regressions are listed in Table 1.
- SWC versus RS
We evaluated the relationships between soil respiration and SWC. In this investigation, when SWC ranged between 30 and 50% (without drought stress), we applied a linear function (Eq. [4]). When SWC ranged between 20 and 50% (with drought stress), we applied a convex-shape function (Eq. [5]).
SWC evidenced significantly negative correlations to the daily mean and hourly RS without drought stress. A SWC linear function accounted for approximately 56% of the RS rate on the daily mean data (P < 0.001, N = 66). With regard to the hourly data, we noted a significant linear relationship (R2= 0.42, P < 0.0001, N = 1,952) with SWC. A convex-shaped function with optimal SWC on all data (with drought stress) explained approximately 32% of the daily mean data (P < 0.001, N = 69) and 16% of the hourly data (P < 0.001, N = 2,050). The relationship between SWC and RS was described using a linear function, rather than a second-order polynomial function for water stress conditions ( Fig. 4 ).
RS rates achieved maximum values at SWC levels of approximately positive or convex-type relationships between SWC and RS rates are rarely observed in actual field conditions, except during short-term drought or saturation conditions(Davidson et al. 1998, 2000) or hot-dry summer climate conditions (Xu and Qi 2001). A second-order polynomial function with optimal SWC explained an additional 16 to 32 % of the temporal variability in RS that could not be accounted for by ST alone. Our results indicated that SWC generates considerable scatter in residuals both above and below the predicted RS line for the function of ST. Kang et al. (2003) indicated that spatial variability in RS is better explained with SWC than ST. Moderate drought may also partially explain the inter-annual variability in the RS rate (Epron et al. 1999). Even in forests with less pronounced drought stress, in which drought occurs either infrequently or moderately, the inclusion of moisture regressions in empirical models will increase the degree to which the models fit the data.
- ST × SWC versus RS
To estimate ST × SWC versus the RS model, we used the functions of Eq. (2) and (3) ( Table 1 ). This model showed that the RS rates increased with ST and were suppressed under dry or wet conditions.
Firstly, we utilized a multiple polynomial regression
Lager Image
Relationship between the observed soil CO2 efflux (RS) and predicted soil CO2 efflux (RS) hourly data in temperature dependent on the exponential function (Eq. [1], a) and two polynomial multiple-regression models (Eq. [2], b and Eq. [3], c). The solid and dash lines represent the linear function and confidence limit (95%).
(Eq. [2]) of our previous study (Lee et al. 2002) using the hourly data recorded on sunny days (without rainfall effect) ( Table 1 ). The data on rainy days and those obtained within 24 hours after the last rainfall were excluded from this calculation, thus eliminating the effects of rainfall on the R S .
Lager Image
The model explained approximately 96% of the variance in the daily soil carbon fluxes calculated from ST and SWC on sunny days during August 1999 to November 2000 (Lee et al. 2002). However, no significant relationship was detected between the observed and predicted RS rates using all datasets. In order to validate the regression model, we used all of the data collected between September 1999 and September 2001. There was considerable scatter in residuals both above and below the predicted RS line of the ST × SWC model (Y = 0.35X + 145, R2 = 0.34) ( Fig. 5 b).
Second, we applied the other multiple polynomial function (with rainfall effect, period April-November 2000 and May-November 2001) (Lee et al. 2005) using the hourly all-dataset according to equation (3) ( Table 1 ).
Lager Image
The ST × SWC model derived for the Eq. (3) function was also a good predictor (Y = 0.78X + 228, R 2 = 0.78), although not as good as a model derived specifically for the function Eq. (2) (Fig. 5c). The addition of a SWC function to the temperature-only model significantly increased the predictive power (Y = 0.83X + 89, R 2 = 0.71 for Eq. (1), Fig. 5a) (R 2 = 0.78 for Eq. (3), Fig. 5 c).
Similarly, Savage and Davidson (2001) demonstrated that an empirical model utilizing both ST and SWC had good explanatory power in the Harvard forest (R2 = 0.64 for temperature-only model and R2 = 0.70 for ST × SWC). Hanson et al. (1993) also demonstrated that the inclusion of SWC increased the predictive power of their model at four of their five sampling sites (R2 from 0.49-0.60 increased to 0.67-0.73).
In the current study, Eq. (3) was best explained by temporal variations, as it provided a more unbiased fit to the data as compared to Eq. (2). The relatively modest increases in R2 values when the SWC functions are added to the temperature-only regression models. The Eq. (2) has many more site-specific parameters than Eq. (3) and was significantly underestimated. The shape of the moisture response curve and the moisture content at which maximum respiration occurs depends on an array of site-specific factors such as soil texture and structure, amount and type of organic matter, and ST (Howard and Howard 1993).
Our values of the coefficients of a and b in Eq. (1) were 59.472 (mg CO2 m?2 h?1) and 0.163. In the Japanese cedar forest, however, the values were 360 (mg CO2 m?2 h?1) and 0.058 (Lee et al. 2008). Similarly, the parametrization of temperature functions is likely to be site-specific and to be confounded by other effects, such as interactions with water content. Temperature functions predicted seasonal variation in RS fairly well, but variations in precipitation and SWC were identified as the key to understanding inter-annual variations. We will need to analyze the data of two consecutive wet years to determine whether or not high respiration rates persist during both years or if the high respiration rates only occur during the first wet year after a dry year. Our on-going measurements may also afford us other opportunities to investigate a broad range of interannual variability.
The author would like to thank Prof. H. Koizumi of Waseda University for his encouragement and Dr. N. Saigusa of the National Institute of Environmental Studies for providing the long-term micrometeorological measurement data. The author also wishes to thank Prof. K. Nakane and Prof. T. Nakatsubo of Hiroshima University for their advice and assistance.
Anderson J.M 1973 Carbon dioxide evolution from two temperate deciduous woodland soils J Appl Ecol 10 361 - 378    DOI : 10.2307/2402287
Boone R.D , Nadelhoffer K.J , Canary J.D , Kaye J.P 1998 Roots exert a strong influence on the temperature sensitivity of soil respiration Nature 396 570 - 572    DOI : 10.1038/25119
Bunnell F.L , Tait D.E.N , Flanagan P.W , Van Cleve K 1977 Microbial respiration and substrate weight loss: I A general model of the influences of abiotic variables Soil Biol Biochem 9 33 - 40    DOI : 10.1016/0038-0717(77)90058-X
Carlyle J.C , Than U.B 1988 Abiotic controls of soil respiration beneath an eighteen-year-old Pinus radiata stand in south-eastern Australia J Ecol 76 654 - 662    DOI : 10.2307/2260565
Davidson E.A , Belk E , Boone R.D 1998 Soil water content and temperature as independent or confounded factors controlling soil respiration in a temperate mixed hardwood forest Global Change Biol 4 217 - 227    DOI : 10.1046/j.1365-2486.1998128.x
Davidson E.A , Verchot L.V , Cattanio J.H , Ackerman I.L , Carvalho J.E.M 2000 Effects of soil water content on soil respiration in forests and cattle pastures of eastern Amazonia Biogeochemistry 48 53 - 69    DOI : 10.1023/A:1006204113917
Epron D , Farque L , Lucot E , Badot P.M 1999 Soil CO2efflux in a beech forest: dependence on soil temperature and soil water content Ann For Sci 56 221 - 226    DOI : 10.1051/forest:19990304
Food and Agriculture Organization 1999 FAO STAT Statistical Database 1997 Food and Agriculture Organization Rome
Hanson P.J , Wullschleger S.D , Bohlman S.A , Todd D.E 1993 Seasonal and topographic patterns of forest floor CO2efflux from an upland oak forest Tree Physiol 13 1 - 15
Holt J.A , Hodgen M.J , Lamb D 1990 Soil respiration in the seasonally dry tropics near Townsville North Queensland Aust J Soil Res 28 737 - 745    DOI : 10.1071/SR9900737
Howard D.M , Howard P.J.A 1993 Relationships between CO2evolution moisture content and temperature for a range of soil types Soil Biol Biochem 25 1537 - 1546    DOI : 10.1016/0038-0717(93)90008-Y
Janssens I.A , Sampson D.A , Cermak J , Meiresonne L , Riguzzi F , Overloop S , Ceulemans R 1999 Above- and below-ground phytomass and carbon storage in a Belgian Scots pine stand Ann For Sci 56 81 - 90    DOI : 10.1051/forest:19990201
Janssens I.A , Dore S , Epron D , Lankreijer H , Buchmann N , Longdoz B , Brossaud J , Montagnani L , Valentini R 2003 Climatic influences on seasonal and spatial differences in soil CO2efflux In: Fluxes of Carbon Water and Energy of European Forests Springer-Verlag New York 235 - 256
Kang S.Y , Doh S , Lee D , Lee D , Jin V.L , Kimball J.S 2003 Topographic and climatic controls on soil respiration in six temperate mixed-hardwood forest slopes Korea Global Change Biol 9 1427 - 1437    DOI : 10.1046/j.1365-2486.2003668.x
Keith H , Jacobsen K.L , Raison R.J 1997 Effects of soil phosphorus availability temperature and moisture on soil respiration in Eucalyptus pauciflora forest Plant Soil 190 127 - 141    DOI : 10.1023/A:1004279300622
Killham K 1994 Soil Ecology Cambridge University Press Cambridge
Kucera C.L , Kirkham D.R 1971 Soil respiration studies in tallgrass prairie in Missouri Ecology 52 912 - 915    DOI : 10.2307/1936043
Lee M.S , Lee J.S , Koizumi H 2008 Temporal variation in CO2efflux from soil and snow surfaces in a Japanese cedar (Cryptomeria japonica) plantation Ecol Res 23 777 - 785    DOI : 10.1007/s11284-007-0439-z
Lee M.S , Mo W.H , Koizumi H 2006 Soil respiration of forest ecosystems in Japan and global implications Ecol Res 21 828 - 839    DOI : 10.1007/s11284-006-0038-4
Lee M.S , Nakane K , Nakatsubo T , Koizumi H 2003 Seasonal changes in the contribution of root respiration to total soil respiration in a cool-temperate deciduous forest Plant Soil 255 311 - 318    DOI : 10.1023/A:1026192607512
Lee M.S , Nakane K , Nakatsubo T , Koizumi H 2005 The importance of root respiration in annual soil carbon fluxes in a cool-temperate deciduous forest Agric For Meteorol 134 95 - 101    DOI : 10.1016/j.agrformet.2005.08.011
Lee M.S , Nakane K , Nakatsubo T , Mo W.H , Koizumi H 2002 Effects of rainfall events on soil CO2flux in a cool temperate deciduous broad-leaved forest Ecol Res 17 401 - 409    DOI : 10.1046/j.1440-1703.2002498.x
Lloyd J , Taylor J.A 1994 Funct Ecol 8 315 - 323    DOI : 10.2307/2389824
Matteucci G , Dore S , Stivanello S , Rebmann C , Buchmann N , Schulze E.D 2000 Soil respiration in beech and spruce forests in Europe: trends controlling factors annual budgets and implications for the ecosystem carbon balance In: Carbon and Nitrogen Cycling in European Forest Ecosystems: Ecological Studies Vol 142 Springer-Verlag Berlin
Mishima S.I 2002 Microbial biomass microbial respiration activity and CO2flux derived by microorganisms in cool temperate secondary forest floor in Japan Proceedingsof the VIII INTECOL International Congress of Ecology Seoul 183 - 184
Mo W , Lee M.S , Uchida M , Inatomi M , Saigusa N , Mariko S , Koizumi H 2005 Seasonal and annual variations in soil respiration in a cool-temperate deciduous broad-leaved forest in Japan Agric For Meteorol 134 81 - 94    DOI : 10.1016/j.agrformet.2005.08.015
Murayama S , Saigusa N , Chan D , Yamamoto S , Kondo H , Eguchi Y 2003 Temporal variations of atmospheric CO2concentration in a temperate deciduous forest in central Japan Tellus B 55 232 - 243    DOI : 10.1034/j.1600-0889.2003061.x
Pinol J , Alcaniz J.M , Roda F 1995 Carbon dioxide efflux and pCO2in soils of three Quercus ilex montane forests Biogeochemistry 30 191 - 215    DOI : 10.1007/BF02186413
Raich J.W , Potter C.S 1995 Global patterns of carbon dioxide emissions from soils Global Biogeochem Cycles 9 23 - 36    DOI : 10.1029/94GB02723
Rout S.K , Gupta S.R 1989 Soil respiration in relation to abiotic factors forest floor litter root biomass and litter quality in forest ecosystems of Siwaliks in northern India Acta Oecol Oecol Plant 10 229 - 244
Savage K.E , Davidson E.A 2001 Interannual variation of soil respiration in two New England forests Global Biogeochem Cycles 15 337 - 350    DOI : 10.1029/1999GB001248
Schlentner R.E , van Cleve K 1985 Relationship between CO2evolution from soil substrate temperature and substrate moisture in four mature forest types in interior Alaska Can J For Res 15 97 - 106    DOI : 10.1139/x85-018
Suh S.U , Chun Y.M , Chae N.Y , Kim J , Lim J.H , Yokozawa M , Lee M.S , Lee J.S 2006 A chamber system with automatic opening and closing for continuously measuring soil respiration based on an open-flow dynamic method Ecol Res 21 405 - 414    DOI : 10.1007/s11284-005-0137-7
Witkamp M 1966 Decomposition of leaf litter in relation to environment microflora and microbial respiration Ecology 47 194 - 201    DOI : 10.2307/1933765
Xu M , Qi Y 2001 Soil surface CO2efflux and its spatial and temporal variations in a young ponderosa pine plantation in northern California Global Change Biol 7 667 - 677    DOI : 10.1046/j.1354-1013.2001435.x
Yamamoto S , Murayama S , Saigusa N , Kondo H 1999 Seasonal and inter-annual variation of CO2flux between a temperate forest and the atmosphere in Japan Tellus B 51 402 - 413    DOI : 10.1034/j.1600-0889.1999020.x