Patchy Distribution Fields: Sampling Distance Unit of a Spiral Survey and Reconstruction Adequacy

A mathematical model is used to examine the effects of choosing various units of sampling distance of a spiral acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and an acoustic survey by a spiral of Archimedes along which a unit of sampling distance is set. For comparison, surveys are imitated by parallel transects along which the same unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). In the case of a spiral survey, the mathematical experiments conducted show that an immovable field can be reconstructed properly (r ² > 0.70) if the ratio of the units of sampling distance to the autocorrelation radius for the field averaged in various directions d/R _av < 1.0–1.5. Regarding immovable fields, the spiral surveys ensure, practically speaking, the same adequacy of the field reconstruction as do surveys by parallel transects with the same unit of sampling distance. In regard to movable fields, a comparison of the results of spiral surveys with those of surveys by parallel transects indicates that the former may ensure even higher adequacy of the field reconstruction than do the latter, provided that the units of sampling distances in these surveys are equal to each other.     

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.     

Patchy distribution fields: A spiral survey design and reconstruction adequacy

A mathematical model was used to examine the effects of a spiral survey design on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and acoustic surveys by a spiral of Archimedes; for comparison, surveys by parallel or zigzag transects are imitated. Adequacy of the reconstructed fields to those originally generated was evaluated by calculating their correlations (r). The mathematical experiments conducted showed that spiral surveys ensure, practically speaking, the same adequacy of field reconstruction (both in cases of immovable or movable patches) as do surveys by parallel or zigzag transects with greater sampling effort (overall path). In the case of a spiral survey, a patchy field can be reconstructed properly (r² = 0.70) if the overall survey path is not less than S/R_av = 20–30, where R_av is the autocorrelation radius averaged for various directions. Thus, the results obtained allow us to conclude that a spiral survey design is expedient in cases that the minimal duration of a survey is a decisive factor for its conduction and there is a priori information that no onshore-offshore gradients of fish density exist in a region under study.     

Measuring differences between overlapping but unpaired spatial surveys using a geographic information system

To compare the results from two overlapping but unpaired synoptic surveys, a method using a geographic information system is proposed. The new method uses a GIS to create (i) minimum convex polygons (MCP) enclosing each set of survey data and (ii) Voronoi tesselations assigning area weights to each datum. The pairs of MCP and Voronoi maps are overlaid to produce a mosaic of polygons each with one datum assigned from each survey. The differences between the pairs weighted by the polygon areas provides the basis for statistical testing. Area-weighted means and variances of paired differences are computed and a z-statistic measures the significance of differences for the whole intersection area. A cross-product autocorrelation statistic provides an assessment of the spatial distribution of differences. Alternate, conventional methods are compared with the new method: analysis of variance (ANOVA), analysis of covariance (ANCOVA), and contour overlays. The four methods were applied to the comparison of two macrobenthic surveys conducted in Lake Erie. Measures for five different taxa were examined. The ANOVA and ANCOVA methods found many significant differences between surveys. Both methods were judged inappropriate as sampling data are not expected to be drawn from normally-distributed populations in spatial surveys. Differences between surveys were detected but were difficult to assess using contour overlays. The new method only found a significant difference in one taxon, an expected difference as members of that taxon were absent from most sites in one survey.     

Sampling of Sparse Species with Probability Proportional to Prediction

In this article, a 3P sampling (Probability Proportional to Prediction) approach is presented for surveying sparse species connected to certain types of substrates. The method uses the surveyors judgement of the probability of finding the species on a substrate as the base for selection of substrates for species inventories. The method is presented together with estimators, variances and variance estimators. The method is first presented for sampling based on all substrates in a study area and then as a subsampling technique in a two-stage design in which plots or strips are selected in a first stage. The presented approach was evaluated as a subsampling technique in a strip survey of calicioid lichen species associated with coarse broadleaf trees. In comparison with a simple random subsampling without replacement, 3P subsampling was in one study area found to be an improvement with 30–55% in terms of standard errors. The improvement was more modest in the other study area, only between 10–16%. The strip survey with 3P subsampling was more cost-efficient than a strip survey without subsampling except in one case. Based on the results in the test, 3P sampling seems to have a potential for sampling sparse species.     

