• Energy Research
• 2. Zero hunger close
• CN close
• GB close
• HK close

Twenty soil cores were collected from a field site in Lincolnshire in March 2011, three weeks after planting and Nitrogen fertiliser addition. Soil cores of 150-180 millimetre (mm) depth, containing approximately 1.6 kilogram soil (dry weight) were extracted in Polyvinyl chloride (PVC) pipes (height 215 mm depth 102 mm) and stored at 4 degrees centigrade for 30 days. A four-treatment factorial experiment was designed using soils un-amended or amended with biochar and un-wetted or wetted with deionised water (5 replicates per treatment). Soil in all the cores was mixed to 7 centimetre (cm) depth. To half of the cores, biochar (less than 2 mm) was mixed into the soil at a rate of 3 percent soil dry weight (approximately 22 tons per hectare (t ha-1)). After allowing for any potential Carbon dioxide (CO2) flush from newly-mixed soil to equilibrate for seven days, the cores were placed at 16 degrees centigrade in the dark. Un-wetted soil cores were maintained at 23 percent Gravimetric moisture content (GMC), whilst the GMC of 'wetted' soil cores was increased to 28 percent GMC at the time zero (t0) of four wetting events on day 17, 46, 67 and 116. These water addition rates were based on mean and maximum monthly soil GMC measured in the field between 2009-2010. Data from an investigation of the effects of biochar application to soil on greenhouse gas emissions using soil from a bioenergy crop (Miscanthus X. giganteus). Data include physical (bulk density) and chemical analyses of the soil (total carbon (C) and nitrogen (N), extractable ammonium and nitrate), and greenhouse gas (GHG) emissions (carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O)) during incubations. Data were collected during two incubation experiments investigating the effects of temperature, soil moisture and soil aeration on biochar induced suppression of GHG emissions. Biochar is a carbon rich substances which is being advocated as a climate mitigation tool to increase carbon sequestration and reduce nitrous oxide emissions.

