Next Article in Journal
Bacterial Canker Disease on Populus × euramericana Caused by Lonsdalea populi in Serbia
Previous Article in Journal
Fine-Scale Microclimate Pattern in Forest-Steppe Habitat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reconstruction of Piñon–Juniper Woodlands in the Sky Islands of the Davis Mountains, Texas, USA

1
College of Forestry, Agriculture and Natural Resources, University of Arkansas at Monticello, Monticello, AR 71655, USA
2
Arthur Temple College of Forestry and Agriculture, Stephen F. Austin State University, Nacogdoches, TX 75965, USA
*
Author to whom correspondence should be addressed.
Forests 2020, 11(10), 1079; https://doi.org/10.3390/f11101079
Submission received: 2 September 2020 / Revised: 1 October 2020 / Accepted: 5 October 2020 / Published: 9 October 2020
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Piñon (Pinus spp. L.)–juniper (Juniperus spp. L.) woodlands’ historical stand structures were recreated to provide reference conditions and document long-term changes in the Sky Islands of the Davis Mountains, Texas. Restoration of these isolated woodlands requires insights into the range of variability in current and historical stand structures, as well as an understanding of the spatiotemporal establishment and recruitment patterns of tree species. With drastic changes in forests and woodlands of the Southwestern United States widely reported, the main objective of this study was to reconstruct woodland tree temporal and spatial establishment patterns. A stratified random sampling approach was used to select two study sites each of 3600 m2 in area. Within each site, all individual trees were mapped, measured, and cored for age determination. Age and tree location data were used to recreate the spatiotemporal patterns of tree species establishment and recruitment. Increments in density of both Mexican piñon (Pinus cembroides var. cembroides Zucc.) and alligator juniper (Juniperus deppeana var. deppeana Steud.) reached 422 trees ha−1 in the 115-year period between 1890 and 2005; a yearly increment of 4 trees ha−1 that did not reflect a rapid rate of change in these piñon–juniper woodlands. Age distributions reflected the multi-cohort nature of these woodlands, and spatial autocorrelation measures were useful in the objective delineation of these cohorts. Temporal and functional niche differentiation of juniper was reflected in the development pattern where alligator juniper served as a pioneer species, exhibited a longer period of substantial recruitment, and had greater recruitment rates than that of Mexican piñon. Recruitment of Mexican piñon and alligator juniper occurred in an episodic fashion, with the majority of recruits being acquired between 1890 and 1949.

1. Introduction

Reconstructive studies enable the assessment of long-term changes and ecosystem dynamics, aid in the interpretation of ecosystem disturbance history, and are essential in guiding conservation and management goals [1,2]. Examination of changes in forest and woodland stand structure over time can be accomplished through the use of permanent plots that are monitored for tree recruitment, growth, and mortality at regular intervals. However, this approach may be limited in application due to inherent costs, as well as changes in land ownership. An alternative approach is to determine age structure retrospectively using dead and live trees [3]. The use of static age structures (at one point in time) in the interpretation of stand dynamics has been scrutinized, however. The absence of trees in past age classes has been posited to be a direct result of high rates of mortality and decomposition in lodgepole pine (Pinus contorta Dougl. ex Loud.)-Engelmann spruce (Picea engelmannii Parry ex Engelm.) stands, rather than lack of recruitment, as is commonly interpreted using static age distributions [4]. To circumvent the effect of a presumably high tree decomposition at early stages of stand development [5], it was recommended to construct age distributions at multiple points in time using fixed intervals [4]. This notion has been challenged, however, and static age distributions of ponderosa pine (Pinus ponderosa Lawson & C. Lawson) were deemed a faithful reflection of actual long-term establishment patterns rather than missing trees that died and decomposed with no trace, especially when standing dead and downed trees were included in these studies [6]. Evaluation of stand dynamics of a Mexican piñon pine (Pinus cembroides var. cembroides Zucc.) forest at the southern edge of its distributional range [7] found a preponderance of small trees (diameter <20 cm and height <6 m). The presence of two major cohorts and the lack of trees in past age classes suggested that recruitment and regeneration patterns were in discrete intervals rather than a continuous manner.
Sky islands are defined as an inland terrain that is made up of a widely spaced sequence of valleys and mountains [8], supporting a diversity of ecosystems and species within a few kilometers and relatively small changes in elevation. These isolated ecosystems, including the piñon–juniper (Pinus spp. L. Juniperus spp. L.) woodlands of the Madrean Archipelago, may be the least studied of all ecosystems in North America [9]. The Madrean Archipelago is characterized by high biological diversity due to the presence of two major floristic realms, the Rocky Mountains and the Sierra Madre Occidental, in addition to the convergence of three major climatic zones (temperate, subtropical, and tropical) [8]. The forest and woodland vegetation of the Madrean Archipelago has been classified [10] into three types: (1) encinal woodlands, (2) piñon–juniper woodlands and (3) high elevation coniferous forests.
There is a growing interest in restoring the forests and woodlands of the Southwestern United States [6,11,12,13,14,15,16]. Restoration efforts will be an increasing challenge with climate change [17,18,19,20]. While Tausch et al. [20] and Miller et al. [21] reported an increase in tree density, accompanied by a decline in herbaceous production, as well as vertical expansion of piñon–juniper communities within the Great Basin, Blackburn and Tueller [22] noted the expansion of piñon and juniper into black sagebrush (Artemisia nova) communities in East-Central Nevada, and Landis and Bailey [15] reported substantial increments in density, basal area, and canopy cover of piñon–juniper communities since 1860. These changes in piñon–juniper communities were mainly attributed to overgrazing, fire exclusion, changing climate, and logging [21,22,23,24]. Restoration of these ecosystems requires the determination of reference conditions that can be used as approximations in guiding conservation plans and restoration treatments and take into account keystone variables capable of reflecting the evolutionary environment of a species or group of species [23,25]. These reference conditions can be defined using the range of historic variability in ecological structures and processes, which reflects the recent evolutionary history of these ecosystems [26]. Unfortunately, some restoration efforts were the product of a misconception in which the entire piñon–juniper type was considered “degraded” [27]. Although changes in structure and distribution have occurred throughout the piñon–juniper’s extensive distribution range [15,20,22], some areas did not witness substantial increments in tree density or spatial extent [27]. For example, Ffolliott and Gottfried [28] reported an increase of approximately 148 trees ha−1 in 52 years (between 1938 and 1991) near Flagstaff, Arizona.
This study was conducted to provide reference conditions for restoration treatments and document long-term changes in stand structures of piñon–juniper woodlands of the Davis Mountains Sky Islands. The main objective of this study was to recreate historical stand structures by reconstructing tree temporal and spatial establishment patterns at five-points in time starting from the year 2005 to the earliest establishment date (i.e., 1985, 1965, 1945, 1925, and 1905).

