Next Article in Journal
Comparing the Quantity and Structure of Deadwood in Selection Managed and Old-Growth Forests in South-East Europe
Next Article in Special Issue
Resistance of a Local Ecotype of Castanea sativa to Dryocosmus kuriphilus (Hymenoptera: Cynipidae) in Southern Italy
Previous Article in Journal
Night Light-Adaptation Strategies for Photosynthetic Apparatus in Yellow-Poplar (Liriodendron tulipifera L.) Exposed to Artificial Night Lighting
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Encroachment Dynamics of Juniperus virginiana L. and Mesic Hardwood Species into Cross Timbers Forests of North-Central Oklahoma, USA

1
Department of Natural Resource Ecology and Management, Oklahoma State University, 008 Agricultural Hall, Stillwater, OK 74078, USA
2
Bureau of Indian Affairs, Southern Plains Region, P.O. Box 368, Anadarko, OK 73005, USA
*
Author to whom correspondence should be addressed.
Forests 2018, 9(2), 75; https://doi.org/10.3390/f9020075
Submission received: 8 January 2018 / Revised: 29 January 2018 / Accepted: 1 February 2018 / Published: 3 February 2018
(This article belongs to the Special Issue Ecology and Management of Invasive Species in Forest Ecosystems)

Abstract

:
Cross Timbers forests, typically dominated by Quercus stellata Wangenh. and Q. marilandica Muenchh., are the transition zone between eastern deciduous forest and prairie in the southern Great Plains. Fire exclusion beginning in the mid-1900s has led to increasing stand density and encroachment of fire-intolerant Juniperus virginiana L. and mesic hardwood. We measured current forest structure and tree ages of 25 stands (130 plots) in north-central Oklahoma to characterize the extent and dynamics of encroachment. The respective basal area and stand density of the overstory (diameter at breast height; dbh > 10 cm) were 19.0 m2 ha−1 and 407 trees ha−1 with Q. stellata comprising 43% of basal area and 42% of stand density. Quercus marilandica represented only 3% of basal area and 4% of overstory density. Juniperus virginiana represented 7% of basal area and 14% of stand density while mesic hardwoods, e.g., Celtis spp., Ulmus spp., Carya spp., 33% of basal area and stand density. The sapling layer was dominated by mesic hardwoods (68%) and J. virginiana (25%) while the seedling layer was dominated by mesic hardwoods (74%). The majority of Quercus recruited into the overstory between 1910–1970, while recruitment of J. virginiana and mesic hardwoods began more recently (post 1950s). Growth rate, based on the relationship between age and dbh, was faster for mesic hardwoods than for J. virginiana and Q. stellata. These results indicate that removal of recurrent surface fire as a disturbance agent has significantly altered forest composition in the Cross Timbers region by allowing encroachment of J. virginiana and fire-intolerant, mesic hardwoods. This increases wildfire risk because J. virginiana is very flammable and will alter how these forests respond to future drought and other disturbance events.

1. Introduction

The Cross Timbers region of the United States is the ecotone between the eastern deciduous forest and the southern Great Plains, stretching from Texas north through Oklahoma into southern Kansas. This region is a mosaic comprising oak forest, oak savanna, and tallgrass prairie. The Cross Timbers region historically covered approximately 8 million hectares, about 4.5 million in Oklahoma [1]. The forest areas in the Cross Timbers are typically dominated by an overstory of Quercus stellata (post oak) and Q. marilandica (blackjack oak) that reach approximately 15 m in height [2].
The species composition and structure of Cross Timbers forests were maintained by frequent surface fire [3]. Measurements of the historical fire return interval pre-European American settlement ranged from 2.9–6.7 years [4,5,6,7,8]. Due to a combination of fire exclusion and favorable recruitment following canopy gaps that developed as a result of widespread drought in the 1950s, the Cross Timbers forests have seen an increase in stand basal area and an increase in fire-intolerant species such as Juniperus virginiana (eastern redcedar), Celtis spp., and Ulmus spp. [9]. This change in forest composition and structure in the Cross Timbers parallels the “mesophication” of much of the rest of the eastern deciduous forest, where fire-tolerant oak species are replaced with less fire-tolerant species due to fire exclusion [10]. In particular, the encroachment of highly flammable, evergreen J. virginiana into Cross Timbers forests is of serious concern. Re-measurement of stands initially sampled from 1953–1957 by Rice and Penfound [11] indicated that J. virginiana increased from 0.05 to 2.71 m2 ha−1 in terms of basal area and from 0.73 to 23.85 stems ha−1, while the overall stand basal area and stand density roughly doubled [9,12].
Juniperus virginiana is typically considered a shade-intolerant species; however, within the forest interior of the Cross Timbers, J. virginiana finds enough light availability in canopy gaps and is physiologically active during the winter/spring oak leaf-off period [13]. Once J. virginiana establishes, it shades out herbaceous vegetation [14] and adds litter to the forest floor that decomposes differently from oak leaves [15], alters fungal communities and litter chemistry [16], and burns much less readily than either oak litter or herbaceous vegetation. This reduces the available fine fuel load around J. virginiana trees and saplings and results in a forest that is incapable of maintaining a consistent, low-intensity surface fire. These low-intensity surface fires are necessary to purge J. virginiana seedlings (and other fire intolerant species), which become difficult to kill with prescribed fire once the trees reach a height of approximately 2 m [17,18].
In contrast to its effects on surface fire, J. virginiana encroachment increases the risk of wildfire. Nowacki and Abrams [10] proposed an alternative stable state for the oak forests of the eastern United States where fire exclusion leads to a more mesic forest that is resistant to all types of fire due to reduced fine fuel loading. However, the encroaching J. virginiana is more xeric, resinous, and burns extremely well under drought/wildfire conditions [19]. Hoff et al. [20] found that J. virginiana encroachment in north-central Oklahoma increased available fuel loading for wildfires by 38% through the inclusion of J. virginiana foliage and small branches. This reduction in ability for fine fuels to support a surface fire combined with the increased risk of crown fire from J. virginiana in the midstory sets the stage for wildfires that will be harder to suppress and more dangerous to human health and property. Increased risk and severity of wildfires also may change the regeneration pattern of Cross Timbers forests by replacing surface fires with stand-replacing fires [21].
Understanding the extent of encroachment of fire-intolerant species into Cross Timbers forests and the timeline of the invasion process is relevant to understanding regional forest dynamics and is important information for land managers. Previous work by Stambaugh et al. [21] and DeSantis et al. [9], based on tree aging, documented an increase in J. virginiana within Cross Timbers forests beginning in the mid 1900s. However, neither study looked at the timeline for invasion of more mesic species such as Celtis occidentalis L., Celtis laevigta Willdenow, or Ulmus americana L. In addition, the majority of studies on Cross Timbers forests have occurred on government controlled and managed lands (e.g., [4,5,7,8,22]) that do not reflect the typical disturbance pattern for land owned by small private landowners, which represents the vast majority of land holdings in the region. By knowing the dynamics of J. virginiana and mesic species invasion into Cross Timbers forests, land managers can prioritize and schedule management activities based on current conditions and projected rates of encroachment to reduce J. virginiana, reduce wildfire risk, and/or restore historical ecosystem composition or structure.
To address the issue of J. virginiana and mesic species encroachment into Cross Timbers forests, we measured current forest structure and condition during the summer of 2016 on Bureau of Indian Affairs (BIA) properties in north-central Oklahoma. In addition, we determined the age structure of the various species and related that to the dynamics of J. virginiana and mesic hardwood encroachment. Our objectives were to (1) quantify the current forest overstory, sapling, and seedling composition, (2) determine the extent of encroachment by fire-intolerant species, (3) establish a sub-decadal timeline of encroachment, and (4) compare growth rates of major species groups.