Prediction of the potentially suitable areas of Leonurus japonicus habitability zones with maxent occurrence points：By sorting out the information of Leonurus japonicus specimens recorded in the Chinese Digital Herbarium (CVH, http://www.cvh.ac.cn/), and combining with the L. japonicus presence points in the Global Biodiversity Information Platform (GBIF, https://www.gbif.org/), the existing distribution positions of L. japonicus were preliminarily obtained, and then the corresponding latitude and longitude coordinates of each distribution point were obtained by Baidu coordinate system. All were used for modeling. environmental variables：Species’ ecoloical niches are affected by climate, topography, biology, and other factors. In consideration of the comprehensiveness and complexity of ecological factors, 34 environmental variables which could reflect species’ ecoloical niches were selected. The list included 19 bioclimatic factors, 14 soil factors and a topographic factor (altitude).The current (1970–2000), 2050s (2041–2060), and 2090s (2081–2100) bioclimatic factor data used in this research were derived from the world climate database Worldclim (http://www.worldclim.Org), and the pixel size of the data was 2.5 arc-minutes (-5 km). The climate data of the 2050s and 2090s were obtained from the Beijing Climate Center-Climate System Model-Medium Resolution (BCC-CSM2-MR), one of the Coupled Model Inter-Comparison Project Phase 6 (CMIP6) datasets, which included three scenarios: sustainable development (SSP126), intermediate development (SSP245) and conventional development (SSP585). SSP scenarios have a high accuracy and separation rate and can integrate local development factors, and so are more convincing than CMIP5 data. The data of soil factors and topographic factors were obtained form the World Soil Database (HWSD) of the FAO (http://www.fao.org/faostat/en/#data), and the provincial national vector map were from China’s Ministry of Natural Resources (http://www.mnr.gov.cn/). The environmental variables is in ASCii format. ASCii can be viewed using standard GIS software such as: environmental variables\\climate\\50126\\bio1.asc Naming convention: Type Variables Description UNITS Bio1 Annual Mean Temperature ℃×10 Bioclimatic Bio2 Mean Diurnal Range ℃×10 Variables Bio3 Isothermality 1 Bio4 Temperature Seasonality 1 Bio5 Max Temperatur ℃×10 Bio6 Min Temperature of Coldest Month ℃×10 Bio7 Temperature Annual Range ℃×10 Bio8 Mean Temperature of Wettest Quarter ℃×10 Bio9 Mean Temperature of Driest Quarter ℃×10 Bio10 Mean Temperature of Warmest Quarter ℃×10 Bio11 Mean Temperature of Coldest Quarter ℃×10 Bio12 Annual Precipitation mm Bio13 Precipitation of Wettest Month mm Bio14 Precipitation of Driest Month mm Bio15 Precipitation Seasonality 1 Bio16 Precipitation of Wettest Quarter mm Bio17 Precipitation of Driest Quarter mm Bio18 Precipitation of Warmest Quarter mm Bio19 Precipitation of Coldest Quarter mm T\_GRAVEL Topsoil Gravel Content %vol. Top Soil Variable T\_SAND Topsoil Sand Fraction % wt. T\_SILT Topsoil Silt Fraction % wt. T\_CLAY Topsoil Clay Fraction % wt. T\_USDA\_TEX\_CLASS Topsoil USDA Texture Classification name T\_REF\_BULK\_DENSITY Topsoil Reference Bulk Density kg/dm3 T\_OC Topsoil Organic Carbon % weight T\_PH\_H2O Topsoil pH (H2O) -log(H+) T\_CEC\_CLAY Topsoil CEC (clay) cmol/kg T\_CEC\_SOIL Topsoil CEC (soil) cmol/kg T\_BS Topsoil Base Saturation % T\_TEB Topsoil TEB cmol/kg T\_ESP Topsoil Sodicity (ESP) % T\_ECE Topsoil Salinity (Elco) dS/m Terrain ELEV Elevation m ENMeval package: To avoid overfitting due to the high complexity of the model constructed with the default parameters, which may cause the predicted distribution of the potential habitat of L. japonicus to deviate too much from the actual situation, this study used the ENMeval package in R 4.3.1, and adjusted the two most important parameters, namely, regularization multiplier (RM) and feature combination (FC), to improve the prediction accuracy of the model. CoordinateCleaner：The R software package ‘CoordinateCleaner’ was used to removing records without coordinate precision and suspected outliers. Based on the ‘subset’ ‘clean\_coordinates’ operation in CoordinateCleaner, we obtained the results of bias corrections on the datasets. SpThin package: Spatial thinning of species occurrence records can help address problems associated with spatial sampling biases. Ideally, thinning removes the fewest records necessary to substantially reduce the effects of sampling bias, while simultaneously retaining the greatest amount of useful information. Spatial thinning can be done manually; however, this is prohibitively time consuming for large datasets. Using a randomization approach, the ‘thin’ function in the spThin R package returns a dataset with the maximum number of records for a given thinning distance, when run for sufficient iterations. geosphere package:The geosphere package in the R environment was used to calculate the centroid range shift distance of L. japonicus under different climate change scenarios. SDMTools：The package in R language was used to calculate the location of centroid in the suitable area of Leonurus japonicus under 6 different economic paths in the current and future periods. VIF package: The usdm package provides a set of functions to support dealing with problematic situations in species distribution modelling (e.g., multicollinearity, positional uncertainty).To detect whether predictor variables are subjected to multicollinearity, you may use vif (variance inflation factor) metric, and some methods implemeted in this package including vifstep or vifcor (a stepwise procedure to identify collinear variables). Leonurus japonicus Houtt. is a traditional Chinese medicinal plant with high medicinal and edible value. Wild L. japonicus resources have been reduced dramatically in recent years. This study predicted the response of distribution range of L. japonicus to climate change in China, which provided the scientific basis for the conservation and utilization. In this study, 489 occurrence points of L. japonicus were selected based on GIS technology and spThin package. The default parameters of the Maxent model were adjusted by using ENMeva1 package of the R environment, and the optimized Maxent model was used to analyze the distribution of L. japonicus. When the feature combination in the model parameters is hing and the regularization multiplier is 1.5, the Maxent model has a higher degree of optimization. With the AUC of 0.830 our model showed a good predictive performance The results showed that L. japonicus was widely distributed in the current period. The maximum temperature of the warmest month, the minimum temperature of the coldest month, the precipitation of the wettest month, the precipitation of the driest month and altitude were the main environmental factors affecting the distribution of L. japonicus. Under the three climate change scenarios, the suitable distribution area of L. japonicus will range-shift to high latitudes, indicating that the distribution of L. japonicus has a strong response to climate change. The regional change rate is the lowest under the SSP126-2090s scenario and the highest under the SSP585-2090s scenario.

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.

Climate trends during maize growing period and their impacts on maize yield in Southern hills was investigated. This dataset contains: 1) information of stations in cultivation region for maize in Southern hills; 2) Trend in temperature and its effect on yield in cultivation region for maize in Southern hills; 3) Trend in radiation and its effect on yield in cultivation region for maize in Southern hills; 4) Trend in precipitation and its effect on yield in cultivation region for maize in Southern hills. Climate trends during maize growing period and their impacts on maize yield in Southern hills was investigated. This dataset contains: 1) information of stations in cultivation region for maize in Southern hills; 2) Trend in temperature and its effect on yield in cultivation region for maize in Southern hills; 3) Trend in radiation and its effect on yield in cultivation region for maize in Southern hills; 4) Trend in precipitation and its effect on yield in cultivation region for maize in Southern hills.

