Next Article in Journal
Network Text Analysis of Visitors’ Perception of Multi-Sensory Interactive Experience in Urban Forest Parks in China
Previous Article in Journal
MTL-FFDET: A Multi-Task Learning-Based Model for Forest Fire Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of Slash Management Practices on Water and Nutrient Dynamics in Longleaf Pine Forests

Clemson University Forest Ecology and Fire Science Laboratory, Department of Forestry and Environmental Conservation, College of Agriculture, Forestry and Life Sciences, Clemson University, Clemson, SC 29634, USA
*
Author to whom correspondence should be addressed.
Forests 2022, 13(9), 1449; https://doi.org/10.3390/f13091449
Submission received: 25 July 2022 / Revised: 6 September 2022 / Accepted: 6 September 2022 / Published: 9 September 2022
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
(1) Silvicultural applications that manipulate woody debris loading and the structural composition of a forest can have both short and long-term effects on biogeochemical cycling. Longleaf pine forests have been the historically dominant community types throughout much of the Southeastern United States. Fire exclusion, hardwood encroachment, and resource exploitation have severely reduced the amount of remaining longleaf pine habitats, making ecological restoration necessary. The silvicultural treatments used to reestablish these communities have been widespread, leading to some skepticism regarding the sustainability of certain restoration practices. (2) This study aimed to understand how overstory manipulation and woody debris management affected soil water retention rates and nutrient availability. Using a randomized complete block design, abiotic responses to biomass harvesting, conventional harvesting, and mastication treatments were measured across a soil moisture gradient in the South Carolina sandhills. (3) Our findings indicate that mastication increased soil moisture retention rates by 37% and 41%, on average, compared to conventional harvesting and biomass harvesting, respectively. (4) Additionally, soil nutrient stocks did not decline following any management practice, indicating that both biomass harvesting and mastication treatments may not necessarily impact site productivity in a negative manner. These findings imply that mastication treatments keep moisture retention high and do not immediately change soil nutrient availability in longleaf pine forests. Long-term vegetation response studies should continue to document successional trends in conjunction with moisture retention rates and long-term nutrient pulsing.

1. Introduction

Developed over various landscapes in the Southeastern United States (e.g., flatwoods, sandhills, and clayhills), longleaf pine (Pinus palustris Mill.) forests are often classified based on soil moisture characteristics [1,2,3]. They are one of the most species-rich and diverse forest types outside of the tropics but are also one of the most endangered in North America, represented by only a small fraction of their historical range [4]. Longleaf pine trees thrive in well-drained sandy soils and are both fire-adapted and shade-intolerant, meaning that fire exclusion promotes successional development and inhibits reproduction [5,6]. Historically, these forests have been maintained by frequent (1–3 years) low-intensity fires that exploit the complex herbaceous layer to create a recurrent disturbance regime [7,8,9,10,11,12,13]. As a result, understory vegetation is frequently restored in conjunction with longleaf pines, as it is critical to seedling survival [14,15]. However, plant establishment and restoration are often limited by light availability when there is a densely established overstory, so treatments that decrease the basal area and create canopy gaps will increase light penetration and encourage the growth of early successional grasses and forbs [16].
Many longleaf pine restoration efforts focus on transitioning loblolly pine (Pinus taeda L.), slash pine (Pinus elliottii Engelm), and hardwood stands to longleaf pine forests, using a combination of prescribed burning, herbicide, and mechanical treatments [15,17,18,19,20,21]. Conventional harvesting applications include significantly reducing the basal area, planting longleaf pine seedlings, and employing release treatments, which often consist of 1–2 herbicide applications followed by prescribed burning [15,22,23]. However, establishing herbaceous vegetation in sites such as the South Carolina sandhills can be challenging, as they are often highly eroded and experience water and nutrient deficiencies [24,25,26,27,28]. For example, Van Eerden (1997) found that bunchgrass seedling establishment was lower on xeric sandy sites compared to sites with higher quantities of loam and silt, suggesting that moisture retention was a limiting factor in this region [25]. As a result, restoration methods that increase moisture availability and maintain nutrient stocks could aid in establishing understory vegetation.
With new market demands, alternative silvicultural treatments, such as biomass harvesting and mastication, are sometimes used in restoration, but concerns about their applicability and long-term sustainability have arisen with increased use [29,30,31,32,33,34,35,36,37]. The removal of harvesting residue (i.e., logging slash) and alterations in the spatial distribution of masticated debris influence soil hydrology and nutrient dynamics, but the effects have been widespread [38,39,40,41,42,43,44,45,46,47]. In sandy soils with low cation exchange capacity, soil nutrient capital directly correlates to the amount of soil organic matter present, which is subsequently affected by the quality, quantity, and distribution of input materials [33,36]. Removing tree components with higher nutrient contents than wood (e.g., leaves, cambium, and root tips) could potentially cause declines in long-term site productivity, but could also provide greater germination success [33]. Conversely, the deposition and configuration of masticated wood can result in short-term periods of nitrogen immobilization, directly influencing plant growth and species composition [43,46,47,48]. However, other studies have found that mastication increases inorganic nitrogen and moisture retention rates, which may stimulate the vegetation response [39,40,45,46,47].
By quantifying the abiotic environmental conditions following varying intensities of woody debris manipulation, we hope to identify the optimal restoration methods, while simultaneously setting the stage for future vegetation response studies within the region. Our specific objectives were as follows: (1) to quantify the changes to woody debris loading following different silvicultural treatments (i.e., conventional harvesting, biomass harvesting, mechanical mastication) that result in various spatial configurations and woody debris quantities, (2) to determine how manipulating woody debris during restoration timber harvest affects soil moisture retention rates and nutrient availability along a soil moisture gradient. We hypothesize that (1) the intensified removal of woody debris during the biomass harvesting treatment will result in a parallel decline in nutrient availability, (2) masticated debris beds will yield greater moisture retention rates that will be most apparent in xeric sites during periods of drought and (3) the addition of masticated woody debris will result in a short period of nitrogen immobilization.

2. Materials and Methods

