Cautions on using the Before-After-Control-Impact design in environmental effects monitoring programs
Abstract
Often the Before-After-Control-Impact (BACI) design is suggested as being a statistically powerful experimental design in environmental impact studies. If the timing and location of the impact are known and adequate pre-data are collected, the BACI design is considered optimal to help isolate the effect of the development from natural variability. This paper presents 9 years of results from a long-term BACI experiment tested using a range of statistical models and post-impact monitoring designs. To explore suboptimal designs that are often utilized in environmental effects monitoring, the same data were also explored assuming either no control system was available (Before-After only), or that no pre-impact data were available (Control-Impact only). The results of the BACI design were robust to the statistical model used, and the BACI design was able to detect effects from the impact that the two suboptimal designs failed to detect. However, the BACI design demonstrated different conclusions depending on the number and configuration of post-impact years included in the analysis. Our results reinforce the idea that caution should be employed when using, or interpreting results from, a BACI design in an environmental impact study, but demonstrate that a well-designed BACI remains one of the best models for environmental effects monitoring programs.
Introduction
Environmental effects monitoring design and analyses have been highlighted in the scientific literature for decades (Green 1979; Stewart-Oaten et al. 1986; Underwood 1991; Chapman 1998). Many challenges surround the ability to separate a human-induced change from natural variability, and addressing this real-world problem has been the focus of many sampling designs and statistical analyses. Green (1979) outlines a hierarchical tree of design options dependent on if the location and timing of the impact are known and/or if a control area is available. If a spatial control is lacking, the effects must be inferred from sampling both Before and After an impact (a design termed Intervention Analysis by Stewart-Oaten and Bence 2001), but this design necessarily assumes that an unaffected site would not have changed in a similar fashion. If the impact occurred before any monitoring began, the Impact site can be compared to a nearby (or upstream) Control site, but this design is suboptimal because it assumes the Control and Impact sites were similar prior to the impact. If the timing and location of the impact are unknown, the Reference Condition Approach has been suggested as a method for quantifying impairment of an affected site, but specific and rigorous guidelines must be met for this approach to be considered sound (Bowman and Somers 2005). If the impact has not yet occurred, but the timing and location are known, the Before-After-Control-Impact (BACI) design was suggested as optimal (Green 1979). However, since the BACI design was first introduced by Green (1979) many criticisms and suggestions for improvement have been made.
The first criticism was from Hurlbert (1984) who suggested that the use of a single Control area and a single Impact area (with replicate samples) analyzed by analysis of variance (ANOVA) only demonstrated significant differences between locations, and not of the impact itself because the treatment could not be randomly assigned to the experimental plots. Stewart-Oaten et al. (1986) refuted Hurlbert’s (1984) claim by suggesting that using sampling times as replicates with the Control and Impact sites sampled (near) simultaneously, and each sampling time represented by only one number (as the difference between the Impact and Control samples for that time), allows for the detection of the mean difference between time periods, using any standard two factor test. They do stress, however, that the assumptions of additivity and independence must be met in the data within a time period (Stewart-Oaten et al. 1986). If differences are not used the equivalent BACI test is the test for an interaction effect in a two-way ANOVA (specifically the BA × CI interaction), but the same assumptions apply (Underwood 1991; Smith et al. 1993).
Reviews of designs, new criticisms, and suggestions for design improvement were presented in a series of papers by Underwood (1991, 1992, 1993, 1994). The biggest concerns raised in these papers include the lack of ability to determine causation by the human impact, even with a significant BACI, and the need for multiple Control sites. An asymmetrical design using one Impact site and multiple Control sites was suggested to alleviate the problem of the need for simultaneous sampling, and for interpreting causation so long as the changes observed at the Impact site are greater than the patterns observed in the set of control locations (Underwood 1991, 1994). Stewart-Oaten and Bence (2001) refute Underwood’s criticisms by clarifying that BACI Controls are not experimental controls, and are not chosen randomly but are instead deliberately chosen to be highly correlated with the Impact site to be useful covariates. As such, the variation among Control sites is irrelevant to assessing the problem, and thus multiple Control sites are unnecessary. In addition, the claim that causation cannot be inferred from the BACI model is refuted because of the fundamental reason that the Impact site is not randomly chosen, and therefore the analyses should not be considered that of a typical experiment (Stewart-Oaten and Bence 2001). Their final conclusion is that the BACI design is sound for use in impact assessments so long as statistical analyses are conducted carefully by checking assumptions, exploring sensitivity, inquiring where the “chance” comes from, checking models, and interpreting parameters and results (Stewart-Oaten and Bence 2001).
In Canada, any permanent alteration or destruction of fish habitat is assessed for impacts to affected fisheries to avoid violation of the revised Fisheries Act. Any residual serious harm to fish is addressed via an offsetting plan that includes a monitoring component with the objective of determining whether the proposed measures will be effective in offsetting the serious harm, and that such measures have stabilized and are functioning as intended (DFO 2013). Recent guidelines for effectiveness monitoring designs recommended the BACI design when possible, with a minimum of 3 years of pre-treatment monitoring and a range of options presented for post-term monitoring (Smokorowski et al. 2015). However, recognizing that there are often practical limitations on the duration, frequency, and intensity of monitoring that can be conducted as part of offsetting a development project, a range of suboptimal designs may be adopted, with potential implications for the scientific soundness of the conclusions arising from such monitoring.
With the requirement for effectiveness monitoring now established in regulations, Fisheries and Oceans Canada (DFO) has an opportunity to guide proponent monitoring programs to be scientifically sound, clearly reported, and to contribute towards the long-term improvement of program effectiveness (Smokorowski et al. 2015). Conversely, poorly designed programs may hinder future decision-making ability by presenting, unknowingly, false conclusions about the ecological impacts of development or offsetting actions. With these important consequences in mind, we used data from a long-term BACI experiment as a case study to explore various design and analysis decisions. For illustration purposes, we used 9 years of total invertebrate and fish community biomass and diversity data (3 years pre and 6 years post) from two rivers, one with a hydroelectric facility that changed its operational regime, and one reference river. We used multiple models to explore the results of the BACI analyses using different combinations of years for the “After” period to test the implications of duration of monitoring and balanced vs. unbalanced designs. We used multiple statistical models to explore the effects of model choice on the results, and we examined the resulting confounding factors that can arise in the analyses and interpretation. Specific predictions tested include: (1) the BACI design results would be robust to the statistical model used; (2) the BACI design results would be robust to changes in continuous post-duration monitoring, but that intermittent sampling may affect results; and (3) that the BA or CI design may not produce the same results as the BACI models. Our intention is to show that caution should be employed when using, or interpreting results from, a BACI design in an environmental impact study, but that a well-designed BACI remains one of the best models for environmental effects monitoring programs.
Methods
Case study
To develop the scientific information necessary to design management tools that address both ecosystem integrity and energy efficiency, the DFO, the Ontario Ministry of Natural Resources and Forestry, Brookfield Renewable, and the University of Waterloo collaborated on a long-term, adaptive management experiment to test whether regulating ramping rates through hydroelectric turbines provided ecological benefits, while at the same time minimizing electricity production losses (Smokorowski et al. 2011). The main purpose of this experiment was to determine if removing all operational constraints on ramping rates from a hydroelectric facility was detrimental to the downstream riverine ecology. Aspects of river ecology that were monitored included hydrology, fluvial geomorphology, chemistry and plankton, benthic macroinvertebrates, fish populations and community structure, and food web ecology (via stable isotopes).
We used a BACI design for this experiment, which in this case involved comparing conditions on a regulated peaking river (Impact river; Magpie) to conditions on an unregulated river (Control river; Batchawana), Before and After implementing a change in ramping rates. The Control river had a natural flow regime and was carefully selected to be similar to the Impact system (Metcalfe et al. 2001). The intention was that this approach should allow detection of a change in the metrics of interest that were caused by the experimental ramping rate changes, as the Control river should reflect any changes in metrics resulting from regional environmental factors. Both rivers were similarly sampled during the entire field study, which included a Before period from 2002 to 2004 and an After period from 2005 onward (the planned final year of monitoring is 2017). The purpose of this paper is not to focus on the results of the Magpie project, but instead to use those results to examine the implications of using different BACI analyses methods and different post-impact monitoring designs. Therefore, detailed methods on field sampling will not be provided, nor will a thorough analysis of results from a biological implications perspective.
Study sites
The experimental site was the Magpie River, Wawa, Ontario (48°0′N; 84°7′W) on the 40 km stretch between Steephill Falls and the Harris waterpower facilities (WPF) (Fig. 1). The reference river was the unregulated Batchawana River (47°0′N; 84°3′W), located approximately 60 km north of Sault Ste. Marie, Ontario. Both rivers were divided into transects spaced 500 m apart to provide a spatial reference and allow for the selection of sampling sites in a stratified random sampling design. The first 2.5 km below the dam on the Magpie River was characterized by large deep plunge pools and blasted channels that are atypical of a natural river system. For the Magpie, the most intensively sampled river segments were those considered the most vulnerable to ramping due to their riffle-run-pool riverine nature (total length 18 km, beginning 2.5 km downstream from the dam at Transect 5, Fig. 1) and these are the results presented here. The spatial and sampling coverage was similar for the Batchawana River.
Fig. 1.
Between 2002 and 2004, data were collected from the regulated Magpie River under the original restricted ramping rate regime: ramping rate could not exceed 1 m3·s−1·h−1 from 10 October to 15 November; 2 m3·s−1·h−1 from 16 November until spring freshet (early May); from May until early October, the dam was restricted to an increase or decrease of 25% of the previous hour’s flow. From 2005 to 2010, data were collected with no restrictions on ramping and while the Steephill Falls plant operated in accordance with water availability and market forces (Fig. 2). During the majority of the study period (with one exception noted in the Results), through all seasons the Steephill Falls WPF could not release a discharge lower than 7.5 m3·s−1 as that was the regulated minimum flow. All sampling from the Batchawana River was done contemporaneously.
Fig. 2.
Field sampling
Hydrological data were obtained both from Brookfield Renewable (for data from the Steephill Falls Generating Station on the Magpie River), and from the Water Survey of Canada gauge on the Batchawana River (02BF001). Invertebrate data are from an annual sampling of the community in the thalweg of the river using rock bags. In each year at each site, five mesh rock bags were randomly placed in a riffle, ensuring a minimum distance of 3 m apart, and at a depth to maintain a sufficient flow over the bags throughout low water periods. The rock bags were constructed out of 5.1 cm net mesh, were 122 cm in circumference and 46 cm in length, and were filled with rocks of representative size found along the shoreline at the site of placement until each reached a weight of 7 kg (± 0.5 kg). The actual number of rocks used, their diameter, and the weight of each bag was recorded, as were the water depth and velocity (Marsh McBirney Flomate 2000 Portable Flow Meter) at each bag location. The bags were left in the river for a period of approximately 60 d (June to August), a sufficient length of time for full colonization to reach fluctuating taxa richness, abundance, and biomass (Mason et al. 1973; Shaw and Minshall 1980). Once bags were retrieved the rocks were cleaned and all invertebrates and debris were preserved in 70% ethanol. The entire sample was subsampled for identification to taxonomic level of family, and enumeration, although in each year a number of samples were identified in their entirety to allow for the calculation of accuracy and precision of subsampling procedure, which were always found to be within acceptable limits (defined as being within 20% of true counts, Elliott 1977).
Backpack electrofishing was used to examine relative fish population abundance, biomass, growth, species richness, diversity, and community composition on both rivers. Electrofishing was conducted in July of each year in the Batchawana River and in August of each year in the Magpie River. At each randomly selected transect, all areas ≤60 cm in depth were sampled using backpack electrofishing covering eight 100 m shoreline segments on each river. These segments were sampled according to habitat type (fast or slow) using a standard back-and-fourth electrofishing technique, and shockers were standardized by power (W), rate (s·m−2), area (m2), and time (s). Fish were identified to species and enumerated by habitat area, which were averaged to the transect level. A random subsample of the captured fish was preserved for accurate length, weight and age data. Fish biomass per unit area (g(100 m2)–1) was calculated for each sampling site.
As backpack electrofishing was the exclusive method of fish capture, the fraction of the fish community vulnerable to capture were smaller individuals of larger species, or small species residing in shallower habitats (mean fish length captured was 50 mm ± 26 mm). The rationale was that relative changes observed in this proportion of the fish community would reflect early life history impacts on larger fish species. In addition, as home range size increases allometrically with body size (Minns 1995), sampling smaller body sizes ensured samples were representative of local conditions. Because the home range of fish 150 mm or smaller is <500 m2 (Minns 1995), and the average distance between adjacent transects was 571 m (min−max = 200–5000 m) on rivers with approximate mean widths of >40 m, fishing transects were considered independent replicates.
Statistical analysis
All statistics were calculated using Statistica version 10 (StatSoft Inc., Tulsa, Oklahoma) software. Results were considered significant at p-values that were adjusted using the sequentially rejective Bonferroni test (or Holm–Bonferroni; Holm 1979), corrected for the number of tests within each data set and family of tests. We chose to use the Holm–Bonferroni test to minimize the family-wise error rate for multiple tests, because it is considered more powerful than the overly conservative Bonferroni method which runs the risk of not detecting real effects. We did not consider the full suite of tests in our corrections given that any one study would not normally conduct all analyses presented here, and we felt that correction would be overly conservative. Simply stated, with the Holm–Bonferroni test p-values for your family of tests (n) are ordered from smallest to largest, with the smallest p-value considered relative to α/n (where α = 0.05), the same adjustment as the Bonferroni method. If the first test is found to be not-significant, the procedure stops there; but if it is significant, you compare the second lowest p-value to α/(n−1), and then continue the procedure until the sequentially larger p-value is found to be not significant, at which point all testing stops. Data were tested for parametric assumptions (normality and homogeneity of variance of the data) and transformed when necessary. Residuals from analyses of transformed data were tested for normality. Additivity was tested on transformed data by testing for a zero slope in the regression of the difference and averages of corresponding points between the Impact and Control rivers, which essentially corresponds to Tukey’s test for additivity (Stewart-Oaten et al. 1986). Mean annual flow for each river was calculated as the mean discharge from January to January of each year, and a correlation analysis was run between rivers to test if changes in flow among years were similar.
Fish species and invertebrate families were used to calculate diversity (probability of interspecific encounter (PIE); Hurlburt 1971). PIE is an unbiased diversity measure that calculates the chance that two individuals drawn at random from a population represent different families:where n is the number of all individuals in the sample, ni is the number of individuals of a family in the sample, and s is the number of families (Hurlbert 1971). PIE was selected over other diversity indices because it provides a statistically and biologically understandable probability (out of 100%; the higher the number the more diverse the community), unlike more traditional diversity measures (Gottelli and Graves 1996).
(1)
Mean annual invertebrate community abundance (numbers per rock bag), and mean annual fish biomass-per-unit-effort (by habitat type), and diversity for both taxon groups, were calculated and plotted to illustrate annual variability. To explore the effect of using different statistical models on the results, multiple approaches were used. First, to assess the effect of differing monitoring durations, two-way ANOVAs were run on the BACI data for fish and invertebrate biomass and diversity, using three options for the “after” period, each representing a different combination of continuous 3 year monitoring blocks: (1) 2005–2007 representing the first 3 years post-ramping change and a balanced design (short-term monitoring); (2) 2005–2010 representing a 6 year longer-term monitoring period but creating an unbalanced design, which is unavoidable with a fixed duration pre-period (complete long-term monitoring); and (3) 2008–2010 representing a long-term response period but keeping the design balanced for statistical purposes (balanced long-term monitoring). The main statistic of interest in a BACI analysis is the interaction term (Before-After × Control-Impact), which would be significant when a change occurs at the impact site but not at the control site. To explore the effect of using an alternate statistical model, as recommended by Stewart-Oaten et al. (1986), differences between the Impact minus Control sites were calculated using means by transect and year, and a t test was run on the Before and After periods (including the three continuous block After durations outlined above) to assess if the BACI was significant. To test if the unbalanced design using the 6 year complete long-term monitoring had an effect on the BACI model validity, we compared the standard ANOVA result with the full data set to the same data using a restricted maximum likelihood estimation (REML) ANOVA, which is considered more robust for unbalanced designs (Robinson 1987; Fletcher and Underwood 2002).
In addition, recognizing that there are often pressures to minimize the duration and frequency of post-monitoring requirements, and that there are often practical limitations on monitoring designs, our data were used to test a variety of possible post-impact sampling timetables, using a range of combinations of intermittent sampling in the BACI analyses. Finally, recognizing that sometimes the incorporation of a control system is not feasible, or that the development will occur without the opportunity to collect pre-treatment data, we examined the results from conducting a Before-After only analysis (BA: no Control), and a Control-Impact only analysis (CI: no Before), to examine the conclusions that would have resulted from these suboptimal designs.
Permitting
Fish for this project were collected under the Ontario Ministry of Natural Resources and Forestry Collection of Fish for Scientific Purposes permit number 1000784. Animal care permits for this project were approved by the Canadian Council for Animal Care certified GLLFAS/NWRI Animal Care Committee (permit number GLLFAS/NWRI 0907). Work was reviewed and permits issued on an annual basis.
Results
All abundance data were ln(x + 1) transformed and PIE was arcsine transformed to meet assumptions of additivity, normality, and homogeneity of variance. All transformed data met assumptions except for the invertebrate diversity measure (arcsinePIE), which did not meet the additivity assumption. In addition, residuals from the BACI analyses were not significantly different from normal with the exception of the fish diversity measure (arcsinePIE) in the slow habitat. We present the analyses for invertebrate and slow fish diversity despite these violations, but interpret the results with caution.
The natural hydrology of the Batchawana resulted in greater peak flows and lower minimum flows relative to the altered Magpie River (Fig. 2). In 2002, when ramping rate was restricted, the dam operated on a reduced peaking cycle. During the week when water supply was high, flow was “perched” on an elevated minimum, and on weekends when demand was low water levels did not reach full turbine flow (Fig. 2a). In 2005, however, full ramping from the maximum turbine discharge to minimum regulated flow occurred at a much greater frequency because the rate of change was unrestricted (Fig. 2b). Mean annual flow of the rivers tracked each other over time confirming the suitable use of the Batchawana River as a regional control for the Magpie River (r = 0.78, p = 0.01, Fig. 3). Coincidentally, 2005, which was the year corresponding to the change in unlimited ramping, also saw the beginning of a 3 year period of drought on both rivers. The return to normal annual flow levels was observed on both rivers in 2008 and 2009, but 2010 was again a relatively dry year (Fig. 3).
Fig. 3.
The average number of invertebrates per rock bag was generally greater for the Magpie than the Batchawana, and the pattern among years was roughly similar (Fig. 4a). Invertebrate diversity was initially greater for the Magpie River, but in 2005 the Batchawana diversity became higher than the Magpie until it switched back again in 2009 (Fig. 4b). The BACI analysis interaction was not significant for all year combinations for abundance (Table 1 and Fig. 5a), and was significant for diversity when both short and complete long-term monitoring were tested, but not when the balanced long-term monitoring was tested (Table 1 and Fig. 5b). Calculating the difference between Impact minus Control means by site and running the two sample t test yielded similar results for diversity, but different results for abundance, which was found to be significant against the adjusted p crit values (Table 2). The results for the REML ANOVA for unbalanced designs were the same as the standard BACI factorial ANOVA. As mentioned previously, invertebrate diversity statistics should be interpreted with caution due to the violation of the additivity assumption.
Fig. 4.
Fig. 5.
Table 1.
| After years included | SS | F (df) | p | p crit | |
|---|---|---|---|---|---|
| Invertebrate lnN BACI | |||||
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | 1.8 | 3.8(1,66) | 0.05 | — |
| 1, 2, 3 | Short-term balanced | 1.1 | 4.3(1,48) | 0.04 | — |
| 4, 5, 6 | Long-term balanced | 1.6 | 4.3(1,40) | 0.04 | 0.0125 |
| Invertebrate arcsinePIE BACI | |||||
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | 412.5 | 6.7(1,66) | 0.01 | 0.01 |
| 1, 2, 3 | Short-term balanced | 824.1 | 13.4(1,48) | <0.001 | 0.008 |
| 4, 5, 6 | Long-term balanced | 6.5 | 0.08(1,40) | 0.78 | — |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, and 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance, where in this case the total number of tests was 6. Bold text indicates a significant difference.
Table 2.
| Source of variation—invertebrates | Mean difference Before ± SE (n) | Mean difference After ± SE (n) | t | p | p crit |
|---|---|---|---|---|---|
| lnN (1, 2, 3, 4, 5, 6) | 1.1 ± 0.13 (13) | 0.43 ± 0.17 (22) | 2.7 | 0.01 | 0.0125 |
| lnN (1, 2, 3) | 1.1 ± 0.13 (13) | 0.51 ± 0.16 (13) | 2.8 | <0.01 | 0.01 |
| lnN (4, 5, 6) | 1.1 ± 0.13 (13) | 0.32 ± 0.35 (9) | 2.4 | 0.03 | 0.025 |
| arcsinePIE (1, 2, 3, 4, 5, 6) | 6.5 ± 3.4 (13) | −3.5 ± 2.1 (22) | 2.6 | 0.01 | 0.017 |
| arcsinePIE (1, 2, 3) | 6.5 ± 3.4 (13) | −9.4 ± 1.2 (13) | 4.4 | <0.001 | 0.008 |
| arcsinePIE (4, 5, 6) | 6.5 ± 3.4 (13) | 5.0 ± 3.2 (9) | 0.32 | 0.75 | — |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, and 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance. Bold text indicates a significant difference from before to after the change in dam operations.
Fish biomass per unit effort (BPUE) followed a similar pattern among years on both rivers, although overall interannual variability of the Batchawana was less than the Magpie (Figs. 6a, 6b). For the fish data analyses the total number of tests considered in the Holm–Bonferroni correction was 12, resulting in no significant results regardless of the years included in the design (Table 3 and Figs. 7a, 7b). Calculating the difference between Impact minus Control biomass means by transect and running the two sample t test found the difference between the Before and After period were not significant for either habitat type in the short- or long-term monitoring, but both habitat types were significant when testing the balanced long-term monitoring design (Table 4). Again, the REML ANOVA results were the same as the BACI ANOVA results.
Fig. 6.
Fig. 7.
Table 3.
| After years included | Design type | Habitat | SS | F (df) | p | p crit |
|---|---|---|---|---|---|---|
| Fish lnBPUE BACI | ||||||
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | Fast | 4.2 | 5.5(1,140) | 0.02 | — |
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | Slow | 3.4 | 3.3(1,132) | 0.07 | — |
| 1, 2, 3 | Short-term balanced | Fast | 1.1 | 2.0(1,92) | 0.16 | — |
| 1, 2, 3 | Short-term balanced | Slow | 0.34 | 0.45(1,88) | 0.50 | — |
| 4, 5, 6 | Long-term balanced | Fast | 6.1 | 7.9(1,92) | 0.006 | — |
| 4, 5, 6 | Long-term balanced | Slow | 7.3 | 7.8(1,84) | 0.006 | 0.004 |
| Fish arcsinePIE BACI | ||||||
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | Fast | 123 | 0.95(1,140) | 0.33 | — |
| 1, 2, 3, 4, 5, 6 | Long-term unbalanced | Slow | 311 | 1.4(1,132) | 0.23 | — |
| 1, 2, 3 | Short-term balanced | Fast | 14.9 | 0.15(1,90) | 0.69 | — |
| 1, 2, 3 | Short-term balanced | Slow | 95 | 0.57(1,88) | 0.45 | — |
| 4, 5, 6 | Long-term balanced | Fast | 237 | 1.7(1,92) | 0.19 | — |
| 4, 5, 6 | Long-term balanced | Slow | 453.8 | 1.7(1,84) | 0.19 | — |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, and 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance, where in this case the total number of tests was 12.
Table 4.
| Source of variation—fish lnBPUE or diversity (arcsinePIE) | Habitat | Mean difference Before ±SE (n) | Mean difference After ±SE (n) | t | p | p crit |
|---|---|---|---|---|---|---|
| lnBPUE (1, 2, 3, 4, 5, 6) | Fast | 0.84 ± 0.20 (24) | 0.12 ± 0.16 (46) | 2.6 | 0.01 | 0.005 |
| Slow | 0.25 ± 0.22 (22) | −0.43 ± 0.20 (46) | 2.0 | 0.05 | — | |
| lnBPUE (1, 2, 3) | Fast | 0.84 ± 0.20 (24) | 0.40 ± 0.19 (24) | 1.6 | 0.12 | — |
| Slow | 0.25 ± 0.22 (22) | 0.005 ± 0.23 (24) | 0.69 | 0.49 | — | |
| lnBPUE (4, 5, 6) | Fast | 0.84 ± 0.20 (24) | −0.17 ± 0.26 (24) | 3.1 | 0.003 | 0.0045 |
| Slow | 0.25 ± 0.23 (22) | −0.90 ± 0.28 (22) | 3.2 | 0.003 | 0.004 | |
| arcsinePIE (1, 2, 3, 4, 5, 6) | Fast | −8.5 ± 2.6 (24) | −12.5 ± 2.6 (48) | 0.96 | 0.34 | — |
| Slow | −5.6 ± 4.0 (22) | −12.1 ± 3.1 (46) | 1.2 | 0.22 | — | |
| arcsinePIE (1, 2, 3) | Fast | −8.5 ± 2.6 (24) | −10.1 ± 3.1 (24) | 0.39 | 0.70 | — |
| Slow | −5.6 ± 4.0 (22) | −9.7 ± 2.5 (24) | 0.88 | 0.38 | — | |
| arcsinePIE (4, 5, 6) | Fast | −8.5 ± 2.6 (24) | −14.8 ± 4.1 (24) | 1.3 | 0.20 | — |
| Slow | −5.6 ± 4.0 (22) | −14.7 ± 5.8 (22) | 1.3 | 0.20 | — |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, and 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance. Bold text indicates a significant difference from before to after the change in dam operations.
Fish diversity in both habitat types was greater for the Batchawana than for the Magpie River (Figs. 6c, 6d). The interannual pattern of fish diversity in fast habitat was similar among years between the two rivers except in 2009. Fish diversity in the slow habitat was relatively consistent among years for the Batchawana River, whereas the interannual variability of the Magpie was great. The interaction term for fish diversity was not significant for any combination of years in either habitat type (Table 3 and Figs. 7c, 7d). The t test results similarly demonstrated that the difference in fish diversity from Before to After was not significant for any year combination (Table 4).
Exploring various intermittent sampling designs demonstrated that some significant impacts were no longer apparent, and very different conclusions could result depending on the years included in the monitoring program (Tables 5 and 6). The use of Before-After Magpie River only design (no Control) was not able to detect a significant change in either the fish or invertebrate community in the impact river, although in most cases there was a decreasing trend in the metric from Before to After the change to unlimited ramping (Table 7). The Control-Impact design (no Before) was similarly not able to detect a difference in abundance between the Batchawana and Magpie Rivers, except that it did capture that both invertebrate and fish diversity was lower in the impact Magpie River relative to the control Batchawana River (Table 7).
Table 5.
| After years included | SS | F (df) | p | p crit | |
|---|---|---|---|---|---|
| Invertebrate lnN BACI | |||||
| 1, 3, 6 | BA × CI | 0.91 | 1.9(1,46) | 0.17 | — |
| 2, 4, 6 | BA × CI | 0.82 | 3.0(1,42) | 0.09 | — |
| 1, 3, 5 | BA × CI | 1.8 | 7.0(1,46) | 0.01 | 0.01 |
| Invertebrate arcsinePIE BACI | |||||
| 1, 3, 6 | BA × CI | 351.9 | 4.4(1,46) | 0.04 | 0.0125 |
| 2, 4, 6 | BA × CI | 104.3 | 1.3(1,42) | 0.26 | — |
| 1, 3, 5 | BA × CI | 556.8 | 7.9(1,46) | 0.007 | 0.008 |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, and 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance. Bold text indicates a significant difference.
Table 6.
| After years included | Habitat | SS | F (df) | p | p crit |
|---|---|---|---|---|---|
| Fish lnBPUE BACI | |||||
| 1, 3, 6 | Fast | 1.7 | 3.2(1,92) | 0.08 | — |
| 1, 3, 6 | Slow | 1.3 | 1.7(1,86) | 0.20 | — |
| 2, 4, 6 | Fast | 2.4 | 3.1(1,92) | 0.08 | — |
| 2, 4, 6 | Slow | 4.0 | 3.6(1,86) | 0.06 | — |
| 1, 3, 5 | Fast | 3.9 | 5.6(1,92) | 0.02 | 0.004 |
| 1, 3, 5 | Slow | 1.5 | 1.9(1,86) | 0.17 | — |
| Fish arcsinePIE BACI | |||||
| 1, 3, 6 | Fast | 25.6 | 0.28(1,92) | 0.60 | — |
| 1, 3, 6 | Slow | 78.7 | 0.38(1,86) | 0.54 | — |
| 2, 4, 6 | Fast | 21.6 | 0.16(1,92) | 0.69 | — |
| 2, 4, 6 | Slow | 190.2 | 0.94(1,86) | 0.34 | — |
| 1, 3, 5 | Fast | 213.5 | 2.1(1,92) | 0.15 | — |
| 1, 3, 5 | Slow | 285.3 | 1.2(1,86) | 0.27 | — |
Note: In all cases the Before period included 2002–2004. Years included in the After period varied as indicated, where 1 = 2005, 2 = 2006, 3 = 2007, 4 = 2008, 5 = 2009, 6 = 2010. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance.
Table 7.
| Test | Habitat | lnN/Fish BPUE F (df), p, (trend) | arcsinePIE F (df), p, (trend) | p crit |
|---|---|---|---|---|
| Invertebrates Before-After (no Control) | 2.6(1,33), 0.12, (↓) | 2.6(1,33), 0.11, (↓) | — | |
| Invertebrates Control-Impact (no Before) | 3.2(1,42), 0.08, (↑) | 4.0(1,42), 0.05, (↓) | 0.0125 | |
| Fish Before-After (no Control) | Fast | 1.9(1,70), 0.18, (↓) | 0.32(1,70), 0.57, (↓) | — |
| Slow | 0.87(1,66), 0.35, (↓) | 0.004(1,66), 0.95, (↔) | — | |
| Fish Control-Impact (no Before) | Fast | 0.41(1,94), 0.52, (↑) | 24.0(1,94), <0.001, (↓) | 0.00625 |
| Slow | 3.5(1,90), 0.06, (↓) | 15.6(1,90), <0.001, (↓) | 0.007 |
Note: The trend arrow indicates the directional difference (increasing or decreasing) between the mean values from Before to After or between Control and Impact. p crit provides the critical p-value as adjusted by the Holm–Bonferroni method for determining significance, where in this case the total number of tests was 4 for invertebrates and 8 for fish.
Discussion
The purpose of this paper is to explore how results and conclusions can change depending on the design and analysis of the monitoring program, and to raise awareness of resulting implications. While it seems clear that the change in ramping rate regime had an effect on the biota of the Magpie River (details of the biotic response are discussed elsewhere), the interpretation of the data and conclusions of the statistical significance of those effects varied depending on the number and configuration of years included in the analysis, and the statistical model used to analyze the data. The computational procedure used to analyze the unbalanced BACI design resulting from continuous long-term monitoring with a fixed baseline phase did not have an effect on the interpretation of the data. However, it did appear that analyzing the difference between the reference and control before and after the impact occurred (using a t test as recommended by Stewart-Oaten et al. 1986), was a more powerful statistical test than the BACI ANOVA, as significant differences were declared despite the higher bar set by the Holm–Bonferroni adjusted p crit values. Intermittent sampling designs often failed to detect effects that were found with continuous sampling, as did the suboptimal designs used when no Control system or Before data were available. We are unaware of any similar exploration of design and analysis options on environmental effects monitoring results and conclusions. Although others have examined optimal study designs for environmental effects or restoration programs (e.g., Benedetti-Cecchi 2001; Munkittrick et al. 2002, 2009; Liermann and Roni 2008), few have had the duration and consistency of data necessary for such analyses.
Using a balanced number of years before and after affecting a perturbation is usually considered important in statistical analyses due to its effect on the power of the test, the implications regarding the assumption of equality in variance, and in how you estimate parameters and partition the sum of squares (SS) to test the hypotheses of interest (Shaw and Mitchell-Olds 1993). Yet due to a myriad of factors, ecologists rarely have completely balanced data and need to face decisions about imposing balance by eliminating values (thus losing information), or to use computational methods that are specifically designed for unbalanced data (Robinson 1987; Shaw and Mitchell-Olds 1993). For any long-term environmental effects monitoring study, the period of pre-treatment or “before” data is often limited, and therefore the longer duration of monitoring results in greater imbalance in the design. Yet long-term monitoring programs, although relatively rare especially in environmental effects assessments, are considered important to capture the full manifestation of management actions as system response times can be on the order of decades (Burt et al. 2008; Roni et al. 2008; Lindenmayer and Likens 2009). Therefore, reconciling the trade-off between long-term monitoring and balanced statistical analyses becomes necessary.
In our example, the original balanced design was termed “short-term” monitoring because it included only 3 years post-impact monitoring. While not having any statistical concerns due to imbalanced design, if sampling ended in 2007 we would have concluded that the unlimited ramping had no significant effect on relative fish biomass or diversity, but did have a significant effect on invertebrate abundance (using the difference tests) and diversity. After 6 years of post-impact monitoring differences in fish biomass were greater but still not significant using the adjusted p crit values. Using a computational procedure specifically designed for unbalanced designs (REML ANOVA, Robinson 1987) did not alter the significance of our BACI interaction term, or the interpretation of our data. If our main concern was achieving statistical balance, this was easily achieved by eliminating data and analyzing 3 year packages, which we refer to as balanced long-term monitoring. Losing years of data would result in detecting a statistically significant difference in fish biomass results, but eliminating significance in invertebrate diversity, and would have altered our conclusions. Therefore, achieving statistical balance by eliminating data should not be considered an option in environmental effects monitoring, preferentially using a computational procedure specific to unbalanced designs (available in most statistical packages) to ensure that statistical bias is not introduced by long-term environmental effects monitoring.
Smokorowski et al. (2015) provided recommendations for effectiveness monitoring designs that included: (1) 3 years before development at the impact site to establish a baseline, (2) 3 years sampling immediately post-change, (3) an additional 3 years of sampling at a later time (e.g., 4–6 years post or some later time), and (4) a revisit 10 years after project impact to capture longer-term changes to the site. These 3 year blocks of time were recommended to allow for the quantification of interannual variability, increasing the likelihood of being able to distinguish project effects from natural variability. However, when negotiating the duration and frequency of environmental effects monitoring programs, resource managers are often faced with proposals suggesting alternate years or multiple skipped years to increase the overall duration while minimizing the costs. For example, it may be suggested that years 1, 3, and 6 post-impact (or post-offsetting) should be monitored to capture a longer duration response but still minimize effort, but in this case no differences would have been found. Various scenarios could be envisioned and a number were tested with our data, which clearly demonstrated that the conclusions from intermittent designs would have been different, and would sometimes have failed to detect an effect that was found using continuous sampling. It is not unusual for long-term datasets to vary in design, with intermittent sampling as a common feature, which can have significant impacts on the ecological measure and change detection ability (Magurran et al. 2010).
When annual data are available, it is important that statistical model results should not be the sole basis for reaching conclusions of impact. A visual examination of the annual trend plots clearly shows why selecting different combinations of post-monitoring years could affect the significance of the statistical result. For example, the flip-flop between the Control and Impact sites for invertebrate diversity emphasize that caution is needed when interpreting the “significant impact” that was found in the long-term continuous design. In 2010, the invertebrate diversity at the Impact site was greater than at the Control site, despite continued unlimited ramping, indicating that invertebrate diversity is likely influenced by factors other than ramping rates, and the statistical significance was still being driven by 2005–2008 data. The violation of the additivity assumption for invertebrate diversity similarly calls any significant impact into question (Stewart-Oaten et al. 1986, 1992). The annual plots also highlight the potential pitfalls of intermittent sampling. For example, monitoring invertebrate abundance during 2006, 2008, and 2010 would have hit high abundance years in the Magpie River and concluded no impact; monitoring 2005, 2007 and 2009 would have hit low abundance years in the Magpie River and would have concluded the ramping change had a significant effect on invertebrate abundance. Thus, long-term monitoring programs that consider 1 year snapshots of data on an intermittent basis may have difficulty arriving at any conclusions, whether statistical or observational.
Invertebrate abundance and diversity were highly variable over the 10 years studied, and no one consistent trend was apparent in either abundance or diversity for the statistical models or range of designs tested. A quantitative review of the literature demonstrated mixed responses of macroinvertebrate abundance to changes in flow magnitude, similarly preventing the development of any robust statistical relationships (Poff and Zimmerman 2010). The same review demonstrated a consistent negative response of fishes to alteration in flow magnitude (Poff and Zimmerman 2010). So perhaps our inconsistent invertebrate statistics are as much a reflection of the true (and inconsistent) invertebrate community response, and that this rapid-responding, ephemeral ecological group is not the best to use to as an indicator of flow alteration in environmental assessments.
The fish community appears to be a more consistent and potentially reliable indicator of environmental effects. Backpack shocking sampled the smaller size classes and species of fish in the wadable portion of the river. Thus, if response to environmental change at the population level is linked to life history duration, sampling these smaller species suggests that change, if it was to occur, would manifest in 3–5 years, or the maximum age reached by the majority of fishes sampled. Fish community biomass generally did not fluctuate greatly year-to-year, and the patterns between the Magpie and Batchawana were similar, with the greatest fluctuations observed in the Magpie close to the dam. This similarity reinforces the Batchawana as a good covariate for the Magpie in this BACI experiment (Stewart-Oaten and Bence 2001). Our BACI statistics indicate that after 3 years a negative impact on fish biomass was detectable; a time period that corresponds to when the majority of species will have completed one generation. The intermittent sampling analysis, however, failed to detect this effect. Fish diversity was never found to be significantly altered by the change to unlimited ramping no matter what combination of years were included in the analyses, but the Control-Impact analysis correctly detected that fish diversity was consistently greater in the Batchawana than the Magpie. If no Before data were available, the implication would have been that the lower diversity was caused by the unlimited ramping regime, highlighting the critical importance of Before data in helping to tease out causation.
Environmental impact assessments conducted in the real world are often limited by temporal and financial constraints, resulting in lack of adequate before data or multiple controls. This paper demonstrates that employing a long-term BACI design remains one of the best ways of detecting an effect from an environmental change, but as cautioned by Stewart-Oaten and Bence (2001), this is true only if statistical analyses are conducted carefully. The use of the BACI design with a minimum of 3 years of before data was suggested as the best scientific framework for assessing the effectiveness of offsetting activities required under the Fisheries Act (2012) in Canada (Smokorowski et al. 2015), and this paper further supports that recommendation. In addition, blocks of continuous monitoring after the impact provide more robust conclusions, and it is to the benefit of the project proponent, the regulator, and the aquatic system that poor designs and misleading statistical tests do not lead to flawed decisions. As far as the ramping rate experiment goes, exploring the mechanisms behind the observed biotic responses, if they continue with the inclusion of additional years of monitoring, requires an examination of environmental characteristics such as the flow and temperature regimes, and this is the subject of ongoing analyses. Also still to be considered is if the statistical significance translates into biological relevance for the impacted system (Arciszewski and Munkittrick 2015), but that question is beyond the scope of this paper.
Acknowledgements
We thank Fisheries and Oceans Canada, Ontario Ministry of Natural Resources and Forestry, and Brookfield Renewable for supporting this long-term monitoring program. We also thank Evan Timusk for assisting with data compilation and validation.
References
Arciszewski TJ, and Munkittrick KR. 2015. Development of an adaptive monitoring framework for long-term programs: an example using indicators of fish health. Integrated Environmental Assessment and Management, 11(4): 701–718.
Benedetti-Cecchi L. 2001. Beyond BACI: optimization of environmental sampling designs through monitoring and simulation. Ecological Applications, 11: 783–799.
Bowman FB, and Somers KM. 2005. Considerations when using the reference condition approach for bioassessment of freshwater ecosystems. Water Quality Research Journal of Canada, 40: 347–360.
Burt TP, Howden NJK, Worrall F, and Whelan MJ. 2008. Importance of long-term monitoring for detecting environmental change: lessons from a lowland river in south east England. Biogeosciences, 5: 1529–1535.
Chapman MG. 1998. Improving sampling designs for measuring restoration in aquatic habitats. Journal of Aquatic Ecosystem Stress and Recovery, 6: 235–251.
DFO. 2013. Fisheries productivity investment policy: a proponent’s guide to offsetting. Catalogue Number: Fs23-596/2013E-PDF. [online]: Available from http://www.dfo-mpo.gc.ca/pnw-ppe/offsetting-guide-compensation/offsetting-guide-compensation-eng.pdf.
Elliott JM. 1977. Some methods for the statistical analysis of samples of benthic invertebrates. 2nd edition. Freshwater Biological Association Scientific Publication No. 25, Ambleside, UK. 159 p.
Fletcher DJ, and Underwood AJ. 2002. How to cope with negative estimates of components of variance in ecological field studies. Journal of Experimental Marine Biology and Ecology, 273: 89–95.
Gottelli NJ, and Graves GR. 1996. Null models in ecology. Smithsonian Institution Press, Herndon, Virginia. 368 p.
Green RH. 1979. Sampling design and statistical methods for environmental biologist. Wiley Interscience, Chichester, UK.
Holm S. 1979. A simple sequentially rejective multiple test procedure. Scandinavian Journal of Statistics, 6(2): 65–70.
Hurlbert SH. 1971. The nonconcept of species diversity: a critique and alternative parameters. Ecology, 52: 577–586.
Hurlbert SH. 1984. Pseudoreplication and the design of ecological field experiments. Ecological Monographs, 54: 187–211.
Liermann M, and Roni P. 2008. More sites or more years? Optimal study design for monitoring fish response to watershed restoration. North American Journal of Fisheries Management, 28(3): 935–943.
Lindenmayer DB, and Likens GE. 2009. Adaptive monitoring: a new paradigm for long-term research and monitoring. Trends in Ecology & Evolution, 24(9): 482–486.
Magurran AE, Baillie SR, Buckland ST, Dick JMCP, Elston DA, Scott EM, et al. 2010. Long-term datasets in biodiversity research and monitoring: assessing change in ecological communities through time. Trends in Ecology & Evolution, 25(10): 574–582.
Mason WT, Weber CI, Lewis PA, and Julian EG. 1973. Factors affecting the performance of basket and multiplate macroinvertebrate samplers. Freshwater Biology, 3: 409–436.
Metcalfe RA, House DA, and Jahncke R. 2001. Waterpower project—physical examination of possible reference study sites for the Magpie river ramping study, Watershed Science Centre Report: WSC.01.1, Ontario Ministry of Natural Resources; 18 p.
Minns CK. 1995. Allometry of home range sizes in lake and river fishes. Canadian Journal of Fisheries and Aquatic Sciences, 52: 1499–1508.
Munkittrick KR, McGeachy A, McMaster ME, and Courtenay SC. 2002. Overview of freshwater fish studies from the pulp and paper environmental effects monitoring program. Water Quality Research Journal of Canada, 37(1): 49–77.
Munkittrick KR, Arens CJ, Lowell RB, and Kaminski GP. 2009. A review of potential methods of determining critical effect size for designing environmental monitoring programs. Environmental Toxicology and Chemistry, 28(7): 1361–1371.
Poff NL, and Zimmerman JKH. 2010. Ecological responses to altered flow regimes: a literature review to inform the science and management of environmental flows. Freshwater Biology, 55: 194–205.
Robinson DL. 1987. Estimation and use of variance components. The Statistician, 36: 3–14.
Roni P, Hanson K, and Beechie T. 2008. Global review of the physical and biological effectiveness of stream habitat rehabilitation techniques. North American Journal of Fisheries Management, 28: 836–890.
Shaw DW, and Minshall GW. 1980. Colonization of an introduced substrate by stream macroinvertebrates. Oikos, 34: 259–271.
Shaw RG, and Mitchell-Olds T. 1993. ANOVA for unbalanced data: an overview. Ecology, 74(6): 1638–1645.
Smith EP, Orvos BW, and Cairns J Jr. 1993. Impact assessment using the Before-After-Control-Impact (BACI) model: concerns and comments. Canadian Journal of Fisheries and Aquatic Sciences, 50: 627–637.
Smokorowski KE, Bradford MJ, Clarke KD, Clément M, Gregory RS, and Randall RG. 2015. Assessing the effectiveness of habitat offset activities in Canada: monitoring design and metrics. Canadian Technical Report of Fisheries and Aquatic Sciences: 3132: vi + 48 p.
Smokorowski KE, Metcalfe RA, Finucan SD, Jones N, Marty J, Power, et al. 2011. Ecosystem level assessment of environmentally-based flow restrictions for maintaining ecosystem integrity: a comparison of a modified peaking versus unaltered river. Ecohydrology, 4: 791–806.
Stewart-Oaten A, and Bence JR. 2001. Temporal and spatial variation in environmental impact assessment. Ecological Monographs, 71: 305–339.
Stewart-Oaten A, Bence JR, and Osenberg CW. 1992. Assessing effects of unreplicated perturbations: no simple solutions. Ecology, 73(4): 1396–1404.
Stewart-Oaten A, Murdoch WW, and Parker KR. 1986. Environmental impact assessment: “Pseudoreplication” in time? Ecology, 67: 929–940.
Underwood AJ. 1991. Beyond BACI: experimental designs for detecting human environmental impacts on temporal variations in natural populations. Australian Journal of Marine and Freshwater Research, 42: 569–587.
Underwood AJ. 1992. Beyond BACI: the detection of environmental impacts on populations in the real, but variable, world. Journal of Experimental Marine Biology and Ecology, 161: 145–178.
Underwood AJ. 1993. The mechanics of spatially replicated sampling programmes to detect environmental impacts in a variable world. Australian Journal of Ecology, 18: 99–116.
Underwood AJ. 1994. On beyond BACI: sampling designs that might reliably detect environmental disturbances. Ecological Applications, 4(1): 3–15.
Information & Authors
Information
Published In
FACETS
Volume 2 • Number 1 • January 2017
Pages: 212 - 232
Editor: Brett Favaro
History
Received: 30 September 2016
Accepted: 5 December 2016
Version of record online: 2 March 2017
Copyright
© 2017 Her Majesty the Queen in right of Canada. This work is licensed under a Creative Commons Attribution 4.0 International License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Data Availability Statement
All relevant data are within the paper.
Key Words
Sections
Subjects
Authors
Author Contributions
All conceived and designed the study.
KES performed the experiments/collected the data.
All analyzed and interpreted the data.
KES contributed resources.
All drafted or revised the manuscript.
Competing Interests
The authors have declared that no competing interests exist.
Metrics & Citations
Metrics
Other Metrics
Citations
Cite As
K.E. Smokorowski and R.G. Randall. 2017. Cautions on using the Before-After-Control-Impact design in environmental effects monitoring programs. FACETS.
2(): 212-232. https://doi.org/10.1139/facets-2016-0058
Export Citations
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.
Cited by
1. Can linear transportation infrastructure verges constitute a habitat and/or a corridor for vascular plants in temperate ecosystems? A systematic review
2. Evaluating ecological restoration in urban ecosystems with acoustic telemetry: marine and freshwater case studies
3. Evaluating the efficacy of ecological restoration of fish habitat in coastal waters of Lake Ontario
4. Reduction of large vessel traffic improves water quality and alters fish habitat-use throughout a large river
5. Long-term outcomes of passive and active restoration approaches following a vessel grounding in Hawaii, USA
6. A multi‐site and hypothesis‐driven approach to identify controls on the bedload transport regime of an anthropised gravel‐bed river
7. Ground-Based Pyrethroid Adulticides Reduce Mosquitoes But Not Nontarget Insects in Central Florida
8. Identifying anthropogenic legacy in freshwater ecosystems
9. Temporal change and impact on air quality of an energy recovery plant using the M-BACI design in Gipuzkoa
10. Spatial and temporal patterns of diffuse attenuation coefficient in Guánica, Puerto Rico: 15 years after development of a watershed management plan
11. Cumulative Effects of Watershed Disturbances and Run-of-river Dams on Mercury Cycling: Case Study and Recommendations for Environmental Managers
12. Comment: A reexamination of , bed-scale impact and recovery of a commercially important intertidal seaweed
13. Assessing wild turkey productivity before and after a 14‐day delay in the start date of the spring hunting season in Tennessee
14. Towards the development of a standardized method to monitor lakeshore hydromorphological restoration measures
15. From DNA to diagnostics: A case study using macroinvertebrate metabarcoding to assess the effectiveness of restoration measures in a Dutch stream
16. Swimmer’s itch control: Timely waterfowl brood relocation significantly reduces an avian schistosome population and human cases on recreational lakes
17. Gravel washing as a lacustrine spawning habitat restoration method for smallmouth bass
18. The effects of gas extraction under intertidal mudflats on sediment and macrozoobenthic communities
19. Grazing capacity change on a small Mediterranean-type South African offshore island following the control of invasive alien European rabbits (
Oryctolagus cuniculus
)
20. Before–after control–impact (BACI) assessment of the effects of a deep-water no-take fishery reserve to recover Norway lobster (
Nephrops norvegicus
) overfished populations and coexisting megafauna
21. Long‐Term
Standardized Sampling as a Tool for Evaluating Management Actions:
BACI
Analysis to Reevaluate Predatory Control of Stunted Crappie with Saugeye
22. Temporal development of fish communities in brackish lagoons of the Baltic Sea considering the invasion of the non-indigenous round goby (Neogobius melanostomus)
23. Towards a better future for biodiversity and people: Modelling Nature Futures
24. Understanding ‘it depends’ in ecology: a guide to hypothesising, visualising and interpreting statistical interactions
25. Immediate impact of fires on roadkilling of wild vertebrates on a highway in southeast Brazil
26. Moderate recreational oyster harvest has variable impacts on biophysical properties of reefs and reef‐associated fauna
27. Controlling common carp (Cyprinus carpio): barriers, biological traits, and selective fragmentation
28. A need to better monitor the effects of coastal defence measures on coastal socio-ecological systems to improve future adaptation solutions
29. Hydrologic restoration of anthropoentically altered springs in the Sheldon National Wildlife Refuge in the Great Basin, USA
30. Potential protection effects in a small marine reserve: a demersal community analysis from stereo baited remote underwater video
31. Attaining freshwater and estuarine-water soil saturation in an ecosystem-scale coastal flooding experiment
32. Effects of Wind Power Development on Reindeer: Global Positioning System Monitoring and Herders’ Experience
33. Grassland restoration benefits declining farmland birds: A landscape-scale before-after-control-impact experiment
34. Spawning Atlantic cod (
Gadus morhua
L.) exposed to noise from seismic airguns do not abandon their spawning site
35. Accelerated marsh erosion following the Deepwater Horizon oil spill confirmed, ameliorated by planting
36. How do changes in flow magnitude due to hydropower operations affect fish abundance and biomass in temperate regions? A systematic review
37. Evaluating the impact of adaptation interventions on vulnerability and livelihood resilience
38. The direct and indirect effects of damming on the Hippopotamus amphibius population abundance and distribution at Bui National Park, Ghana
39. Whole-lake nanosilver additions reduce northern pike (Esox lucius) growth
40. Can we incorrectly link armouring to damming? A need to promote hypothesis-driven rather than expert-based approaches in fluvial geomorphology
41. Small Hydropower—Small Ecological Footprint? A Multi-Annual Environmental Impact Analysis Using Aquatic Macroinvertebrates as Bioindicators. Part 1: Effects on Community Structure
42. Climate change mitigation and improvement of water quality from the restoration of a subtropical coastal wetland
43. Looking for a Simple Assessment Tool for a Complex Task: Short-Term Evaluation of Changes in Fisheries Management Measures in the Pomo/Jabuka Pits Area (Central Adriatic Sea)
44. The effects of population management on wild ungulates: A systematic map of evidence for UK species
45. Do headwater lakes moderate downstream temperature response to forest harvesting? Illustrating opportunities and obstacles associated with virtual experiments
46. Impacts of coastal realignment on biodiversity. A systematic review and meta-analysis
47. Assessing the short-term response of fish assemblages to damming of an Amazonian river
48. Changes in carabid diversity indicate successful restoration of riparian habitats
49. Are natural floods accelerators for streambank vegetation development in floodplain restoration?
50. Design considerations for rapid biodiversity reconnaissance surveys and long‐term monitoring to assess the impact of wildfire
51. Can reductions in water residence time be used to disrupt seasonal stratification and control internal loading in a eutrophic monomictic lake?
52. Lobster reserves as a management tool in coastal waters: Two decades of experience in Norway
53. A spatial causal network approach for multi-stressor risk analysis and mapping for environmental impact assessments
54. What’s in My Toolkit? A Review of Technologies for Assessing Changes in Habitats Caused by Marine Energy Development
55. The effects of endangered freshwater pearl mussels on channel morphology and flow in a low-gradient sandy river
56. Spatio-temporal impacts of aerial adulticide applications on populations of West Nile virus vector mosquitoes
57. Bird Displacement by Wind Turbines: Assessing Current Knowledge and Recommendations for Future Studies
58. Editorial: Before-After Control-Impact (BACI) Studies in the Ocean
59. Assessing the effectiveness of global protected areas based on the difference in differences model
60. Effects of Livestock Exclusion on Stream Banks and Riparian Vegetation in Washington and Oregon
61. Recent Invaders in Small Mediterranean Islands: Wild Boars Impact Snakes in Port-Cros National Park
62. Distance-Based Sampling Methods for Assessing the Ecological Effects of Offshore Wind Farms: Synthesis and Application to Fisheries Resource Studies
63. Assessment of coastal anthropo-ecological system dynamics in response to a tsunami catastrophe of an unprecedented magnitude encountered in Japan
64. Differential effects of a catastrophic wildfire on downstream fish assemblages in an aridland river
65. A Critical Review of the IUCN UK Peatland Programme’s “Burning and Peatlands” Position Statement
66. Ecological impacts of the LED-streetlight retrofit on insectivorous bats in Singapore
67. Nutrient enrichment effects are conditional on upstream nutrient concentrations: Implications for bioassessment in multi-use catchments
68. Differential reporting of biodiversity in two citizen science platforms during COVID-19 lockdown in Colombia
69. Wood Removals from Lakes May Not Necessarily Elicit Fish Population Responses
70. Early Ecosystem Development Varies With Elevation and Pre-Restoration Land Use/Land Cover in a Pacific Northwest Tidal Wetland Restoration Project
71. Quantifying and addressing the prevalence and bias of study designs in the environmental and social sciences
72. Evidence of the impact of noise pollution on biodiversity: a systematic map
73. Collateral benefits of targeted supplementary feeding on demography and growth rate of a threatened population
74. Sharing the neighbourhood: assessing the impact of kelp harvest on foraging behaviour of the European shag
75. Trends in water quality in a subtropical Australian river-estuary system: Responses to damming, climate variability and wastewater discharges
76. Using Assessment, Inventory, and Monitoring data for evaluating rangeland treatment effects in Northern New Mexico
77. Monitoring fisheries resources at offshore wind farms: BACI vs. BAG designs
78. Effects of crowding due to habitat loss on species assemblage patterns
79. Assessing macroinvertebrate community response to restoration of Big Spring Run: Expanded analysis of before‐after‐control‐impact sampling designs
80. The effectiveness of spawning habitat creation or enhancement for substrate-spawning temperate fish: a systematic review
81. Do static and dynamic marine protected areas that restrict pelagic fishing achieve ecological objectives?
82. Simple study designs in ecology produce inaccurate estimates of biodiversity responses
83. Research questions to facilitate the future development of European long-term ecosystem research infrastructures: A horizon scanning exercise
84. Estimated Sediment Protection Efficiences for Increasing Levels of Best Management Practices on Forest Harvests in the Piedmont, USA
85. Timber harvest alters mercury bioaccumulation and food web structure in headwater streams
86. Endangered species protection and evidence-based decision-making: Case study of a quarry proposal in endangered turtle habitat
87. Multi-Decadal Change in Reef-Scale Production and Calcification Associated With Recent Disturbances on a Lizard Island Reef Flat
88. The Vulnerability of Microbial Ecosystems in a Changing Climate: Potential Impact in Shark Bay
89. Population‐level consequences of seismic surveys on fishes: An interdisciplinary challenge
90. The Effect of Hydrological Connectivity on Fish Assemblages in a Floodplain System From the South-East Gulf of California, Mexico
91. Impacts of freshwater aquaculture on fish communities: A whole‐ecosystem experimental approach
92. Rangifer within areas of human influence: understanding effects in relation to spatiotemporal scales
93. Spatial distribution shifts in two temperate fish species associated to a newly-introduced tropical seaweed invasion
94. A Review of Approaches for Monitoring the Effectiveness of Regional River Habitat Restoration Programs
95. Does restoration work? It depends on how we measure success
96. An integrated approach to assessing marine seismic impacts: Lessons learnt from the Gippsland Marine Environmental Monitoring project
97. Salty stories, fresh spaces: Lessons for aquatic protected areas from marine and freshwater experiences
98. Evidence‐based restoration in the Anthropocene—from acting with purpose to acting for impact
99. Wildfire affects space use and movement of white-tailed deer in a tropical pyric landscape
100. Large-Scale Removal of Invasive Honeysuckle Decreases Mosquito and Avian Host Abundance
101. Increased size and relative abundance of migratory fishes observed near the Athabasca oil sands
102. How experimental biology and ecology can support evidence-based decision-making in conservation: avoiding pitfalls and enabling application