Sampling design for long-term regional trends in marine rocky intertidal communities

Probability-based designs reduce bias and allow inference of results to the pool of sites from which they were chosen. We developed and tested probability-based designs for monitoring marine rocky intertidal assemblages at Glacier Bay National Park and Preserve (GLBA), Alaska. A multilevel design was used that varied in scale and inference. The levels included aerial surveys, extensive sampling of 25 sites, and more intensive sampling of 6 sites. Aerial surveys of a subset of intertidal habitat indicated that the original target habitat of bedrock-dominated sites with slope ≤30° was rare. This unexpected finding illustrated one value of probability-based surveys and led to a shift in the target habitat type to include steeper, more mixed rocky habitat. Subsequently, we evaluated the statistical power of different sampling methods and sampling strategies to detect changes in the abundances of the predominant sessile intertidal taxa: barnacles Balanomorpha, the mussel Mytilus trossulus, and the rockweed Fucus distichus subsp. evanescens. There was greatest power to detect trends in Mytilus and lesser power for barnacles and Fucus. Because of its greater power, the extensive, coarse-grained sampling scheme was adopted in subsequent years over the intensive, fine-grained scheme. The sampling attributes that had the largest effects on power included sampling of “vertical” line transects (vs. horizontal line transects or quadrats) and increasing the number of sites. We also evaluated the power of several management-set parameters. Given equal sampling effort, sampling more sites fewer times had greater power. The information gained through intertidal monitoring is likely to be useful in assessing changes due to climate, including ocean acidification; invasive species; trampling effects; and oil spills.     

Design to monitor trend in abundance and presence of American beaver (Castor canadensis) at the national forest scale

Wildlife conservationists design monitoring programs to assess population dynamics, project future population states, and evaluate the impacts of management actions on populations. Because agency mandates and conservation laws call for monitoring data to elicit management responses, it is imperative to design programs that match the administrative scale for which management decisions are made. We describe a program to monitor population trends in American beaver (Castor canadensis) on the US Department of Agriculture, Black Hills National Forest (BHNF) in southwestern South Dakota and northeastern Wyoming, USA. Beaver have been designated as a management indicator species on the BHNF because of their association with riparian and aquatic habitats and its status as a keystone species. We designed our program to monitor the density of beaver food caches (abundance) within sampling units with beaver and the proportion of sampling units with beavers present at the scale of a national forest. We designated watersheds as sampling units in a stratified random sampling design that we developed based on habitat modeling results. Habitat modeling indicated that the most suitable beaver habitat was near perennial water, near aspen (Populus tremuloides) and willow (Salix spp.), and in low gradient streams at lower elevations. Results from the initial monitoring period in October 2007 allowed us to assess costs and logistical considerations, validate our habitat model, and conduct power analyses to assess whether our sampling design could detect the level of declines in beaver stated in the monitoring objectives. Beaver food caches were located in 20 of 52 sampled watersheds. Monitoring 20 to 25 watersheds with beaver should provide sufficient power to detect 15–40% declines in the beaver food cache index as well as a twofold decline in the odds of beaver being present in watersheds. Indices of abundance, such as the beaver food cache index, provide a practical measure of population status to conduct long-term monitoring across broad landscapes such as national forests.     

Sampling The Soil In Long-Term Forest Plots: The Implications Of Spatial Variation

Long-term monitoring of forest soils as part of a pan-European network to detect environmental change depends on an accurate determination of the mean of the soil properties at each monitoring event. Forest soil is known to be very variable spatially, however. A study was undertaken to explore and quantify this variability at three forest monitoring plots in Britain. Detailed soil sampling was carried out, and the data from the chemical analyses were analysed by classical statistics and geostatistics. An analysis of variance showed that there were no consistent effects from the sample sites in relation to the position of the trees. The variogram analysis showed that there was spatial dependence at each site for several variables and some varied in an apparently periodic way. An optimal sampling analysis based on the multivariate variogram for each site suggested that a bulked sample from 36 cores would reduce error to an acceptable level. Future sampling should be designed so that it neither targets nor avoids trees and disturbed ground. This can be achieved best by using a stratified random sampling design.     

Assessment of wadeable stream resources in the driftless area ecoregion in Western Wisconsin using a probabilistic sampling design