The Hardscramble tract is 305 has and is located in the city of Camden in Kershaw County, SC, USA (coordinates: 34.26732, −80.65668) and exhibits a humid, subtropical climate, with year-round rainfall and long hot summers. Precipitation averages about 1100 mm annually, with average minimum and maximum temperatures ranging from 10.3 °C to 24 °C [49]. The underlying geology of the site is comprised of Cretaceous sediments of the Middendorf Formation, which, in Kershaw county, range from 0 to 107 m in depth and sit on top of crystalline rocks from the Paleozoic age [50]. Thick beds of Cretaceous clays form the first visible layer and are located directly beneath deep layers of highly permeable quartz sand [27]. The property retains fourteen different soil series, all of which classify as either Ultisols or Entisols, and is characterized by various habitats, ranging from xeric mature longleaf pine stands to mesic bottomland hardwood stands. Within the property, the study area is a 28.7 ha subsection in the northeast corner (Figure 1). While mature longleaf pine was present, the majority of the study area was classified as an uneven-aged mixed pine and oak forest, with a loblolly overstory and a densely stratified midstory of hardwood species. Intense shading, caused by a prolonged lack of forest management, inhibited understory development and longleaf pine regeneration for much of the study area.
Pre-treatment data were collected from August 2019 to August 2020, and post-treatment data were collected from May 2021 to September 2021. In Fall 2019, a landscape ecological classification (LEC) system, developed for the Hilly Coastal Plain province of South Carolina, was used to classify sites along a soil moisture gradient, with the following four classifications ranging from the wettest sites to the driest sites: mesic, submesic, subxeric, xeric [3]. This designation was used to create a randomized complete block design that organized the 28.7 ha into four distinct soil moisture blocks. Each soil moisture block was subdivided into three compartments, where we implemented our treatments. Three random sampling points were placed within each treatment using ArcGIS mapping applications (36 total sampling points) (Figure 2a). A seed-tree timber harvest was conducted from November to December of 2020 over the entire study area, using mechanized harvesting equipment (i.e., feller buncher, grapple skidder and knuckle-boom loader). Post-harvest woody debris manipulation simulated the site conditions of the following three different silvicultural treatments: conventional harvesting, biomass harvesting, and fuel mastication. Harvesting slash was scattered across the site for conventional harvesting treatments, removed for biomass harvesting treatments, and masticated for the remaining treatments.
Prior to the implementation of our silvicultural treatments, woody debris measurements and water retention rates were collected at each point using the line intersect method and a series of tensiometers (Soil Measurement Systems, Tucson, AZ, USA) placed at each location [51,52]. Thirty-six 15.24 cm (6 in) tensiometers were used to measure biweekly soil water potential (kPa) from 20 July 2020 to 27 August 2020. Woody debris, with diameters larger than 7.62 cm, were labeled as coarse woody debris (CWD) and volumes were calculated using the conic paraboloid equation (Aa = cross-sectional area at the upper end; Ab = the cross-sectional area at the lower end; L = length) [53].
Conic-paraboloid volume = (L/12) × (5Aa + 5Ab + (AaAb)1/2)
Woody debris, with diameters less than 7.62 cm, were labeled as fine woody debris (FWD) and subdivided into three categories based on diameter size (size class 1 = 0–0.64 cm, size class 2 = 0.64–2.54 cm, and size class 3 = 2.54–7.62 cm). FWD volume was estimated using Huber’s volume formula (Am = cross-sectional area at the midpoint; L = length) [53].
Huber’s volume = L × Am
Estimated per unit values were found using Waddell’s (2002) equations based on DeVries’ (1973) formulas (L= total length of transect line; Vm = volume in cubic meters of individual piece; li = length of individual piece; f = 10,000 m2/ha; SpG = specific gravity of wood; DCR = decay class reduction factor) [51,54,55].
Cubic meters per ha = (π/2 L) × (Vm/li) × f
Kilograms per ha = (cubic meters per ha) × 1000 kg/m3 × SpG × DCR
For FWD biomass estimates, volumes were multiplied by the average bulk density for Quercus spp. (579.7 kg m−3) and an assumed decay class reduction factor of 0.80 [56]. Average oak bulk density was used to avoid underestimation of FWD, which often happens when using average pine bulk densities, especially in mixed stands [56]. Carbon estimates were found by multiplying the per unit biomass weight by a conversion factor (0.491 for hardwoods and 0.521 for softwoods) [54].
Post-treatment, each tensiometer was replaced and measured weekly from May 2021 to September 2021. A 4 m × 20 m strip plot was positioned approximately 2.5 m from the sampling point at the same azimuth as the pre-harvest woody debris transects [38]. Within the post-treatment strip plot, four destructive FWD samples were taken at each plot corner, four soil samples were taken at the midpoints of each plot boundary line, and three bulk density samples were taken at 5 m intervals in the center of each CWD plot (Figure 2b).
The length, large end diameter, and small end diameter for CWD logs were recorded if the length was larger than or equal to 1 m, if the large end diameter was greater than or equal to 15 cm, and if the small end diameter was greater than or equal to 7.62 cm [54,57]. Litter and duff were combined and FWD was separated into diameter size classes. All samples were washed, air-dried for 48-h before being weighed, dried in an air-dry oven at 105 degrees C for 72 h, and finally reweighed. Dry weights were converted to per-unit values (kg ha−1) using sampling estimations and bulk densities provided by the Forest Inventory and Analysis (FIA) program [56]. Per-unit values (kg ha−1) were obtained by scaling up our strip plots, using the combined weight of the woody debris inventoried within the boundaries of each plot.
All the soil samples were sandy and coarse-textured, with a clay layer (B-horizon) at a depth of at least 101.6 cm. The soils were tested for pH, buffer pH, extractable elements (P, K, Ca, Mg, Zn, Mn, Cu, B, Na, S), organic matter, nitrate-nitrogen ( NO 3 N ), and soluble salts. Calculations for the cation exchange capacity (CEC), acidity, and percent base saturation were also included in the analysis. Soil sample concentrations were classified as either sufficient or insufficient based on the Clemson University Agricultural Service Laboratory Soil Test Rating System, which defines soil fertility status based on a nutrient element concentration range [58]. In order to supplement our soil sample findings, weekly soil solution samples were obtained from twelve 30 cm porous cup suction lysimeters (Hanna Instruments HI89300-30) from June 2021 to August 2021. One lysimeter was placed 15.24 cm into the soil in each treatment. Using spectrophotometry, ammonium, nitrate, and phosphate concentrations were measured using the Environmental Protection Agency (EPA) compliant FIA-012 method, the FIA-026 cadmium reduction method (EPA 353.2, SM 4500-NO3), and the FIA-073 sequential flow injection method (EPA 365.1), respectively [59].
All data transformations and analyses were conducted using Microsoft Excel [60] and R [61]. The treatment and moisture block variables were converted into factor data, measures of central tendency and dispersion for our dependent variables were determined, and finally we tested for normal distribution using normal Q-Q plots and the Shapiro–Wilk test on the ANOVA residuals. Levene’s test was used to measure the homogeneity of variances. Data that were not normally distributed were log-transformed to meet the assumptions of normality. In order to better understand the effects of our treatments on soil moisture characteristics and nutrient availability, we used the data compiled from 2020 to 2021 to investigate the different trends. Soil nutrient availability and downed woody debris were analyzed using a two-way analysis of variance (ANOVA), and time-integrated measures (matric water potential) were analyzed with a repeated-measures ANOVA using a Bonferroni correction. Tukey’s HSD post-hoc tests were performed for pairwise comparisons and, for any significant interactions, we developed an interaction plot using least-squares means tests for different factor combinations.

3. Results

3.1. Woody Debris

A preharvest inventory of CWD found that there were no significant differences in the mean amounts of biomass (kg ha1) or the CWD predicted weight of carbon (kg ha1) across the soil moisture blocks or treatments. There was an average of 168.3 kg ha1 (min = 11.6 kg ha1; max = 3303.9 kg ha1) across the entire study area, with only two samples indicating more than 1000 kg ha1. The post-treatment CWD analysis found that there were significant differences in the means between treatments for CWD biomass (p = 0.030) and CWD carbon estimates (p = 0.028) (Figure 3). However, for both woody biomass and carbon, there were no significant findings between soil moisture blocks (p = 0.188; p = 0.284, respectively) and there were no interaction effects between the independent variables (p = 0.056; p = 0.052, respectively). The amount of CWD biomass increased across the treatments from our masticated sites to our conventional harvest sites (Figure 3). The pairwise comparisons found that, on average, conventional harvesting produced 7349 kg ha1 (142%) more CWD and 3645.8 kg ha1 (141%) more CWD carbon compared to our mastication treatment (Figure 3).
The average pre-treatment FWD across all sites was 2802.4 kg ha1 and there were no significant mean differences across the soil moisture blocks (p = 0.080) or across treatments (p = 0.775). Post-treatment destructive sampling found that the mean dry weight for litter and duff was 4046.7 kg ha1 (SD = 2432.6 kg ha1) across the entire study area and that the means across treatments (p = 0.541) and soil moisture blocks (p = 0.535) were not significantly different (Figure 4A). Post-treatment duff and litter depths were not measured. However, using average pre-treatment duff and litter depths, we found that post-treatment bulk density averaged about 16.2 kg m3.
Post-treatment measurements found significant differences between our treatments for size classes 1, 2, and 3 (p = <0.05; p = 0.014; p = 0.009, respectively) (Figure 4). On average, the mastication treatment produced 6642.4 kg ha1 and 6116.8 kg ha1 more woody debris in size class 1 and 7316.9 kg ha1 and 7236.0 kg ha1 in size class 2, compared to the biomass and conventional harvest treatments, respectively (Figure 4B,C). This trend reverses with our largest FWD size class, producing the least amount of woody debris per unit area in masticated treatments (Figure 4D). The conventional harvest treatment produced the largest amount of woody debris in size class 3, and significantly more than the mastication treatment (75% more kg ha1) (Figure 4D). For all moisture blocks and treatments, the harvest generated an average of 5328.5 kg ha1, 11,978.5 kg ha1, and 16,398.5 kg ha1 for our size classes in ascending order.

3.2. Soil Moisture Availability

