The size and form of sampling units-SU have always been variables considered in planning and structuring forest inventories, being performed in forests or in plantations. The experimental work outlined to deal with the problem was conducted in an area of araucaria forest, in Sao Joao de Triunfo, PR, Brazil. The forms of sampling units circle, square and rectangular were evaluated, whose areas ranged from 200 m2 to 1000 m2. Time was recorded using a stopwatch and computed separately for locomotion and measurement. The power model was used to adjust the relations of the times of locomotion and measurement, as well as a hyperbolic one for the coefficient of variation, all taking as function of the SUs sizes. To achieve the analytical solution for the optimum size of the SU, it was necessary simulating the behavior of the three functions until the size of 10,000 m2. By taking the derivative of the combined function it was found the maximum point, which allowed optimizing the size of the SU in 600 m2 for the structured experimental conditions. This result proved the formulated hypothesis, performing even a critical analysis of the inclusion of other relevant variables, such as size of the area to be inventoried, number of SUs performed in one day of work, average distance between the SUs, average speeds for locomotion between SUs and for taking all respective measurements.
Forest Inventories were and still are being informative instruments essential for the planning and forest man- agement, being performed in plantations or in native forests. The information is usually obtained by sampling and this is taken in a sample space previously defined.
The sampling unit-SU, as it is specified a priori in its shape and size, characterizes a selection with probability proportional to the area and frequency of individuals who comprise it, and is called a method of Fixed Area.
It will be considered three relevant aspects that directly affect the size and shape of SUs: 1) The degree of va- riability of the experimental variable used to express the SU, usually the volume; 2) The costs or time spent to obtain the information in the SUs and 3). The degree of effectiveness to accomplish the work of data collection in minimal time.
It is known, of several published works in the field of forest inventory, that the shape and size of SUs has been decided much more by convenience and the operability of their location and demarcation in the field, than by analytical criterion (Péllico Netto & Brena, 1997; Silva, Xavier, Leite, & Pires, 2003). Whilst it is extremely valid worrying about the more appropriate conditions to sample forests or plantations, it will be equally relevant worrying about optimizing the time and the cost of the measurements as well as when conducting other activities in the sampling procedures.
Even the researchers dedicated to this theme, attempted to establish a close relationship between cost and time spent in activities of sampling to reach sizes and most appropriate ways of SUs, ceased to evaluate the influences of biological variables and observe physical reality in each forest, such as size of the population used as a Production Unit (PU), logistic conditions in inaccessible areas and planning the activities to be carried out by the team of measurement, especially when their permanence in the field work is of limited durations. Although these latter variables are difficult to assess, they influence the final decision of this optimization process.
Pearce (1953), referring specifically to this problem, says there is a contradictory situation to reach an opti- mum size of the SUs, because the smaller units provide economy of time and larger ones provide economy of labor.
In fact, rectangular units have mostly been used in sampling the native forests. In addition to the better control of marginal trees, the detection of species diversity becomes a relevant aspect to this decision (Péllico Netto & Brena, 1997).
In forest inventories conducted in Europe, especially in Scandinavia, there is a strong preference for using small circular plots, whose decision fits into the context of Pearce (1953), i.e., given the small number of forest species occurring in the northern hemisphere and practicality in sampling with circular plots, whose radius do not exceed 3.5 meters, because, in these circumstances, you can easily use an aluminum rod to control border- line trees. In this sense, experiments carried out in Germany confirm that the circular plots up to size of 0.1 ha are more efficient to inventory native forests; when the plots, for some reason should be of larger sizes, it has been recommended the use of rectangular or square plots (Richter & Grossmann, 1951).
Messavage and Grosenbaugh (1956), working in field experiments in the north of Arkansas, in the United States, have used systematic sampling with plots arranged in rectangular structures (1:4), quadrangular (1:1) and random sampling with 16 plots of 1.6 acres in a population of 25.6 acres and held, in addition to the evaluation of accuracy of sampling, the evaluation of time spent for the measurements and for walking between the plots, aiming to compare finally the costs of these activities in the two structures of sampling.
To analyze the results of this work, the authors stated: 1) the efficiency of the sampling is the result of the joint evaluation of the cost and accuracy of the different sampling procedures used, i.e. it is a function of the re- ciprocal of the product of the square of the sampling error, given in percentage, by the cost. This is considered an optimal when the result is maximum; 2) citing Cochran (1963), they concluded that the systematic structure (1:1) resulted, for the same cost obtained, that the random sampling has an accuracy of the order of (1∕2)−0.5, and regarding the cost, for the same accuracy between both, the random sampling would cost twice as much as what it would be spend if used the systematic sampling; 3) confirmation of decrease of the coefficient of variation and sampling errors estimates with the increase of plot sizes; 4) efficiency was diverse, sometimes occurred to units of smaller sizes, sometimes to larger ones, which permitted the conclusion that is relevant introducing other va- riables sensitive to population variations, to make possible the obtainment of the optimum plot size for each for- est reality.
Considering the important results from experiments reported by Messavage and Grosenbaugh (1956), Freese (1962) suggested that the sampling efficiency could be measured as the inverse of the product of the square of the coefficient of variation by the cost
The circle shows, under the conditions mentioned above, minimization of occurrences of borderline trees, given this, among other geometric shapes of the same size, the one which has the smaller perimeter (Prodan, 1965), constituting this a first analytical contribution, for such circumstance, to decide the shape of the SU.
Smith (1938) was the precursor for the various scientific studies about size and shape of SUs, demonstrating that there is a negative exponential relationship between the plot size and the mean squared error obtained for the variable used to express it. This relationship, also known as an indicator of “Maximum Curvature”, does not reach a critical point to the choice of the plot size. Usually it is expressed using the coefficient of variation as a function of the size of the SU.
Many researchers have tried to explain scientifically this exponential reduction, because if the coefficient of variation is the ratio between the standard deviation and the mean of the estimated
It was only with the work presented by Loetsch, Zöhrer and Haller (1973) that such reduction theoretically was explained. This relationship is expressed by
Note that
Péllico Netto (1979), considering the suggestions made by Messavage and Grosenbaugh (1956), has proposed to increase the number of variables to achieve an optimum size for the SU, including: 1) Size of the area to be inventoried—Af; 2) Time for locomotion or displacement between SUs—T1; 3) Measuring time of the SUs—T2 and 4) The number of hours per day worked by the measurement team—T3. This proposal is based on the con- cept of Efficiency per Day of Work-EDT, which is based on the equation used in physics to obtain speed as a function of space and time and is strictly related to the size of the SU, expressed by:
where:
Note that the concept of EDT is equivalent to the time T3, i.e., the time spent in locomotion and in measure- ments on a day of work. If this working time, optimistically, is always close to 8 hours actually worked in a day of work, then aforciori the inventory as a whole will be optimized.
Zaide (1980) also considered time as important variable for plot size optimization and stated that the optimum achievement is reached when the time is minimum for the execution of activities of installation and measure- ment of the SUs with accuracy established for the variable considered in the sampling (usually the volume). Therefore, any experiment performed to achieve optimization of the SUs depends on the total time in function of the size of the SUs and consider such variations in a relationship that will allow minimizing the total time of measurement.
Therefore the total time for the day of work can be expressed by:
where:
As proposed by Péllico Netto (1979), it can be observed that the time of locomotion
Or:
The same way, the time of measurement may be taken as:
Given the theoretical content, the objective of this work was to integrate such knowledge experimentally to test the following hypothesis: It is possible to derive an analytical solution for a more appropriate size of the sample unit-SU and effective to be used in all types of forest inventories.
The data used to illustrate the present process of obtaining the size of SUs were collected in a native forest of Araucaria angustifolia, located in Sao Joao do Triunfo, PR, at an altitude of 780 m (25˚34'18" Latitude and 50˚05'56" Longitude W), belonging to the Federal University of Parana, with 32 ha (Tello, 1980) and will be mentioned later as araucaria experiment—AE.
The demarcation of the experimental area was performed considering the topography and edaphic conditions the most uniform possible, to allow comparisons between shapes, sizes and time spent for measurements of SUs, so that they could be performed with maximum similarity.
The experimental area was delimited in nine fields, each with one hectare, in which were laid down the varia- tions of sample structure more commonly used in forest inventories, i.e.: circle, square and rectangular SUs, with variations from 200 m2 up to 1000 m2, as is shown in
After demarcations of the SUs, DBH, total and commercial heights of all trees were collected, being that for Araucária (DAP
The equations used for obtaining the volume per species are presented in
The sample which was applied to each size of SUs came from a pilot inventory to have the preliminary esti- mators of mean and variance of various shapes and sizes of SUs, to what it has been established an error limit of 10% of the average volume with 95% probability. The results for a sampling area of 9 ha are presented in
The time spent for each one of the activities carried out in the SUs were obtained using a stopwatch and computed individually, as suggested by Stöhr (1978), to prevent the operations of subtraction performed to ob- tain the liquid times by activities, when using continuous timing. The result of this operation is shown in
Tello (1980) using such obtained estimators and having them applied in Equation (1) proposed by Freese (1962) concluded that the maximum efficiency in all forms occurred with the size of 1000 m2 with a small ad- vantage for the circular shape, which minimized the total time for this size. This conclusion clearly shows that
. Shapes and sizes of SUs used in the experiment.
Sizes | Circles (Rays) (m) | Square Shapes (Sides) (m) | Rectangular Shapes (Sides) (m) |
---|---|---|---|
200 m2 | 7.98 | 14.14 × 14.14 | 10.00 × 20.00 |
400 m2 | 11.28 | 20.00 × 20.00 | 14.14 × 28.28 |
600 m2 | 13.82 | 24.49 × 24.49 | 17.32 × 34.64 |
800 m2 | 15.96 | 28.28 × 28.28 | 20.00 × 40.00 |
1000 m2 | 17.84 | 31.62 × 31.62 | 22.36 × 44.72 |
. Equations adjusted for calculations of tree volumes in the SUs.
Species | Equations |
---|---|
Araucaria (Total Volume) | . Equations adjusted for calculations of tree volumes in the SUs. |
Araucaria (Commercial Volume) | . Equations adjusted for calculations of tree volumes in the SUs. |
Other species (Commercial Volume) | . Equations adjusted for calculations of tree volumes in the SUs. |
. Sample sizes demanded for sampling and average distance (D) between the SUs.
Size (m2) | Circle Shape | Square Shape | Rectangular Shape | |||
---|---|---|---|---|---|---|
n | D (m) | n | D (m) | n | D (m) | |
200 | 135 | 25.82 | 120 | 27.38 | 166 | 23.28 |
400 | 69 | 36.12 | 62 | 38.10 | 62 | 38.10 |
600 | 43 | 45.75 | 46 | 44.23 | 42 | 46.29 |
800 | 29 | 55.71 | 32 | 53.03 | 30 | 54.77 |
1000 | 24 | 61.24 | 23 | 62.55 | 21 | 65.47 |
Note: n = sample size, D = average distance between SUs.
. Average speeds for locomotion, measurements and total in different sampling schemes.
Form | Size (m2) | Locomotion (Minutes) | Measurement (Minutes) | Total (Minutes) | Average Speeds | nd | |
---|---|---|---|---|---|---|---|
v1 (m∙h−1) | v2 (m2∙h−1) | ||||||
Circle | 200 | 08.65 | 09.92 | 18.57 | 179.10 | 1209.67 | 25 |
400 | 10.01 | 20.54 | 30.55 | 216.50 | 1168.45 | 15 | |
600 | 11.13 | 33.54 | 44.67 | 246.63 | 1073.35 | 10 | |
800 | 13.70 | 45.55 | 59.25 | 243.99 | 1053.79 | 08 | |
1000 | 14.20 | 53.45 | 67.65 | 258.76 | 1122.54 | 07 | |
Average | 1125.56 | ||||||
Square | 200 | 08.99 | 16.50 | 25.49 | 182.74 | 727.27 | 18 |
400 | 10.85 | 28.21 | 39.06 | 210.69 | 850.76 | 12 | |
600 | 11.03 | 43.95 | 54.98 | 240.60 | 819.11 | 08 | |
800 | 13.10 | 58.19 | 71.29 | 242.88 | 824.88 | 06 | |
1000 | 14.30 | 67.55 | 81.85 | 262.45 | 888.23 | 05 | |
Average | 822.05 | ||||||
Rectangular | 200 | 08.32 | 15.50 | 23.82 | 167.88 | 774.19 | 20 |
400 | 10.85 | 27.61 | 38.46 | 210.69 | 869.25 | 12 | |
600 | 11.16 | 42.95 | 54.11 | 248.87 | 838.18 | 08 | |
800 | 13.25 | 53.59 | 66.84 | 248.02 | 895.69 | 07 | |
1000 | 14.45 | 68.85 | 83.30 | 271.85 | 871.46 | 05 | |
Average | 849.75 |
the experiment could not cover broad conditions in integration of the functions, hence the need to use simulation to achieve more appropriate conditions in the evaluation of SUs of sizes larger than 1000 m2. Still, the incorpo- ration of other variables in this process, as has been suggested by Messavage and Grosenbaugh (1956) was cri- tical to ensure greater consistency in the final decision of the size and shape of the SU considering those expe- rimental conditions.
The simulation was structured using the mathematical models: power for the times and a hyperbole for the coefficient of variation, as is shown in
The estimate of the total time was calculated by the sum of the estimators of the times of locomotion and measurement, as presented in Equation (7):
where:
The estimate of relative efficiency was calculated using the Equation (1). If
The derivative of
where:
Therefore:
To obtain the maximum point of the function the derivative was equaled to zero.
To simplify the equation and create an independent term, all terms were divided by
The value of the area
The equations adjusted and the statistics of the regression for the times of locomotion, measurements and coefficient of variation, obtained with the adjustment of the models shown in
estimate of the total time was calculated by the sum of these estimators using Equation (7), and is presented in
. Mathematical models used.
Variable | Model used |
---|---|
Time of locomotion | . Mathematical models used. |
Time of taking measurement | . Mathematical models used. |
Coefficient of Variation | . Mathematical models used. |
Note: where:
. Adjusted Equations and statistics for T1, T2 and CV.
Variable | Adjusted Equation | F | . Adjusted Equations and statistics for T1, T2 and CV. | Syx% |
---|---|---|---|---|
Locomotion Time | . Adjusted Equations and statistics for T1, T2 and CV. | 209.92** | 93.72% | 1.92% |
Measurement Time | . Adjusted Equations and statistics for T1, T2 and CV. | 2.14** | 92.43% | 4.60% |
Coefficient of Variation | . Adjusted Equations and statistics for T1, T2 and CV. | 5.92** | 87.15% | 8.24% |
Note: where:
Estimate of the times for locomotion, measurements and total
lation to the size of sample units. The dots represent the values observed for the sizes of plots of 200, 400, 600,
800 and 1000 m2 for shapes: circle, square and rectangular. The curve represents the values estimated for the sizes up to 10,000 m2.
The relative efficiency
The graph in
Substituting in Equation (13) the values of the coefficients obtained by the adjustment of the models listed in
It was used to solve the function Matlab software, as is shown in
The form of SUs was analyzed considering the minimum total time spent to locate and measure them. As they are presented in
For sizes larger than 0.1 ha is advisable to use the rectangular units, because its efficiency increases due to
Coefficient of variation
Relative efficiency
Analytical solution for the optimum size of WATER obtained by Matlab software
better control of edge effects, as well as favoring the detection of species diversity in native forests. In the case of forest plantations, due to the regularity of the lines of plantations and the spacing between trees, the research- ers prefer to use the rectangular shape of SU, for ease of installation and demarcation of the SU, mainly when applying successive sampling in continuous forest inventories (Cesario, Engel, Finger, & Schneider, 1994; Simplício, Muniz, Aquino, & Soares, 1996; Strand, 1957).
The maximum amount of sample units per day of work
The graph in
Sampling units per day of work
measured in 8 hours (one day of work), for this data set, is 7000 m2.
The simulations using the integration of the assessed variables permitted the optimization, by means of an analytical process to obtain the most appropriated plot size for the sampling conditions used in the forest popu- lation, in this case 600 m2, which proves the formulated hypothesis. Such procedure was possible, however, due to some relevant aspects detected in structuring this process: 1) it was not possible to derive Equation (8) using the power function for expressing the relationship of CV with the plot sizes, although it has been used in many scientific works since its proposition by Lessman and Atkins (1963). It was due to the proposed use of a hyper- bolic function by the authors which became possible the functionality of the analytical process until obtaining the optimum plot size; 2) Although this analytic solution has already been formalized by Zaide (1980), using the power functions for the three variables of the process (times of locomotion and measurements, and coefficient of variation), he succeeded because he determined the power values of such functions. When we obtained them by adjusting a linear regression model, the derivative process did not work.
Applying the experimental results in Equation (2) it was possible to get the value of
Cesaro et al. (1994) compared fixed area plots with the Bitterlich and Prodan SUs in Rio Grande do Sul, hav- ing assessed carefully the times of location, installation and measurements of dbh and height of the trees in an area of 4.6 ha for the three sampling methods. They used a fixed area plot of 600 m2, which allowed compari- sons with the results obtained in this work, in spite of their experiment have been conducted in stands of Pinus sp., with 28 years of age, which will be referenced later as pine experiment—PE.
In the PE, after having performed the calculation of the estimators needed for comparative purposes, the av- erage distance between the units of fixed area sampled resulted in
In several forest inventories of companies carried out in pine plantations have been reported that such effi- ciency is possible, however it is worth noting that the conditions of access and the level of experience of the field teams should be very good. In most adverse conditions of access the average velocities, both for locomo- tion or for measurement, tend to be smaller. Still, as discussed in the AE, such number of SUs per day of work was obtained without the incorporation of the variable coefficient of variation.
This type of experiment can be repeated both in natural forest, as well as in plantations to define the appropri- ate plot size to be used in the forest inventory, however, from a practical point of view, it can be assumed, at the light of the result obtained in the AE, the constraints that will allow approaching to the optimum plot size, mainly in the case of plantations.
Faced with the premise that such evaluation can be performed with the same effectiveness by the sampling team in one day of work (EDT), it is not relevant anymore to focus on planning of the forest inventory as a whole, and what is important to be known is summarized in
The area of the company to be measured is 35 ha, with four compartments, which has two or three plots in each one of them. It will be assumed that the company has performed a timing to see about the times spent for locomotion between the plots and also to perform the measurements of all dbh and heights of a percentage of trees. With these data at his disposal it is possible to calculate the speeds for locomotion and measurement and can also clearly specify in its planning that the team must work 8 hours per day performing the sampling activi- ties.
Using the result obtained previously
Replacing it in the data of
As noted, when it was expanded the number of variables in the decision-making process to decide what is the optimum size of SU, the measurement of 18 SUs per day allows you to measure units of only 471.14 m2 which can be rounded up to 470 m2. However, if the company believes that the unit of 600 m2 should be maintained as the best option, they may reduce the number of units to be measured in a day of work
. Preliminary data collected by the company for definition of the optimum size of the SU.
Fields | No. of plots | Areas (ha) | Average distance between plots (m) | Time Locom. (minutes) | Time Measur. (minutes) | Speed. v1 (m∙h-1) | Speed. v2 (m2∙h-1) |
---|---|---|---|---|---|---|---|
T1 | 2 | 05 | 158.11 | 11.83 | ̶ | 801.91 | ̶ |
T2 | 2 | 08 | 200.00 | 14.97 | ̶ | 801.60 | ̶ |
T3 | 3 | 12 | 200.00 | 14.97 | ̶ | 801.60 | ̶ |
T4 | 2 | 10 | 223.61 | 16.73 | ̶ | 801.95 | ̶ |
Total | 9 | 35 | 197.20 | 14.75 | 19.94 | 802.18 | 1805.42 |
The hyperbolic function proposed by the authors to express the coefficient of variation as a function of plot size became the possible functionality of the analytical process to obtain the optimum plot size;
The simulation using the integration of the assessed variables permitted to get the optimization, and the most appropriated plot size for the sampling conditions of that forest population was 600 m2;
The use of additional variables as size of the area to be inventoried
Applying the combination of the two procedures: optimization process and EDT evaluation, it was possible to obtain the most appropriate number of plots of 600 m2 to be measured in one day of field work, which in the case were 15 plots.