High Rainfall Inhibited Soil Respiration in an Asian Monsoon Forest in Taiwan

: Soil respiration represents the second largest carbon ﬂux, next to photosynthesis of the terrestrial biosphere, and thus plays a dual role in regional and global carbon cycles. However, soil respiration in Asian monsoon forests with high rainfall has rarely been studied. In this study, we continuously measured soil respiration using a 12-channel automated chamber system in a 61-year-old Japanese cedar forest in central Taiwan with annual rainfall greater than 2500 mm. A 4-year (2011–2014) continuous half-hourly dataset was used to quantify the inﬂuences of soil temperature and moisture, especially rainfall events, on both total soil respiration (Rs) and heterotrophic respiration (Rh). The annual mean Rs was approximately 10.8 t C ha − 1 (ranging from 10.7 to 10.9) t C ha − 1 , with Rh contributing approximately 74.6% (ranging from 71.7% to 80.2%). Large seasonal variations in both Rs and Rh were primarily controlled by soil temperature. Over 45.8% of total annual rainfall amounts were provided by strong rainfall events (over 50 mm), and over 40% of rainfall events occurred during summers between 2012 and 2014. These strong rainfall events caused rainwater to enter soil pores and cover the soil surface, which resulted in limited soil microorganism activity and, consequently, restricted CO 2 production. The mean Q 10 values were 2.38 (ranging from 1.77 to 2.65) and 2.02 (ranging from 1.71 to 2.34) for Rs and Rh, respectively. The Q 10 values in this study, which were lower than in global forest ecosystems, may imply that the interannual Rs values observed in this study that were caused by high rainfall were less temperature-dependent than the Rs levels in global forest ecosystems. Both Rs and Rh were negatively correlated with soil moisture, which indicated that the soil moisture levels in the studied forest were usually under saturated conditions. These results also provide the lack of data for respiration in the Asian monsoon region under high-rainfall conditions.