The pre-harvest measures of soil water matric potential yielded significant differences in the means across dates (p = <0.05), but not across moisture gradients (p = 0.438) or treatments (p = 0.177). Following the implementation of our treatments, we found significant differences between soil moisture blocks (p = <0.05), treatments (p = <0.05), and dates (p = <0.05), indicating that, in addition to sporadic rainfall events, our treatments influenced soil water availability for vegetation (Figure 5). Our statistical analysis produced significant interaction effects between our moisture gradient and our treatments (p = <0.05), as well as between our treatments and recording dates (p = <0.05) (Figure 5). Our masticated treatment remained consistently high in water availability across dates, and the greatest differences between treatments were observed in our xeric sites during dry periods (Figure 5). Soil tension increased along our moisture gradient from mesic to xeric, but did not differ significantly between biomass and conventional harvesting treatments. However, for five sampling periods during the growing season, moisture became a limiting (<−30 kPa) factor or near-limiting factor for the biomass and conventional harvest treatment, while soil water in the masticated treatments never dropped below −20 kPa (Figure 5). For these five sampling dates, soil moisture was significantly higher in the masticated treatments compared to the other treatments (Figure 5).

3.3. Soil Quality

Across all treatments, soil pH was strongly acidic (M = 4.9, SD = 0.2) and did not exhibit a significant difference in the means across treatments (p = 0.607) or between soil moisture blocks (p = 0.338) (Table 1). There was a similar trend for the amount of extractable potassium (K), where low K levels were present across the study area (M = 20.8 ppm, SD = 8.7 ppm) and did not show significant differences between the treatments or soil moisture blocks (p = 0.234; p = 0.593, respectively) (Table 1). Nitrate nitrogen ( NO 3 N ) was low across all sites, with 80.6% of all samples having 0 ppm, 16.7% having 1 ppm, and only one sample showing 5 ppm (Table 1). The only elements that were sufficiently available were manganese (Mn) in 25% of the samples and magnesium (Mg) in 5.6% of the samples (Table 1). Extractable phosphorus (P) was available in insufficient amounts as well (M = 4.3 ppm, SD = 3.3 ppm), exhibiting low amounts of P (0–15 ppm) in 97% of our samples (Table 1). However, there was a significant difference in the means between treatments (p = 0.005) and soil moisture blocks (p = 0.005), with our submesic biomass harvest exhibiting the greatest amount of P per unit area. Pairwise comparisons found that, on average, P levels were 113% greater in biomass harvest treatments compared to our conventional treatments, and 210% higher in submesic sites compared to subxeric sites. However, there were no significant interaction effects (p = 0.248).
Lysimeter readings found that there were significant nitrate ( NO 3 ) concentration differences across soil moisture blocks (p = 0.002), with pairwise comparisons showing higher average concentrations for the subxeric sites. The highest maximum concentration was 6.3 ppm and was located in a subxeric moisture block (Table 2). However, there were no significant differences or interactions between treatments (p = 0.988; p = 0.317, respectively). Additionally, we did not find significant mean differences in the Gaussian peak predictions for our treatments (p = 0.769) or soil moisture blocks (p = 0.063).
Comparatively, mean ammonium ( NH 4 + ) concentrations were low but were significantly different between soil moisture blocks (p = 0.003) and treatments (p = <0.001). On average, our biomass treatment produced significantly lower ammonium concentrations compared to all other treatments, and our mesic site produced significantly higher ammonium concentrations than all other soil moisture blocks (Table 2). We also found significant interaction effects between our independent variables, with conventional harvest sites on mesic soil producing very high concentrations (p = 0.002). However, there were no significant findings for the Gaussian peak predictions across or between the independent variables.   PO 4 3   concentrations showed no significant differences between the means or between the mean Gaussian peak predictions across or between independent variables. Maximum concentrations ranged from 1.34 to 44.94 ppm and the highest predicted mean Gaussian peak prediction was 1.14 ppm, with a majority of the results having less than 0.5 ppm (Table 2).
Soil bulk density (g cm3) increased, on average, as the soil moisture gradient became drier and sandier. Significant differences were found between soil moisture blocks (p = 0.001) but not between treatments (p = 0.900), and there were no significant interaction effects (p = 0.752). Xeric and subxeric soil bulk densities were 15.6% and 13.3% higher, respectively, compared to mesic sites. Soil porosity was lower in the xeric sites, ranging from 43 to 45%, on average, across treatments, suggesting stronger matric potentials. Soil organic matter averaged 1.92% of the overall soil composition by weight across all sites. There were no significant mean differences between treatments (p = 0.166) or soil moisture blocks (p = 0.066). However, there was a significant interaction effect between our independent variables (p = 0.012). Biomass harvesting in the subxeric moisture block exhibited the greatest amount of soil organic matter and significantly higher amounts of organic matter compared to conventional harvesting treatments in the submesic (+1.55%) and subxeric sites (+1.41%), and biomass harvesting treatments in the submesic (+1.87%) and xeric sites (1.61%).

4. Discussion

4.1. Downed Woody Debris

The pre-treatment measurements of downed woody debris found that CWD and FWD loadings were similar across the study area, confirming that our site conditions before the harvest were relatively uniform. Post-treatment FWD and CWD were, on average, different across treatments. Masticated treatments had significantly more woody debris in size classes 1 and 2 compared to the other treatments, which showed no difference. Biomass harvest treatments exhibited the lowest amount of downed woody debris. While masticated treatment should have produced the same amount of woody debris as the conventional treatment, the latter was estimated to produce over 2000 more kg ha1 than mastication. This finding could be a result of our sampling methods, which assumed a consistent shape for the CWD logs. Additionally, we used a plot-based method for measuring all woody debris, which is less consistent than a hybrid methodology [38]. Due to the irregular shape of masticated wood particles, it is recommended that quantification of woody debris should come in the form of a hybrid methodology, in which masticated material is measured using plot-based measurements and larger, coarse woody debris is calculated using the standard planar intercept method [38,62].

4.2. Soil Moisture Availability

Soil water matric potential was consistent across all sites for the pre-treatment findings. Our most significant fluctuations were observed between dates that did not experience precipitation events close to our observation period. Expectedly, the matted debris produced by our mastication treatment had higher and more consistent moisture retention rates even during dry conditions and this was the only treatment where water was never a limiting factor for plant growth. Kreye et al. (2012) found that fuelbeds composed of 0–0.64 cm and 0.64–2.54 cm particles increased water retention, producing response drying times for fuels that ranged between 40 and 69 h [39]. Masticated treatments have not been found to influence drying rates for individual particles when compared to non-fragmented particles of the same size class [39,42]. However, the compacted nature of the bed slows the overall drying rate, potentially resulting in positive vegetation responses and pulses of available nutrients, known as the Assart effect [42,63]. A major concern is that the slower overall drying rates, the decrease in bare soil, and the structural configuration of the woody debris will limit prescribed fire capabilities and nitrogen availability in the short-term [39,40,47,64].
Our differences were magnified in the xeric sites, showing an increase in matric water tension in conventional and biomass treatments, likely due to mineral soil exposure. Loose sandy soils can produce additional heat through reradiation, which subsequently increases the evapotranspiration of surface moisture [65]. Masticated fuelbeds can reduce the rate of surface water evaporation, stabilize soil temperature, and slow the rate of infiltration [42,44,66]. These effects can positively influence vegetation responses, but can also influence when and to what extent certain nutrients become available [44].

4.3. Soil Quality

