Long-Term Effects of Repeated Prescribed Fire and Fire Surrogate Treatments on Forest Soil Chemistry in the Southern Appalachian Mountains (USA)

: From 2001–2018, a series of fuel reduction and ecosystem restoration treatments were implemented in the southern Appalachian Mountains near Asheville, North Carolina, USA. Treatments consisted of prescribed fire (four burns), mechanical cutting of understory shrubs and mid-story trees (two cuttings), and a combination of both cutting and prescribed fire (two cuts + four burns). Soils were sampled in 2018 to determine potential treatment impacts for O horizon and mineral soil (0–10 cm depth) carbon (C) and nitrogen (N) and mineral soil calcium (Ca), magnesium (Mg), phosphorus (P), potassium (K), and pH. Results suggested that mean changes in O horizon C and N and mineral soil C, N, C:N, Ca, and P from 2001–2018 differed between the treatments, but only mineral soil C, N, C:N and Ca displayed differences between at least one fuel reduction treatment and the untreated control. One soils-related restoration objective was mineral soil N reduction and the cut + burn treatment best achieved this result. Increased organic matter recalcitrance was another priority, but this was not obtained with any treatment. When paired with previously reported fuels and vegetation results from this site, it appeared that continued use of the cut + burn treatment may best achieve long-term management objectives for this site and other locations being managed for similar long-term restoration and fuels management objectives.


Introduction
Wildland fires may impact soil physical, biological, and chemical properties in a multitude of ways [1][2][3][4][5]. These impacts may be related to factors such as fire intensity, severity, and frequency; vegetative species composition; soil sampling period (i.e. time since fire occurred); soil property investigated (i.e. O horizon properties vs. mineral soil properties); and many other factors [4][5][6]. These variables differ greatly between most wildfires and prescribed fires in a given ecosystem, therefore, soil responses to wildfires and prescribed fires may also differ [2,6,7]. Generally, lowintensity, low-severity prescribed surface fires impact soils less than high-intensity, high-severity wildfires [3,7,8]. Likewise, harvesting practices have also been shown to impact soils variably depending upon factors such as harvesting method, amount of material harvested, amount of mineral soil exposed, and topography [9][10][11].
On the Green River Game Lands near Asheville, North Carolina, USA, in the southern Appalachian Mountains, fuel reduction and woodland restoration treatments consisting of prescribed fire only, mechanical cutting of mid-story and understory vegetation, and a combination of these two treatments were first implemented in 2001 as part of the National Fire and Fire Surrogate Study (FFSS) [12][13][14]. Historic forests of the Appalachian region, dominated by upland oaks (Quercus spp.), hickories (Carya spp.), and pines (Pinus spp.), were nitrogen (N)-limited [15,16]. These ectomycorrhizal species thrived in soils with low inorganic N content and highly recalcitrant organic matter [17]. Many locations within the Appalachian region would now be considered N-saturated as a result of industrialization and atmospheric N deposition [18]. Therefore, low total ecosystem N and recalcitrant organic matter were desired, soils-related, treatment outcomes [15].
O horizon and mineral soil (0-10 cm depth) carbon (C), N, and C:N and mineral soil calcium (Ca), magnesium (Mg), phosphorus (P), potassium (K), and pH were evaluated 1 to 2 years following the first implementation of these treatments on the Green River Game Lands [13]. O horizon C:N and mineral soil Ca differed between treatments. O horizon C:N was lower on the cut + burn treatment than on all others and mineral soil Ca differed between the cut-only and cut + burn treatments and an untreated control. Two to four years post-treatment, however, differences between treatments for these variables were not present. Soil pH differed between the cut-only and untreated control treatments, but that difference was only 0.12. Therefore, the investigators suggested that treatmentrelated soil responses were time-sensitive and one treatment implementation could be expected to have little biological significance. Structural and functional restoration of this ecosystem was expected to require repeated treatments every 3-8 years and potentially higher intensity and severity of prescribed burns [15].
Since 2001-2002, additional treatments have been conducted on the Green River Game Lands: four total burns for the burn-only treatment; two total cuts for the cut-only treatment; and two total cuts + 4 burns for the cut + burn treatment [14]. Recent fuels and vegetation research following the repeated use of these treatments at this site suggested that: 1) all fuel reduction and restoration treatments reduced one or more fuel types, 2) the burn-only and cut + burn treatments reduced most fuels and most likely decreased wildfire hazard best, 3) desired upland oak and pine regeneration and graminoid density and cover increased following cut + burn treatments, and 4) cut + burn treatments induced forest structural changes most consistent with restoration goals, but no treatment regime had fully achieved the desired woodland structure [12,14]. No additional soils research from this site has been published, however.
Given the long-term and repeated use of these fuel reduction and land management strategies in the southern Appalachian Mountains, we sought to determine potential long-term treatmentinduced soil chemical changes at this site. In summer 2018, O horizon and mineral soil samples were collected. O horizon and mineral soil C, N, and C:N and mineral soil Ca, K, Mg, P, and pH were evaluated to determine changes in values of these properties since initial sampling was conducted in 2001. Based upon the additional treatments that have been conducted since 2001-2002, our hypotheses were: 1) mean differences in 2001 and 2018 O horizon and mineral soil C, N, and C:N values would differ most for the burn and cut + burn treatments (suggesting that the repeated use of these treatments has led to more N-limited soil conditions and more recalcitrant organic matter) and 2) mean differences in 2001 and 2018 mineral soil Ca, Mg, P, K, and pH values would be greatest for the cut + burn treatment.