The Wisconsin Department of Natural Resources (WDNR), with support from the U.S. EPA, conducted an assessment of wadeable streams in the Driftless Area ecoregion in western Wisconsin using a probabilistic sampling design. This ecoregion encompasses 20% of Wisconsin’s land area and contains 8,800 miles of perennial streams. Randomly-selected stream sites (n = 60) equally distributed among stream orders 1–4 were sampled. Watershed land use, riparian and in-stream habitat, water chemistry, macroinvertebrate, and fish assemblage data were collected at each true random site and an associated “modified-random” site on each stream that was accessed via a road crossing nearest to the true random site. Targeted least-disturbed reference sites (n = 22) were also sampled to develop reference conditions for various physical, chemical, and biological measures. Cumulative distribution function plots of various measures collected at the true random sites evaluated with reference condition thresholds, indicate that high proportions of the random sites (and by inference the entire Driftless Area wadeable stream population) show some level of degradation. Study results show no statistically significant differences between the true random and modified-random sample sites for any of the nine physical habitat, 11 water chemistry, seven macroinvertebrate, or eight fish metrics analyzed. In Wisconsin’s Driftless Area, 79% of wadeable stream lengths were accessible via road crossings. While further evaluation of the statistical rigor of using a modified-random sampling design is warranted, sampling randomly-selected stream sites accessed via the nearest road crossing may provide a more economical way to apply probabilistic sampling in stream monitoring programs.     

Sampling Scale Effects in Random Fields and Implications for Environmental Monitoring

The concept of a sampling scale triplet of spacing, extent and support is used to define the spatial dimensions of a monitoring network or a field study. The spacing is the average distance between samples, the extent is the size of the domain sampled and the support is the averaging area of one sample. The aim of this paper is to examine what is the bias and the random error (uncertainty) introduced by the sampling scale triplet into estimates of the mean, the spatial variance and the integral scale of a variable in a landscape. The integral scale is a measure of the average distance over which a variable is correlated in space. A large number of two dimensional random fields are generated from which hypothetical samples, conforming to a certain sampling scale triplet, are drawn which in turn are used to estimate the sample mean, spatial variance and integral scale. The results indicate that the biases can be up to two orders of magnitude. The bias of the integral scale is positively related to the magnitude of any of the components of the scale triplet while the bias of the spatial variance is different for different components of the scale triplet. All sampling scale effects are relative to the underlying correlation length of the variable of interest which is closely related to the integral scale. The integral scale can hence be used for sampling design and data interpretation. Suggestions are given on how to adjust a monitoring network to the scales of the variables of interest and how to interpret sampling scale effects in environmental data.     

The BIOTA Biodiversity Observatories in Africa—a standardized framework for large-scale environmental monitoring

The international, interdisciplinary biodiversity research project BIOTA AFRICA initiated a standardized biodiversity monitoring network along climatic gradients across the African continent. Due to an identified lack of adequate monitoring designs, BIOTA AFRICA developed and implemented the standardized BIOTA Biodiversity Observatories, that meet the following criteria (a) enable long-term monitoring of biodiversity, potential driving factors, and relevant indicators with adequate spatial and temporal resolution, (b) facilitate comparability of data generated within different ecosystems, (c) allow integration of many disciplines, (d) allow spatial up-scaling, and (e) be applicable within a network approach. A BIOTA Observatory encompasses an area of 1?km² and is subdivided into 100 1-ha plots. For meeting the needs of sampling of different organism groups, the hectare plot is again subdivided into standardized subplots, whose sizes follow a geometric series. To allow for different sampling intensities but at the same time to characterize the whole square kilometer, the number of hectare plots to be sampled depends on the requirements of the respective discipline. A hierarchical ranking of the hectare plots ensures that all disciplines monitor as many hectare plots jointly as possible. The BIOTA Observatory design assures repeated, multidisciplinary standardized inventories of biodiversity and its environmental drivers, including options for spatial up- and downscaling and different sampling intensities. BIOTA Observatories have been installed along climatic and landscape gradients in Morocco, West Africa, and southern Africa. In regions with varying land use, several BIOTA Observatories are situated close to each other to analyze management effects.     

Selecting the optimum plot size for a California design-based stream and wetland mapping program

