Towards the harmonization between National Forest Inventory and Forest Condition Monitoring. Consistency of plot allocation and effect of tree selection methods on sample statistics in Italy

In the frame of a process aiming at harmonizing National Forest Inventory (NFI) and ICP Forests Level I Forest Condition Monitoring (FCM) in Italy, we investigated (a) the long-term consistency between FCM sample points (a subsample of the first NFI, 1985, NFI_1) and recent forest area estimates (after the second NFI, 2005, NFI_2) and (b) the effect of tree selection method (tree-based or plot-based) on sample composition and defoliation statistics. The two investigations were carried out on 261 and 252 FCM sites, respectively. Results show that some individual forest categories (larch and stone pine, Norway spruce, other coniferous, beech, temperate oaks and cork oak forests) are over-represented and others (hornbeam and hophornbeam, other deciduous broadleaved and holm oak forests) are under-represented in the FCM sample. This is probably due to a change in forest cover, which has increased by 1,559,200 ha from 1985 to 2005. In case of shift from a tree-based to a plot-based selection method, 3,130 (46.7 %) of the original 6,703 sample trees will be abandoned, and 1,473 new trees will be selected. The balance between exclusion of former sample trees and inclusion of new ones will be particularly unfavourable for conifers (with only 16.4 % of excluded trees replaced by new ones) and less for deciduous broadleaves (with 63.5 % of excluded trees replaced). The total number of tree species surveyed will not be impacted, while the number of trees per species will, and the resulting (plot-based) sample composition will have a much larger frequency of deciduous broadleaved trees. The newly selected trees have—in general—smaller diameter at breast height (DBH) and defoliation scores. Given the larger rate of turnover, the deciduous broadleaved part of the sample will be more impacted. Our results suggest that both a revision of FCM network to account for forest area change and a plot-based approach to permit statistical inference and avoid bias in the tree sample composition in terms of DBH (and likely age and structure) are desirable in Italy. As the adoption of a plot-based approach will keep a large share of the trees formerly selected, direct tree-by-tree comparison will remain possible, thus limiting the impact on the time series comparability. In addition, the plot-based design will favour the integration with NFI_2.     

Three-Dimensional Modeling of an Urban Park and Trees by Combined Airborne and Portable On-Ground Scanning LIDAR Remote Sensing

In this study, we confirmed the utility of airborne and portable on-ground scanning light detection and ranging (LIDARs) for three-dimensional visualization of an urban park and quantification of biophysical variables of trees in the park. The digital canopy height model (DCHM) and digital terrain model generated from airborne scanning LIDAR data provided precise images of the ground surface and individual tree canopies. The heights of 166 coniferous and broadleaf trees of 11 species in the park were estimated from the DCHM images with slight underestimation (mean error = −0.14 m, RMSE = 0.30 m). Portable on-ground scanning LIDAR provided images of individual trees with detailed features. Tree height and trunk diameter were estimated to be within 0.31 m and 1 cm, respectively, from the on-ground LIDAR images. We combined airborne and on-ground LIDAR images to overcome blind regions and created a complete three-dimensional model of three standing trees. The model allowed not only visual assessment from all viewpoints but also quantitative estimation of canopy volume, trunk volume, and canopy cross-sectional area.     

Measurements of stem diameter: implications for individual- and stand-level errors

