Next Article in Journal
Effects of CO2 and Soil Moisture Treatments on Morphological and Allometric Trait Variation in Coppiced Seedlings: A Study of Four Early-Successional Deciduous Species
Next Article in Special Issue
Range-Wide Assessment of Recent Longleaf Pine (Pinus palustris Mill.) Area and Regeneration Trends
Previous Article in Journal
Relative Expression of Genes Elicited by Clonostachys rosea in Pinus radiata Induces Systemic Resistance
Previous Article in Special Issue
Differentiating Historical Open Forests and Current Closed Forests of the Coastal Plain, Southeastern USA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Currently Available Site Index Equations That Use On-Site Tree Measurements for Naturally Regenerated Longleaf Pine in Its Historical Range

by
Curtis L. VanderSchaaf
Central Mississippi Research and Extension Center, Mississippi State University, Raymond, MS 39154, USA
Forests 2024, 15(5), 855; https://doi.org/10.3390/f15050855
Submission received: 1 March 2024 / Revised: 2 May 2024 / Accepted: 6 May 2024 / Published: 14 May 2024
(This article belongs to the Special Issue Longleaf Pine Ecology, Restoration, and Management)

Abstract

:
The longleaf pine (Pinus palustris P. Mill.) forest type occupied a much greater area in the United States prior to the arrival of Europeans, estimated to be around 37.2 million hectares. This area has been greatly reduced, and these ecosystems now occupy only about 1.2 to 1.6 million hectares. However, there has been a great interest in the restoration of this forest type mainly due to concerns about the loss of ecosystem services associated with these forests; the improved seedling quality and yield potentials bolster those efforts. Beyond that, existing stands are actively managed through different types of practices, including thinnings, prescribed burns often to manipulate the vegetation of other species, and the various timings of clearcuts. Thus, managers need tools to estimate site quality and ultimately productivity. A commonly used measure of site quality is site index, or the height of some defined dominant portion of the stand at a standardized base or index age. The primary objectives are to summarize the 16 existing equations to estimate site index and dominant height in naturally regenerated longleaf pine stands and to examine and visually compare their predicted behavior across a range of site quality and age conditions. Important considerations when using site index of anamorphism and polymorphism as well as base-age invariance are reviewed. Biologically, polymorphism is often considered advantageous since for many species differences in site quality not only result in different asymptotic dominant heights, but also varying rates in their approach to the asymptote. Of the 16 equations examined, only nine of them were polymorphic in nature, but all equations were base-age invariant. There is not an individual equation that is clearly superior because, for instance, it is either anamorphic in nature, is polymorphic but developed based on anamorphic curves, fit using data obtained from temporary plots, or it is limited geographically. Given these limitations, others can use this publication as a reference to determine which equation they feel is best for their particular situation.

1. Introduction

The longleaf pine (Pinus palustris P. Mill.) forest type was one of the most extensive in North America prior to the arrival of Europeans. It has been estimated that these ecosystems were dominant on nearly 37.2 million hectares [1]. However, it is thought that they now occupy only about 1.3 million hectares [1], but estimates vary to some extent (e.g., Table 1). The species was historically found in a range from southern Virginia to East Texas (Figure 1) and occupied a significant amount of land in the southern half of the Gulf Coastal Plain [2,3]. Fire suppression and exclusion, replanting with faster-growing species such as loblolly pine (Pinus taeda L.) and slash pine (Pinus elliottii Engelm.) following the initial cutover of sites with westward expansion in the late 19th and early 20th centuries [4] and continuing today, and denser stands to increase timber production have transformed the southeastern United States for centuries contributing to the substantial decrease in longleaf pine-dominated forests [4,5].
In its native area, longleaf pine can be a viable alternative to other southern yellow pines, such as loblolly pine and slash pine, because it is resistant to many pathogens and insects and has relatively low site requirements [2]. With proper management, more valuable poles can be produced [6,7,8]. Additionally, longleaf may incur less damage during hurricanes [9,10]. Concerns about the lost ecosystem services associated with a reduction in longleaf pine ecosystems [4,5] have prompted restoration efforts. Over the past 35 years or so, and continuing today, there has been a strong interest in restoring many longleaf pine ecosystems, including the use of natural regeneration approaches such as through shelterwood and seed tree operations. Thus, estimates of site quality and growth and yield are needed to help make optimal management decisions to meet landowner objectives [6]. According to the USDA Forest Service Forest Inventory and Analysis (FIA) program [11], there are around 893,000 naturally regenerated hectares in the southeastern United States (Table 1).
Table 1. Amount of USDA Forest Service Forest Inventory and Analysis (FIA) forestland estimated hectares by state in the southeastern United States of the longleaf pine forest type by stand origin [11]. The amount of hectares established using seeding practices is not reported by the FIA program.
Table 1. Amount of USDA Forest Service Forest Inventory and Analysis (FIA) forestland estimated hectares by state in the southeastern United States of the longleaf pine forest type by stand origin [11]. The amount of hectares established using seeding practices is not reported by the FIA program.
StateNaturalPlantedTotal
Alabama148,034157,861305,895
Florida289,130115,188404,318
Georgia73,392181,887255,279
Louisiana48,95713,39162,348
Mississippi101,55323,998125,551
North Carolina106,77452,210158,985
South Carolina109,278108,641217,919
Texas15,479547620,955
Virginia041744174
Total892,598662,8261,555,424
The primary objective of this publication is to summarize the currently available site index and dominant height curves and equations for naturally regenerated longleaf pine stands that predict site index and/or dominant height as a direct function of on-site tree measurements. In addition, predictions from these equations are visually compared to see their behavior across a range of site conditions.

1.1. Site Quality Estimation

Often site quality estimates are used to directly identify the applicability of management actions and to determine the ability of a site to meet certain management objectives, and they are often included in growth and yield models for the prediction of future yields. A commonly used measure of site quality is site index (e.g., [12]), or the height of some defined dominant portion of the stand at a standardized base or index age ([13,14] and [15] (pp. 304–305)). Site index is used across the world [16,17,18]. Conceptually, for a particular species, site index is a collective influence of soil factors and climatic conditions, and when excluding extremes, site index is thought to be independent of stand density ([13,14] and [15] (p. 303)). Site index is both advantageous and non-advantageous because it is a function of the existing trees—thus, the existing genetics and management practices of the current rotation and previous rotations. However, it is often non-advantageous because it does not provide a direct explanation of site growing conditions. Although site index does not perfectly differentiate productive capabilities among sites ([15] (pp. 311–312) and [19,20]), it is a simple-to-understand tool and it is widely accepted as a means of differentiating the ability of sites to produce a particular product class.

