Data from: Species-specific responses cushion fish community body size against concurrent global change pressures
Data files
Jul 31, 2026 version files 248.98 KB
-
LMM_Oder_Jul_26.csv
160.20 KB
-
open_script_Oder_July__31_26.R
11.51 KB
-
README.md
16.42 KB
-
TaxaSpecific_Oder_Jul_26.csv
60.85 KB
Abstract
Multiple climate warming pressures can simultaneously affect biological communities and catalyze future disturbance. However, understanding how multiple pressures interact to affect ecological structure is challenging. Here, we evaluate the relationship between climate warming pressures (warming, hydrologic anomalies, and eutrophication), and fish community body size before and after a toxic algal bloom causing a large fish kill of a large-river ecosystem using 26 years of ecological data. We find that warming was associated with decreased community-level fish body size, yet changes in temperature and flow did not interact to affect body size. Furthermore, following the timing of an algal bloom, community-level body size increased for a given temperature range, caused by strong size increases of certain species. Thus, we find that at the community-level, deleterious effects of warming were not compounded by concurrent pressures, in part because of species-specific responses. Ultimately, our study provides evidence that warming decreases community-level body size, but such effects are not necessarily amplified by concurrent pressures within natural ecosystems. Data are organized in comma-delimited files containing the sizes of fish captures, taxon-specific relative biomass and community-level geometric mean body size, site-level and environmental covariates regarding temperature and discharge, and Principal Component Site scores.
This README provides data context for included data files by describing the organization of the dataset, the individual data files, and variables.
Reference Information
Provenance for this README
- File name: README.md
- Authors: Benjamin B. Tumolo
- Other contributors: Thomas Mehner, Paula Torre Zaffaroni, Sascha Starck, Christian Wolter
- Date created: 2026-07-28
Dataset Attribution and Usage
- Dataset Title: Data for the article "Species-specific responses cushion fish community body size against concurrent global change pressures"
- Dataset Contributors:
- Creators: Benjamin B. Tumolo, Thomas Mehner, Paula Torre Zaffaroni, Sascha Starck, Christian Wolter
- Date of Issue: 2026-07-28
- Suggested Citations:
-
Dataset citation:
Benjamin B. Tumolo, Thomas Mehner, Paula Torre Zaffaroni, Sascha Starck, Christian Wolter. Data for the article " Species-specific responses cushion fish community body size against concurrent global change pressures" 2026. , Dryad, Dataset
-
Contact Information
- Name: Benjamin B. Tumolo
- Affiliations: Department of Biology, Northern Michigan University, Marquette, MI 49855 USA
- Email: bbtumolo@gmail.com
Dates and Locations
- Dates of data collection: Data was collected from 1998 to 2024.
- Geographic locations of data collection: Main channel of the Oder River Germany.
Methodological Information
Methods of data collection/generation (see manuscript for more details):
We used fish mass (g) as the basis for patterns in body size.
We converted individual masses to community-level geometric mean mass (i.e., central tendency of mean body size) to examine patterns of body size at the community-level.
At each site, IGB quantitatively sampled fish from mid-channel river habitats by trawling. Trawling was performed using a bottom otter trawl with 8 x 1.5 m (width x height) net opening and 10 mm mesh size in the cod end pulled downstream at a speed of 8 km h-1.
Each haul length was 1000-2000 m. Collected fish were identified, their total length measured to the nearest mm for fish up to 10 cm length and nearest 0.5 cm for larger individuals.
Data and File Overview
Summary Metrics
- File count: 3
- Total file size: 319.55 KB
- File formats: .csv
Table of Contents
- LMM_Oder_Jul_26.csv
- TaxaSpecific_Oder_Jul_26.csv
- open_script_Oder_July__31_26.R
Setup
- Unpacking instructions: n/a
- Minumum Recommended software/tools: RStudio 2021.09.2; R version 3.6.3
File/Folder Details
Details for: LMM_Oder_Jul_26.csv
- Description: a comma-delimited file containing the relative biomass of fish species, community-level geometric body size, and environmental variables. Data are organized based sample site, sample date, and sample id.
- Format(s): .csv
- Size(s): 160.2 KB
- Dimensions: 361 rows x 65 columns
Variables:
- SITE_Event_ID: Unique sample site identification for each site on a specific date. The site event ID is composed of 2 parts: the site number or location and a specific date of sampling. The site location is a three-digit river km number followed by a month as 1-12 and then a four-digit year.
The sampling sites river KM are as follows: River KM 561, 568, 570, 572, 586, 590, 594, 597, 599, 601, 607, 609, 678, 682, 684, 686, 688, 693, 695, 697, 699. - Site: Sample site identification code based on river km and GPS coordinates latitude and longitude.
Sites and associated latitude/longitude are as follows. Note that coordinate data for endangered species are not included in this list.
KM 561 = 52.20451 14.69962
KM 568 = 52.25723 14.68590
KM 570 = 52.26245 14.65746
KM 572 = 52.26752 14.62962
KM 586 = 52.36412 14.55262
KM 590 = 52.39723 14.53432
KM 594 = 52.43160 14.54788
KM 597 = 52.44531 14.58282
KM 599 = 52.45944 14.60073
KM 601 = 52.47513 14.61227
KM 607 = 52.51745 14.61088
KM 609 = 52.53501 14.60525
KM 678 = 52.96815 14.16498
KM 699 = 53.09841 14.37011
- Stream_KM: numerical stream km site corresponding with location along the main channel of the Oder River. The above sample id's provide specific geographic latitude and longitude coordinates for each one of these values.
- Anguilla.anguilla_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Leuciscus.idus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Rhodeus.amarus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Barbus.barbus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Abramis.brama_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Squalius.cephalus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Ctenopharyngodon.idella_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Carassius.gibelio_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Perca.fluviatilis_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Gobio.gobio_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Blicca.bjoerkna_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Leuciscus.leuciscus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Esox.lucius_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Cyprinus.carpio_rel_biomass: relative biomass as a percentage or this specific taxa used for species matrix for PCA
- Gymnocephalus.cernua_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Coregonus.maraena_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Rutilus.rutilus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Osmerus.eperlanus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Leuciscus.aspius_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Romanogobio.belingi_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Alburnus.alburnus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Silurus.glanis_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Vimba.vimba_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Sander.lucioperca_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Ballerus.ballerus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Barbatula.barbatula_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Neogobius.melanostomus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Cobitis.taenia_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Tinca.tinca_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Gasterosteus.aculeatus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Scardinius.erythrophthalmus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Lota.lota_rel_biomass: relative biomass for this specific taxa used for species matrix for PCA
- Acipenser.oxyrinchus_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Salmo.trutta_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- Acipenser.baerii_rel_biomass: relative biomass as a percentage for this specific taxa used for species matrix for PCA
- geom_mean_comm: Geometric mean body mass (in grams) of fish community for a given sampling.
- Month: Numeric value of indicating the month of sampling event (4-12) April(4) through December(12) no samples are present from the months of Jan(1)-March(3).
- date_ymd: full date of sampling event from June-1-1998 to April-1-2024.
- Year: Year of sampling event from 1998 to 2024
- Qstation: Station name of discharge data collection. Daily discharge data starting in 1950 were obtained from the Global Runoff Data Centre and aggregated to a monthly resolution. Two gauging stations were used to represent sites upstream and downstream of the Warta River. River discharge at Gozdowice was used as representative value of the study section encompassed within Lower Oder River, downstream of the Warta River. This includes the bulk of the samples included in the analysis (n=317 out of 344, or 88% of the samples). A second gauging station at Eisenhüttenstadt was to characterize these more upstream sites (n=44, or 12% of observations).
- LT_average_Q: Long-term average discharge value in m3 s-1 from 1950 to 2024
- LT_month_average_Q: Long-term average discharge value for the given month in m3 s-1 from 1950 to 2024
- LT_season_average_Q: Long-term average discharge value for the given season in m3 s-1 from 1950 to 2024
- longtermAnomaly_Q: Long-term anomaly discharge value for the time series in m3 s-1 from 1950 to 2024. Anomalies were calculated as the difference between each value and the long-term mean for its corresponding month.
- monthAnomaly_Q: Long-term anomaly discharge value for the time series for a given month in m3 s-1 from 1950 to 2024. Anomalies were calculated as the difference between each value and the long-term mean for its corresponding month.
- LT_averageTemp: Long-term average surface temperature in (°C) Land surface temperature data was derived from the ERA5- Land reanalysis gridded at a 5 km resolution, covering the period 1950-2025
- LT_month_averageTemp: long-term average temperature value for the given month in (°C)
- longtermAnomalyTemp: Long-term anomaly temperature value for the time series in (°C), calculated as the difference between each value and the long-term mean for its corresponding month.
- monthAnomalyTemp: Long-term anomaly temperature value in (°C) for the time series for a given month calculated as the difference between each value and the long-term mean for its corresponding month.
- Season: Season of data collection
- Summer = June, July, August
- Fall =September, October, November
- Winter = December
- Spring =April, May
- Season_Range: Seasonal range of data collection where Fall and Winter were combined and Spring and Summer were combined
- FallWinter = September, October, November, December
- SpringSummer = April, May, June, July, August
- fish_n: The number value of fish collected on event in nets
- mean_mass: Average mass in grams of all fish collected on event
- mean_length: The avearge total length in cm of fish collected
- total_mass: summed total mass in grams of fish mass for event
- fished_length: estimate for fishing effort in haul length distance length in m. How far the trawl net was hauled and sampling for.
- fished_area: the estimated for fishing effort area in taking account distance of trawl and trawl net size.
- CPUE: No units. Estimated catch per unit effort of fishing event
- KM_Range2: Range of KM to distinguish upper and lower sections of the study area based on the Warta River Confluence.
- higher = above or upstream of the Warta confluence and includes KM 561-KM 601
- lower = below or downstream of the Warta confluence and includes KM 607-699
- Bloom: Timing of the large-scale toxic algal bloom that resulted in large-scale fish kill in August of 2022.
- pre = time period prior to the bloom including June 1 1998- May 1 2022
- post = time period after the bloom including November 1 2022- April 1 2024
- PC1: Principal component 1 specific coordinated extracted from PCA analysis
- PC2: Principal component 2 specific coordinated extracted from PCA analysis
Missing data codes: missing data or data that was not collected is indicated with NA
Details for: TaxaSpecific_Oder_Jul_26.csv
- Description: a comma-delimited file containing the measured relative abundances, relative biomass' and average mass for specific focal taxa.These taxa were focal taxa because they were identified by other analyses as being important to community composition.
- Format(s): .csv
- Size(s): 59.4 KB
- Dimensions: 640 rows x 12 columns
Variables:
- SITE_Event_ID: Unique sample site identification for each site on a specific date. Unique sample site identification for each site on a specific date. The site event ID is composed of 2 parts: the site number or location and a specific date of sampling. The site location is a three-digit river km number followed by a month as 1-12 and then a four-digit year.
The sampling sites river KM are as follows: River KM 561, 568, 570, 572, 586, 590, 594, 597, 599, 601, 607, 609, 678, 682, 684, 686, 688, 693, 695, 697, 699. - Site: Sample site identification code based on river km and GPS coordinates latitude and longitude.
Sites and associated latitude/longitude are as follows. Note that coordinate data for endangered species are not included in this list.
KM 561 = 52.20451 14.69962
KM 568 = 52.25723 14.68590
KM 570 = 52.26245 14.65746
KM 572 = 52.26752 14.62962
KM 586 = 52.36412 14.55262
KM 590 = 52.39723 14.53432
KM 594 = 52.43160 14.54788
KM 597 = 52.44531 14.58282
KM 599 = 52.45944 14.60073
KM 601 = 52.47513 14.61227
KM 607 = 52.51745 14.61088
KM 609 = 52.53501 14.60525
KM 678 = 52.96815 14.16498
KM 699 = 53.09841 14.37011
- Bloom: Timing of the large-scale toxic algal bloom that resulted in large-scale fish kill in August of 2022.
pre = time period prior to the bloom including June 1 1998- May 1 2022
post = time period after the bloom including November 1 2022- April 1 2024 - Taxa: the taxonomic names or classification of organisms that were collected during a sampling event
- rel_abundance: Measured relative abundance as a percentage of a particular taxa on for a given sampling event
- avg_mass: Measured average mass in g of a particular taxa on for a given sampling event
- rel_biomass: Measured relative biomass as a percentage of a particular taxa on for a given sampling event
- Month: Numeric value indicating the month of sampling event (4-12) April (4) through December(12); no samples are present from the months of Jan(1)-March(3).
- date_ymd: full date of sampling event from June-1-1998 to April-1-2024.
- Year: Year of sampling event from 1998 to 2024
- Season: Season of data collection
- Summer= June, July, August
- Fall=September, October, November
- Winter= December
Spring=April, May
- Season_Range: Seasonal range of data collection where Fall and Winter were combined and Spring and Summer were combined
- FallWinter = September, October, November, December
- SpringSummer= April, May, June, July, August
Missing data codes: missing data or data that was not collected is indicated with NA
Details for: open_script_Oder_July_31_26.R
- Description: a File containing code documentation used for analysis of the above comma-delimited files.
- Format(s): .R
- Size(s): 98.5 KB
- Dimensions: 291 lines
- There are Season-specific analyses that require subsetting of the larger dataset "LMM_Oder_Jul_26.csv" into "Season_Range" categories. Seasonal range of data collection where Fall and Winter were combined and Spring and Summer were combined
I have not provide these additional subseted dataframes, but they can be delineated accordingly.
FallWinter = September, October, November, December;
SpringSummer= April, May, June, July, August
END OF README