1. Predicted changes in the frequency and intensity of extreme rainfall events in the UK have the potential to disrupt terrestrial ecosystem function. However, responses of different trophic levels to these changes in rainfall patterns, and the underlying mechanisms, are not well characterised. 2. This study aimed to investigate how changes in both the quantity and frequency of rainfall events will affect the outcome of interactions between plants, insect herbivores (above- and below- ground) and natural enemies. 3. Hordeum vulgare L. plants were grown in controlled conditions and in the field, and subjected to three precipitation scenarios: ambient (based on a local 10 year average rainfall); continuous drought (40% reduction compared to ambient); drought/ deluge (40% reduction compared to ambient at a reduced frequency). The effects of these watering regimes and wireworm (Agriotes species) root herbivory on the performance of the plants, aphid herbivores above-ground (Sitobion avenae, Metapolophium dirhodum and Rhopalosiphum padi), and natural enemies of aphids including ladybirds (Harmonia axyridis) were assessed from measurements of plant growth, insect abundance and mass, and assays of feeding behaviour. 4. Continuous drought decreased plant biomass, whereas reducing the frequency of watering events did not affect plant biomass but did alter plant chemical composition. In controlled conditions, continuous drought ameliorated the negative impact of wireworms on plant biomass. 5. Compared to the ambient treatment, aphid mass was increased by 15% when feeding on plants subjected to drought/ deluge; and ladybirds were 66% heavier when feeding on these aphids but this did not affect ladybird prey choice. In field conditions, wireworms feeding below-ground reduced the number of shoot-feeding aphids under ambient and continuous drought conditions but not under drought/ deluge. 6. Predicted changes in both the frequency and intensity of precipitation events under climate change have the potential to limit plant growth, but reduce wireworm herbivory, while simultaneously promoting above-ground aphid numbers and mass, with these effects transferring to the third trophic level. Understanding the effect of future changes in precipitation on species interactions is critical for determining their potential impact on ecosystem functioning and constructing accurate predictions under global change scenarios. Controlled environment and field experimental dataData file containing all data reported in the paper including plant, soil and insect data from controlled environment and field experiments. First spreadsheet in the data file contains a key to explain all abbreviations used throughout the file.Experimental data.xlsx