2. Materials and Methods

2.1. Study Area

We used 25 properties managed in trust for the Pawnee Tribe by the BIA, Southern Plains Region (Pawnee and Payne Counties Oklahoma, USA). These are typically landholdings ~65 ha in size for which the federal government holds legal title but the beneficial interest and management remain with the individual tribal member. We did not identify recent fires on any of the 25 properties. This subset of the BIA managed properties in the region (Figure 1) was randomly selected and had to meet two criteria; they had to be in excess of 16 ha and have approximately 50% oak forest canopy cover during leaf-on conditions. These 25 properties fell on a north-south swath from the Arkansas River in the north to the Cimarron River in the south (Figure 2). The eastern extent was limited by the availability of BIA properties in the Southern Plains region while the western extent was limited by the 50% oak forest requirement as Cross Timbers forests transitions into prairie. This area in Oklahoma averages 100 cm of annual precipitation. The average daytime high temperature is 15.5 °C, ranging from 34.0 °C in July to 8.8 °C in January. The growing season averages 197 days. Winds from the south and southeast dominate, averaging 14.5 km hr−1 (OK Climatological Survey 2016). Soils are predominantly silt loam complexes and feature prominent interspersed rock outcroppings [23].
Once the forest areas were identified on each of the 25 properties [20], five field measurement plots per property (one larger property had ten field measurement plots) were randomly located within the forest interior using ArcMAP (Esri, Redlands, CA, USA). Plot locations were found in the field using the onboard GPS capability of an iPad Air (Apple, Cupertino, CA, USA) with an accuracy of roughly 5 m. Field plot centers were marked with a metal stake and tagged. Plots were visited from May to August 2016, and data were collected using modified FIREMON protocol [25] designed for long-term fire effects monitoring.

2.2. Data Collection

Overstory hardwood tree data were collected using a circular fixed radius plot of 11.3 m (0.04 ha). All overstory trees, defined as having a diameter at breast height (dbh) greater than 10 cm, were measured 1.37 m above the ground to the nearest 0.5 cm and identified by species. On every 5th tree, starting with the most dominant tree (determined by relative size and crown dominance) on the plot and moving clockwise, we measured tree height to the nearest 0.1 m using a laser hypsometer (Truepulse 200b, Laser Technology Inc., Centennial, CO, USA) or a telescopic measuring rod (Hastings, Hastings, MI, USA). Cores (5.15 mm diameter) were taken at breast height using an increment corer (Haglof, Långsele, Sweden) on this subset of trees. We cored a minimum of three overstory hardwood trees per plot. If a plot had an insufficient number of mature stems to sample three stems by choosing every fifth, an alternative sampling strategy was used. For example, a plot with nine trees would have every third tree sampled and for a plot with only three trees, all were sampled.
Hardwood sapling data were collected within a smaller fixed radius subplot of 3.57 m (0.004 ha). Saplings were defined by a dbh between 0.1 and 10 cm. We identified the species of each stem, measured height to the nearest 0.1 m and recorded dbh as one of five size classes, i.e., 0–2, 2–4, 4–6, 6–8, 8–10 cm. We quantified seedling density and size within a subplot with radius of 1.78 m (0.001 ha) by recording a count by species and height class. Height class intervals were 0–0.25, 0.25–0.5, 0.5–0.75, 0.75–1.0 and 1.0–1.37 m.
All J. virginiana overstory trees, saplings, and seedlings were sampled on the same 11.3 m (0.04 ha) fixed radius plot used for mature hardwood trees. All J. virginiana seedlings were counted into the same height class intervals as hardwood seedlings. Saplings over 1.37 m and under 3 m tall were recorded by count, and all J. virginiana larger than 3 m tall were measured for dbh to the nearest 0.5 cm and height to the nearest 0.1 m. Every fifth J. virginiana with dbh > 0.1 cm beginning with the most dominant and moving clockwise was selected for ageing. A core was taken if the tree was of sufficient size (approximately 10 cm dbh); otherwise, a stem disk was taken approximately 15 cm above ground level using a bow saw. An alternative numerical sampling strategy was used as described above if necessary to age a minimum of three J. virginiana stems per plot; however, because not all plots had J. virginiana present, some plots had less than three trees sampled.
We used two different methods to record standing dead tree (snag) density. We initially assumed snags would be rare and recorded snags by count and species on a circular plot with a radius of 25 m (0.19 ha). Because snag density was higher than expected, we began to collect snag data within the overstory tree plot (radius 11.3 m) and included dbh measurements from each snag and species when possible (n = 53 plots). In this manuscript, snag relative densities are reported from this subset of plots with more rigorous data collection.

2.3. Age

Juniperus virginiana will produce false rings. Because of this and the extremely narrow growth rings of some slower growing hardwoods, all cores were aged using a variable magnification sliding stage scope after progressive sanding to 400 grit using a mechanical sander. Additional hand sanding was undertaken if necessary using 3000 grit sandpaper. Cores usually included the pith. When the pith was not included in the core, the core was either discarded or the age to the pith was estimated based on growth trends for that core. This was only done when the core was judged to be close to the center (~3 rings). Stem disks were sanded and aged using a 12× hand lens or microscope. Linear regression analysis for the relationship between age (ground level) and height of J. virginiana between 1.37 and 3 m tall in our data set determined that it took approximately 9 years (8.87 year) for J. virginiana to reach a height of 1.37 m. Accordingly, we adjusted the J. virginiana ages for those measured at breast height by adding 9 years. When comparing to hardwoods, unadjusted ages of cored J. virginiana were used for comparison to hardwood ages because those were based on cores taken at breast height.

2.4. Analyses

To simplify the analyses of species relative density, we grouped similar species together. Quercus stellata, Q. marilandica and J. virginiana were analyzed as stand-alone species. The two Carya species, C. texana Buckley and C. illinoinensis (Wangenh.) K. Koch, were grouped together due to the small sample size of C. illinoinensis and similar fire tolerance. Celtis leavigata Willdenow and C. occidentalis L., which frequently hybridize and share high fire sensitivity, were grouped together. We encountered one Ulmus rubra Muhl. and included it with U. americana. We encountered both Q. velutina Lam. and Q. shumardii Buckland and grouped them together in a Red Oak group. Species that typically never reach the 10 cm dbh threshold were classified as Shrubs and were analyzed separately for relative frequency. All remaining tree species were considered Other, and had a relative density of no more than 2% in the overstory and sapling categories or 5% in the seedling category. Five-year precipitation averages were calculated using data available from the National Oceanic and Atmospheric Administration, National Centers for Environmental Information for the north central region of Oklahoma [26].
Slopes and least squares means (LS means) for the relationship between age and dbh were tested to determine if there were differences between J. virginiana, Q. stellata, and a combined Mesic species group comprising Carya spp., Celtis spp. and Ulmus spp. using Analysis of Covariance (ANCOVA) using mixed model analyses (PROC MIXED) of SAS version 9.3 (SAS Institute, Cary, NC, USA). The age vs. height and age vs. dbh relationships using ages corrected to ground level were determined using regression analysis (PROC REG) in SAS version 9.3 to calculate the growth rate of J. virginiana.

