- Research Paper
- Published:

# Structure area curves in Eastern Hardwoods: implications for minimum plot sizes to capture spatially explicit structure indices

*Annals of Forest Science*
**volume 78**, Article number: 16 (2021)

## Abstract

### • *Key message*

*Key message*

**Sampling needs differ by forest type for timber inventory and structural complexity metrics. We demonstrate in a typical mixed Eastern Hardwoods forest that optimal sampling of timber inventory metrics and spatially explicit structure indices may be achieved in one large plot plus a cruise for large diameter trees, but accurately capturing inventory metrics may not be possible with sparse large-scale sampling.**

### • *Context*

*Context*

Managing forest stand structures for multiple objectives require accurate and precise estimates of structural features that may be best estimated at different scales.

### • *Aims*

*Aims*

We document minimum necessary plot sizes for structural metrics and spatially explicit indices to characterize structure in a mature North American Eastern hardwoods forest.

### • *Methods*

*Methods*

Metrics and indices (Index of Aggregation, Diameter Differentiation Index, Dissimilarity Coefficient, Structural Complexity Index) were calculated within 0.05–1.75-ha plots for 1000 iterations of random placement in two 2.0-ha macroplots. Estimation adequacy required (1) precision (varied < 10% among plots) and (2) accuracy (within 10% of the 2.0-ha value at 5th and 95th percentiles).

### • *Results*

*Results*

Minimum single plot sizes to achieve estimation adequacy were 0.25–0.75 ha for spatially explicit indices and 0.5–2 ha for stand metrics. A minimum of five 0.10-ha subplots would be needed for most indices and 6–25 for most metrics, but an untenable 375+ for the density of large diameter trees.

### • *Conclusion*

*Conclusion*

Estimation adequacy for structural complexity requires no greater sampling intensity than for timber metrics, except for density of large trees. A single large plot may be most cost-effective. National inventories in Eastern hardwoods may not estimate structural complexity well due to inadequate sampling intensity.

## 1 Introduction

The objective of forest inventory is optimized parameter estimation: accurately characterizing the population of interest while minimizing the resources required to do so. Consequently, there is a long history of optimizing plot sizes, shapes, and layouts to ensure a swift and adequately realistic assessment of the timber parameters (e.g., density, basal area) that inform stocking charts (Bormann 1953; Freese 1967; Zeide 1980; Kenkel and Podani 1991; Avery and Burkhart 2001). Recognizing the importance of incorporating local variability into the optimal design of a sampling scheme (Bormann 1953; Reich and Arvanitis 1992; Avery and Burkhart 2001), most inventory protocols for timber resources employ relatively large numbers of small, widely distributed sampling plots (e.g., Forest Inventory and Analysis (FIA), Gormanson et al. 2018). The efficiency of small, often variable radius, plots for characterizing traditional timber metrics persists across forest types and management histories (Berrill and O’Hara 2012; Du et al. 2018).

Management objectives, however, have evolved, and with them the parameters of interest to inventory on both managed lands and strict reserves. As interest in the complex stand structures integral to biodiversity has grown (Crow et al. 2002; Larson 2007), existing inventory protocols have struggled to keep pace. Although the data from inventories based on many small plots can be used to plot diameter distributions and summarize tree species richness and the density of large diameter trees (e.g., Brown et al. 1997; Crow et al. 2002), both theory and practice indicate that special features often have population characteristics that differ from those of density and tree diameter (Kenkel and Podini 1991). The inherent variability of some features (Franklin and Van Pelt 2004; Spies 2004) often dictates more intensive sampling effort to achieve robust estimates (Nagel et al. 2007; Zenner and Peck 2009; Král et al. 2010). Further, while many traditional metrics correlate well between field measurements and remotely sensed (e.g., LiDAR-derived) estimates, predictive models based on remotely obtained textural data are still often inadequate for structural features such as the density of large diameter trees (Kane et al. 2010) or variation in tree size (Mura et al. 2015; Meng et al. 2016), and poorer than expected (Kekunda et al. 2019) or even demonstrably poor for metrics incorporating spatial arrangement (unless additionally drawing on more costly spectral data; Meng et al. 2016; Kandare 2017).

Because the spatial distribution of features within a stand determines their probability of inclusion in a sampling frame, influencing statistic power, field inventory protocols using many small plots are challenged by the incorporation of structural features typical of older, unmanaged forests, which are highly variable in frequency, abundance, and spatial arrangement (Spies and Franklin 1991; Reich and Aravanitis 1992; Gray 2003). Features that are rare and/or unevenly distributed on the landscape (e.g., large trees) are analogous to rare species, which are better captured in fewer large plots than in more small plots (McCune and Lesica 1992). Traditional protocols can be expanded by tacking on larger supplementary plots (e.g., lichen survey plots, FIA, Gormanson et al. 2018), but the root causes of this bias—the influence of spatial pattern on variability—is largely unaddressed when field protocols prioritize the number of parameters over spatial extent and resolution (Proulx 2007).