Experimental Design and Treatment Descriptions
Treatments for this study consisted of prescribed burn only, mechanical cutting only (cutting of all mountain laurel and rhododendron stems and all other woody stems > 1.8 m tall and < 10 cm diameter at breast height (DBH) to create a vertical fuel break), a combination of mechanical cutting + prescribed burning, and an untreated control [12]. The burn-only and cut-only treatments were designed to reduce fuels with 80% overstory survivability in the event of a wildfire occurring on an 80 th percentile weather day [19]. For this study, a randomized complete block design was selected (4 treatments in 3 replicate blocks with blocks representative of separate, but similar geographic locations) ( Figure 1). All replicates consisted of a contiguous block of land, which helped control variability among sites. Each stand in the 3 replicates was approximately 14 ha in size and basal area ranged from 21-31 m 2 ha -1 prior to treatment. Stands selected for this study were approximately 80 to 120 years old and had at least 5 years without fire and 10 years without cutting prior to initial treatment installation [13].
The original cutting treatments were conducted for the cut-only and cut + burn treatments December 2001-February 2002. A second cutting operation for both treatments occurred January-February 2012 [9]. As of 2018, 4 dormant season prescribed burns had occurred in the burn-only and cut + burn units: February-March 2003, and 2015. These burns were ignited using aerial ignition (spot fire) or strip head fire techniques [14].

Soil Sampling, Processing, and Analyses
Ten randomly oriented 20 m × 50 m plots were established within each treatment-replicate block, resulting in 120 total treatment-replicate plots. These plots were permanently marked with rebar in each of the plot corners. Latitude and longitude coordinates were logged for each plot using GPS units. Within each plot, ten 10 m × 10 m subplots were created ( Figure 2). In 2001, O horizon samples were taken from six 10 m × 10 m subplots within each of the 120 plots using a shovel (Figure 3a). Six mineral soil samples (0-10 cm depth) were collected directly below locations where the O horizon samples were obtained to evaluate mineral soil chemistry [13]. These plots were visited often and sampled for additional ecosystem components, therefore, the soil sampling locations within each plot were oriented to minimize potential disturbance from other sampling activities ( Figure 2). All 120 plots were re-visited May-July 2018. Three O horizon grab samples were collected within 1 m of the 10 m x 10 m subplots in the 120 treatment-replicate plots. Three mineral soil samples (0-10 cm depth) were obtained in these exact locations using an Oakfield Model H soil probe (Oakfield Apparatus, Fond du Lac, WI, USA; inner diameter 2.06 cm) ( Figure 3b).  potassium (K) concentrations were derived using Mehlich III methodology [21] and resulting analyses for each element were conducted using ICP-Optical Emission Spectrometry [22]. Soil pH was determined using a 1:1 soil to water solution [23].