3. Results

3.1. Forest Composition

Overstory Q. stellata averaged approximately 23 cm dbh and 12.0 m in height (Table 1). It was the most frequently encountered tree species representing 171 trees ha−1 and 8.2 m2 ha−1 of basal area per plot and occurred on 92% of properties and 78% of plots. In contrast, Q. marilandica represented only 16 trees ha−1 and 0.7 m2 ha−1 of basal area (Table 2) and occurred on 52% of properties and 21% of plots (Table 3). Quercus marilandica were slightly smaller than Q. stellata, averaging 21 cm dbh and 10.7 m in height (Table 1). Quercus velutina were relatively large, i.e., 32 cm dbh and 16.3 m tall (Table 1), but only represented 2.6 m2 ha−1 basal area and 26 stems ha−1 (Table 2) occurring on 72% of properties and 35% of plots (Table 3). Carya texana (52% of properties, 22% of plots), C. occidentalis (68% of properties, 35% of plots), C. leavigata (60% of properties, 22% of plots), and U. americana (88% of properties, 38% of plots) made up the remainder of the commonly encountered hardwood tree species in the overstory (Table 3). Juniperus virginiana comprised 57 trees ha−1 and 1.3 m2 ha−1 basal area (Table 2) and occurred on 96% of properties and 53% of plots (Table 3) but were smaller than hardwood trees with an average dbh of 16 cm and average height of 8.6 m (Table 1). The average number of snags was 126.2 ha−1 with an average dbh of 23 cm.
Juniperus virginiana was common in the sapling layer (100% of properties and 84% of plots; Table 3) with an average density of 232 saplings ha−1 (Table 2). Carya texana (67 stems ha−1, 40% of properties and 16% of plots), C. occidentalis (271 stems ha−1, 56% of properties and 30% of plots) and Sapindus drummondii (67 stems ha−1, 16% of properties and 4% of plots) were also tree species that were major contributors to the sapling layer (Table 2 and Table 3). The shrub species Cercis canadensis and Cornus drummondii were frequently encountered sapling-sized stems with densities of 205 and 325 stems ha−1, respectively (Table 2). Conspicuously absent from the sapling size class were the traditionally dominant overstory tree species, Q. marilandica and Q. stellata, which together accounted for only 62 stems ha−1. Ulmus americana was also present in relatively small quantity of 52 stems ha−1. In the seedling category, U. americana accounted for 3331 ha−1 (100% of properties and 61% of plots) more than any other species. For the seedling size class, Celtis spp. were all classified as C. occidentalis due to difficulties in identification. This group averaged 2600 seedlings ha−1 (92% of properties and 61% of plots). Juniperus virginiana only had 171 stems ha−1 in the seedling layer.

3.2. Relative Dominance and Density

Most tree species had relative densities similar in magnitude to relative basal area. Quercus stellata was the dominant overstory species comprising 43% and 42% of the basal area and stem density distribution, respectively (Figure 3). Quercus marilandica, a historically abundant overstory tree in Cross Timbers forests, represented only 3% of the basal area and 4% of the relative stem density of overstory trees (Figure 3). The Red Oak group, dominated by Q. velutina and commonly associated with stream courses in the Cross Timbers, represented an outsized proportion of basal area (14%) compared to its relative density of 7% (Figure 3) due to large average dbh. Juniperus virginiana accounted for 14% of the relative density in the overstory tree size class but only 7% of the relative basal area of (Figure 3) due to smaller than average dbh. The more mesic species, i.e., Cetlis spp., Ulmus spp., and Carya spp. (we lumped C. texana with the ‘mesic’ species as it has a similar fire sensitivity and litter dynamics as the more mesic species) combined, accounted for 22% of the relative basal area and 23% of relative stem density in the overstory (Figure 3). The species distribution of snags closely mirrored the mature tree stem distribution (Figure 3). The exception was Q. marilandica, which had relative density in the snag distribution almost twice as great as its relative density in the overstory. While occurrence of Biscogniauxia spp. (formerly Hypoxylon) canker was not formally recorded, almost every Q. marilandica snag displayed signs of infection, as did some live stems. No other species were noticed to have signs of Biscogniauxia spp. canker on live trees. Quercus stellata snags frequently displayed some signs of Biscogniauxia spp., but it was not as prevalent as on Q. marilandica snags.
Juniperus virginiana stems represented 25% of the sapling relative density and the more mesic species, i.e., Cetlis spp., Ulmus spp., and Carya spp., made up an additional 48% of the saplings (Figure 3). In contrast, all Quercus species combined made up only 13% of tree saplings (Figure 3). Combined, Cetlis spp., Ulmus spp., and Carya spp. made up 60% of the seedlings (Figure 3). Juniperus virginianas relative density in the seedling size class was 2%, which was much lower than for the sapling and overstory size classes. Combined, the Quercus species comprised 24% of seedlings.
Cornus drumundii (53%) and C. canadensis (34%) comprised the majority of shrubs in the sapling size class. Prunus mexicana composed 8% and the ‘other’ category an additional 4%. Rhus copallina is an early successional species that was commonly encountered on forest edges, but rare within the forest interior (1%) where inventory plots were located. In total, all shrubs made up 623 stems ha−1 of 1542 total stems ha−1 in the sapling layer (40.4%). Within the seedling layer, shrub species comprised 2323 stems ha−1 out of 13,202 stems ha−1 total (17.5%) (Table 2).

3.3. Age Distribution