Choosing an optimal spatial extent (i.e., plot size), however, is itself challenged by the continuous—and nearly functional (cf. Král et al. 2010)—decrease in variation of most metrics with increasing plot size (Busing and White 1993; Zenner 2005; Zenner and Peck 2009; Berrill and O’Hara 2012; Guillemette et al. 2012; Zenner et al. 2015, 2019; Lombardi et al. 2015; Du et al. 2018; Kekunda et al. 2019). In an effect very like the flattening of the species area curve with sampling effort, which can be captured by plotting “structure area curves” of variation against spatial scale (Zenner 2005), parameter estimation for spatially non-random features often requires larger plot sizes (Kenkel and Podini 1991), especially as forest heterogeneity increases (Kekunda et al. 2019). As spatial extent increases and plots increasingly incorporate diverse features by absorbing multiple patches (Zenner et al. 2019), within-plot variance increases at the expense of between-plot variance (Scott 1998). The heterogeneity evident at small scales is averaged across at larger scales; the resulting homogenization at large scales, known as spatial smoothing (e.g., Zenner and Peck 2018), renders only small gains in estimation efficiency with the addition of more large plots (Kenkel and Podani 1991).

Thus, the need to efficiently incorporate structural features with variable spatial patterns apparently dictates sampling protocols using fewer, larger plot sizes. Further, because plot size influences the assessment of spatial pattern (Zenner and Peck 2009; Fonton et al. 2011; Carrer et al. 2018), the recognition of spatial pattern as a parameter in and of itself (e.g., Aldrich et al. 2003) has different sampling requirements from traditional inventory parameters (Kenkel and Podani 1991). Yet there remains no consensus on a minimum standard plot size for inventorying spatially dependent structural features, even within forest type, because virtually no information is available on the structure area curves of most spatially explicit indices of structural complexity (but see Maleki and Kiviste 2015; Kekunda et al. 2019). Although 1.0-ha sample plots are no longer uncommon when structural features are of interest (e.g., Guillemette et al. 2012; Lombardi et al. 2015; Grotti et al. 2019; Kekunda et al. 2019; Zenner et al. 2019), most inventory protocols continue to use relatively small fixed or variable radius plots (e.g., McGee et al. 1999; Crow et al. 2002; Gormanson et al. 2018), and only a minority of national inventory programs quantify spatially explicit structural features (Winter et al. 2008). As land managers look toward revising these protocols to incorporate structural features, the question remains as to how large a plot is necessary.

The objectives of the current study, therefore, were to (i) derive a structure area curve for a mature Eastern hardwoods forest and (ii) determine the minimum acceptable single large plot size and/or small (0.1 ha) plot sample sizes necessary to estimate structural parameters, including spatially explicit indices, with adequate precision and accuracy.

## 2 Material and methods

### 2.1 Sampling

Two adjacent stands (RS1, SV3) were sampled within a typical Eastern hardwoods forest (Fig. 1) on Penn State Stone Valley Forest in central Pennsylvania, USA (40° 37′ 37–54″ N, 77° 54′ 5–9″ W). The stands are on well-drained channery loam soils, classified as medium sites with even-aged stratified mixed hardwoods dominated by oak (*Quercus*) and maple (*Acer*) with a component of white pine (*Pinus strobus* L.). Both stands originated naturally following clearcutting ca. 1920, received minor improvement cuts mid-century, and had been unmanaged since at least 1985.

Each stand was sampled using a ca. 2.25-ha macroplot (82 × 276 m and 108 × 208 m, respectively) in the summer of 2007, each shaped to best capture the individual stand. In each macroplot, the position of all live trees ≥ 5 cm in diameter at breast height (DBH) was stem-mapped and slope-corrected distances and azimuths were converted to Cartesian coordinates. For each tree, species and DBH were recorded (Peck and Zenner 2021). Both macroplots exhibited a reverse-J diameter distribution.

### 2.2 Parameters

Structure area curves for each macroplot were derived by calculating each parameter at several different scales within each macroplot by overlaying plots of 500, 1000, 2500, 5000, 7500, 10,000, 12,500, 15,000, and 17,500 m^{2} onto each stem map. Plot shapes were kept identical to the 2-ha macroplot core (i.e., rectangular), but rescaled to 0.05 to 1.75 ha in size. We simulated “sampling” by randomly placing 1000 of these variously sized plots within each macroplot. To ensure consistency across plot size/scale, 1000 [*X*,*Y*] coordinate positions were randomly selected in each 2-ha macroplot core and then used as the plot centre for sampling at each spatial scale. Torroidal edge correction was used to ensure equal sample probability of all trees while correcting for edge effects (Griffith 1983; Boots and Getis 1988; Gray 2003). This adjustment joins the opposite ends of a mapped area, creating a continuous surface for random plot placement throughout the mapped area (Gray 2003).