Introduction
Forest ecosystems significantly impact the global carbon cycle and serve as carbon sinks for the atmosphere [1]. The amount of biomass stored in forest ecosystems is estimated to be approximately 80% of the global terrestrial biomass [2]. Global forests contain approximately 1146 Pg of carbon, which represents 69% of the total carbon stored in forest soil and litter [3]. Greenhouse gas release rates of forest soil also affect global climate change. Soil carbon cycling processes have received much attention from ecological scientists and policy makers because of their possibility of influencing the soil carbon cycle through land use management. Therefore, accounting for forest carbon to clarify the response of carbon cycles to control mechanisms and to climate change is necessary to predict carbon magnitudes and identify possible mitigation and adaptation methods. Soil respiration (Rs) is mainly derived from autotrophic respiration (Rr) and heterotrophic respiration (Rh) in soil and is one of the largest carbon flux components in terrestrial ecosystems [4][5][6][7]. Of the total ecosystem respiration in temperate forests, 60-90% is contributed by Rs [8,9]. Rs is affected by biotic factors, such as root biomass and soil microorganism activity [10,11], and abiotic factors, such as soil temperature and soil moisture [12][13][14][15][16]. Rs generally has a positive relationship with soil temperature and can be described by a simple exponential function [12,[17][18][19][20][21] and by the exponential "Q 10 " value in temperate and tropical regions. Recently, however, more studies have focused on subtropical regions to fill the gaps in Rs data [19,[21][22][23][24][25][26]. On the other hand, Rs shows an inconsistent relationship with soil moisture. The inconsistent relationship between Rs and soil moisture is due to the effects of unusually low or high soil moisture levels, which suppress Rs production [15]. Liang et al. [12] found that Rs in a 45-year-old Japanese larch plantation in northern Japan had low correlation coefficients (R 2 < 0.05) with soil moisture. Adachi et al. [27] also demonstrated a weak relationship between soil moisture and Rs in primary and secondary tropical forests. However, Kosugi et al. [28] found that Rs had a strong linear relationship with soil moisture in a Southeast Asian tropical rainforest, while Adachi et al. [29] demonstrated a polynomial relationship between Rs and soil moisture in a dry tropical forest in Thailand.
Seasonal variations in Rs and their relationships with environmental factors, such as soil temperature and soil moisture, have been documented [15,26,30]. Although previous reports have explained the interannual variability in Rs as being due to spring soil temperatures or summer rainfall and soil moisture [15], the inconsistency of the results suggests that the cause of interannual Rs variability may depend on climate. The terrestrial carbon cycle is strongly influenced by climate [31]. Extreme climate events, such as droughts or storms, lead to disturbances in the carbon cycle of Southeast Asia, but scientific understanding of the impact on the carbon cycle is still obscure [31]. Extreme climate events intensify typhoon damage and affect Rs, Rh, and Rh/Rs. It is necessary to be able to forecast Rs and Rh and their relationships with climate change. However, Rs and Rh in Asian monsoon forests with high rainfall have rarely been reported. Taiwan is struck by typhoons at a rate of approximately three times per year, which also brings high rainfall [32]. Soil moisture is high during the typhoon season, especially in forest ecosystems that receive high rainfall amounts. Conducting field studies in subtropical forests is challenging because of typhoon damage and high rainfall, which make the quantification of soil carbon magnitudes difficult. Therefore, long-term Rs and Rh monitoring is necessary to clarify the dominant cause of variations and how this cause changes over years and seasons.
Cryptomeria japonica (L.) D. Don (Japanese cedar), an evergreen coniferous tree, is naturally distributed in Japan and China and is one of the most important and dominant plantation species in Japan and Taiwan. Japanese cedar plantations cover approximately 4.5 million ha in Japan and comprise approximately 45% of the total forest plantation area [33]. Japanese cedar was first introduced in Xitou, central Taiwan, in 1910 and was widely grown in the mountains of northern and central Taiwan between 1950 and 1980 and represented 10% of the total forest plantation area [34]. Investigating the intra-and interannual responses of Rs in these Japanese forests can help explain how forest ecosystems respond to climate change. Therefore, long-term Rs data are necessary to support such analyses. We installed a multichannel automated chamber Rs and Rh observation system in a Japanese cedar forest in Xitou in late 2011 and continuously measured Rs and Rh for 4 years. The objectives of this study were to (1) estimate Rs magnitudes at interannual timescales, (2) determine the relationships between Rs and environmental factors, and (3) quantify the influences of soil temperature and moisture, especially rainfall events, on both Rs and Rh.

Site Description
The Xitou Long-Term Ecosystem Research Site (XT, 23 • 39 N, 120 • 47 E) was established in 1950 in the Experimental Forest, National Taiwan University (EXFO, NTU), central Taiwan. In 2009, a flux tower was built at XT to observe carbon fluxes. This evergreen coniferous forest site is located at an altitude of 800-2000 m a.s.l. and the average slope is 15 • . The dominant species, C. japonica, was planted over 60 years ago. The stand density is 721 stems ha −1 , and the total basal area is 20.6 m 2 ha -1 in 2011. The soil is an Inceptisol, and the bedrock consists of lithic sandstone. Soil depths range from 40 to 50 cm. Soil pH is below 5.0 and contains low levels of base cations above the O horizon because of the high rainfall and subsequent leaching. Long-term weather data that were obtained from the Xitou agricultural weather station (270 m away from the studied plot) between 1941 and 2012 showed that the mean annual temperature (MAT) was 16.6 • C, with a mean maximum of 20.8 • C in July and mean minimum of 12 • C in January and the mean annual rainfall (MAP) was 2635 mm (Figure 1).

Site Description
The Xitou Long-Term Ecosystem Research Site (XT, 23°39′ N, 120°47′ E) was established in 1950 in the Experimental Forest, National Taiwan University (EXFO, NTU), central Taiwan. In 2009, a flux tower was built at XT to observe carbon fluxes. This evergreen coniferous forest site is located at an altitude of 800-2000 m a.s.l. and the average slope is 15°. The dominant species, C. japonica, was planted over 60 years ago. The stand density is 721 stems ha −1 , and the total basal area is 20.6 m 2 ha -1 in 2011. The soil is an Inceptisol, and the bedrock consists of lithic sandstone. Soil depths range from 40 to 50 cm. Soil pH is below 5.0 and contains low levels of base cations above the O horizon because of the high rainfall and subsequent leaching. Long-term weather data that were obtained from the Xitou agricultural weather station (270 m away from the studied plot) between 1941 and 2012 showed that the mean annual temperature (MAT) was 16.6 °C, with a mean maximum of 20.8 °C in July and mean minimum of 12 °C in January and the mean annual rainfall (MAP) was 2635 mm (Figure 1).

Automated Chamber Installation and Measurement
Rs is usually measured by the chamber-based method [12]. In this study, Rs was measured using a multichannel, automated chamber system [6,12]. Briefly, the main components consisted of a control unit (i.e., a control panel with infrared gas analyzer, IRGA, LI-840, LI-COR, Lincoln, NE, USA), data logger (CR1000, Campbell Scientific Inc., Logan, UT, USA), gas sampler and air compressor. There were 12 automated chambers (65 cm × 50 cm × 50 cm) that were constructed with aluminum bars and PVC plates as the frame and walls, respectively. The automated chambers were also equipped with two microfans for mixing air and with sampling tubes for controlling the chamber lids. The chamber lids are controlled by two pneumatic cylinders at a pressure of 0.2 MPa that are driven by a microcompressor (M-10, Hitachi Ltd., Tokyo, Japan). The chamber air was made to flow through the IRGA by a microdiaphragm pump (5 Lmin −1 ; CH-50, Enomoto Ltd., Tokyo, Japan), and the flow rate was 700 mL min −1 when the IRGA was measuring CO2 concentrations. During the measurements, the chamber air was mixed by two microfans when the lid was closed. Between measurements, the chamber lid was raised by air pressure. The chambers were closed sequentially, and Rs data were collected. The sampling period for each chamber was 150 s. Finally, datasets from each chamber were collected every halfhour. In September 2011, 12 chambers were installed randomly on the forest floor within a circular area of 20 m diameter near the flux tower, and all understory plants inside the chambers were cut every two weeks. From January 2012 to 2014, all chambers were divided into total Rs and heterotrophic groups, each with six chambers. Around the chambers used for measuring Rh, a 5 cm wide × 50 cm deep trench was dug along the chamber boundary using a sharp shovel to cut any roots. PVC plates (200 cm long × 50 cm wide × 0.4 cm thick) were then inserted into the trench to prevent root growth into the chambers.

Automated Chamber Installation and Measurement
Rs is usually measured by the chamber-based method [12]. In this study, Rs was measured using a multichannel, automated chamber system [6,12]. Briefly, the main components consisted of a control unit (i.e., a control panel with infrared gas analyzer, IRGA, LI-840, LI-COR, Lincoln, NE, USA), data logger (CR1000, Campbell Scientific Inc., Logan, UT, USA), gas sampler and air compressor. There were 12 automated chambers (65 cm × 50 cm × 50 cm) that were constructed with aluminum bars and PVC plates as the frame and walls, respectively. The automated chambers were also equipped with two microfans for mixing air and with sampling tubes for controlling the chamber lids. The chamber lids are controlled by two pneumatic cylinders at a pressure of 0.2 MPa that are driven by a microcompressor (M-10, Hitachi Ltd., Tokyo, Japan). The chamber air was made to flow through the IRGA by a microdiaphragm pump (5 Lmin −1 ; CH-50, Enomoto Ltd., Tokyo, Japan), and the flow rate was 700 mL min −1 when the IRGA was measuring CO 2 concentrations. During the measurements, the chamber air was mixed by two microfans when the lid was closed. Between measurements, the chamber lid was raised by air pressure. The chambers were closed sequentially, and Rs data were collected. The sampling period for each chamber was 150 s. Finally, datasets from each chamber were collected every half-hour. In September 2011, 12 chambers were installed randomly on the forest floor within a circular area of 20 m diameter near the flux tower, and all understory plants inside the chambers were cut every two weeks. From January 2012 to 2014, all chambers were divided into total Rs and heterotrophic groups, each with six chambers. Around the chambers used for measuring Rh, a 5 cm wide × 50 cm deep trench was dug along the chamber boundary using a sharp shovel to cut any roots. PVC plates (200 cm long × 50 cm wide × 0.4 cm thick) were then inserted into the trench to prevent root growth into the chambers.
Soil temperatures at 5-cm depth were monitored in each chamber using a thermocouple probe that was connected to a data logger. The volumetric soil water contents at 10-cm depth were measured around the chambers with a soil moisture sensor (SM300, Delta-T Device Ltd., Cambridge, UK). This soil moisture sensor was also connected to the data logger for data collection.

Data Analysis
Soil respiration rates (Rs, µmol CO 2 m −2 s −1 ) that were obtained from the chamber measurements (six chambers per treatment) were used to calculate mean soil respiration values, and the equation used is as follows where V is the chamber volume (m 3 ), P is the atmospheric pressure (Pa), R is the ideal gas constant (8.314 Pa m 3 K -1 mol -1 ), S is the soil surface area in the chamber (m 2 ), T is the air temperature in the chamber (K), and δC/δt is the slope of the CO 2 ratio (µmol mol −1 s −1 ). The chamber volume and soil surface area were 0.1625 m 3 and 0.325 m 2 , respectively. To analyze the influence of soil temperature on soil respiration in each chamber, the equation was transformed as follows (Lloyd and Taylor, 1994) where T soil is the soil temperature at a depth of 5 cm, a is the soil respiration at 0 • C, and b is the temperature sensitivity constant. The parameters of the different regression equations for the control and trench treatments were compared using analysis of covariance, and statistical significance was defined as p < 0.05. The value of b was also used to calculate the temperature sensitivity index (Q 10 coefficient), which is the relative increase in Rs with a 10 • C increase in soil temperature To analyze the response of Rs to soil moisture, we subtracted the simulated values of Rs by using Equation (2), and the residual values (RRs, the temperature-normalized Rs, µmol CO 2 m −2 s −1 ) were used for analysis. The relationship was fitted with the following concave-down regression equation where RRs (µmol CO 2 m −2 s −1 ) is the temperature-normalized Rs, θ is the volumetric soil moisture content (%), and c 1 , c 2 , and c 3 are curve-fitting parameters. The parameters of the different regression equations for the control and trench treatments were compared using analysis of covariance. Missing Rs and Rh data (gaps) in each chamber were estimated based on the annual temperature-response equation derived from Equation (2). The number of gap days totaled 102 during the 1487 days of measurements from September 2011 to December 2014, which corresponded to 6.8% of the entire measurement period. The largest gap was from 13 March to 7 April 2013 due to an IRGA problem, and this gap was the only gap longer than 30 days during the four years of observations. Finally, regression analysis and variance analysis of the continuous data were performed using Sigmaplot 14.0 software (Systat Software Inc., San Jose, CA, USA).

Environmental Factor Magnitudes on Interannual Timescales
Daily soil temperatures at 5-cm depths exhibited a typical seasonal pattern with a range of 8.9-20.3 • C and peaks in July and August ( Figure 2). The mean annual soil temperatures exhibited no differences from 2012 to 2014 (p > 0.05, Figure 3a). Daily soil moisture varied between 14.8% and 46.9%, with peaks generally occurring on days with high rainfall. Soil moisture gradually decreased from 40% to 20% because of a lack of rainfall from September to November and from January to April ( Figure 2). Mean annual soil moisture showed significant differences from 2012 to 2014 (p < 0.05, Figure 3b). From January 2011 to December 2014, the climatic conditions exhibited patterns that were typical for central Taiwan at the study site. The annual rainfall amounts were 3267 mm (24% higher than MAP), 2403 mm (9% lower than MAP) and 2062 mm (22% lower than MAP) from 2012 to 2014, respectively (Figures 2 and 3c). High-intensity rainfall events generally occurred from April to August (Figure 2). The characteristics of rainfall events in the study area are shown in Table 1. peratures exhibited no differences from 2012 to 2014 (p > 0.05, Figure 3a). Daily soil moisture varied between 14.8% and 46.9%, with peaks generally occurring on days with high rainfall. Soil moisture gradually decreased from 40% to 20% because of a lack of rainfall from September to November and from January to April (Figure 2). Mean annual soil moisture showed significant differences from 2012 to 2014 (p < 0.05, Figure 3b). From January 2011 to December 2014, the climatic conditions exhibited patterns that were typical for central Taiwan at the study site. The annual rainfall amounts were 3267 mm (24% higher than MAP), 2403 mm (9% lower than MAP) and 2062 mm (22% lower than MAP) from 2012 to 2014, respectively (Figures 2 and 3c). High-intensity rainfall events generally occurred from April to August (Figure 2). The characteristics of rainfall events in the study area are shown in Table 1.  peratures exhibited no differences from 2012 to 2014 (p > 0.05, Figure 3a). Daily soil moisture varied between 14.8% and 46.9%, with peaks generally occurring on days with high rainfall. Soil moisture gradually decreased from 40% to 20% because of a lack of rainfall from September to November and from January to April (Figure 2). Mean annual soil moisture showed significant differences from 2012 to 2014 (p < 0.05, Figure 3b). From January 2011 to December 2014, the climatic conditions exhibited patterns that were typical for central Taiwan at the study site. The annual rainfall amounts were 3267 mm (24% higher than MAP), 2403 mm (9% lower than MAP) and 2062 mm (22% lower than MAP) from 2012 to 2014, respectively (Figures 2 and 3c). High-intensity rainfall events generally occurred from April to August ( Figure 2). The characteristics of rainfall events in the study area are shown in Table 1.   Annual rainfall events with levels greater than 50 mm comprised 45.8% to 57.4% of strong rainfall events between 2012 and 2014 (Table 1). Rainfall event frequencies of 59%, 49%, and 40% were concentrated in the summer (June-August) from 2012 to 2014, respectively. Six-day periods with over 100 mm of daily rainfall occurred in April-August, and the highest recorded daily rainfalls were 770 mm in 2012, 269 mm in 2013, and 322 mm in 2014, due to typhoon movement. In contrast to the strong rainfall events, the total accumulation from small rainfall events (<10 mm) comprised 9.6% to 16.4% of annual rainfall between 2011 and 2014 (Table 1). Short-term drought occurred from January to March 2013, especially in March (less than 64% of the mean monthly rainfall). Minimum soil moisture occurred in late autumn and early spring when rainfall stopped.   Significant relationships between temperature-normalized soil respiration and volumetric soil moisture were observed over the 4 years (p < 0.001, Figure 5). The regression curves of the relationships with volumetric soil moisture at 10 cm depth of Rs exhibited significant differences to those of Rh in each year (p < 0.05, Figure 5b

Response of Rainfall Intensity to Rs and Rh Magnitudes in the Summer Season
Pulse increases in Rs occurred after light rainfall events, and these increases continued for 2 to 3 days (Figure 6a,b). The soil moisture also increased immediately after rainfall events and decreased when rainfall stopped until the next rainfall event (Figure 6b). However, pulse decreases in Rs occurred immediately after strong rainfall events for which the hourly rainfall amounts exceeded 12 mm (Figure 6c,d).

Response of Rainfall Intensity to Rs and Rh Magnitudes in the Summer Season
Pulse increases in Rs occurred after light rainfall events, and these increases continued for 2 to 3 days (Figure 6a,b). The soil moisture also increased immediately after rainfall events and decreased when rainfall stopped until the next rainfall event (Figure 6b). However, pulse decreases in Rs occurred immediately after strong rainfall events for which the hourly rainfall amounts exceeded 12 mm (Figure 6c,d).

Variations of Rs and Rh with Soil Temperature, Temperature Sensitivity, and Soil Moisture in Summer
Rs and Rh showed significant exponential relationships with soil temperature at 5 cm depths in summer for each year (p < 0.001, Figure 7). The regression curves of the relationships with soil temperature at 5 cm depth of Rs exhibited significant differences to those of Rh in summer in each year (p < 0.05, Figure 7a  The relationships between temperature-normalized soil respiration and volumetric soil moisture in summer over the 3-year period were calculated (p < 0.001, Figure 8). The regression curves of the relationships with volumetric soil moisture at 10 cm depth of Rs exhibited significant differences with those of Rh in each year (p < 0.05, Figure 8a

Environmental Factors, Rs, and Rh Magnitude in Summer
The mean summer value of soil temperature in 2014 was significantly higher than those in 2012 and 2013 (p < 0.05, Figure 9a). Furthermore, the mean summer value of volumetric soil moisture in 2012 was significantly higher than those in 2013 and 2014 (p < 0.05, Figure 9c). The sum of summer rainfall decreased with increasing years.

Environmental Factors, Rs, and Rh Magnitude in Summer
The mean summer value of soil temperature in 2014 was significantly higher than those in 2012 and 2013 (p < 0.05, Figure 9a). Furthermore, the mean summer value of volumetric soil moisture in 2012 was significantly higher than those in 2013 and 2014 (p < 0.05, Figure 9c). The sum of summer rainfall decreased with increasing years.

Annual Rs and Rh Magnitudes for the Forest Ecosystem Carbon Magnitudes
The annual mean Rs values were approximately 10.8 (ranging from 10.7 to 10.9) t C ha −1 , with Rh contributing approximately 74.6% (ranging from 71.7% to 80.2%) during the study period (2012-2014). The annual Rs values were within the range of other reported estimates from subtropical forests in Southwest China (12.5 t C ha −1 y −1 [26], South China (10.2 t C ha −1 y −1 [35], and Taiwan (8.93 t C ha −1 y −1 ) [36].

Effects of Environmental Factors (e.g., Soil Temperature, Soil Moisture, and Rainfall) on Interannual Changes in Rs and Rh
In general, interannual changes in Rs are predominantly controlled by soil temperature [14,15,[36][37][38]. However, high rainfall occurred from 2012-2013, which resulted in well-watered conditions in the Xitou area ( Figure 2). Therefore, a weak relationship between Rs and soil temperature was observed in summer due to the high rainfall frequency and well-watered conditions (Figures 2,7 and 9 and Table 1).
Soil moisture exhibited a concave-down relationship with temperature-normalized soil respiration ( Figure 5). The relationships of Rs and Rh with soil moisture were weaker than those with soil temperature in terms of interannual changes. Cook et al. [39] reported that microbial decomposition was promoted in the early summer rainy season and was inhibited in late summer drought conditions. A shortened dry spring occurred in 2013, with the lowest mean soil moisture in March (28%, versus an annual mean of 36%, Figure  2). The relatively lower soil moisture may have reduced microorganism activity and root respiration, especially in the initial growing season. Therefore, soil temperatures showed a weaker relationship with Rs and Rh in 2013 than in the other years ( Figure 4).

Annual Rs and Rh Magnitudes for the Forest Ecosystem Carbon Magnitudes
The annual mean Rs values were approximately 10.8 (ranging from 10.7 to 10.9) t C ha −1 , with Rh contributing approximately 74.6% (ranging from 71.7% to 80.2%) during the study period (2012-2014). The annual Rs values were within the range of other reported estimates from subtropical forests in Southwest China (12.5 t C ha −1 y −1 [26], South China (10.2 t C ha −1 y −1 [35], and Taiwan (8.93 t C ha −1 y −1 ) [36].

Effects of Environmental Factors (e.g., Soil Temperature, Soil Moisture, and Rainfall) on Interannual Changes in Rs and Rh
In general, interannual changes in Rs are predominantly controlled by soil temperature [14,15,[36][37][38]. However, high rainfall occurred from 2012-2013, which resulted in well-watered conditions in the Xitou area ( Figure 2). Therefore, a weak relationship between Rs and soil temperature was observed in summer due to the high rainfall frequency and well-watered conditions (Figures 2, 7 and 9 and Table 1).
Soil moisture exhibited a concave-down relationship with temperature-normalized soil respiration ( Figure 5). The relationships of Rs and Rh with soil moisture were weaker than those with soil temperature in terms of interannual changes. Cook et al. [39] reported that microbial decomposition was promoted in the early summer rainy season and was inhibited in late summer drought conditions. A shortened dry spring occurred in 2013, with the lowest mean soil moisture in March (28%, versus an annual mean of 36%, Figure 2). The relatively lower soil moisture may have reduced microorganism activity and root respiration, especially in the initial growing season. Therefore, soil temperatures showed a weaker relationship with Rs and Rh in 2013 than in the other years ( Figure 4).
Pulse increases in Rs after rainfall events have been demonstrated in incubation [40], temperate forests [5,12,[41][42][43][44], semiarid steppes [45], and subtropical forests [46]. Rainfall events stimulate the mineralization of soil organic carbon by microorganisms, which lead to rapid CO 2 release (Figure 6a,b). On the other hand, when a rainfall event occurred, rainwater filled the soil pores and forced CO 2 from the soil pores into the atmosphere, which resulted in increased soil respiration. However, high rainfall occurred and brought strong rainfall events due to the monsoons in this area during the research period (Table 1); therefore, high soil moisture is the threshold to restrict Rs and Rh (Figure 8c). The strong rainfall events caused rainwater to not only enter the soil pores, but also to cover soil surfaces, which resulted in limited oxygen and CO 2 diffusion. Therefore, oxygen deficiency in soil pores decreased soil microorganism activity and consequently restricted CO 2 production, as shown in Figure 6c,d. On the other hand, the mean summer values of both Rs and Rh significantly increased with increasing years, which resulted in decreased rainfall sums and soil moisture in summer ( Figure 9). Annual rainfall amounts and strong rainfall events tended to decrease (Table 1), which indicated a decrease in water cover of the soil surface and an increase in Rs and Rh releases from the soil. However, both Rs and Rh were negatively correlated with soil moisture, which indicated that soil moisture in the studied forest usually represented saturated conditions (Figure 8). The reason for this finding may be that the higher MAP in this study area when compared to those in previous subtropical forest studies (i.e., MAP 1243 mm, Yang et al. [47]; MAP 1,840 mm, Tan et al. [26]), and the subsequent high soil moisture, suppressed Rs and Rh.

Trench Treatment and Interannual Changes in Rh
Rh accounted for more than 70% of Rs (Figure 3d), and the contribution of Rr was less than 30% of Rs. The contribution of Rh to Rs in the study was in the upper range of those in previous reports [48]. The trench method is widely used to measure Rh [5]. The initial Rh was high due to the disturbances and decomposition of the dying roots following installation of the trench treatment (Figures 2 and 3d). Ohashi et al. [49] reported that the acceleration of microbial respiration offset the decline in root respiration in the first year in a 10-year-old Japanese cedar forest. The rapid decomposition of dead roots may have offset the decrease in root respiration at this study site. Makita et al. [30] indicated that Rh is sensitive to seasonal rainfall patterns because soil moisture strongly affects microbial activity. However, the studied forest is mature (64 years old); at this forest age, the biomass of fine roots tends to decrease, and the biomass of coarse roots tends to increase. Zhang and Wang [50] reported that the decomposition rates of fine and coarse roots were different among conifer species. Nevertheless, the fine roots rapidly decomposed in the trench treatment, and the coarse roots consequently decomposed and offset the decrease in root respiration. Diaz-Pines et al. [51] suggested that the influence of dead root decomposition continued for several years after the installation of a trench treatment. In our study, the interannual value of Rh gradually decreased from 2012 (6.6 t C ha −1 y −1 ) to 2014 (7.8 t C ha −1 y −1 ), which resulted in a decreased biomass of dead roots (Figure 9d).

Interannual Changes in Temperature Sensitivity
Seasonal patterns in temperature sensitivity may be related to the decomposition of soil C and microbial metabolic activity [30,52]. Lee et al. [53] found similar results: Rs was lower in spring and early summer and higher in autumn. This difference may be explained by seasonal changes in soil heat transport and CO 2 fluxes [30]. Rs contributions from root and microbial activity are reduced due to the low temperatures and soil moisture in deeper soil layers in spring. In contrast, Rs increased in response to the warming of deeper soil layers in late summer and autumn, which resulted in the higher Q 10 observed in summer (Figure 7). Soil heat increases and decreases from the top downward in spring and in autumn, respectively. Therefore, the variations in CO 2 production with soil depth during the growing season may affect heat-transport-based hysteresis [30]. In this study, the mean Q 10 values were 2.38 (ranging from 1.77 to 2.65) and 2.02 (ranging from 1.71 to 2.34) for Rs and Rh, respectively, which were similar to the mean values found for coniferous forests (1.69) [54], Japanese cedar plantations (1.91) [40,53] and subtropical China (1.82) [40] but were lower than the mean value estimated for global forest ecosystems (3.0) [55]. The Q 10 values observed in this study, which were lower than those in global forest ecosystems, may imply that the interannual Rs caused by high rainfall in this study is less temperaturedependent than the Rs values in global forest ecosystems.

Conclusions
The present study provides important information for understanding the dominant cause of variations in soil respiration and how the dominant cause changes among years. The Rs and Rh values in this Japanese cedar plantation showed high interannual variations. Based on the 4-year continuous monitoring results, the primary environmental factor that controlled interannual variations in both Rs and Rh was soil temperature. With frequent, high rainfall and minimal drought in this study area, soil moisture may be the limiting factor for annual Rs and Rh. However, the temperature sensitivity of interannual Rs and Rh was lower than the mean value estimated for global forest ecosystems and may have been affected by the high rainfall in the study area. We did find that a significant effect of soil moisture in the studied forest usually occurred under saturated conditions, and high rainfall events inhibited both Rs and Rh, especially in summer. In conclusion, our results suggest that the annual mean Q 10 values under high rainfall in this study were lower than those of global forest ecosystems. These results also alleviate the lack of data on respiration in the Asian monsoon region under high rainfall conditions.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.