The species distribution for hardwood tree age reported in bi-decadal increments based on pith age at breast height was dominated by Q. stellata, particularly at older ages (Figure 4). The oldest tree sampled was a 233-year-old Q. stellata, which reached a height of 1.37 m between 1780 and 1785. The average Q. stellata sampled for age was 75 years old, with 83% of stems sampled between age 50 and 100 years. Quercus marilandica and Red Oaks, with a mean age of 63 and 73% stems between 40 and 90 years of age, were on average slightly younger than Q. stellata, but older than the more mesic species. The Cetlis spp., Ulmus spp., and Carya spp. stems were younger and began to reach the sapling size class after the 1950s with only six out of 92 stems sampled over 70 years of age. This group averaged 44 years of age, and 76% of stems sampled were between 20 and 60 years old. This parallels the results from the relative frequency data that seem to indicate recent and ongoing recruitment of the mesic species given their abundance in the sapling and seedling size classes. Measured at breast height, J. virginiania began to reach the sapling size class beginning in the 1950s.
The J. virginiana age distribution of all measured trees based on pith age at ground level is presented in bi-decadal increments (Figure 5). The oldest stem dated to the last years of the 1920s. The majority of J. virginiana recruitment began in the 1950s, and the rate of recruitment increased through the early 2000s. The rate of increase between 1940 and 2005 was mostly linear with some variation associated with differences in 5-year precipitation averages. For instance, less recruitment during the 1955–1960 and 1975–1980 periods corresponded to periods of below average precipitation. Recruitment during the 2005–2010 and 2010+ intervals was lower because we did not age trees shorter than 1.4 m and many trees that recently established were thus not sampled as it takes approximately 9 years for J. virginiana to reach a height of 1.4 m.
The relationship between age at ground-level and dbh for J. virginiana was best characterized by a power function (Figure 6). This equation form was superior to a linear relationship because it provided a better fit for stems with small dbh. A stem at 20 years of age added approximately 0.30 cm dbh per year, while a stem age 50 added approximately 0.45 cm dbh annually. The relationship between age and height was adequately modeled with a linear regression forced through the origin. Juniperus virginiana grew approximately 16 cm annually in height.

3.4. Age Size Relationships between Groups

The relationship between dbh and age measured at breast height for J. virginiana, Q. stellata, and a combined grouping of Cetlis spp., Ulmus spp., and Carya spp. is presented in Figure 7. There was a C. occidentalis tree (confirmed) with a dbh of 127 cm that was 60 cm greater than the next largest tree measured for age. If the ‘outlier’ tree was removed from the analysis, the slopes were not significantly different from one another (p = 0.12). However, analysis of LS Means indicated the dbh of the mesic species grouping (28.28 ± 1.16 s.e.) was greater than Q. stellata (21.12 ± 0.82 s.e.) and J. virginiana (22.47 ± 1.16 s.e.) (p < 0.0001), indicating greater dbh for the mesic species at a given age (Figure 7). When including this ‘outlier’ tree, the slopes of the relationship between dbh and age were significantly greater for the mesic species as were the LS Means. Pairwise comparisons indicated that J. virginiana and Q. stellata slopes between age and dbh were not different (p = 0.66), but the slope of the Cetlis spp., Ulmus spp., and Carya spp. grouping was significantly greater than Q. stellata (p = 0.0001) and J. virginiana (p = 0.001). Analysis of LS Means indicated that Q. stellata (21.26 ± 0.89 s.e.) and J. virginiana (21.13 ± 1.23 s.e.) were not significantly different (p = 0.21), but that both were lower than the mesic species grouping (30.02 ± 1.13 s.e.) (p < 0.0001).

4. Discussion

The areas of privately managed Cross Timbers forests in north-central Oklahoma we measured are undergoing encroachment by J. virginiana, densification (increase in tree density and basal area), and mesophication (increase in Celtis spp., Ulmus spp., and other mesic, fire-intolerant species). The age structure as well as relative abundance of species in the overstory, sapling, and seedlings size classes indicate that this is a fairly recent and ongoing transformation. These changes have important implications for the future forest condition and successional processes related to drought and fire as well as management options for applying prescribed fire and restoration activities.
Our overstory age data agree with previous studies from Cross Timbers forests [9,21] reporting that J. virginiana recruitment accelerated after the 1950s. This coincided with fire suppression in the post WWII era and droughts during the 1950s, which opened canopy gaps due to tree mortality. Likewise, the relative density distributions among the overstory, sapling, and seedling size classes from our study support a mesophication and densification narrative, as the fire-intolerant tree species are most abundant in the sapling and seedling size classes. Three stands in north-central Oklahoma first measured in the middle of 1950s [11] and later remeasured in the late 2000s [9] exhibited an increase in basal area from 10.4 to 22.4 m2 ha−1, an increase in overstory density from 531 to 754 trees ha−1, and an increase in overstory J. virginiana basal area from 0 to 0.7 m2 ha−1 and density from 0 to 82 trees ha−1 over a period of 50 years. The stand conditions we measured were similar to those measured in the 2000s by DeSantis et al. [9] with the exception that we had approximately twice the basal area of J. virginiana, which could be related to an additional decade for J. virginana growth between the two studies.
While Q. stellata dominated the overstory, the most abundant species in the sapling size class were fire-intolerant trees and shrubs, e.g., C. drummondii, C. occidentallis, J. virginiana, C. canadensis, S. drummondii, and U. americana which indicates a potential transition away from a Quercus-dominated overstory if fire continues to be excluded in future. Similarly, the most abundant seedling species were C. occidentalis and U. americana. Although scarce in the sapling size class, Quercus spp. likely will become relatively more abundant if fire is reintroduced. Quercus stellata and Q. marilandica seedlings can survive in the understory for decades, allocating much of their acquired carbon to growing root systems which allows them to resprout multiple times following topkill from fire [27].
The faster diameter development of the fire-intolerant, mesic grouping of the Celtis spp., Ulmus spp., Carya spp. further indicates that once established and free from disturbance, the rate of meshophication may accelerate. In contrast, the dbh of more xeric J. virginiana grew approximately at the same rate as did Q. stellata. In regards to wildfire fuels, the growth of J. virginiana is particularly important. On these same sites, Hoff et al. [20] measured that J. virginiana composed 18.4% of canopy cover and added 6.3 Mg ha−1 to available wildfire fuel. Given the large number of J. virginiana in the sapling size class, these values are in the process of increasing further. Height and dbh growth rates of J. virginiana from our study were similar to another report of J. virginiana growing in a forest [28], but 20% less in height growth rate and 56% less in diameter growth rate than those reported for J. virginiana in prairie environments [18]. The relatively slow growth of J. virginiana growing in the forest interior indicates that it would take approximately 12 years to reach 2 m tall at which time they are more difficult to kill with prescribed fire. This is helpful for planning prescribed fire and indicates it may be possible to burn less frequently than the historical fire return interval (3–6 years) to accomplish management objectives related to J. virginiana reduction.
There were fewer J. virginiana seedlings than might have been expected from their high relative density in the sapling size class. The dips in J. virginiana recruitment associated with below average 5-year precipitation suggests that drought may reduce seedling/sapling establishment and survival. The recent period ending in 2015 was the driest since 1975–1980, which may help explain the low numbers of J. virginiana in the current seedling size class. Juniperus virginiana transpires water all year in the Cross Timbers [13]. While extremely drought-tolerant, J. virginiana seedlings may be more susceptible to drought than larger J. virginiana trees because seedlings have smaller, shallower root systems. Cones of J. virginina typically disperse February to March and seeds germinate the following spring such that an increase in seedlings germination from seed stored in the seedbed should not be expected post drought. Alternative explanations for the relatively low frequency of J. virginiana seedlings could be that survival of J. virginiana seedlings might be greater such that fewer seedlings are needed to achieve successful recruitment into the sapling size class or, the rate of new J. virginiana recruitment may be beginning to wane if forest density has increased the point of hindering seedling establishment.
Cross Timbers forests are considered to be dominated by Q. stellata and Q. marilandica, yet we encountered relatively few Q. marilandica in the overstory, sapling, or seedling classes. In the 1950s Rice and Penfound [11] measured that Q. marilandica had 55% and 73% as much basal area stem density as did Q. stellata. When remeasured in the 2000s by DeSantis et al. [12], Q. marilandica had only 26% and 41% as much basal area and stem density, respectively, as did Q. stellata. The decline of the relative dominance and density of Q. marilandica was due to both a decrease in Q. marilandica as well as an increase in Q. stellata [12]. While no historical data are available for the stands measured in our study, greater relative snag density for Q. marilandica probably indicates greater mortality of mature trees for this species, which reduces dominance. Drought is a predisposing factor for many pathogens including Biscogniuxia spp. infections on Quercus spp. [29]. In Cross Timbers forests, Q. marilandica seemed predisposed to infection. The years of 2011 and 2012 were among the driest years in north-central Oklahoma since the late 1970s [30]. Our field observations agree with DeSantis et al. [12] and Masters and Waymire [31], who reported high occurrence of canker on Q. marilandica in response to drought and competition. Without significant disturbance events, particularly fire to stimulate regeneration and reduce competition, Q. marilandica are disappearing and will likely continue to decline.
The Cross Timbers region frequently undergoes moderate to severe droughts, which may be exacerbated by increased competition for water within current, denser forest conditions. The increasingly more common mesic species may experience greater drought mortality during extreme events as they lack the drought tolerance of the formerly dominant oak species. This potential for increased drought mortality may further increase encroachment of J. virginiana, as it has a high physiological drought tolerance [32,33]. Therefore, J. virginiana will likely survive drought better than the mesic species [34] and fill any canopy gaps resulting from forest dieback. Future fire regimes also will be affected by ongoing changes within the Cross Timbers forests, particularly by J. virginiana encroachment. Increased fuel loading for wildfires results from this encroachment [20] and has the potential to cause stand replacing wildfires or widespread damage to overstory oaks [21]. Stand replacing fires do not have an analog in the known fire history of the Cross Timbers forests. Given the high flammability of J. virginiana, the mesophication process of Cross Timbers forests produces a different potential outcome than further east. In the Eastern Deciduous Forest, mesophication causes a loss of fire as a potential disturbance altogether. However, the encroachment of J. virginiana into Cross Timbers forests may reduce land manager’s ability to use prescribed fire, but increase the likelihood of stand replacing, crown fires.

