Edge effects confirmed at the clear-cut area of Korean red pine forest in Uljin, eastern Korea

Abstract

Background

Forest edges create distinctive ecological space as adjacent constituents, which distinguish between different ecosystems or land use types. These edges are made by anthropogenic or natural disturbance and affects both abiotic and biotic factors gradually. This study was carried out to assess edge effects on disturbed landscape at the pine-dominated clear-cut area in a genetic resources reserve in Uljin-gun, eastern Korea. This study aims to estimate the distance of edge influence by analyzing changes of abiotic and biotic factors along the distance from forest edge. Further, we recommend forest management strategy for sustaining healthy forest landscapes by reducing effects of deforestation.

Results

Distance of edge effect based on the abiotic factors varied from 8.2 to 33.0 m. The distances were the longest in Mg²⁺ content and total nitrogen, K⁺, Ca²⁺ contents, canopy openness, light intensity, air humidity, Na⁺ content, and soil temperature followed. The result based on biotic factors varied from 6.8 to 29.5 m, coverage of tree species in the herb layer showed the longest distance and coverage of shrub plant in the herb layer, evenness, species diversity, total coverage of herb layer, and species richness followed. As the result of calculation of edge effect by synthesizing 26 factors measured in this study, the effect was shown from 11.0 m of the forest interior to 22.4 m of the open space. In the result of stand ordination, Rhododendron mucronulatum, R. schlippenbachii, and Fraxinus sieboldiana dominated arrangement of forest interior sites and Quercus mongolica, Vitis amurensis, and Rubus crataegifolius dominated spatial distribution of the open area plots.

Conclusions

Forest interior habitat lies within the influence of both abiotic and biotic edge effects. Therefore, we need a forest management strategy to sustain the stability of the plant and further animal communities that depend on its stable conditions. For protecting forest interior, we recommend selective logging as a harvesting method for minimizing edge effects by anthropogenic disturbance. In fact, it was known that selective logging contributes to control light availability and wind regime, which are key factors affecting microclimate. In addition, ecological restoration applying protective planting for the remaining forest in the clear-cut area could contribute to prevent continuous disturbance in forest interior.

Background

Forest edges create distinctive ecological space as adjacent constituents, which distinguish between different ecosystems or land use types. These edges are made by anthropogenic or natural disturbance and affects both abiotic and biotic factors gradually (Harper et al. 2004). Edge effects can have serious impacts on species diversity and composition, community dynamics, and ecosystem functioning (Saunders et al. 1991, Laurance et al. 2006). Many edge effects are variable in space and time (Ewers and Didham 2006).

The influence of the adjacent non-forest environment on forest structure and species composition at created edges is now widely recognized. The altered habitat may be contributing to forest degradation and the loss of biodiversity in fragmented landscapes (Saunders et al. 1991, Gascon et al. 2000, Laurance et al. 2002).

Forest edges are becoming more abundant in many regions around the world because of the loss of forest arising from human activity, including settlement, agriculture, resource extraction, and timber harvesting. Most of all, clear cutting is a controversial forest management tool. It can create early successional habitats and edge areas at a landscape scale, which are preferred by many species (Pykälä 2004), while creating also a large portion of the landscape, which experiences edge effects (Keenan and Kimmins 1993). The edge effect has been a major topic of interest in studies of the landscape patterns and processes associated with edge creation and fragmentation during the last few decades (Harper et al., 2005).

Negative effects of edge creation have become apparent, including structural damage (Laurance et al. 1998) and depressed breeding success of songbirds (Gates and Gysel 1978) at forest edges. More recently, studies on changes at forest edges have revealed that edge effect can lead to the degradation of forest fragments (Gascon et al. 2000, Laurance et al. 2002). These negative consequences have fostered much interest in edges and fragmentation in conservation biology (Harper et al., 2005).

The strength of edge effects diminishes as one moves to deeper inside forests, but in addition, many edge phenomena vary markedly even within the same habitat fragment or landscape. Factors that might promote edge-effect variability include the age of habitat edges (Matlack 1993), edge aspect (Turton and Freiburger 1997), fragment size (Ewers et al. 2007), the structure of the adjoining matrix vegetation (Pohlman et al. 2007), seasonality (Young and Mitchell 1994), extreme weather events (Laurance et al. 2001), and fires (Cochrane and Laurance 2002).