Patchy Distribution Fields: Sampling Distance Unit of an Interleaved Survey and Reconstruction Adequacy

A mathematical model was used to examine the effects of choosing various units of sampling distance of an interleaved (two-pass) acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and an interleaved survey by parallel or zigzag transects, along which a unit of sampling distance is set. The efficiency of a survey is determined by the adequacy of a reconstructed field to that originally generated, which is evaluated by calculating their correlations (r). Regarding immovable fields, the experiments conducted indicate that a patchy field can be reconstructed properly (r ² > 0.70) if the distance between transects D < (1.5–2.0)R and the unit of sampling distance d < (1.5– 2.0)R _p. As, for regular surveys d < (1.0–1.5)R _p, it may be concluded that interleaved surveys are more efficient than regular ones because of the factor studied. In regard to movable fields, a comparison of the results of interleaved surveys with those of regular surveys directly indicates that the former may ensure a more adequate field reconstruction than the latter do. This fact confirms the previous conclusion that an interleaved survey is expedient in cases where there is no preference regarding the position of a vessel for further work.     

Using variance components to estimate power in a hierarchically nested sampling design

We used variance components to assess allocation of sampling effort in a hierarchically nested sampling design for ongoing monitoring of early life history stages of the federally endangered Devils Hole pupfish (DHP) (Cyprinodon diabolis). Sampling design for larval DHP included surveys (5 days each spring 2007–2009), events, and plots. Each survey was comprised of three counting events, where DHP larvae on nine plots were counted plot by plot. Statistical analysis of larval abundance included three components: (1) evaluation of power from various sample size combinations, (2) comparison of power in fixed and random plot designs, and (3) assessment of yearly differences in the power of the survey. Results indicated that increasing the sample size at the lowest level of sampling represented the most realistic option to increase the survey’s power, fixed plot designs had greater power than random plot designs, and the power of the larval survey varied by year. This study provides an example of how monitoring efforts may benefit from coupling variance components estimation with power analysis to assess sampling design.     

Aerial Monitoring of Marine Waterfowl in the Alaskan Beaufort Sea

The purpose of this study was to design and test a monitoring protocol for marine waterfowl in the central Alaskan Beaufort Sea. The study provides an important case-study of how a long-term monitoring program may be affected by unanticipated human disturbances.Because of its overwhelming and widespread abundance, relatively sedentary behavior, ease in counting, and the extensive historical database, the long-tailed duck (Clangula hyemalis) was selected as the focal species. Two null hypotheses were formulated concerning potential changes in the numbers and distribution of long-tailed ducks in relation to disturbance in an industrial study area, compared to a reference study area located about 50 km to the east.A 9-year historical database (1977–1984, 1989) of long-tailed duck densities and other important data recorded during systematic aerial surveys was analyzed retrospectively using multiple regression techniques. The retrospective analyses determined which of several predictor variables recorded were significantly related to long-tailed duck density. Separate analyses were conducted for two periods: (1) the overall period when long-tailed ducks were present in the lagoon study areas, and (2) the shorter adult male molt period. The results of the two analyses indicated that 57% and 68%, respectively, of the total variation in long-tailed duck density during the two periods could be explained by variables recorded during the surveys. Predictor variables representing habitat, day of the year, time of day, amount of ice, and wave height recorded on-transect during surveys were most closely associated with long-tailed duck density. Measurement error during the surveys, and influences outside the study area such as nesting success in tundra habitats and mortality during migration and in over-wintering areas likely also had strong influences on the results, but these factors were not measurable in our study.Based on results of the retrospective analyses, a long-term monitoring protocol consisting of a program of systematic aerial surveys and an analyses of variance and covariance (ANOVA and ANCOVA) statistical procedure was designed and initially tested in 1990 and 1991. This 2-year testing phase resulted in several revisions to the monitoring protocol. Refinements were made to the original sampling procedures, to the survey schedule, and to the recommended statistical analysis procedures. Results of the ANOVA and ANCOVA indicated that there was no evidence of a change in long-tailed duck densities that could be attributable to disturbance (from any source) in the industrial study area relative to a reference area with no industrial development. Other analyses indicated that the sampling and analysis procedures would be adequate to detect long-term trends in long-tailed duck density and localized disturbance effects, but that the monitoring program should be continued well beyond two years to detect statistically significant changes. As a result, additional aerial surveys of both study areas were conducted again during 1999–2001.Results of the revised ANOVA and ANCOVA of the 1990–1991 and 1999–2001 survey data indicated that the density of long-tailed ducks had significantly declined in coastal lagoons along the central Alaskan Beaufort Sea coast during the study period. In addition, disturbances throughout the barrier island-lagoon systems used by these ducks, including both the industrial and the reference study areas, had significantly increased over the same period. However, because unanticipated disturbances from a variety of anthropogenic sources, and not just industry sources, increased in both study areas, the reference study area was not an effective statistical control. As a result, the decline in long-tailed duck density in both study areas was not attributable to industry-related activities. Although the monitoring protocol described here is an effective method to detect statistically significant changes in long-tailed duck distribution and abundance in the nearshore Alaskan Beaufort Sea, many more years of sampling would be necessary to attribute observed changes to industry-related disturbances.     