Statistical Analyses
Sample sizes differed between treatments as a result of failed collections for a few plots in either sampling year (Tables 1, 2, and 3). Outlier analyses were conducted, but no plot values were removed from the statistical analyses described in this section. Values for each treatment-replicate plot were entered into JMP Pro 14 (SAS Institute Inc., Cary, NC, USA). An analysis of variance (ANOVA) was conducted to determine differences between treatments during each sampling year. Differences between the 2018 and 2001 values were determined for each available treatment-replicate plot by subtracting the 2001 values from the 2018 values (subsequently referred to as ∆ values). An ANOVA was conducted to determine the significance of mean ∆ values between treatments. Additionally, mean differences between the 2018 and 2001 values were evaluated between the replications to evaluate spatial heterogeneity in soil chemical values. For each of these tests, least square means were determined and a Tukey's HSD test was used to separate means when differences were detected. Differences were declared statistically significant at α = 0.05.
The authors noted that mineral soil C, N, C:N, Ca, and Mg ∆ values were not normally distributed when evaluated with the Kolmogorov-Smirnov test. Multiple data transformations (i.e log, ln, square root, etc.) were not successful in resolving this issue. Non-parametric analyses were then conducted for these variables, but difference determinations between treatments were the same as those estimated using a normal distribution. Therefore, all results reported assume a normal distribution.

O horizon Carbon, Nitrogen, and C:N
Pre-treatment O horizon C and N differed between the assigned treatments (p < 0.0001) ( Table  1), indicating the high degree of soil spatial heterogeneity present within this landscape. In 2018, differences were detected between O horizon N treatment means (p = 0.0001), but not O horizon C treatment means (p = 0.0695). Tukey's HSD indicated O horizon N was greatest for the cut-only treatment (8.69 g kg -1 ) and differed only with the cut + burn (6.28 g kg -1 ) treatment, not the untreated control (7.96 g kg -1 ) or burn-only (7.12 g kg -1 ) treatments. The O horizon C:N mean ∆ values were positive for all of the fuel reduction and restoration treatments and negative for the untreated control. Although these were not statistically significant, this trend may be noteworthy as an indicator of potential increased organic matter recalcitrance. Mean ∆ values differed between the treatment replications for each of the O horizon variables (C p = 0.0002; N p = 0.0032; C: N p = 0.0122) ( Table 2), further highlighting the spatial variability present for these soil properties.

Mineral Soil Carbon, Nitrogen, and C:N
Pre-treatment mineral soil C (p < 0.0001), N (p = 0.0057), and C:N (p < 0.0001) differed between the assigned treatments in 2001 (Table 1), again indicating the high degree of spatial heterogeneity present within this landscape. In 2018, mineral soil C (p = 0.0074) and N (p = 0.0412) differed between the treatments, but C: N did not (p = 0.6152). Mean mineral soil C in 2018 was highest for the cut-only treatment (32.68 g kg -1 ) and untreated control (32.49 g kg -1 ) and these values differed from the cut + burn (24.81 g kg -1 ) values, but not the burn-only (27.54 g kg -1 ) values. Mean mineral soil N was also highest for the untreated control (1.22 g kg -1 ) and this differed from the cut + burn (0.92 g kg -1 ), but not the burn-only (1.04 g kg -1 ) or cut-only (1.16 g kg -1 ) means.
Mean ∆ values differed between treatments for C, N, and C:N (p < 0.0001, p = 0.0040, and p = 0.0111, respectively) ( Table 1). Tukey's HSD indicated the cut + burn values differed from both the untreated control and burn-only values for each variable, but they did not differ from the cut-only values. Mineral soil C mean ∆ values were positive for the untreated control and burn-only treatments and the cut + burn and cut-only mean ∆ values were negative. Similarly, the burn-only N mean ∆ value was positive and the cut + burn and cut-only N mean ∆ values were negative. The untreated control N mean ∆ value was 0.0 g kg -1 , however. The 2018 mean mineral soil N range was 0.92-1.22 g kg -1 N, a difference of only 0.30 g kg -1 . Despite this narrow range and potentially high variability in inherent soil properties at individual sampling locations, the overall 0.42 g kg -1 N reduction from 2001-2018 for the cut + burn treatment was highest between the treatments (a 31% reduction in N between 2001-2018). Mean ∆ values for mineral soil C, N, and C:N differed between the treatment replications (C and N p = < 0.0001; C:N p = 0.0254) ( Table 2), further highlighting the spatial variability present for these soil properties.