1.2. Site Index Construction and Configuration

Initially, site index curves were developed using interpolative graphical methods (e.g., [21,22]). Some measure of site trees, “dominant height,” was conducted across a range of sites, age was determined, and the paired height measurements were plotted over age. Site index curves were then visually interpolated. However, as statistics became more common and more sophisticated computing ability became available, empirically based site index curves were developed from some definition of height and age. For empirically based site index curves, there are two generally accepted forms, or types [23]. The first is anamorphic curves and the second is polymorphic curves.

1.2.1. Anamorphic or Guide Curve

Anamorphic curves are often referred to as “guide curve” curves. This approach statistically fits a “guide curve” through the data (paired height–age observations) based on some model fitting criteria—often the minimization of squared deviations. A regionally, and silviculturally, applicable species base age is then selected and the “guide curve” site index is determined. Applicable site indexes are then graphically or mathematically obtained by relating the “guide curve” site index to desired site indexes. As noted by many, site index curves using this approach achieve different asymptotic (maximum) heights, but a constant proportional shape of curves is assumed regardless of site quality. This statement is somewhat confusing though. Often, the curves of different site qualities will not have the exact same shape through time visually or graphically. This is an important clarification, since a constant shape across time is often not observed graphically because the proportional rate among site indexes is constant across time, but the rate of change in height for the guide curve is not constant across time. As ([21], p. 319) state in their work:
“To draw the curves in this manner is equivalent to an assumption that differences in site result in similar differences in the rate of height growth at all ages”.
Thus, the constant proportion of sites to the guide curve is multiplied by the rate of change in height, and therefore the shape is not always constant across time among all the site index classes. Thus, anamorphic curves can often have a polymorphic curve appearance, making it difficult sometimes to differentiate between anamorphic and polymorphic curves graphically.
When applying site index curves, it is important to take into account the range of ages used in construction. As noted by many studies, most likely the majority of the predicted dominant heights at older ages are biased low, or biased down (e.g., [22,24,25,26]). Stated differently, height is underpredicted. Due to the long-term nature of forestry, it could be stated with some degree of confidence that all site index equations have been developed from operational landbases, not experimentally designed research studies. Thus, often very few older stands exist on the landscape from which measurements can be made because the majority of the stands are harvested at younger ages. A confounding factor is that usually the older stands that do exist on the landscape are of poorer quality. This usually results from factors such as harvesting costs are not justified on lower-quality sites or silvicultural investment being made in other more inherently productive stands [24,26]. An example can be seen in Table 1 of the publication by [27] where only two sites were used in the equation fitting of ages 71–90 years for the site indexes of 76–85 feet (23.2 m to 25.9 m, base age 50 years), and none were used for the site indexes of 86–95 feet (26.2 m to 29.0 m). There were a total of 38 observations for ages 71–90 years, and hence the “section” of curves for those ages are from lower-quality sites. Similar behavior can be seen in ([28] (Table 96, p. 116)) and [29]. Therefore, since the guide curve is fit through the middle of the data, and older ages are only represented by poorer sites, when site index is desired on better quality sites at older ages the estimated height–age patterns will likely be biased low.

1.2.2. Polymorphic

Polymorphic site index curve approaches were developed based on the findings that for most species the rate of height change through time is not constant across site qualities and that growth rates are not proportional among site qualities, and height curves for a species differ in shape depending upon the site [22,23]. Higher site qualities not only possess a greater height asymptote (maximum attainable height), but they approach and ultimately reach that asymptote faster, and hence sooner, at earlier ages. Lower site qualities have a lower height asymptote and most often a slower approach to that asymptote. Bull (1931, [30]) is often credited as being the first to develop polymorphic site index equations regardless of species. Clutter et al. (1983, [23]) define two types of polymorphic site index curves. The first being polymorphic-disjoint and the second being polymorphic-nondisjoint. For polymorphic-disjoint curves, unlike anamorphic curves, the constant proportional relationship among site qualities does not necessarily hold, but, similar to anamorphic curves, the curves still do not cross. However, not only does the constant proportionality not generally hold for polymorphic-nondisjoint curves, but the curves of site qualities may also actually cross one another at some point in time.

1.3. Base-Age Invariant

An issue with many initially developed polymorphic site index curves was that site index needed to be determined at a preselected base age [31,32]. This was problematic because often the predicted dominant height development would differ, for the same dataset, based on what base age was preselected [31]. Hence, the polymorphic curves were specific to the preselected base age [32]. Bailey and Clutter (1974, [32]) were one of the first to introduce mathematical and statistical approaches to avoid having to use a preselected base age, rather the equation form would allow the user to select any desired base age. For southern yellow pine species, this has been extremely important since in many applications base ages have changed from 50 years to 25 years. Base-age invariant equations are advantageous because any base age can be selected. Because of the constant proportionality assumption, and the manner in which anamorphic curves are created, there are usually no issues associated with base-age invariance for anamorphic curves.

2. Available Site Index Equations for Naturally Regenerated Longleaf Pine Forests

