• Energy Research
• Open Access close
• AT close
• TH close
• Aurora Universities Network close

Complete replication data and code for article "A Meta-Analysis of Country Level Studies on Environmental Change and Migration". The rdata file contains both the meta and country level data. The data is also saved separately as xlsx files.

The response of the single-celled ciliates to increased temperature during global warming is critical for the structure and functioning of freshwater food webs. I conducted a meta-analysis of the literature from field studies and experimental evidence to assess the parameters characterising the thermal response of freshwater ciliates. The shape of the thermal performance curve predicts the ciliates’ survival at supraoptimal temperatures (i.e., the width of the thermal safety margin, TSM). The ciliates’ typical TSM is ~5°C. One-third of the freshwater ciliates dwelling permanently or occasionally in the pelagial cannot survive at temperatures exceeding 30°C. Likewise, cold-stenothermic species, which represent a significant fraction of euplanktonic ciliates, cannot survive by evolutionary adaptation to rapidly warming environments. The statistical analysis revealed that the ciliates’ thermal performance is affected by their planktonic lifestyle (euplanktonic versus tychoplanktonic), ability to form cysts, and nutritional ecology. Bactivorous ciliates have the widest temperature niche, and algivorous ciliates have the narrowest temperature niche. Phenotypic plasticity and genetic variation, favouring the selection of pre-adapted species in a new environment, are widespread among freshwater ciliates. However, the lack of evidence for the temperature optima and imprecisely defined tolerance limits of most species hamper the present analysis. The extent of acclimation and adaptation requires further research with more ciliate species than the few chosen thus far. Recent eco-evolutionary experimental work and modelling approaches demonstrated that the ciliates’ thermal responses follow general trends predicted by the metabolic theory of ecology and mechanistic functions inherent in enzyme kinetics. The present analysis identified current knowledge gaps and avenues for future research that may serve as a model study for other biota. Thermal adaptation may conflict with adaptation to other stressors (predators, food availability, pH), making general predictions on the future role of freshwater ciliates in a warmer environment difficult, if not impossible, at the moment. # Data from: Thermal response of freshwater ciliates: can they survive at elevated lake temperatures? [https://doi.org/10.5061/dryad.jdfn2z3jr](https://doi.org/10.5061/dryad.jdfn2z3jr) The dataset results from a meta-analysis to assess the parameters characterising the thermal response of freshwater ciliates (i.e., minimum and maximum temperature tolerated, temperature niche breadth). Cyst formation, the nutritional type, and the planktonic lifestyle were considered as factors affecting the ciliates’ thermal performance. ## Description of the data and file structure The main dataset reporting ciliate species and synonyms, taxonomic affiliation, minimum and maximum temperature and the temperature range tolerated, cysts formation, mixotrophic nutrition, food type, and planktonic lifestyle are reported in the 'Dataset_v4.xlsx' file. This is the main document. Taxonomic affiliation (i.e., order) following Adl et al. (2019, reference [65]J, the GBIF Backbone Taxonomy, and Lynn (2008; reference [66]). Details on the references - i.e., authors, publication year, title, journal/book, volume, and page/article numbers used to compile this dataset and some comments can be found in 'References.xlsx'. Empty cells mean that information is unavailable. References A-E are the main sources of the dataset, i.e., comprehensive review articles published by W. Foissner and colleagues in the 1990s. References 1-64 are case studies, published mainly after 1999. References 65 and 66 refer to the taxonomic affiliation of the ciliate species. More details about each column of the main document can be found in the 'Units_table.xlsx' file. ## Sharing/Access information Data was derived from the following sources: * ISI Web of Science (All Data Bases) * Google Scholar ## Code/Software R statistical software (v 4.0.5, R Core Team 2021) with the packages lme4, lmtest, multcomp, AICcmodavg. WebPlotDigitizer (Version 4.6) for data extraction from figures ## Version changes **06-aug-2024**: Taxonomic affiliation (order) corrected according to GBIF. Genus *Tintinnidium* is now in the order Oligotrichida. I scrutinised the detailed literature compilations by Foissner and colleagues published in the 1990s; these references are listed as primary sources A-E in the Dataset, see References.xlsx and README.txt) to obtain an overview of the thermal performance, resting cyst formation, and nutritional ecology of planktonic freshwater ciliates. I then searched the ISI Web of Science (All Data Bases) for updates and cross-references of Foissner’s works and further temperature records from (mainly) field studies. Search terms (in all fields) for the latter were ciliate* AND temperature NOT marine NOT ocean NOT soil NOT parasit* (1,339 hits). I followed the PRISMA guidelines in combination with EndNote 20 to filter out the records eligible for screening and analysis. Temperature data for assessing the minimum (Tmin) and maximum temperature (Tmax) of occurrence were eventually extracted from 68 publications. However, because Foissner’s works present extensive reviews, the actual number of publications used for the analysis is much higher. The final dataset obtained from field studies comprised 206 ciliate species. Next, I searched the ISI Web of Science for experimental results, using ciliate* AND temperature AND growth rate* NOT marine as search terms (218 records). Removing results from unsuitable research areas (mainly from medical research) reduced the records to 71 publications, which were screened. The combination of ciliate* AND numerical response NOT marine yielded 40 studies, ciliate* AND thermal performance 21 hits. I checked the selected articles for citations and cross-references using Google Scholar to identify any publications that might have slipped my attention. Eventually, I picked experimental results from 18 studies. If the literature data were only shown in figures, I extracted the data from the plots with WebPlotDigitizer (Version 4.6). I analysed the dataset with the R Statistical Software using the packages lme4, lmerTest, stats, multcomp, AICcmodavg and car.