Stem diameter is one of the most common measurements made to assess the growth of woody vegetation, and the commercial and environmental benefits that it provides (e.g. wood or biomass products, carbon sequestration, landscape remediation). Yet inconsistency in its measurement is a continuing source of error in estimates of stand-scale measures such as basal area, biomass, and volume. Here we assessed errors in stem diameter measurement through repeated measurements of individual trees and shrubs of varying size and form (i.e. single- and multi-stemmed) across a range of contrasting stands, from complex mixed-species plantings to commercial single-species plantations. We compared a standard diameter tape with a Stepped Diameter Gauge (SDG) for time efficiency and measurement error. Measurement errors in diameter were slightly (but significantly) influenced by size and form of the tree or shrub, and stem height at which the measurement was made. Compared to standard tape measurement, the mean systematic error with SDG measurement was only ?0.17 cm, but varied between ?0.10 and ?0.52 cm. Similarly, random error was relatively large, with standard deviations (and percentage coefficients of variation) averaging only 0.36 cm (and 3.8%), but varying between 0.14 and 0.61 cm (and 1.9 and 7.1%). However, at the stand scale, sampling errors (i.e. how well individual trees or shrubs selected for measurement of diameter represented the true stand population in terms of the average and distribution of diameter) generally had at least a tenfold greater influence on random errors in basal area estimates than errors in diameter measurements. This supports the use of diameter measurement tools that have high efficiency, such as the SDG. Use of the SDG almost halved the time required for measurements compared to the diameter tape. Based on these findings, recommendations include the following: (i) use of a tape to maximise accuracy when developing allometric models, or when monitoring relatively small changes in permanent sample plots (e.g. National Forest Inventories), noting that care is required in irregular-shaped, large-single-stemmed individuals, and (ii) use of a SDG to maximise efficiency when using inventory methods to assess basal area, and hence biomass or wood volume, at the stand scale (i.e. in studies of impacts of management or site quality) where there are budgetary constraints, noting the importance of sufficient sample sizes to ensure that the population sampled represents the true population.     

Crown density assessments, control surveys and reproducibility

This study aims to evaluate the quality of crown density data, based on independent, pairwise tree assessments. The data originates from monitoring of forest health (crown condition) in Norway; 250 plots, comprising 12 000 individual trees of Norway spruce, have been reassessed by a single observer during 1990–95. Of the trees, 2300 were controlled more than twice, providing the possibility of evaluating the quality of assessed temporal changes of crown density. True errors (standard deviation) are estimated to be about 10% for single trees and 5% for plot means, while the real standard deviation of the differences were slightly higher. The errors of the temporal changes of crown density were of similar magnitude. Systematic differences in crown density were found between sites and plot types, partly resulting from observer bias. However, the results suggest that observer bias is really the result of each observer's personal style in assessment.     

Relation between individual tree mortality and tree characteristics in a polluted and non-polluted environment

Data on individual tree mortality in relatively healthy (Berezinskiy biosphere reserve) and damaged (surroundings of the mineral fertilizer plant 'Achema') even-aged Scots pine (Pinus sylvestris L.) stands are presented. Tree size and competition intensity were found to be the most significant predictors of individual tree mortality in all age groups of the relatively healthy Scots pine stands, however, an essential reduction in the closeness of relations between the tree mortality rate and these variables was determined with the aging of stands. An exponential decrease in tree mortality probability with an increase of tree size is characteristic for trees suffering different competition intensity, however, this decrease becomes much more pronounced as the competition pressure increases. The relations of different tree and stand variables with tree mortality probability have been found to become much weaker in the polluted environment. An exponential increase in tree mortality probability with an increase of crown defoliation was characteristic of damaged Scots pine stands, however, the rate of the increase was different in different age and diameter classes. The impact of crown defoliation on tree mortality rate increased with the aging of stands. At the same defoliation level, individual tree mortality probability was much higher for smaller suppressed trees, however, a relative increase in tree mortality probability along with an increase of crown defoliation was more pronounced for dominant trees. Conclusion: a higher mortality of damaged (defoliated) trees should be considered while assessing losses in forest productivity in a polluted environment.     

Relation Between Individual Tree Mortality and Tree Characteristics in a Polluted and Non-Polluted Environmentx

Data on individual tree mortality in relatively healthy (Berezinskiy biosphere reserve) and damaged (surroundings of the mineral fertilizer plant ‘Achema’) even-aged Scots pine (Pinus sylvestris L.) stands are presented. Tree size and competition intensity were found to be the most significant predictors of individual tree mortality in all age groups of the relatively healthy Scots pine stands, however, an essential reduction in the closeness of relations between the tree mortality rate and these variables was determined with the aging of stands. An exponential decrease in tree mortality probability with an increase of tree size is characteristic for trees suffering different competition intensity, however, this decrease becomes much more pronounced as the competition pressure increases. The relations of different tree and stand variables with tree mortality probability have been found to become much weaker in the polluted environment. An exponential increase in tree mortality probability with an increase of crown defoliation was characteristic of damaged Scots pine stands, however, the rate of the increase was different in different age and diameter classes. The impact of crown defoliation on tree mortality rate increased with the aging of stands. At the same defoliation level, individual tree mortality probability was much higher for smaller suppressed trees, however, a relative increase in tree mortality probability along with an increase of crown defoliation was more pronounced for dominant trees. Conclusion: a higher mortality of damaged (defoliated) trees should be considered while assessing losses in forest productivity in a polluted environment.     