Overall, soil composition was characteristic of South Carolina sandhill ecosystems, exhibiting high acidity with extremely well-drained sandy soils across the entire site. We found similar findings to Brendemuehl (1967) and Faust (1976), who noted that sandhills’ fall line soils average less than 2% organic matter, with a soil pH range from 4.0 to 5.6 [65,67]. Low pH values can cause deficiencies in soil nutrients, especially nitrogen and phosphorus. Nitrification is highly sensitive to changes in soil pH, litter composition, and soil water content, essentially ceasing at values less than 5.0 pH [68,69,70]. Additionally, inorganic phosphorus is frequently deficient at pH values below 5.5, and these effects are often exacerbated by a lack of organic matter [71].
Soil bulk density varies greatly by soil type and depth, with porous soils containing large amounts of organic matter, silt, or clay having lower bulk densities (0.5–1.6 g cm3) compared to sandy soils with less total pore space (1.3–1.7 g cm3). Deeper sandy soil layers can have even higher rates of compaction (> 2.0 g cm3) [72,73]. For sands and loamy sands, such as those present in the South Carolina sandhills, the USDA-NRCS identifies ideal bulk densities for plant growth as <1.60 g cm3, with bulk densities > 1.80 g cm3 restricting root growth. High bulk densities could, thus, affect long-term vegetation responses in xeric sites by either restricting plant growth or slowing the rate of percolation and enhancing water retention capacities. We found no significant differences between treatments, although bulk density increased in the xeric block. Long-term monitoring should be conducted across the study area to quantify the effects of decomposition on bulk density, especially in masticated sites, where organic material may be more likely to become incorporated into the mineral soil over time.
Due to high soil temperatures, organic matter in this region oxidizes quickly and releases nutrients that are readily leached from the upper soil levels by frequent rainfall events [27,65,74]. Of the 17 essential nutrients and the 6 primary nutrients (C, H, O, N, P, K) required for plant growth, nitrogen (N), sulfur (S), and phosphorus (P) are highly influenced by the presence of organic matter and are often limited in sandhill communities [71,75]. Moreover, following disturbances, subsurface pools of N and P are key indicators for ecosystem recovery and nutrient availability [59]. Longleaf ecosystems in sandhill regions are typically nutrient deficient, and while restorative methods that remove soil carbon (e.g., prescribed burning and biomass harvesting) might result in short-term changes to nutrient availability, they may have little long-term effect on soil quality [76,77]. For example, after only one year following prescribed fire in a longleaf pine stand in the Florida sandhills, C and N pools recovered to 67% and 76% of pre-treatment levels [78].
Inorganic nitrogen ( NO 3 and NH 4 + ) is the most commonly limiting nutrient in longleaf sandhills, due to the low litter quality (high C:N ratio) and high soil acidity. High C:N ratios result in ammonium immobilization, while nitrification is inhibited at low pH values [43,48]. Inorganic phosphorus ( H 2 PO 4 , HPO 4 2 , and PO 4 3 ) is essential for plant growth, energy transfer (ADP and ATP), reproduction, photosynthesis, and several other plant functions, frequently showing deficiencies in soils with pH values below 5.5 or above 6.5 [71]. However, increases in soil organic matter have been found to reduce the strength of P adsorption and the maximum phosphate buffering capacity [79]. Similarly, inorganic sulfur ( SO 4 2 ) is also considerably influenced by soil organic matter and is most important in the formation of chlorophyll. Sulfur is often deficient in coarse sandy soils with good drainage [70,71,80].
We observed that   NO 3 concentrations were considerably higher than   NH 4 + concentrations but were still low (<1 ppm) in most of our soil tests and solutions. This finding was interesting, considering other longleaf restoration studies found that   NH 4 + can supplement longleaf growth during growth-limiting environmental periods that exhibit low   NO 3 concentrations [69]. While we did observe a post-harvest spike of NO 3 in our subxeric sites, there were no discernable trends for nutrient availability across our treatments. Wilson et al. (2002) found that soil temperature increases were the principal influence in producing large pools of inorganic nitrogen in xeric sites. However, this does not explain why there was not a similar pulse of nitrogen compared to our subxeric results [70]. We also noted a parallel spike in soil organic matter in our subxeric site, but it did not seem to influence   NH 4 + or   PO 4 3 . Considering these soils exhibit a low CEC with very little organic matter, nutrients may have leached below the depth of our lysimeters.
Without an observed period of immobilization in the mastication treatment or a heightened deficiency in the biomass treatment, it could be argued that the initial disturbance response to these treatments is negligible. Coates et al. (2010) found similar results, in which fuel reduction treatments had little to no effect on soil properties 2–4 years following their study in the Southern Appalachians [81]. In an already nutrient deficient environment, we found that both biomass harvesting and mastication are essentially equally sustainable with regard to soil quality. Without nutrient deficiencies, biomass harvesting exhibited similar soil characteristics compared to conventional harvesting methods. Conversely, mastication treatments showed similar trends to other studies, which noted that the addition of masticated materials did not result in short-term changes to nutrient availability [47,82,83]. However, the additional and consistent water retention rates produced by that treatment could improve the long-term vegetation responses. Water availability often constrains the establishment and performance of native understory plants in longleaf pine forests [24]. Thus, any treatment that enhances water availability during the growing season could be a viable option for managers seeking to restore the understory community.
Due to the fragmented and compacted nature of masticated woody debris, decomposition typically happens at faster rates than coarse woody debris, decreasing the period of nitrogen immobilization and resulting in faster mineralization [46,47]. In turn, this could potentially influence long-term nitrogen dynamics. For example, Rhoades et al. (2012) noted an initial decline in plant-available nitrogen following mastication, but then observed a rapid rate of mineralization, showing an elevated level of available nitrogen (+32%) three to five years following the treatment [46]. Young et al. (2013) found similar results and concluded that mastication could improve pine seedling growth by increasing the available nitrogen and moisture retention [47].
It is also possible that we are currently experiencing a period of immobilization compared to pre-harvest concentration levels and that, because the C:N ratio would increase similarly across all treatments just in different spatial arrangements, our initial findings do not exhibit any differences between treatments. This could lead to a delayed mineralization process that manifests itself several years following the initiation of our treatments [46]. Pre-treatment concentrations were not evaluated, but long-term studies on mineralization could still indicate significant fluctuations in soil organic matter and nutrient availability. For example, nitrogen limitations could constrain the initial establishment of pioneer vegetation and influence successional dynamics as these plant communities develop. Long-term studies that monitor changes in soil chemistry and vegetation responses could identify successional trends directly correlated to nutrient stock fluctuations.

5. Conclusions

Our study indicates that the intense removal or addition of biomass caused by our treatments did not measurably impact the nutrient pools. These findings suggest that both biomass harvesting and mastication could be utilized as sustainable alternatives for longleaf pine restoration in the region, compared to conventional harvesting methods. Mastication reduced woody debris loading for larger size classes, while simultaneously retaining soil moisture at significantly higher rates compared to other treatments. By exhibiting consistent soil moisture retention rates, even during the driest parts of the growing season, masticated fuelbeds could potentially result in faster vegetation responses and/or a different successional trajectory compared to other treatments. Over time, this will likely result in more rapid nutrient mineralization rates and organic matter incorporation. However, masticated fuelbeds often smolder when prescribed fire is employed, meaning that if there is not a positive vegetation response (e.g., the establishment of fine live fuels), fire implementation could be limited. Without a robust herbaceous understory to carry low-intensity surface fires, masticated sites could experience encroachment by fire-sensitive species. To quantify these changes and understand the influence of silvicultural treatments with regard to restoration, additional management applications and long-term studies need to be applied to comprehend the ecological trends in the region.

Author Contributions

Conceptualization, J.M., D.H. and P.H.; methodology, P.H.; software, P.H.; validation, D.H., P.H. and R.B.; formal analysis, J.M., P.H. and D.H.; investigation, J.M.; resources, P.H., D.H. and R.B.; data curation, J.M.; writing—original draft preparation, J.M.; writing—review and editing, J.M. and D.H.; visualization, J.M.; supervision, D.H.; project administration, J.M.; funding acquisition, R.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Margaret H. Lloyd-SmartState Endowment through Clemson University and received no external funding.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to acknowledge members of the Margaret H. Lloyd Foundation and Clemson’s Department of Forestry and Environmental Conservation, including numerous staff members, faculty, and students who provided insight and direction during the research and publication process. This includes, but is not limited to, Wayne Chao, Phil Gaines, Leoncia Cruz, and Ellie Johnson. Additionally, Coy Myers and Forest Land Management Inc. were instrumental in our logistical planning process, particularly as it related to scheduling silvicultural applications. We also acknowledge the contributions of Jesse Kreye, who provided guidance on our woody debris sampling methodologies.

Conflicts of Interest

The funders were involved in the design of the study and were invested in the improvement of the study site. Clemson University received both the Hardscramble property and a monetary endowment, which was used to fund this research.