The expansion of irrigated agriculture has increased global crop production but resulted in widespread stress to freshwater resources. Ensuring that increases in irrigated production only occur in places where water is relatively abundant is a key objective of sustainable agriculture, and knowledge of how irrigated land has evolved is important for measuring progress towards water sustainability. Yet a spatially detailed understanding of the evolution of global area equipped for irrigation (AEI) is missing. Here we utilize the latest sub-national irrigation statistics (covering 17298 administrative units) from various official sources to develop a gridded (5 arc-min resolution) global product of AEI for the years 2000, 2005, 2010, and 2015. We find that AEI increased by 11% from 2000 (297 Mha) to 2015 (330 Mha) with locations of both substantial expansion (e.g., northwest India, northeast China) and decline (e.g., Russia). Combining these outputs with information on green (i.e., rainfall) and blue (i.e., surface and ground) water stress, we also examine to what extent irrigation has expanded unsustainably (i.e., in places already experiencing water stress). We find that more than half (52%) of irrigation expansion has taken place in regions that were already water stressed, with India alone accounting for 36% of global unsustainable expansion. These findings provide new insights into the evolving patterns of global irrigation with important implications for global water sustainability and food security. Recommended citation: Mehta, P., Siebert, S., Kummu, M. et al. Half of twenty-first century global irrigation expansion has been in water-stressed regions. Nat Water (2024). https://doi.org/10.1038/s44221-024-00206-9 Open-access peer reviewed publication available at https://www.nature.com/articles/s44221-024-00206-9 Files G_AEI_*.ASC were produced using the GMIA dataset[https://data.apps.fao.org/catalog/iso/f79213a0-88fd-11da-a88f-000d939bc5d8]. Files MEIER_G_AEI_*.ASC were produced using Meier et al. (2018) dataset [https://doi.pangaea.de/10.1594/PANGAEA.884744].

Climate trends during maize growing period and their impacts on maize yield in Northwest China was investigated. This dataset contains: 1) information of stations in cultivation region for maize in Northwest China; 2) Trend in temperature and its effect on yield in cultivation region for maize in Northwest China; 3) Trend in radiation and its effect on yield in cultivation region for maize in Northwest China; 4) Trend in precipitation and its effect on yield in cultivation region for maize in Northwest China. Climate trends during maize growing period and their impacts on maize yield in Northwest China was investigated. This dataset contains: 1) information of stations in cultivation region for maize in Northwest China; 2) Trend in temperature and its effect on yield in cultivation region for maize in Northwest China; 3) Trend in radiation and its effect on yield in cultivation region for maize in Northwest China; 4) Trend in precipitation and its effect on yield in cultivation region for maize in Northwest China.