The structure parameters assessed in each plot at each scale included five simple metrics (the mean and standard deviation (STD) of DBH (cm), basal area (m^{2}/ha), tree density (trees/ha), and density of large trees (DBH ≥ 50 cm) (large trees/ha)) and six spatially explicit structure indices (see Zenner et al. 2015 for formulas). The index of aggregation (*R*, Clark and Evans 1954) was used to capture horizontal structure, ranging from 0 (maximum aggregation/clustering of trees) to 2.1491 (a regular hexagonal arrangement) with a value of 1 when the spatial distribution is random. The dissimilarity coefficient (*DC*) of Hagner and Nyquist (1998) was used to quantify vertical structure as size differences among pairs of neighboring trees, ranging from 0 to 1 with a value of 0.5 when tree sizes are drawn independently from an exponential distribution. Vertical structure was also measured as the difference in tree size among four nearest neighbors (*T*, the diameter differentiation index of Füldner 1995, ranging from 0 to 1). Indices were also calculated based on neighborhoods identified after connecting trees to form a triangular surface that, when extended across the entire sampling area, forms a triangulated irregular network (TIN) of non-overlapping triangles (Fraser and van den Driessche 1972). The average difference in tree sizes within these three-tree triangles (*DBHdiff3*), which was also reported scaled from 0 to 1 (*Dd3*), was also used to quantify vertical structure. Finally, both vertical and horizontal structures were assessed using the structural complexity index (*SCI*; Zenner and Hibbs 2000), in which trees are represented as irregularly spaced three-dimensional data points (*x*, *y* = spatial coordinates, *z* = tree DBH). The *SCI* was calculated as the sum of the surface areas of the TINs for a plot divided by the projected ground area of all triangles in the plot.

### 2.3 Estimation

The mean across the 1000 simulated samples was used as an unbiased measure of central tendency for each parameter (except *R*; see below) at each scale. The estimation error (i.e., bias) for each parameter at each scale was presented as the relative deviation of the simulated sample estimate (i.e., the mean across the 1000 iterations at each scale) from the “true” value, defined as the known value for the 2-ha reference macroplot core (Table 1):

where *xij* is the value of the parameter from a plot of size *i* in macroplot *j* and *μj* is the 2-ha value in macroplot *j* (after Gray 2003).

The mean and the 5th and 95th percentiles of the distribution of the estimation error were calculated for the 1000 simulated plots of each parameter in each macroplot. Variation among scales, and therefore adequacy of sampling, was evaluated using two standards. First, precision was determined by identifying the minimum scale at which the parameter estimates varied 10% or less among the iterations for a given scale (i.e., coefficient of variation ≤ 10%). Second, accuracy was determined by identifying the minimum scale at which the parameter estimates at both the 5th and 95th percentiles of the iteration distribution for a given scale were within 10% of the 2.0-ha value (i.e., comparable to an effect size of 10% at a two-sided alpha of 0.05).

To evaluate the trade-off between using a single large plot and a larger number of smaller plots, we used estimates of variance from the simulated plots to calculate the necessary sample size if multiple “subplots” were sampled at different spatial scales to obtain satisfactory estimates for each parameter. Using the average variance across both macroplots (*s*
^{2}), an effect size (*E*) of 10% of the average 2-ha value, and a desired alpha of 0.05, the sample size needed to estimate the mean was calculated as,

(Freese 1967).

Although frequently reported as a measure of horizontal structural complexity, the index of aggregation (*R*) presents a unique statistical and interpretive challenge. While the index provides output in the form of a continuous dataset, the interpretation is nearer to that of a class variable: values near 0 indicate aggregation/clustering, near 1 a random spatial distribution, and those significantly > 1 a regular/dispersed arrangement. Consequently, averaging across these values can smooth a combination of aggregated and dispersed plots to give the impression of spatial randomness. Thus, rather than using the mean across the 1000 iterations to calculate sample size, we report instead the proportion of plots that would be considered aggregated/clustered or dispersed as opposed to random.

All calculations and simulations were performed in Matlab V. 8.2.0 (Mathworks Inc.).

## 3 Results

Due to relatively low variation within simulated plots of a given scale, precision was obtained at smaller single plot extents than accuracy (Table 2). The minimum scale at which estimates were both precise and accurate varied among parameters from 0.25 ha for some structure indices to the full 2.0 ha for the density of large trees. While a single large plot of 0.5 ha would be adequate to sample most spatially explicit structure indices, most simple metrics required at least a 1-ha plot. Despite the proximity and similarity of history for the two stands, larger plot sizes were necessary for RS1 than for SV3 to achieve precision and accuracy for several parameters.