Contrasting physiological mortality with predator-induced mortality is of tremendous importance for the population dynamics of many organisms but is difficult to assess. I performed a meta-analysis using planktonic ciliates as model organisms to estimate the maximum physiological mortality rates (δmax) across pelagic ecosystems in relation to environmental and biotic factors. Data were compiled from published numerical response (NR) experiments and experimentally determined rates of decline (ROD). Variables reported are ciliate species and order, ciliate specific growth rates (rmax), prey species, temperature, habitat (marine vs freshwater), the coefficients of the numerical response experiments, and reported or calculated ciliate mortality rates. The median δmax of planktonic ciliates was 0.62 d−1 and did not differ between marine and freshwater species. Maximum ciliate mortality rates were species-specific and affected by their rmax, cell volume, and ability to encyst. Cyst-forming species had, on average, higher δmax than species unable to encyst. Maximum mortality rates of ciliates were positively related to rmax but appeared unaffected by temperature. I conclude that (i) in the ocean, physiological mortality is more critical for controlling ciliate population size than ciliate losses imposed by microcrustacean predation, but (ii) in many lakes, the opposite holds; (iii) cyst-formation is an effective ciliate trait to cope with the high mortality of motile cells upon starvation. The lack of a temperature effect on δmax deserves further study; if correct, planktonic ciliates may take advantage of rising ocean and lake temperatures, with important implications for the pelagic food web. I used ISI Web of Science and Google Scholar to search for experiments that measured growth and mortality rates of ciliates as a function of prey concentration (i.e. numerical responses). The search terms were “growth (rate)” or “numerical response” in combination with “ciliate*” to search for numerical response experiments and “starvation” or “starved” in combination with “ciliate*” to search for mortality experiments. In addition, I searched the literature cited in these publications for further datasets. I considered only planktonic ciliates. When studies did not report all parameters of the NR curve, the data were extracted from figures with DataThief III or WebPlotDigitizer (Version 4.6) and fitted with a modified Michaelis-Menten equation that included the threshold prey concentration (P’) as an additional parameter. Mortality rates obtained by ROD experiments used the δmax reported in the respective study or calculated δmax from the maximum rate of decline after digitizing the data from the original curves, as described above. The literature search yielded δmax reported from 41 studies investigating 56 species or strains in 81 NR experiments and 19 ROD experiments. The final dataset (n = 77) included 37 studies and 48 species. I analyzed the dataset using the R Statistical Software using the packages lme4, lmerTest, AICcmodavg, and MuMIn. # Physiological mortality rates of planktonic ciliates ## Description of the Data and file structure I used ISI Web of Science and Google Scholar to search for experiments that measured growth and mortality rates of ciliates as a function of prey concentration (i.e. numerical responses). The main dataset containing available experimental studies reporting ciliate species, experimental temperature, prey species, ciliate maximum growth rates, ciliate cell volumes, habitat of ciliate isolation, method of study and reported or calculated ciliate mortality rates are reported in the 'Dataset_v2.xlsx' file. This is the main document. Missing data codes: N.A. = not available; n/a = not applicable. More details about each column of the main document can be found in the 'Units_table.xlsx' file. Details on the references - i.e. authors, publication year, title, journal/book, volume and page/article numbers - used to compile this dataset can be found in 'References.xlsx'. ## Sharing/access Information The individual data were derived mainly from the ISI Web of Science. The data compilation is novel. Excel, R