2. Materials and Methods

The study was located within The Texas Nature Conservancy (TNC) Davis Mountains Preserve in Jeff Davis County, Texas. Climate is cool temperate with warm summers. Average annual rainfall is approximately 406 mm, with summer monsoons comprising two-thirds of the annual rainfall. Mean summer temperature is 27 °C, whereas winter average temperatures range between −1 and 10 °C. Soils of the study area belong mainly to the Puerta–Madrone association (Alfic Lithic Argiustolls), with Loghouse association soils (Udic Haplustalfs), along stream channels of the study area [29].
Overstory vegetation includes: ponderosa pine, Mexican piñon, Emory oak (Quercus emoryi Torr.), gray oak (Quercus grisea Libem.), silverleaf oak (Quercus hypoleucoides A. Camus), and alligator juniper (Juniperus deppeana var. deppeana Steud.). Common understory herbaceous vegetation includes sideoats grama (Bouteloua curtipendula (Michx.) Torr.), bulb panicgrass (Panicum bulbosum Kunth), pinyon ricegrass (Piptochaetium fimbriatum (Kunth) Hitchc.), and longleaf cologania (Cologania angustifolia Kunth). Catclaw mimosa (Mimosa aculeaticarpa Ortgea) is often abundant in the understory layer [16].
The TNC Davis Mountains Preserve was stratified subjectively but without preconceived bias [30] into homogeneous areas based on vegetation, slope, elevation, and aspect. Stratification was achieved by merging 1 m Digital Orthophoto Quadrangles and 30 m Digital Elevation Models using a suite of tools in the ArcToolbox of ArcGIS 9.0 (ESRI, Redlands, CA, USA). To ensure that representative piñon–juniper communities were selected, areas less than 50 ha in size were excluded. From this sampling frame, two study sites (60 × 60 m) of similar vegetation and topography were randomly selected.
For each site, all overstory (height ≥ 1.4 m) vegetation was identified to the species level, and diameter at breast height (1.37 m; dbh) measured to the nearest 0.1 cm. All trees greater than 15.0 cm in dbh were cored at the base using an increment borer, and root collar diameter (RCD) and total height were measured to the nearest 0.1 cm and 0.1 m, respectively. Due to the high density of small-diameter (dbh ≤ 15.0 cm) trees, a sub-sample consisting of 15 trees (1 tree per diameter class) per species was randomly selected for age determination. Within this sub-sample, trees >5.0 cm dbh were cored at the base, whereas trees ≤5.0 cm in dbh were cut and cross-sections collected. Dead, standing and downed, trees as well as stumps were included, regardless of decay stage, in the sample to extend age structure and minimize bias toward younger trees [4,6]. Dead wood accounted for 2% (8 trees ha−1) and 7% (67 trees ha−1) for site 1 and site 2, respectively. All trees within a site were positioned using a global positioning system receiver (Trimble GeoXT®, Sunnyvale, CA, USA). For trees along the perimeter of the site (perimeter trees), distance to the nearest neighbor was measured to the nearest 0.1 m when the nearest neighbor was outside the perimeter of the plot [31]. For multi-stemmed trees, only the largest diameter (dbh) stem was considered for sampling. All seedlings (0.1 m ≤ height < 1.4 m) within each site were tallied and identified to the species level.
Tree cores were air dried, mounted, and sanded using progressively finer sand paper (50, 240, and 320 grit) following standard methodology [32]; cross-sections were processed in the same manner. Cores from the same species and same site were cross-dated (matching ring pattern) to ensure identification of missing and false rings by graphing ring-width measurements against the number of years [32]. For dead trees and stumps, death or cut dates were determined as the date of the outermost ring [33,34]. Decay was not sufficiently severe to prevent direct measurement and estimation of most death dates. Annual growth increments (ring-width) were measured to the nearest 0.001 mm using a Velmex TA system (Velmex Inc., Bloomfield, NY, USA). The quality of cross-dating was checked using COFECHA software [35,36].
Dbh and age frequency distributions were constructed and tree density (trees ha−1) and basal area (m2 ha−1) calculated. Age structure (frequency distribution of tree age data) was constructed at multiple points in time using 20-year intervals. Estimates of tree density (trees ha−1) at multiple points in time were derived as the number of established trees that were still alive at that point in time.
To examine contemporary and historical tree spatial pattern, distance to the nearest neighbor for each tree within each site was calculated using an Arc Macro Language (AML) program in ArcInfo Workstation (ESRI, Redlands, CA, USA). Nearest neighbor distances were replaced with field measurements for perimeter trees, whose nearest neighbor lay outside the site. Nearest neighbor distances were pooled to calculate the observed mean nearest neighbor distance for each site [31]. The observed and expected mean nearest neighbor distance as well as the R statistic (the ratio of the observed mean nearest neighbor distance to the expected mean nearest neighbor distance) and the standard normal deviate (Z score) were calculated using Average Nearest Neighbor in ArcToolbox of ArcGIS 9.0 (ESRI, Redlands, CA, USA). The null hypothesis “the observed tree spatial pattern is random across the study site” was tested using an alpha level of 0.05 (−1.96 ≤ ZCritical ≤ 1.96).
For each site, a total of 20 adjacency matrices (adjacent tree pairs represented by 1 and non-adjacent tree pairs represented by 0) were constructed. Adjacency rule differed for each adjacency matrix and was defined in 3 m distance classes [37]. Adjacent tree pairs were defined as those pairs that are 0.1–3.0, 3.1–6.0, 6.1–9.0, and 57.1–60.0 m apart (one adjacency rule for each adjacency matrix). Adjacency matrices were used as the weighting matrices in calculating Moran’s I spatial autocorrelation statistic for each distance class (d) for the variable tree age. Moran’s I values were calculated using Spatial Autocorrelation in ArcToolbox of ArcGIS 9.0 (ESRI, Redlands, CA, USA). To test the null hypothesis “the observed spatial distribution of tree ages for each d is random”, the standard normal deviate of that distance class (Z (d)) was compared to the critical standard normal deviate (−1.96 ≤ ZCritical ≤ 1.96) using an alpha level of 0.05. Correlograms were constructed by plotting the values of Z (d) on the vertical axis Y-axis) against d on the horizontal axis (X-axis) to explore the presence of even and uneven aged tree patches [37]. To identify even and uneven aged tree patches, Anselin Local Moran’s I was calculated using Cluster and Outlier Analysis in ArcToolbox of ArcGIS 9.0 (ESRI, Redlands, CA, USA).