Deviation decreased with increasing scale for all parameters (Fig. 2). Although this overall pattern of typical spatial smoothing was observed in both stands, RS1 had greater overall tree density and mean tree diameters while SV3 had notably more large trees (and somewhat higher basal area). While all metrics and indices, except the density of large trees, deviated on average ≤ 10% from true by 0.75 ha, the rate of spatial smoothing varied among metrics and indices. Very few iterations of most spatially explicit structure indices deviated ≥ 10% by 0.5 ha (Appendix Table 5), such that their means reached ≤ 5% by that same spatial scale (Fig. 2). In contrast, a comparably low number of iterations deviating ≤ 10% was not observed for most simple metrics until 1 ha in size (Appendix Table 5), with means stabilizing at ≤ 5% deviation at 0.75 ha for basal area and the standard deviation of tree diameters, 1.25 ha for mean tree diameter, and 1.5 ha for tree density (Fig. 2).

Consequently, the minimum sampling extent required for adequate precision and accuracy also varied by metric and index (Table 3). Although a single 0.25-ha plot would be sufficient to estimate *Dd3* and *T*, the remaining spatially explicit structure indices and simple metrics required larger extents to achieved desired levels of precision and accuracy: one 0.50 ha plot for *DC*; one 0.75 ha plot for the standard deviation of tree diameters, basal area, and the *SCI*; 1.0 ha for *DBHdiff3*; 1.25 ha for mean tree diameter; and 1.5 ha for the density of all trees. The density of large trees could not be robustly estimated in less than the full 2.0 ha plot. Conversely, if multiple smaller sampling extents were an option, then fifteen 0.25-ha subplots would suffice for all metrics and indices except the density of large trees.

Spatial arrangement varied slightly with scale (Table 4). The proportion of plot iterations with a non-random spatial arrangement declined from 15% at 0.05 ha to 0.05% by 1 ha, with both stands exhibiting a random tree arrangement at 2 ha. However, these non-random arrangements were predominantly aggregated/clustered at the smaller scales in both stands, although by 0.75 ha they were predominantly dispersed in the SV3 stand (at which scale all iterations were already random in the RS1 stand).

## 4 Discussion

Sampling design is invariably faced with a dilemma: the required sampling effort depends upon the variability and spatial distribution of structural attributes within a stand, which are rarely known before an inventory is conducted. Previous research can provide some guidance, given observations that simple metrics such as mean tree diameter and density can require single plot extents on the order of 0.35 ha in beech (Král et al. 2010) and 0.50 in red pine (Zenner and Peck 2009) and old-growth *Tsuga*-mixed hardwoods (Busing and White 1993). Stands of similar management history and forest type provide the best analogues, hinting that, as in Northern hardwoods (Guillemette et al. 2012), a single plot of 0.5 ha in extent may be required to estimate basal area or of 1.0 ha to assess growth parameters. Likewise, metrics such as basal area may be most efficiently captured in mixed hardwoods using numerous small (often < 0.1 ha) dispersed subplots (e.g., Becker and Nichols 2011).

In fact, adequate estimation of simple timber metrics in the mature Eastern hardwoods forest considered here would indeed require a large single sample plot (upwards of 1.0 ha) or numerous smaller fixed area subplots (> 10 0.1 ha). Deviation of the traditional inventory metrics was still high at even large plot sizes, stabilized only at the largest extents, and most metrics could not be adequately estimated from a single plot covering less than at least 50% of the macroplot area. These results confirm the need for a high sampling intensity in even-aged stratified mixed stands, which are typically characterized by a reverse-J diameter distribution (Ashton and Kelty 2018) due to the high density of small-diameter trees, and are consistent with the large sampling extents needed to reliably capture diameter distributions (Rubin et al. 2006) and improve the precision of tree density estimates in similarly structured selection forests (Jazbec et al. 2011).

However, structural features subject to strong species area curve-like trends, such as tree species composition (Busing and White 1993) or the density of large trees (Lombardi et al. 2015), often require even greater sampling effort. As is often the case (Gray 2003; Zenner and Peck 2009), large diameter trees were relatively scarce in these stands, yet are an important parameter for the restoration of Eastern forests (McGee et al. 1999; Crow et al. 2002) because they are thought to have once contributed considerably more biomass in old-growth forests than currently in young managed forests (Brown et al. 1997). Due to the high variability in capturing large trees in individual sample plots, more than ten times as many 0.1-ha subplots would be required to adequately estimate the density of large diameter trees than would have even fit within the macroplots being assessed. Alternatively, adequately estimating this parameter with a single large plot would have required the full extent of the area to be sampled (as was also seen in red pine, Zenner and Peck 2009). In the current study, even supplementary sampling through a nested plot design (e.g., a larger plot around each subplot just for large diameter trees, such as the FIA macroplot; Gormanson et al. 2018) would not have achieved sufficient accuracy and precision if it were less than the full macroplot area. Rather than untenably increasing the number of subplots or plot extent, however, such rare large trees may instead require a separate/additional sampling protocol (Thompson and Burnham 2004). Features that are known to be relatively rare within a stand may be best estimated using an entirely different sampling frame, such as the combination of fixed and variable radius plots that improves sampling precision for tree density (Packard and Radtke 2007), or point (Ritter and Saborowski 2014) or line transects (Bate et al. 1999) that can be most efficiency sampled while traveling between plots (Johnson et al. 2008; Bäuerle et al. 2009).