5. Conclusions

Our data provide important insight into mesophication related to fire exclusion of Cross Timbers forests. The encroachment of J. virginiana in the tree and sapling size classes and the increase of mesic, fire-intolerant hardwood species such as Celtis spp., Carya spp. and Ulmus spp. in the sapling and seedling classes provide evidence that the composition of Cross Timbers forests is changing. The linear increase in J. virginiana recruitment into the overstory indicates that J. virginiana likely will continue to increase in dominance as the younger trees become older and larger. In addition, the presence of mesic, fire-intolerant species is expected to increase in the absence of fire due to their abundance in the sapling and seedling size classes as well as their faster growth relative to Quercus spp. or J. virginiana. These data and other reports [9,12,21] provide compelling evidence that removal of recurrent surface fire as a disturbance agent is significantly altering the trajectory of forest composition in the Cross Timbers. This has significant implications for wildfire risk and how these forests will respond to future disturbance events.

Acknowledgments

Special thanks to: OSU interns Giovanni Serreau and Caren Mendonica, Bureau of Indian Affairs Pathways students Ryan Morrison and Joe Rule for their hard work in the field, and Chris Clay for facilitating access to BIA trust lands. This work was supported by the Bureau of Indian Affairs project A16AC00016, the Oklahoma Agricultural Experiment Station, and McIntire Stennis project OKL0 2929.

Author Contributions