Accurate estimates of the extent and distribution of wetlands and streams are the foundation of wetland monitoring, management, restoration, and regulatory programs. Traditionally, these estimates have relied on comprehensive mapping. However, this approach is prohibitively resource-intensive over large areas, making it both impractical and statistically unreliable. Probabilistic (design-based) approaches to evaluating status and trends provide a more cost-effective alternative because, compared with comprehensive mapping, overall extent is inferred from mapping a statistically representative, randomly selected subset of the target area. In this type of design, the size of sample plots has a significant impact on program costs and on statistical precision and accuracy; however, no consensus exists on the appropriate plot size for remote monitoring of stream and wetland extent. This study utilized simulated sampling to assess the performance of four plot sizes (1, 4, 9, and 16 km²) for three geographic regions of California. Simulation results showed smaller plot sizes (1 and 4 km²) were most efficient for achieving desired levels of statistical accuracy and precision. However, larger plot sizes were more likely to contain rare and spatially limited wetland subtypes. Balancing these considerations led to selection of 4 km² for the California status and trends program.     

Review of environmental monitoring methods: survey designs

During the past decade and a half, environmental monitoring programshave increased in numberand importance. Large scale environmental monitoring programs often present design difficulties because they tend to measure many(sometimes hundreds) of parameters through space and time. Thispaper reviewed and summarized one important component of environmental monitoring programs, the statistical survey design. Survey designsused for long-term monitoring programs lasting multiple (3)occasions were reviewed, paying specialattention to those published after 1985. During this review, two key components of the overall survey design were identified. The first key component wasthe membership design. Groups of population units sampled the same occasionwere called panels here, and the membership design specified which units were members of which panels. The second component was the revisit design that specified when panels were to be revisited. Membership designs varied, butsome form of simple random or systematic design was popular.Among revisit designs, four basicpatterns were found in the literature and their relative strengths andweaknesses were summarized. To efficiently discuss revisit designs, anew unified short-hand notation was proposed and adopted.     

Towards a sampling strategy for the assessment of forest condition at European level: combining country estimates

A correct characterization of the status and trend of forest condition is essential to support reporting processes at national and international level. An international forest condition monitoring has been implemented in Europe since 1987 under the auspices of the International Co-operative Programme on Assessment and Monitoring of Air Pollution Effects on Forests (ICP Forests). The monitoring is based on harmonized methodologies, with individual countries being responsible for its implementation. Due to inconsistencies and problems in sampling design, however, the ICP Forests network is not able to produce reliable quantitative estimates of forest condition at European and sometimes at country level. This paper proposes (1) a set of requirements for status and change assessment and (2) a harmonized sampling strategy able to provide unbiased and consistent estimators of forest condition parameters and of their changes at both country and European level. Under the assumption that a common definition of forest holds among European countries, monitoring objectives, parameters of concern and accuracy indexes are stated. On the basis of fixed-area plot sampling performed independently in each country, an unbiased and consistent estimator of forest defoliation indexes is obtained at both country and European level, together with conservative estimators of their sampling variance and power in the detection of changes. The strategy adopts a probabilistic sampling scheme based on fixed-area plots selected by means of systematic or stratified schemes. Operative guidelines for its application are provided.     

Integrated response plot designs for indicators of desertification

The improvement of land management practices on lands susceptible to desertification requires information on the status and condition of the existing resources as well as any change occurring in the resource condition over time. The Environmental Monitoring and Assessment Program (EMAP) of the U.S. Environmental Protection Agency has developed a statistical survey design for monitoring the condition of ecological resources on large spatial scales. EMAP-Rangelands used a uniformity sampling study in 1993 to evaluate response plot designs for three categories of indicators (soils, vegetation, and spectral reflectance) to be used for monitoring ecological condition of a site. The response plot design study was developed to integrate on-site measurements for the three indicator categories. The study was conducted on the Colorado Plateau in southern Utah in three rangeland resource classes (grassland, desertscrub, and conifer woodland) of differing productivity levels in an attempt to develop a common plot design for all three resource classes. Basic measurement units were developed to facilitate integration of data collection. Preliminary spatial analysis of the sampling study found considerable differences in variation patterns among the study sites and measurement categories for the indicator classes used by EMAP-Rangelands. Evidence of substantial trends in the indicator measurements on monitoring sites relative to regional trends leads to the conclusion that nonstationary spatial models for biological processes on a monitoring site may be needed to fulfill the requirements for developing plot designs and indicator criteria.The U.S. Environmental Protection Agency, through the Office of Research and Development, funded the research described here. This paper has been subjected to the Agency's peer and administrative review and has been approved as an EPA publication. The U.S. Government has the right to retain a nonexclusive, royalty-free license in and to any copyright covering this article.     