In contrast to the simple metrics, most of the spatially explicit structural indices could be adequately estimated in a single plot of less than one hectare in extent or in fewer than ten 0.1-ha subplots. Spatial arrangement (*R*) stabilized by 0.5 ha and spatially explicit vertical structure was adequately captured in single plots of 0.25 ha (*T*) or 0.5 ha (*DC*) in size, in keeping with previous observations in old silver birch (Maleki and Kiviste 2015) and mature red pine (Zenner and Peck 2009). However, although not as extensive as that required in savanna woodland, montane conifers or northern hardwoods (up to 1 ha, Fonton et al. 2011; Guillemette et al. 2012; Carrer et al. 2018), a single plot of 0.75 ha would be required to capture both horizontal and vertical structure (*SCI*) in these mature Eastern hardwoods stands. The predominance of a random spatial distribution of trees (almost 90% at 0.1 ha) in these stands contributed to the more rapid stabilization of the structure area curves for most spatially explicit measures of structural complexity than for the simple metrics. Rapid spatial smoothing has also been observed for size class abundances and thus diameter distribution forms, which stabilized by ~ 0.15 ha in old-growth Douglas-fir (Zenner and Peck 2018) and beech (Zenner et al. 2018). Likewise, the pattern of dominance by all-sized tree neighborhoods became clear by 0.1 ha in old-growth beech (Zenner et al. 2019; 2020) and that of subsequently assigned development phases by 0.125 ha (Zenner et al. 2020). In fact, neighborhood-level tree size differences may be most clearly expressed at the fine scales (e.g., < 0.1 ha; Zenner et al. 2019) capturing individual tree processes, which coalesce at larger scales into patterns of tree size distribution (Zenner et al. 2015) (a transition in perspective across scale that is only observable using structure area curves, Zenner 2005).

As a consequence of the fine-scale structural complexity in these mature Eastern hardwoods stands, spatially explicit structural complexity indices did not necessarily require an increase in either sampling extent (for single plots) or intensity (for subplots) over what would already be required to achieve adequate estimates of the simple metrics: i.e., by the time a sufficient number of subplots was sampled for traditional timber metrics, the minimum number required for spatially explicit indices would already be met. On the one hand, this indicates that structural complexity could be tacked on to sampling protocols intended for estimating simple metrics, such as basal area, without requiring additional subplots or even a change in sampling frame from many small subplots to one large plot. On the other hand, however, greater sampling effort is nonetheless required, due to both the necessity of stem-mapping and of sampling smaller trees than is often typical (e.g., smaller than the 12.7 cm cutoff for FIA subplots, Gormanson et al. 2018). Given that the intensive sampling of neighboring trees is more suited to fixed area plots (Berrill and O’Hara 2012) and the greater efficiency of sampling one large plot over many smaller subplots (Jazbec et al. 2011), inventories assessing structural complexity in Eastern hardwoods may be optimally achieved using a single large plot (e.g., 1 ha; Grotti et al. 2019)—although an additional transect or other supplemental approach may be needed for the density of large trees when present.

The implication of these results for large national inventories relying on a small number of sparsely distributed, relatively small subplots (e.g., FIA with four spatially linked 0.07 ha plots) is that they are best suited to what they were designed for: broad trends in the “extent, condition, volume, growth, and use of trees” (Gormanson et al. 2018). Regardless of whether simple metrics or spatially explicit structural indices are involved, the high size variability in the Eastern hardwoods forest type necessitates a high degree of sampling effort when precision and accuracy are desired. Further, meaningful efforts to estimate structural complexity from data collected using inventory protocols like FIA would be futile in Eastern hardwoods, not only because complete stem-mapped data are generally lacking, but because sampling intensities are simply inadequate. This may explain why the recommended applications of FIA data only include assessments of forest structure when additional (e.g., remote sensing) data are available (Tinkham et al. 2018). Thus, the cost-effective estimation of spatially dependent forest structures in Eastern hardwoods is likely still some years away, as currently inadequate remote sensing technologies (Kekunda et al. 2019) continue to evolve (e.g., Meng et al. 2016; Kandare 2017).