3. Results

Diameter distributions for all trees within each site were positively skewed, indicating an abundance (92%) of small-diameter (dbh < 20.0 cm) trees. Tree density and basal area were 711 trees/ha and 12 m2 ha−1, respectively. Alligator juniper was the most dominant species with mean basal area of 6 m2 ha−1 and tree density of 283 trees/ha. Gray oak and Mexican piñon were ranked second in abundance with mean basal areas of 3 m2 ha−1 and tree densities of about 200 trees/ha. Emory oak was a minor component but dominated the seedling strata.
Combined age distribution for piñon and juniper was unimodal, with a substantial establishment occurring from 1925 to 1945 (Figure 1). Except for the period 1900–1929, recruitment rates of alligator juniper were always greater than piñon, with both species showing substantial increases in tree densities over time. Alligator juniper density increased from 33 trees/ha in the 1890s to 274 trees/ha in 2005, whereas Mexican piñon density increased from 7 trees/ha in the 1890s to 188 trees/ha in 2005 (Figure 1).
Contemporary spatial patterns of Mexican piñon, alligator juniper, oak, and all trees combined were clustered and significantly deviated from a random distribution. Historical and contemporary piñon–juniper joint patterns were also clustered (Table 1; Figure 2 and Figure 3). Mexican piñon spatial patterns were significantly clustered, except in 1905, where the pattern was random (Table 1); alligator juniper spatial patterns were always significantly clustered.
Positive spatial autocorrelation for Mexican piñon tree age was shown for 2005, 1985, and 1965 for both sites (Table 2). Positive autocorrelation indicated that tree clusters identified using nearest neighbor distances were of relatively similar age. Although not statistically significant, age spatial pattern had a negative autocorrelation for the periods 1925 for site 1 and 1905 for site 2, indicating that early establishment had trees with relatively contrasting ages during these times. Age data did not extend beyond 1905 for either site or species. Alligator juniper age spatial autocorrelation had a similar pattern to that of Mexican piñon with significant positive autocorrelation for 2005, 1985, and 1965 (Table 2). No significant spatial autocorrelation was revealed for alligator juniper age patterns at any time period for site 2.
Correlograms, using adjacency matrices to define tree neighborhood, showed significant positive spatial autocorrelation at distances of 0.1–3.0, 6.1–9.0, and 18.1–21.0 m for 2005 and 1985 for Mexican piñon tree ages at site 1 and significant negative spatial autocorrelation in the same years at distances of 36.1–39.0 and 42.1–45.0 m (Figure 4). In 1965, significant positive spatial autocorrelation was revealed at distances of 0.1–3.0 and 18.1–21.0 m, whereas negative spatial autocorrelation was revealed at distances of 36.1–39.0 and 42.1–45.0 m. In 1945, significant positive autocorrelation was revealed at a distance of 21.1–24.0 m. Within site 2, significant positive spatial autocorrelation was revealed for Mexican piñon age data in 2005 at a distance of 0.1–3.0 m (Figure 5). In 1905, significant positive spatial autocorrelation was revealed at a distance of 24.1–27.0 m; however, only three neighbors were used in the calculation of Moran’s I for this distance class.
Alligator juniper correlograms showed a significant positive spatial autocorrelation in 1985 and 2005 at distances of 0.1–3.0 and 57.1–60.0 m for site 1; eleven and nine neighbors were used in the calculation of Moran’s I for these distance classes, respectively (data not shown). In 1985, spatial autocorrelation for alligator juniper age data was similar to that of 2005. In 1965, the first distance class (0.1–3.0 m) had significant positive spatial autocorrelation, but only 10 neighbors were used in Moran’s I calculation. Negative autocorrelation (14 neighbors) was shown for 1945 at a distance of 15.1–18.0 m. In 1925, significant positive (at 42.1–45.0 m) and negative (at 45.1–48.0 and 51.1–54.0 m) spatial autocorrelation were revealed, and in 1905, significant positive (at 36.1–39.0 and 42.1–45.0 m) spatial autocorrelation was revealed. Within site 2, alligator juniper age data of 2005 and 1985 were significantly positively autocorrelated only at 57.1–60.0 m (data not shown). In 1925, significant negative spatial autocorrelation was revealed at distances of 3.1–6.0 and 54.1–57.0 m. In 1905, significant positive spatial autocorrelation was revealed at a distance of 51.1–54.0 m. Mapping of Anselin Local Moran’s I values indicated the presence of local clusters of relatively similar age trees for 2005, 1985, and 1965 (Figure 6 and Figure 7).