To assess edge effects on disturbed landscape, several studies have estimated the magnitude of edge influence along the distance from forest edge (Harper and Macdonald 2011, Harper et al. 2015) and developed statistical methodology for calculating the distance of edge influence (Wales 1972, MacQuarrie and Lacroix 2003).

The majority of current researches on edge effect focused on proving differences along an interior-edge-exterior gradient and changes of species composition after clear cutting into early successional stage in Korea (Park et al. 2010, Ming et al. 2013, Kim 2010, Kim 2014). These researches designate arbitrary blocks and thus exclude a series of continuous environmental gradients. Any studies on the distance of edge influence, which a given variable is significantly different from the interior forest, were not carried out to date.

First of all, we aim to analyze changes of abiotic and biotic factors along the distance from forest edge in the forest area disturbed by clear-cut. Second, we aim at estimating the distance of edge influence by synthesizing those data. Based on our finding, we recommend forest management strategy for reducing effects of deforestation and thus recovering healthy forest landscapes.

Methods

Site description and data collection

The study was conducted in a harvest area of forest genetic resource reserves in Uljin-gun, eastern Korea. Mean annual precipitation and temperature of this region are 1119 mm and 12.6 °C, respectively (Korea Meteorological Administration 2011). The geology is Precambrian granite gneiss (Park and Yoon 1968). Plants of northern and southern provinces appear together due to topographical features connected to the high mountains called “Baekdudaegan.” The forest landscape of this study area is dominated by Korean red pine (Pinus densiflora for. Erecta) forest and Mongolian oak (Quercus mongolica) forest. The former is dominated by Korean red pine and Acer pseudo-sieboldianum, Rhododendron schlippenbachii, Rhododendron mucronulatum, and so on, appear in undergrowth. The latter is dominated by Mongolian oak and Styrax obassia, Lespedeza cyrtobotrya, Lindera obtusiloba, and so on, grow as undergrowth.

The study area is located in 37° 1′ 56.22″ N and 129° 12′ 4.65″ E geographically and at an elevation of approximately 690 m above sea level on NE aspects topographically (Fig. 1). Harvesting was conducted by applying a clear-cut method in 2013. Areal size that clear-cut is executed amounts to 1.54 ha. We carried out a field survey from May to September 2016. Spatial range for field survey includes the forest edge and interior and exterior of the forest.

Fig. 1

A map showing the geographic location of study area

We established five 80-m-long belt transects with 5-m breadth perpendicular to forest edge, which are extended to 40 m toward both directions of the forest interior and the harvest area. Belt transects were placed next to each other. Along each transect, we established 5 m × 5 m plots at the following distances from the edge: −38 m (C5), −25 m (C4), −16 m (C3), −9 m (C2), −3 m (C1), 3 m (O1), 9 m (O2), 16 m (O3), 25 m (O4) and 38 m (O5) (negative distances represent them within forest interior) (Fig. 1).

Measurement of abiotic factors

Microclimate

We installed meteorological sensors (Model HOBO Pro v2 U23–001, U23–004, Onsetcompany, Bourne, USA) on 1 m above ground and under 15 cm below ground at each distance to measure air temperature, air humidity, and soil temperature during the growing season. Each sensor was protected from direct sunlight by a protective case and observations were made at 30-min intervals. Soil humidity was measured during the same period with the other microclimatic factors and observed at 10-min intervals with a moisture probe (Pico Trime HD2, AGEON, Germany).

Canopy openness and light availability

Canopy openness and light availability were calculated using imaging software (GAP Light Analyzer 2.0) to extract canopy structural parameters and gap light transmission indices from true-color fisheye photographs (Nikon D80 digital camera, Sigma 4.5, F 2.8 EX DC Circular Fisheye HSM lens). All hemispherical fisheye photographs were taken from 1 m above ground at the center of plots.

Physicochemical properties of soil

Soil samples were collected from plots at all distance. In each distance, soil samples were collected in June 2016 from the top 10 cm at five random points in each plot, pooled, air-dried under the shaded condition, and sieved though 2-mm mesh. pH was measured in a 1:5 mixture of soil and distilled water then measuring after agitation using a pH meter (MAPA 1994). Total nitrogen and available phosphate were measured according to the Walkley and Black method (Walkely and Black 1934) and Lancaster method (Murphy and Riley 1962). Cation exchange capacity (CEC) and exchangeable bases (K⁺, Na⁺, Ca²⁺, and Mg²⁺) were determined by the ammonium acetate method (Chapman 1965) and atomic absorption spectrophotometry (MAPA 1994).