## 5 Conclusion

It is often assumed that assessment of spatially explicit measures of structural complexity requires larger plot extents than traditional timber metrics. The results of the current study indicate that adequate estimation of simple timber metrics in spatially random mature Eastern hardwoods would actually require an even larger single plot extent than the spatially explicit indices—particularly to estimate the density of large diameter trees. Although it could be concluded that a trade-off is inevitable and sampling designs must focus on optimizing some metrics over others, a more flexible approach may allow efficient characterization of all parameters of interest through the employment of a mixed inventory sampling scheme. Some combination approaches are already employed, such as a mixture of fixed and variable-radius plots (Packard and Radtke 2007; Gormanson et al. 2018). The current study demonstrates that basal area could be adequately estimated from ~ 10 0.1-ha subplots and that the spatially explicit indices could also be estimated if half of them were also stem-mapped. If the same level of accuracy and precision was desired for the estimation of the density of large trees, they might be most efficiently estimated through line transects connecting these subplots. Alternatively, capturing an array of forest metrics may be most cost-effectively achieved through a combination of sampling a single plot spanning 50% of the stand with an additional cruise of the remaining area for large-diameter trees. Finally, the high variability observed in these even-aged stratified mixed stands indicates that structural complexity may not be well estimated from the few subplots of small extent typically employed by national inventory programs designed to capture broad trends in forest extent and condition.

## Data availability

The raw data are available through the Penn State Data Commons at https://doi.org/10.26208/wx8r-fb53.

## References

Aldrich PR, Parker GR, Ward JS, Michler CH (2003) Spatial dispersion of trees in an old-growth temperate hardwood forest over 60 years of succession. For Ecol Manage 180:475–491

Ashton MS, Kelty MJ (2018) The Practice of Silviculture: Applied Forest Ecology, 10th edn. John Wiley and Sons, Hoboken, NJ

Avery TE, Burkhart HE (2001) Forest measurements, 5th edn. McGraw-Hill, New York

Bate LJ, Garton EO, Wisdom MJ (1999) Estimating snag and large tree densities and distributions on a landscape for wildlife management. USDA For Ser PNW-GTR-425. Portland, OR

Bäuerle H, Nothdurft A, Kändler G, Bauhus J (2009) Monitoring habitat trees and coarse woody debris based on sampling schemes. Allg Forst u Jagdztg 180:249–260

Becker P, Nichols T (2011) Effects of basal area factor and plot size on precision and accuracy of forest inventory estimates. N J Appl For 28:152–156

Berrill JP, O'Hara KL (2012) Influence of tree spatial pattern and sample plot type and size on inventory. In: Standiford RB, Weller TJ, Piirto DD, Stuart JD (cords) Proceedings of coast redwood forests in a changing California: a symposium for scientists and managers. USDA For Ser PSW-GTR-238. Albany, CA, pp 485–497

Boots BN, Getis A (1988) Point pattern analysis. SAGE Publications, Newbury Park, CA

Bormann FH (1953) The statistical efficiency of sample plot size and shape in forest ecology. Ecology 34:474–487

Brown S, Schroeder P, Birdsey R (1997) Aboveground biomass distribution of US eastern hardwood forests and the use of large trees as an indicator of forest development. For Ecol Manage 96:37–47

Busing RT, White PS (1993) Effects of area on old-growth forest attributes: implications for the equilibrium landscape concept. Lands Ecol 8:119–126

Carrer M, Castagneri D, Popa I, Pividori M, Lingua E (2018) Tree spatial patterns and stand attributes in temperate forests: the importance of plot size, sampling design, and null model. For Ecol Manage 407:125–134

Clark PJ, Evans FC (1954) Distance to nearest neighbors: a measure of spatial relationships in populations. Ecology 35:445–453

Crow TR, Buckley DS, Nauertz EA, Zasada JC (2002) Effects of management on the composition and structure of northern hardwood forest in Upper Michigan. For Sci 48:129–145

Du J, Zhao W, He Z, Yang J, Chen L, Zhu X (2018) Characterizing stand structure in a spruce forests: effects of sampling protocols. Sci Cold Arid Reg 7:245–256

Fonton NH, Atindogbe G, Hounkonnou NM, Dohou RO (2011) Plot size for modelling the spatial structure of Sudanian woodland trees. Ann For Sci 68:1315–1321

Franklin JF, Van Pelt R (2004) Spatial aspects of structural complexity in old-growth forests. J For 102:22–28

Fraser AR, van den Driessche P (1972) Triangles, density, and pattern in point populations. Proceeding of 3rd Conference of Advisory Group of Forest Statisticians International Union of Forest Research Organisations. INRA, Jouy-en-Josas, pp 277–286