4. Discussion

Using diameter as a surrogate for age, reverse “J” shape diameter distributions are considered an indication of uneven-age, all-age, or multi-cohort stands [5]. Hump-like diameter distribution patterns may indicate the presence of at least two cohorts: an old cohort represented by large diameter trees and a young cohort represented by small-diameter trees.
Contemporary age structure indicated the presence of at least two cohorts: a remnant or initial alligator juniper cohort and a mixed piñon–juniper secondary cohort. The descriptors remnant and initial refer to two possible scenarios that may describe the development of this cohort. The initial scenario suggests the encroachment of alligator juniper into an open area, which may be attributed to comparatively greater seed longevity [38], greater tolerance of photosynthetic apparatus to water stress [39], greater drought tolerance, and better use of shallow soil moisture than piñon pine [40]. These pioneer members may have served as nurse trees for subsequent recruitment of piñon and juniper and may have also served as seed sources [38]. The remnant scenario, on the other hand, suggests opening of the stand by means of disturbance (e.g., firewood harvesting and/or fire) prior to 1890 that resulted in survival of this cohort. This scenario is consistent with the establishment of nearby Fort Davis in 1854. The secondary cohort corresponded to substantial increments in piñon–juniper recruitment over several decades (1890–1949). Substantial recruitment rates of Mexican piñon coincided with that of alligator juniper, although alligator juniper recruitment rates extended a little longer (1890–1969). The idea that junipers persist in the seed banks and contribute a continuous input into the seed bank supports this notion [38]. The period of substantial recruitment rates could be a result of favorable establishment and survival conditions and/or a result of response to disturbance in the form of fire exclusion and/or overgrazing [20,21,22,23,24]. A portion (1935–1949) of the period of substantial recruitment rates for piñon–juniper exhibited an above average annual precipitation, and the period was preceded by the establishment of Fort Davis and the Euro-American settlement of the area. This establishment period coincides with establishment periods for piñon–juniper in Northern Arizona [15] and the Great Basin [20].
Piñon–juniper age distributions between 1905 and 1945 suggest a stand initiation stage which occurred over several decades. The distribution of Mexican piñon indicated less variation in age as compared to alligator juniper, and this may be attributed to shorter recruitment periods of Mexican piñon. In addition, age distributions suggested a longer lifespan for alligator juniper, which may also have served as a pioneer species in structural development.
Greater recruitment rates and density estimates of alligator juniper between 1890 and 2005 could be explained by the same factors potentially responsible for the establishment of junipers in canopy interspaces and open areas [38,39,40].
Significant spatial clustering may be explained by seed dispersal mechanisms in which piñons are mainly dispersed by birds (e.g., piñon jays (Gymnorhinus cyanocephalus)) and rodents (e.g., piñon mice (Peromyscus truei)),whereas junipers are mainly dispersed by birds (e.g., bluebirds (Sialia mexicana and S. currucoides)), rodents, as well as mammals (e.g., mule deer (Odocoileus hemionus) and coyote (Canis latrans)), after passing through the digestive system [38]. The spatial clustering may also have created habitats for this variety of animals, further enhancing seed dispersal by these species. Animal-dependent seed dispersal (scattering in multi-seed caches), seed dispersal by gravitational forces, and the spatial heterogeneity of the environment may have contributed to the observed spatial pattern.
The observed spatial pattern may be confounded by the inability to capture spatial arrangement of seedlings that might have occurred in great numbers and experienced high rates of mortality and decay [5]. Temporal spatial patterns indicated that Mexican piñon and alligator juniper were equally successful in inter-canopy and open area establishment as well as establishment within existent clusters, in contrast to Everett et al. [41] who found that piñon seedlings were more abundant under tree canopies than in inter-canopy spaces in Utah. In the Davis Mountains, however, summer monsoons may have lessened the need for nurse-tree micro-environments that are necessary for establishment and survival in the Great Basin [38]. In addition, differences among species (Mexican piñon and alligator juniper versus piñon pine and Utah juniper) in terms of light, moisture, competition tolerance, and nutrient levels needed for germination and seedling survival and development may have attributed to the observed differences. Climate change would certainly impact these development patterns, and under future climate projections, nurse-tree micro-environments may become more important for establishment and survival of piñon and juniper in the Davis Mountains [17,18,19].
Significant positive spatial autocorrelation of Mexican piñon in 2005, 1985, and 1965 indicated the presence of piñon clusters of relatively similar age; this similarity may be attributed to the lack of piñon density variation among these periods. In 1945, 1925, and 1905, the spatial distribution of ages was similar to that expected from an independent random process. The presence of clusters of relatively similar age piñons and clusters of relatively dissimilar age piñons reflected the multi-cohort nature of Mexican piñons.
Under the assumption that juniper regeneration occurred in a clustered pattern, differences in the spatial distribution of successful recruitment between 1925 and 1945 may be attributed to higher rates of seedling mortality (e.g., low intensity surface fire or a more pronounced within-cluster competition) prior to 1925. In general, the shift from negative to positive spatial autocorrelation emphasized the clustered recruitment patterns of piñon and juniper, and contemporary spatial distributions of piñon–juniper ages reflected the irregular (different size clusters) multi-cohort nature of these stands.

5. Conclusions