Soil respiration

Soil respiration was measured by applying the closed dynamic chamber method (Bekku et al. 1997) using a portable closed chamber infrared gas analysis system (SRC-1 with EGM-4, PP-Systems, Hitchin, Herts, UK). We installed six soil collars in plots of each distance. The CO₂ efflux was calculated from the concentration change within the chamber headspace and measured at 2-s intervals for 2 min by applying the following equation (Eq. 1). Soil respiration was measured periodically in every month for 3 months from June to August. Measurement was carried out at four times duplication from 10:00 to 13:00 in all plots of each distance.

$$ \mathrm{Soil}\ \mathrm{respiration}\ \left({\mathrm{m}\mathrm{gCO}}_2\kern0.24em {\mathrm{m}}^{-2}{\mathrm{h}}^{-1}\right)=3600\times {10}^{-6}\ a\ \rho\ V\ {S}^{-1} $$

(1)

$$ a:\mathrm{change}\ \mathrm{of}\ {\mathrm{CO}}_2\kern0.24em \mathrm{concentration}\ \left(\upmu \mathrm{mol}\ {\mathrm{mol}}^{-1}{\mathrm{s}}^{-1}\right) $$

$$ \rho :{\mathrm{CO}}_2\kern0.24em \mathrm{density}\ \left(\mathrm{mg}\ {\mathrm{m}}^{-3}\right) $$

$$ V:\mathrm{volume}\ \mathrm{of}\ \mathrm{the}\ \mathrm{chamber}\ \left({\mathrm{m}}^3\right) $$

$$ S:\mathrm{surface}\ \mathrm{area}\ \left({\mathrm{m}}^2\right) $$

Measurement of biotic factors

Species diversity

Shannon diversity index (H′) (Shannon and Weaver 1949, Brower and Zar 1984), richness (R) (Margalef 1958), and evenness (J) (Pielou 1969) were calculated to evaluate changes of species diversity along the distance from forest edge using PC-Ord 4.0 (McCune and Mefford 1999). Diversity indices were calculated by applying the following equations.