Study site This is a dataset of soil physiochemical properties, bacterial and fungal abundance, and above and belowground plant and invertebrate biomass, sampled at 40 soil plots in the Hengill geothermal valley, Iceland, from 15th to 22nd August 2018. The plots, measuring approximately 1 m2, evenly span a temperature gradient of 10-35°C. The dataset also includes data on the decomposition rate of soil organic matter, which was sampled at 60 plots in the Hengill valley from May to July 2015 (see Robinson et al. 2021 for plot details and sampling regime). Soil properties Soil temperature was measured at 5 cm depth at each plot on 15th, 18th, and 22nd August, and a mean plot temperature calculated. Soil physiochemical properties were analysed from 3 soil cores of 3 cm in diameter, taken from the upper 10 cm soil stratum at each plot; one quarter of each subsample was pooled to obtain an estimate per plot. Aboveground plant matter, excluding roots, were removed from each core. Percentage soil moisture was calculated by measuring the weight of one pooled soil sample before and after drying for 24 h in a 70°C drying oven. Soil pH was obtained from 20 g of the dry soil by adding 100 ml distilled water, shaking for 5 min on 150 rpm, letting the sample stand for 2 h, and measuring soil pH from the water layer using an InoLab pH 720 (WTW) probe. Soil PO4, NH4, and NO3 concentrations were analysed from a second pooled soil; 60 g of fresh soil was extracted in 100 ml distilled water, filtered through a GF/C (1.2μm) glass microfiber filter (Whatman, GE Healthcare Europe GmbH), and analysed using a Lachat QuikChem 8000 analyser (Zallweger Analytics, Inc., Lachat Instruments Division, USA). Total mineral N was calculated as the sum of NH4 and NO3. Soil organic matter content (excluding dry root biomass) was calculated as the weight lost from an oven dried (105°C for 24 hours) soil sample after heating at 550 °C for 5 h. Decomposition rate of soil organic matter was measured using the Cotton-strip Assay method (Tiegs et al. 2013) by placing a 2.5 cm x 8 cm strip of Fredrix-brand unprimed 12-oz. heavyweight cotton fabric (Style #548) 5 cm belowground at 60 plots, concurrently with a Maxim Integrated DS1921G Thermocron iButton temperature logger, on 13th May 2015. The strips were collected on 3rd July, rinsed with stream water to remove residual soil, soaked in 96% ethanol for 30 seconds to kill bacteria and halt decomposition, and dried at 60 °C for 12 h. Using a universal testing machine (Instron 5866 with 500 kN tensile holding clamps), maximum tensile strench of each cotton strip was measured. % tensile loss (proxy for decomposition) was calculated as (C-T) / C x 100, where T is the maximum tensile strength for each strip collected from the field, and C is the mean tensile strength of seven control strips, which had not been placed in the ground. See Robinson et al. 2021 for detailed description of plots sampled in 2015. Microbial abundance Bacterial and fungal abundance was estimated from additional soil cores of 3 cm in diameter taken from the upper 4 cm soil stratum (including the litter layer) at each plot. DNA was extracted using the PowerSoil DNA Isolation Kit (Qiagen, Germany). DNA was quantified using the high-sensitivity Qubit assay (Thermo Fisher Scientific, Switzerland). Relative abundances of bacterial and fungal communities were determined by quantitative PCR (qPCR) on an ABI7500 Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). PCR amplification of partial bacterial small-subunit ribosomal RNA genes (region V1–V3 of 16S; primers 27F and 512R) and fungal ribosomal internal transcribed spacers (region ITS2; primers IT3 and ITS4) was performed as described previously (Frey et al. 2020, Frey et al. 2021). For qPCR analyses, 2.5 ng DNA in a total volume of 6.6 µL and 8.4 µL GoTaq qPCRMaster Mix (Promega, Switzerland), containing 1.8 mM of each primer and 0.2 mg mL-1 of BSA, were used. The PCR conditions consisted of an initial denaturation at 95 ºC for 10 min, 40 cycles of denaturation at 95 ºC for 40 s, annealing at 58 ºC for 40 s and elongation at 72 ºC for 60 s followed by the final data acquisition step at 80 ºC for 60 s. The specificity of the amplification products was confirmed by melting-curve analysis. Three standard curves per target region (correlations ≥0.997) were obtained using tenfold serial dilutions (10-1 to 10-9 copies) of plasmids generated from cloned targets (Frey et al. 2020). Data were converted to represent the average copy number of targets per μg DNA and per g soil. Vegetation properties Vascular plant biomass was measured from a randomly placed 30 x 30 cm quadrat at each plot. To measure aboveground biomass (AGB) of plants, the aboveground layer of vegetation was cut and removed, dried at 70 °C for 24 h and weighed to obtain biomass per unit area. AGB was estimated as the biomass of graminoids plus forbs; total biomass of mosses was also estimated. Graminoid leaf N concentration was analysed from dried and ground leaf material using a LECO CNS-2000 analyser (LECO Corporation, Saint Joseph, MI, USA). Belowground biomass (BGB) of vascular plants was estimated from a soil core of 3 cm in diameter taken from the 10 cm upper soil stratum (excluding aboveground plant material) at each quadrat. Roots were extracted from the soil cores by rinsing in water using a 250-μm sieve, dried at 70 °C for 24 hours and weighed to obtain biomass per unit area. Root to shoot ratio was calculated as dry weight of BGB per cm2 divided by dry weight of AGB per cm2, and the total vascular plant biomass as the sum of AGB and BGB. Invertebrate community Enchytraied and nematode biomass was estimated from 3 soil cores of 3 cm in diameter taken from the upper 4 cm soil stratum (including litter layer) at each plot. Enchytraieds were extracted using wet funnels (O'Connor 1962) from a pooled sample of one half of each of the three soil cores, counted live, and classified into size classes (length 0-2, 2.1-4, 4.1-6, 6.1-8, 8.1-10, 10.1-12 or >12 mm) and their biomass was calculated according to Abrahamsen (1973). Nematodes were also extracted using wet funnels (Sohlenius 1979) from a pooled sample of a quarter of each of the three soil cores, counted live and preserved in 70% ethanol. Fifty individuals from each sample were identified and classified by trophic group (bacterivore, fungivoe, herbivore, omnivore, predator; Yeates et al. 1993). Soil micro-arthropods were extracted using a modified high-gradient-extractor (MacFayden 1961) from soil cores of 5.4 cm in diameter, taken from the upper 4 cm soil straum (including litter layer) at each plot. Total micro-arthropod biomass was calculated as the sum of all individual species' biomasses, obtained using length-weight regressions (see Robinson et al. 2021), and abundance of individual trophic groups (microbivore/detritivore, herbivore, omnivore, predator) calculated. Epigeal invertebrates were sampled by deploying five pitfall traps in each plot. White plastic cups of 7 cm in diameter and 8.5 cm in depth were filled with 10 ml of ethylene glycol and 30 ml of stream water, and left for 48 h before collection. Samples from the five traps at each plot were combined into a 250-μm sieve and stored in 70% ethanol. Invertebrate activity density (abundance) was estimate as the total number of individuals in the five traps, and total biomass as the sum of all individual species' biomasses. Invertebrates were identified to species level where possible and split into trophic groups, exluding adult Diptera, Hymenoptera, and Lepidoptera. Further details of sampling and collection of epigeal invertebrates are detailed in Robinson et al. (2018). References: Abrahamsen G. (1973) Studies on body-volume, body-surface area, density, and live weight of enchytraeidae (Oligochaeta). Pedobiologia 13: 6–15. Frey B, Carnol M, Dharmarajah A, Brunner I, Schleppi P. (2020) Only minor changes in the soil microbiome of a sub-alpine forest after 20 years of moderately increased nitrogen loads. Frontiers in Forests and Global Change 3: 77. Frey B, Walthert L, Perez-Mon C, Stierli B, Köchli R, Dharmarajah A, Brunner I (2021) Deep soil layers of drough-exposed forests harbor poorly known bacterial and fungal communities. Frontiers in Microbiology 12: 1061. MacFayden A. (1961) Improved funnel-type extractors for soil arthropods. Journal of Animal Ecology 30: 171–184. O’Connor FB. (1962) The extraction of Enchytraeidae from soil. In: P. W. Murphy (Ed.) Progress in soil zoology. Butterworth, London, UK; 279–285. Robinson SI, McLaughlin ÓB, Marteinsdóttir B, O'Gorman EJ. (2018) Soil temperature effects on the structure and diversity of plant and invertebrate communities in a natural warming experiment. Journal of Animal Ecology 87: 634–46. Robinson SI, Mikola J, Ovaskainen O, O’Gorman EJ. (2021) Temperature effects on the temporal dynamics of a subarctic invertebrate community. Journal of Animal Ecology 90: 1217-1227. Sohlenius B. (1979) A carbon budget for nematodes, rotifers and tardigrades in a Swedish coniferous forest soil. Holarctic Ecology 2: 30–40. Tiegs SD, Clapcott JE, Griffiths NA, Boulton AJ. (2013) A standardized cotton-strip assay for measuring organic-matter decomposition in streams. Ecological Indicators 32: 131–139. Yeates GW, Bongers T, De Goede RGM, Freckman DW, Georgieva SS. (1993) Feeding habits in soil nematode families and genera—an outline for soil ecologists. Journal of Nematology 25: 315–331. This is a dataset of soil physiochemical properties, bacterial and fungal abundance, and above and belowground plant and invertebrate biomass, sampled at 40 plots in the Hengill geothermal valley, Iceland, from 15th to 22nd August 2018. The plots span a temperature gradient of 10-35 °C over the sampling period, and this temperature gradient is consistent over time. The dataset also includes data on the decomposition rate of soil organic matter, which was sampled at 60 plots in the Hengill valley from May to July 2015. See README_Robinson_Hengill2018.txt