Double sampling for stratification for the monitoring of sparse tree populations: the example of Populus euphratica Oliv. forests at the lower reaches of Tarim River, Southern Xinjiang, China

Desertification is a pressing issue in the dry Tarim River basin, which is under anthropogenic stresses. In this study, double sampling for stratification (DSS) is employed to inventory Populus euphratica Oliv. forests in the lower reaches of the Tarim River Basin in Xinjiang, China. The two objectives were evaluating DSS as a sampling technique for monitoring desertification and generating baseline information for permanent observation. Here, DSS consists of two phases: in phase 1, crown cover is observed on a large sample of plots on a high resolution satellite image, and these photo-plots are stratified into five crown cover strata. Phase 2 is a stratified random sample from these photo-plots and the sampled plots are field observed. Approximately 32% of the study area is without P. euphratica trees. As expected, estimated mean poplar tree density and basal area increase with crown cover. DSS takes advantages of stratification (fieldwork efficiency and statistical precision) without the need for a priori strata delineation. It proves feasible for inventory the sparse poplar population and holds promise for the assessment of trees outside the forest, where density varies considerably and pre-stratification is intractable. It can be integrated into permanent observation systems for monitoring vegetation changes.     

Watershed-Based Survey Designs

Watershed-based sampling design and assessment tools help serve the multiple goals for water quality monitoring required under the Clean Water Act, including assessment of regional conditions to meet Section 305(b), identification of impaired water bodies or watersheds to meet Section 303(d), and development of empirical relationships between causes or sources of impairment and biological responses. Creation of GIS databases for hydrography, hydrologically corrected digital elevation models, and hydrologic derivatives such as watershed boundaries and upstream–downstream topology of subcatchments would provide a consistent seamless nationwide framework for these designs. The elements of a watershed-based sample framework can be represented either as a continuous infinite set defined by points along a linear stream network, or as a discrete set of watershed polygons. Watershed-based designs can be developed with existing probabilistic survey methods, including the use of unequal probability weighting, stratification, and two-stage frames for sampling. Case studies for monitoring of Atlantic Coastal Plain streams, West Virginia wadeable streams, and coastal Oregon streams illustrate three different approaches for selecting sites for watershed-based survey designs.     

Multiscale satellite and spatial information and analysis framework in support of a large-area forest monitoring and inventory update

Many countries undertake a national forest inventory to enable statistically valid monitoring in support of national and international reporting of forest conditions and change. Canada’s National Forest Inventory (NFI) program is designed to operate on a 10-year remeasurement cycle, with an interim report produced at the 5-year mid-point. The NFI is a sample-based inventory, with approximately 18,850 2 ×2-km photo plots across the country, distributed on a 20×20-km grid of sample points; these photo plots are the primary data source for the NFI. Capacity to provide annual monitoring information is required to keep policy and decision makers apprised of current forest conditions. In this study, we implemented a multistage monitoring framework and used a Moderate Resolution Imaging Spectroradiometer (MODIS) change product to successfully identify 78% of the changes in forest cover area that were captured with a Landsat change detection approach. Of the NFI photo plots that were identified by both the Landsat and MODIS approaches as having changes in forest cover, the proportion of change area within the plots was similar (R ²?=?0.78). Approximately 70% of the Landsat-derived change events occupied less than 40% of a single MODIS pixel, and more than 90% of the change events of this size were successfully detected with the MODIS product. Finally, MODIS estimates of the proportion of forest cover change at the NFI photo plot level were comparable to change estimates for the ecoregions as a whole (R ²?=?0.95). High-temporal, low-spatial resolution imagery such as MODIS, in combination with other remotely sensed data sources, can provide information on disturbance events within a national forest inventory remeasurement cycle, thereby satisfying the interim information needs of policy and decision makers as well as the requirements of national and international reporting commitments.     

Optimising the sampling effort in riparian surveys