$$ {H}^{\hbox{'}}=-{\sum}_i{P}_i\ln \left({P}_i\right) $$

(2)

$$ R=\left({n}_i-1\right)/\ln (N) $$

(3)

$$ J={H}^{\hbox{'}}/{H}^{\hbox{'}}\max $$

(4)

$$ i:1,2,3,\dots S $$

$$ {P}_i:\mathrm{proportion}\ \mathrm{of}\ \mathrm{in}\mathrm{dividuals}\ \mathrm{in}\ \mathrm{the}\;i\mathrm{th}\ \mathrm{species} $$

$$ {n}_i:\mathrm{number}\ \mathrm{of}\ \mathrm{in}\mathrm{dividuals}\ \mathrm{in}\ \mathrm{the}\;i\mathrm{th}\ \mathrm{species} $$

$$ N:\mathrm{total}\ \mathrm{number}\ \mathrm{of}\ \mathrm{individuals} $$

Species composition

A crown projection diagram of herb layer was drawn for the whole study plot, which is extended over 80 m × 25 m divided with a 5 m × 5 m grid. All plant species appeared in each subplot were identified, following Lee (1985) and spatial range of crown that each plant covers the ground was measured using tape ruler in four directions of east, west, south, and north. The coverage of herb layer was calculated as a real size of land covered by the crown using ArcGIS program (ver. 10.0).

Vegetation survey was carried out Based on the frequency and dominance (Braun-Blanquet 1964), in 1 m × 1 m subplots arranged in diagonal direction in each plot. Differences in species composition were analyzed with non-metric multidimensional scaling (NMDS) (Kruskal 1964) stand ordination based on Euclidean distance, which was performed using the function “vegan” of the R statistical package. The importance value of each species was fed in a matrix for NMDS ordination.

Aboveground biomass

Estimation of aboveground biomass in the herb layer was made on each subplot applying direct sampling method after vegetation survey. Aboveground of all plants appeared in each subplot was harvested, oven dried at 80 °C to constant mass (Drying Oven, Daeil Eng. DDO102), and weighted using an electronic scale (SHIMADZU EB-3200HU).

Edge effect

The functions “dist” and “hclust” of the R statistical package (ver. 3.3.3) were used to determine inner forest plot as reference information.

Magnitude of edge influence (abbreviated as MEI hereafter) for each response variable was calculated by applying the following equation (Harper and Macdonald 2011):

$$ \mathrm{MEI}=\left|\frac{\overline{e}-\overline{r}}{\overline{e}+\overline{r}}\right| $$

$$ \overline{e}:\mathrm{value}\ \mathrm{of}\ \mathrm{the}\ \mathrm{parameter}\ \mathrm{at}\ \mathrm{edge} $$

$$ \overline{r}:\mathrm{value}\ \mathrm{of}\ \mathrm{the}\ \mathrm{parameter}\ \mathrm{in}\ \mathrm{the}\ \mathrm{in}\mathrm{terior} $$

Distance of edge influence (abbreviated as DEI hereafter) is defined as the distance that a given response variable is significantly different from the interior forest. Distance of edge influence can be calculated based on change of MEI along the distance. We calculated DEI with piecewise regression models, which make two or more lines joined at unknown points, called “breakpoints.” Breakpoints can be used as estimates of thresholds (Toms and Lesperance 2003). Breakpoints are calculated using the function “segmented” of the R statistical package and displayed as dashed gray lines.

On the other hand, DEI of each environmental factor was calculated based on the actual measurement value of the factor.

Results

Determination of the range of forest interior

As the result of hierarchical clustering, the forest interior plot (C5), which is the farthest from the forest edge formed by clear cutting and the other plots of the forest interior were divided into different groups from each other (Fig. 2). Based on the result, ecological information of the forest interior plots except for C5 plot was used for calculating MEI and DEI to exclude the effects of the other environmental factors beside the edge effect.

Fig. 2

A result of hierarchical clustering based on the abiotic and biotic factors (C5: −38 m, C4: −25 m, C3: −16 m, C2: −9 m, C1: −3 m, O1: 3 m, O2: 9 m, O3: 16 m, O4: 25 m, O5: 38 m). The optimal solution with k = 4 is indicated by different coloring

Magnitude of edge effect based on the abiotic factors

As the result of calculation for the distance of edge effect by applying a regression analysis method (Fig. 3), soil temperature showed significant edge effect from 7.6 to 15.8 m in the open space by lumbering. Humidity showed significant edge effect from 1.6 to 13.6 m of the forest interior. Air temperature and soil moisture did not show any significant edge effect. Canopy openness showed significant edge effect from the borderline of the lumbering to 1.5 and 11.3 m toward the forest interior and the open area, respectively (Fig. 4). Light intensity showed edge effect from 4.7 m of the forest interior to 12.1 m of the open space (Fig. 4).

Fig. 3

Changes of soil temperature and air humidity along the distance from forest edge (soil temperature: r ² = 0.7318, air humidity: r ² = 0.9954, black thick line: piecewise regression, dotted line: Two break points)

Fig. 4

Changes of canopy openness and light intensity along the distance from forest edge (canopy openness: r ² = 0.9623, transmitted light: r ² = 0.9497, black thick line: piecewise regression, dotted line: two break points)

As the results of calculation based on physicochemical properties, total nitrogen, K⁺, Mg²⁺, and Na⁺ showed significant edge effects from 10.0 m, 14.1 m, 12.9 m, and 4.4 m of forest interior, to 21.4 m, 14.0 m, 20.1 m, and 4.0 m of the open space, respectively (Fig. 5). Ca²⁺ showed an edge effect from 2.5 m to 20.0 m of the open space. On the other hand, pH, available phosphorous, cation exchange capacity, and soil respiration did not show any significant effect.

Fig. 5

A change of physicochemical properties of soil along the distance from forest edge (total nitrogen: r ² = 0.6148, K⁺: r ² = 0.5490, Na⁺: r ² = 0.7622, Ca²⁺: r ² = 0.6207, Mg²⁺: r ² = 0.5942, black thick line: piecewise regression, dotted line: two break points)

Distance of edge effect based on the biotic factors

Species diversity, evenness, and richness showed the edge effects from 8.2 to 23.0 m in the open space, from 7.8 to 22.9 m in the open space, and from 3.2 m of the forest interior to 3.6 m of the open space (Fig. 6). Coverages of tree and shrub species in the herb layer showed edge effects from 7.6 m of the forest interior to 12.0 m of the forest interior (Fig. 7). On the other hand, total coverage of the herb layer showed edge effect from 1.4 to 11.9 m of the forest interior (Table 1). But coverage of subtree, vine, and herbaceous plants and aboveground biomass did not show any significant effect.

Fig. 6

Changes of species diversity, evenness and richness indices along the distance from forest edge (species diversity: r ² = 0.8442, species evenness: r ² = 0.5901, species richness: r ² = 0.9371, black thick line: piecewise regression, dotted line: two break points)

Fig. 7

A change of coverage of undergrowth along the distance from forest edge (tree species in herbal layer: r ² = 0.6984, shrub species in herbal layer: r ² = 0.6751, total vegetation area in herbal layer: r ² = 0.754, black thick line: piecewise regression, dotted line: two break points)

Table 1 Estimates of the edge effect distance and approximate 95% confidence intervals (CI) for abiotic and biotic factors (r ² > 0.6)

Distance of synthetic edge effect

As the result of calculation of edge effect by synthesizing 26 factors measured in this study, the effect was shown from 11.0 m of the forest interior to 22.4 m of the open space. In the result of calculation based on 15 factors, which showed significant effect, the effect was shown from 3.3 m of the forest interior to 4.0 m of the lumbered area (Fig. 8, Table 2).

Fig. 8

A change of magnitude of edge effect according to the distance from forest edge (left: with all factors, r ² = 0.9860; right: with significant factors, r ² = 0.9828, black thick line: piecewise regression, dotted line: two break points)

Table 2 Estimates of the distance of edge influence and approximate 95% confidence intervals (CI) (r ² > 0.6)

Response of species composition

As the result of stand ordination based on vegetation data of the herb layer collected in the forest interior, edge, and lumbered area, plots of the forest interior and the open area were divided into the left and the right parts on the AXISI, respectively, and sites of the edge were arranged between them. In the result of stand ordination, Rhododendron mucronulatum, R. schlippenbachii, and Fraxinus sieboldiana dominated arrangement of forest interior sites and Quercus mongolica, Vitis amurensis, and Rubus crataegifolius dominated spatial distribution of the open area plots (Fig. 9).

Fig. 9

NMDS biplot of sites (stress = 0.1724) and environmental variables. All species variables had significant impacts on community differences (based on 1000 permutations, p = 0.001)

Discussion

Edge effect on abiotic factors

Primary responses of ecological process emerge immediately due to changes in forest canopy structure. The direct effects of edge creation led to primary responses of biophysical process, which are significantly different between interior and exterior of forest (Harper and Macdonald 2001).

Among abiotic factors, air humidity and soil temperature stand as significant factors in edge effect. Air humidity began to decrease from 14 m within the forest and thereby reduced more than 3% in the forest edge compared with that of the forest interior. Air humidity was as low as 76% in the forest exterior compared with that of the forest interior. Soil temperature showed the opposite trend and thus increased continuously to 16 m in the open area (Fig. 3).

Air humidity changed rapidly but soil temperature showed a slightly different trend in edge area (Fig. 3). Although microclimate of forest interior and exterior has been compared previously (Geiger 1965, Lee 1978), it is important to consider the edge as both a separate microclimate and a climatic mediator between forest and clear-cut area. Forest edge has distinguished temperature and moisture regimes resulting from wind. As near the edge winds are weaker than in clear-cut area, this produces relative stable air, which allows more extreme temperature and humidity (Chen et al. 1993).

The pattern of microclimate apparently mirrors that of light exposure because solar radiation causes soil heating (Davies-colley et al. 2000). Our finding of edge effect on canopy openness and light intensity is closely comparable with that of Delado et al. (2007) who reported that the edge effects on canopy openness and light intensity may penetrate as much as 10 and 5 m into Canary island pine (Pinus canariensis) stand. The edge effects on canopy openness and light intensity are more influential in forest interior than in forest exterior, showing complete change over about 10 m exterior forest (Fig. 4).

Total nitrogen contents in forest interior and edge were similar to previous research, which was carried out in a Pinus densiflora–Quercus mongolica community of Korea (Korea Forest Service 2007). There was little difference in the average value between forest interior and edge, but the edge area showed a greater variation than forest interior. Clear-cutting influences nutrient contents, which are released from decomposing logging residues (Palviainen et al. 2004), soil microorganisms, and soil structure, which are occupied by tree roots (Keenan and Kimmins 1993).

Ca²⁺, Mg²⁺, and K⁺ contents showed a notable difference between forest interior and exterior, but no any difference between forest interior and edge. This result due to the lowered pH in forest interior because of pine trees composing of the forest (Hur and Joo 2002) was similar to that of Lippok et al. (2014) (Fig. 5).

Soil of forest interior have higher exchangeable Na⁺ than that of forest exterior. But other studies showed that there was no significant difference between forest interior and exterior (Johnson et al. 1991). Therefore, further studies are required in order to make a conclusive determination for edge effect of this element.

Edge effect on biotic factors

Species richness, diversity, and evenness in herb layer rapidly increased in edge area (Fig. 6). Edge effect on plant community tends to promote the establishment of generalist species and light-demanding species (Alignier et al. 2014) and consequently reconstruction of species composition due to abiotic changes (van Oijen et al. 2005). Although species diversity and evenness in herb layer decreased in forest exterior, species richness was usually higher in forest exterior. Canopy open increases open habitat species and consequently it contributed to increasing species richness with regeneration of dominant species (Fig. 6).

Coverage of tree and shrub in undergrowth increased due to regeneration of woody plant after clear-cutting (Wendel 1975, Johnson 1977). Total coverage of undergrowth was relatively in a steady state in the forest interior, whereas it increased when it was becoming far from the forest edge in the harvest area. Total coverage of undergrowth tended to increase rapidly after lumbering (Marozas et al. 2005, Fig. 7).

Forest management strategy depending on magnitude of edge effect

The responses of species to ecotones and correspondingly to edge effects occur on both sides of an edge located on forest interior and exterior (Zurita et al. 2012). Overall estimates of DEI toward forest interior and exterior was 11 and 22 m, respectively.

Species composition in forest edge showed a feature that forest interior and exterior are mixed. But there was a difference between forest edge and exterior as sprouts of woody plant (Q. mongolica) were generated from stump and thereby occupied high coverage on forest floor after clear-cutting. In addition, vines (R. crataegifolius and V. amurensis) preferring open and disturbed areas had significant impacts on change of species composition (Suzuki 1987, Krestov et al. 2015). Although other studies (Harper et al., 2005, McDonald and Urban 2006) showed that creating edge area induced invasion of exotic species, we did not find any exotic species in this study site as it is protected as a forest genetic resource reserve.

Forest interior habitat lies within the influence of both abiotic and biotic edge effects. Therefore, we need forest management to sustain the viability of the plant and animal communities that depend on stable conditions. We recommend selective logging as a harvesting method within 20 m from the remaining forest based on DEI obtained from this study for minimizing edge effects (Fig. 8). Selective logging in edge area toward forest exterior contributes to control light availability (Matlack 1993) and wind regime, which are key factors affecting microclimate (Davies-Colley et al. 2000, Fig. 8, Table 2).

In addition, protective planting could be applied as a measure to prevent continuous disturbance in forest interior in the clear-cut area.

Conclusions

First of all, this study was carried out to evaluate the magnitude of edge effect that lumbering caused. Further, we recommend a forest management strategy to sustain stable forest ecosystem by minimizing the effect. Based on the results of this study, lumber harvest changed abiotic factors such as physicochemical properties of soil, atmospheric humidity, and light intensity and thereby induced edge effect. Further, it also caused changes of biotic factors including species composition from forest margin into forest interior. The result of evaluation on the edge effect by synthesizing both abiotic and biotic factors clarified that the magnitude of edge effect reaches up to 11 m of forest interior. The result of stand ordination based on vegetation data represented that shade-intolerant species tended to increase on both forest edge and harvested sites. Selective cutting rather than clear-cutting was recommended as a harvesting method to reduce edge effect in the area bordered on the remaining forest. Further, ecological restoration applying protective planting was suggested as another forest management strategy to prevent continuous disturbance like invasion and expansion of exotic species following forest edge formation into the forest interior.