Can tree species diversity be assessed with Landsat data in a temperate forest?

The diversity of forest trees as an indicator of ecosystem health can be assessed using the spectral characteristics of plant communities through remote sensing data. The objectives of this study were to investigate alpha and beta tree diversity using Landsat data for six dates in the Gönen dam watershed of Turkey. We used richness and the Shannon and Simpson diversity indices to calculate tree alpha diversity. We also represented the relationship between beta diversity and remotely sensed data using species composition similarity and spectral distance similarity of sampling plots via quantile regression. A total of 99 sampling units, each 20 m × 20 m, were selected using geographically stratified random sampling method. Within each plot, the tree species were identified, and all of the trees with a diameter at breast height (dbh) larger than 7 cm were measured. Presence/absence and abundance data (tree species number and tree species basal area) of tree species were used to determine the relationship between richness and the Shannon and Simpson diversity indices, which were computed with ground field data, and spectral variables derived (2 × 2 pixels and 3 × 3 pixels) from Landsat 8 OLI data. The Shannon-Weiner index had the highest correlation. For all six dates, NDVI (normalized difference vegetation index) was the spectral variable most strongly correlated with the Shannon index and the tree diversity variables. The Ratio of green to red (VI) was the spectral variable least correlated with the tree diversity variables and the Shannon basal area. In both beta diversity curves, the slope of the OLS regression was low, while in the upper quantile, it was approximately twice the lower quantiles. The Jaccard index is closed to one with little difference in both two beta diversity approaches. This result is due to increasing the similarity between the sampling plots when they are located close to each other. The intercept differences between two investigated beta diversity were strongly related to the development stage of a number of sampling plots in the tree species basal area method. To obtain beta diversity, the tree basal area method indicates better result than the tree species number method at representing similarity of regions which are located close together. In conclusion, NDVI is helpful for estimating the alpha diversity of trees over large areas when the vegetation is at the maximum growing season. Beta diversity could be obtained with the spectral heterogeneity of Landsat data. Future tree diversity studies using remote sensing data should select data sets when vegetation is at the maximum growing season. Also, forest tree diversity investigations can be identified by using higher-resolution remote sensing data such as ESA Sentinel 2 data which is freely available since June 2015.     

Spatial pattern analysis in Persian oak (Quercus brantii var. persica) forests on B&W aerial photographs

The purpose of this investigation was to develop a method to determine the spatial pattern of trees as a robust indicator to monitor changes from B&W aerial photographs in Persian oak forests of Zagros, Iran. A 500 x 600 m study area was selected in Servak forests next to Yasuj city in Kohgiluyeh-Va-BuyerAhmad Province. All the trees were tagged in the study area and the point map of stems were prepared. The spatial distribution of trees was determined as "dispersed" using nearest neighbour technique. Then the index of "C" calculated by T-square sampling method was applied to the point map of the study area in 30 systematic sample points in a 100 x 100 m network. Comparing the results of this method with the true spatial pattern of the study area showed that "C" can detect the spatial arrangement of trees. Thereafter the index was used on the air photo of the study area that was made of B&W aerial photographs. The method suggested in this study provides a suitable approach for detecting the spatial pattern of trees in Zagros forests on B&W air photos.     

Evaluation of a simultaneous sampling method of PAHs, PCDD/Fs and dl-PCBs in ambient air