Many site index and dominant height curves and equations requiring on-site tree measurements have been developed for naturally regenerated longleaf pine forests (Table 2, Table 3 and Table 4). Kush et al. (2006, [6]) provided a brief review of longleaf pine site index curves. The first site index curve is from [28]. It is a series of graphical curves only, and no equation was provided, with site qualities ranging from 45 to 125 feet (13.7 m to 38.1 m) for a base age of 50 years. The curves were developed using data from across the southeastern US, and most likely entirely from naturally regenerated stands. The measurement plots were established in fully stocked stands, or in “normal” stands. Farrar (1973, [33]) provided equations depicting the graphically only site index curves found in [28]. Kush et al. (2006, [6]) state that the conditioning of Farrar’s equation resulted in it failing to adequately represent the behavior of the [28] curves for ages 15 years and younger. The equation found in [33] was used as the site index equation in [34], where [34] defined site index to be dominants and co-dominants. Carmean et al. (1989, [22]) took interpreted measurements from the site index curves found in [28] and fit an equation capable of polymorphic behavior, base age 50 years. As stated in [35], the USDA Forest Service Forest Vegetation Simulator (FVS) uses [22]’s Equation (2), while the equation found in [36], Equation (3), is used by the USDA Forest Service Forest Inventory and Analysis (FIA) program [37].
The equations found in [28,36] were presented within overall growth and yield model systems. In some ways, it could be argued that [40] presented site index curves for longleaf stands based on his two-point concept. Only the equations produced by [22,39] exhibit polymorphic behavior. All equations are base-age invariant, and although in many ways anamorphic curves do not adequately represent the biological height patterns of naturally regenerated longleaf stands, anamorphic curves have the advantage of being base-age invariant.
Equations (5)–(13) and (16) were developed using the same naturally regenerated dataset established in the states of Alabama, Florida, Georgia, and Mississippi [27], commonly referred to as the Regional Longleaf Pine Growth Study (RLGS). Farrar (1981, [27]) produced initial curves using this dataset with a 4th-degree polynomial equation form, fit using weighted least squares with the inverse of age2 as the weight. It does not appear that serial correlation was accounted for. When using additional data beyond that used by [27], Rayamajhi et al. (1999) [29] updated his fourth-degree polynomial equation, and fit a second-degree polynomial equation and a Chapman–Richards equation. Ordinary least squares was used, and serial correlation was not accounted for. Rayamajhi et al. (1999, [29]) ultimately suggested their 2nd-degree polynomial equation. Lauer and Kush (2010, [39]) used the most recently updated dataset. Lauer and Kush (2010, [39]), using additional data beyond that used by [29], developed a dynamic site index equation, using as the base the Hossfeld model form, capable of producing polymorphic and even site-specific site index curves. Of all the equations discussed in this paper, only Equations (11) and (13) account for serial correlation when estimating parameters, and Equation (12) does not account for serial correlation. Lauer and Kush (2010, [39]) recommended Equation (11), and it was subsequently used as the site index equation within an even-aged, naturally regenerated growth and yield model system presented in [41]. For the RLGS dataset, different equations not only resulted from the updated datasets, but also the selection of different functional forms. Carmean et al. (1989, [22]) took the equation from [27] and produced a polymorphic site index equation of the Chapman–Richards form. A different weighting approach from [27] was used by [22] when estimating the parameters for the site index equation.
It should be noted that there is one potential equation for naturally regenerated stands not included in Table 2, Table 3 and Table 4 which was presented by [42]. The publication is hard to find, but predictions for a site index of 60 feet (18.3 m, base age 50 years) were presented by [40]. Zeide (1999, [40]) never commented on whether the data was obtained from naturally regenerated sites or plantations. Popham et al. (1979, [43]) presented the equation form used by [42] without the parameter estimates, but they also did not mention the data source.
All site index curves and equations referenced in Table 2, Table 3 and Table 4 were developed using data obtained from even-aged, naturally regenerated stands. The RLGS (Equations (5)–(13) and (16)) measured longleaf pine trees in permanently established 0.1- and 0.2-acre plots (0.0405 hectares to 0.0809 hectares, [27]), and plots were added as the study progressed [38,39]. Plots were mostly measured at 5-year intervals [38,39]. As previously mentioned, only the [39] Equations (11) and (13) were developed taking into account serial correlation resulting from repeated measurements. Schumacher and Coile (1960, [36]) measured trees in 0.2-acre (0.0809 hectares) temporary plots. Specific plot measurement approaches were never clearly specified in [28], as noted also by [22] since the number of dominant trees per plot was not provided. Only Equations (5)–(13) and (16) were fit using repeated measurement data which is considered superior to temporary plots when developing site index and dominant height equations [23] (p. 46).
Due to issues with determining plot age in naturally regenerated even-aged longleaf pine stands, and also varying sample trees per plot, [39] modeled individual trees rather than plot averages [23] (pp. 39–40). They stated that the use of plot averaged age and height when using Equation (11) does not produce substantial errors, particularly for stands over 15 years of age, but technically individual tree height projections should be obtained first, and then these estimates should be averaged to obtain plot-level site index estimates. The users of Equations (11)–(13) should be aware of this issue when applying these equations.

3. Dominant Height Predictions