Abstract. Reliable quantification of the sources and sinks of greenhouse gases, together with trends and uncertainties, is essential to monitoring the progress in mitigating anthropogenic emissions under the Paris Agreement. This study provides a consolidated synthesis of CH4 and N2O emissions with consistently derived state-of-the-art bottom-up (BU) and top-down (TD) data sources for the European Union and UK (EU27+UK). We integrate recent emission inventory data, ecosystem process-based model results, and inverse modelling estimates over the period 1990–2018. BU and TD products are compared with European National GHG Inventories (NGHGI) reported to the UN climate convention secretariat UNFCCC in 2019. For uncertainties, we used for NGHGI the standard deviation obtained by varying parameters of inventory calculations, reported by the Member States following the IPCC guidelines recommendations. For atmospheric inversion models (TD) or other inventory datasets (BU), we defined uncertainties from the spread between different model estimates or model specific uncertainties when reported. In comparing NGHGI with other approaches, a key source of bias is the activities included, e.g. anthropogenic versus anthropogenic plus natural fluxes. In inversions, the separation between anthropogenic and natural emissions is sensitive to the geospatial prior distribution of emissions. Over the 2011–2015 period, which is the common denominator of data availability between all sources, the anthropogenic BU approaches are directly comparable, reporting mean emissions of 20.8 Tg CH4 yr−1 (EDGAR v5.0) and 19.0 Tg CH4 yr−1 (GAINS), consistent with the NGHGI estimates of 18.9 ± 1.7 Tg CH4 yr−1. TD total inversions estimates give higher emission estimates, as they also include natural emissions. Over the same period regional TD inversions with higher resolution atmospheric transport models give a mean emission of 28.8 Tg CH4 yr−1. Coarser resolution global TD inversions are consistent with regional TD inversions, for global inversions with GOSAT satellite data (23.3 Tg CH4yr−1) and surface network (24.4 Tg CH4 yr−1). The magnitude of natural peatland emissions from the JSBACH-HIMMELI model, natural rivers and lakes emissions and geological sources together account for the gap between NGHGI and inversions and account for 5.2 Tg CH4 yr−1. For N2O emissions, over the 2011–2015 period, both BU approaches (EDGAR v5.0 and GAINS) give a mean value of anthropogenic emissions of 0.8 and 0.9 Tg N2O yr−1 respectively, agreeing with the NGHGI data (0.9 ± 0.6 Tg N2O yr−1). Over the same period, the average of the three total TD global and regional inversions was 1.3 ± 0.4 and 1.3 ± 0.1 Tg N2O yr−1 respectively, compared to 0.9 Tg N2O yr−1 from the BU data. The TU and BU comparison method defined in this study can be operationalized for future yearly updates for the calculation of CH4 and N2O budgets both at EU+UK scale and at national scale. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.4288969 (Petrescu et al., 2020).

Abstract. Regional land carbon budgets provide insights on the spatial distribution of the land uptake of atmospheric carbon dioxide, and can be used to evaluate carbon cycle models and to define baselines for land-based additional mitigation efforts. The scientific community has been involved in providing observation-based estimates of regional carbon budgets either by downscaling atmospheric CO2 observations into surface fluxes with atmospheric inversions, by using inventories of carbon stock changes in terrestrial ecosystems, by upscaling local field observations such as flux towers with gridded climate and remote sensing fields or by integrating data-driven or process-oriented terrestrial carbon cycle models. The first coordinated attempt to collect regional carbon budgets for nine regions covering the entire globe in the RECCAP-1 project has delivered estimates for the decade 2000–2009, but these budgets were not comparable between regions, due to different definitions and component fluxes reported or omitted. The recent recognition of lateral fluxes of carbon by human activities and rivers, that connect CO2 uptake in one area with its release in another also requires better definition and protocols to reach harmonized regional budgets that can be summed up to the globe and compared with the atmospheric CO2 growth rate and inversion results. In this study, for the international initiative RECCAP-2 coordinated by the Global Carbon Project, which aims as an update of regional carbon budgets over the last two decades based on observations, for 10 regions covering the globe, with a better harmonization that the precursor project, we provide recommendations for using atmospheric inversions results to match bottom-up carbon accounting and models, and we define the different component fluxes of the net land atmosphere carbon exchange that should be reported by each research group in charge of each region. Special attention is given to lateral fluxes, inland water fluxes and land use fluxes.