Stand structure reconstruction in this piñon–juniper woodland revealed substantial increments in density, similar to those reported in the Great Basin and Northern Arizona [15,20]. However, when viewed in the context of the elapsed time, annual increments were similar to those reported by Ffolliott and Gottfried [28] and did not reflect a rapid rate of change in these piñon–juniper communities. Age distributions reflected the multi-cohort nature of these communities. Temporal and functional niche differentiation of alligator juniper was reflected in the development pattern where alligator juniper served as a pioneer species, exhibited a longer period of substantial recruitment, and had greater recruitment rates than that of Mexican piñon, which may have been influenced by climate conditions that favored germination and establishment of woody species in this region. Recruitment of Mexican piñon and alligator juniper occurred in an episodic fashion with the majority of recruits being acquired between 1890 and 1949. Examination of spatial pattern and spatial autocorrelation measures suggested that recruitment occurred in different-size clusters within the canopies of existent clusters as well as among inter-canopy spaces.
Restoration efforts should consider other lines of evidence (e.g., historical photographs and early historical accounts), in conjunction with stand reconstruction data, when establishing reference or desired conditions for restoration treatments. Although stand structure reconstruction may have provided insight into the development pattern of these piñon–juniper woodlands, an understanding of the development pattern is far from being complete without consideration of its contextual past climate and disturbance history [42]. When stand structure reconstruction is coupled with the evaluation of disturbance history, especially the reconstruction of fire frequency, and tree species distribution climate change projections a better understanding of the development patterns and restoration of resilient woodlands would be attained.

Author Contributions

Conceptualizations, M.M.B. and B.P.O.; methodology, M.M.B., B.P.O., H.M.W., D.R.U., and I.-K.H.; formal analysis, M.M.B., D.R.U., and I.-K.H.; investigation, M.M.B.; data curation, M.M.B. and B.P.O.; writing—original draft preparation, M.M.B. and B.P.O.; writing—review and editing, H.M.W., D.R.U., and I.-K.H.; supervision, B.P.O.; project administration and funds acquisition, B.P.O. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for this research was provided to the Arthur Temple College of Forestry and Agriculture at Stephen F. Austin State University by the USFS Rocky Mountain Research Station (Grant 00-JV-11221615-108) and the McIntire-Stennis Cooperative Forest Research Program (Program TEXY158).

Acknowledgments

The authors also acknowledge the support of The Nature Conservancy for access to the Davis Mountain Preserve and for local logistical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analysis, interpretation of data, the writing of the manuscript, or the decision to publish the results.

