Regeneration Traits of Four Dominant Species in a Cool-Temperate Conifer-Hardwood Mixed Forest , Northern Japan

Regeneration traits of four dominant species were studied during a decade in a cool-temperate conifer-hardwood mixed forest, northern Japan. Dominant species were three canopy species (Tilia japonica, Acer mono and Abies sachalinensis) and a subcanopy species Prunus ssiori. Regeneration traits differed among the four dominant species. The regeneration of a conifer Abies sachalinensis largely depends on major disturbances because its size structure was a bell-shaped pattern. The growth rate of the sub-canopy species Prunus ssiori increased faster with tree size than the other three species, and therefore, this species can reach reproductive stage faster. Although the number of recruits of Tilia japonica was less than the other three dominant species, about one-fourth of Tilia japonica regenerated by the sprouts. Thus, vegetative reproduction is important for the regeneration of Tilia japonica. On the contrary, many recruits were observed in Acer mono, and its recruits and saplings concentrated in canopy gaps, suggesting that the regeneration of Acer mono restricted to canopy gaps. Thus, clear differences in the regeneration traits were recognized in the four dominant species. This study suggests that the species-specific regeneration traits contribute to the species coexistence of the four dominant species through different regeneration niches.


Introduction
Conifer-hardwood mixed forests widely distribute in the middle latitude of How to cite this paper: Takahashi, K. (2017) Regeneration Traits of Four Dominant Species in a Cool-Temperate Conifer-Hardwood Mixed Forest, Northern Northern Hemisphere [1] [2] [3].Forest structure and dynamics of this type of forest ecosystem have often been examined in relation to major disturbances [4] [5].For example, several researchers have shown discontinuous age and size distributions, and discussed the importance of major disturbances for species coexistence [6] [7] [8] [9].However, all species in a conifer-hardwood mixed forest do not always depend on major disturbances for their regeneration.Return times for major disturbances are unpredictable, and sometimes exceed more than a millennium [10] [11] [12].Therefore, coexistence mechanism in conifer-hardwood mixed forests cannot be adequately explained by major disturbances alone.
In equilibrium conditions without major disturbances, tree competition is important role in structuring force of community through interspecific differences in density-dependent demographic rates (recruitment, growth and mortality), i.e., higher growth and recruitment rates and lower mortality in less crowded conditions, and vice versa in more crowded conditions [13] [14] [15] [16] [17].Kobe et al. [18] has examined the interspecific variation in sapling mortality in relation to the growth rate, and showed a positive correlation between growth rate in high light conditions and mortality in low light conditions.Thus, the potential growth rate of shade-intolerant species is high, but their regeneration is often restricted to canopy gaps [19] [20].Therefore, tree competition (or canopy cover) largely influences the regeneration dynamics through the density-dependent demographic rates.
A conifer Abies sachalinensis Masters and three deciduous hardwood species Acer mono Maxim., Tilia japonica (Miq.)Simonkai and Prunus ssiori Fr.Schm.dominate in a cool-temperate conifer-hardwood mixed forest in northern Japan [21].This study analyzes the population structure and dynamics of the four most dominant species based on size structure, gap regeneration and tree competition by using the 10-year data of permanent plots to examine how regeneration traits differ among the four most dominant species.

Study Site
This study was conducted in a cool-temperate conifer-hardwood mixed forest in Hokkaido, northern Japan (43˚1'N, 143˚38'E, 200 m above sea level).Annual precipitation, recorded at Honbetsu Weather Station (ca. 10 km from the study site), was 837 mm during 1993 to 2004.Mean monthly temperatures of the coldest month of January and the hottest month of August were -8.7˚C and 19.4˚C, respectively, with 6.2˚C of annual mean temperature.Forest floor was patchily covered with Sasa nipponica Makino.The forest was well reserved, and therefore, there were no anthropogenic effects on the vegetation.Takahashi [21] described the vegetation and stand structure of this forest.

Field Methods
Three plots were established in a conifer-hardwood mixed forest in 1994.Plot size was 40 × 100 m, 30 × 50 m and 60 × 100 m (total 1.15 ha).Each plot was sub-divided into 5 × 5 m quadrats.I determined canopy gaps in each 5 × 5 m quadrat in 1994; quadrats with less than 50% of canopy coverage were defined as canopy gaps in this study.All trees ≥ 1.5 cm in diameter at breast height (DBH) were tagged, and their DBH was measured in 1994.Census was conducted also in 1995, 1996, 1998, 1999 and 2004.Dead trees and newly recruited trees, growing to ≥1.5 cm DBH, during successive censuses were recorded.