Mineral Soil Macronutrients and pH
Pre-treatment mineral soil P (p = 0.0006), K (p = 0.0115), and Ca (p = 0.0245) differed between the assigned treatments (Table 3). Differences between the 2018 treatment means were not detected for P, K, Ca, Mg, or pH (all p > 0.1688) ( Table 3). All mean ∆ values for these variables were negative. The reductions were substantial, accounting for approximately 50%-57% P reductions and 41%-73% Ca reductions, for example. Differences in sieving methods between 2001 (Sawyer mill) and 2018 (handsieved) may be related to both the direction and magnitude of these differences. Mean ∆ values between treatments were detected for mineral soil Ca (p = 0.0116) and P (p = 0.0115) ( Table 3). Tukey's HSD suggested the Ca untreated control values differed from the cut + burn and cut-only values. The cut + burn Ca reduction was only 37% of the untreated control Ca reduction. Phosphorus mean ∆ values differed between the burn-only and cut-only treatments with the burnonly accounting for the greatest reduction. No mean ∆ value differences were detected between treatments for Mg (p = 0.6269), K (p = 0.4787), or pH (p = 0.3537) ( Table 3). Magnesium (p = 0.0186) and K (p = 0.0003) mean ∆ values differed between the treatment replications ( Table 2), evidence of the spatial variability present for these soil macronutrients.

