Basal area has shortcomings as a measure of stand density, but it is often preferred for operational assessments because it is easy to measure. Previous work has demonstrated that an additive version of Reineke's stand density index can be estimated by a simple tree count using a modified horizontal point sampling technique. We show that this technique can be extended further to estimate a mixed-species density measure that has been developed for complex stands in the northeastern United States, using wood specific gravity to harmonize the density contributions of different species. The sampling technique provides design-unbiased estimates of stand density from a weighted tree count, where the weights depend on specific gravity but not on diameter. Rounding the specific gravity values for different species in the calculation of estimates introduces a trivial amount of bias but streamlines the procedure for rapid use in the
field.
1. Introduction
The ability to measure forest stand density is fundamental for both description and prescription in applied ecology including silviculture, forest health, and wildlife habitat management [1, 2]. Basal area is the most commonly used density metric in practice, primarily because it is very fast and simple to measure in the field [3–5]. However, the basal area at which key ecological behaviors occur (such as crown closure, the onset of self-thinning, and asymptotic density) depends on the size of the trees in the stand, the species of those trees, and also potentially on the age and site quality of the stand. Other metrics, such as Reineke’s [6] Stand Density Index (SDI) and related metrics such as those of Curtis [7] and Long and Daniel [8], offer improvements over basal area for monospecific stands. Reineke’s [6] SDI is
SDI=N(Dq25)1.6,
where N is number of trees per hectare, and Dq is stand quadratic mean diameter. SDI has units of trees/ha. Zeide ([2, page 5]) writes that SDI is “still the best [metric] we have” for assessing stand density.
Recent research efforts have attempted to extend Reineke’s [6] basic work to more complex stands than the even-aged, monospecific subjects of the original SDI. Long and Daniel [7] developed an additive version of SDI closely related to an index explored by Curtis [8]:
ASDI=ΣiNi(Di25)1.6,
where Ni and Di are the number of trees per hectare and diameter of size classes of trees, and the summation is taken over an exhaustive, mutually-exclusive set of classes. (For continuous diameter distributions, the summation in (2) can be replaced by an integral and Ni can be replaced by the product of N and an appropriate probability density function; see [9, 10], e.g.) Several authors [1, 7, 11] have suggested using ASDI for multicohort stands. More recently, Woodall et al. [12, 13], building on earlier mechanistic work [14], suggested that the maximum ASDI of mixed species stands could be estimated as a function of the mean specific gravity of the wood of trees on the plots and proposed a modified density measure. Ducey and Knapp [15], building on [12], suggested a relative density metric of the form
RD=ΣijNij(b0+b1SGj)(Di25)1.6,
where Nij is the number of trees per hectare in the ith diameter class of the jth species, b0 and b1 are parameters to be estimated from plot data, and SGj is the specific gravity of wood for the jth species.
A challenge for any stand density metric more complex than basal area is that of obtaining a fast and simple estimate in the field. That challenge seems more and more difficult as stand density metrics such as that in (3) attempt to encompass stands that are more complex in diameter distribution and species mixture. Rapid field estimates are particularly valuable for diagnosis and description of stand conditions when a density metric will be used in situ in conjunction with expert opinion, and for controlling the marking of thinnings, shelterwoods, single-tree selection, and other partial harvests. Basal area would seem to have a distinct advantage in these situations, because a simple count of trees using a prism or other angle gauge, multiplied by the basal area factor of that gauge, provides a direct estimate of basal area per hectare [16, 17]. However, a modified prism technique can provide estimates of ASDI through a simple count as well [18].
The goal of this study is to show how the modified prism technique [18] can be modified to provide direct estimates of the Ducey and Knapp [15] relative density for mixed species stands. We provide tables and instructions for implementing the approach and illustrate its application with data from a mixed-species stand in southeastern New Hampshire.
2. Theory2.1. Mixed-Species Density Measure
Ducey and Knapp [15] proposed RD as a stand density metric for the mixed-species forests of New England and New York in the USA. They used data from the Forest Inventory and Analysis program of the U.S. Forest Service to estimate b0=0.00015 and b1=0.00218, using specific gravity at 12% moisture content for each species. RD is a unitless fraction, scaled so that RD=0 represents a completely empty stand, and RD=1 represents a stocking level consistent with “normal” or A-line stocking using conventional stocking guides. For example, if RD=0.5 for a particular stand, that stand has 50% of the stocking level implied by “normal” or A-line stocking.
Now, consider (3). It can be rewritten as
RD=Σj(b0+b1SGj)[ΣiNij(Di25)1.6]=Σj(b0+b1SGj)ASDIj
which makes clear that RD is simply a weighted form of ASDI, where the weights for different species are determined by their specific gravity.
Ducey and Valentine [18] developed a technique analogous to sampling for basal area with a prism, in which tree count times a factor gives ASDI. We will show that it is possible to modify that technique to give RD as a sum of weights for the tallied trees. First, we review the modified prism technique for ASDI.
2.2. Direct Sampling Technique
Recall that in ordinary horizontal point sampling (HPS) with a prism, relascope, or other angle gauge, each tallied tree represents the same amount of basal area per hectare [16]. That quantity is called the basal area factor or BAF, and that quantity is constant for a given gauge because the area of the inclusion zone (the area within which a specific tree will be tallied) is directly proportional to the basal area of the individual tree. This result follows from the geometry of the gauge: the radius of the inclusion zone (also called the limiting distance) equals the gauge constant (or ratio of length to width) times the diameter of the tree. Mathematical proofs are provided by [19], while [20] provides a thorough discussion of HPS in its many applications.
If the inclusion zone radius could be made proportional to D0.8, rather than proportional to D, then inclusion zone area would be proportional to D1.6 and each tallied tree would represent a constant amount of ASDI [18]. This quantity is called the stand density factor (SDF), by analogy to the BAF of ordinary HPS. Recall that ASDI has units of trees/ha; so SDF also has units of trees/ha. If the method of [18] is used, and if m trees are tallied at n points, then the estimate of ASDI is simply
ASDI=(mn)SDF.
If the SDF is chosen with reasonable care, then the limiting distances for a given SDF are very close to those for the BAF that would ordinarily be used in similar stand conditions [18] (Figure 1). The slight departure, depending only on tree diameter, could easily be adjusted in the field by “pushing the point” or “pulling the point.” Pushing the point [21] consists of holding the prism not directly over the sample point but extended toward the tree of interest. In ordinary HPS pushing the point is a common field error, and it creates biased estimates. Pulling the point consists of holding the prism back and sighting across the sample point to determine whether a tree should be tallied. A simple table of the distance to push or pull the point, in relation to D for individual trees, is adequate to provide a sample with a constant SDF [18]. In practice, it is not usually necessary to know the exact D of an individual tree to determine whether it should be included in the sample. Because the degree of push or pull changes slowly with D, an initial ocular estimate is usually sufficient. D should be measured on borderline trees, and if an estimate of ASDI is desired for different diameter classes, then D can be measured on the tallied trees as well.
Area of the inclusion zone and limiting distance (or inclusion zone radius), for ordinary HPS with a BAF 4 m2/ha gauge and direct sampling for ASDI with a stand density factor of 50.
Now suppose that we wish to modify the technique to estimate RD, building on (4). Let mj be the number of trees of the jth species tallied, using n sample points as before. Then
ASDIj=(mjn)SDF
provides an estimate of the ASDI of species j, and
(b0+b1SGj)ASDIj=(b0+b1SGj)(mjn)SDF
provides a natural estimator of the summand in (4). Rearranging slightly, and substituting, we have
RD=Σj(mjn)[(b0+b1SGj)SDF]
as a natural estimator for RD. In other words, RD can be estimated as a weighted tree count, in which each tree at a single sample point represents [(b0+b1SGj)SDF] units of RD. Recall that RD is a unitless fraction; therefore the weight or contribution of an individual tree [(b0+b1SGj)SDF] is also a unitless fraction. Proofs of unbiasedness of (5), (6), and (8) are given in the appendix.
In practice, if “exact” (or at least, unrounded) values of SGj are used, the quantity [(b0+b1SGj)SDF] will likely be a cumbersome fraction. This is not a problem if, for example, field computers are used to provide nearly instantaneous feedback on sample results. However, for the cruiser or timber marker who desires a simple, rapid field procedure, it may be more convenient to group species and use an approximate value so that [(b0+b1SGj)SDF] is always a nearly round percentage. This requires advance tabulation of species into groups, which we undertake hereinafter. Furthermore, it does introduce a small amount of bias due to rounding error. As we show, this is likely to be negligible for samples taken at typical operational intensities.
3. Application
One of the first steps in direct sampling for ASDI is to choose a BAF and corresponding SDF [18]. In the northeastern United States, a BAF of 2.3 m2/ha is a common choice but can lead to nondetection errors; so a larger BAF is desirable [22, 23]. Prisms with a BAF of 4 m2/ha are widely available and fall within the range recommended for the region [24]. For a given BAF, the SDF should be chosen such that the limiting distances are similar. This allows the push or pull to be done quickly and accurately. A good match to BAF 4 m2/ha is an SDF of 50 trees/ha: the limiting distances differ by less than 2 m over a broad range of DBH, and “pulling the point” is only necessary for rare trees with DBH>85 cm. The limiting distances and push/pull for an SDF 50 using a BAF 4 are shown in Table 1.
Critical distances and required “push’’ for sampling with a stand density factor (SDF) of 50, using a basal area factor (BAF) 4 m2/ha prism. When the “push’’ is shown as negative, the prism must be pulled back across the point.
DBH, cm
BAF=4 Critical Distance, m
SDF=50 Critical Distance, m
Required Push, m
5
1.25
2.20
0.95
10
2.50
3.83
1.33
15
3.75
5.30
1.55
20
5.00
6.67
1.67
25
6.25
7.98
1.73
30
7.50
9.23
1.73
35
8.75
10.44
1.69
40
10.00
11.62
1.62
45
11.25
12.77
1.52
50
12.50
13.89
1.39
55
13.75
14.99
1.24
60
15.00
16.07
1.07
65
16.25
17.14
0.89
70
17.50
18.18
0.68
75
18.75
19.21
0.46
80
20.00
20.23
0.23
85
21.25
21.24
-0.01
90
22.50
22.23
-0.27
95
23.75
23.21
-0.54
100
25.00
24.19
-0.81
Ducey and Knapp [15] tabulate SG for all the species found in New York and New England in the USFS FIA database. However, for fast application in the field it is useful to tabulate rounded values of [(b0+b1SGj)SDF] as weights. This is done in Table 2, with weights multiplied by 100 so that they reflect percents rather than fractions.
Rounded weights (values of 100×[(b0+b1SGj)SDF]) for common species in New York and New England, using an SDF of 50.
Acer platanoides, Acer saccharum, Betula allegheniensis, Betula lenta, Betula nigra, Carya laciniosa, Crataegus spp., Fagus grandifolia, Gleditsia tricanthos, Malus spp., Quercus alba and the other white oaks (except bicolor), Quercus rubra and the other red oaks (except velutina and ilicifolia), Robinia pseudoacacia, Ulmus thomasii
Chamaecyparis thyoides, Populus balsamifera, Thuja occidentalis
0.31–0.34
4
Armed with Tables 1 and 2, it is straightforward to conduct field sampling and to calculate estimates. For example, suppose that we have used our prism and Table 1 to determine which trees to sample at a single point. We tallied 7 trees: 3 Betula papyrifera, 1 Picea rubens, 2 Abies balsamea, and 1 Fagus grandifolia. We can immediately determine that the ASDI at this point is 7×50=350, and RD is 3×7+1×5+2×5+1×8=44%, or 0.44.
4. Field Demonstration
To illustrate the calculations associated with the method, we conducted sampling in a 4.3 ha mixed-species stand located in Durham, NH. 12 sample points were located systematically on a 60 m by 60 m grid. We used a BAF 4 m2/ha prism, and the push-pull values from Table 1, to obtain a sample of trees with an SDF of 50 trees/ha. Both the exact weights as calculated from the specific gravities tabulated in [15], and rounded weights from Table 2, were used to estimate RD for comparison.
A total of 133 trees were tallied on the 12 plots. Of these, 49 were Tsuga canadensis, 33 were Pinus strobus, 31 were Acer rubrum, 10 were Quercus rubra, 5 were Fagus grandifolia, 4 were Betula allegheniensis, and 1 was Quercus alba. RD as calculated using the exact weights was 65.6%, with a standard error of 7.6%; the coefficient of variation of the sample points was 40%. RD as calculated using the rounded weights was 63.3%, with a standard error of 7.7%; the coefficient of variation of the sample points was 42%. In this stand at least, the difference caused by rounding the weights is trivial in comparison with the sampling variability of RD.
5. Discussion and Conclusions
Metrics of stand density can be viewed in two ways: as hypotheses about or models of stand dynamics and as operational aids to silviculture. From the first perspective, stand basal area has long been rejected. Even in single-species, single-cohort stands, the inadequacy of basal area for comparing density in stands of different tree size and species was known long before Reineke’s [6] work and provided much of the impetus for his efforts. Even in those stands, Reineke’s work represents a simplification, obscuring potentially important interspecific differences in allometry and growth [25]. Extension of Reineke’s ideas to multicohort stands and other stands with complex structure has been an important area of recent work [7, 26–28]. While empirical tests of those efforts in comparison to previous measures have been relatively rare [9, 29], some progress has been made at least in understanding how one might approach the design of experiments [10].
A variety of studies have explored adapations of SDI to mixed-species stands [12, 13, 30–36]. However, the Ducey and Knapp [15] formulation has two properties that are both required for the straightforward use of the sampling method developed here. The first property is additivity [9]. An index is said to be additive if, for two stand components A and B, the density of A and B considered together equals the density of A considered alone, plus the density of B considered alone. Additivity is required for unbiased estimation of any stand density metric using traditional methods (such as fixed-area plot sampling or HPS, or any method that can be identified as using a Horvitz-Thompson [37] estimator) [38, 39]. Among previous mixed-species indices, only that of Stout and colleagues [31, 32] is additive; that of Woodall et al. [12, 13] meets the criterion of separability but not additivity [9]. By analogy with site index curves, we may also say that a density measure is anamorphic if, for two different species, the density contribution of an individual of species X divided by that of an individual of the same D but of species Y is a constant that does not depend on D. The density measure of Stout and colleagues [12, 13] is not anamorphic. If a density measure is not anamorphic, then a different “push” or “pull” would be required for each species or species group. Thus, the Ducey and Knapp [15] formulation has certain practical advantages from a measurement standpoint. We caution, however, that simplicity and practicality may be advantageous but do not imply biological correctness [9]. Both additivity and anamorphism are likely incorrect from a strictly biological and mechanistic standpoint [9, 25, 40]. On the other hand, detecting departures from additivity and anamorphism implies the availability of an additive, anamorphic measurement from which to construct a baseline.
At its heart, a stand density index is a model, hypothesis, or conjecture about the relationship between measurable quantities and ecological behaviors. Statistician G.E.P. Box [41, page 2] wrote “Models, of course, are never true, but fortunately it is only necessary that they be useful.” The durability of basal area as a practical measure of stand density points to the importance of rapid, simple measurement as a component of utility. Box [42, page 792] also wrote “Since all models are wrong the scientist must be alert to what is importantly wrong. It is inappropriate to be concerned about mice when there are tigers abroad.” One tiger abroad in stand density measurement is sampling variability. Fast, repeatable measurement with relatively low sampling variability has been a reality for basal area since the seminal work of Bitterlich [16] on HPS. However, the recent generalization of HPS for ASDI [18], and the approach for RD as explored here, may close the gap somewhat.
In theory, the modified HPS approach [18] and this method could be used as a general sampling approach, if tree diameters are measured. However, the need to push (or pull) the point does open the door to field errors. The results from ordinary HPS should be more precise and less prone to operator error, even when a final variable of interest is RD. Rather, the strength of the approach we present here is for near-instantaneous estimation of stand density in three primary situations.
(1) The first situation is to allow quick assessment and description of stand conditions. Fast calculation of RD would allow more effective coupling between quantitative measurement and subjective, expert opinion [18]. Density indices and related diagrams are best viewed as complementary to qualitative silvicultural expertise [43]; a rapid measurement technique can make quantitative information available to the forester or ecologist on-site.
(2) The second one is to enable control of residual density while marking partial harvests such as thinnings and shelterwoods.
(3) The third one is to provide a simple point density estimate for use in ecological studies that would be available immediately for stratification or subsampling purposes.
Equation [7] provides an estimating equation for RD when trees are sampled using the modified point sampling technique. In this appendix, we present a formal proof of design-unbiasedness for this estimating equation.
Suppose that we are interested in the population of trees occurring in a tract of area A. Let the trees be indexed by k=1,…,K, where K is the total number of trees. K is fixed but unknown, as are the diameters Dk of the trees at the outset of sampling. Without loss of generality, let us locate a single point at random and with uniform density within the tract. This point will be our sample point, and from it we will include all trees located at a distance no greater than cD0.8, where c is a constant chosen in advance. Denote the inclusion area of the kth tree by ak; ak is the area where the sample point may fall and the kth tree will be included. From basic geometry we have
ak=π(cD0.8)2=πc2Dk1.6,
and the probability pk that the kth tree will be included in the sample (ignoring boundary slopover, which can easily be corrected using the walkthrough method [44]) is
pk=akA.
We can rewrite (2) as
ASDI=(1A)Σk(Dk25)1.6.
Now let dk be an indicator variable taking the value 1 if the kth tree is included in the sample, and 0 otherwise; the expectation is E[dk]=pk. The total number of trees included in the sample is
m=Σkdk,
and its expectation is
E[m]=E[Σkdk]=Σkpk=(πc2A)ΣkDk1.6=[πc225-1.6]ASDI
from which it follows directly that if SDF=251.6/(πc2), then mSDF is a design-unbiased estimator of ASDI, since
E[mSDF]=SDFE[m]=ASDI,
and if we collect data at n such sample points, the sample mean (m/nSDF) will also provide a design-unbiased estimate. This proves the unbiasedness of the estimator in (5). Furthermore, if we restrict the population to consist only of trees of species j and replace m with mj and ASDI with ASDIj in the above equations, the equations still hold; this proves the unbiasedness of the estimator in (6).
In general, if there is an attribute yk associated with the kth tree, an estimator of the form
Y=Σkdkykpk
is a Horvitz-Thompson [37] estimator and provides a design-unbiased estimate of the population total Σkyk. Now, consider (4), rewritten as a population total:
RD=(1A)(b0+b1SGk)(Dk25)1.6,
where SGk is the specific gravity of the kth tree, so that
yk=(1A)(b0+b1SGk)(Dk25)1.6.
Now, if we substitute this into the basic Horvitz-Thompson estimator, we obtain
Y=Σkdk[(1/A)(b0+b1SGk)(Dk/25)1.6]pk=Σkdk[(1A)(b0+b1SGk)(Dk25)1.6]/[πc2Dk1.6/A]=Σkdk[(b0+b1SGk)251.6πc2]=SDFΣkdk(b0+b1SGk).
Now, recognizing that all trees of the jth species have the same specific gravity for our purposes, and letting SGj represent the specific gravity of the jth species as before, we see that
Y=Σjmj[(b0+b1SGj)SDF]
is the design-unbiased Horvitz-Thompson estimator associated with a single sample point, and the mean of several sample points is identical to (8), which must also be design-unbiased.
Acknowledgments
Colleen Didas assisted with the field work. Development of the relative density index was supported under a grant from the National Science Foundation, “Collaborative Research: The northeastern carbon sink: enhanced growth, regrowth, or both?” to M. J. Ducey. This paper is a scientific contribution of the New Hampshire Agricultural Experiment Station.
LongJ. N.A practical approach to density management19856112327ZeideB.How to measure stand density200519111410.1007/s00468-004-0343-x2-s2.0-12544252732SpurrS. H.1952New York, NY, USARonald PressBickfordC. A.BakerF. S.WilsonF. G.Stocking, normality, and measurement of stand density195755299104AveryT. E.BurkhartH. E.20025thNew York, NY, USAJohn Wiley & SonsReinekeL. H.Perfecting a stand-density index for even-aged forests1933467627638CurtisR. O.A tree area power function and related density measures for Douglas-fir1971172146159LongJ. N.DanielT. W.Assessment of growing stock in uneven-aged stands1990539396DuceyM. J.LarsonB. C.Is there a correct stand density index? An alternate interpretation20031831791842-s2.0-0042665477DuceyM. J.The ratio of additive and traditional stand density indices20092415102-s2.0-63349108847FiedlerC. E.CullyJ. E.A silvicultural approach to develop Mexican spotted owl habitat in Southwest forests1995104144148WoodallC. W.MilesP. D.VissageJ. S.Determining maximum stand density index in mixed species stands for strategic-scale stocking assessments20052161–336737710.1016/j.foreco.2005.05.0502-s2.0-24944572233WoodallC. W.PerryC. H.MilesP. D.The relative density of forests in the United States20062261–336837210.1016/j.foreco.2006.01.0322-s2.0-33645863129DeanT. J.BaldwinV. C.Jr.The relationship between Reineke's stand-density index and physical stem mechanics1996811–325342-s2.0-0002450090DuceyM. J.KnappR. A.Utilizing quantile regression to develop a stand density index for mixed species stands in the northeastern United StatesCanadian Journal of Forest Research. In reviewBitterlichW.Die Winkelzählprobe1948591-245GrosenbaughL. R.Plotless timber estimates—new, fast, easy19525013237DuceyM. J.ValentineH. T.Direct sampling for stand density index200823278822-s2.0-43149110210GrosenbaughL. R.Point sampling and line sampling: probability theory, geometric implications, synthesis1957160Asheville, NC, USAUSDA Forest Service, Southern Forest Experiment StationBitterlichW.1984Slough, UKCommonwealth Agricultural Bureaux, Farnham RoyalOderwaldR. G.GregoireT. G.Overstated volumes from ‘pushing the point’1995194162165WenselL. C.LevitanJ.BarberK.Selection of basal area factor in point sampling19807828384WiantH. V.Jr.YandleD. O.AndreasR.Is BAF 10 a good choice for point sampling?1984122325DuceyM. J.Pre-cruise planningWorkshop Proceedings: Forest Measurements for Natural Resource Professionals. II. Getting the Most from Your Cruise2001Durham, NH, USAUniversity of New Hampshire Cooperative ExtensionPretzschH.BiberP.A re-evaluation of Reineke's rule and stand density index20055143043202-s2.0-23244438998LongJ. N.O'HaraK.Using stand density index to regulate stocking in uneven-aged stands199556Montana Forest and Conservation Experiment Station110122LongJ. N.A technique for the control of stocking in two-storied stands199611259612-s2.0-0002003070ShawJ. D.Application of stand density index to irregularly structured stands200015140422-s2.0-0002896591WoodallC. W.FiedlerC. E.MilnerK. S.Stand density index in uneven-aged ponderosa pine stands20033319610010.1139/x02-1682-s2.0-0037282357StoutS. L.NylandR. D.Role of species composition in relative density measurement in Allegheny hardwoods19861635745792-s2.0-0022834031StoutS. L.MarquisD. A.ErnstR. L.A relative density measure for mixed-species stands198785745472-s2.0-0023528605PuettmannK. J.HibbsD. E.HannD. W.The dynamics of mixed stands of Alnus rubra and
Pseudotsuga menziesii: extension of size-density analysis to species mixture19928034494582-s2.0-0027063366PuettmannK. J.HannD. W.HibbsD. E.Evaluation of the size-density relationships for pure red alder and Douglas-fir stands1993391727CochranP. H.GeistJ. M.ClemensD. L.ClausnitzerR. R.PowellD. C.Suggested stocking levels for forest stands in northeastern Oregon and southeastern Washington1994PNW-RN-513Portland, Ore, USAU.S. Department of Agriculture Forest Service, Pacific Northwest Research StationTorres-RojoJ. M.MartinezA.Relative stand density index for mixed even-aged stands20004497507WilliamsR. A.Use of stand density index as an alternative to stocking percent in upland hardwoods20032031371422-s2.0-4644335507HorvitzD. G.ThompsonD. J.A generalization of sampling without replacement from a finite universe1952472608996DuceyM. J.LarsonB. C.Thinning decisions using stand density indices: the influence of uncertainty199712389922-s2.0-0042111485DuceyM. J.mjducey@cisunix.unh.eduLarsonB. C.bruce.larson@.yale.eduAccounting for bias and uncertainty in nonlinear stand density indices19994534524572-s2.0-0028978761PretzschH.Species-specific allometric scaling under self-thinning: evidence from long-term plots in forest stands2006146457258310.1007/s00442-005-0126-02-s2.0-29144535281BoxG. E. P.Some problems of statistics and everyday life19727436514BoxG. E. P.Science and statistics197671356791799LongJ. N.ShawJ. D.A density management diagram for even-aged ponderosa pine stands20052042052152-s2.0-27644487669DuceyM. J.GoveJ. H.ValentineH. T.A walkthrough solution to the boundary overlap problem20045044274352-s2.0-4644340892