Riparian condition is commonly measured as part of stream health monitoring programs as riparian vegetation provides an intricate linkage between the terrestrial and aquatic ecosystems. Field surveys of a riparian zone provide comprehensive riparian attribute data but can be considerably intensive and onerous on resources and workers. Our objective was to assess the impact of reducing the sampling effort on the variation in key riparian health indicators. Subsequently, we developed a non-parametric approach to calculate an information retained (IR) statistic for comparing several constrained systematic sampling schemes to the original survey. The IR statistic is used to select a scheme that reduces the time taken to undertake riparian surveys (and thus potentially the costs) whilst maximising the IR from the original survey. Approximate bootstrap confidence intervals were calculated to improve the inferential capability of the IR statistic. The approach is demonstrated using riparian vegetation indicators collected as part of an aquatic ecosystem health monitoring program in Queensland, Australia. Of the nine alternative sampling designs considered, the sampling design that reduced the sampling intensity per site by sixfold without significantly comprising the quality of the IR, results in halving the time taken to complete a riparian survey at a site. This approach could also be applied to reducing sampling effort involved in monitoring other ecosystem health indicators, where an intensive systematic sampling scheme was initially employed.     

Sampling for landscape elements—a case study from Lower Saxony,Germany

The estimation of coverage, i.e., the proportion of the total area in a study region covered by a given target class, is essential to many aspects of environmental monitoring. We analyze and compare the efficiency of different sample-based approaches for the estimation of coverage of different land cover classes from aerial imagery in a case study in Lower Saxony, Germany on the basis of the estimated standard errors. A complete delineation of vegetation classes in n?=?279 aerial photo plots of 400 × 400 m thrown onto the study region of 1,117.7 km² in accordance with a systematic grid is compared to different configurations of line intercept sampling and clusters of points. The observation designs under study are characterized by different complexity and total size of the observation units and therefore also to the efforts related to yield a single observation. Especially for those classes that cover a relatively large proportion of the sampling frame, our results show that difference in performance between the different designs are negligible. A cluster of four transects of 200 m each allows estimating the area of land cover classes with high coverage with nearly similar precision as a complete mapping of fixed area plots of 16 ha each. Clusters of points show unexpected high precision for the estimated coverage of land cover classes with relatively high coverage.     

Integrating probabalistic and targeted compliance monitoring for comprehensive watershed assessment

Environmental monitoring typically falls into one of two broad categories. Targeted designs, utilizing fixed stations, focus on describing and quantifying impacts, tracking trends, and assessing compliance with regulatory guidelines or limits. Probabilistic designs, in contrast, draw sampling stations at random from an area or region, and the stations are used to describe conditions in the region of interest based on a subpopulation of sites. These two design approaches are usually viewed as mutually exclusive, with randomized designs used for broader regional assessments of overall ambient condition and targeted designs for demonstrating regulatory compliance and/or characterizing specific, localized impacts. Combining elements of both approaches into a single design provides benefits not available from either design alone. Embedding targeted monitoring within the framework of a probabilistic design enables data from targeted stations to be viewed in a more accurate regional context and provides a consistent background against which to identify characteristic regional patterns of contamination and impact. We use the San Gabriel River Regional Monitoring Program, recently implemented in southern California, to illustrate the structure of a hybrid design and how it enables data analyses and assessments that provide a more complete picture of conditions in the watershed. For example, the hybrid design showed that approximately 80% of the metals levels at compliance sites were below the 25th percentile of the overall watershed condition as indicated by the probabilistic sampling.     

Effect of diel period and season on seining effort required to detect changes in Lake Erie beach fish assemblages

Using beach seining data from Lake Erie beaches, we examined the effect of diel period and season on two common study objectives considered sensitive to sampling effort: (1) has fish assemblage composition changed between two sampling time periods?, and (2) how many samples are necessary to provide a representative sample of the fish assemblage? Across a range of effect sizes and power levels, sample sizes were estimated for species richness, total fish abundance, round goby (Neogobius melanostomus) abundance and emerald shiner (Notropis atherinoides) abundance. For most sampling periods and effect sizes, several hundred to several thousand seine hauls were estimated to be required for detecting large-scale (e.g. basin-wide) changes in round goby and emerald shiner abundance. Depending on the indicator of interest, night seining (either spring or fall) was interpreted as the most efficient approach. For both seasons, sampling effort beyond 40 to 50 night seine hauls results in few additional species being detected. Both species accumulation curves and estimated sample sizes indicate that less effort is required to detect species richness changes than fish abundance. For most effect sizes and power levels, estimates of the number of seine hauls required to detect changes in fish abundance are unrealistically large for most monitoring programs. More practical monitoring strategies could be achieved by adopting alternative indicators (e.g. guilds), a more liberal significance level (p=0.1), a paired-site sampling design, or including reference sites.     