The Stockholm Convention on Persistent Organic Pollutants (POPs) was signed in May 2001 by 127 countries. Currently, 12 substances are regulated by the convention, and the work on finding new candidate chemicals to the convention has started. Among these 12 substances, dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs) and polychlorinated biphenyls (PCBs) are the objects of this study. There are no European standard methods for the simultaneous determination of these microorganic compounds, together with polycyclic aromatic hydrocarbons (PAHs) in ambient air--they must be referred to three different methods developed by the U.S. Environmental Protection Agency (EPA). The quali-quantitative analysis of these microorganic pollutants is an important challenge due to the low concentrations at which they may be present. In this study, the development of a simplified, alternative, fast and affordable sampling method for the determination of PAHs, PCDDs, PCDFs and PCBs in ambient air was performed. The sampling time was extended from 24 h to 7 d in order to enrich the sample, to fall within the instrumental limits of detection and to reduce the number of samples to be processed and, therefore, errors that may arise. First of all, experiments with labelled standards were conducted in the research area of Montelibretti (rural station, which is sited about 20 km northeast of Rome), with the purpose of optimizing sampling efficiency. Finally, the method was applied to the analysis of these compounds in the air of a suburban area with small industrial plants in order to evaluate the feasibility of the proposed sampling method system, by comparing concentrations of native compounds acquired during simultaneous daily and weekly sampling.     

Accuracy assessment of a small-area method for estimating the spatial distribution of the degree of tree damage

Tree damage, gauged by the amount of defoliation, is one of the basic criteria used to determine treatments for protected and economic forests. Monitoring should include an assessment of the degree of tree damage in different spatial scales. Therefore, in addition to the commonly applied large-area methods, small-area methods should be used. The aim of the paper is to present the results of the accuracy assessment of a small-area method, proposed by Podlaski (2005) [Podlaski, R. (2005). Inventory of the degree of tree defoliation in small areas. Forest Ecology and Management, 215, 361–377], for monitoring the degree of tree damage. The degree of tree damage was shown in sub-blocks P₃ of the system of information on natural environment (SINUS). To estimate the spatial distribution of the degree of tree defoliation, survey sampling, based on simple random sampling with replacement (SRSWR), was used. The degree of damage to fir (Abies alba Mill.) and beech (Fagus sylvatica L.) was analysed in the Święty Krzyż forest section in the Świętokrzyski National Park. The maximum total estimation errors for the proportion of trees with a degree zero of damage, and with second and third degrees of damage together (for α = 0.05) were at most 30.8% for fir and 24.3% for beech trees. For standard, small-area evaluations, these are satisfactory values. In the Święty Krzyż forest section, the number of P₃ sub-blocks with 0.00–5.00% of undamaged trees and with 80.01–100.00% of moderately- or severely-damaged trees was significantly greater for fir than for beech. These results indicate that the fir population was unhealthier than the beech group in the study area. P₃ sub-blocks of the SINUS system, in which the proportion of the healthiest trees was highest, were situated at the forest margin, bordering on meadows and arable fields (in the case of fir) and forming dense patches consisting of several sub-blocks, or occurring singly in the whole study area (in the case of beech). The results show the significant differentiation of forest tree health in small areas.     

Mapping trees outside forests using high-resolution aerial imagery: a comparison of pixel- and object-based classification approaches

Discrete trees and small groups of trees in nonforest settings are considered an essential resource around the world and are collectively referred to as trees outside forests (ToF). ToF provide important functions across the landscape, such as protecting soil and water resources, providing wildlife habitat, and improving farmstead energy efficiency and aesthetics. Despite the significance of ToF, forest and other natural resource inventory programs and geospatial land cover datasets that are available at a national scale do not include comprehensive information regarding ToF in the United States. Additional ground-based data collection and acquisition of specialized imagery to inventory these resources are expensive alternatives. As a potential solution, we identified two remote sensing-based approaches that use free high-resolution aerial imagery from the National Agriculture Imagery Program (NAIP) to map all tree cover in an agriculturally dominant landscape. We compared the results obtained using an unsupervised per-pixel classifier (independent component analysis—[ICA]) and an object-based image analysis (OBIA) procedure in Steele County, Minnesota, USA. Three types of accuracy assessments were used to evaluate how each method performed in terms of: (1) producing a county-level estimate of total tree-covered area, (2) correctly locating tree cover on the ground, and (3) how tree cover patch metrics computed from the classified outputs compared to those delineated by a human photo interpreter. Both approaches were found to be viable for mapping tree cover over a broad spatial extent and could serve to supplement ground-based inventory data. The ICA approach produced an estimate of total tree cover more similar to the photo-interpreted result, but the output from the OBIA method was more realistic in terms of describing the actual observed spatial pattern of tree cover.     