Temporal and Spatial Variability of Soil Properties and Fire Effects
Mean ∆ value differences between treatments were observed for O horizon C and N and mineral soil C, N, C: N, Ca, and P ( Table 1; Table 3). However, these differences existed only between some of the fuel reduction treatments, not the untreated control, for O horizon C and N (  Table 3). The cut + burn mean ∆ values differed between the untreated control and the cut-only treatment for mineral soil C and the untreated control and the burn-only treatment for both mineral soil N and C: N (Table 1). Both the cut + burn and cut-only Ca mean ∆ values differed from the untreated control (Table 3). The results for these variables suggested that changes induced by the periodic and repeated use of these treatments were greater than changes measured in the untreated control during this time period. When the 2018 ANOVA results were evaluated for these variables (Table 1), only mineral soil C and N for the cut + burn treatment differed significantly from the untreated control. Therefore, mineral soil C and N were significantly changed by the repeated use of the cut + burn treatment and these changes produced significantly altered mean values that differed from the untreated control mean values. With this in mind, the cut + burn treatment enacted the most effective change toward one of the desired fuel reduction and ecosystem restoration goals: N limitation [15]. Not only was the change in mineral soil N from 2001-2018 greatest for the cut + burn treatment, but the resulting 2018 mineral soil N mean was also significantly different from the 2018 untreated control N mean. The reduction in mineral soil N for the cut + burn treatment did not differ from the cut-only or burn-only treatments, however. A second soils-related restoration priority for the Green River Game Lands was increased organic matter recalcitrance [15]. No treatment successfully obtained this result, however, as O horizon C:N did not differ significantly between treatments (Table 1). However, all fuel reduction mean ∆ values were positive and the untreated control mean ∆ value was negative, suggesting that O horizon C:N was higher in 2018 than 2001 for all treatments. Perhaps that is a non-significant, yet noteworthy trend. Mineral soil C:N mean ∆ values differed between the cut + burn treatment and the untreated control and burn-only treatments (Table 1). However, when the 2018 mean values were compared, no differences existed between treatments and the range in C:N values was 27.35-28.78 or a difference of only 1.43 from highest to lowest. In this regard, any treatment-related C:N change most likely had little biological significance in the field and may be more related to inherent soil variability as opposed to actual treatment differences.
As noted in previous studies, soil physical, chemical, and biological properties vary spatially regardless of land uses or management regimes [2][3][4][5]. This was noted in our study with pre-treatment differences for many variables based upon the assigned treatment designations (Tables 1 and 3) and mean ∆ value differences between replications ( Table 2). These inherent differences were one primary reason mean ∆ values were selected for part of the statistical analyses included in this study. In order to accurately assess true, treatment-induced soil changes, pre-existing heterogeneity had to be addressed. Soil order (ultisols), species composition, landform, aspect, elevation, and disturbance history were similar for the forests selected as part of this study [12][13][14]. Despite this similarity, however, some soil chemical properties differed within the Green River Game Lands based upon spatial arrangement, independent of treatment-related effects.
Additionally, treatment effects on the Green River Game Lands were potentially impacted by time since disturbance. One to two years following the first implementation of these treatments, O horizon C:N (p = 0.0472) and mineral soil Ca (p = 0.0270) differed between treatments [13]. O horizon C:N was significantly lower than all other treatments and mineral soil Ca was lower on the cut-only and cut + burn treatments than the untreated control. Two to four years post-treatment, however, these variables did not differ between treatments, but soil pH was lower on the cut-only treatment than the untreated control (p = 0.0457). The 2018 samples were taken 3 years since the last prescribed burns occurred and 6 years since the last cutting operations were conducted. In 2018, mean O horizon (p = 0.0042) and mineral soil N (p = 0.0412) and mineral soil C (p = 0.0074) differed between treatments (Table 1). Of these variables in 2018, mineral soil C and N were the only ones that displayed differences between the untreated control and one of the fuel reduction treatments: the cut + burn treatment.
Inherently, however, the 2018 samples were also influenced not only by time since disturbance, but also treatment frequency. Four prescribed burns and two cutting operations have been conducted since 2001. Repeated treatment implementation is therefore reflected in the 2018 values and the subsequent mean ∆ values. Soils-related treatment effects vary not only as a result of time since treatment, but also as a result of treatment frequency [1,2]. In the case of wildland fire, they may also vary as a result of other fire behavior (i.e. fuels, weather, and topography) or fire regime (i.e. seasonality) factors [14][15][16]. This is particularly noteworthy when considering wildfire and prescribed fire impacts.
Most recently, C and N have received increased research attention as a result of increasing wildfire intensity and frequency in many locations globally. In a recent global meta-analysis, Pelligrini et al. [24] suggested that increased fire frequency, regardless of ignition source, may detrimentally impair soil C and N. These authors suggested that impacts in needleleaf forests may not agree with these trends, however. This would particularly apply to pine forests in the southern United States, an area heavily influenced by prescribed fires as opposed to wildfires. The Green River results and those from other prescribed fire studies (15,25,26) bear consideration in the interpretation of fire-related soil impacts, as well. Soil N may not be a limitation for all managed vegetative systems and structures in all locations globally [16,17,25,26]. Some management objectives, such as wildlife habitat management and hazardous fuel reduction, for example, necessitate different forest structures and species compositions [15,27]. Altered fire regimes, that may include more frequent and persistent prescribed fire ignitions to achieve desired management goals, may directly affect organisms as a primary soil-forming factor [28,29]. Therefore, soil chemical properties may also be altered as a product of management intent [29]. Furthermore, atmospheric N deposition and N saturation in some geographic locations may actually hamper management objectives because anthropogenic N enrichment may stifle any potential competitive advantages pyrophytic vegetation exhibit [15][16][17][18]25,26]. In a recent study of legumes in coastal longleaf pine (Pinus palustris L.) savannas in the southeastern United States, it was determined that leguminous biological N fixation did not adequately account for N losses as a result of repeated prescribed burning [25]. However, the authors found that soil N mineralization remained high and tree productivity was unrelated to N availability. This led the authors to question N limitation in these soils and also consider the impact of anthropogenic N sources that might enhance soil N. When soil results are presented as a result of anthropogenic natural resource management or natural disturbance, one should consider many factors, including but not limited to 1) the soil property in question, 2) soil type, 3) time of sampling (since treatment implementation), 4) treatment intensity, 5) method of soil analysis, 6) sampling depth, 7) primary soil forming factors, 8) spatial variability of sampling, and 9) overall management objectives. Caution should be taken when comparing wildfire and prescribed fire soil effects because fire intensity and severity may vary greatly between these fire types in a given ecosystem [30].
The Green River results offered meaningful comparisons with other prescribed fire studies globally. Where we observed decreased mineral soil C and N as a result of the cut + burn treatment, others have reported no treatment effects on nutrient cycling following cutting and burning in the southern Appalachian Mountains [31]. A lack of change in burn-only mineral soil C:N has been documented by others in oak systems [26,32,33]. In Coastal Plain pine forests, McKee [8] determined that mineral soil Ca increased at the 0-8 cm depth in Alabama (biennial, winter burns) and Florida (periodic and annual winter burns) as a result of long-term prescribed fire use. Similarly, Lavoie et al. [34] reported increased Ca, Mg, and K in mineral soils to the 0-5 cm depth 1-3 years after a single prescribed fire in a longleaf (Pinus palustris L.) and slash pine (Pinus elliottii Engelm.) forest in north central Florida. The authors suggested that increased Ca, Mg, and K in mineral soils might be related to leaching from the organic horizon as they noted decreased levels of these nutrients for the forest floor.
Soil P mean ∆ value differences noted at Green River were significant between the treatments and the largest loss was found in the burn-only treatment, which differed from the cut-only treatment ( Table 2). Similar results were found as a result of burning in an Arkansas tallgrass prairie [35]. In southern Ohio oak forests, P losses were noted 4 years following prescribed burning, but the authors were not certain if this reduction was due to increased plant uptake post-fire or a change in the supply of organic matter and parent material [36]. Phosphorus has a lower relative volatilization temperature (774°C) than other soil nutrients, such as Ca (1484°C) [19]. Burning in Coastal Plain pine forests had little effect on long-term mineral soil P values in several research studies [8,37]. These varying results strengthen the discussion that pre-treatment soil conditions and factors such as initial base saturation, fire intensity, fire residence time, and time since fire can affect potential mineral soil P [5,6]. Mineral soil K, Mg, and pH did not differ between treatments at Green River (Table 2). This agreed with results from other studies in upland oak systems [32] and tallgrass prairies [35].