D.L.H., R.E.W., C.B.Z., and N.D.L. conceived and designed the study; D.L.H. collected the data; D.L.H. and R.E.W. analyzed the data; D.L.H, R.E.W., C.B.Z., and N.D.L. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kuchler, A.W. Potential Natural Vegetation of the Conterminous United States; Special Publication No. 36; American Geographical Society: New York, NY, USA, 1964. [Google Scholar]
  2. Therrell, M.D.; Stahle, D.W. A predictive model to locate ancient forests in the Cross Timbers of Osage County, Oklahoma. J. Biogeogr. 1998, 25, 847–854. [Google Scholar] [CrossRef]
  3. Hallgren, S.W.; DeSantis, R.D.; Burton, J.A. Fire and vegetation dynamics in the Cross Timbers forests of South-Central North America. In Proceedings of the 4th Fire in Eastern Oak Forests Conference, Springfield, MO, USA, 17–19 May 2011; Dey, D.C., Stambaugh, M.C., Clark, S.L., Schweitzer, C.J., Eds.; USDA Forest Service, Northern Research Station: Newtown Square, PA, USA, 2012. [Google Scholar]
  4. Clark, S.L.; Hallgren, S.W.; Stahle, D.W.; Lynch, T.B. Characteristics of the Keystone Ancient Forest Preserve, an old-growth forest in the Cross Timbers of Oklahoma. Nat. Area J. 2005, 25, 165–175. [Google Scholar]
  5. Stambaugh, M.C.; Guyette, R.P.; Godfrey, R.; McMurry, E.R.; Marschall, J.M. Fire, drought, and human history near the Western terminus of the Cross Timbers, Wichita Mountains, Oklahoma, USA. Fire Ecol. 2009, 5, 51–65. [Google Scholar] [CrossRef]
  6. DeSantis, R.D.; Hallgren, S.W.; Stahle, D.W. Historic fire regime of an upland oak forest in south-central North America. Fire Ecol. 2010, 6, 45–61. [Google Scholar] [CrossRef]
  7. Allen, M.S.; Palmer, M.W. Fire history of a prairie/forest boundary: More than 250 years of frequent fire in a North American tallgrass prairie. J. Veg. Sci. 2011, 22, 436–444. [Google Scholar] [CrossRef]
  8. Stambaugh, M.C.; Sparks, J.C.; Guyette, R.P.; Willson, G. Fire history of a relict oak woodland in northeast Texas. Rangel. Ecol. Manag. 2011, 64, 419–423. [Google Scholar] [CrossRef]
  9. DeSantis, R.D.; Hallgren, S.W.; Stahle, D.W. Drought and fire suppression lead to rapid forest composition change in a forest-prairie ecotone. For. Ecol. Manag. 2011, 261, 1183–1840. [Google Scholar] [CrossRef]
  10. Nowacki, G.J.; Abrams, M.D. The demise of fire and “mesophication” of forests in the Eastern United States. BioScience 2008, 58, 123–138. [Google Scholar] [CrossRef]
  11. Rice, E.L.; Penfound, W.T. The upland forests of Oklahoma. Ecology 1959, 40, 593–608. [Google Scholar] [CrossRef]
  12. DeSantis, R.D.; Hallgren, S.W.; Lynch, T.B.; Burton, J.A.; Palmer, M.W. Long-term directional changes in upland Quercus forests throughout Oklahoma, USA. J. Veg. Sci. 2010, 21, 606–615. [Google Scholar] [CrossRef]
  13. Caterina, G.L.; Will, R.E.; Turton, D.J.; Wilson, D.S.; Zou, C.B. Water use of Juniperus virginiana trees encroached into mesic prairies in Oklahoma, USA. Ecohydrology 2014, 7, 1124–1134. [Google Scholar]
  14. Van Els, P.; Will, R.E.; Palmer, M.W.; Hickman, K.R. Changes in forest understory associated with Juniperus encroachment in Oklahoma, USA. Appl. Veg. Sci. 2010, 13, 356–368. [Google Scholar]
  15. Nunes-Biral, V.C. Impacts of Tree Encroachment into a Grassland on Soil Carbon and Inputs in Northcentral Oklahoma. Master’s Thesis, Oklahoma State University, Stillwater, OK, USA, 2017. [Google Scholar]
  16. Williams, R.J.; Hallgren, S.W.; Wilson, G.W.T.; Palmer, M.W. Juniperus virginiana encroachment into upland oak forests alters arbuscular mycorrhizal abundance and litter chemistry. Appl. Soil Ecol. 2013, 65, 23–30. [Google Scholar] [CrossRef]
  17. Ansley, R.J.; Wiedemann, H.T. Reversing the woodland steady state: Vegetation responses during restoration of Juniperus-dominated grasslands with chaining and fire. In Western North American Juniperus Communities: A Dynamic Vegetation Type; Van Auken, O.W., Ed.; Springer: New York, NY, USA, 2008; pp. 272–290. [Google Scholar]
  18. Owensby, C.E.; Blan, K.R.; Eaton, B.J.; Russ, O.G. Evaluation of eastern redcedar infestations in the Northern Kansas Flint Hills. J. Range Manag. 1973, 26, 256–260. [Google Scholar] [CrossRef]
  19. Twidwell, D.; Fuhlendorf, S.D.; Taylor, C.A., Jr.; Rogers, W.E. Refining thresholds in coupled fire-vegetation models to improve management of encroaching woody plants in grasslands. J. Appl. Ecol. 2013, 50, 603–613. [Google Scholar] [CrossRef]
  20. Hoff, D.; Will, R.E.; Zou, C.B.; Weir, J.; Gregory, M.; Lillie, N. Estimating increased fuel loading within the Cross Timbers forest matrix of Oklahoma, USA due to an encroaching conifer, Juniperus virginiana, using leaf-off imagery. For. Ecol. Manag. 2018, 409, 215–224. [Google Scholar] [CrossRef]
  21. Stambaugh, M.C.; Marschall, J.M.; Guyette, R.P. Linking fire history to successional changes of xeric oak woodlands. For. Ecol. Manag. 2014, 320, 83–95. [Google Scholar] [CrossRef]
  22. Bragg, D.C.; Stahle, D.W.; Cerny, K.C. Structural attributes of two old-growth Cross Timbers stands in western Arkansas. Am. Midl. Nat. 2012, 167, 40–55. [Google Scholar] [CrossRef]
  23. Scott, G.F.; Boyett, J.B.; Dennis, J.; White, J.D. Supplement to the Soil Survey of Pawnee County, Oklahoma; USDA Natural Resource Conservation Service: Washington, DC, USA, 2006. [Google Scholar]
  24. Olson, D.M.; Dinerstein, E.; Wikramanayake, E.D.; Burgess, N.D.; Powell, G.V.N.; Underwood, E.C.; D’Amico, J.A.; Itoua, I.; Strand, H.E.; Morrison, J.C.; et al. Terrestrial Ecoregions of the world: A new map of life on earth. BioScience 2001, 51, 933–938. [Google Scholar] [CrossRef]
  25. Lutes, D.C.; Keane, R.E.; Caratti, J.F.; Key, C.H.; Benson, N.C.; Sutherland, S.; Gangi, L.J. FIREMON: Fire Effects Monitoring and Inventory System; USDA Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2006. [Google Scholar]
  26. Southern Climate Impacts Planning Program. Climate Trends for Oklahoma North-Central Region. Available online: http://charts.srcc.lsu.edu/trends/ (accessed on 2 October 2017).
  27. Clark, S.L.; Hallgren, S.W. Dynamics of oak (Quercus marilandica and Q. stellata) reproduction in an old-growth Cross Timbers forest. Southeast. Nat. 2003, 2, 559–574. [Google Scholar]
  28. Lykin, R.W. Estimation of Aboveground Eastern Redcedar Biomass. Master’s Thesis, Oklahoma State University, Stillwater, OK, USA, 1995. [Google Scholar]
  29. Despres-Loustau, M.L.; Marcais, B.; Nageleisen, L.M.; Piou, D.; Vannini, A. Interactive effects of drought and pathogens in forest trees. Ann. For. Sci. 2006, 63, 597–612. [Google Scholar] [CrossRef]
  30. Oklahoma Climatological Survey. The Climate of Pawnee County. Available online: http://climate.mesonet.org/county_climate/Products/County_Climatologies/county_climate_pawnee.pdf (accessed on 26 May 2017).
  31. Masters, R.E.; Waymire, J.R. Oak savanna restoration: Oak response to fire and thinning through 28 years. In Proceedings of the 4th Fire in Eastern Oak Forests Conference, Springfield, MO, USA, 17–19 May 2011; Dey, D.C., Stambaugh, M.C., Clark, S.L., Schweitzer, C.J., Eds.; USDA Forest Service, Northern Research Station: Newtown Square, PA, USA, 2012. [Google Scholar]
  32. Sperry, J.S.; Tyree, M.T. Water-stress-induced xylem embolism in three species of conifers. Plant Cell Environ. 1990, 13, 427–436. [Google Scholar] [CrossRef]
  33. Willson, C.J.; Manos, P.S.; Jackson, R.B. Hydraulic traits are influenced by phylogenetic history in the drought-resistant, invasive genus Juniperus (Cupressaceae). Am. J. Bot. 2008, 95, 299–314. [Google Scholar] [CrossRef] [PubMed]
  34. Volder, A.; Tjoelker, M.G.; Briske, D.D. Contrasting physiological responsiveness of establishing trees and a C4 grass to rainfall events, intensified summer drought, and warming in oak savanna. Glob. Chang. Biol. 2010, 16, 3349–3362. [Google Scholar] [CrossRef]