Inventory methods for trees in nonforest areas in the great plains states

The US Forest Service’s Forest Inventory and Analysis (FIA) program collects information on trees in areas that meet its definition of forest. However, the inventory excludes trees in areas that do not meet this definition, such as those found in urban areas, in isolated patches, in areas with sparse or predominantly herbaceous vegetation, in narrow strips (e.g., shelterbelts), or in riparian areas. In the Great Plains States, little is known about the tree resource in these noninventoried, nonforest areas, and there is a great deal of concern about the potential impact of invasive pests, such as the emerald ash borer. To address this knowledge gap, FIA’s National Inventory and Monitoring Applications Center has partnered with state cooperators and others in a project called the Great Plains Initiative to design and implement an inventory of trees in nonforest areas. The goal of the inventory is to characterize the nonforest tree resource using methods compatible with those of FIA so a holistic understanding of the resource can be obtained by integrating the two surveys. The goal of this paper is to describe the process of designing and implementing the survey, including plot and sample design, and to present some example results from a reporting tool we developed.     

Assessment of vehicular pollution in Kolkata, India, using CALINE 4 model

Michigan water quality standards for public bathing beaches require local health departments to collect and analyze a minimum of three water samples for Escherichia coli during each sampling event. The geometric mean number of E. coli colonies is then compared to the 300 colonies per 100 ml standard to determine compliance. This article compares the results of the currently mandated procedure to a composite sampling method, whereby the three samples are mixed in equal volumes and analyzed once. This effectively replaces the geometric mean of the individual sample results with an arithmetic mean. Although arithmetic means are more affected by outliers, this sensitivity to high concentrations is more health conservative than the geometric mean. During the 2007 sampling season, nine bathing beaches were monitored once each week. Three individual point samples and a composite sample were analyzed for each sampling event. No statistically significant differences in bacteria concentrations were found between composite sample analysis and the arithmetic mean of individual point sample analyses. No violations were detected in the 2007 sampling season, so using historical data, a retrospective analysis was performed on samples gathered at nine bathing beaches in Kalamazoo County, Michigan during the years 2001–2007. The arithmetic mean of the three samples taken at each site served as a surrogate composite sample. The benefits of compositing the three samples were investigated assuming a 2/3 reduction in analytical costs. In the traditional sampling method, three individual samples were obtained and analyzed once in every 3-week period during the summer season, whereas compositing was simulated by taking the arithmetic mean of each week’s results. The results of this retrospective cost analysis indicates that ten to 14 violations would have been missed using the less frequent traditional sampling and analysis methodology. Composite sampling is a cost-saving alternative to traditional sampling techniques that can be more protective of public health, particularly when the savings are applied to increased numbers of samples in time or space.     

Effects of snow-reflected light levels on human visual comfort

The intensity of the sunlight reflected by the snow-covered surfaces is so high that it may disturb humans many times. This study aims to determine the reflected sunlight intensities from snow covered areas at points near (at a distance of 2 m) and under an individual tree and among trees (in the forest area) by accepting the open area as control; the reducing effects of the plant materials on reflected sunlight in percentage by comparing with the values of the open (control) area; and critical reflected sunlight threshold values for human visual comfort. The study was carried out over 22 clear and calm, i.e. sky was cloudless and wind was calm, days between the 1st and 31st days of January 2004, at 8:30 in the morning, at 12:30 at noon and at 14:30 in the afternoon in Erzurum. In order to determine the discomforting light intensity levels, 25 females and 26 male (totally 51) student subjects whose mean age was 20 and who had no visual disorders were selected. Considering the open area as control, mean reflected sunlight reducing effects were found to be 19.0, 66.0 and 82.7% for the 2 m near a tree, under a tree, and forest area, respectively. According to the responses of 51 subjects in the study, visually "very comfortable" range is between 5,000 and 8,000 lx; "comfortable" range is between 11,000 and 75,000 lx (mostly at 12,000 lx); "uncomfortable" condition is above the light intensity value of 43,000 lx and "very uncomfortable" condition is above the intensity of 80,000 lx. Great majority of the subjects (91%) found the value of 103,000 lx to be "very uncomfortable." As it is not an applicable way to use the great and dense tree masses in the cities, at least individual trees should be used along the main pedestrian axels in the cities having the same features with Erzurum to prevent the natural light pollution and discomforting effects of the snow-reflected sunlight.     