References

  1. Foster, D.R.; Otwig, D.A.; McLachlan, J.S. Ecological and conservation insights from. Trends Ecol. Evol. 1996, 11, 419–424. [Google Scholar] [CrossRef]
  2. Swetnam, T.W.; Allen, C.D.; Betancourt, J.L. Applied historical ecology: Using the past to manage for the future. Ecol. Appl. 1999, 9, 1189–1206. [Google Scholar] [CrossRef]
  3. Clark, S.L.; Hallgren, S.W. Age estimation of Quercus marilandica and Quercus stellata: Applications for interpreting stand dynamics. Can. J. For. Res. 2004, 34, 1353–1358. [Google Scholar] [CrossRef]
  4. Johnson, E.A.; Miyanishi, K.; Kleb, H. The hazards of interpretation of static age structures as shown by stand reconstructions in a Pinus contorta—Picea engelmannii forest. J. Ecol. 1994, 82, 923–931. [Google Scholar] [CrossRef]
  5. Oliver, C.D.; Larson, B.C. Forest Stand Dynamics; McGraw-Hill, Inc.: New York, NY, USA, 1996; ISBN 0471138339. [Google Scholar]
  6. Mast, J.N.; Fule, P.Z.; Moore, M.M.; Covington, W.W.; Waltz, A.E.M. Restoration of presettlement age structure of an Arizona Ponderosa pine forest. Ecol. Appl. 1999, 9, 228–239. [Google Scholar] [CrossRef]
  7. Segura, G.; Snook, L.C. Stand dynamics and regeneration patterns of a pinyon pine forest in east central Mexico. For. Ecol. Manag. 1992, 47, 175–194. [Google Scholar] [CrossRef]
  8. Warshall, P. The Madrean Sky Island Archipelago: A planetary overview. In Proceedings of the Biodiversity and Management of the Madrean Archipelago: The Sky Islands of Southwestern United States and Northwestern Mexico, Tucson, AZ, USA, 19–23 September 1994; DeBano, L.F., Ffolliott, P.F., Ortega-Rubio, A., Gottfried, G.J., Hamre, R.H., Edminster, C.B., Eds.; USDA, Forest Service, Rocky Mountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1995; pp. 6–18. [Google Scholar]
  9. DeBano, L.F.; Ffolliott, P.F. The sky island conference: Looking back, looking ahead. In Proceedings of the Biodiversity and Management of the Madrean Archipelago: The Sky Islands of Southwestern United States and Northwestern Mexico, Tucson, AZ, USA, 19–23 September 1994; DeBano, L.F., Ffolliott, P.F., Ortega-Rubio, A., Gottfried, G.J., Hamre, R.H., Edminster, C.B., Eds.; USDA, Forest Service, Rocky Mountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1995; pp. 1–5. [Google Scholar]
  10. Gottfried, G.J.; Ffolliott, P.F.; DeBano, L.F. Forests and woodlands of the sky islands: Stand characteristics and silvicultural prescriptions. In Proceedings of the Biodiversity and Management of the Madrean Archipelago: The Sky Islands of Southwestern United States and Northwestern Mexico, Tucson, AZ, USA, 19–23 September 1994; DeBano, L.F., Ffolliott, P.F., Ortega-Rubio, A., Gottfried, G.J., Hamre, R.H., Edminster, C.B., Eds.; USDA, Forest Service, Rocky Mountain Forest and Range Experiment Station: Fort Collins, CO, USA, 1995; pp. 152–164. [Google Scholar]
  11. Howell, J.J. Pinon and juniper woodlands of the southwest. J. For. 1941, 39, 542–545. [Google Scholar] [CrossRef]
  12. Cooper, C.F. Changes in vegetation structure and growth of southwestern Ponderosa pine forests since white settlement. Ecol. Monogr. 1960, 30, 129–164. [Google Scholar] [CrossRef]
  13. Covington, W.W.; Moore, W.W. Postsettlement changes in natural fire regimes and forest structure: Ecological restoration of old-growth ponderosa pine forests. J. Sustain. For. 1994, 2, 153–181. [Google Scholar] [CrossRef]
  14. Moore, M.M.; Huffman, D.W.; Fulé, P.Z.; Covington, W.W.; Crouse, J.E. Comparison of historical and contemporary forest structure and composition on permanent plots in southwestern ponderosa pine forests. For. Sci. 2004, 50, 162–176. [Google Scholar] [CrossRef]
  15. Landis, A.G.; Bailey, J.D. Reconstruction of age structure and spatial arrangement of piñon-juniper woodlands and savannas of Anderson Mesa, northern Arizona. For. Ecol. Manag. 2005, 204, 221–236. [Google Scholar] [CrossRef]
  16. Bataineh, M.M.; Oswald, B.P.; Bataineh, A.L.; Farrish, K.W.; Coble, D.W.; Edminster, C.B. Plant communities associated with Pinus ponderosa forests in the sky islands of the Davis Mountains, Texas. J. Torrey Bot. Soc. 2007, 134, 468–478. [Google Scholar] [CrossRef]
  17. Friggens, M.M.; Warwell, M.V.; Chambers, J.C.; Kitchen, S.G. Modeling and predicting vegetation response of western USA grasslands, shrublands, and deserts to climate change. In Climate Change in Grasslands, Shrublands, and Deserts of the Interior American West: A Review and Needs Assessment; Gen. Tech. Rep. RMRS-GTR-285; Finch, D.M., Ed.; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2012; pp. 1–20. [Google Scholar]
  18. Rehfeldt, G.E.; Crookston, N.L.; Sáenz-Romero, C.; Campbell, E.M. North American vegetation model for land-use planning in a changing climate: A solution to large classification problems. Ecol. Appl. 2012, 22, 119–141. [Google Scholar] [CrossRef] [PubMed]
  19. Thorne, J.H.; Choe, H.; Stine, P.A.; Chambers, J.C.; Holguin, A.; Kerr, A.C.; Schwartz, M.W. Climate change vulnerability assessment of forests in the Southwest USA. Clim. Chang. 2018, 148, 387–402. [Google Scholar] [CrossRef] [Green Version]
  20. Tausch, R.J.; West, N.E.; Nabi, A.A. Tree age and dominance patterns in Great Basin pinyon-juniper woodlands. J. Range Manag. 1981, 34, 259–264. [Google Scholar] [CrossRef]
  21. Miller, R.F.; Tausch, R.J.; MacArthur, D.; Johnson, D.D.; Sanderson, S.C. Age Structure and Expansion of Piñon-Juniper Woodlands: A Regional Perspective in the Intermountain West; Res. Pap. RMRS-RP-69; Rocky Mountain Research Station: Fort Collins, CO, USA, 2008.
  22. Blackburn, W.H.; Tueller, P.T. Pinyon and juniper invasion in black sagebrush communities in east-central Nevada. Ecology 1970, 51, 841–848. [Google Scholar] [CrossRef]
  23. White, P.S.; Walker, J.L. Approximating nature’s variation: Selecting and using reference information in restoration ecology. Restor. Ecol. 1997, 5, 338–349. [Google Scholar] [CrossRef] [Green Version]
  24. Miller, R.F.; Tausch, R.J. The role of fire in juniper and piñon woodlands: A descriptive analysis. In Proceedings of the Invasive Species Workshop: The Role of Fire in the Control and Spread of Invasive Species—Fire Conference 2000: The First National Congress on Fire Ecology, Prevention, and Management, San Diego, CA, USA, 18 November–1 December 2000; Galley, K.E.M., Wilson, T.P., Eds.; Tall Timbers Research Station: Tallahassee, FL, USA, 2001; pp. 15–30. [Google Scholar]
  25. Moore, M.M.; Covington, W.W.; Fule, P.Z. Reference conditions and ecological restoration: A southwestern ponderosa pine perspective. Ecol. Appl. 1999, 9, 1266–1277. [Google Scholar] [CrossRef]
  26. Fulé, P.Z.; Covington, W.W.; Moore, M.M. Determining reference conditions for ecosystem management of southwestern ponderosa pine forests. Ecol. Appl. 1997, 7, 895–908. [Google Scholar] [CrossRef]
  27. Romme, W.H.; Floyd-Hanna, L.; Hanna, D.D. Ancient pinon–juniper forests of Mesa Verde and the West: A cautionary note for forest restoration programs. In Proceedings of the Conference on Fire, Fuel Teatments, and Ecological Restoration, Fort Collins, CO, USA, 16–18 April 2002; Omni, P.N., Joyce, L.A., Eds.; USDA Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2003; pp. 335–350. [Google Scholar]
  28. Ffolliott, P.F.; Gottfried, G.J. Dynamics of a Pinyon-Juniper Stand in Northern Arizona: A Half-Century History; Res. Pap. RMRS-RP-35; Rocky Mountain Research Station: Fort Collins, CO, USA, 2002.
  29. Turner, A.J. Soil Survey of Jeff Davis County, Texas; USDA, Soil Conservation Service: College Station, TX, USA, 1977. [Google Scholar]
  30. Mueller-Dombois, D.; Ellenberg, H. Aims and Methods of Vegetation Ecology; John Wiley and Sons, Inc.: New York, NY, USA, 1974. [Google Scholar]
  31. Clark, P.J.; Evans, F.C. Distance to nearest neighbor as a measure of spatial relationships in populations. Ecology 2016, 35, 445–453. [Google Scholar] [CrossRef]
  32. Stokes, M.A.; Smiley, T.L. An Introduction to Tree-Ring Dating; The University of Arizona Press: Tucson, AZ, USA, 1996. [Google Scholar]
  33. Bataineh, M.M.; Daniels, L.D. An objective classification of largewood in streams. For. Ecol. Manag. 2014, 313, 1–9. [Google Scholar] [CrossRef]
  34. Daniels, L.D.; Dobry, J.; Klinka, K.; Feller, M.C. Determining year of death of logs and snags of Thuja plicata in southwestern coastal British Columbia. Can. J. For. Res. 1997, 27, 1132–1141. [Google Scholar] [CrossRef]
  35. Grissino-Mayer, H.D. Evaluating crossdating accuracy: A manual and tutorial for the computer program COFECHA. Tree-Ring Res. 2001, 57, 205–221. [Google Scholar]
  36. Holmes, R.L. Computer-assisted quality control in tree-ring dating and measurement. Tree-Ring Bull. 1983, 43, 69–78. [Google Scholar]
  37. Duncan, R.P.; Stewart, G.H. The temporal and spatial analysis of tree age distributions. Can. J. For. Res. 1991, 21, 1703–1710. [Google Scholar] [CrossRef]
  38. Chambers, J.C.; Schupp, E.W.; Vander Wall, S.B. Seed dispersal and seedling establishment of piñon and juniper species within the piñon-juniper woodland. In Proceedings of the Ecology and Management of Piñon-Juniper Communities Within the Interior West, Provo, UT, USA, 15–18 September 1997; Monsen, S.B., Stevens, R., Eds.; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Ogden, UT, USA, 1999; pp. 29–34. [Google Scholar]
  39. Nowak, R.S.; Moore, D.J.; Tausch, R.J. Ecophysiological patterns of piñon and juniper. In Proceedings of the Ecology and Management of Piñon-Juniper Communities Within the Interior West, Provo, UT, USA, 15–18 September 1997; Monsen, S.B., Stevens, R., Eds.; U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station: Ogden, UT, USA, 1999; pp. 35–46. [Google Scholar]
  40. Breshears, D.D.; Myers, O.B.; Johnson, S.R.; Meyer, C.W.; Martens, S.N. Differential use of spatially heterogeneous soil moisture by two semiarid woody species: Pinus edulis and Juniperus Monosperma. J. Ecol. 1997, 85, 289–299. [Google Scholar] [CrossRef]
  41. Everett, R.L.; Koniak, S.; Budy, J.D. Pinyon Seedling Distribution among Soil Surface Microsites; US Department of Agriculture, Forest Service, Intermountain Research Station: Ogden, UT, USA, 1986; Volume INT-363.
  42. Harris, J.A.; Hobbs, R.J.; Higgs, E.; Aronson, J. Ecological restoration and global climate change. Restor. Ecol. 2006, 14, 170–176. [Google Scholar] [CrossRef]