Freese F (1967) Elementary statistic methods for foresters. USDA Forest Service Agricultural Handbook 317, Madison WI. https://www.fpl.fs.fed.us/documnts/usda/ah317.pdf. Accessed 18.08.20

Füldner K (1995) Zur Strukturbeschreibung in Mischbeständen Forstarchiv 66:235–240

Gormanson DD, Pugh SA, Barnett CJ, Miles PD, Morin RS, Sowers PA, Westfall JA (2018) Statistics and quality assurance for the Northern Research Station Forest Inventory and Analysis Program. USDA For Ser NRS-GTR-178. Newtown Square, PA

Griffith DA (1983) The boundary value problem in spatial statistical analysis. J Reg Sci 23:377–387

Gray A (2003) Monitoring stand structure in mature coastal Douglas-fir forests: effect of plot size. For Ecol Manage 175:1–16

Grotti M, Chianucci F, Puletti N, Fardusi MJ, Castaldi C, Corona P (2019) Spatio-temporal variability in structure and diversity in a semi-natural mixed oak-hornbeam floodplain forest. Ecol Indic 104:576–587

Guillemette F, Lambert MC, Bédard S (2012) Sampling design and precision of basal area growth and stand structure in uneven-aged northern hardwoods. For Chron 88:30–39

Hagner M, Nyquist H (1998) A coefficient for describing size variation among neighboring trees. J Agr Biol Env Stat 3:62–74

Jazbec A, Vedriš M, Božić M, Goršić E (2011) Efficiency of inventory in uneven-aged forests on sample plots with different radii. Croat J For Eng 32:311–312

Johnson SE, Mudrak EL, Beever EA, Sanders S, Waller DM (2008) Comparing power among three sampling methods for monitoring forest vegetation. Can J For Res 38:143–156

Kane VR, McGaughey RJ, Bakker JD, Gersonde RF, Lutz JA, Franklin JF (2010) Comparisons between field-and LiDAR-based measures of stand structural complexity. Can J For Res 40:761–773

Kukunda CB, Beckschäfer P, Magdon P, Schall P, Wirth C, Kleinn C (2019) Scale-guided mapping of forest stand structural heterogeneity from airborne LiDAR. Ecol Indic 102:410–425

Kenkel NC, Podani J (1991) Plot size and estimation efficiency in plant community studies. J Veg Sci 2:539–544

Kandare K (2017) Fusion of airborne laser scanning and hyperspectral data for predicting forest characteristics at differernt spatial scales. Dissertation, Norwegian University of Life Sciences, Ås, Norway. https://static02.nmbu.no/mina/forskning/drgrader/2017-Kandare.pdf. Accessed 18.08.20

Král K, Janík D, Vrška T, Adam D, Hort L, Unar P, Šamonil P (2010) Local variability of stand structural features in beech dominated natural forests of Central Europe: Implications for sampling. For Ecol Manage 260:2196–2203

Larson BC (2007) Thoughts on the development of new, appropriate measures to complexity. In: Complex stand structures and associated dynamics: measurement indices and modelling approaches. Ontario For Res Inst For Res Info Paper No. 67, Sault Ste. Marie, Ontario, pp 54–55

Lombardi F, Marchetti M, Corona P, Merlini P, Chirici G, Tognetti R, Burrascano S, Alivernini A, Puletti N (2015) Quantifying the effect of sampling plot size on the estimation of structural indicators in old-growth forest stands. For Ecol Manage 346:89–97

Maleki K, Kiviste A (2015) Effect of sample plot size and shape on estimates of structural indices: a case study in mature silver birch (

*Betula pendula*Roth) dominating stand in Järvselja. For Studies 63:130–150McCune B, Lesica P (1992) The trade-off between species capture and quantitative accuracy in ecological inventory of lichens and bryophytes in forests in Montana. Bryol 95:296–304

McGee GG, Leopold DJ, Nyland RD (1999) Structural characteristics of old-growth, maturing, and partially cut northern hardwood forests. Ecol Appl 9:1316–1329

Meng J, Li S, Wang W, Liu Q, Xie S, Ma W (2016) Estimation of forest structural diversity using the spectral and textural information derived from SPOT-5 satellite images. Rem Sens 8:125–149

Mura M, McRoberts RE, Chirici G, Marchetti M (2015) Estimating and mapping forest structural diversity using airborne laser scanning data. Rem Sens Env 170:133–142

Nagel LM, Janowiak MK, Webster CR (2007) Spatial scale affects diameter distribution shape in uneven-aged northern hardwoods. In: Complex stand structures and associated dynamics: measurement indices and modelling approaches. Ont For Res Inst For Res Info Paper No. 67, Sault Ste. Marie, Ontario, pp 41–43

Packard KC, Radtke PJ (2007) Forest sampling combining fixed-and variable-radius sample plots. Can J For Res 37:1460–1471