Assessing Spatial, Temporal, and Analytical Variation of Groundwater Chemistry in a Large Nuclear Complex, USA

Statistical analyses were applied at the Hanford Site, USA, to assess groundwater contamination problems that included (1) determining local backgrounds to ascertain whether a facility is affecting the groundwater quality and (2) determining a ‘pre-Hanford' groundwater background to allow formulation of background-based cleanup standards. The primary purpose of this paper is to extend the random effects models for (1) assessing the spatial, temporal, and analytical variability of groundwater background measurements; (2) demonstrating that the usual variance estimate s ², which ignores the variance components, is a biased estimator; (3) providing formulas for calculating the amount of bias; and (4) recommending monitoring strategies to reduce the uncertainty in estimating the average background concentrations. A case study is provided. Results indicate that (1) without considering spatial and temporal variability, there is a high probability of false positives, resulting in unnecessary remediation and/or monitoring expenses; (2) the most effective way to reduce the uncertainty in estimating the average background, and enhance the power of the statistical tests in general, is to increase the number of background wells; and (3) background for a specific constituent should be considered as a statistical distribution, not as a single value or threshold. The methods and the related analysis of variance tables discussed in this paper can be used as diagnostic tools in documenting the extent of inherent spatial and/or temporal variation and to help select an appropriate statistical method for testing purposes.     

Extensive tree health monitoring networks are useful in revealing the impacts of widespread biotic damage in boreal forests

We surveyed the regional distribution of conifer defoliation in Finland with an extensive monitoring network during 1995–2006 (EU Forest Focus Level I). The average defoliation in the whole Finland was 10.3% in pine and 19.9% in spruce. The sharp changes were often related to abiotic and biotic factors. The mean age of the stand explained more than one half of the between-plot variance in defoliation. In a variance component analysis, the main effect of years was negligible, while most of the random variation was due to plot main effect and plot × year interaction. About one fifth of the defoliation could be attributed to abiotic or biotic damage, and there were strong local correlations, e.g., between the changes in defoliation and degree of pine sawfly (Diprionidae) damage. There were clear temporal and spatial patterns in the incidence of the most important causes [Scots pine: Scleroderris canker (Gremmeniella abietina), pine shoot beetles (Tomicus sp.), and pine sawflies (Diprion pini, Neodiprion sertifer); Norway spruce: rust fungi (primarily Chrysomyxa ledi)]. Our results suggest that extensive monitoring networks can reveal useful information about the widespread outbreaks of pest organisms (insects and fungi) already in their increase phases, giving some time for management decisions. In a changing climate, large-scale, regular monitoring of tree health, including abiotic and biotic causes, is more important than ever before.     

Using Spatial Pattern to Quantify Relationship Between Samples, Surroundings, and Populations

The need for accurate carbon budgeting, climate change modelling, and sustainable resource management has lead to an increase in the number of large area forest monitoring programs. Large area forest monitoring programs often utilize field and remotely sensed data sources. Sampling, via field or photo plots, enables the collection of data with the desired level of categorical detail in a timely and efficient manner. When sampling, the aim is to collect representative detailed data enabling the statistical reporting upon the characteristics of larger areas. As a consequence, approaches for investigating how well sample data represent larger areas (i.e., the sample neighbourhood and the population) are desired. Presented in this communication is a quantitative approach for assessing the nature of sampled areas in relation to surrounding areas and the overall population of interest. Classified Landsat data is converted to forest/non-forest categories to provide a consistent and uniform data set over a 130,000 km² study region in central British Columbia, Canada. From this larger study area 322 2 × 2 km photo plots on a 20 × 20 km systematic grid are populated with composition and configuration information for comparison to non-sampled areas. Results indicate that typically, within the study area, the spatial pattern of forest within a photo plot is representative of the forest patterns found within primary and secondary neighbourhoods and over the entire population of the study. These methods have implications for understanding the nature of data used in monitoring programs worldwide. The ability to audit photo and field plot information promotes an increased understanding of the results developed from sampling and provides tools identifying locations of possible bias.