Figure 2 shows substantial differences in predicted “dominant” height by Equations (1), (3), (5), (7)–(16) across time for the stands with the site indices of 70 and 90 feet (21.3 m and 27.4 m, base age 50 years). The only polymorphic equations were produced by [22] and [39]. For both site indices 70 and 90 feet (21.3 m and 27.4 m, base age 50 years), at age 80 years, Equation (7) (found in [38]) predicted the shortest heights at 77.6 feet and 99.7 feet (23.7 m and 30.4 m), while Equation (1) (found in [33]) had the greatest heights of 84.8 feet and 109.0 feet (25.8 m and 33.2 m), respectively. Carmean et al. (1989, [22]) presented two separate equations when estimating the curves of [27,28,36], one to estimate dominant height from site index (used in Figure 2 and presented in Table 4) and one to estimate site index from dominant height (presented in Table 3). For Figure 2, since dominant height is being predicted, the parameter estimates from the dominant height equations were used (Table 4).
The recommended equation by [39], Equation (11), predicted the values of 84.4 feet and 103.4 feet (25.7 m and 31.5 m) for the site indexes of 70 feet and 90 feet (21.3 m and 27.4 m, base age 50 years) at age 80 years, respectively. Equation (11), found in [39], was fit using an updated RLGS dataset (originally described by [27]) relative to [38], and they used different equation forms (Table 3). For a site index of 70 feet (21.3 m, base age 50), at age 80 years, Equations (5), (12), (13), and (16) all predict similarly to Equation (11), but Equations (8)–(10) generally predict slightly shorter heights. Predictions for Equations (8), (9), and (10) were 80.3 feet (24.5 m), 77.7 feet (23.7 m), and 77.6 feet (23.7 m), respectively.
Similarly, for a site index of 90 feet (27.4 m, base age 50 years), Equations (9) and (10) predict shorter heights; the heights were 99.9 feet (30.4 m) and 99.8 feet (30.4 m), respectively. However, Equation (8) predicted height to be 103.2 feet (31.5 m) which is essentially the same as Equation (11), and height predictions from Equations (5) and (16), which were 107.6 feet (32.8 m) and 106.3 feet (32.4 m), respectively, exceeded that of Equation (11).
The differences in the estimates, and the relative behavior of the equations across site qualities, could be due to the anamorphic nature of the equations presented by [27,29,38], different model forms, as well as the addition of older aged data in RLGS when equations were fit by [39]. Although Equation (16), found in [22], is polymorphic in nature, data were obtained from predictions using an anamorphic equation—Equation (5) from [27].
Since Equation (14), which is presented by [22] and based upon the published curves of [28], is used by the FVS growth and yield model simulator for longleaf pine forests, and Equation (11) is recommended by [39], and since Equation (3) as reported by [36] is used by FIA, Figure 3 presents predictions only from these equations. In addition to that, Equation (5) presented by [27] was included since it is the original equation from the RLGS. Carmean et al. (1989, [22]) presented two separate equations when estimating the curves of [28], one to estimate dominant height from site index (used in Figure 2 and Figure 3 and presented in Table 4) and one to estimate site index from dominant height (presented in Table 3). For Figure 3, since dominant height is being predicted, the parameter estimates from the dominant height equation were used.
For the lower site quality (site index 70 (21.3 m), base age 50 years), at age 80 years, predictions do not differ substantially among the four equations (Figure 3). The equation by [22] (Equation (14)) had the greatest predicted height of 84.6 feet (25.8 m) while the equation by [36] (Equation (3)) had the shortest predicted height of 80.1 feet (24.4 m). Equation (5), as found in [27], exhibits odd behavior for both site indexes for ages younger than 15 years (as noted by [27]), and it is the curve that has the greatest predictions from ages 20 to 40. Most likely this erratic behavior is due to the use of a polynomial equation form (a 4th-degree polynomial). Often polynomial equations provide very good model fits of the data, but do not extrapolate well beyond the range of conditions used to fit the equation.
For the higher site quality (site index 90 (27.4 m), base age 50 years), at age 80 years, the dominant heights do not differ substantially (5.8 feet, 1.8 m). Similar to the lower quality site, the equation by [22] had the greatest height (108.8 feet, 33.2 m), and the equation by [36] (Equation (3)) had the shortest predicted height of 102.9 feet (31.4 m). Remember the equation from [22] (Equation (14)) is a polymorphic equation developed based on anamorphic dominant height–age paired data from [28]. This could lead to some inflexibility in the predictions across age, and perhaps result in greater predictions at older ages. Equations (5) and (11) were produced using the same source of data; however, the equation by [39] (Equation (11)) was developed using a substantially updated dataset. For both site qualities, at younger ages, [36] (Equation (3)) produces shorter trees relative to equations by [22] (Equation (14)) and [39] (Equation (11)). Schumacher and Coile (1960, [36]) developed the equation using stands in the Atlantic Coastal Plain of the southeastern US while those produced by [27] (Equation (5)) and [39] (Equation (11)) are from the East Gulf Region. Schumacher and Coile (1960—Equation (3), [36]) may differ in part from [22] (Equation (14)) and [39] (Equation (11)) because of its anamorphic nature.
It is hard to know specifically the ages used to produce Equation (3) [36]. Although stands as young as 11 years of age were used by [27], the dataset is dominated by stands 25 years of age and older.

4. Predictions of Site Index

Site index predictions (base age 50 years) for a 15-year-old stand with a dominant height of 35 feet (10.7 m, condition 1) differ substantially among equations while condition 2, a 35-year-old stand with a dominant height of 58 feet (17.7 m), has far less variability (Table 5). For condition 1, the site index values average 96.6 feet (29.4 m) and range from 88.7 (27.0 m) to 105.8 feet (32.2 m), a 19.3% difference. The site indices have a standard deviation of 5.49 feet (1.7 m). For condition 2, the site index values average 69.8 feet (21.3 m) and range from 66.9 (20.4 m) to 71.7 feet (21.9 m), a 7.1% difference. The site indices have a standard deviation of 1.67 feet (0.5 m).
Condition 1 is on the fringe of the majority of modeling fitting data ages (Table 2), while condition 2 is well within the range of the model fitting data ages, and thus this may explain in part why predictions for condition 2 have less variability among the equations. In relation to condition 1, it is known that site index estimates obtained when using younger data are often inconsistent and thus there are often significant errors.
Despite their site index curves, Equations (5)–(13), coming from nearly the same model fitting RLGS dataset, for condition 1, clearly model form resulted in different predictions. For the RLGS dataset, there was some variability among the anamorphic equations, but all anamorphic equations overpredicted site index compared to the polymorphic equations of both [22,39]. It could very well be that, when using younger data such as age 15 dominant heights, an anamorphic equation form is not flexible enough to account for reduced height growth near the ages of 50 years and beyond [39].
Farrar (1981, [27]) compared predictions from his equation (Equation (5)) with predictions from [28] and [36] (Equation (3)). He concluded for the RLGS dataset, based on his selected approach, that Equation (5) produced slightly better estimates of site index (base age 50 years). However, he only recommended his equation for the RLGS geographic region—the East Gulf. Rayamajhi et al. (1999, [29]) concluded that their second-degree polynomial Equation (8) produced better estimates relative to the fourth-degree polynomial Equation (5) from [27]. Beyond that, since Equation (8) is a second-degree polynomial equation it is more parsimonious, which is advantageous. When comparing the estimated site indexes (base age 50 years) of RLGS measurement plots, [39] concluded that Equation (11) improved the estimation of site index in young stands relative to Equation (5) from [27]. They noted that the site index estimates from Equation (11) when using younger data were still subject to error though, which is not surprising given it is well known that the use of young ages can produce large errors. For a 15-year-old stand with a dominant height of 35 feet (10.7 m, condition 1 in Table 5), Equation (5) overpredicted by 12.9 feet (3.9 m) in site index (base age 50 years) relative to Equation (11). For other ages, as shown in Figure 2 and Figure 3 as well, [39] found that both equations predict similarly, but that with age Equation (5) tends to overpredict relative to Equation (11).
No previous analysis has compared predictions from Equations (8)–(11). For conditions 1 and 2 in Table 5, Equation (8) overpredicts site index (base age 50 years) relative to Equation (11) for the 15-year-old stand (condition 1) and underpredicts site index (base age 50 years) for the 35-year-old stand (condition 2). Although [39] never compared the two equations, and the results in Table 5 are not a validation analysis, since an updated dataset was used to produce Equation (11) and it is a dynamic equation capable of polymorphic behavior, while Equation (8) is a polynomial anamorphic equation, it is understood that Equation (11) is preferable to Equation (8) for applicable longleaf pine populations.