Figure 1. Establishment rates and density estimates per decade for piñon and juniper from two study sites within the Nature Conservancy Davis Mountains Preserve, Jefferson County, Texas.
Figure 1. Establishment rates and density estimates per decade for piñon and juniper from two study sites within the Nature Conservancy Davis Mountains Preserve, Jefferson County, Texas.
Forests 11 01079 g001
Figure 2. Spatial pattern of Mexican piñon (Pinus cembroides var. cembroides Zucc.), alligator juniper (Juniperus deppeana var. deppeana Steud.), and oak species (Quercus spp. L.) within study site 1 in 1905 through 2005.
Figure 2. Spatial pattern of Mexican piñon (Pinus cembroides var. cembroides Zucc.), alligator juniper (Juniperus deppeana var. deppeana Steud.), and oak species (Quercus spp. L.) within study site 1 in 1905 through 2005.
Forests 11 01079 g002
Figure 3. Spatial pattern of Mexican piñon (Pinus cembroides var. cembroides Zucc.), alligator juniper (Juniperus deppeana var. deppeana Steud.), and oak species (Quercus spp. L.) within study site 2 in 1905 through 2005.
Figure 3. Spatial pattern of Mexican piñon (Pinus cembroides var. cembroides Zucc.), alligator juniper (Juniperus deppeana var. deppeana Steud.), and oak species (Quercus spp. L.) within study site 2 in 1905 through 2005.
Forests 11 01079 g003
Figure 4. Contemporary and historical spatial autocorrelation correlograms for Mexican piñon (Pinus cembroides var. cembroides Zucc.) of study site 1. Dotted lines in the middle of each graph indicate the 0.05 significance level, and stars (★) indicate fewer than 20 neighbors for the specified distance class. The 1905 correlogram is left intentionally empty to emphasize that sampling depth did not extend into 1905.
Figure 4. Contemporary and historical spatial autocorrelation correlograms for Mexican piñon (Pinus cembroides var. cembroides Zucc.) of study site 1. Dotted lines in the middle of each graph indicate the 0.05 significance level, and stars (★) indicate fewer than 20 neighbors for the specified distance class. The 1905 correlogram is left intentionally empty to emphasize that sampling depth did not extend into 1905.
Forests 11 01079 g004
Figure 5. Contemporary and historical spatial autocorrelation correlograms for Mexican piñon (Pinus cembroides var. cembroides Zucc.) of study site 2. Dotted lines in the middle of each graph indicate the 0.05 significance level and stars (★) indicate fewer than 20 neighbors for the specified distance class.
Figure 5. Contemporary and historical spatial autocorrelation correlograms for Mexican piñon (Pinus cembroides var. cembroides Zucc.) of study site 2. Dotted lines in the middle of each graph indicate the 0.05 significance level and stars (★) indicate fewer than 20 neighbors for the specified distance class.
Forests 11 01079 g005
Figure 6. Anselin Local Moran’s I values for Mexican piñon (Pinus cembroides var. cembroides Zucc.) (top left), alligator juniper (Juniperus deppeana var. deppeana Steud.) (top right), and piñon–juniper (bottom middle) age data for 1965, 1985, and 2005 within study site 1. Trees of relatively similar age are presented by black circles, whereas trees of relatively dissimilar age are presented by open circles.
Figure 6. Anselin Local Moran’s I values for Mexican piñon (Pinus cembroides var. cembroides Zucc.) (top left), alligator juniper (Juniperus deppeana var. deppeana Steud.) (top right), and piñon–juniper (bottom middle) age data for 1965, 1985, and 2005 within study site 1. Trees of relatively similar age are presented by black circles, whereas trees of relatively dissimilar age are presented by open circles.
Forests 11 01079 g006
Figure 7. Anselin Local Moran’s I values for Mexican piñon (Pinus cembroides var. cembroides Zucc.) (top left), alligator juniper (Juniperus deppeana var. deppeana Steud.) (top right), and piñon–juniper (bottom middle) age data for 1965, 1985, and 2005 within study site 2. Trees of relatively similar age are presented by black circles, whereas trees of relatively dissimilar age are presented by open circles.
Figure 7. Anselin Local Moran’s I values for Mexican piñon (Pinus cembroides var. cembroides Zucc.) (top left), alligator juniper (Juniperus deppeana var. deppeana Steud.) (top right), and piñon–juniper (bottom middle) age data for 1965, 1985, and 2005 within study site 2. Trees of relatively similar age are presented by black circles, whereas trees of relatively dissimilar age are presented by open circles.
Forests 11 01079 g007
Table 1. The observed, expected, R, and the Z score for contemporary and historical mean nearest neighbor distance of Mexican piñon (Pinus cembroides var. cembroides Zucc.) and alligator juniper (Juniperus deppeana var. deppeana Steud.) jointly and separately.
Table 1. The observed, expected, R, and the Z score for contemporary and historical mean nearest neighbor distance of Mexican piñon (Pinus cembroides var. cembroides Zucc.) and alligator juniper (Juniperus deppeana var. deppeana Steud.) jointly and separately.
SpeciesYearNumber of Trees (n)Observed MeanExpected MeanRZ Score
Piñon–juniper–oak20055121.5120.070.08−40.06 *
Piñon–juniper20053601.8623.930.08−33.47 *
19853361.9224.770.08−32.35 *
19653231.9025.270.08−31.80 *
19452472.3128.830.08−27.65 *
1925924.1047.100.09−16.75 *
1905397.4370.390.11−10.69 *
Mexican piñon20051563.0636.160.08−21.87 *
19851393.2238.300.08−20.66 *
19651373.2238.580.08−20.52 *
19451183.5341.480.09−19.01 *
1925405.4969.370.08−11.14 *
1905117.997.111.120.78
Alligator juniper20052042.4931.650.08−25.18 *
19851972.5532.210.08−24.73 *
19651862.5133.140.08−24.11 *
19451293.2439.790.08−19.96 *
1925526.1562.650.10−12.44 *
1905288.9983.080.11−9.03 *
* Indicates a calculated Z outside the range of critical Z values (−1.96 ≤ ZCritical ≤ 1.96) for an alpha level of 0.05 and a spatial pattern that is significantly different from that expected by an independent random process.
Table 2. Observed, expected, variance, and Z scores of contemporary and historical Moran’s I spatial autocorrelation statistic for Mexican piñon (Pinus cembroides var. cembroides Zucc.) and alligator juniper (Juniperus deppeana var. deppeana Steud.). Inverse distance was used as the weighting function in Moran’s I calculation.
Table 2. Observed, expected, variance, and Z scores of contemporary and historical Moran’s I spatial autocorrelation statistic for Mexican piñon (Pinus cembroides var. cembroides Zucc.) and alligator juniper (Juniperus deppeana var. deppeana Steud.). Inverse distance was used as the weighting function in Moran’s I calculation.
SiteYearNumber of Trees (n)Observed Moran’s I bExpected Moran’s IVarianceZ Score
Mexican piñon
12005570.13−0.020.00372.47 *
1985570.13−0.020.00372.47 *
1965550.16−0.020.00432.8 *
1945400.07−0.030.00591.25
19255−0.13−0.250.50180.17
19050NA aNANANA
22005820.04−0.010.00062.01 *
1985820.04−0.010.00061.99 *
1965820.04−0.010.00061.99 *
1945780.04−0.010.0071.85
1925350.00−0.030.00260.57
190511−0.13−0.100.0074−0.34
Alligator juniper
12005380.15−0.030.00402.80 *
1985380.15−0.030.00402.80 *
1965322.05−0.030.00652.94 *
1945200.00−0.050.01790.40
192515−0.05−0.070.03980.11
1905110.15−0.100.08440.85
220051590.01−0.010.00021.18
19851590.01−0.010.00021.18
19651540.01−0.010.00021.24
19451090.00−0.010.00030.55
192537−0.06−0.030.0051−0.42
190517−0.10−0.060.0030−0.67
* Indicates a calculated Z outside the range of critical Z values (−1.96 ≤ ZCritical ≤ 1.96) for an alpha level of 0.05 and a spatial autocorrelation that is significantly different from that expected by an independent random process. a Not available; b Positive values indicate positive spatial autocorrelation and negative values indicate negative spatial autocorrelation.