References

  1. Abrahamson, W.G.; Hartnett, D.C. Pine flatwoods and dry prairies. In Ecosystems of Florida; University of Central Florida Press: Orlando, FL, USA, 1990; pp. 103–149. [Google Scholar]
  2. Myers, R.L. Scrub and high pine. In Ecosystems of Florida; Myers, R.L., Ewel, J.J., Eds.; The University of Central Florida Press: Gainesveille, FL, USA, 1990; pp. 150–193. [Google Scholar]
  3. Van Lear, D.H.; Jones, S.M. An Example of Site Classification in the Southeastern Coastal Plain Based on Vegetation and Land Type¹. South. J. Appl. For. 1987, 11, 23–28. [Google Scholar] [CrossRef]
  4. Peet, R.K.; Allard, D.J. Longleaf pine vegetation of the southern Atlantic and eastern Gulf Coast regions: A preliminary clas-sification. In Proceedings of the Tall Timbers Fire Ecology Conference, Tallahassee, FL, USA, 30 May 1993; Volume 18, pp. 45–81. [Google Scholar]
  5. Frost, C.C. Four centuries of changing landscape patterns in the longleaf pine ecosystem. In Proceedings of the Tall Timbers Fire Ecology Conference, Tallahassee, FL, USA, 30 May 1993; Volume 18, pp. 17–43. [Google Scholar]
  6. Longleaf Pine Ecosystem Fact Sheet–U.S. Fish and Wildlife Service 2003. Available online: http://www.southeast.fws.gov/pfwpine.html (accessed on 5 November 2021).
  7. Christensen, N.L. Fire regimes in southeastern ecosystems. Fire Regimes Ecosyst. Prop. 1981, 112, 136. [Google Scholar]
  8. Cox, A.; Gordon, D.; Slapcinsky, J.; Seamon, G. Understory restoration in longleaf pine Sandhills. Nat. Areas J. 2004, 24, 4–14. [Google Scholar]
  9. Engstrom, R.T. Characteristic mammals and birds of longleaf pine forests. In Proceedings of the Tall Timbers Fire Ecology Conference, Tallahassee, FL, USA, 30 May 1993; Volume 18, pp. 127–138. [Google Scholar]
  10. Stambaugh, M.C.; Guyette, R.P.; Marschall, J.M. Longleaf pine (Pinus palustris Mill.) fire scars reveal new details of a frequent fire regime. J. Veg. Sci. 2011, 22, 1094–1104. [Google Scholar] [CrossRef]
  11. Reinhart, K.O.; Menges, E.S. Effects of re-introducing fire to a central Florida sandhill community. Appl. Veg. Sci. 2004, 7, 141–150. [Google Scholar] [CrossRef]
  12. Wade, D.D.; Lundsford, J. Fire as a forest management tool: Prescribed burning in the southern United States. Unasylva 1990, 41, 28–38. [Google Scholar]
  13. Walker, J.; Peet, R. Composition and species diversity of pine-wiregrass savannas of the Green Swamp, North Carolina. Plant Ecol. 1984, 55, 163–179. [Google Scholar] [CrossRef]
  14. Walker, J.L.; Silletti, A.M. Restoring the ground layer of longleaf pine ecosystems. In The Longleaf Pine Ecosystem; Jose, S., Jokela, E.J., Miller, D.L., Eds.; Springer: New York, NY, USA, 2007; pp. 297–333. [Google Scholar] [CrossRef]
  15. Brockway, D.G. Restoration of Longleaf Pine Ecosystems; USDA Forest Service, Southern Research Station: Gainesville, FL, USA, 2005; Volume 83. [Google Scholar]
  16. Sharma, A.; Jose, S.; Bohn, K.K.; Andreu, M.G. Effects of reproduction methods and overstory species composition on un-derstory light availability in longleaf pine–slash pine ecosystems. For. Ecol. Manag. 2012, 284, 23–33. [Google Scholar] [CrossRef]
  17. Brockway, D.G.; Outcalt, K.W.; Estes, B.L.; Rummer, R. Vegetation response to midstorey mulching and prescribed burning for wildfire hazard reduction and longleaf pine (Pinus palustris Mill.) ecosystem restoration. For. Int. J. For. Res. 2009, 82, 299–314. [Google Scholar] [CrossRef]
  18. Outcalt, K.W. Factors affecting wiregrass (Aristida stricta Michx.) cover on uncut and site prepared sandhills areas in Central Florida. Ecol. Eng. 1992, 1, 245–251. [Google Scholar] [CrossRef]
  19. Outcalt, K.W.; Lewis, C.E. Response of wiregrass (Aristida stricta) to mechanical site preparation. In Wiregrass Biology and Management Symposium Proceedings; KBN Engineering and Applied Science: Gainesville, FL, USA, 1990. [Google Scholar]
  20. Provencher, L.; Herring, B.J.; Gordon, D.R.; Rodgers, H.L.; Galley, K.E.M.; Tanner, G.W.; Hardesty, J.L.; Brennan, L.A. Effects of Hardwood Reduction Techniques on Longleaf Pine Sandhill Vegetation in Northwest Florida. Restor. Ecol. 2001, 9, 13–27. [Google Scholar] [CrossRef]
  21. Walker, J.L.; van Eerden, B.P.; Robinson, D.; Hausch, M. Burning and chopping for woodpeckers and wiregrass? In Red-Cockaded Woodpecker: Road to Recovery; Costa, R., Daniels, S.J., Eds.; Hancock House Publishers: Blaine, WA, USA, 2004; pp. 683–686. [Google Scholar]
  22. Knapp, B.O.; Wang, G.G.; Walker, J.L.; Cohen, S. Effects of site preparation treatments on early growth and survival of planted longleaf pine (Pinus palustris Mill.) seedlings in North Carolina. For. Ecol. Manag. 2006, 226, 122–128. [Google Scholar] [CrossRef]
  23. Knapp, B.O.; Walker, J.L.; Wang, G.G.; Hu, H.; Addington, R.N. Effects of overstory retention, herbicides, and fertilization on sub-canopy vegetation structure and functional group composition in loblolly pine forests restored to longleaf pine. For. Ecol. Manag. 2014, 320, 149–160. [Google Scholar] [CrossRef]
  24. Hagan, D.L.; Jose, S.; Thetford, M.; Bohn, K. Production physiology of three native shrubs intercropped in a young longleaf pine plantation. Agrofor. Syst. 2009, 76, 283–294. [Google Scholar] [CrossRef]
  25. Van Eerden, B.P. Studies on the Reproductive Biology of Wiregrass (Aristida stricta Michaux) in the Carolina Sandhills. Doc-toral dissertation, University of Georgia, Athens, GA, USA, 1997. [Google Scholar]
  26. MLRA 137: Carolina and Georgia Sand Hills. Available online: https://www.nrcs.usda.gov/wps/portal/nrcs/detail/ga/soils/surveys/?cid=nrcs144p2_021884 (accessed on 21 October 2021).
  27. South Carolina Encyclopedia—Sandhills. Available online: https://www.scencyclopedia.org/sce/entries/Sandhills/ (accessed on 8 July 2021).
  28. USDA-NRCS. Land resource regions and major land resource areas of the United States, the Caribbean, and the Pacific Ba-sin. In US Department of Agriculture Handbook; USDA: Washington, DC, USA, 2006; p. 296. [Google Scholar]
  29. Aguilar, F.X.; Mirzaee, A.; McGarvey, R.G.; Shifley, S.R.; Burtraw, D. Expansion of US wood pellet industry points to positive trends but the need for continued monitoring. Sci. Rep. 2020, 10, 1–17. [Google Scholar] [CrossRef]
  30. Costanza, J.K.; Abt, R.C.; McKerrow, A.J.; Collazo, J.A. Linking state-and-transition simulation and timber supply models for forest biomass production scenarios. AIMS Environ. Sci. 2015, 2, 180–202. [Google Scholar] [CrossRef]
  31. North, B.W.; Pienaar, E.F. Continued obstacles to wood-based biomass production in the southeastern United States. GCB Bioenergy 2021, 13, 1043–1053. [Google Scholar] [CrossRef]
  32. U.S. Department of Energy. 2016 Billion-Ton Report: Advancing Domestic Resources for a Thriving Bioeconomy, Volume 2: Environmental Sustainability Effects of Select Scenarios from Volume 1; ORNL/TM-2016/727; Efroymson, R.A., Langholtz, M.H., Johnson, K.E., Stokes, B.J., Eds.; Oak Ridge National Laboratory: Oak Ridge, TN, USA, 2017; p. 640. [Google Scholar] [CrossRef]
  33. Janowiak, M.K.; Webster, C.R. Promoting ecological sustainability in woody biomass harvesting. J. For. 2010, 108, 16–23. [Google Scholar] [CrossRef]
  34. Richardson, J.; Björheden, R.; Hakkila, P.; Lowe, A.T.; Smith, C.T. Bioenergy from Sustainable Forestry: Guiding Principles and Practice; Springer Science & Business Media: Berlin, Germany, 2006; Volume 71. [Google Scholar]
  35. Vance, E.D.; Aust, W.M.; Strahm, B.D.; Froese, R.E.; Harrison, R.B.; Morris, L.A. Biomass Harvesting and Soil Productivity: Is the Science Meeting our Policy Needs? Soil Sci. Soc. Am. J. 2014, 78, S95–S104. [Google Scholar] [CrossRef]
  36. Walker, T.; Cardellichio, P.; Colnes, A.; Gunn, J.; Kittler, B.; Perschel, B.; Initiative, N.C. Biomass Sustainability and Carbon Policy Study; Manomet Center for Conservation Sciences: Plymouth, MA, USA, 2010. [Google Scholar]
  37. Wang, J.; LeDoux, C.B.; Edwards, P.; Jones, M. Soil bulk density changes caused by mechanized harvesting: A case study in central Appalachia. For. Prod. J. 2005, 55, 37–40. [Google Scholar]
  38. Kane, J.M.; Varner, J.M.; Knapp, E.E. Novel fuelbed characteristics associated with mechanical mastication treatments in northern California and south-western Oregon, USA. Int. J. Wildland Fire 2009, 18, 686–697. [Google Scholar] [CrossRef]
  39. Kreye, J.K.; Varner, J.M.; Knapp, E.E. Moisture desorption in mechanically masticated fuels: Effects of particle fracturing and fuelbed compaction. Int. J. Wildland Fire 2012, 21, 894–904. [Google Scholar] [CrossRef]
  40. Kreye, J.K.; Kobziar, L.N.; Camp, J.M. Immediate and short-term response of understory fuels following mechanical mastication in a pine flatwoods site of Florida, USA. For. Ecol. Manag. 2014, 313, 340–354. [Google Scholar] [CrossRef]
  41. Hood, S.; Wu, R. Estimating fuel bed loadings in masticated areas. In Proceedings of the Fuels Management-How to Measure Success Conference, Portland, OR, USA, 28–30 March 2006; Andrews, P.L., Butler, B.W., Eds.; US Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2006; Volume 41, pp. 333–340. [Google Scholar]
  42. Kreye, J.; Varner, J.M. Moisture dynamics in masticated fuelbeds: A preliminary analysis. In Proceedings of the fire environment--innovations, management, and policy conference, Destin, FL, USA, 26–30 March 2007; Butler, B.W., Cook, W., Eds.; US Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2007; Volume 46, pp. 173–186. [Google Scholar]
  43. Overby, S.; Gottfries, G. Does mastication residue alter soil nitrogen dynamics in woodlands of Southwest Colorado? In Proceedings of the North American Forest Ecology Workshop conference, Logan, UT, USA, 22–26 June 2009. [Google Scholar]
  44. Jain, T.; Sikkink, P.; Keefe, R.; Byrne, J. To masticate or not: Useful tips for treating forest, woodland, and shrubland vegetation. Gen. Technol. Rep. 2018, 381, 5–32. [Google Scholar] [CrossRef]
  45. Kreye, J.K.; Varner, J.M.; Knapp, E.E. Effects of particle fracturing and moisture content on fire behaviour in masticated fuelbeds burned in a laboratory. Int. J. Wildland Fire 2011, 20, 308–317. [Google Scholar] [CrossRef]
  46. Rhoades, C.; Battaglia, M.; Rocca, M.; Ryan, M. Short- and medium-term effects of fuel reduction mulch treatments on soil nitrogen availability in Colorado conifer forests. For. Ecol. Manag. 2012, 276, 231–238. [Google Scholar] [CrossRef]
  47. Young, K.; Roundy, B.A.; Eggett, D.L. Plant Establishment in Masticated Utah Juniper Woodlands. Rangel. Ecol. Manag. 2013, 66, 597–607. [Google Scholar] [CrossRef]
  48. Sahrawat, K.L. Factors affecting nitrification in soils. Commun. Soil Sci. Plant Anal. 2008, 39, 1436–1446. [Google Scholar] [CrossRef]
  49. Kershaw County Climate—South Carolina State Climatology Office 2022. Available online: https://www.dnr.sc.gov/climate/sco/ClimateData/countyData/county_kershaw.php (accessed on 29 June 2022).
  50. Benjamin, J.G.; Seymour, R.S.; Meacham, E.; Wilson, J. Impact of whole-tree and cut-to-length harvesting on postharvest condition and logging costs for early commercial thinning in maine. North. J. Appl. For. 2013, 30, 149–155. [Google Scholar] [CrossRef]
  51. Newcome, R. Ground-Water Resources of Kershaw County, South Carolina; South Carolina State Documents Depository: Columbia, SC, USA, 2002. [Google Scholar]
  52. Briedis, J.I.; Wilson, J.S.; Benjamin, J.G.; Wagner, R.G. Logging residue volumes and characteristics following integrated roundwood and energy-wood whole-tree harvesting in central maine. North. J. Appl. For. 2011, 28, 66–71. [Google Scholar] [CrossRef]
  53. Fraver, S.; Ringvall, A.; Jonsson, B.G. Refining volume estimates of down woody debris. Can. J. For. Res. 2007, 37, 627–633. [Google Scholar] [CrossRef]
  54. Waddell, K.L. Sampling coarse woody debris for multiple attributes in extensive resource inventories. Ecol. Indic. 2002, 1, 139–153. [Google Scholar] [CrossRef]
  55. De Vries, P.G. A General Theory on Line Intersect Sampling: With Application to Logging Residue Inventory; Forest Mensura-tion Department, Agricultural University: Wageningen, NL, USA, 1973; No. 73-11. [Google Scholar]
  56. Woodall, C.W.; Monleon, V.J. Sampling, Estimation, and Analysis Procedures for the Down Woody Materials Indicator; USDA Forest Service General Technical Report NRS-22 Northern Research Station: Newtown Square, PA, USA, 2008. [Google Scholar]
  57. Waldrop, T.; Phillips, R.A.; Simon, D.A. Fuels and predicted fire behavior in the southern Appalachian Mountains after fire and fire surrogate treatments. For. Sci. 2010, 56, 32–45. [Google Scholar] [CrossRef]
  58. Agricultural Service Laboratory Soil Testing Results—Clemson University Agricultural Service Laboratory. 2014. Available online: http://psaweb.clemson.edu/soils/htdocs/understandingreport.pdf (accessed on 12 October 2021).
  59. Klimas, K.B. Prescribed Fire Effects on Water Quality Variables in the Southern Appalachian Region. Doctoral Dissertation, Clemson University, Clemson, SC, USA, 2020. [Google Scholar]
  60. Microsoft Corporation. Microsoft Excel. 2018. Available online: https://office.microsoft.com/excel (accessed on 11 June 2021).
  61. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: https://www.R-project.org/ (accessed on 22 September 2021).
  62. Brown, J.K. Handbook for Inventorying Downed Woody Material; US Department of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station: Ogden, UT, USA, 1974; p. 16. [Google Scholar]
  63. Wiebe, S.; Morris, D.; Luckai, N.; Reid, D. Coarse Woody Debris Dynamics Following Biomass Harvesting: Tracking Carbon and Nitrogen Patterns During Early Stand Development in Upland Black Spruce Ecosystems. Int. J. For. Eng. 2012, 23, 25–32. [Google Scholar] [CrossRef]
  64. Rhoades, C.; Battaglia, M.; Ryan, M.G.; Rocca, M.E. Woody Mulch Effects on Soil Nitrogen Availability in Mechanical Fuel Reduction Treatments. In Proceedings of the North American Forest Ecology Workshop conference, Logan, UT, USA, 22–26 June 2009. [Google Scholar]
  65. Faust, W.Z. A vegetation analysis of the Georgia fall-line Sandhills. Rhodora 1976, 78, 525–531. [Google Scholar]
  66. Qin, S.; Li, S.; Kang, S.; Du, T.; Tong, L.; Ding, R. Can the drip irrigation under film mulch reduce crop evapotranspiration and save water under the sufficient irrigation condition? Agric. Water Manag. 2016, 177, 128–137. [Google Scholar] [CrossRef]
  67. Brendemuehl, R.H. Research progress in the use of fertilizers to increase pine growth on the Florida Sandhills. For. Fertil. -Zation-Theory Pract. 1968, 1, 191–196. [Google Scholar]
  68. Linn, D.M.; Doran, J.W. Effect of Water-Filled Pore Space on Carbon Dioxide and Nitrous Oxide Production in Tilled and Nontilled Soils. Soil Sci. Soc. Am. J. 1984, 48, 1267–1272. [Google Scholar] [CrossRef]
  69. McCaskill, G.L.; Jose, S.; Chauhan, A.; Ogram, A.V. Soil nitrogen dynamics as an indicator for longleaf pine restoration. Restor. Ecol. 2017, 26, 264–274. [Google Scholar] [CrossRef]
  70. Wilson, C.A.; Mitchell, R.J.; Boring, L.R.; Hendricks, J.J. Soil nitrogen dynamics in a fire-maintained forest ecosystem: Results over a 3-year burn interval. Soil Biol. Biochem. 2002, 34, 679–689. [Google Scholar] [CrossRef]
  71. Nutrients Plants Require for Growth—University of Idaho College of Agricultural and Life Sciences 2004. Available online: https://www.extension.uidaho.edu/publishing/pdf/CIS/CIS1124.pdf (accessed on 17 August 2021).
  72. McKenzie, N.; Coughlan, K.; Cresswell, H. Soil Physical Measurement and Interpretation for Land Evaluation; Csiro Publishing: Clayton, Australia, 2002; Volume 5. [Google Scholar]
  73. Bulk Density Measurement—Fact Sheets. Available online: https://soilquality.org.au/factsheets/bulk-density-measurement (accessed on 14 September 2021).
  74. Burns, R.M.; Hebb, E.A. Site Preparation and Reforestation of Droughty, Acid Sands; US Department of Agriculture, Forest Ser-vice: Washington, DC, USA, 1972. [Google Scholar]
  75. Donovan, L.A.; West, J.B.; McLeod, K.W. Quercus species differ in water and nutrient characteristics in a resource-limited fall-line sandhill habitat. Tree Physiol. 2000, 20, 929–936. [Google Scholar] [CrossRef] [PubMed]
  76. Lavoie, M.; Mack, M.C.; Hiers, J.K.; Pokswinski, S.; Barnett, A.; Provencher, L. Effects of Restoration Techniques on Soil Carbon and Nitrogen Dynamics in Florida Longleaf Pine (Pinus palustris) Sandhill Forests. Forests 2014, 5, 498–517. [Google Scholar] [CrossRef]
  77. Raison, R.J.; Khanna, P.K.; Jacobsen, K.L. Effects of Fire on Forest Nutrient Cycles Joan Romanya and Isabel Serrasolses. In Fire Effects on Soils and Restoration Strategies, 1st ed.; CRC Press: Boca Raton, FL, USA, 2009; pp. 241–272. [Google Scholar]
  78. Lavoie, M.; Starr, G.; Mack, M.C.; Martin, T.A.; Gholz, H.L. Effects of a Prescribed Fire on Understory Vegetation, Carbon Pools, and Soil Nutrients in a Longleaf Pine-Slash Pine Forest in Florida. Nat. Areas J. 2010, 30, 82–94. [Google Scholar] [CrossRef]
  79. Yang, X.; Chen, X.; Yang, X. Effect of organic matter on phosphorus adsorption and desorption in a black soil from Northeast China. Soil Tillage Res. 2018, 187, 85–91. [Google Scholar] [CrossRef]
  80. Sulfur—The 4th Major Nutrient. 2010. Available online: www.ipni.net/pnt (accessed on 27 September 2021).
  81. Coates, T.A.; Shelburne, V.B.; Waldrop, T.A.; Smith, B.R.; Hoke, S.; Hill, J.; Simon, D.M. Forest soil response to fuel reduction treatments in the southern Appalachian Mountains. In Proceedings of the 14th Biennial Southern Silvicultural Research Conference, Asheville, NC, USA, 26 February–1 March 2010; Stanturf, J.A., Ed.; US Department of Agriculture, Forest Service, Southern Research Station: Asheville, NC, USA, 2010; Volume 121, pp. 283–287. [Google Scholar]
  82. Kobziar, L.N.; Long, A.J.; Zipperer, W.C.; Kreye, J.K. Characterization of Masticated Fuelbeds and Fuel Treatment Effectiveness in Southeastern US Pine Ecosystems. Jt. Fire Sci. Program, UNL Digital Commons 2014. Available online: https://digitalcommons.unl.edu/jfspresearch/35/ (accessed on 4 September 2021).
  83. Overby, S.T.; Gottfried, G.J. Microbial and nitrogen pool response to fuel treatments in Pinyon-Juniper woodlands of the southwestern USA. For. Ecol. Manag. 2017, 406, 138–146. [Google Scholar] [CrossRef]