5. Factors Impacting Site Index and Dominant Height Predictions

Several factors can impact the behavior of dominant tree predictions, and hence greatly impact an estimate of site quality. These factors, then, should be considered when selecting what site index curve is to be used. One factor is the range of ages and likely even the stand densities of the stands used in model fitting. Yet another factor would certainly be whether anamorphic or polymorphic curves are used and yet another would be the functional form, or the mathematical equation itself being used to represent dominant height. For instance, whether a Schumacher model [36], Chapman–Richards equation [22], or Hossfeld [39] model form is used. Yet another important consideration would be whether the equation form allows essentially site-specific predictions (dynamic) such as those of [39]. Lauer and Kush (2010, [39]) obtain site-specific, and hence polymorphic behavior, through the use of a dynamic site index equation approach (GADA—generalized algebraic difference approach). Beyond that, only Equations (5)–(13) and (16) (from the RLGS dataset) were fit using repeated measurement data which often better depicts the long-term trends in height growth than model-fitting data obtained from temporary plots [23] (p. 46). All of the factors above can particularly impact predictions at the extremes of stand age, very young stands and older/mature stands. It should be noted that for the equations from [39], the selection of the years to breast height age (G in Table 3), can have a meaningful impact on dominant height and site index predictions. Beyond that, for [22], the use of either the dominant height (Table 4) or site index (Table 3) parameter estimates can impact predictions.

6. Conclusions

Predictions for conditions that are extrapolations beyond the values used in model fitting may result in odd behavior. For these naturally regenerated longleaf stands, a more older-age-centric dataset may suffice since most managers are aware of erratic height behavior for naturally regenerated longleaf pine stands up to around age 20 years, and sometimes likely even to age 30 years. As found for several other species, and demonstrated here, data across the entire range of ages needs to be combined with a polymorphic equation form which allows for the flexibility needed for faster growth rates at younger ages and slower, more asymptotically related, growth rates at older ages [39]. For naturally regenerated longleaf pine stands, when estimating either dominant height or site index, based on the literature as well as the examinations presented here and others not reported, it likely will be best to use Equation (11) as produced by [39]. However, no one equation examined here is superior across all the properties of concern. For example, although Equation (11) is a dynamic and polymorphic equation fit using repeated measurement data, Equation (11) is geographically limited; technically it should not be used outside of the states of Alabama, Florida, Georgia, and Mississippi. Thus, each user will need to decide for themselves which equation they feel is most appropriate for their particular situation.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflicts of interest.

