Enhancing Overland Flow Inﬁltration through Sustainable Well-Managed Thinning: Contour-Aligned Felled Log Placement in a Chamaecyparis obtusa Plantation

: Contour felling is a restoration method used to decrease overland ﬂow (OF) and soil erosion in the world. However, the impact of thinning and the placement of felled logs on OF remains inconclusive. Low ground cover and soil permeability promote OF in Chamaecyparis obtusa (Siebold et Zucc.) Endl plantations, making thinning a method for reducing runoff. We examined the relationship between OF and ground cover in a C. obtusa plantation in Japan. Event-based runoff was monitored in three plots from 2016 to 2021, with 40% thinning conducted in 2019. In plot T1, logs were randomly scattered, and, in T2, logs followed contour lines, while control plots stayed the same. After thinning, both treatment plots showed lower OF than the control plot. The ANCOVA test shows a signiﬁcant slope reduction in treatment plots compared to the control plot from pre-thinning to post-thinning (T1: 0.67 to 0.26, T2: 0.66 to 0.12, p < 0.001, Tukey HSD test). However, in plot T2, OF remained stable for two years post-thinning, afﬁrming the enduring effectiveness of contour-aligned log placement. This study backs the notion that aligning fallen logs with contour lines boosts long-term OF inﬁltration, supporting sustainable forest and soil management.