Share and Cite

MDPI and ACS Style

Bataineh, M.M.; Oswald, B.P.; Williams, H.M.; Unger, D.R.; Hung, I.-K. Reconstruction of Piñon–Juniper Woodlands in the Sky Islands of the Davis Mountains, Texas, USA. Forests 2020, 11, 1079. https://doi.org/10.3390/f11101079

AMA Style

Bataineh MM, Oswald BP, Williams HM, Unger DR, Hung I-K. Reconstruction of Piñon–Juniper Woodlands in the Sky Islands of the Davis Mountains, Texas, USA. Forests. 2020; 11(10):1079. https://doi.org/10.3390/f11101079

Chicago/Turabian Style

Bataineh, Mohammad M., Brian P. Oswald, Hans M. Williams, Daniel R. Unger, and I-Kuai Hung. 2020. "Reconstruction of Piñon–Juniper Woodlands in the Sky Islands of the Davis Mountains, Texas, USA" Forests 11, no. 10: 1079. https://doi.org/10.3390/f11101079

APA Style

Bataineh, M. M., Oswald, B. P., Williams, H. M., Unger, D. R., & Hung, I. -K. (2020). Reconstruction of Piñon–Juniper Woodlands in the Sky Islands of the Davis Mountains, Texas, USA. Forests, 11(10), 1079. https://doi.org/10.3390/f11101079

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