Figure 1. Study area represented by the black rectangle. The lightly shaded area is the Terrestrial Ecoregions of the World [24] designation for the transition between central forest and grasslands. The dark area represents the Cross Timbers from the Potential Natural Vegetation Type for the Continental U.S. [1]. This ecotone is commonly referred to as the Cross Timbers in Texas, Oklahoma and Kansas.
Figure 1. Study area represented by the black rectangle. The lightly shaded area is the Terrestrial Ecoregions of the World [24] designation for the transition between central forest and grasslands. The dark area represents the Cross Timbers from the Potential Natural Vegetation Type for the Continental U.S. [1]. This ecotone is commonly referred to as the Cross Timbers in Texas, Oklahoma and Kansas.
Forests 09 00075 g001
Figure 2. BIA-managed properties in this area of the Southern Plains region are outlined in black. The 25 properties randomly selected for field measurements in this study are circled in red.
Figure 2. BIA-managed properties in this area of the Southern Plains region are outlined in black. The 25 properties randomly selected for field measurements in this study are circled in red.
Forests 09 00075 g002
Figure 3. Tree relative dominance based on basal area (Tree RBA) as well as overstory tree, sapling, seedling, and snag relative density (RD) measured in Cross Timbers forests of north-central Oklahoma. Categories are either individual species or groupings of related species.
Figure 3. Tree relative dominance based on basal area (Tree RBA) as well as overstory tree, sapling, seedling, and snag relative density (RD) measured in Cross Timbers forests of north-central Oklahoma. Categories are either individual species or groupings of related species.
Forests 09 00075 g003
Figure 4. Major tree species associations measured in Cross Timbers forests of north-central Oklahoma presented in bi-decadal increments. All data were collected using increment cores at breast height such that each data point represents the 5-year period when the stem was recruited into the sapling size class. Juniperus virginiana were sampled separately from the other species and are offset for clarity of display.
Figure 4. Major tree species associations measured in Cross Timbers forests of north-central Oklahoma presented in bi-decadal increments. All data were collected using increment cores at breast height such that each data point represents the 5-year period when the stem was recruited into the sapling size class. Juniperus virginiana were sampled separately from the other species and are offset for clarity of display.
Forests 09 00075 g004
Figure 5. All Juniperus virginiana stems by bi-decadal age class with 5-year precipitation average. Data were based on stem disks at ground level or from cores at breast height that were corrected to year of establishment. Young J. virginiana with establishment dates after 2007 are not expected to have reached breast height and were not included in the sampling.
Figure 5. All Juniperus virginiana stems by bi-decadal age class with 5-year precipitation average. Data were based on stem disks at ground level or from cores at breast height that were corrected to year of establishment. Young J. virginiana with establishment dates after 2007 are not expected to have reached breast height and were not included in the sampling.
Forests 09 00075 g005
Figure 6. Relationship between age and dbh and relationship between age and height for Juniperus virginiana stems. Data were based on stem disks at ground level or from cores at breast height that were corrected to year of establishment. Trees shorter than 1.37 m were not sampled for age.
Figure 6. Relationship between age and dbh and relationship between age and height for Juniperus virginiana stems. Data were based on stem disks at ground level or from cores at breast height that were corrected to year of establishment. Trees shorter than 1.37 m were not sampled for age.
Forests 09 00075 g006
Figure 7. Relationship between age and dbh for Juniperus virginiana (J.v.), Quercus stellata (Q.s.), and a combined Mesic group of Carya spp., Celtis spp., and Ulmus spp. (CCU). All ages were taken at dbh and not adjusted. The circled outlier for the mesic species grouping was removed for calculation of relationship between age and dbh.
Figure 7. Relationship between age and dbh for Juniperus virginiana (J.v.), Quercus stellata (Q.s.), and a combined Mesic group of Carya spp., Celtis spp., and Ulmus spp. (CCU). All ages were taken at dbh and not adjusted. The circled outlier for the mesic species grouping was removed for calculation of relationship between age and dbh.
Forests 09 00075 g007
Table 1. Number encountered (N) and average dbh and height of trees (dbh > 10 cm), sapling (dbh 0.1 to 10 cm) and seedlings (shorter than 1.37 m) from 130 plots measured in Cross Timbers forests of north-central Oklahoma. Std. = standard deviation. Trees were measured in plots 0.04 ha in size (5.2 ha total). For hardwoods, saplings were measured in plots 0.004 ha in size (0.52 ha total) and seedlings in plots 0.001 ha in size (0.13 ha total). Juniperus virginiana L. overstory, saplings and seedlings were measured in the 0.04 ha plots (5.2 ha total). For seedlings, all Celtis spp. are lumped as Celtis occidentalis L.
Table 1. Number encountered (N) and average dbh and height of trees (dbh > 10 cm), sapling (dbh 0.1 to 10 cm) and seedlings (shorter than 1.37 m) from 130 plots measured in Cross Timbers forests of north-central Oklahoma. Std. = standard deviation. Trees were measured in plots 0.04 ha in size (5.2 ha total). For hardwoods, saplings were measured in plots 0.004 ha in size (0.52 ha total) and seedlings in plots 0.001 ha in size (0.13 ha total). Juniperus virginiana L. overstory, saplings and seedlings were measured in the 0.04 ha plots (5.2 ha total). For seedlings, all Celtis spp. are lumped as Celtis occidentalis L.
Overstory TreesSaplingsSeedlings
dbh (cm)Height (m)dbh (cm)Height (m)Height (m)
Genus SpeciesNMeanStd.MeanStd.NMeanStd.MeanStd.NMeanStd.
Acer negundo L.221.014.1--54.62.64.91.70--
Carya illinoinensis (Wangenh.) K.Koch3727.814.113.35.544.53.05.22.880.40.2
Carya texana Buckley9018.38.511.94.4364.32.84.22.3620.40.3
Celtis laevigata Willdenow5019.19.711.14.5212.82.03.61.90--
Celtis occidentalis L.16420.814.111.34.71412.42.13.11.63380.40.3
Cercis canadensis L.410.91.0--1072.42.12.91.2690.50.3
Cornus drummondii C.A.Mey.0----1491.10.72.00.91580.80.4
Diospyros virginiana L.1018.07.710.55.184.02.44.21.580.60.3
Fraxinus spp. L.3628.613.815.74.663.02.53.52.150.80.3
Gleditsia triacanthos L.112.5---13.0-4.4-30.50.1
Gymnocladus dioicus (L.) K. Koch621.822.212.95.50----0--
Juglans nigra L.2323.816.010.34.936.33.14.71.70--
Juniperus virginiana L.29616.46.28.62.112103.02.43.21.58910.70.4
Morus rubra L.2919.59.412.03.70----50.50.2
Platanus occidentalis L.1136.713.918.70.10----0--
Prunus angustifolia Marsh.211.01.4--0----0--
Prunus mexicana S.Watson812.92.46.42.5261.91.74.61.9540.80.4
Prunus serotina Ehrh.0----0----20.40.0
Quercus macrocarpa Michx.1024.215.110.72.10----20.40.0
Quercus marilandica Muenchh.8221.19.310.72.795.23.12.40.9300.20.2
Quercus muehlenbergii Engelm.1827.013.810.32.8112.52.44.41.8570.30.3
Quercus shumardii Buckland231.00.713.1-0----0--
Quercus stellata Wangenh.88822.99.312.03.2234.73.33.32.31570.30.2
Quercus velutina Lam.13731.716.216.34.3162.92.61.80.51610.30.3
Rhus copallina L.0----41.00.02.10.4180.40.3
Robinia pseudoacacia L.0----41.51.03.21.521.00.2
Salix nigra Marshall337.55.812.4-0----0--
Sapindus drummondii L.1415.94.19.31.9352.31.83.61.6310.30.3
Sideroxylon languinosum Michx.4117.86.09.52.592.12.32.81.1820.40.3
Ulmus americana L.15018.810.711.03.5273.52.93.81.44330.30.2
Ulmus rubra Muhl.334.226.520.1-43.02.32.80.260.10.0
Viburnum rufidulum Raf.0----31.71.21.60.130.30.1
Zanthoxylum americanum Mill.0----101.00.02.91.60--
Table 2. Average basal area (BA) and density for overstory trees (dbh >10 cm) and average density for sapling (dbh 0.1 to 10 cm) and seedlings (shorter than 1.37 m) measured in 130 plots in Cross Timbers forests of north-central Oklahoma. Std. = standard deviation. For seedlings, all Celtis spp. are considered Celtis occidentalis.
Table 2. Average basal area (BA) and density for overstory trees (dbh >10 cm) and average density for sapling (dbh 0.1 to 10 cm) and seedlings (shorter than 1.37 m) measured in 130 plots in Cross Timbers forests of north-central Oklahoma. Std. = standard deviation. For seedlings, all Celtis spp. are considered Celtis occidentalis.
Overstory TreesSaplingsSeedlings
BA (m2 ha−1)Stems (ha−1)Stems (ha−1)Stems (ha−1)
Genus SpeciesMeanStd.MeanStd.MeanStd.MeanStd.
Acer negundo0.00.20.43.19.68.10.00.0
Carya illinoinensis0.52.27.125.07.76.961.57.5
Carya texana0.61.917.341.667.323.4476.933.9
Celtis laevigata0.31.19.627.040.418.50.00.0
Celtis occidentalis1.64.131.568.8271.269.92600.0111.6
Cercis canadensis0.00.10.85.3205.8109.4530.840.1
Cornus drummondii0.00.00.00.0325.0140.21215.4164.7
Diospyros virginiana0.10.51.914.115.413.361.57.5
Fraxinus spp.0.62.26.925.311.56.838.56.5
Gleditsia triacanthos0.00.00.22.21.92.423.16.6
Gymnocladus dioicus0.10.81.28.20.00.00.00.0
Juglans nigra0.31.64.415.75.84.20.00.0
Juniperus virginiana1.33.356.5117.2232.7429.1171.4484.0
Morus rubra0.20.75.616.00.00.038.55.7
Platanus occidentalis0.32.02.119.90.00.00.00.0
Prunus angustifolia0.00.00.44.40.00.00.00.0
Prunus mexicana0.00.11.56.050.019.2415.439.4
Prunus serotina0.00.00.00.00.00.015.43.1
Quercus macrocarpa0.10.71.99.70.00.015.44.4
Quercus marilandica0.72.415.848.317.311.0230.817.5
Quercus muehlenbergii0.31.03.510.221.212.9438.536.4
Quercus shumardii0.00.20.43.10.00.00.00.0
Quercus stellata8.28.1170.8174.444.221.41207.771.3
Quercus velutina2.65.326.445.330.811.91238.569.2
Rhus copallina0.00.00.00.07.79.8138.517.5
Robinia pseudoacacia0.00.00.00.07.77.715.43.1
Salix nigra0.10.70.66.60.00.00.00.0
Sapindus drummondii0.10.52.721.667.345.0238.541.6
Sideroxylon lanuginosum0.20.77.921.817.37.0630.829.1
Ulmus americana1.12.728.964.751.920.23330.8198.4
Ulmus rubra0.10.70.63.87.75.946.29.3
Viburnum rufidulum0.00.00.00.05.84.223.13.8
Zanthoxylum americanum0.00.00.00.019.224.40.00.0
Total19.0-407.0-1542.0-13,202.0-
Table 3. Number of times that species were identified on properties (Prop; 25) or plots (130, 5 plots per property with 1 property having 10) for overstory trees (dbh >10 cm), saplings (dbh 0.1 to 10 cm), and seedlings (shorter than 1.37 m), expressed as a percent. Overstory trees were measured in plots 0.04 ha in size. For hardwoods, saplings were measured in plots 0.004 ha in size and seedlings were measured in plots 0.001 ha in size. Juniperus virginiana saplings and seedlings were measured in the 0.04 ha plots. For seedlings, all Celtis spp. are lumped as Celtis occidentalis.
Table 3. Number of times that species were identified on properties (Prop; 25) or plots (130, 5 plots per property with 1 property having 10) for overstory trees (dbh >10 cm), saplings (dbh 0.1 to 10 cm), and seedlings (shorter than 1.37 m), expressed as a percent. Overstory trees were measured in plots 0.04 ha in size. For hardwoods, saplings were measured in plots 0.004 ha in size and seedlings were measured in plots 0.001 ha in size. Juniperus virginiana saplings and seedlings were measured in the 0.04 ha plots. For seedlings, all Celtis spp. are lumped as Celtis occidentalis.
OverstorySaplingSeedling
Genus SpeciesPropPlotPropPlotPropPlot
Acer negundo828200
Carya illinoinensis241082245
Carya texana522240166021
Celtis laevigata602228800
Celtis occidentalis683556309261
Cercis canadensis8256274817
Cornus drummondii0048186018
Diospyros virginiana12382245
Fraxinus spp.40128442
Gleditsia triacanthos414141
Gymnocladus dioicus1220000
Juglans nigra2898200
Juniperus virginiana9653100849246
Morus rubra521400163
Platanus occidentalis1220000
Prunus angustifolia410000
Prunus mexicana28624104816
Prunus serotina000042
Quercus macrocarpa1650041
Quercus marilandica52211244414
Quercus muehlenbergii44121244417
Quercus shumardii820000
Quercus stellata92782068834
Quercus velutina723536106833
Rhus copallina0041246
Robinia pseudoacacia004282
Salix nigra410000
Sapindus drummondii82164165
Sideroxylon lanuginosum60182878834
Ulmus americana8838281210061
Ulmus rubra828282
Viburnum rufidulum0082122
Zanthoxylum americanum004100