References

  1. Oswalt, C.M.; Cooper, J.A.; Brockway, D.G.; Brooks, H.W.; Walker, J.L.; Connor, K.F.; Oswalt, S.N.; Conner, R.C. History and Current Condition of Longleaf Pine in the Southern United States; USDA Forest Service, General Technical Report SRS-166; Southern Research Station: Asheville, NC, USA, 2012; 51p.
  2. Boyer, W.D. Longleaf pine: Pinus palustris Mill. In Silvics of North America; Agriculture Handbook 654; U.S. Department of Agriculture Forest Service: Washington, DC, USA, 1990; pp. 405–412. [Google Scholar]
  3. Little, E.L., Jr. Atlas of United States Trees. Volume 1. Conifers and Important Hardwoods; Miscellaneous Publication 1146; Department of Agriculture, Forest Service: Washington, DC, USA, 1971; 9p.
  4. Hanberry, B.B.; Stober, J.M.; Bragg, D.C. Documenting two centuries of change in longleaf pine (Pinus palustris) forests of the Coastal Plain Province, Southeastern USA. Forests 2023, 14, 1938. [Google Scholar] [CrossRef]
  5. Guldin, J.M. Restoration of native fire-adapted southern pine-dominated forest ecosystems: Diversifying the tools in the silvicultural toolbox. For. Sci. 2019, 65, 508–518. [Google Scholar] [CrossRef]
  6. Kush, J.S.; Goelz, J.C.G.; Williams, R.A.; Carter, D.R.; Linehan, P.E. Longleaf Pine Growth and Yield. In The Longleaf Pine Ecosystem: Ecology, Silviculture, and Restoration; Jose, S., Jokela, E., Miller, D., Eds.; Springer: New York, NY, USA, 2006; pp. 9–42. [Google Scholar]
  7. Stottlemyer, A.; Dougal, E.; Zhang, X.; Zehnder, R.; Taylor, E.; Carraway, B. Harvest Trends 2019; Texas A&M Forest Service: College Station, TX, USA, 2021; 24p. [Google Scholar]
  8. Measells, M. Mississippi Timber Price Report—4th Quarter 2022; Starkville, M.S., Ed.; Mississippi State University Extension: Starkville, MS, USA, 2022; 4p. [Google Scholar]
  9. Johnsen, K.H.; Butnor, J.R.; Kush, J.S.; Schmidtling, R.C.; Nelson, C.D. Hurricane Katrina winds damaged longleaf pine less than loblolly pine. South. J. Appl. For. 2009, 33, 178–181. [Google Scholar] [CrossRef]
  10. Haywood, J.D.; Sayer, M.A.S.; Sung, S.-J.S. Comparison of planted loblolly, longleaf, and slash pine development through 10 growing seasons in central Louisiana-an argument for longleaf pine. In Proceedings of the 17th Biennial Southern Silvicultural Research Conference, Shreveport, LA, USA, 5–7 March 2013; Holley, G.A., Connor, K.F., Haywood, J.D., Eds.; 2015; pp. 383–390. [Google Scholar]
  11. USDA Forest Service, Forest Inventory and Analysis Program, Mon, 12 February 2024 21:01:01 GMT; Forest Inventory EVALIDator web-application Version 2.1.0.; U.S. Department of Agriculture, Forest Service, Northern Research Station: St. Paul, MN, USA, 2024. Available online: https://apps.fs.usda.gov/fiadb-api/evalidator (accessed on 12 February 2024).
  12. King, J.E. Site Index Curves for Douglas-fir in the Pacific Northwest; Weyerhaeuser Forest Paper 8; Weyerhaeuser Company, Weyerhaeuser Forestry Research Center: Centralia, WA, USA, 1966; 49p. [Google Scholar]
  13. Burkhart, H.E.; Tennent, R.B. Site index equations for radiata pine in New Zealand. N. Z. J. For. Sci. 1977, 7, 408–416. [Google Scholar]
  14. Ritchie, M.; Zhang, J.; Hamilton, T. Effects of stand density on top height estimation for ponderosa pine. West. J. Appl. For. 2012, 27, 18–24. [Google Scholar] [CrossRef]
  15. Burkhart, H.E.; Avery, T.E.; Bullock, B.P. Forest Measurements, 6th ed.; Waveland Press, Inc.: Long Grove, IL, USA, 2019; 434p. [Google Scholar]
  16. Grey, D.C. Site index—A review. South Afr. For. J. 1989, 148, 28–32. [Google Scholar] [CrossRef]
  17. Kędziora, W.; Tomusiak, R.; Borecki, T. Site index research: A literature review. For. Res. Pap. 2020, 81, 91–98. [Google Scholar] [CrossRef]
  18. Socha, J.; Tymińska-Czabańska, L.; Grabska, E.; Orzeł, S. Site index models for main forest-forming tree species in Poland. Forests 2020, 11, 301. [Google Scholar] [CrossRef]
  19. Amateis, R.L.; Burkhart, H.E. Simulating the effects of site index variation within loblolly pine plantations using an individual tree growth and yield model. In Proceedings of the 18th Biennial Southern Silvicultural Research Conference, Knoxville, TN, USA, 2–5 March 2015; Schweitzer, C.J., Clatterbuck, W.K., Oswalt, C.M., Eds.; 2016; pp. 245–251. [Google Scholar]
  20. VanderSchaaf, C.L. Implications of incorrectly determining site index on stand-level management activities and financial returns in older generation loblolly pine plantations. Reforesta 2023, 16, 1–15. [Google Scholar] [CrossRef]
  21. Bruce, D.; Schumacher, F.X. Forest Mensuration; McGraw-Hill Book Company, Inc.: New York, NY, USA, 1935; 376p. [Google Scholar]
  22. Carmean, W.H.; Hahn, J.T.; Jacobs, R.D. Site Index Curves for Forest Tree Species in the Eastern United States; General Technical Report NC-128; U.S. Department of Agriculture Forest Service, North Central Forest Experiment Station: St. Paul, MN, USA, 1989; 142p.
  23. Clutter, J.L.; Fortson, J.C.; Pienaar, L.V.; Brister, G.H.; Bailey, R.L. Timber Management: A Quantitative Approach; John Wiley and Sons: New York, NY, USA, 1983; 333p. [Google Scholar]
  24. Beck, D.E.; Trousdell, K.B. Site Index: Accuracy of Prediction; Research paper SE-108; Southeastern Forest and Range Experiment Station, U.S. Forest Service: Asheville, NC, USA, 1973; 8p.
  25. Smith, V.G. Asymptotic site-index curves, fact or artifact? For. Chron. 1984, 60, 150–156. [Google Scholar] [CrossRef]
  26. Walters, D.K.; Gregoire, T.G.; Burkhart, H.E. Consistent estimation of site index curves fitted to temporary plot data. Biometrics 1989, 45, 23–33. [Google Scholar] [CrossRef]
  27. Farrar, R.M., Jr. A site index function for naturally regenerated longleaf pine in the East Gulf area. South. J. Appl. For. 1981, 5, 150–153. [Google Scholar] [CrossRef]
  28. USDA Forest Service. Volume, Yield, and Stand Tables for Second-Growth Southern Pines; USDA Misc. Publ. 50; (rev. 1976); USDA Forest Service: Washington, DC, USA, 1929; 202p.
  29. Rayamajhi, J.N.; Kush, J.S.; Meldahl, R.S. An updated site index equation for naturally regenerated longleaf pine stands. In Proceedings of the Tenth Biennial Southern Silvicultural Research Conference, Asheville, NC, USA, 16–18 February 1999; Haywood, J.D., Ed.; 1999; pp. 542–545. [Google Scholar]
  30. Bull, H. The use of polymorphic curves in determining site quality in young red pine plantation. J. Agric. Res. 1931, 43, 1–28. [Google Scholar]
  31. Heger, L. Effect of index age on the precision of site index. Canad. J. For. Res. 1973, 3, 1–6. [Google Scholar] [CrossRef]
  32. Bailey, R.L.; Clutter, J.L. Base-age invariant polymorphic site curves. For. Sci. 1974, 20, 155–159. [Google Scholar]
  33. Farrar, R.M., Jr. Southern pine site index equations. J. For. 1973, 71, 696–697. [Google Scholar]
  34. Lohrey, R.E.; Bailey, R.L. Yield Tables and Stand Structure for Unthinned Longleaf Pine Plantations in Louisiana and Texas; Research Paper SO-133; U.S. Department of Agriculture, Forest Service, Southern Forest Experiment Station: New Orleans, LA, USA, 1977; 52p.
  35. FVS Staff. Southern (SN) Variant Overview—Forest Vegetation Simulator; Internal Report; (revised April 2024); U. S. Department of Agriculture, Forest Service, Forest Management Service Center: Fort Collins, CO, USA, 2008; 89p.
  36. Schumacher, F.X.; Coile, T.S. Growth and Yields of Natural Stands of the Southern Pines; T.S. Coile, Inc.: Durham, NC, USA, 1960; 115p. [Google Scholar]
  37. Vissage, J.S.; Greer, T.R., Jr.; Brandeis, T.J. Site Class and Site Index: Two Estimates of Site Quality for the Southern Research Station Forest Inventory and Analysis Program; e-Res. Note SRS-025; U.S. Department of Agriculture Forest Service, Southern Research Station: Asheville, NC, USA, 2019; 13p.
  38. Farrar, R.M., Jr.; Matney, T.G. A dual growth simulator for natural even-aged stands of longleaf pine in the South’s East Gulf region. South. J. Appl. For. 1994, 18, 147–155. [Google Scholar] [CrossRef]
  39. Lauer, D.K.; Kush, J.S. Dynamic site index equation for thinned stands of even-aged natural longleaf pine. South. J. Appl. For. 2010, 34, 28–37. [Google Scholar] [CrossRef]
  40. Zeide, B. Pattern of height growth for southern pine species. For. Ecol. Manage. 1999, 118, 183–196. [Google Scholar] [CrossRef]
  41. Lauer, D.K.; Kush, J.S. A Variable Density Stand Level Growth and Yield Model for Even-aged Natural Longleaf Pine; Special Report No. 10; Alabama Agricultural Experiment Station: Auburn, AL, USA, 2011; 16p. [Google Scholar]
  42. Larson, E.H.; Moehring, D.M. Site Index Curves for Longleaf Pine in East Texas; Research Notes 1; Department of Forest Science, Texas A&M University: College Station, TX, USA, 1972; 3p. [Google Scholar]
  43. Popham, T.W.; Feduccia, D.P.; Dell, T.R.; Mann, W.F., Jr.; Campbell, T.E. Site Index for Loblolly Plantations on Cutover Sites in the West Gulf Coastal Plain; USDA Forest Service, Research Note SO-250; Southern Forest Experiment Station: New Orleans, LA, USA, 1979; 7p.
