Data and code from: Invertebrate responses one year after beaver reintroduction
Data files
Nov 03, 2026 version files 1.20 MB
-
alldata.csv
8.47 KB
-
allraw.csv
4.44 KB
-
CompleteBeaverAnalysis.R
66.58 KB
-
InvertebrateDataManipulation.Rmd
32.23 KB
-
inverts.csv
506.50 KB
-
panTrap.csv
274.19 KB
-
panTrapalldata.csv
2.84 KB
-
pitfall.csv
43.18 KB
-
pitfalldata.csv
2.80 KB
-
README.md
6.56 KB
-
SurberSamples.csv
248.25 KB
Abstract
Reintroductions of Eurasian beavers (Castor fiber) are gathering pace across the UK after a 400-year local extinction. Amid substantial global invertebrate biodiversity declines, the reintroduction of beavers poses opportunities for gains both for invertebrate and beaver populations, since beavers are “ecosystem engineers”: they create substantial changes in hydrology and therefore instream and riparian habitats. Previous studies on beavers (both C. canadensis in N. America and C. fiber in Europe) have primarily used space-for-time substitutions, due to the lack of opportunity to examine a controlled before/after approach. Other studies have focused on North America or mainland Europe and have not typically included both freshwater and terrestrial invertebrates, or quantified changes in invertebrate communities using biomass. To fill these knowledge gaps, we examined invertebrate communities, algal concentrations, and water quality one year before and after the reintroduction of beavers to enclosures in 2023 at Spains Hall Estate, Essex. We found that beaver introduction increased terrestrial invertebrate abundance and species richness, but had no significant impacts on freshwater invertebrate richness, abundance, or biomass at the community level. On a finer scale, beavers significantly impacted freshwater indicator taxa abundance and biomass and shifted functional feeding group biomass. This was mediated by abiotic shifts, including an increase in dissolved oxygen content throughout the stream and increases in nitrate concentrations in the water near the dams. Downstream, nitrate was nearly completely filtered out from the water following beaver introduction. Our study shows that impacts are felt at least a kilometer downstream of beaver dams, and that beavers’ benefits to water quality and invertebrate biodiversity start to accumulate within only a year. These findings highlight the ecological value of beaver reintroductions and support their consideration in UK conservation planning.
Dataset DOI: 10.5061/dryad.15dv41pc0
Description of the data and file structure
We collected a series of water quality, algal concentration, and invertebrate samples 1 year before and 1 year after a beaver reintroduction trial in Spains Hall Estate, Essex, UK. The water quality and algal samples can be seen in allraw.csv, and have been processed into alldata.csv according to InvertebrateDataManipulation.Rmd. The other .csv files contain information from the invertebrate sampling. Aquatic invertebrates were collected using Surber sampling, while terrestrial invertebrates were sampled using pan and pitfall traps, left out for 24 hours.
CompleteBeaverAnalysis.R contains all GLMs and graphs.
InvertebrateDataManipulation.R is the data cleaning script on the raw data
SurberSamples.csv, allraw.csv, panTrap.cs,v and pitfall.csv are raw files.
alldata.csv, inverts.csv, panTrapalldata.csv and pitfalldata.csv are all processed, cleaned files.
NOTE: In some of the raw files, some cells are empty. This is because some sites had differing numbers of replicates for different abiotic variables, or that invertebrate identification was not done to that fine a specification. The R code in InvertebrateDataManipulation.Rmd manipulates these into NA values, so running this script remedies any blank cells.
Files and variables
File: CompleteBeaverAnalysis.R
Description: produces all GLMs and graphs used in the paper
File: InvertebrateDataManipulation.Rmd
Description: the data cleaning script on the raw data, including turning blank cells into NA where necessary.
File: SurberSamples.csv
Description: long-format raw file of the invertebrates sampled using Surber sampling
Variables
- Before/After: before or after beavers were reintroduced
- Site: Upstream, North (Dammed Enclosure), South (Undammed Enclosure), Downstream
- Sample number: which of the 5 repeated samples each invertebrate is contained in.
- Phylum/Class/Order/Family/Genus/Species: taxonomic ID to the finest practical resolution according to keys from Craig Macadam (2010), Brock, Paul D. (2021), Dobson (2012), Hildrew and Edington (1995), Wallace (2003), and Elliott (2010) [see paper citations].
- Body Length (mm)/Body Width (mm)/Head Width (mm): body measurements calculated using mm graph paper below the Petri dish.
- Notes: any further ID notes or notes on the condition of samples.
File: alldata.csv
Description: processed data file containing data compiled for each site
Variables
- time: before or after beaver introduction
- Site: Upstream, North (Dammed Enclosure), South (Undammed Enclosure), Downstream
- Date.Time: [not used for analysis] specific time stamp of sampling in the Before period according to the BenthoTorch
- Latitude/Longitude: coordinates of repeat samples for algal sampling as read by the BenthoTorch
- Sample: repeat sample number
- Cyano.g.cm./Green.Algae.g.cm./Diatoms.g.cm./Total.Conc..g.cm.: algal concentrations as read by the BenthoTorch
- oxygen..ppm./temperature/conductivity.S.cm./pH/TDS..ppm./nitrate..ppm./phosphate..mg.L.: abiotic variables measuring stream water quality
- avgPhosphate: [unused] mean phosphate across any repeats within a site.
- oxygen....: [unused] oxygenation percentage, for Before sites only.
- invertBiomass/evenness: NA
- Site_Time: Combination of site and time for the sake of the code, e.g. afterSouth or beforeNorth
- Site_numeric: 1 (Upstream), 2 (North/Dammed Enclosure), 3 (South/Undammed Enclosure), 4 (Downstream)
- speciesDiv/genusDiv/familyDiv/orderDiv/classDiv/phylumDiv: richness of each taxon
- richness.y: invertebrate taxonomic richness, calculated as the sum of the number of unique entries in each taxonomic level. Where a unique entry at a finer taxonomic resolution is used, its coarser resolutions were not used, so as not to double-count (e.g. counting Coleoptera as well as Arthropoda in the same measurement)
- shannon_index: [unused] Shannon index of diversity
- EPTrichness: richness as calculated above, specifically of members of Ephemeroptera, Plecoptera & Trichoptera
- EPTshannon_index: [unused] Shannon index of diversity specifically of members of Ephemeroptera, Plecoptera & Trichoptera
- abundance: number of individuals in the sample
- EPTabundance: number of EPT individuals in the sample
- totalBiomass: biomass of the sample, as calculated by O’Gorman et al. 2017; Benke et al. 1999 (see paper)
- EPTtotalBiomass: biomass of the EPT individuals in the sample, as calculated by O’Gorman et al. 2017; Benke et al. 1999 (see paper)
File: allraw.csv
Description: raw data which is processed to make the alldata.csv file
Variables
- See alldata.csv notes
Files: panTrapalldata.csv, pitfalldata.csv
Description: short-form data compiled from pan trap and pitfall trap data.
Variables
- See alldata.csv notes
Files: pitfall.csv, panTrap.csv
Description: raw long-form data of individuals identified in the pitfall & pan traps
Variables
- Before.After: before or after beavers were introduced
- Site: Upstream, North (Dammed Enclosure), South (Undammed Enclosure), Downstream
- Sample.number: 1,3,5 (since only 3 out of the 5 samples collected could be identified within the timeframe). The individual trap within which the individuals were identified.
- P: [unused] the original number of the trap before being reassigned 1,3,5 for data handling purposes.
- Phylum/Class/Order/Family/Genus/Species: taxonomic ID
- Notes: any further ID notes or notes on the condition of samples.
- finest_taxonomic_resolution: finest taxonomic resolution of each sample
- family_taxonomic_resolution: the family of each sample
- order_taxonomic_resolution: the order of each sample
File: inverts.csv
Description: raw long-form data of individuals identified in the Surber samples
Variables
- See pitfall.csv, panTrap.csv and surberSamples.csv
Code/software
All visualisation and analysis were completed in R (v 4.3.2, www.R-project.org/). Packages ggplot, cowplot, dplyr, grid, forcats, and purrr were used for data handling and presentation, and emmeans, goeveg, lme4, and easystats were used for modelling and statistical analysis.
