Data from: Earthworms do not increase greenhouse gas emissions (CO2 and N2O) in an ecotron experiment simulating a realistic three-crop rotation system
Data files
Nov 08, 2023 version files 1.74 MB
-
ew_extraction_biomass.csv
4.64 KB
-
ew_ghg_gram_day.csv
1.73 MB
-
README.md
5.41 KB
Dec 01, 2023 version files 2.74 MB
-
ew_extraction_biomass.csv
4.64 KB
-
ew_ghg_gram_day.csv
2.73 MB
-
README.md
5.41 KB
Dec 18, 2023 version files 2.74 MB
-
ew_extraction_biomass.csv
4.69 KB
-
ew_ghg_gram_day.csv
2.73 MB
-
README.md
5.41 KB
Abstract
Earthworms are known to stimulate soil greenhouse gas (GHG) emissions, but the majority of previous studies have used simplified model systems or lacked continuous high-frequency measurements. To address this, we conducted a two-year study using large lysimeters (5 m2 area and 1.5 m soil depth) in an ecotron facility, continuously measuring ecosystem-level CO2, N2O, and H2O fluxes. We investigated the impact of endogeic and anecic earthworms on GHG emissions and ecosystem water use efficiency (WUE) in a simulated agricultural setting. Although we observed transient stimulations of carbon fluxes in the presence of earthworms, cumulative fluxes over the study indicated no significant increase in CO2 emissions. Endogeic earthworms reduced N2O emissions during the wheat culture (-44.6%), but this effect was not sustained throughout the experiment. No consistent effects on ecosystem evapotranspiration or WUE were found. Our study suggests that earthworms do not significantly contribute to GHG emissions over a two-year period in experimental conditions that mimic an agricultural setting. These findings highlight the need for realistic experiments and continuous GHG measurements.
README: Earthworms greenhouse gas emissions in an ecotron experiment datasets
https://doi.org/10.5061/dryad.mgqnk9955
The experimental setup involved applying an earthworm treatment with three levels (i) anecic earthworms, (ii) endogeic earthworms, and (iii) a control with very low earthworm biomass (as the total elimination of earthworm is unfeasible) in large lysimeters. The experiment simulated a three plant species crop rotation consisting of a succession of Triticum aestivum (winter wheat) - Sinapsis alba (winter white mustard) - Zeamays (maize). The experiment was conducted in the European Ecotron of Montpellier (Montferrier-sur-Lez, France, www.ecotron.cnrs.fr). It utilized the Macrocosms experimental platform consisting in twelve identical and independent units, each consisting of an 30 m³ aboveground compartment enclosed by a highly transparent material to light and UV radiation (250 μm thick Teflon-FEP film, DuPont, USA). Additionally, it included a belowground compartment containing a 5 m2 stainless steel lysimeter hosting 14 t of soil (volume of 7.5 m3).
Each treatment level was replicated four times in four separate domes.
Description of the data and file structure
Two datasets are used in this experiment.
(i) The file "ew_ghg_gram_day.csv" contains the measurements for all green house gases and the fresh earthworm biomass measured at the end of each crop growing season and used as a covariable in the mixed model. The columns are described in the table below :
Variable | Levels_number | Description | Units |
---|---|---|---|
mac | 12 | macrocosm identity | number; 1-12 |
date | 664 | day of the measurements | yyyy-mm-dd |
Week | 96 | week of the year | number; 1-96 |
EW_T | 3 | earthworm treatment levels: Control, Anecic, Endogeic | unitless |
crop | 5 | Cultivated crop and intercrop: wheat, mustard, maize, intercrop_1, intercrop_2 | unitless |
month | 12 | month of the measurements | Abbreviated month name |
EW_BM | 36 | eartworm fresh weight biomass | g m-2 |
N2O | 7018 | concentration in N2O | mg N-N2O m-2 day-1 |
NEE | 7297 | net ecosystem exchange; concentration in CO2 | g C-CO2 m-2 day-1 |
Reco | 7304 | ecosystem respiration; concentration in CO2 | g C-CO2 m-2 day-1 |
GPP | 7304 | gross primary production; concentration in CO2 | g C-CO2 m-2 day-1 |
ET | 1208 | evapotranspiration | kg H2O m-2day-1 |
WUE | 7940 | water use efficiency | g C-CO2 kg H2O m-2 day-1 |
AirHR | 1189 | air relative humidity | % |
AirTemp | 1188 | air temperature | Celsius degrees |
Rg | 100 | global radiation | W m-2 |
SWC_13 | 1189 | soil water content at 13cm depth | % |
SWC_25 | 1185 | soil water content at 25cm depth | % |
SWC_80 | 1189 | soil water content at 80cm depth | % |
SoilTemp_13 | 1189 | soil temperature at 13cm depth | Celsius degrees |
SoilTemp_25 | 1189 | soil temperature at 25cm depth | Celsius degrees |
SoilTemp_80 | 1189 | soil temperature at 80cm depth | Celsius degrees |
VPD | 1189 | vapor pressure deficit | kPa |
NA values for all variables are missing values due to recording device malfunctioning.
(ii) the file "ew_extraction_biomass.csv" contains the measurements of fresh earthworm biomass as they were measured, extracted or added at several periods during the experiment. The columns are described in the table below :
Variable | Levels_number | Description | Units |
---|---|---|---|
crop | 4 | Cultivated crop and intercrop: wheat, mustard, maize, fallow | unitless |
EW_T | 3 | earthworm treatment levels: Control, Anecic, Endogeic | unitless |
mac | 12 | macrocosm identity | number; 1-12 |
bm | 56 | eartworm fresh weight biomass | g m-2 |
date | 4 | month of the measurements | Month name yyyy |
status | 2 | Status of the action on the earthworm fresh weight biomass: either measured, added or removed | NA |
Code/Software
We provide here the script ("EW_GHG_statistical_analysis_final.rmd") and the html version ("EW_GHG_statistical_analysis_final.html") of its reporting that contains the statistical analyses, plots and model output tables for the paper article titled "Earthworms do not increase greenhouse gas emissions (CO2 and N2O) in an ecotron experiment simulating a realistic three-crop rotation system".
The analysis was conducted using R (version 4.2.1 (2022-06-23 ucrt)) and RStudio (version 2023.06.2 Build 561) and the following R packages: kableExtra (version 1.3), nlme (version 3.1 - 157), MASS (version 7.3-57), plotly (version 4.10.1), patchwork (version 1.1.2), plotrix (version 3.8-2), broom.mixed (version 0.2.9.4), multilevelmod (version 1.0.0), suncalc (version 0.5.1), ggtext (version 0.1.2), multcomp (version 1.4-22), gginnards (version 0.1.1), lemon (version 0.4.6), tidyverse (version 2.0.0).
This script is structured in such a way that it enables the writing of the article using the html report produced by knitting the script without the need to re run parts of the script.
Methods
The CO2 concentration was measured at the inlet and outlet of each dome (every 12 min) using a multiplexer system coupled with two LI-7000 CO2/H2O analysers (LI-COR Biosciences, Lincoln, NE, USA). The Reichstein et al. (2005) C flux partitioning algorithm was used to estimate the daytime ecosystem respiration based on an exponential regression model (Lloyd et al., 1994). This allowed for the estimation of ecosystem respiration over 24 h (Reco = Reco_night + Reco_day) and gross primary production (GPP = NEE_day – Reco_day).
Ecosystem-level N2O fluxes were measured continuously as an open system using a TILDAS Compact Single analyser (N2O Aerodyne Research, Inc., USA). The analyzer was coupled to a multiplexer system allowing N2O fluxes measurement every 72 min for each Macrocosm. Evapotranspiration (ET) was computed as the lysimeter weight difference between two consecutive days. Four shear beam load cells per lysimeter (CMI-C3, Precia-Molen, Privas CEDEx France), with an accuracy of ±200 g, were used to measure the changes in mass. Ecosystem WUE was estimated as the ratio of GPP to ET derived from measurements by lysimeter weight changes over 24h.
Earthworms were sampled by utilizing the non-invasive octet electric method (Schmidt et al., 2001), which enabled sampling of a 1 m² surface area for a duration of 50 minutes using a customized version of the octet device manufactured by Electrotechnik Schuller (Darmstadt, Germany). To achieve coverage of the entire 5 m² surface area, five devices were simultaneously deployed.
References:
Reichstein, M. et al. On the separation of net ecosystem exchange into assimilation and ecosystem respiration: review and improved algorithm. Glob. Change Biol. 11, 1424–1439 (2005).
Lloyd, J. & Taylor, J. A. On the Temperature Dependence of Soil Respiration. Funct. Ecol. 8, 315–323 (1994).
Schmidt, O. Appraisal of the electrical octet method for estimating earthworm populations in arable land. Ann. Appl. Biol. 138, 231–241 (2001).