Data Analyses
The population dynamics of the four dominant species (Abies sachalinensis, Acer mono, Prunus ssiori and Tilia japonica) was analyzed.I analyzed the effects of tree competition on DBH growth rates of target trees, applying the following equation [9] [22] [23] [24] [25] [26], where ADGR is the absolute diameter growth rate (mm⋅year −1 ) of a target tree during 1994 to 2004, and BA L is the cumulative basal area (cm 2 •m −2 ) of neighboring trees larger than the target tree.Coefficient a 0 , a 1 and a 2 are constants.Coefficient a 0 is a growth rate at DBH 0 cm without suppression from neighboring trees, a 1 means the rate of increase of growth rate with increasing DBH, and a 2 means the rate of decrease of growth rate with increasing neighboring trees if the coefficient is a minus value.These coefficients were compared among the four species to examine the growth traits.Target individual trees were divided into five arbitrary DBH classes chosen to arrange similar vertical strata (1.5 -4.9, 5.0 -9.9, 10.0 -19.9, 20.0 -39.9, ≥40.0 cm DBH).I defined the neighboring trees if they satisfied the following two conditions: trees belonging to the same or larger size class of the target individual tree because the primary interest of this analysis was to express competition for light; trees in the same 10 × 10 m quadrat in which the target tree was located [9] [22] [23] [24] [25] [26].Local crowding was calculated for each target individual tree.Multiple linear regression analysis was conducted for the ADGR equation.ADGR was ln-transformed before the analysis to improve homoscedasticity.
Tree mortality is a discrete event.A datum can have only the value 0 (live) or 1 (dead).Therefore, I used the logistic equation to model the probability of tree mortality of the four dominant species as: where M is the annual mortality during 1994 to 2004 (year −1 ), b X ′ is a linear combination of parameter b and independent variable X, and t is the observation period (t years).M ranges between 0 and 1, so percent mortality (% year -1 ) can be obtained by multiplying M by 100.This study examined the effects of DBH, growth rate and competition on the tree mortality.Size-dependent mortality sometimes shows a U-shaped pattern with higher mortality at saplings and large canopy trees because of suppression of saplings by tall canopy trees and because of senescence of large canopy trees [26] [27].Such a U-shaped mortality can be expressed by a quadric function of size.Degree of local crowding by larger neighboring trees on target dead trees can be expressed as the cumulative basal area of neighbors within the same 10 × 10 m quadrat in which the target dead tree was located.Several researchers also showed that tree mortality is higher for trees with lower growth rates, and therefore, growth rate can be used as an integrated measure of whole-plant carbon gain [18] [27] [28] [29] [30].Combining these effects, the hypothesized mortality model is: in 1994, and BA L is the total basal area (cm 2 m -2 ) of larger neighbors than the target dead tree.I estimated the parameters of the logistic equation using maximum likelihood methods.Backward stepwise regression was also used.The criterion for including and removing an independent variable was set as P = 0.1 [31].
In this paper, two life phases were defined: sapling (1.5 ≤ DBH < 10 cm) and canopy tree (DBH ≥ 10 cm).To examine the canopy-gap dependency of saplings and recruits (during 1994 to 2004) of the four dominant species, the percentage of gap saplings (or recruits) was determined.Data in 1994 were used for saplings.I conducted a randomization test to determine significant differences in the observed percentage of gap saplings (or recruits) from random distribution.The data set of saplings (or recruits) was re-assorted randomly.I calculated the percentage of gap saplings (or recruits) for 1000 random re-assortment of saplings (or recruits), keeping the spatial distribution of canopy gaps fixed.If the observed percentage of gap saplings (or recruits) fell in either of the 2.5% tail of the calculated-value-distribution, the observed percentage of gap saplings (or recruits) was assumed to be significantly different from random.This calculation was conducted on the basis of 5 × 5 m quadrats.

Results
Total 30 hardwood species and a conifer Abies sachalinensis, larger than 1.5 cm DBH, were recorded in the three plots.The four species (Abies sachalinensis, Acer mono, Prunus ssiori and Tilia japonica) occupied 72% of the total density (1677 trees•ha −1 ) in 1994.The four species varied in the maximum DBH, the DBH frequency distribution and percentage of sprouts.24.4% of trunks of Tilia japonica was sprouts, while the percentage of sprouts was only ca.1% (Acer mono and Prunus ssiori) and 0% (Abies sachalinensis).Abies sachalinensis showed a bell-shaped DBH frequency distribution, while the three other species showed L-shaped patterns (Figure 1).Canopy gaps occupied 19.4% of the plots in 1994.Saplings and recruits of Abies sachalinensis significantly distributed more under closed canopies than in canopy gaps (Table 1).The number of recruits was quite low in Prunus ssiori and Tilia japonica during the census period, and their saplings were significantly   concentrated in canopy gaps (Table 1).The number of recruits of Acer mono was greater than the other three species, and they were concentrated in canopy gaps (Table 1).Saplings of Acer mono were also more concentrated in canopy gaps than the other three species (Table 1).
Absolute diameter growth rates (ADGR) of the four dominant species increased with DBH (Figure 2), and were decreased by the increase of local crowding of neighbors (Figure 3).The degree of the reduction of ADGR due to the increase of local crowding (a 2 ) did not differ among the four dominant species because the 95% confidence interval overlapped among the four species (Figure 3).However, the two other parameters (a 0 and a 1 ) varied among the four  species.Initial ADGR at DBH 0 cm (a 0 ) was lower in Acer mono and Prunus ssiori than Abies sachalinensis and Tilia japonica (Figure 3).The increase of ADGR with DBH (a 1 ) was greater in the subcanopy Prunus ssiori than the three other canopy species, while that of Acer mono was lowest among the four species (Figure 2 and Figure 3).
Mortality of Acer mono was a U-shaped pattern against DBH, while that of the three other species monotonically decreased with DBH (Figure 4).Although the mortality of the four species showed no relationship with local crowding of larger neighbors (BA L ), the mortality of the three species, except for Acer mono, increased with decreasing ADGR (Table 2).