Management Implications
As noted previously, multiple applications of these treatments have now been completed on the Green River Game Lands [14]. Until 2018, however, no additional soil analyses had been completed as additional treatment installations occurred. Other properties have been evaluated, however, such as fuels, forest structure, and species composition. Waldrop et al. [12] investigated long-term fuel and vegetation dynamics on the Green River Game Lands following the repeated use of these treatments (two cuttings for the cut-only and cut + burn treatments; three burns for the burn-only and cut + burn treatments). Their investigation suggested that the burn-only treatment reduced canopy cover slightly and removed the shrub layer, while the cut + burn treatment produced a stand structure most like that of the desired open woodlands. This was achieved through delayed overstory mortality and a consistent top-killing of most shrub stems as repeated treatment implementations occurred. The burn-only and cut + burn treatments also led to increased oak regeneration and the reduction of most fuel components throughout the study. This fuel reduction in the burn-only and cut + burn treatments would most likely decrease the risk of future wildfire ignitions if the treatments were continued longterm. Additionally, the authors concluded that the burn-only and cut + burn treatments might also eventually produce the desired long-term structural conditions, but many repeated fires and/or fires occurring in different seasons would be necessary. In another investigation of vegetative dynamics at this location (two cuttings for the cut-only and cut + burn treatments; four burns for the burn-only and cut + burn treatments), Oakman et al. [14] determined that structural changes induced by the burn-only and cut + burn treatments resulted in the most favorable conditions for oak and pine regeneration. When comparing the two, the cut + burn treatment was most desired for oak and pine regeneration, however, because the burn-only treatment failed to significantly reduce undesired hardwood competition. Additionally, graminoid density and cover as additional components of the desired, open woodland structure increased most as a result of the cut + burn treatment.
As noted previously in this Discussion, mineral soil C and N results would agree with this assertion: the cut + burn treatment may best promote the desired treatment effects (Table 1). However, this is not uniform for every soil property investigated. For example, mean ∆ values for O horizon C and N were negative for all of the treatments, including the untreated control (i.e. the 2018 values were less than the 2001 values) ( Table 1). Decreases in both were least for the cut-only treatment and most for the burn-only treatment. This may be most attributed to the natural process of combustion and volatilization of organic matter in the burn-only treatment [9]. Prescribed burns in this treatment predominantly consumed surface fuels as opposed to larger coarse woody debris that may have been present in the cut + burn treatment following the cutting practices [38]. Likewise, greater retention of C and N in the cut-only treatment might also be related to the increased fuel loads generated through the cutting practices [39,40]. These changes did not correlate with significant changes in O horizon C: N, however, suggesting that forest floor organic matter quality has remained relatively unchanged since this treatment regime was established.
It is additionally only speculative to suggest correlations between these soils results and the recent vegetation results mentioned above [12,14]. However, perhaps the mineral soil N reduction observed in the cut + burn treatment is related to a desired, increased presence of ectomycorrhizal trees and plants [15]. Without sampling both soils and vegetation during the same sampling periods and locations, this is only hypothesized, but this correlation could be meaningful. Therefore, future research at the Green River Game Lands should include statistical correlations between soils, fuels, and vegetation data collected at the same time. This comprehensive view may more fully inform managers and practitioners about the achievement of long-term restoration and fuel reduction objectives in light of complex ecosystem properties and responses. It should be noted that these results were collected in one geographic region in one predominant soil order. However, this soil order is generally thought to represent approximately 8%-9% of Earth's ice-free land surface, supporting approximately 18% of the world's population [41] and the management practices evaluated have only increased in interest globally since this study began in 2001. Therefore, this longterm evaluation of soil responses to repeated treatments may serve as a valuable source of information for managers, practitioners, and researchers globally.
Additional treatment implementations may be needed to obtain and maintain the open woodland structure, species composition, soil chemistry, and hazardous wildfire fuel reduction desired at this site. A treatment frequency of 3-8 years has been recommended [15]. Fire behavior alterations may be needed to accelerate this process or potentially improve treatment success [12]. Such alterations might include continued dormant season burns at lower, acceptable relative humidity or increased days since precipitation, for example, to potentially increase fire intensity and severity [42]. Growing season burning might also promote enhanced fire behavior, but evidence for this is not conclusive at this time [43,44].