Spatial scales of variation in lichens: implications for sampling design in biomonitoring surveys

The variability of biological data is a main constraint affecting the quality and reliability of lichen biomonitoring surveys for estimation of the effects of atmospheric pollution. Although most epiphytic lichen bioindication surveys focus on between-site differences at the landscape level, associated with the large scale effects of atmospheric pollution, current protocols are based on multilevel sampling, thus adding further sources of variation and affecting the error budget. We test the hypothesis that assemblages of lichen communities vary at each spatial scale examined, in order to determine what scales should be included in future monitoring studies. We compared four sites in Italy, along gradients of atmospheric pollution and climate, to test the partitioning of the variance components of lichen diversity across spatial scales (from trunks to landscapes). Despite environmental heterogeneity, we observed comparable spatial variance. However, residuals often overcame between-plot variability, leading to biased estimation of atmospheric pollution effects.     

Assessing differences in macrofaunal assemblages as a factor of sieve mesh size,distance between samples,and time of sampling

With increasing cascading effects of climate change on the marine environment, as well as pollution and anthropogenic utilization of the seafloor, there is increasing interest in tracking changes to benthic communities. Macrofaunal surveys are traditionally conducted as part of pre-incident environmental assessment studies and post-incident monitoring studies when there is a potential impact to the seafloor. These surveys usually characterize the structure and/or spatiotemporal distribution of macrofaunal assemblages collected with sediment cores; however, many different sampling protocols have been used. An assessment of the comparability of past and current survey methods was in need to facilitate future surveys and comparisons. This was the aim of the present study, conducted off the Oregon coast in waters 25–35 m deep. Our results show that the use of a sieve with a 1.0-mm mesh size gives results for community structure comparable to results obtained from a 0.5-mm mesh size, which allows reliable comparisons of recent and past spatiotemporal surveys of macroinfauna. In addition to our primary objective of comparing methods, we also found interacting effects of seasons and depths of collection. Seasonal differences (summer and fall) were seen in infaunal assemblages in the wave-induced sediment motion zone but not deeper. Thus, studies where wave-induced sediment motion can structure the benthic communities, especially during the winter months, should consider this effect when making temporal comparisons. In addition, some macrofauna taxa-like polychaetes and amphipods show high interannual variabilities, so spatiotemporal studies should make sure to cover several years before drawing any conclusions.     

A Multilevel Modeling Approach to Assessing Regional and Local Landscape Features for Lake Classification and Assessment of Fish Growth Rates

The ecoregion and watershed frameworks are landscape-based classifications that have been used to group waterbodies with respect to measures of community structure; however, they have yet to be evaluated for grouping lakes for demographic characteristics of fish populations. We used a multilevel modeling approach to determine if variability in mean fish length at age could be partitioned by ecoregions and watersheds. For the ecoregions analysis, we then examined if within-ecoregion variability could be explained by local water quality and lake morphometry characteristics. We used data from agency surveys conducted during 1974-1984 for age 2 and 3 fish of seven common warm and coolwater fish species. Variance in mean length at age between ecoregions for all species was not significant, and between-watershed variance estimates were only significant in 3 out of 14 analyses; however, the total amount of variation between watersheds was very small (ranging from 1.8% to 3.7% of the total variance), indicating that ecoregions and watersheds were ineffective in partitioning variability in mean length at age. Within ecoregions, water quality and lake morphometric characteristics accounted for 2%-23% of the variation in mean length at age. Measures of lake productivity were the most common significant covariates, with mean length at age increasing with increasing lake productivity. Much of the variability in mean length at age was not accounted for, suggesting that other local factors such as biotic interactions, fish density, and exploitation are important. The results indicate that the development of an effective regional framework for managing inland lakes will require a substantial effort to understand sources of demographic variability and that managers should not rely solely on ecoregions or watersheds for grouping lakes with similar growth rates.     

Comparison between electrofishing and snorkeling surveys to describe fish assemblages in Laurentian streams

We evaluated the relative performance of electrofishing and visual surveys (snorkeling) for estimating the abundance of combinations of fish species and size classes in rivers. We also assessed the effect of environmental conditions on potential differences between the results obtained using these two sampling methods. Sampling sites were distributed in the Laurentian region of Québec. Both methods were used while sections were blocked. Three snorkelers swam the river sections upstream while identifying and counting fish of each species and size classes. Three-pass electrofishing was performed in the same sites and abundances were estimated with a maximum likelihood depletion model. Greater abundances of fish were observed by snorkeling than by electrofishing at all sites. Snorkeling species richness was higher or equal to electrofishing richness in, respectively, 60 % and 40 % of sampled sites. Differences in the fish communities observed by both sampling methods were not related to environmental conditions. The results of our work are therefore contrary to that of most published studies that suggested the use of electrofishing over visual surveys. This study highlights that conclusions derived from previous work on sampling gear comparisons may not be generalisable; rather survey methods might benefit from being selected based on fish community composition.     

Geochemistry in the modern soil survey program

Elemental analysis has played an important role in the characterization of soils since inception of the soil survey in the US. Recent efforts in analysis of trace and major elements (geochemistry) have provided necessary data to soil survey users in a variety of areas. The first part of this paper provides a brief overview of elemental sources, forms, mobility, and bioavailability; critical aspects important to users of soil survey geochemical data for appropriate use and interpretations. Examples are provided based on data gathered as part of the US soil survey program. The second part addresses the organization of sample collection in soil survey and how soil surveys are ideally suited as a sampling strategy for soil geochemical studies. Geochemistry is functional in characterization of soil types, determining soil processes, ecological evaluation, or issues related to soil quality and health, such as evaluating suitability of soils for urban or agricultural land use. Applications of geochemistry are on-going across the US and are documented herein. This analytical direction of soil survey complements historic efforts of the National Cooperative Soil Survey Program and addresses the increasing need of soil survey users for data that assists in understanding the influence of human activities on soil properties.     

Long-Term Biomonitoring with Lichens: Comparing Data from Different Sampling Procedures

The results of 2 biomonitoring surveys carried out over a span of 7 years (1996–2003) in the Province of Ancona (central Italy) using the diversity of epiphytic lichens are reported. According to a recent standardisation process, in 2003 a different protocol was adopted, raising the problem of data comparability with the previous study. For this reason, a ‘pilot study' was carried out in which the two sampling tactics were both used. The high correlation between biodiversity values proved by the results of this study, let to compare the two lichen diversity data set and to identify the impacts on lichen communities of temporal variations of atmospheric quality in the study area, showing a reduction of stations classified as ‘altered' and an increase in stations classified as ‘semi-altered'. Furthermore, possible reduction of the sampling effort was detected performing a study on the sampling density of the survey area, indicating in a 6 × 6 km grid the best density to adopt in order to save financial resources for the future.     

Application of statistical modeling to optimize a coastal water quality monitoring program

The long-term water quality monitoring program implemented by the Massachusetts Water Resources Authority in 1992 is extensive and has provide substantial understanding of the seasonality of the waters in both Boston Harbor and Massachusetts Bay and the response to improvements in effluent quality and offshore transfer of the effluent in September 2000. The monitoring program was designed with limited knowledge of spatial and temporal variability and long-term trends within the system. This led to an extensive spatial and temporal sampling program. The data through 2003 showed high correlation within physical parameters measured (e.g., salinity, dissolved oxygen) and in biological measures such as chlorophyll fluorescence. To address the potential sampling redundancies in the measurement program, an assessment of the impact of reduced levels of monitoring on the ability to make water quality decisions was completed. The optimization was conducted by applying statistical models that took into account whether there was evidence of a seasonal pattern in the data. The optimization used model survey average readings to identify temporal fixed effects, model survey-average-corrected individual station readings to identify spatial fixed effects, corrected the individual station readings for temporal and spatial fixed effects and derived a correlation model for the corrected data, and applied the correlation model to characterize the correlation of annual average readings from reduced monitoring programs with true parameter levels. Reductions in the number of sampling stations were found less detrimental to the quality of the data for annual decision-making than reductions in the number of surveys per year, although there is less of a difference in this regard for dissolved oxygen than there is for chlorophyll. The analysis led to recommendations for a substantially lower monitoring effort with minimal loss of information. The recommendation supported an annual budget savings of approximately $183,000. Most of the savings was from fewer surveys as approximately $21,000 came from the reduction in the number of stations monitored from 21 to 7 and associated laboratory analytical costs.     

Improving the precision of sample-based forest damage inventories through two-phase sampling and post-stratification using remotely sensed auxiliary information

Many countries have a national forest inventory (NFI) designed to produce statistically sound estimates of forest parameters. However, this type of inventory may not provide reliable results for forest damage which usually affects only small parts of the forest in a country. For this reason, specially designed forest damage inventories are performed in many countries, sometimes in coordination with the NFIs. In this study, we evaluated a new approach for damage inventory where existing NFI data form the basis for two-phase sampling for stratification and remotely sensed auxiliary data are applied for further improvement of precision through post-stratification. We applied Monte Carlo sampling simulation to evaluate different sampling strategies linked to different damage scenarios. The use of existing NFI data in a two-phase sampling for stratification design resulted in a relative efficiency of 50 % or lower, i.e., the variance was at least halved compared to a simple random sample of the same size. With post-stratification based on simulated remotely sensed auxiliary data, there was additional improvement, which depended on the accuracy of the auxiliary data and the properties of the forest damage. In many cases, the relative efficiency was further reduced by as much as one-half. In conclusion, the results show that substantial gains in precision can be obtained by utilizing auxiliary information in forest damage surveys, through two-phase sampling, through post-stratification, and through the combination of these two approaches, i.e., post-stratified two-phase sampling for stratification.     

Evaluation of significant sources influencing the variation of water quality of Kandla creek, Gulf of Katchchh, using PCA

To evaluate the significant sources contributing to water quality parameters, we used principal component analysis (PCA) for the interpretation of a large complex data matrix obtained from the Kandla creek environmental monitoring program. The data set consists of analytical results of a seasonal sampling survey conducted over 2 years at four stations. PCA indicates five principal components to be responsible for the data structure and explains 76% of the total variance of the data set. The study stresses the need to include new parameters in the analysis in order to make the interpretation of principal components more meaningful. The PCA could be applied as a useful tool to eliminate multi-collinearity problems and to remove the indirect effect of parameters.     

Choosing an appropriate soil-gas survey method

Three types of soil-gas surveys were performed over the same location at an industrial facility. The surveys performed were (1) a passive technique using buried collectors, (2) an active technique in which samples were collected in the field by pumping with subsequent analysis at a field laboratory, and (3) a screening technique that employed a hand-held survey instrument. Each of the surveys had desirable features, but there was considerable variability in the quality of data obtained, cost, and ease of performance. The passive and active surveys provided excellent detection limits, but the passive technique does not provide data in real time. The active survey was also time-consuming and required an analytical chemist in the field. Neither the active nor the passive survey were able to delineate the plume of contaminated groundwater. Consequently, the results from the rapid and inexpensive reconnaissance technique provided functionally equivalent data.Publication No. 3672, Environmental Sciences Division, ORNL.Operated by Martin Marietta Energy Systems, Inc. for the U.S. Department of Energy under Contract No. DE-AC05-840R21400.