Jonathan Lefcheck
Department of Biological Sciences
Virginia Institute of Mariner Science
The College of William & Mary
Gloucester Point, VA USA
E-mail: [email protected]

J. Emmett Duffy
Tennenbaum Marine Observatories Network
Smithsonian Institution
Washington, D.C. USA
E-mail: [email protected]

initial_species_data.csv (MD5: a18ff9ecba3416b181591c3e94277cf0)

final_species_data.csv (MD5: 37134e48f4d7b3de20bb25aaaecf5b4e)

final_afdm_data.csv (MD5: fc558a1f8f7a9db215b67126a94e8bd2)

trait_data.csv (MD5: bc4d4f19b0fb4377c70ec8f1d72afde3)

Lefcheck_mesocosm_analysis.R (MD5: cbcfe9d5f8cb4a457332f4279e115374)

dbFD_new.R (MD5: dd2b3639948332c1751a5e997c509739)

These supplements should allow the reader to exactly reproduce every analysis in the manuscript, including the generalized linear mixed effects models, analysis of individual traits, and structural equation models, as well as most graphics.

The data set is broken up into four files: (1) initial_species_data.csv contains the stocked abundances of all organisms used in the experiment; (2) final_species_data.csv contains the final abundances for stocked and recruiting invertebrates in their respective sieve size classes; (3) final_afdm_data.csv contains the final biomasses for stocked predators and Gracilaria spp., as well as recruiting algae and invertebrates; and (4) trait_data.csv contains the traits used in the analysis of functional diversity, derived from literature values and field observations.

The main analysis is in the Lefcheck_mesocosm_analysis.R file. The dbFD_new.R file is sourced from the main script and includes minor modifications to the dbFD() function in the FD package to account for missing cells when calculating functional diversity.

Column metadata

initial_species_data.csv

id: the number of species (mono = 1, 3, 6, or poly = 9), the functional diversity treatment (low or high), and finally the replicate number (1-6).

treatment: the level of functional diversity (individual species, low, high, or poly). ‘Low’ and ‘high’ values are followed by the number of species (3 or 6).

block: experimental block.

initial.fd: initial functional richness (min, low, high, or max).

initial.gracilaria.biomass: initial wet weight of Gracilaria spp.

initial.ampithoe.abund: initial abundance of Ampithoid spp.

initial.bittium.abund: initial abundance of Bittolium varium

initial.callinectes.abund: initial abundance of Callinectes sapidus

initial.erichsonella.abund: initial abundance of Erichsonella attenuata

initial.fundulus.abund: initial abundance of Fundulus heteroclitus

initial.gammarus.abund: initial abundance of Gammarus mucronatus

initial.hippolyte.abund: initial abundance of Hippolyte pleuracanthus

initial.palaemonetes.abund: initial abundance of Palaemonetes pugio

initial.syngnathus.abund: initial abundance of Syngnathus spp.

initial.callinectes.length: initial length of Callinectes sapidus (in mm)

initial.fundulus.length: initial length of Fundulus heteroclitus (in mm)

initial.hippolyte.biomass: initial biomass of Hippolyte pleuracanthus (in g)

initial.palaemonetes.biomass: initial biomass of Palaemonetes pugio (in g)

initial.syngnathus.length: initial length of Syngnathus spp. (in mm)

 

final_species_data.csv

date: date the sample was processed and weighed

id: the number of species (mono = 1, 3, 6, or poly = 9), the functional diversity treatment (low or high), and finally the replicate number (1-6).

group: the taxonomic group of the species (e.g., “Molluscs”, “Crustacean”, “Polychaete”,”Caprellids”) to be used in the calculation of ash-free dry mass from Edgar equations (see R script).

Species: species whose dry mass whose abundance was recorded.

8.0_mm: the number of individuals in the 8 mm sieve.

5.6_mm: the number of individuals in the 5.6 mm sieve.

4.0_mm: the number of individuals in the 4 mm sieve.

2.8_mm: the number of individuals in the 2.8 mm sieve.

2.0_mm: the number of individuals in the 2.0 mm sieve.

1.4_mm: the number of individuals in the 1.4 mm sieve.

1.0_mm: the number of individuals in the 1 mm sieve.

0.71_mm: the number of individuals in the 0.71 mm sieve.

0.5_mm: the number of individuals in the 0.5 mm sieve.

 

final_afdm_data.csv

date: date the sample was processed and weighed

id: the number of species (mono = 1, 3, 6, or poly = 9), the functional diversity treatment (low or high), and finally the replicate number (1-6).

species: species whose weight was recorded.

number.length: the number of individuals or length (in mm) of the species in the sample.

Tin: weight of drying tin by itself (g).

Tin_Dry: weight of drying tin plus dried sample (g).

Tin_Ash: weight of trying tin plus ashed sample (g).

AFDM: ash-free dry mass (g) = (Tin_Dry – Tin) – (Tin_Ash – Tin)

Notes: additional notes.

 

trait_data.csv

species: the species used in the experiment.

armor: whether the organisms is defended. Categorical: None, shell (chitin, calcium carbonate).

body.plan: the body plan of the organism. Categorical: Articulate (laterally-, ventrally-compressed, subcylindrical), shelled conic, filiform, fusiform).

trophic.level: The general trophic level of the organism. Categorical: Grazer, omnivore, predator.

max.biomass: Maximum observed biomass in field surveys of local eelgrass bed (mg).

mean.biomass: Average observed biomass in field surveys of local eelgrass bed (mg).

mobility: The average mobility of the organism. Categorical: Swimmer (low, high), tube-builder, crawler.

reproductive.mode: The reproductive mode of the organism. Categorical: Direct, planktotrophic, ovoviviparous, oviparous.

month.max.abund: The month of maximum abundance as observed in field surveys of local eelgrass bed (February – November).