Figure 1. Historical range of longleaf pine in the southeastern United States based on Little (1971, [3]). Obtained from https://databasin.org/galleries/5e449f3c91304f498a96299a9d5460a3/#expand=12211, accessed on 30 April 2024.
Figure 1. Historical range of longleaf pine in the southeastern United States based on Little (1971, [3]). Obtained from https://databasin.org/galleries/5e449f3c91304f498a96299a9d5460a3/#expand=12211, accessed on 30 April 2024.
Forests 15 00855 g001
Figure 2. Predicted “dominant height” across time by equation for site indices of 70 feet (21.3 m, left) and 90 feet (27.4 m, right), base age 50 years. Curves from [22] used parameter estimates obtained when dominant height (Table 4), and not site index, was the dependent variable. Thus, Equations (14)–(16) were used.
Figure 2. Predicted “dominant height” across time by equation for site indices of 70 feet (21.3 m, left) and 90 feet (27.4 m, right), base age 50 years. Curves from [22] used parameter estimates obtained when dominant height (Table 4), and not site index, was the dependent variable. Thus, Equations (14)–(16) were used.
Forests 15 00855 g002
Figure 3. Predicted “dominant height” growth for four site index equations. The curves on the left are for a site index of 70 feet (21.3 m) and those on the right are for a site index of 90 feet (27.4 m), both base age 50 years. Equation (14) is from [22] based upon [28], Equation (3) is from [36], Equation (5) is from [27], and Equation (11) is from [39]. The curve from [22] was produced using parameter estimates obtained when dominant height (Table 4), and not site index, was the dependent variable.
Figure 3. Predicted “dominant height” growth for four site index equations. The curves on the left are for a site index of 70 feet (21.3 m) and those on the right are for a site index of 90 feet (27.4 m), both base age 50 years. Equation (14) is from [22] based upon [28], Equation (3) is from [36], Equation (5) is from [27], and Equation (11) is from [39]. The curve from [22] was produced using parameter estimates obtained when dominant height (Table 4), and not site index, was the dependent variable.
Forests 15 00855 g003
Table 2. Equations available for naturally regenerated longleaf pine forests to estimate site index and dominant height. Carmean et al. (1989, [22]) fit polymorphic equations from previously published anamorphic site index curves, where Eqn—equation, SEUS—southeastern US, and ACP—Atlantic Coastal Plain.
Table 2. Equations available for naturally regenerated longleaf pine forests to estimate site index and dominant height. Carmean et al. (1989, [22]) fit polymorphic equations from previously published anamorphic site index curves, where Eqn—equation, SEUS—southeastern US, and ACP—Atlantic Coastal Plain.
EqnAuthorSite StatesAge of Data (Years)Base Age (Years)TypeThinned
[28] *Old-field and cutoverSEUS10–11050AnamorphicUnknown
(1)[33] from [28]Old-field and cutoverSEUS10–11050 (invariant)AnamorphicUnknown
(2), (14)[22] from [28]Old-field and cutoverSEUS10–110InvariantPolymorphicUnknown
(3)[36] **UnknownACPUnknown50 (invariant)AnamorphicUnthinned/Both
(4), (15)[22] from [36] **UnknownACPUnknownInvariantPolymorphicUnthinned/Both
(5)[27]UnknownAL, FL, GA, and MS11–9050 (invariant)AnamorphicBoth
(6), (16)[22] from [27] UnknownAL, FL, GA, and MS11–90InvariantPolymorphicBoth
(7)[38]UnknownAL, FL, GA, and MS15–95InvariantAnamorphicBoth
(8)[29]UnknownAL, FL, GA, and MS20–11050 (invariant)AnamorphicBoth
(9)[29]UnknownAL, FL, GA, and MS20–11050 (invariant)AnamorphicBoth
(10)[29]UnknownAL, FL, GA, and MS20–110InvariantAnamorphicBoth
(11)[39]UnknownAL, FL, GA, and MS~20–120InvariantPolymorphicBoth
(12)[39]UnknownAL, FL, GA, and MS~20–120InvariantPolymorphicBoth
(13)[39]UnknownAL, FL, GA, MS~20–120InvariantPolymorphicBoth
* USDA Forest Service (1929, [28]) only presents published curves and no associated equation. ** Within Schumacher and Coile (1960, [36]) it is difficult to tell definitively if thinned stands were included in the modeling dataset. It is stated, “No plot was taken if past cutting had any apparent effect on stand structure or crown canopy…”.
Table 3. Available site index (SI) equations for naturally regenerated longleaf pine forests where site quality is estimated based on some quantification of on-site tree measurements. Carmean et al. (1989, [22]) fit polymorphic site index curves from previously published anamorphic site index curves, where Eqn—equation and HtD—dominant height.
Table 3. Available site index (SI) equations for naturally regenerated longleaf pine forests where site quality is estimated based on some quantification of on-site tree measurements. Carmean et al. (1989, [22]) fit polymorphic site index curves from previously published anamorphic site index curves, where Eqn—equation and HtD—dominant height.
Eqn from Table 2AuthorEquation
(1)[33] from [28] SI = Ht D 10 11.8701 * 1 50 1 Age + 1263.79 1 50 3 1 Age 3 12 , 409.5 * 1 50 4 1 Age 4
(2)[22] from [28] SI = 0.6333 Ht D 1.0310 1 e 0.0278 Age 1.6207 Ht D 0.0878
(3)[36] SI = 10 L o g 10 Ht D + 6 . 33645 1 Age 1 50 + 44.550 1 Age 2 1 50 2
(4)[22] from [36] SI = 0.7948 H t D 1.0171 1 e 0.0470 A g e 2.1964 H t D 0.0864
(5)[27] SI = Ht D 10 29.468 1 50 1 Age + 938.97 1 50 2 1 Age 2 16 , 102 1 50 3 1 Age 3 + 88 , 775 1 50 4 1 Age 4
(6)[22] from [27] SI = 0.6398 Ht D 1.0401 1 e 0.0345 Age 2.6970 Ht D 0.1749
(7)[38] SI = Ht D e 71.60 Age 3 2 50 3 2
(8)[29] SI = Ht D 10 6.74634 1 50 1 Age 35.95195 1 50 2 1 Age 2
(9)[29] SI = Ht D 10 14 . 67550 1 50 1 Age 1005.43358 1 50 2 1 Age 2 + 17,194 1 50 3 1 Age 3 103 , 876 1 50 4 1 Age 4
(10)[29] SI = Ht D 1 e 0.0568896 Age 2 1 e 0.0568896 Age 1 2.095444
(11)[39] SI n + G =   4.5 + 77.080 + X 0 1 + 1723.39 X 0 * ( B n + G G ) 1.235
where:
X0 = 0.5( Ht D 4.5 − 77.080
+ ( Ht D     4.5 77.080 ) 2 + 4 1723.39 Ht D 4.5 Age G 1.235 ), and
G = age at which trees reach 4.5 feet (breast height).
(12)[39] SI n + G = 4.5 + 73.327 + X 0 1 + 1877.54 X 0 ( B n + G G ) 1.226
where:
X0 = 0.5( Ht D 4.5 − 73.327
+ ( Ht D 4.5 73.327 ) 2 + 4 1877.54 Ht D 4.5 Age G 1.226 ), and
G = age at which trees reach 4.5 feet (breast height).
(13)[39] SI n + G = 4.5 + 77.385 + X 0 1 + 1708.07 X 0 ( B n + G G ) 1.236
where:
X0 = 0.5( Ht D − 4.5 − 77.385
+ ( Ht D 4.5 77.385 ) 2 + 4 1708.07 Ht D 4.5 Age G 1.236 ), and
G = age at which trees reach 4.5 feet (breast height).
Table 4. Available dominant height (HtD) equations produced by [22] for naturally regenerated longleaf pine forests where site quality is estimated based on some quantification of on-site tree measurements. Carmean et al. (1989, [22]) fit polymorphic dominant height curves from previously published anamorphic site index curves, where Eqn—equation and SI—site index.
Table 4. Available dominant height (HtD) equations produced by [22] for naturally regenerated longleaf pine forests where site quality is estimated based on some quantification of on-site tree measurements. Carmean et al. (1989, [22]) fit polymorphic dominant height curves from previously published anamorphic site index curves, where Eqn—equation and SI—site index.
Eqn from Table 2AuthorEquation
(14)[22] from [28] Ht D = 1.4210 SI 0.9947 1 e 0.0269 Age 1.1344 SI 0.0109
(15)[22] from [36] Ht D = 1.1672 SI 1.0010 1 e 0.0459 Age 1.3460 SI 0.0204
(16)[22] from[27] Ht D = 1.3196 SI 1.0000 1 e 0.0356 Age 1.4271 SI 0.000017
Table 5. Site index predictions (base age 50 years) in meters (m) and feet (ft) for two naturally regenerated stand conditions. Condition 1 is for a 15-year-old naturally regenerated stand with a dominant height of 35 feet (10.7 m). Condition 2 is for a 35-year-old naturally regenerated stand with a dominant height of 58 feet (17.7 m). Estimates from [22] were produced using parameter estimates obtained when site index is the dependent variable (Table 3), where Eqn—equation, SEUS—southeastern US, and ACP—Atlantic Coastal Plain.
Table 5. Site index predictions (base age 50 years) in meters (m) and feet (ft) for two naturally regenerated stand conditions. Condition 1 is for a 15-year-old naturally regenerated stand with a dominant height of 35 feet (10.7 m). Condition 2 is for a 35-year-old naturally regenerated stand with a dominant height of 58 feet (17.7 m). Estimates from [22] were produced using parameter estimates obtained when site index is the dependent variable (Table 3), where Eqn—equation, SEUS—southeastern US, and ACP—Atlantic Coastal Plain.
EqnAuthorStatesAge of Data
(Years)
TypeCondition 1Condition 2
ftmftm
(1)[33] from [28]SEUS10–110Anamorphic94.828.971.121.7
(2)[22] from [28]SEUS10–110Polymorphic88.727.071.421.8
(3)[36]ACPUnknownAnamorphic104.731.968.620.9
(4)[22] from [36]ACPUnknownPolymorphic88.927.168.821.0
(5)[27]AL, FL, GA, and MS11–90Anamorphic105.832.267.120.5
(6)[22] from [27]AL, FL, GA, and MS11–90Polymorphic96.029.269.921.3
(7)[38]AL, FL, GA, and MS15–95Anamorphic98.029.966.920.4
(8)[29]AL, FL, GA, and MS20–110 yearsAnamorphic101.030.868.620.9
(9)[29]AL, FL, GA, and MS20–110 yearsAnamorphic100.930.870.021.3
(10)[29]AL, FL, GA, and MS20–110 yearsAnamorphic98.830.169.621.2
(11)[39]AL, FL, GA, and MS~20–120 yearsPolymorphic92.928.371.721.8
(12)[39]AL, FL, GA, and MS~20–120 yearsPolymorphic92.928.371.521.8
(13)[39]AL, FL, GA, and MS~20–120 yearsPolymorphic92.928.371.721.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

VanderSchaaf, C.L. Currently Available Site Index Equations That Use On-Site Tree Measurements for Naturally Regenerated Longleaf Pine in Its Historical Range. Forests 2024, 15, 855. https://doi.org/10.3390/f15050855

AMA Style

VanderSchaaf CL. Currently Available Site Index Equations That Use On-Site Tree Measurements for Naturally Regenerated Longleaf Pine in Its Historical Range. Forests. 2024; 15(5):855. https://doi.org/10.3390/f15050855

Chicago/Turabian Style

VanderSchaaf, Curtis L. 2024. "Currently Available Site Index Equations That Use On-Site Tree Measurements for Naturally Regenerated Longleaf Pine in Its Historical Range" Forests 15, no. 5: 855. https://doi.org/10.3390/f15050855

APA Style

VanderSchaaf, C. L. (2024). Currently Available Site Index Equations That Use On-Site Tree Measurements for Naturally Regenerated Longleaf Pine in Its Historical Range. Forests, 15(5), 855. https://doi.org/10.3390/f15050855

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