Share and Cite

MDPI and ACS Style

Hoff, D.L.; Will, R.E.; Zou, C.B.; Lillie, N.D. Encroachment Dynamics of Juniperus virginiana L. and Mesic Hardwood Species into Cross Timbers Forests of North-Central Oklahoma, USA. Forests 2018, 9, 75. https://doi.org/10.3390/f9020075

AMA Style

Hoff DL, Will RE, Zou CB, Lillie ND. Encroachment Dynamics of Juniperus virginiana L. and Mesic Hardwood Species into Cross Timbers Forests of North-Central Oklahoma, USA. Forests. 2018; 9(2):75. https://doi.org/10.3390/f9020075

Chicago/Turabian Style

Hoff, Daniel L., Rodney E. Will, Chris B. Zou, and Nathan D. Lillie. 2018. "Encroachment Dynamics of Juniperus virginiana L. and Mesic Hardwood Species into Cross Timbers Forests of North-Central Oklahoma, USA" Forests 9, no. 2: 75. https://doi.org/10.3390/f9020075

APA Style

Hoff, D. L., Will, R. E., Zou, C. B., & Lillie, N. D. (2018). Encroachment Dynamics of Juniperus virginiana L. and Mesic Hardwood Species into Cross Timbers Forests of North-Central Oklahoma, USA. Forests, 9(2), 75. https://doi.org/10.3390/f9020075

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