Monocultural rubber plantations have replaced tropical forest, causing biodiversity loss. While protecting intact or semi-intact biodiverse forest is paramount, improving biodiversity value within the 11.4 million hectares of existing rubber plantations could offer important conservation benefits, if yields are also maintained. Some farmers practice agroforestry with high-yielding clonal rubber varieties to increase and diversify incomes. Here, we ask whether such rubber agroforestry improves biodiversity value or affects rubber yields relative to monoculture. We surveyed birds, fruit-feeding butterflies and reptiles in 25 monocultural and 39 agroforest smallholder rubber plots in Thailand, the world’s biggest rubber producer. Management and vegetation structure data were collected from each plot, and landscape composition around plots was quantified. Rubber yield data were collected for a separate set of 34 monocultural and 47 agroforest rubber plots in the same region. Reported rubber yields did not differ between agroforests and monocultures, meaning adoption of agroforestry in this context should not increase land demand for natural rubber. Butterfly richness was greater in agroforests, where richness increased with greater natural forest extent in the landscape. Bird and reptile richness were similar between agroforests and monocultures, but bird richness increased with the height of herbaceous vegetation inside rubber plots. Species composition of butterflies differed between agroforests and monocultures, and in response to natural forest extent, while bird composition was influenced by herbaceous vegetation height within plots, the density of non-rubber trees within plots (representing agroforestry complexity), and natural forest extent in the landscape. Reptile composition was influenced by canopy cover and open habitat extent in the landscape. Conservation priority and forest-dependent birds were not supported within rubber. Synthesis and applications. Rubber agroforestry using clonal varieties provides modest biodiversity benefits relative to monocultures, without compromising yields. Agroforests may also generate ecosystem service and livelihood benefits. Management of monocultural rubber production to increase inter-row vegetation height and complexity may further benefit biodiversity. However, biodiversity losses from encroachment of rubber onto forests will not be offset by rubber agroforestry or rubber plot management. This evidence is important for developing guidelines around biodiversity-friendly rubber and sustainable supply chains, and for farmers interested in diversifying rubber production. The accompanying ReadMe.txt file explains the contents of each .csv file, including definitions of each column. Sampling protocols are outlined in the paper in Journal of Applied Ecology.

Abstract. Dynamic global vegetation models are used to predict the response of vegetation to climate change. They are essential for planning ecosystem management, understanding carbon cycle–climate feedbacks, and evaluating the potential impacts of climate change on global ecosystems. JULES (the Joint UK Land Environment Simulator) represents terrestrial processes in the UK Hadley Centre family of models and in the first generation UK Earth System Model. Previously, JULES represented five plant functional types (PFTs): broadleaf trees, needle-leaf trees, C3 and C4 grasses, and shrubs. This study addresses three developments in JULES. First, trees and shrubs were split into deciduous and evergreen PFTs to better represent the range of leaf life spans and metabolic capacities that exists in nature. Second, we distinguished between temperate and tropical broadleaf evergreen trees. These first two changes result in a new set of nine PFTs: tropical and temperate broadleaf evergreen trees, broadleaf deciduous trees, needle-leaf evergreen and deciduous trees, C3 and C4 grasses, and evergreen and deciduous shrubs. Third, using data from the TRY database, we updated the relationship between leaf nitrogen and the maximum rate of carboxylation of Rubisco (Vcmax), and updated the leaf turnover and growth rates to include a trade-off between leaf life span and leaf mass per unit area.Overall, the simulation of gross and net primary productivity (GPP and NPP, respectively) is improved with the nine PFTs when compared to FLUXNET sites, a global GPP data set based on FLUXNET, and MODIS NPP. Compared to the standard five PFTs, the new nine PFTs simulate a higher GPP and NPP, with the exception of C3 grasses in cold environments and C4 grasses that were previously over-productive. On a biome scale, GPP is improved for all eight biomes evaluated and NPP is improved for most biomes – the exceptions being the tropical forests, savannahs, and extratropical mixed forests where simulated NPP is too high. With the new PFTs, the global present-day GPP and NPP are 128 and 62 Pg C year−1, respectively. We conclude that the inclusion of trait-based data and the evergreen/deciduous distinction has substantially improved productivity fluxes in JULES, in particular the representation of GPP. These developments increase the realism of JULES, enabling higher confidence in simulations of vegetation dynamics and carbon storage.