Figure 1. Map of Kershaw county, SC, USA (A), the Hardscramble location imbedded in Kershaw county (B), and the 28.7 ha longleaf pine restoration study site in Camden, SC, USA (C).
Figure 1. Map of Kershaw county, SC, USA (A), the Hardscramble location imbedded in Kershaw county (B), and the 28.7 ha longleaf pine restoration study site in Camden, SC, USA (C).
Forests 13 01449 g001
Figure 2. (a) Map of the longleaf pine restoration study site in Camden, SC, USA, including soil moisture blocks, treatments, and sampling points; (b) post-treatment sampling scheme, which includes plot-based sampling points for FWD, CWD, bulk density, soil samples, and the tensiometer location.
Figure 2. (a) Map of the longleaf pine restoration study site in Camden, SC, USA, including soil moisture blocks, treatments, and sampling points; (b) post-treatment sampling scheme, which includes plot-based sampling points for FWD, CWD, bulk density, soil samples, and the tensiometer location.
Forests 13 01449 g002
Figure 3. Bar plot (mean ± SD) of CWD (kg ha1) comparisons between pre- and post-treatment findings for biomass, mastication, and conventional treatments in the longleaf restoration site in Camden, SC, USA. Significant mean differences before and after treatment implementation are illustrated by lowercase and capital letters (α < 0.05). Different capitalized letters indicate statistically significant differences between treatments (α < 0.05).
Figure 3. Bar plot (mean ± SD) of CWD (kg ha1) comparisons between pre- and post-treatment findings for biomass, mastication, and conventional treatments in the longleaf restoration site in Camden, SC, USA. Significant mean differences before and after treatment implementation are illustrated by lowercase and capital letters (α < 0.05). Different capitalized letters indicate statistically significant differences between treatments (α < 0.05).
Forests 13 01449 g003
Figure 4. Post-treatment bar plots with standard error bars for mean weights (kg ha1) of litter and duff (A), woody debris in size class 1 (0–0.64 cm) (B), woody debris in size class 2 (0.64–2.54 cm) (C), and woody debris in size class 3 (2.54–7.62 cm) (D) across soil moisture blocks and treatments. Statistically different means are represented by an asterisk located next to the legend for between-treatment differences (α < 0.05).
Figure 4. Post-treatment bar plots with standard error bars for mean weights (kg ha1) of litter and duff (A), woody debris in size class 1 (0–0.64 cm) (B), woody debris in size class 2 (0.64–2.54 cm) (C), and woody debris in size class 3 (2.54–7.62 cm) (D) across soil moisture blocks and treatments. Statistically different means are represented by an asterisk located next to the legend for between-treatment differences (α < 0.05).
Forests 13 01449 g004
Figure 5. Average soil matric water potential (kPa) readings between treatments for all moisture blocks taken weekly from 24 May 2021 to 15 September 2021 in the Hardscramble longleaf pine restoration site in Camden, SC. Dates marked by an asterisk have significantly lower soil water tension in the mastication treatment compared to both other treatments (α < 0.05).
Figure 5. Average soil matric water potential (kPa) readings between treatments for all moisture blocks taken weekly from 24 May 2021 to 15 September 2021 in the Hardscramble longleaf pine restoration site in Camden, SC. Dates marked by an asterisk have significantly lower soil water tension in the mastication treatment compared to both other treatments (α < 0.05).
Forests 13 01449 g005
Table 1. Soil composition derived from soil samples taken from the Hardscramble property longleaf pine restoration site in Camden, SC, USA. For each soil moisture block and treatment, soil pH, extractable elements (P, K, Ca, Mg, Zn, Mn, Cu, B, Na, S), cation exchange capacity (CEC), nitrate-nitrogen ( NO 3 N ), and organic matter were measured.
Table 1. Soil composition derived from soil samples taken from the Hardscramble property longleaf pine restoration site in Camden, SC, USA. For each soil moisture block and treatment, soil pH, extractable elements (P, K, Ca, Mg, Zn, Mn, Cu, B, Na, S), cation exchange capacity (CEC), nitrate-nitrogen ( NO 3 N ), and organic matter were measured.
Moisture BlockTreatmentSoil pHP (kg/ha)K (kg/ha)Ca (kg/ha)Mg (kg/ha)Zn (kg/ha)Mn (kg/ha)Cu (kg/ha)B (kg/ha)Na (kg/ha)S (kg/ha)CEC (mmhos/cm) NO 3 N   ppm OM %
MesicConventional511.228.089.712.31.213.50.20.25.67.80.0301.12
4.93.443.7100.928.00.711.20.10.36.712.30.0301.83
53.461.6244.335.91.522.40.20.311.211.20.0512.28
Biomass515.726.996.412.30.710.10.20.24.57.80.0311.96
4.93.430.368.415.70.87.80.10.26.713.50.0501.39
5.27.896.4358.753.81.942.60.40.410.112.30.0502.98
Mastication5.24.542.6107.619.10.415.70.20.36.79.00.0412.17
5.26.740.497.523.50.944.80.30.39.09.00.0301.74
4.513.542.6154.726.91.69.00.20.25.67.80.0502.84
SubmesicConventional513.519.1103.110.10.815.70.30.24.59.00.0301.26
4.717.944.8130.023.51.319.10.30.37.89.00.0601.52
4.95.637.060.524.70.65.60.30.39.038.10.0401.78
Biomass4.611.226.970.612.30.75.60.20.27.86.70.0501.21
4.839.237.049.310.10.86.70.20.25.66.70.0501.19
4.820.229.194.214.60.814.60.20.26.711.20.0401.2
Mastication5.15.657.2191.729.11.610.10.10.35.612.30.0402.01
511.267.3233.133.62.526.90.30.39.010.10.0501.79
4.84.542.6178.226.91.224.70.10.24.59.00.0312.38
SubxericConventional4.82.216.834.77.80.75.60.20.25.614.60.0301.66
4.73.423.557.216.80.87.80.20.413.514.60.0401.68
4.93.441.5146.831.41.115.70.30.313.514.60.0401.63
Biomass4.65.691.9520.178.53.052.70.30.417.911.20.0703.89
4.85.681.8378.874.02.646.00.30.417.913.50.0702.92
4.75.644.8134.519.11.328.00.30.36.714.60.0952.39
Mastication5.35.650.4274.648.21.032.50.20.310.111.20.0402.17
55.647.1182.732.51.542.60.20.311.210.10.0301.40
4.64.566.1130.026.90.913.50.20.37.89.00.0502.15
XericConventional4.93.449.368.419.10.66.70.20.411.215.70.0511.35
4.75.660.5171.531.41.821.30.20.35.611.20.0502.58
4.83.443.7108.725.81.012.30.20.26.79.00.0401.88
Biomass4.922.458.3288.147.11.726.90.30.35.69.00.0401.7
517.933.695.316.80.931.40.40.36.79.00.0301.32
5.17.830.3156.922.40.926.90.20.35.612.30.0301.34
Mastication4.64.532.580.714.60.914.60.30.27.810.10.0502.5
4.812.348.2193.925.81.151.60.30.311.214.60.0401.92
5.17.881.8404.678.52.058.30.30.410.110.10.0612.05
Table 2. Mean parameters of soil solution variable nutrient concentrations (ppm) at a 15.24 cm soil depth between woody debris manipulation treatments and soil moisture blocks in Camden, SC, USA. Parenthetical values are standard deviation.
Table 2. Mean parameters of soil solution variable nutrient concentrations (ppm) at a 15.24 cm soil depth between woody debris manipulation treatments and soil moisture blocks in Camden, SC, USA. Parenthetical values are standard deviation.
Moisture GroupTreatment N O 3 N H 4 + P O 4 3
Predicted PeakAverageMaxPredicted PeakAverageMaxPredicted PeakAverageMax
MesicMastication0.56 (0.03)0.43 (0.17)0.580.21 (0.34)0.17 (0.24)0.610.24 (0.16)1.11 (1.84)6.35
Biomass0.56 (0.01)0.52 (0.10)0.620.03 (0.02)0.04 (0.06)0.191.14 (1.57)5.51 (13.40)44.94
Conventional0.74 (0.05)0.65 (0.12)0.790.19 (0.18)0.29 (0.18)0.660.23 (0.04)0.57 (1.06)2.83
SubmesicMastication0.91 (0.35)0.87 (0.40)1.270.04 (0.02)0.05 (0.03)0.130.63 (0.75)3.68 (9.75)33.55
Biomass0.67 (0.19)0.59 (0.18)0.890.01 (0.01)0.01 (0.04)0.140.38 (0.35)1.20 (2.84)10.01
Conventional0.65 (0.13)0.52 (0.18)0.80.01 (0.01)0.03 (0.05)0.160.17 (0.06)0.33 (0.36)1.34
SubxericMastication3.35 (2.88)3.33 (2.59)6.320.01 (0.01)0.03 (0.06)0.190.26 (0.14)0.77 (1.97)5.01
Biomass2.13 (1.23)2.03 (1.15)3.160.01 (0.01)0.04 (0.07)0.240.18 (0.08)0.94 (2.12)7.48
Conventional2.17 (2.56)2.07 (2.35)5.120.05 (0.04)0.05 (0.04)0.120.41 (0.46)0.92 (2.36)8.08
XericMastication0.57 (0.02)0.52 (0.11)0.640.07 (0.08)0.08 (0.1)0.350.27 (0.21)1.71 (4.47)15.74
Biomass0.75 (0.32)0.65 (0.31)1.110.03 (0.04)0.01 (0.04)0.080.19 (0.08)0.39 (1.50)4.57
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Murray, J.; Hagan, D.; Hiesl, P.; Baldwin, R. The Influence of Slash Management Practices on Water and Nutrient Dynamics in Longleaf Pine Forests. Forests 2022, 13, 1449. https://doi.org/10.3390/f13091449

AMA Style

Murray J, Hagan D, Hiesl P, Baldwin R. The Influence of Slash Management Practices on Water and Nutrient Dynamics in Longleaf Pine Forests. Forests. 2022; 13(9):1449. https://doi.org/10.3390/f13091449

Chicago/Turabian Style

Murray, Jacob, Donald Hagan, Patrick Hiesl, and Robert Baldwin. 2022. "The Influence of Slash Management Practices on Water and Nutrient Dynamics in Longleaf Pine Forests" Forests 13, no. 9: 1449. https://doi.org/10.3390/f13091449

APA Style

Murray, J., Hagan, D., Hiesl, P., & Baldwin, R. (2022). The Influence of Slash Management Practices on Water and Nutrient Dynamics in Longleaf Pine Forests. Forests, 13(9), 1449. https://doi.org/10.3390/f13091449

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop