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.     

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.     

Patchy distribution fields: sampling distance unit of a zigzag survey and reconstruction adequacy

A mathematical model was used to examine the effects of choosingvarious units of sampling distance of a zigzag survey on the adequacy of reconstructing patchy distribution fields. The modelsimulates fish or plankton patches (or gaps) of different shapesand spatial orientations, and an acoustic survey by zigzag or parallel transects along which a unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). A priori information on the autocorrelation radii for the field in the directions of the survey (R_s) and perpendiculardirection (R_p) allows optimisation of the survey design andthe algorithm of data analysis. A field can be reconstructed properly (r² > 0.70) if the distance between transects D < (1.0–1.5)R_s and the unit of sampling distance d < (1.0–1.5)R_p. A posteriori determination of patchorientation allows reconstruction of the best field attainable onthe basis of the survey data. In cases of field movement, if thedimension of patches in the direction of movement exceeds that ofa surveyed area, a survey in the opposite direction gives bestresults; in contrast, if the dimension of moving patches is smaller than that of a surveyed area, it is reasonable to carryout a survey in the same direction. The criterion remains validwhen a survey is carried out by zigzag transects and a unit ofsampling distance is set along them. The results obtained indicate that, for a fixed transect spacing and a given number of sampling points on each full transect, zigzag pattern allowsless adequate reconstruction of an original distribution field (in cases of both immovable and movable fields) than corresponding parallel pattern.     

Patchy distribution fields: a zigzag survey design and reconstruction adequacy

A mathematical model was used to examine the effects of a zigzagsurvey 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 an acoustic survey by a set of transects forming a zigzag. Adequacyof the reconstructed fields to those originally generated was evaluated by calculating their correlations (r). A prioriinformation on the autocorrelation radius for the field in thedirection of a survey (R) allows optimization of survey design and the algorithm of data analysis. A patchy field can be reconstructed properly (r > 0.80) if the distance between transects D < (1.0–1.5)R. If a priori information onthe field is not available, the autocorrelation radius should bedetermined when reconstructing the field, i.e. a posteriori. In cases of field movement, the criterion for choosing a survey direction is based on the relationship betweenthe dimension of moving patches in the direction of movement andthat of the surveyed area. The results obtained indicate that, for a fixed transect spacing, zigzag pattern allow less adequatereconstruction of an original distribution field (in cases of both immovable and movable fields) than corresponding parallel pattern.     

Patchy Distribution Fields: An Interleaved Survey Design and Reconstruction Adequacy

A mathematical model was used to compare the effects of a regular (one-pass) or interleaved (two-pass) acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches of different shapes and spatial orientations, and a set of parallel or zigzag transects forming a regular or interleaved acoustic survey. 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. Regarding the immovable fields, the efficiency of a regular or interleaved acoustic survey was tested with the following two alternative assumptions: (1) the entire survey was completed; (2) the survey was interrupted for some reason at the moment when one transect remained non-accomplished. In the former case, the efficiencies of both acoustic surveys were nearly the same; in the latter case, the efficiency of an interleaved survey was superior to that of a regular one. With respect to movable fields, the efficiency of the completed interleaved surveys was even higher than that of the regular ones. Thus, the results obtained allow us to conclude that an interleaved survey is expedient in cases where there is no preference regarding the position of a vessel for further work.     

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.     

Evaluation of Terrestrial and Streamside Salamander Monitoring Techniques at Shenandoah National Park

In response to concerns about amphibian declines, a study evaluating and validating amphibian monitoring techniques was initiated in Shenandoah and Big Bend National Parks in the spring of 1998. We evaluate precision, bias, and efficiency of several sampling methods for terrestrial and streamside salamanders in Shenandoah National Park and assess salamander abundance in relation to environmental variables, notably soil and water pH. Terrestrial salamanders, primarily redback salamanders (Plethodon cinereus), were sampled by searching under cover objects during the day in square plots (10 to 35 m²). We compared population indices (mean daily and total counts) with adjusted population estimates from capture-recapture. Analyses suggested that the proportion of salamanders detected (p) during sampling varied among plots, necessitating the use of adjusted population estimates. However, adjusted population estimates were less precise than population indices, and may not be efficient in relating salamander populations to environmental variables. In future sampling, strategic use of capture-recapture to verify consistency of p's among sites may be a reasonable compromise between the possibility of bias in estimation of population size and deficiencies due to inefficiency associated with the estimation of p. The streamside two-lined salamander (Eurycea bislineata) was surveyed using four methods: leaf litter refugia bags, 1 m² quadrats, 50 × 1 m visual encounter transects, and electric shocking. Comparison of survey methods at nine streams revealed congruent patterns of abundance among sites, suggesting that relative bias among the methods is similar, and that choice of survey method should be based on precision and logistical efficiency. Redback and two-lined salamander abundance were not significantly related to soil or water pH, respectively.     

Judging Survey Quality: Local Variances

The present study was undertaken to explore possibilities to judge survey quality on basis of a limited and restricted number of a-priori observations. Here, quality is defined as the ratio between survey and local variance (signal-to-noiseratio). The results indicate that the presented surveys do not permit such judgement. Furthermore, the discussion alsosuggests that the 5-fold local sampling strategies do not merit any sound judgement in about 10% of all cases.The results further imply that surveys will benefit from procedures, controls and approaches in sampling and sample handling, to assess both average, variance and the nature of the distribution of elemental concentrations in local sites. This reasoning is compatible with the idea of the site as a basic homogeneous survey unit, which is implicitly and conceptually underlying any survey performed.     

Feeding of western gray whales during a seismic survey near Sakhalin Island, Russia

Exxon Neftegas Limited, as operator of the Sakhalin-1 consortium, is developing oil and gas reserves on the continental shelf off northeast Sakhalin Island, Russia. DalMorNefteGeofizika (DMNG) on behalf of the Sakhalin-1 consortium conducted a 3-D seismic survey of the Odoptu license area during 17 August–9 September 2001. A portion of the primary feeding area of the endangered western gray whale (Eschrichtius robustus) is located in the vicinity of the seismic survey. This paper presents data to assess whether western gray whale bottom feeding activity, as indicated by visible mud plumes, was affected by seismic operations. The mitigation and monitoring program associated with the seismic survey included aerial surveys during 19 July–19 November 2001. These aerial surveys documented the local and regional distribution, abundance, and bottom feeding activity of western gray whales. Data on gray whale feeding activity before, during and after the seismic survey were collected, with the whales assumed to be feeding on the benthos if mud plumes were observed on the surface. The data were used to assess the influence of seismic survey and other factors (including environmental) on feeding activity of western gray whales. A stepwise multiple regression analysis failed to find a statistically significant effect (α = 0.05) of the seismic survey on frequency of occurrence of mud plumes of western gray whales used as a proxy to evaluate bottom feeding activity in Piltun feeding area. The regression indicated that transect number (a proxy for water depth, related to distance from shore) and swell height (a proxy for sea state) were the only variables that had a significant effect on frequency of whale mud plumes. It is concluded that the 2001 seismic survey had no measurable effect (α = 0.05) on bottom feeding activity of western gray whales off Sakhalin Island. M.W. Newcomer, deceased     

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.     

Assessment of Temporal and Spatial Variation of Nitrate Removal in Riparian Zones

The aim of this study was to monitor long-term temporal and spatial groundwater NO₃^- removal efficiencies in different riparian zones via a limited number of sampling wells. Groundwater NO₃^- concentrations were measured fortnightly or monthly over a period of two years using transects of ground water sampling wells. Depending on the level of the NO₃^- load (up to 120mgNL^-1 at the input side of the riparian zone a distance of 10 to 30m was needed to remove NO₃^- from the groundwater below 11.3mgNL^-1. Considering all seasons, the mixed vegetation and grass riparian site succeeded to remove groundwater NO₃^- efficiently (92—100% within a distance of 30m. The forested riparian zone removed 72—90% of the total NO₃^- input within a distance of 30m. Evidence emerged that NO₃^- could also be removed actively at depths up to 2m, due to the presence of organically enriched layers of alluvial deposits or roots. Our four dimensional approach (three dimensional space and time), in combination with a limited number of sampling wells, was shown to be a useful monitoring tool to assess the variability of NO₃^- removal in riparian zones.     

Assessment of river water quality during snowmelt and base flow periods in two catchment areas with different land use

River water quality was evaluated with respect to eutrophication and land use during spring snowmelt and summer base flow periods in Abashiri (mixed cropland-livestock farming) and Okoppe (grassland-based dairy cattle farming), eastern Hokkaido, Japan. Water from rivers and tributaries was sampled during snowmelt and summer base flow periods in 2005, and river flow was measured. Total N (TN), NO₃–N, and Si concentrations were determined using standard methods. Total catchment and upland areas for each sampling site were determined with ArcGIS hydrology modeling software and 1:25,000-scale digital topographic maps. Specific discharge was significantly higher during snowmelt than during base flow. In both areas, TN concentrations increased, whereas Si concentrations decreased, with increased specific discharge, and were significantly higher during snowmelt. The Si:TN mole ratio decreased to below or close to the threshold value for eutrophication (2.7) in one-third of sites during snowmelt. River NO₃–N concentrations during base flow were significantly and positively correlated with the proportion of upland fields in the catchment in both the Abashiri (r = 0.88, P < 0.001) and Okoppe (r = 0.43, P < 0.01) areas. However, the regression slope, defined as the impact factor (IF) of water quality, was much higher in Abashiri (0.025) than in Okoppe (0.0094). The correlations were also significantly positive during snowmelt in both areas, but IF was four to eight times higher during snowmelt than during base flow. Higher discharge of N from upland fields and grasslands during snowmelt and the resulting eutrophication in estuaries suggest that nutrient discharge during snowmelt should be taken into account when assessing and monitoring the annual loss of nutrients from agricultural fields.     

A More Cost-Effective Emap Benthic Macrofaunal Sampling Protocol

Benthic macrofaunal sampling protocols in the U.S. Environmental Protection Agency's Environmental Monitoring and Assessment Program (EMAP) are to collect 30 to 50 random benthic macrofauna [defined as animals retained on a 0.5 mm (East and Gulf Coasts, USA) or a 1.0 mm mesh sieve (West Coast, USA)] samples per reporting unit using a 0.044 m² (East and Gulf Coasts) or 0.1 m² (West Coast) grab. Benthic macrofaunal community conditions in the reporting unit are characterized by cumulative distribution functions (CDFs) on end points of interest, such as number of species (S), abundance (A), and Shannon--Wiener diversity (H′). An EMAP and a companion field study were conducted concurrently in Tillamook Bay (Oregon, USA) to compare the cost effectiveness of benthic macrofauna samples collected using the EMAP West Coast (0.1 m² × ≥7 cm deep, 1.0 mm mesh), a 0.01 m² × 5 cm deep, 1.0 mm mesh, and a 0.01 m² × 5 cm deep, 0.5 mm mesh sampling protocol. Cost was estimated in relative laboratory sample-processing time. Sampling protocols were judged equally effective for EMAP purposes if, after linear transformation to adjust for scale changes in end point distributions, their S, A, and H′ CDFs were not significantly different. The 0.01 m² × 5 cm deep, 1.0 mm mesh sampling protocol was the most cost effective.     

Terrestrial vertebrate fauna surveys for the preparation of environmental impact assessments; how can we do it better? A Western Australian example

The Western Australian Environmental Protection Authority (EPA) in 2002 released Position Statement, No. 3, Terrestrial Biological Surveys as an Element of Biodiversity Protection outlining how terrestrial fauna survey data are to be used and interpreted in the preparation of environmental impact assessments (EIA). In 2004, the EPA released its Guidance for the Assessment of Environmental Factors, Terrestrial Fauna Surveys for Environmental Impact Assessment in Western Australia, No. 56. This paper briefly assesses the adequacy of recent terrestrial fauna surveys undertaken to support publicly released EIAs and indicates that the EPA is not always adhering to its own position and guidance statements. This paper argues that the current fauna survey guidelines are in need of improvement. The approach and requirements of some other Australian states are briefly assessed to identify similarities and where improvements can be made to the Western Australian (WA) guidelines. This paper concludes with suggestions on how the process and the guidelines in WA can be revised to more adequately assess the impact of developments on terrestrial vertebrate biodiversity and ecosystem function. These suggestions may have relevance for other areas where fauna surveys are undertaken to support EIAs.     

Within-and among-clutch variation of organochlorine residues in eggs of black-crowned night-herons

Within-clutch variability of DDE and PCB residues in eggs from 62 clutches of black-crowned night-herons (Nycticorax nycticorax) was small (12% and 17%) compared to among-clutch variability (88% and 83%). Significant correlations between concentrations of DDE (median r=0.8885) and of PCBs (median r=0.8244) occurred when 501 correlations were run on two randomly selected eggs from within the same clutch; no significant correlation occurred for either concentrations of DDE (median r=0.0353) or PCBs (median r=–0.0843) when eggs were not restricted to the same clutch but were restricted to the same colony. The probability of finding infrequently detected organochlorine contaminants (e.g., DDT, cis-chlordane) in eggs from the same clutch varied from 43–96% and increased as the chemical became more prevalent and the number of eggs per clutch became smaller. These results further support one of the basic assumptions of the sample egg technique, that the chemical residues in one egg in a clutch accurately reflect residues in the remaining eggs of the clutch.     

Predicting Water Quality Impaired Stream Segments using Landscape-Scale Data and a Regional Geostatistical Model: A Case Study in Maryland

In the United States, probability-based water quality surveys are typically used to meet the requirements of Section 305(b) of the Clean Water Act. The survey design allows an inference to be generated concerning regional stream condition, but it cannot be used to identify water quality impaired stream segments. Therefore, a rapid and cost-efficient method is needed to locate potentially impaired stream segments throughout large areas. We fit a set of geostatistical models to 312 samples of dissolved organic carbon (DOC) collected in 1996 for the Maryland Biological Stream Survey using coarse-scale watershed characteristics. The models were developed using two distance measures, straight-line distance (SLD) and weighted asymmetric hydrologic distance (WAHD). We used the Corrected Spatial Akaike Information Criterion and the mean square prediction error to compare models. The SLD models predicted more variability in DOC than models based on WAHD for every autocovariance model except the spherical model. The SLD model based on the Mariah autocovariance model showed the best fit (r² = 0.72). DOC demonstrated a positive relationship with the watershed attributes percent water, percent wetlands, and mean minimum temperature, but was negatively correlated to percent felsic rock type. We used universal kriging to generate predictions and prediction variances for 3083 stream segments throughout Maryland. The model predicted that 90.2% of stream kilometers had DOC values less than 5 mg/l, 6.7% were between 5 and 8 mg/l, and 3.1% of streams produced values greater than 8 mg/l. The geostatistical model generated more accurate DOC predictions than previous models, but did not fit the data equally well throughout the state. Consequently, it may be necessary to develop more than one geostatistical model to predict stream DOC throughout Maryland. Our methodology is an improvement over previous methods because additional field sampling is not necessary, inferences about regional stream condition can be made, and it can be used to locate potentially impaired stream segments. Further, the model results can be displayed visually, which allows results to be presented to a wide variety of audiences easily.     

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.     

The monitoring results of electromagnetic radiation of 110-kV high-voltage lines in one urban location in Chongqing P.R. China

To understand electromagnetic radiation field strength and its influencing factors of certain 110-kV high-voltage lines in one urban area of Chongqing by measuring 110-kV high-voltage line’s electromagnetic radiation level. According to the methodology as determined by the National Hygienic Standards, we selected certain adjacent residential buildings, high-voltage lines along a specific street and selected different distances around its vertical projection point as monitoring points. The levels of electromagnetic radiations were measured respectively. In this investigation within the frequency of 5–1,000 Hz both the electric field strength and magnetic field strength of each monitoring sites were lower than the public exposure standards as determined by the International Commission on Non-Ionizing Radiation Protection. However, the electrical field strength on the roof adjacent to the high-voltage lines was significantly higher than that as measured on the other floors in the same buildings (p < 0.05). The electromagnetic radiation measurements of different monitoring points, under the same high-voltage lines, showed the location which is nearer the high-voltage line maintain a consistently higher level of radiation than the more distant locations (p < 0.05). Electromagnetic radiation generated by high-voltage lines decreases proportionally to the distance from the lines. The buildings can to some extent shield (or absorb) the electric fields generated by high-voltage lines nearby. The electromagnetic radiation intensity near high-voltage lines may be mitigated or intensified by the manner in which the high-voltage lines are set up, and it merits attention for the potential impact on human health.     

Estimating Field Volatility of Soil Fumigants Using CHAIN_2D: Mitigation Methods and Comparison Against Chloropicrin and 1,3-Dichloropropene Field Observations

Academic, government, and industrial field researchers have generated a significant database of field studies of the volatility of soil applied fumigants. However, limited work exists in validating physical models against field volatility data sets and fully exploring the volatility parametric response surface. Field studies quantifying atmospheric flux for soil fumigants 1,3-dichloropropene and chloropicrin are validated against the United States Department of Agriculture (USDA Salinity Laboratory) soil physics model CHAIN_2D that was modified specifically for agronomic uses of soil fumigants. Comparison between model predictions and field observations for six unique field trials in five different states indicate that CHAIN_2D effectively captures the magnitude and duration of fumigant emission from soil observed experimentally with r ²  ∼ 0.14–0.96 (avg. 0.66) for peak emission, and r ²  ∼ 0.76–1.0 (avg. 0.91) for cumulative emissions (r equals the Pearson moment correlation coefficient). Correct prediction magnitudes suggest that CHAIN_2D is a useful tool for extrapolating flux predictions to diverse scenarios not addressed by field trials. Examples of mitigation strategies such as the use of agricultural films (tarps), increased soil injection depth, and management of soil water content, under near semi-infinite parameter combinations of soil, meteorological, and agronomic properties are discussed.     

Optimum macrobenthic sampling protocol for detecting pollution impacts in the Southern California Bight

The optimum macrobenthic sampling protocol [sampling unit, sieve mesh size, and sample size (n)] was determined for detecting ecologically important pollution impacts in the Southern California Bight, U.S.A. Cost, in laboratory processing time, was determined for samples obtained using fourteen sampling units (0.005–0.1 m² surface area) and two sieve mesh sizes (1.0 and 0.5 mm). Statistical power analyses for t-tests of means were performed to estimate the minimum sample size (n _min) needed to reliably (=0.05, 1–0.95) reject the null hypothesis of no difference between a reference and both a stimulated and a degraded station on twelve measures of community structure. The optimum sampling protocol for detecting impacts was determined as that with the lowest total cost ×n _min on most measures.Five replicate, 0.02 m²×5 cm deep, 1.0 mm mesh samples per station could reliably distinguish reference from impacted conditions on nine or ten measures of community structure at less than one quarter of the cost of the standard sampling protocol of 5 replicate, 0.1 m², 1.0 mm mesh samples per station. About 5 replicate, small (<0.1 m²), 1.0 mm mesh samples per station may often be optimal for detecting important structural changes in macrobenthic communities with naturally high species richness and abundance.