Introduction
Overland flow (OF) generated on steep mountains leads to sediment transport, landslides, and flooding during intense rainfall over a short period.OF occurs when the following three events occur together [1]: (1) saturation of the entire soil profile during a rainfall event, (2) more rainfall intensity than the infiltration capacity of the soil surface, and (3) delay in water infiltration into the soil because of severe soil water repellency.Vegetation covers, such as forests, reduce OF because they have a high soil water infiltration capacity [1].However, high OF has been reported even in forestlands in China [2] and in dense Chamaecyparis obtusa (Siebold et Zucc.)Endl (C.obtusa) plantations in Japan [3,4].C. obtusa is an evergreen coniferous tree with scaly leaves [5] and shallow root systems [6].Its litterfall occurs between October and November (fall).The scaly leaves of C. obtusa are fragile and easily fragmented into small pieces after falling to the ground.Therefore, the litter layer in C. obtusa plantations is thin because the fallen fragmented leaves are easily eroded by wind and OF [7].
Most mountainous areas are covered with forests in Japan, accounting for 67%, with C. obtusa plantations occupying 10% of the forest areas [8].Headwater forests play a critical role in regulating water resources in Japan.Previous studies have suggested several reasons for the high OF in C. obtusa plantations, such as localized redistributed rainwater (i.e., throughfall and stemflow; [9]), sparse understory due to poor light conditions [10,11], a thin litter layer [7], soil crusting [12], and soil water repellency [13].
Ecohydrological processes in forest ecosystems are driven by various environmental and artificial factors [14], such as forest thinning [15], which affects ecohydrological processes by reducing the leaf area index and canopy interception loss [16][17][18].The selective removal of trees, also known as thinning, is a common forest management practice used to improve tree growth [19] and develop understory vegetation [20].Thinning is performed at different intensities for various purposes [21].Commercial thinning aims to enhance timber quality, while non-commercial thinning seeks to improve the ecological and hydrological conditions within the forest [14].In dense Japanese C. obtusa plantations, thinning is often conducted non-commercially to improve tree growth and poor light conditions on the forest floor; subsequently, the felled logs are left on the forest floor.
In recent years, the practice of strategically placing felled logs parallel to contour lines in forest plantations has gained prominence in Japan as an innovative approach aimed at promoting sustainability.This method aims to curtail OF and mitigate soil erosion by bolstering ground resistance and augmenting water infiltration into the soil [10,[22][23][24].Notably, while contour felling has been recognized as a standard rehabilitation technique in European and USA forests following fire events [25][26][27], as well as in heavily disturbed soils due to machinery operations in Iran [28], its adoption in Japan remains relatively novel, particularly in the context of post-thinning management [22].Moreover, the outcomes of previous research efforts [10,[22][23][24] concerning the combined impact of thinning practices and the strategic placement of felled logs on OF have yielded inconclusive results.
In light of these knowledge gaps and the pressing need for sustainable land management strategies, our present study endeavors to shed light on the intricate relationship between thinning practices, the placement of felled logs, and their collective influence on OF dynamics.Our central hypothesis posits that conducting thinning operations while aligning the placement of felled logs parallel to contour lines-a practice akin to wellmanaged thinning-will yield more effective results in curtailing OF.To rigorously test this hypothesis, we conducted a meticulous comparative analysis, evaluating the OF characteristics of two distinct thinned plots characterized by differing felled log arrangements, in comparison to an adjacent control plot.Furthermore, we delve into the mechanisms underpinning the observed changes in OF subsequent to thinning and the strategic placement of felled logs, thus contributing to a comprehensive understanding of sustainable forest and soil management practices.

Study Site
The study was conducted in a C. obtusa plantation (catchment area: 1.45 ha) located in the Obora experimental forest site of Toyota City, Aichi Prefecture, Japan (35 • 16 N, 137 • 15 E; 585 m a.s.l., Figure 1a), from January 2016 to December 2021.The climate of the study area can be characterized by extremely hot and humid summers and moderately cold winters, and it is categorized as wet, warm, and temperate [29].The climate class is warm temperate and fully humid with hot summers (Cfa), according to the Köppen-Geiger climate classification [30].The mean annual precipitation and temperature are 1470 mm and 15.3 • C (1991-2020), respectively [31].Rainfall occurs mainly in May-June (rainy season) and September-October (typhoon season).The C. obtusa plantations were established in 1989.The geology of the study area is a Mesozoic late cretaceous layer with granite bedrock [32].The soil is brown forest soil [33] equivalent to cambisol [34] and covered with almost no understory (Figure 2), a less than 1 cm litter layer, and a ca. 3 cm root mat.A northwest-facing (N5°W) slope (37°) was selected for this study (Figure 1b).

Experimental Design
Three study plots (each 4 m width, 10 m length) adjacent to each other were deployed along the same slope (Figures 1b and 2).A four-meter buffer zone was placed between each plot.Plots T1 (area: 30.8 m 2 ) and T2 (area: 31.7 m 2 ) were established as treatment plots, and plot C (area: 32.1 m 2 ) as a control (Figures 1b and 2).Each plot was delimited by inserting the plastic sheet into the ground (Figure 2).Overland flow (OF) was guided to a tipping bucket rain gauge (Uizin UIZ-TB200, 200 mm/tip) through a gutter at the bottom boundary of each plot (Figure 2).A plastic sheet was installed on the ground to prevent direct rainfall input to the gutter (Figure 2).
Rainfall data were recorded via a tipping bucket rain gauge (Ota Keiki OW-34-BP, 0.5 mm/tip) at a nearby weather station (ca.500 m from the studied plots, Figure 1b).A rain event was divided into 6 h of the rain-free period with a minimum 5.5 mm amount of rainfall (i.e., the minimum amount of rainfall that generated OF in all three studied plots).The tipping times of OF and rainfall data were recorded using the same type of data logger (Onset HOBO Pendant ® Event Data Logger UA-003-64).The area of plots T1 and T2 was thinned (ca.40% tree ha -1 ) in early 2019, and plot C remained intact as a control The C. obtusa plantations were established in 1989.The geology of the study area is a Mesozoic late cretaceous layer with granite bedrock [32].The soil is brown forest soil [33] equivalent to cambisol [34] and covered with almost no understory (Figure 2), a less than 1 cm litter layer, and a ca. 3 cm root mat.A northwest-facing (N5 • W) slope (37 • ) was selected for this study (Figure 1b).mm and 15.3 °C (1991-2020), respectively [31].Rainfall occurs mainly in May-June (rainy season) and September-October (typhoon season).The C. obtusa plantations were established in 1989.The geology of the study area is a Mesozoic late cretaceous layer with granite bedrock [32].The soil is brown forest soil [33] equivalent to cambisol [34] and covered with almost no understory (Figure 2), a less than 1 cm litter layer, and a ca. 3 cm root mat.A northwest-facing (N5°W) slope (37°) was selected for this study (Figure 1b).

Experimental Design
Three study plots (each 4 m width, 10 m length) adjacent to each other were deployed along the same slope (Figures 1b and 2).A four-meter buffer zone was placed between each plot.Plots T1 (area: 30.8 m 2 ) and T2 (area: 31.7 m 2 ) were established as treatment plots, and plot C (area: 32.1 m 2 ) as a control (Figures 1b and 2).Each plot was delimited by inserting the plastic sheet into the ground (Figure 2).Overland flow (OF) was guided to a tipping bucket rain gauge (Uizin UIZ-TB200, 200 mm/tip) through a gutter at the bottom boundary of each plot (Figure 2).A plastic sheet was installed on the ground to prevent direct rainfall input to the gutter (Figure 2).
Rainfall data were recorded via a tipping bucket rain gauge (Ota Keiki OW-34-BP, 0.5 mm/tip) at a nearby weather station (ca.500 m from the studied plots, Figure 1b).A rain event was divided into 6 h of the rain-free period with a minimum 5.5 mm amount of rainfall (i.e., the minimum amount of rainfall that generated OF in all three studied plots).The tipping times of OF and rainfall data were recorded using the same type of data logger (Onset HOBO Pendant ® Event Data Logger UA-003-64).The area of plots T1 and T2 was thinned (ca.40% tree ha -1 ) in early 2019, and plot C remained intact as a control

Experimental Design
Three study plots (each 4 m width, 10 m length) adjacent to each other were deployed along the same slope (Figures 1b and 2).A four-meter buffer zone was placed between each plot.Plots T1 (area: 30.8 m 2 ) and T2 (area: 31.7 m 2 ) were established as treatment plots, and plot C (area: 32.1 m 2 ) as a control (Figures 1b and 2).Each plot was delimited by inserting the plastic sheet into the ground (Figure 2).Overland flow (OF) was guided to a tipping bucket rain gauge (Uizin UIZ-TB200, 200 mm/tip) through a gutter at the bottom boundary of each plot (Figure 2).A plastic sheet was installed on the ground to prevent direct rainfall input to the gutter (Figure 2).
Rainfall data were recorded via a tipping bucket rain gauge (Ota Keiki OW-34-BP, 0.5 mm/tip) at a nearby weather station (ca.500 m from the studied plots, Figure 1b).A rain event was divided into 6 h of the rain-free period with a minimum 5.5 mm amount of rainfall (i.e., the minimum amount of rainfall that generated OF in all three studied plots).The tipping times of OF and rainfall data were recorded using the same type of data logger (Onset HOBO Pendant ® Event Data Logger UA-003-64).The area of plots T1 and T2 was thinned (ca.40% tree ha -1 ) in early 2019, and plot C remained intact as a control plot.Thinning was performed using a chainsaw, minimizing any disruption to the soil surface within each treatment plot.As a result, the stand density of studied plots changed from 2018 to 2020 as follows: plot T1 changed from 1623 tree ha -1 to 974 tree ha -1 , plot T2 changed from 1299 tree ha -1 to 649 tree ha -1 , and the control plot remained unchanged (974 tree ha -1 ) (Table 1).In addition, five felled logs (ca. 2 m in length and 30 cm in diameter) were left randomly on the slope (i.e., not parallel to contour lines) in plot T1, and another five felled logs were placed parallel to contour lines in plot T2 (Figure 2).Felled logs covered 5% of the ground in each treatment plot. 1 Field survey for tree height, DBH, and stand density was conducted in 2018 (pre-thinning) and 2020 (postthinning). 2OF and SD denote overland flow and standard deviation, respectively.

Statistical Analyses
Only OFs that were generated in all three studied plots were included in the data analyses.In this way, we can compare the event-based OF of treatment plots with the control plot.The data loss of OF during the thinning operation was excluded from the analyses (Figure 3).The relationship between the event-based OF and rainfall characteristics was investigated by performing a Pearson's correlation (r) analysis to understand the OF generation mechanism in each plot in the pre-and post-thinning periods.Rainfall amount, intensity, duration, and antecedent precipitation index or soil moisture conditions were selected as rainfall characteristics.
Antecedent precipitation index (API) was calculated for the rain events to understand the surrogate antecedent soil moisture conditions of the studied site prior to each OF event during the study period (Figure S1).We defined the API for a particular day as follows: where API d is the antecedent precipitation index for the day d, k is an empirical decay factor less than one, and P d is rainfall for the day d.k is usually between 0.85 and 0.98.A constant value of 0.95 was recommended by Hill et al. [35].In this study, we used the API for the 7th and 14th days before each OF event.We used these data for Pearson's correlation (r) analysis.Further, API7 and API14 represented the soil surface and subsurface moisture conditions, respectively.We performed analysis of covariance (ANCOVA) to detect differences in even-based OF between the thinned and control plots in the pre-and post-thinning periods.Post hoc comparisons with Tukey's honestly significant difference (HSD) tests were performed to determine whether the slopes of the regression lines differed between years (in the preand post-thinning periods).Finally, the slopes of the regression lines were compared to detect differences between the thinned and control plots.The ANCOVA results showed the impact of thinning and felled logs on OF in the thinned plots compared with the control plot.The R software ver.3.4.2(R Development Core Team, Vienna, Austria) was used for the statistical analyses.Antecedent precipitation index (API) was calculated for the rain events to understand the surrogate antecedent soil moisture conditions of the studied site prior to each OF event during the study period (Figure S1).We defined the API for a particular day as follows: where APId is the antecedent precipitation index for the day d, k is an empirical decay factor less than one, and Pd is rainfall for the day d.k is usually between 0.85 and 0.98.A constant value of 0.95 was recommended by Hill et al. [35].In this study, we used the API for the 7th and 14th days before each OF event.We used these data for Pearson's correlation (r) analysis.Further, API7 and API14 represented the soil surface and subsurface moisture conditions, respectively.We performed analysis of covariance (ANCOVA) to detect differences in even-based OF between the thinned and control plots in the pre-and post-thinning periods.Post hoc comparisons with Tukey's honestly significant difference (HSD) tests were performed to determine whether the slopes of the regression lines differed between years (in the preand post-thinning periods).Finally, the slopes of the regression lines were compared to detect differences between the thinned and control plots.The ANCOVA results showed the impact of thinning and felled logs on OF in the thinned plots compared with the control plot.The R software ver.3.4.2(R Development Core Team, Vienna, Austria) was used for the statistical analyses.

Overall OF
In total, 203 rainfall storms resulted in OF.In the pre-thinning period, 96 rainfall storms occurred (36,35, and 25 in 2016, 2017, and 2018, respectively), and 107 occurred during the post-thinning period (27, 45, and 35 in 2019, 2020, and 2021, respectively).During the pre-thinning period, 41%, 32%, and 27% of the rainfall event amounts were <25, 25-50 mm, and >50 mm, respectively.In the post-thinning period, 36%, 36%, and 28% of the rainfall event amounts were <25, 25-50, and >50 mm, respectively.During the pre-thinning period, the monthly OF was almost identical among the three plots (Figure 3).Contrastingly, in the post-thinning periods, the monthly OF in plot C was higher than that in plots T1 and T2 (Figure 3).Plot T2 had the lowest monthly and total OF in the post-thinning period compared to plots T1 and C (Figure 3 and Table 1).The OF frequency also reduced from the pre-to post-thinning periods in both thinned plots (Table 1).
Figure 4 illustrates the general trend in OF changes across three study plots from the pre-thinning to the post-thinning periods.During the pre-thinning period, all studied plots had similar flow-duration curves (FDCs, Figure 4, total, 2016, 2017).However, in the first year after thinning (2019), the FDC of plot C was higher than that of plots T1 and T2 (Figure 4, 2019).The FDCs were distinctly separated in the second year after thinning (2020), with plots T2, T1, and C having the lowest, moderate, and the highest FDCs, respectively (Figure 4, 2020).In the third year after thinning (2021), the FDC of plot T2 remained the lowest, whereas that of plot T1 shifted to plot C (Figure 4, 2021).

Pre-Thinning Period
The rainfall amount was positively correlated with OF in all plots (Table 2).Additionally, rainfall duration was positively correlated with OF in plot C (Table 2).Contrastingly, API7 and API14 were negatively correlated with OF in plot T1 (Table 2 and Figure S1).
plots had similar flow-duration curves (FDCs, Figure 4, total, 2016Figure 4, total, , 2017)).However, in the first year after thinning (2019), the FDC of plot C was higher than that of plots T1 and T2 (Figure 4, 2019).The FDCs were distinctly separated in the second year after thinning (2020), with plots T2, T1, and C having the lowest, moderate, and the highest FDCs, respectively (Figure 4, 2020).In the third year after thinning (2021), the FDC of plot T2 remained the lowest, whereas that of plot T1 shifted to plot C (Figure 4, 2021).

Pre-Thinning Period
The rainfall amount was positively correlated with OF in all plots (Table 2).Additionally, rainfall duration was positively correlated with OF in plot C (Table 2).Contrastingly, API7 and API14 were negatively correlated with OF in plot T1 (Table 2 and Figure S1).

Post-Thinning Period
During the post-thinning period, the rainfall amount was positively correlated with OF in all plots, with a higher correlation coefficient (r) value than in the pre-thinning period (Table 2).Further, in this period, APIs and OF were not correlated in plot T1 (Table 2).Additionally, rainfall intensity was positively correlated with OF in plot T2 (Table 2).Further, the relationship between rainfall characteristics and OF remained unchanged in plot C in this period compared to that in the pre-thinning period (Table 2).

Impact of Thinning and Placement of Felled Logs on OF
The relationship between OF in thinned plots and pre-and post-thinning in plot C is shown in Figure 5.The ANCOVA results showed that in both thinned plots, the slope of the regression line in the post-thinning period and in the first year after thinning (2019) was lower than that in the pre-thinning period (Figure 5c,d and Table 3).Differences in OF between the thinned plots began in the second year after thinning (Figure 5c,d and Table 3).The regression line for the second year post-thinning (2020) was higher than that for the first year post-thinning (2019) in plot T1, and remained unchanged from the second to the third year post-thinning (Figure 5c and Table 3).Contrastingly, in plot T2, the slope of the regression lines remained unchanged from the first year post-thinning until the end of the study period (Figure 5d and Table 3).Table 3. Results of covariance (ANCOVA 1 ) analysis between the slope of regression plots in pre-thinning and post-thinning periods.

Discussion
The study findings showed that OF was affected by the thinning and placement of felled logs.Previous studies have demonstrated the importance of ground cover in increasing ground resistance and reducing OF in forests [36,37].Their results indicated that OF increased with the removal of ground litter layers [36,37].Furthermore, dos Santos et al. [38] observed that after thinning a native dry tropical forest in Brazil, a welldeveloped herbaceous layer acted as a barrier, reducing surface runoff.Wagenbrenner et al. [26] reported the effectiveness of contour felling in reducing the peak runoff rates in a post-fire rehabilitation treatment in the Colorado Front Range.Our results showed that in the thinned C. obtusa plantations with low ground cover, the remaining felled logs played a vital role in maintaining low OF by increasing water residence time and providing an opportunity for the generated OF to infiltrate the soil.Furthermore, differences in microtopography could have caused differences in the OF generation mechanism in each plot during the pre-thinning period (Table 2).Sato et al. [39] experimented with small traps (for sediment, litter, and OF measurements) in the same C. obtusa catchment, and suggested that the high OF variability in a C. obtusa plantation could be affected by the microtopography or internal heterogeneity between the studied plots.
The thinned plots responded to OF reduction by changing the relationship between rainfall characteristics and OF from the pre-thinning to the post-thinning period (Table 2).Nevertheless, the limited correlation observed between rainfall characteristics and OF suggests the presence of additional factors influencing OF generation within the treatment plots (Table 2).In plot T1, the random placement of felled logs significantly affected the OF generation mechanism by eliminating the impact of API7 and API14 in the postthinning period (Table 2).These results indicated that soil moisture was an influential factor in generating OF in plot T1 prior to thinning.Contrastingly, the OF generation mechanism changed when the felled logs were randomly placed on the ground during the post-thinning period.In plot T2, rainfall intensity was a new parameter influencing OF during the post-thinning period (Table 2), thus supporting our hypothesis that the parallel placement of felled logs on contour lines changes the OF generation mechanism in the post-thinning period (Table 2).Kuraji et al. [24] also mentioned that the parallel placement of felled logs decreased the peak OF and increased subsurface runoff after the thinning of a C. obtusa stand compared to that in a grassland area.In our study, thinning and placing felled logs decreased OF in the treatment plots compared to the control plot with no thinning.Cheng et al. [40] reported that surface runoff in a heavily thinned Larix principisrupprechtii forest in China depended on the rainfall amount and intensity and litter interception.In plot C, the ground cover remained unchanged from the pre-thinning to post-thinning periods (Figure 2); therefore, the OF generation mechanism in this plot was similar during the post-thinning period (Table 2).
The FDCs and ANCOVA results showed the effectiveness of thinning in reducing OF following thinning in 2019 (Figures 4 and 5, and Table 3).The FDCs demonstrated the changes in the flow frequency (i.e., OF occurrences) of thinned plots compared to that of plot C from the pre-thinning to post-thinning periods (Figure 4, total).Additionally, a shift in FDC in plot T1 toward plot C in the second year after thinning suggested that the impact of thinning was longer in plot T2, with the felled logs placed parallel to the contour lines in the post-thinning periods (Figure 4, 2020).The differences in the FDCs between plots T1 and T2 were related to the direction of the felled log.Placing felled logs (or any barrier) on a slope helps to infiltrate water into the soil and stop sediments from being carried by OF [41].In the first year after thinning, felled logs in both thinned plots caused sediment deposits on the upper sides of the felled logs.The sediment deposited from parallel-felled logs in plot T2 enhanced OF infiltration into the soil.In plot T1, the random placement of felled logs negatively impacted OF reduction from the second year after thinning because the deposited sediments directed the OF downslope rather than infiltrating into the soil.Jourgholami et al. [28] reported a decrease in runoff and sediment yield when felled logs were aligned with contour lines in a natural broad-leaved forest in northern Iran disturbed by a heavy operating machine.In contrast, Fernández and Vega [25] observed that erosion barriers (i.e., contour felling with a trench on its upslope) were effective for a relatively short period (a few months) in controlling sediment and OF after severe wildfires in Spain.The discrepancy between these results and our results could be related to the high level of disturbance after the wildfire compared to the well-managed thinning in our study.
The ANCOVA results also showed that OF remained low in plot T2 in the second and third years after thinning, whereas in plot T1, OF increased in the second year after thinning (Figure 5c,d).Kubota et al. [42] also reported an increase in runoff in the second year after the 50% thinning of a small catchment of C. obtusa, where felled logs were left randomly.
Our findings underscore the effectiveness of well-managed thinning coupled with the strategic placement of felled logs parallel to contour lines in controlling OF within plantations, particularly in the context of C. obtusa where the understory is underdeveloped.By interrupting the movement of water across the soil surface, this approach demonstrated its capacity to efficiently manage OF dynamics.Furthermore, the sustained alteration of soil surface hydrology, achieved through the parallel arrangement of felled logs following thinning, highlights the potential for enduring impacts on water movement patterns.
While our study provides compelling support for the positive outcomes of thinning and felled log placement in reducing OF, it is imperative to acknowledge the long-term considerations inherent in sustainable land management.For a comprehensive grasp of the broader implications, it is recommended that future research delves into the influence of decomposing felled logs on soil hydraulic properties, and their role in rainfall-runoff processes over extended periods.Additionally, there exists a valuable opportunity to disentangle the distinct impacts of thinning practices and felled log placement by incorporating a thinned plot devoid of felled log interventions.Such multifaceted investigations will contribute to refining our understanding of sustainable forestry practices aimed at minimizing overland flow.

Conclusions
This study investigated the effects of the thinning and placement of felled logs on OF generation in a C. obtusa plantation in Japan.This study is the first to compare the impact of thinning on OF by considering the placement direction of felled logs (i.e., parallel versus random).We hypothesized that the thinning and parallel placement of felled logs would efficiently reduce OF.The results showed lower total and monthly and event-based OF and runoff frequency in the plot with a parallel placement of felled logs than in the plot with a random placement of felled logs in the post-thinning period, thus supporting our hypothesis.Our findings showed the importance of ground cover (i.e., felled logs) and its direction (i.e., random placement or parallel placement to the contour line) in controlling OF by altering soil surface hydrology after thinning a C. obtusa plantation.However, future research is necessary to improve the OF reduction mechanism using parallel-felled log placement.In the future, the impacts of thinning and felled logs on OF could be separately assessed by comparing the OF between thinned plots with and without felled

Figure 1 .
Figure 1.Map of the study site (a) and plots (b) in the Obora experimental forest site in Toyota City, Aichi Prefecture, Japan.

Figure 2 .
Figure 2. Condition of plot T1: felled logs placed in random directions on the slope; plot T2: felled logs placed parallel to contour lines on the slope; and plot C: control plot after thinning in 2019.

Figure 1 .
Figure 1.Map of the study site (a) and plots (b) in the Obora experimental forest site in Toyota City, Aichi Prefecture, Japan.

Figure 1 .
Figure 1.Map of the study site (a) and plots (b) in the Obora experimental forest site in Toyota City, Aichi Prefecture, Japan.

Figure 2 .
Figure 2. Condition of plot T1: felled logs placed in random directions on the slope; plot T2: felled logs placed parallel to contour lines on the slope; and plot C: control plot after thinning in 2019.

Figure 2 .
Figure 2. Condition of plot T1: felled logs placed in random directions on the slope; plot T2: felled logs placed parallel to contour lines on the slope; and plot C: control plot after thinning in 2019.

Sustainability 2023 , 12 Figure 3 .
Figure 3. Monthly precipitation and overland flow (OF) in plots T1, T2, and C in the pre-and postthinning periods.

Figure 3 .
Figure 3. Monthly precipitation and overland flow (OF) in plots T1, T2, and C in the pre-and post-thinning periods.

Figure 4 .
Figure 4. Flow-duration curve for plots T1, T2, and C in the pre-and post-thinning periods.

Figure 5 .
Figure 5. OF in plot T1 vs. plot C (OF T1 vs. OF C ) (a), plot T2 vs. plot C (OF T2 vs. OF C ) (b) in the pre-and post-thinning periods and OF T1 vs. OF C (c) and OF T2 vs. OF C (d) in pre-thinning and the years 2019, 2020, and 2021.

Table 1 .
General statistics of plots T1, T2, and C in pre-and post-thinning periods at the Obora experimental forest site.

Correlation Coefficient p-Value Correlation Coefficient p-Value
Significant correlations at a 95% confidence level.

Table 3 .
Results of covariance (ANCOVA 1 ) analysis between the slope of regression lines of studied plots in pre-thinning and post-thinning periods.