Peck JE, Zenner EK (2021) Structure Area Curves in Eastern Hardwoods: implications for minimum plot sizes to capture spatially-explicit structure indices. [dataset]. V1. Penn State Data Commons repository. https://doi.org/10.26208/wx8r-fb53

Proulx R (2007) Ecological complexity for unifying ecological theory across scales: a field ecologist’s perspective. Ecol Compl 4:85–92

Reich RM, Arvanitis LG (1992) Sampling unit, spatial distribution of trees, and precision. N J Appl For 9:3–6

Ritter T, Saborowski J (2014) Efficient integration of a deadwood inventory into an existing forest inventory carried out as two-phase sampling for stratification. Forestry 87:571–581

Rubin BD, Manion PD, Faber-Langendoen D (2006) Diameter distributions and structural sustainability in forests. For Ecol Manage 222:427–438

Scott C (1998) Sampling methods for estimating change in forest resources. Ecol Appl 8:228–233

Spies TA (2004) Ecological concepts and diversity of old-growth forests. J For 102:14–20

Spies TA, Franklin JF (1991) The structure of natural young, mature, and old-growth Douglas-fir forests in Oregon and Washington. In: Spies TA, Franklin JF (eds) Wildlife and vegetation of unmanaged Douglas-fir forests, USDA Forest Service PNW-GTR-285. Portland OR, pp 91–109

Tinkham WT, Mahoney PR, Hudak AT, Domke GM, Falkowski MJ, Woodall CW, Smith AM (2018) Applications of the United States Forest Inventory and Analysis dataset: a review and future directions. Can J For Res 48:1251–1268

Thompson WL, Burnham KP (2004) Sampling rare or elusive species: concepts, designs, and techniques for estimating population parameters. Island Press, Washington DC

Winter S, Chirici G, McRoberts RE, Hauk E, Tomppo E (2008) Possibilities for harmonizing national forest inventory data for use in forest biodiversity assessments. Forestry 81:33–44

Zeide B (1980) Plot size optimization. For Sci 26:251–257

Zenner EK (2005) Investigating scale-dependent stand heterogeneity with structure-area-curves. For Ecol Manage 209:87–100

Zenner EK, Hibbs D (2000) A new method for modeling the heterogeneity of forest structure. For Ecol Manage 129:75–87

Zenner EK, Peck JE (2009) Characterizing structural conditions in mature managed red pine: spatial dependency of metrics and adequacy of plot size. For Ecol Manage 257:311–320

Zenner EK, Peck JE (2018) Floating neighborhoods reveal contribution of individual trees to high sub-stand scale heterogeneity. For Ecol Manage 412:29–40

Zenner EK, Peck JE, Hobi ML, Commarmot B (2015) The dynamics of structure across a primeval European beech stand. Forestry 88:180–189

Zenner EK, Peck JE, Sagheb-Talebi K (2018) One shape fits all, but only in the aggregate: diversity in sub-stand scale diameter distributions. J Veg Sci 29:501–510

Zenner EK, Peck JE, Sagheb-Talebi K (2019) Patchiness in old-growth oriental beech forests across development stages at multiple neighborhood scales. Eur J For Res 138:739–752

Zenner EK, Peck JE, Hobi ML (2020) Development phase convergence across scale in a primeval European beech (

*Fagus sylvatica*L.) forest. For Ecol Manage 460: 117889

## Acknowledgements

Assistance in the field and with data entry was provided by Dan Heggenstaller, Darren Wolfgang, and Jeff Watson, and Joe Harding provided site information.

## Funding

Funding was provided by the Department of Ecosystem Science and Management, the College of Agricultural Sciences (CAS), a CAS Seed Grant from The Pennsylvania State University, and USDA National Institute of Food and Agriculture Hatch Appropriations [#PEN04639, Accession #1015105].

## Author information

### Authors and Affiliations

### Corresponding author

## Ethics declarations

### Conflict of interest

The authors declare that they have no conflict of interest.

## Additional information

**Handling Editor: **Andreas Bolte

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Contributions of the co-authors **Eric Zenner: Conceptualization

Jeri Peck, Eric Zenner: Methodology

Eric Zenner, Jeri Peck: Formal Analysis

Jeri Peck, Eric Zenner: Writing

## Appendix

### Appendix

## Rights and permissions

## About this article

### Cite this article

Peck, J., Zenner, E. Structure area curves in Eastern Hardwoods: implications for minimum plot sizes to capture spatially explicit structure indices.
*Annals of Forest Science* **78**, 16 (2021). https://doi.org/10.1007/s13595-021-01036-5

Received:

Accepted:

Published:

DOI: https://doi.org/10.1007/s13595-021-01036-5

### Keywords

- Sampling
- Rarity
- Estimation error
- Oak
- Quercus
- Structural complexity