Support,shape and number of replicate samples for tree foliage analysis

Many fundamental features of a sampling program are determined by the heterogeneity of the object under study and the settings for the error (alpha), the power (beta), the effect size (ES), the number of replicate samples, and sample support, which is a feature that is often overlooked. The number of replicates, alpha, beta, ES, and sample support are interconnected. The effect of the sample support and its shape on the required number of replicate samples was investigated by means of a resampling method. The method was applied to a simulated distribution of Cd in the crown of a Salix fragilis L. tree. Increasing the dimensions of the sample support results in a decrease in the variance of the element concentration under study. Analysis of the variance is often the foundation of statistical tests, therefore, valid statistical testing requires the use of a fixed sample support during the experiment. This requirement might be difficult to meet in time-series analyses and long-term monitoring programs. Sample supports have their largest dimension in the direction with the largest heterogeneity, i.e. the direction representing the crown height, and this will give more accurate results than supports with other shapes. Taking the relationships between the sample support and the variance of the element concentrations in tree crowns into account provides guidelines for sampling efficiency in terms of precision and costs. In terms of time, the optimal support to test whether the average Cd concentration of the crown exceeds a threshold value is 0.405 m3 (alpha = 0.05, beta = 0.20, ES = 1.0 mg kg(-1) dry mass). The average weight of this support is 23 g dry mass, and 11 replicate samples need to be taken. It should be noted that in this case the optimal support applies to Cd under conditions similar to those of the simulation, but not necessarily all the examinations for this tree species, element, and hypothesis test.     

超低沸点挥发性有机物测定中的要点分析

采用大气浓缩仪-气相色谱法测定超低沸点挥发性有机物乙烷、乙烯、丙烷、丙烯、乙炔,选用PLOT-Q柱分析,控制真空罐内相对湿度约50%,通过优化测定条件使方法线性良好,5种低沸点化合物的检出限为0.08 nmol/mol～0.2 nmol/mol,样品保存时间不超过10 d。用该方法测定某地区环境空气中VOCs,结果检出率均为100%。     

Determining the agricultural ammonia immission using bark bio-monitoring: comparison with passive sampler measurements

Bark samples of spruce, pine and oak trees were collected at two sites in southern Bavaria which are characterized by high agricultural ammonia emissions. The samples were taken using a recently developed bark sampling device which removes a defined layer of the bark. The bark was then analysed for ammonium concentration in order to reflect the environmental ammonia immission. The measured bark concentrations decreased with rising distance between the sample trees and the ammonia source. This applied (i) to measurements inside a closed forest stand ranging from forest edge with high immission to forest interior with much lower immission, and (ii) to the open field where single-standing trees were sampled. Comparing the ammonium concentrations among the three different tree species revealed significant correlations. Thus, it could be shown that old spruce trees are as usable for bark bio-monitoring as the traditionally used pine and oak trees. The ammonium concentrations of the bark were significantly correlated to measurements taken by ammonia passive samplers at the same locations. These results indicate that bark samples may be used for a standardised monitoring of airborne ammonia load. A major advantage of the technique is the determination of the long-term accumulative ammonia load using a single measurement.     

Diameter distribution of some subalpine fir stands in central British Columbia

Diameter distributions of subalpine fir (Abies lasiocarpa (Hook.) Nutt.) in central British Columbia were investigated. Nine fire-originated, old-growth stands were selected with maximum tree age at breast height of about 300 years. The stage of stand development was determined by testing how well a negative exponential function matched the cumulative diameter distribution. As a comparison, the negative exponential function was also fitted on the frequency distributions. Eight out of the nine distributions showed fewer medium-sized trees and more large-sized trees than predicted by the negative exponential function. One stand showed relatively good fit, but failed the lack-of-fit test. Although the diameter distributions are close to a balanced stage, stand establishment patterns are still evident 300 years after disturbance.