Conclusions
We studied southern Appalachian Mountain forest soils at the Green River Game Lands (North Carolina, USA) to determine if long-term and repeated applications of fuel reduction and ecosystem restoration treatments had substantive effects on forest soil chemistry from 2001-2018. O horizon and mineral soil (0-10 cm depth) C, N, C:N and mineral soil Ca, K, Mg, P, and pH were sampled and assessed. Changes in O horizon C and N and mineral soil C, N, C:N, Ca, and P over this time period differed between treatments. These long-term results did not align with treatment differences detected 1-2 and 2-4 years post-treatment for some variables, highlighting the time-sensitive nature of some soil responses. The change in mineral soil N from 2001-2018 was greatest for the cut + burn treatment and this result may most closely align with one of the restoration goals of promoting more N-limited mineral soils to support ectomycorrhizal, pyrophytic vegetation. However, changes in other soil parameters as a result of the treatments did not fully align with these results. When these soil results were combined with the results of other studies focused on the vegetation and fuels of the Green River Game Lands, it appeared that the continued use of these fuel reduction treatments every 3-8 years may be necessary to fully achieve this study site's desired structure and function. Future fires of higher intensity might also accelerate the achievement of the specific restoration objectives that were targeted at this site and at similar locations throughout the Appalachian Region.