Discussions
The four dominant species showed different regeneration traits.Abies sachalinensis is shade tolerant species, and forms sapling bank under closed canopy in Figure 4. Relationship between mortality and DBH at five DBH classes (1.5 -4.9, 5.0 -9.9, 10.0 -19.9, 20.0 -39.9 and over 40 cm) for the four dominant species.subalpine coniferous forests [32].In this stand, the saplings and recruits of Abies sachalinensis were more concentrated under closed-canopy conditions than in canopy gaps.This suggests that Abies sachalinensis did not depend on canopy gaps for its recruitment.However, the bell-shaped frequency distribution of DBH suggests that many Abies sachalinensis regenerated in a certain period and that regeneration traits of Abies sachalinensis are different between this conifer-hardwood mixed forest and subalpine coniferous forests.It is reported that the discontinuous age structure of Abies sachalinensis in a conifer-hardwood mixed forest in northern Japan [4].They suggested the importance of catastrophic windblown disturbances such as typhoon on the regeneration of Abies sachalinensis.Thus, Abies sachalinensis can grow and survive under closed-canopy conditions, but sporadic major disturbances accelerate the recruitment in conifer-hardwood forests like this study site.
The increase of ADGR with increasing DBH was greater in the subcanopy species Prunus ssiori than the three other dominant canopy species.Higher growth rate enables Prunus ssiori to reach the reproductive stage faster, compared with the other species.Leaf emergence and leaf fall of Prunus ssiori are considerably earlier than those of other deciduous broad-leaved species [33].
Thus, Prunus ssiori can receive much light before the leaf emergence of the other tall deciduous species, which increases the carbon gain and growth rate of Prunus ssiori in the understory.
Although both Tilia japonica and Acer mono are deciduous canopy tree species, the two species differed in the growth and recruitment rates, and mode of reproduction.The ADGR of Tilia japonica was higher than that of Acer mono at any DBH.This difference is partly due to the difference in wood density.Wood density of Tilia japonica (0.325 g•cm −3 ) is much smaller than that of Acer mono (0.632 g•cm −3 ) (K. Takahashi, unpublished data).Lighter wood of Tilia japonica would contribute to its higher growth rate per unit dry mass, compared with Acer mono [34] [35].On the contrary, the number of recruits was less in Tilia japonica than Acer mono (3 versus 24 recruits during 1994 to 2004).However, one-fourth of Tilia japonica regenerated by vegetative reproduction of sprouts.
An advantage of sprouting is long-term persistence [36] [37] [38].Thus, the population of Tilia japonica is maintained by vegetative reproduction.Many recruits were observed in Acer mono and its recruits and saplings regenerated more in canopy gaps, compared with the other species.However, the gap-dependent regeneration does not mean that Acer mono is a shade-intolerant species.Shade-intolerant species often show a bell-shaped or bimodal size structure, while shade-tolerant species show a L-shaped size structure [39].The size structure of Acer mono was a typical L-shaped pattern in this stand, indicating a shade-tolerant species.Acer mono increases the height growth rate in small canopy gaps through plastic changes of crown architecture [40], suggesting that Acer mono is a small-gap specialist.
The four dominant species showed species-specific regeneration strategies, i.e., disturbance-dependent regeneration (Abies sachalinensis), fast growth rate (Prunus ssiori), vegetative reproduction (Tilia japonica), and regeneration in canopy gaps (Acer mono).It is suggested that the species-specific regeneration traits contribute to the coexistence of the four dominant species through different regeneration niches.
where coefficient b 0 , b 1 , b 2 , b 3 and b 4 are constants, DBH is the initial DBH (cm)

Figure 1 .
Figure 1.Frequency distributions of DBH for Abies sachalinensis, Acer mono, Tilia japonica and Prunus ssiori in a conifer-hardwood mixed forest, northern Japan, in 1994.Solid bars indicate the number of sprouts.

Table 1 .
Number of recruits (passing the minimum boundary of 1.5 cm DBH during 1994 to 2004) and saplings (1.5 ≤ DBH < 10 cm) in canopy gaps and under closed canopies in a cool-temperate conifer-hardwood mixed forest in northern Japan.The gap area ratio of the plots was 19.4%.

Table 2 .
Estimated parameters for the mortality equation of the four dominant species.Parameters selected by the backward stepwise regression were listed.